A generally applicable metadynamics scheme for predicting the free-energy profile of ligand binding to G-protein coupled receptors (GPCRs) is described. A common and effective collective variable (CV) has been defined using the ideally placed and highly conserved Trp6.48 as a reference point for ligand-GPCR distance measurement and the common orientation of GPCRs in the cell membrane. Using this single CV together with well-tempered multiple-walker 2 metadynamics with a funnel-like boundary allows an efficient exploration of the entire ligandbinding path from the extracellular medium to the orthosteric binding site, including vestibule and intermediate sites. The protocol can be used with X-ray structures or high-quality homology models for the receptor and is universally applicable to agonists, antagonists, partial and reverse agonists. The root mean square error (RMSE) in predicted binding free energies for 12 diverse ligands in five receptors (a total of 23 data points) is surprisingly small (less than 1 kcal mol − 1 ).The RMSEs for simulations that use receptor X-ray structures and homology models are very similar.

An index of the activation of Class A G-protein-coupled receptors (GPCRs) has been trained using interhelix distances from a series of microsecond molecular-dynamics simulations and tested for 268 published X-ray structures. In a three-class model that includes intermediate structures, 63% of the active structures are classified in agreement with the experimental assignment, 81% of the intermediate structures, and 89% of the inactives. An alternative two-state model classifies 94% of the actives and 99% of the inactives correctly. The intermediate structures are distributed 2:1 between actives and inactives. X-ray structures with protein nanobodies give good agreement between the assigned activation state and the predictions of the model, whereby many active nanobody structures are predicted to be weakly active. The five interhelix C α −C α distances that occur in the model relate clearly to the established activation mechanism. The model is available as a Python script or via an interactive web page. It can thus be used to classify both experimental and computational GPCR structures.

No abstract

Protein nanobodies have been used successfully as surrogates for unstable G-proteins in order to crystallize G-protein-coupled receptors (GPCRs) in their active states. We used molecular dynamics (MD) simulations, including metadynamics enhanced sampling, to investigate the similarities and differences between GPCR-agonist ternary complexes with the α-subunits of the appropriate G-proteins and those with the protein nanobodies (intracellular binding partners, IBPs) used for crystallization. In two of the three receptors considered, the agonist-binding mode differs significantly between the two alternative ternary complexes. The ternary-complex model of GPCR activation entails enhancement of ligand binding by bound IBPs: Our results show that IBP-specific changes can alter the agonist binding modes and thus also the criteria for designing GPCR agonists.

We present a generally applicable metadynamics protocol for characterizing the activation free-energy profiles of class A G-protein coupled receptors and a proof-of-principle study for the 5HT1A-receptor. The almost universal A100 activation index, which depends on five inter-helix distances, is used as the single collective variable in well-tempered multiple-walker metadynamics simulations. Here, we show free-energy profiles for the serotonin receptor as binary (apo-receptor + G-protein-α-subunit and receptor + ligand) and ternary complexes with two prototypical orthosteric ligands: the full agonist serotonin and the partial agonist aripiprazole. Our results are not only compatible with previously reported experimental and computational data, but they also allow differences between active and inactive conformations to be determined in unprecedented atomic detail, and with respect to the so-called microswitches that have been suggested as determinants of activation, giving insight into their role in the activation mechanism.