The dynamic process of ß(2)-adrenergic receptor activation.

G-protein-coupled receptors (GPCRs) can modulate diverse signaling pathways, often in a ligand-specific manner. The full range of functionally relevant GPCR conformations is poorly understood. Here, we use NMR spectroscopy to characterize the conformational dynamics of the transmembrane core of the ß(2)-adrenergic receptor (ß(2)AR), a prototypical GPCR. We labeled ß(2)AR with (13)CH(3)ε-methionine and obtained HSQC spectra of unliganded receptor as well as receptor bound to an inverse agonist, an agonist, and a G-protein-mimetic nanobody. These studies provide evidence for conformational states not observed in crystal structures, as well as substantial conformational heterogeneity in agonist- and inverse-agonist-bound preparations. They also show that for ß(2)AR, unlike rhodopsin, an agonist alone does not stabilize a fully active conformation, suggesting that the conformational link between the agonist-binding pocket and the G-protein-coupling surface is not rigid. The observed heterogeneity may be important for ß(2)AR's ability to engage multiple signaling and regulatory proteins.

Agonist-independent GPCR activity regulates anterior-posterior targeting of olfactory sensory neurons.

G-protein-coupled receptors (GPCRs) are known to possess two different conformations, active and inactive, and they spontaneously alternate between the two in the absence of ligands. Here, we analyzed the agonist-independent GPCR activity for its possible role in receptor-instructed axonal projection. We generated transgenic mice expressing activity mutants of the ß2-adrenergic receptor, a well-characterized GPCR with the highest homology to odorant receptors (ORs). We found that mutants with altered agonist-independent activity changed the transcription levels of axon-targeting molecules--e.g., Neuropilin-1 and Plexin-A1--but not of glomerular segregation molecules--e.g., Kirrel2 and Kirrel3--thus causing shifts in glomerular locations along the anterior-posterior (A-P) axis. Knockout and in vitro experiments demonstrated that Gs, but not Golf, is responsible for mediating the agonist-independent GPCR activity. We conclude that the equilibrium of conformational transitions set by each OR is the major determinant of expression levels of A-P-targeting molecules.

Single-molecule analysis of ligand efficacy in ß2AR-G-protein activation.

G-protein-coupled receptor (GPCR)-mediated signal transduction is central to human physiology and disease intervention, yet the molecular mechanisms responsible for ligand-dependent signalling responses remain poorly understood. In class A GPCRs, receptor activation and G-protein coupling entail outward movements of transmembrane helix 6 (TM6). Here, using single-molecule fluorescence resonance energy transfer imaging, we examine TM6 movements in the ß2 adrenergic receptor (ß2AR) upon exposure to orthosteric ligands with different efficacies, in the absence and presence of the Gs heterotrimer. We show that partial and full agonists differentially affect TM6 motions to regulate the rate at which GDP-bound ß2AR-Gs complexes are formed and the efficiency of nucleotide exchange leading to Gs activation. These data also reveal transient nucleotide-bound ß2AR-Gs species that are distinct from known structures, and provide single-molecule perspectives on the allosteric link between ligand- and nucleotide-binding pockets that shed new light on the G-protein activation mechanism.

Allosteric coupling from G protein to the agonist-binding pocket in GPCRs.

G-protein-coupled receptors (GPCRs) remain the primary conduit by which cells detect environmental stimuli and communicate with each other. Upon activation by extracellular agonists, these seven-transmembrane-domain-containing receptors interact with heterotrimeric G proteins to regulate downstream second messenger and/or protein kinase cascades. Crystallographic evidence from a prototypic GPCR, the ß2-adrenergic receptor (ß2AR), in complex with its cognate G protein, Gs, has provided a model for how agonist binding promotes conformational changes that propagate through the GPCR and into the nucleotide-binding pocket of the G protein α-subunit to catalyse GDP release, the key step required for GTP binding and activation of G proteins. The structure also offers hints about how G-protein binding may, in turn, allosterically influence ligand binding. Here we provide functional evidence that G-protein coupling to the ß2AR stabilizes a 'closed' receptor conformation characterized by restricted access to and egress from the hormone-binding site. Surprisingly, the effects of G protein on the hormone-binding site can be observed in the absence of a bound agonist, where G-protein coupling driven by basal receptor activity impedes the association of agonists, partial agonists, antagonists and inverse agonists. The ability of bound ligands to dissociate from the receptor is also hindered, providing a structural explanation for the G-protein-mediated enhancement of agonist affinity, which has been observed for many GPCR­G-protein pairs. Our data also indicate that, in contrast to agonist binding alone, coupling of a G protein in the absence of an agonist stabilizes large structural changes in a GPCR. The effects of nucleotide-free G protein on ligand-binding kinetics are shared by other members of the superfamily of GPCRs, suggesting that a common mechanism may underlie G-protein-mediated enhancement of agonist affinity.

Membrane-Facilitated Receptor Access and Binding Mechanisms of Long-Acting ß2-Adrenergic Receptor Agonists.

The drugs salmeterol, formoterol, and salbutamol constitute the frontline treatment of asthma and other chronic pulmonary diseases. These drugs activate the ß2-adrenergic receptors (ß2-AR), a class A G protein-coupled receptor (GPCR), and differ significantly in their clinical onset and duration of actions. According to the microkinetic model, the long duration of action of salmeterol and formoterol compared with salbutamol were attributed, at least in part, to their high lipophilicity and increased local concentrations in the membrane near the receptor. However, the structural and molecular bases of how the lipophilic drugs reach the binding site of the receptor from the surrounding membrane remain unknown. Using a variety of classic and enhanced molecular dynamics simulation techniques, we investigated the membrane partitioning characteristics, binding, and unbinding mechanisms of the ligands. The obtained results offer remarkable insight into the functional role of membrane lipids in the ligand association process. Strikingly, salmeterol entered the binding site from the bilayer through transmembrane helices 1 and 7. The entry was preceded by membrane-facilitated rearrangement and presentation of its phenyl-alkoxy-alkyl tail as a passkey to an access route gated by F193, a residue known to be critical for salmeterol's affinity. Formoterol's access is through the aqueous path shared by other ß2-AR agents. We observed a novel secondary path for salbutamol that is distinct from its primary route. Our study offers a mechanistic description for the membrane-facilitated access and binding of ligands to a membrane protein and establishes a groundwork for recognizing membrane lipids as an integral component in the molecular recognition process. SIGNIFICANCE STATEMENT: The cell membrane's functional role behind the duration of action of long-acting ß2-adrenergic receptor (ß2-AR) agonists such as salmeterol has been a subject of debate for a long time. This study investigated the binding and unbinding mechanisms of the three commonly used ß2-AR agonists, salmeterol, formoterol, and salbutamol, using advanced simulation techniques. The obtained results offer unprecedented insights into the active role of membrane lipids in facilitating access and binding of the ligands, affecting the molecular recognition process and thus their pharmacology.

Energy Landscapes Reveal Agonist Control of G Protein-Coupled Receptor Activation via Microswitches.

Agonist binding to G protein-coupled receptors (GPCRs) leads to conformational changes in the transmembrane region that activate cytosolic signaling pathways. Although high-resolution structures of different receptor states are available, atomistic details of allosteric signaling across the membrane remain elusive. We calculated free energy landscapes of ß2 adrenergic receptor activation using atomistic molecular dynamics simulations in an optimized string of swarms framework, which shed new light on how microswitches govern the equilibrium between conformational states. Contraction of the extracellular binding site in the presence of the agonist BI-167107 is obligatorily coupled to conformational changes in a connector motif located in the core of the transmembrane region. The connector is probabilistically coupled to the conformation of the intracellular region. An active connector promotes desolvation of a buried cavity, a twist of the conserved NPxxY motif, and an interaction between two conserved tyrosines in transmembrane helices 5 and 7 (Y-Y motif), which lead to a larger population of active-like states at the G protein binding site. This coupling is augmented by protonation of the strongly conserved Asp792.50. The agonist binding site hence communicates with the intracellular region via a cascade of locally connected microswitches. Characterization of these can be used to understand how ligands stabilize distinct receptor states and contribute to development drugs with specific signaling properties. The developed simulation protocol can likely be transferred to other class A GPCRs.

Comparative Docking to Distinct G Protein-Coupled Receptor Conformations Exclusively Yields Ligands with Agonist Efficacy.

G protein-coupled receptors exist in a whole spectrum of conformations that are stabilized by the binding of ligands with different efficacy or intracellular effector proteins. Here, we investigate whether three-dimensional structures of receptor conformations in different states of activation can be used to enrich ligands with agonist behavior in prospective docking calculations. We focused on the ß 2-adrenergic receptor, as it is currently the receptor with the highest number of active-state crystal structures. Comparative docking calculations to distinct conformations of the receptor were used for the in silico prediction of ligands with agonist efficacy. The pharmacology of molecules selected based on these predictions was characterized experimentally, resulting in a hit rate of 37% ligands, all of which were agonists. The ligands furthermore contain a pyrazole moiety that has previously not been described for ß 2-adrenergic receptor ligands, and one of them shows an intrinsic efficacy comparable to salbutamol. SIGNIFICANCE STATEMENT: Structure-based ligand design for G protein-coupled receptors crucially depends on receptor conformation and, hence, their activation state. We explored the influence of using multiple active-conformation X-ray structures on the hit rate of docking calculations to find novel agonists, and how to predict the most fruitful strategy to apply. The results suggest that aggregating the ranks of molecules across docking calculations to more than one active-state structure exclusively yields agonists.

Cell-based assays and molecular dynamics analysis of a boron-containing agonist with different profiles of binding to human and guinea pig beta2 adrenoceptors.

The design of beta2 adrenoceptor (ß2AR) agonists is attractive because of their wide-ranging applications in medicine, and the details of agonist interactions with ß2AR are interesting because it is considered a prototype for G-protein coupled receptors. Preclinical studies for agonist development have involved biological assays with guinea pigs due to a similar physiology to humans. Boron-containing Albuterol derivatives (BCADs) designed as bronchodilators have improved potency and efficacy compared with their boron-free precursor on guinea pig ß2ARs (gpß2ARs), and two of the BCADs (BR-AEA and boronterol) conserve these features on cells expressing human ß2ARs (hß2ARs). The aim of this study was to test the BCAD Politerol on gpß2ARs and hß2ARs in vitro and in silico. Politerol displayed higher potency and efficacy on gpß2AR than on hß2AR in experimental assays, possible explanations are provided based on molecular modeling, and molecular dynamics simulations of about 0.25 µs were performed for the free and bound states adding up to 2 µs in total. There were slight differences, particularly in the role of the boron atom, in the interactions of Politerol with gpß2ARs and hß2ARs, affecting movements of transmembrane domains 5-7, known to be pivotal in receptor activation. These findings could be instrumental in the design of compounds selective for hß2ARs.

Design, synthesis and biological evaluation of 5-(2-amino-1-hydroxyethyl)-8-hydroxyquinolin-2(1H)-one derivatives as potent ß2-adrenoceptor agonists.

A series of novel ß2-adrenoceptor agonists with a 5-(2-amino-1-hydroxyethyl)-8-hydroxyquinolin-2(1H)-one moiety was designed, synthesized and evaluated for biological activity in human embryonic kidney 293 cells and isolated guinea pig trachea. Compounds 9g and (R)-18c exhibited the most excellent ß2-adrenoceptor agonistic effects and high ß2/ß1-selectivity with EC50 values of 36 pM for 9g and 21 pM for (R)-18c. They produced potent airway smooth muscle relaxant effects with fast onset of action and long duration of action in an in vitro guinea pig trachea model of bronchodilation. These results support further development of the two compounds into drug candidates.

Long-Acting ß2-Adrenoceptor Agonists Enhance Glucocorticoid Receptor (GR)-Mediated Transcription by Gene-Specific Mechanisms Rather Than Generic Effects via GR.

In asthma, the clinical efficacy of inhaled corticosteroids (ICSs) is enhanced by long-acting ß2-adrenoceptor agonists (LABAs). ICSs, or more accurately, glucocorticoids, promote therapeutically relevant changes in gene expression, and, in primary human bronchial epithelial cells (pHBECs) and airway smooth muscle cells, this genomic effect can be enhanced by a LABA. Modeling this interaction in human bronchial airway epithelial BEAS-2B cells transfected with a 2× glucocorticoid response element (2×GRE)-driven luciferase reporter showed glucocorticoid-induced transcription to be enhanced 2- to 3-fold by LABA. This glucocorticoid receptor (GR; NR3C1)-dependent effect occurred rapidly, was insensitive to protein synthesis inhibition, and was maximal when glucocorticoid and LABA were added concurrently. The ability of LABA to enhance GR-mediated transcription was not associated with changes in GR expression, serine (Ser203, Ser211, Ser226) phosphorylation, ligand affinity, or nuclear translocation. Chromatin immunoprecipitation demonstrated that glucocorticoid-induced recruitment of GR to the integrated 2×GRE reporter and multiple gene loci, whose mRNAs were unaffected or enhanced by LABA, was also unchanged by LABA. Transcriptomic analysis revealed glucocorticoid-induced mRNAs were variably enhanced, unaffected, or repressed by LABA. Thus, events leading to GR binding at target genes are not the primary explanation for how LABAs modulate GR-mediated transcription. As many glucocorticoid-induced genes are independently induced by LABA, gene-specific control by GR- and LABA-activated transcription factors may explain these observations. Because LABAs promote similar effects in pHBECs, therapeutic relevance is likely. These data illustrate the need to understand gene function(s), and the mechanisms leading to gene-specific induction, if existing ICS/LABA combination therapies are to be improved.

Identifying Functional Hotspot Residues for Biased Ligand Design in G-Protein-Coupled Receptors.

G-protein-coupled receptors (GPCRs) mediate multiple signaling pathways in the cell, depending on the agonist that activates the receptor and multiple cellular factors. Agonists that show higher potency to specific signaling pathways over others are known as "biased agonists" and have been shown to have better therapeutic index. Although biased agonists are desirable, their design poses several challenges to date. The number of assays to identify biased agonists seems expensive and tedious. Therefore, computational methods that can reliably calculate the possible bias of various ligands ahead of experiments and provide guidance, will be both cost and time effective. In this work, using the mechanism of allosteric communication from the extracellular region to the intracellular transducer protein coupling region in GPCRs, we have developed a computational method to calculate ligand bias ahead of experiments. We have validated the method for several ß-arrestin-biased agonists in ß2-adrenergic receptor (ß2AR), serotonin receptors 5-HT1B and 5-HT2B and for G-protein-biased agonists in the κ-opioid receptor. Using this computational method, we also performed a blind prediction followed by experimental testing and showed that the agonist carmoterol is ß-arrestin-biased in ß2AR. Additionally, we have identified amino acid residues in the biased agonist binding site in both ß2AR and κ-opioid receptors that are involved in potentiating the ligand bias. We call these residues functional hotspots, and they can be used to derive pharmacophores to design biased agonists in GPCRs.

Small-Molecule Positive Allosteric Modulators of the ß2-Adrenoceptor Isolated from DNA-Encoded Libraries.

Conventional drug discovery efforts at the ß2-adrenoceptor (ß2AR) have led to the development of ligands that bind almost exclusively to the receptor's hormone-binding orthosteric site. However, targeting the largely unexplored and evolutionarily unique allosteric sites has potential for developing more specific drugs with fewer side effects than orthosteric ligands. Using our recently developed approach for screening G protein-coupled receptors (GPCRs) with DNA-encoded small-molecule libraries, we have discovered and characterized the first ß2AR small-molecule positive allosteric modulators (PAMs)-compound (Cmpd)-6 [(R)-N-(4-amino-1-(4-(tert-butyl)phenyl)-4-oxobutan-2-yl)-5-(N-isopropyl-N-methylsulfamoyl)-2-((4-methoxyphenyl)thio)benzamide] and its analogs. We used purified human ß2ARs, occupied by a high-affinity agonist, for the affinity-based screening of over 500 million distinct library compounds, which yielded Cmpd-6. It exhibits a low micro-molar affinity for the agonist-occupied ß2AR and displays positive cooperativity with orthosteric agonists, thereby enhancing their binding to the receptor and ability to stabilize its active state. Cmpd-6 is cooperative with G protein and ß-arrestin1 (a.k.a. arrestin2) to stabilize high-affinity, agonist-bound active states of the ß2AR and potentiates downstream cAMP production and receptor recruitment of ß-arrestin2 (a.k.a. arrestin3). Cmpd-6 is specific for the ß2AR compared with the closely related ß1AR. Structure-activity studies of select Cmpd-6 analogs defined the chemical groups that are critical for its biologic activity. We thus introduce the first small-molecule PAMs for the ß2AR, which may serve as a lead molecule for the development of novel therapeutics. The approach described in this work establishes a broadly applicable proof-of-concept strategy for affinity-based discovery of small-molecule allosteric compounds targeting unique conformational states of GPCRs.

Allosteric regulation of G protein-coupled receptor activity by phospholipids.

Lipids are emerging as key regulators of membrane protein structure and activity. These effects can be attributed either to the modification of bilayer properties (thickness, curvature and surface tension) or to the binding of specific lipids to the protein surface. For G protein-coupled receptors (GPCRs), the effects of phospholipids on receptor structure and activity remain poorly understood. Here we reconstituted purified ß2-adrenergic receptor (ß2R) in high-density lipoparticles to systematically characterize the effect of biologically relevant phospholipids on receptor activity. We observed that the lipid headgroup type affected ligand binding (agonist and antagonist) and receptor activation. Specifically, phosphatidylgycerol markedly favored agonist binding and facilitated receptor activation, whereas phosphatidylethanolamine favored antagonist binding and stabilized the inactive state of the receptor. We then showed that these effects could be recapitulated with detergent-solubilized lipids, demonstrating that the functional modulation occurred in the absence of a bilayer. Our data suggest that phospholipids act as direct allosteric modulators of GPCR activity.

Catecholic Isoquinolines from Portulaca oleracea and Their Anti-inflammatory and ß2-Adrenergic Receptor Agonist Activity.

Isoquinoline alkaloids possess a wide range of structural features and pharmaceutical activities and are promising drug candidates. Ten water-soluble catecholic isoquinolines were isolated from the medicinal plant Portulaca oleracea, including three new (1-3) and seven known compounds (4-10), along with the known catecholamines 11 and 12 and four other known compounds (13-16). A method of polyamide column chromatography using EtOAc-MeOH as the mobile phase was developed for the isolation of catecholic isoquinolines. Alkaloids 1-12 exhibited anti-inflammatory activities (EC50 = 18.0-497.7 µM) through inhibition of NO production in lipopolysaccharide-induced murine macrophage RAW 264.7 cells. Among these compounds, 11, 2, 5, 4, and 8 were more potent than was the positive control, 3,4-dihydroxybenzohydroxamic acid (EC50 = 82.4 µM), with EC50 values of 18.0, 18.1, 35.4, 36.3, and 58.7 µM, respectively. Additionally, at 100 µM, compounds 1-12 showed different degrees of ß2-adrenergic receptor (ß2-AR) agonist activity in the CHO-K1/GA15 cell line which stably expressed ß2-AR as detected by a calcium assay. The EC50 values of 2 and 10 were 5.1 µM and 87.9 nM, respectively.

Crystal structure of the ß2 adrenergic receptor-Gs protein complex.

G protein-coupled receptors (GPCRs) are responsible for the majority of cellular responses to hormones and neurotransmitters as well as the senses of sight, olfaction and taste. The paradigm of GPCR signalling is the activation of a heterotrimeric GTP binding protein (G protein) by an agonist-occupied receptor. The ß(2) adrenergic receptor (ß(2)AR) activation of Gs, the stimulatory G protein for adenylyl cyclase, has long been a model system for GPCR signalling. Here we present the crystal structure of the active state ternary complex composed of agonist-occupied monomeric ß(2)AR and nucleotide-free Gs heterotrimer. The principal interactions between the ß(2)AR and Gs involve the amino- and carboxy-terminal α-helices of Gs, with conformational changes propagating to the nucleotide-binding pocket. The largest conformational changes in the ß(2)AR include a 14 Å outward movement at the cytoplasmic end of transmembrane segment 6 (TM6) and an α-helical extension of the cytoplasmic end of TM5. The most surprising observation is a major displacement of the α-helical domain of Gαs relative to the Ras-like GTPase domain. This crystal structure represents the first high-resolution view of transmembrane signalling by a GPCR.

Conformational changes in the G protein Gs induced by the ß2 adrenergic receptor.

G protein-coupled receptors represent the largest family of membrane receptors that instigate signalling through nucleotide exchange on heterotrimeric G proteins. Nucleotide exchange, or more precisely, GDP dissociation from the G protein α-subunit, is the key step towards G protein activation and initiation of downstream signalling cascades. Despite a wealth of biochemical and biophysical studies on inactive and active conformations of several heterotrimeric G proteins, the molecular underpinnings of G protein activation remain elusive. To characterize this mechanism, we applied peptide amide hydrogen-deuterium exchange mass spectrometry to probe changes in the structure of the heterotrimeric bovine G protein, Gs (the stimulatory G protein for adenylyl cyclase) on formation of a complex with agonist-bound human ß(2) adrenergic receptor (ß(2)AR). Here we report structural links between the receptor-binding surface and the nucleotide-binding pocket of Gs that undergo higher levels of hydrogen-deuterium exchange than would be predicted from the crystal structure of the ß(2)AR-Gs complex. Together with X-ray crystallographic and electron microscopic data of the ß(2)AR-Gs complex (from refs 2, 3), we provide a rationale for a mechanism of nucleotide exchange, whereby the receptor perturbs the structure of the amino-terminal region of the α-subunit of Gs and consequently alters the 'P-loop' that binds the ß-phosphate in GDP. As with the Ras family of small-molecular-weight G proteins, P-loop stabilization and ß-phosphate coordination are key determinants of GDP (and GTP) binding affinity.

Structure and function of an irreversible agonist-ß(2) adrenoceptor complex.

G-protein-coupled receptors (GPCRs) are eukaryotic integral membrane proteins that modulate biological function by initiating cellular signalling in response to chemically diverse agonists. Despite recent progress in the structural biology of GPCRs, the molecular basis for agonist binding and allosteric modulation of these proteins is poorly understood. Structural knowledge of agonist-bound states is essential for deciphering the mechanism of receptor activation, and for structure-guided design and optimization of ligands. However, the crystallization of agonist-bound GPCRs has been hampered by modest affinities and rapid off-rates of available agonists. Using the inactive structure of the human ß(2) adrenergic receptor (ß(2)AR) as a guide, we designed a ß(2)AR agonist that can be covalently tethered to a specific site on the receptor through a disulphide bond. The covalent ß(2)AR-agonist complex forms efficiently, and is capable of activating a heterotrimeric G protein. We crystallized a covalent agonist-bound ß(2)AR-T4L fusion protein in lipid bilayers through the use of the lipidic mesophase method, and determined its structure at 3.5 Å resolution. A comparison to the inactive structure and an antibody-stabilized active structure (companion paper) shows how binding events at both the extracellular and intracellular surfaces are required to stabilize an active conformation of the receptor. The structures are in agreement with long-timescale (up to 30 µs) molecular dynamics simulations showing that an agonist-bound active conformation spontaneously relaxes to an inactive-like conformation in the absence of a G protein or stabilizing antibody.

Structure of a nanobody-stabilized active state of the ß(2) adrenoceptor.

G protein coupled receptors (GPCRs) exhibit a spectrum of functional behaviours in response to natural and synthetic ligands. Recent crystal structures provide insights into inactive states of several GPCRs. Efforts to obtain an agonist-bound active-state GPCR structure have proven difficult due to the inherent instability of this state in the absence of a G protein. We generated a camelid antibody fragment (nanobody) to the human ß(2) adrenergic receptor (ß(2)AR) that exhibits G protein-like behaviour, and obtained an agonist-bound, active-state crystal structure of the receptor-nanobody complex. Comparison with the inactive ß(2)AR structure reveals subtle changes in the binding pocket; however, these small changes are associated with an 11 Å outward movement of the cytoplasmic end of transmembrane segment 6, and rearrangements of transmembrane segments 5 and 7 that are remarkably similar to those observed in opsin, an active form of rhodopsin. This structure provides insights into the process of agonist binding and activation.

Differences between G-Protein-Stabilized Agonist-GPCR Complexes and their Nanobody-Stabilized Equivalents.

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.

Conformational and Thermodynamic Landscape of GPCR Activation from Theory and Computation.

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.