Abstract
The activation of G protein-coupled receptors (GPCRs) is a complex multibody multievent process involving agonist binding, receptor activation, G protein coupling, and subsequent G protein activation. The order and energetics of these events, though crucial for the rational design of selective GPCR drugs, are challenging to characterize and remain largely underexplored. Here, we employed molecular dynamics simulations and our recently developed traveling salesman-based automated path searching (TAPS) algorithm to efficiently locate the minimum free-energy paths for the coupling events of the CXC chemokine receptor 4 (CXCR4) with its endogenous ligand CXC chemokine ligand 12 (CXCL12) and Gi protein. We show that, after overcoming three low energy barriers (3.24-6.89 kcal/mol), Gi alone can precouple with CXCR4 even without CXCL12, consistent with previous reports on the existence of the apo CXCR4-Gi complex and our NanoBiT experiments. The highest barrier of 6.89 kcal/mol in this process corresponds to the packing of the proline-isoleucine-phenylalanine (PIF) motif of CXCR4. Interestingly, without Gi, CXCL12 alone cannot activate CXCR4 (high barrier of 18.89 kcal/mol). Instead, it can enhance Gi coupling by circumventing the energy barrier of PIF packing. Shedding new light on the activation mechanism of CXCR4, these results present TAPS as a promising tool for uncovering complete activation pathways of GPCRs and the corresponding agonist design.