Abstract:We present a hybrid computational methodology to predict multiple energetically accessible conformations for G protein-coupled receptors (GPCRs) that might play a role in binding to ligands and different signaling partners. To our knowledge, this method, termed ActiveGEnSeMBLE, enables the first quantitative energy profile for GPCR activation that is consistent with the qualitative profile deduced from experiments. ActiveGEnSeMBLE starts with a systematic coarse grid sampling of helix tilts/rotations (∼13 trillion transmembrane-domain conformations) and selects the conformational landscape based on energy. This profile identifies multiple potential active-state energy wells, with the TM3-TM6 intracellular distance as an approximate activation coordinate. These energy wells are then sampled locally using a finer grid to find locally minimized conformation in each energy well. We validate this strategy using the inactive and active experimental structures of β2 adrenergic receptor (hβ2AR) and M2 muscarinic acetylcholine receptor. Structures of membrane-embedded hβ2AR along its activation coordinate are subjected to molecular-dynamics simulations for relaxation and interaction energy analysis to generate a quantitative energy landscape for hβ2AR activation. This landscape reveals several metastable states along this coordinate, indicating that for hβ2AR, the agonist alone is not enough to stabilize the active state and that the G protein is necessary, consistent with experimental observations. The method's application to somatostatin receptor SSTR5 (no experimental structure available) shows that to predict an active conformation it is better to start from an inactive structure template based on a close homolog than to start from an active template based on a distant homolog. The energy landscape for hSSTR5 activation is consistent with hβ2AR in the role of the G protein. These results demonstrate the utility of the ActiveGEnSeMBLE method for predicting multiple conformations along the pathways for activating GPCRs and the corresponding energy landscapes, thereby providing detailed structural insights into the initial molecular events of GPCR function that are not easily accessible by experiments.

The Cubic Ternary Complex ReceptorOccupancy Model I. Model Description

The Evolution of Drug-Receptor Models: The Cubic Ternary Complex Model for G Protein-Coupled Receptors

Theoretical Aspects of GPCR-Ligand Complex Pharmacology

A Thermodynamic Model for Receptor Clustering

A 19F-qNMR-Guided Mathematical Model for G Protein-Coupled Receptor Signaling

Toward a Structural Understanding of Class B GPCR Peptide Binding and Activation.

Do pharmacological methods for the quantification of agonists work when the ternary complex mechanism operates?

Derangement model of ligand-receptor binding

A mutation-induced activated state of the beta 2-adrenergic receptor. Extending the ternary complex model.

Application of an Allosteric Ternary Complex Model to the Technique of Pharmacological Resultant Analysis

Receptor Surface Models. 2. Application to Quantitative Structure-Activity Relationships Studies

Cosmological Solutions of Type II Superstrings

Multiple Ligand-Specific Conformations of the Β2-Adrenergic Receptor

Molecular mechanisms for the assembly of the T cell receptor–CD3 complex

Experimental and Computational Models for Side Chain Discrimination in Peptide-Protein Interactions

Structure of a Class C GPCR Metabotropic Glutamate Receptor 1 Bound to an Allosteric Modulator

Seven-helix receptors: structure and modelling

Monitoring the dynamics of ligand-receptor complexes on model membranes

Conformational and Thermodynamic Landscape of GPCR Activation from Theory and Computation

Control of G protein-coupled receptor function via membrane-interacting intrinsically disordered C-terminal domains

Mechanism of psychoactive drug action in the brain: simulation modeling of GABAA receptor interactions at non-equilibrium conditions