Helix 3 acts as a conformational hinge in Class A GPCR activation: An analysis of interhelical interaction energies in crystal structures

Magic angle spinning NMR of G protein-coupled receptors

G protein-coupled receptors (GPCRs) have a simple seven transmembrane helix architecture which has evolved to recognize a diverse number of chemical signals. The more than 800 GPCRs encoded in the human genome function as receptors for vision, smell and taste, and mediate key physiological processes. Consequently, these receptors are a major target for pharmaceuticals. Protein crystallography and electron cryo-microscopy have provided high resolution structures of many GPCRs in both active and inactive conformations. However, these structures have not sparked a surge in rational drug design, in part because GPCRs are inherently dynamic and the structural changes induced by ligand or drug binding to stabilize inactive or active conformations are often subtle rearrangements in packing or hydrogen-bonding interactions. NMR spectroscopy provides a sensitive probe of local structure and dynamics at specific sites within these receptors as well as global changes in receptor structure and dynamics. These methods can also capture intermediate states and conformations with low populations that provide insights into the activation pathways. We review the use of solid-state magic angle spinning NMR to address the structure and activation mechanisms of GPCRs. The focus is on the large and diverse class A family of receptors. We highlight three specific class A GPCRs in order to illustrate how solid-state, as well as solution-state, NMR spectroscopy can answer questions in the field involving how different GPCR classes and subfamilies are activated by their associated ligands, and how small molecule drugs can modulate GPCR activation.

Deconstructing the transmembrane core of class A G protein-coupled receptors

Class A G protein-coupled receptors have evolved to recognize ligands ranging from small-molecule odorants to proteins. Although they are among the most diverse membrane receptors in eukaryotic organisms, they possess a highly conserved core within their seven-transmembrane helix framework. The conservation of the transmembrane core has led to the idea of a common mechanism by which ligand binding is coupled to the outward rotation of helix H6, the hallmark of an active receptor. Nevertheless, there is still no consensus on the mechanism of coupling or on the roles of specific residues within the core. Recent insights from crystallography and NMR spectroscopy provide a way to decompose the core into its essential structural and functional elements that shed new light on this important region.

Structural Assessment of Agonist Efficacy in the μ-Opioid Receptor: Morphine and Fentanyl Elicit Different Activation Patterns

Over the past two decades, the opioid epidemic in the United States and Canada has evidenced the need for a better understanding of the molecular mechanisms of medications used to fight pain. Morphine and fentanyl are widely used in opiate-mediated analgesia for the treatment of chronic pain. These compounds target the μ-opioid receptor (MOR), a class A G protein-coupled receptor (GPCR). In light of described higher efficacy of fentanyl with respect to morphine, we have performed independent μs-length unbiased molecular dynamics (MD) simulations of MOR complexes with each of these ligands, including the MOR antagonist naltrexone as a negative control. Consequently, MD simulations totaling 58 μs have been conducted to elucidate at the atomic level ligand-specific receptor activity and signal transmission in the MOR. In particular, we have identified stable binding poses of morphine and fentanyl, which interact differently with the MOR. Different ligand-receptor interaction landscapes directly induce sidechain conformational changes of orthosteric pocket residues: Asp149^3.32, Tyr150^3.33, Gln126^2.60, and Lys235^5.39. The induced conformations determine Asp149^3.32-Tyr328^7.43 sidechain-sidechain interactions and Trp295^6.48-Ala242^5.46 sidechain-backbone H-bond formations, as well as Met153^3.36 conformational changes. In addition to differences in ligand binding, different intracellular receptor conformational changes are observed as morphine preferentially activates transmembrane (TM) helices: TM3 and TM5, while fentanyl preferentially activates TM6 and TM7. As conformational changes in TM6 and TM7 are widely described as being the most crucial aspect in GPCR activation, this may contribute to the greater efficacy of fentanyl over morphine. These computationally observed functional differences between fentanyl and morphine may provide new avenues for the design of safer but not weaker opioid drugs because it is desirable to increase the safety of medicines without sacrificing their efficacy.

Exploring the activation pathway and Gi-coupling specificity of the μ-opioid receptor

Understanding the activation mechanism of the μ-opioid receptor (μ-OR) and its selective coupling to the inhibitory G protein (Gi) is vital for pharmaceutical research aimed at finding treatments for the opioid overdose crisis. Many attempts have been made to understand the mechanism of the μ-OR activation, following the elucidation of new crystal structures such as the antagonist- and agonist-bound μ-OR. However, the focus has not been placed on the underlying energetics and specificity of the activation process. An energy-based picture would not only help to explain this coupling but also help to explore why other possible options are not common. For example, one would like to understand why μ-OR is more selective to Gi than a stimulatory G protein (Gs). Our study used homology modeling and a coarse-grained model to generate all of the possible "end states" of the thermodynamic cycle of the activation of μ-OR. The end points were further used to generate reasonable intermediate structures of the receptor and the Gi to calculate two-dimensional free energy landscapes. The results of the landscape calculations helped to propose a plausible sequence of conformational changes in the μ-OR and Gi system and for exploring the path that leads to its activation. Furthermore, in silico alanine scanning calculations of the last 21 residues of the C terminals of Gi and Gs were performed to shed light on the selective binding of Gi to μ-OR. Overall, the present work appears to demonstrate the potential of multiscale modeling in exploring the action of G protein-coupled receptors.

Insights into adenosine A2A receptor activation through cooperative modulation of agonist and allosteric lipid interactions

The activation process of G protein-coupled receptors (GPCRs) has been extensively studied, both experimentally and computationally. In particular, Molecular Dynamics (MD) simulations have proven useful in exploring GPCR conformational space. The typical behaviour of class A GPCRs, when subjected to unbiased MD simulations from their crystallized inactive state, is to fluctuate between inactive and intermediate(s) conformations, even with bound agonist. Fully active conformation(s) are rarely stabilized unless a G protein is also bound. Despite several crystal structures of the adenosine A2a receptor (A2aR) having been resolved in complex with co-crystallized agonists and Gs protein, its agonist-mediated activation process is still not completely understood. In order to thoroughly examine the conformational landscape of A2aR activation, we performed unbiased microsecond-length MD simulations in quadruplicate, starting from the inactive conformation either in apo or with bound agonists: endogenous adenosine or synthetic NECA, embedded in two homogeneous phospholipid membranes: 1,2-dioleoyl-sn-glycerol-3-phosphoglycerol (DOPG) or 1,2-dioleoyl-sn-glycerol-3-phosphocholine (DOPC). In DOPC with bound adenosine or NECA, we observe transition to an intermediate receptor conformation consistent with the known adenosine-bound crystal state. In apo state in DOPG, two different intermediate conformations are obtained. One is similar to that observed with bound adenosine in DOPC, while the other is closer to the active state but not yet fully active. Exclusively, in DOPG with bound adenosine or NECA, we reproducibly identify receptor conformations with fully active features, which are able to dock Gs protein. These different receptor conformations can be attributed to the action/absence of agonist and phospholipid-mediated allosteric effects on the intracellular side of the receptor.

Exploring the Activation Mechanism of the mGlu5 Transmembrane Domain

As a class C GPCR and regulator of synaptic activity, mGlu5 is an attractive drug target, potentially offering treatment for several neurologic and psychiatric disorders. As little is known about the activation mechanism of mGlu5 at a structural level, potential of mean force calculations linked to molecular dynamics simulations were performed on the mGlu5 transmembrane domain crystal structure to explore various internal mechanisms responsible for its activation. Our results suggest that the hydrophilic interactions between intracellular loop 1 and the intracellular side of TM6 have to be disrupted to reach a theoretically active-like conformation. In addition, interactions between residues that are key for mGlu5 activation (Tyr659^3.44 and Ile751^5.51) and mGlu5 inactivation (Tyr659^3.44 and Ser809^7.39) have been identified. Inasmuch as mGlu5 receptor signaling is poorly understood, potentially showing a complex network of micro-switches and subtle structure-activity relationships, the present study represents a step forward in the understanding of mGlu5 transmembrane domain activation.

Common activation mechanism of class A GPCRs

Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.

Toward Understanding the Impact of Dimerization Interfaces in Angiotensin II Type 1 Receptor

Angiotensin II type 1 receptor (AT1R) is a prototypical class A G protein-coupled receptor (GPCR) that has an important role in cardiovascular pathologies and blood pressure regulation as well as in the central nervous system. GPCRs may exist and function as monomers; however, they can assemble to form higher order structures, and as a result of oligomerization, their function and signaling profiles can be altered. In the case of AT1R, the classical Gα_q/11 pathway is initiated with endogenous agonist angiotensin II binding. A variety of cardiovascular pathologies such as heart failure, diabetic nephropathy, atherosclerosis, and hypertension are associated with this pathway. Recent findings reveal that AT1R can form homodimers and activate the noncanonical (β-arrestin-mediated) pathway. Nevertheless, the exact dimerization interface and atomic details of AT1R homodimerization have not been still elucidated. Here, six different symmetrical dimer interfaces of AT1R are considered, and homodimers were constructed using other published GPCR crystal dimer interfaces as template structures. These AT1R homodimers were then inserted into the model membrane bilayers and subjected to all-atom molecular dynamics simulations. Our simulation results along with the principal component analysis and water pathway analysis suggest four different interfaces as the most plausible: symmetrical transmembrane (TM)1,2,8; TM5; TM4; and TM4,5 AT1R dimer interfaces that consist of one inactive and one active protomer. Moreover, we identified ILE238^6.33 as a hub residue in the stabilization of the inactive state of AT1R.

Conformational Complexity and Dynamics in a Muscarinic Receptor Revealed by NMR Spectroscopy

The M2 muscarinic acetylcholine receptor (M2R) is a prototypical GPCR that plays important roles in regulating heart rate and CNS functions. Crystal structures provide snapshots of the M2R in inactive and active states, but the allosteric link between the ligand binding pocket and cytoplasmic surface remains poorly understood. Here we used solution NMR to examine the structure and dynamics of the M2R labeled with ¹³CH₃-ε-methionine upon binding to various orthosteric and allosteric ligands having a range of efficacy for both G protein activation and arrestin recruitment. We observed ligand-specific changes in the NMR spectra of ¹³CH₃-ε-methionine probes in the M2R extracellular domain, transmembrane core, and cytoplasmic surface, allowing us to correlate ligand structure with changes in receptor structure and dynamics. We show that the M2R has a complex energy landscape in which ligands with different efficacy profiles stabilize distinct receptor conformations.

Structure-based discovery of selective positive allosteric modulators of antagonists for the M2 muscarinic acetylcholine receptor

The orthosteric binding sites of the five muscarinic acetylcholine receptor (mAChR) subtypes are highly conserved, making the development of selective antagonists challenging. The allosteric sites of these receptors are more variable, allowing one to imagine allosteric modulators that confer subtype selectivity, which would reduce the major off-target effects of muscarinic antagonists. Accordingly, a large library docking campaign was prosecuted seeking unique positive allosteric modulators (PAMs) for antagonists, ultimately revealing a PAM that substantially potentiates antagonist binding leading to subtype selectivity at the M₂ mAChR. This study supports the feasibility of discovering PAMs that can convert an armamentarium of potent but nonselective G-protein–coupled receptor (GPCR) antagonist drugs into subtype-selective reagents.

Subtype-selective antagonists for muscarinic acetylcholine receptors (mAChRs) have long been elusive, owing to the highly conserved orthosteric binding site. However, allosteric sites of these receptors are less conserved, motivating the search for allosteric ligands that modulate agonists or antagonists to confer subtype selectivity. Accordingly, a 4.6 million-molecule library was docked against the structure of the prototypical M₂ mAChR, seeking molecules that specifically stabilized antagonist binding. This led us to identify a positive allosteric modulator (PAM) that potentiated the antagonist N-methyl scopolamine (NMS). Structure-based optimization led to compound ’628, which enhanced binding of NMS, and the drug scopolamine itself, with a cooperativity factor (α) of 5.5 and a K_B of 1.1 μM, while sparing the endogenous agonist acetylcholine. NMR spectral changes determined for methionine residues reflected changes in the allosteric network. Moreover, ’628 slowed the dissociation rate of NMS from the M₂ mAChR by 50-fold, an effect not observed at the other four mAChR subtypes. The specific PAM effect of ’628 on NMS antagonism was conserved in functional assays, including agonist stimulation of [³⁵S]GTPγS binding and ERK 1/2 phosphorylation. Importantly, the selective allostery between ’628 and NMS was retained in membranes from adult rat hypothalamus and in neonatal rat cardiomyocytes, supporting the physiological relevance of this PAM/antagonist approach. This study supports the feasibility of discovering PAMs that confer subtype selectivity to antagonists; molecules like ’628 can convert an armamentarium of potent but nonselective GPCR antagonist drugs into subtype-selective reagents, thus reducing their off-target effects.

Analysis of positive and negative allosteric modulation in metabotropic glutamate receptors 4 and 5 with a dual ligand

As class C GPCRs and regulators of synaptic activity, human metabotropic glutamate receptors (mGluRs) 4 and 5 are prime targets for allosteric modulation, with mGlu5 inhibition or mGlu4 stimulation potentially treating conditions like chronic pain and Parkinson’s disease. As an allosteric modulator that can bind both receptors, 2-Methyl-6-(phenylethynyl)pyridine (MPEP) is able to negatively modulate mGlu5 or positively modulate mGlu4. At a structural level, how it elicits these responses and how mGluRs undergo activation is unclear. Here, we employ homology modelling and 30 µs of atomistic molecular dynamics (MD) simulations to probe allosteric conformational change in mGlu4 and mGlu5, with and without docked MPEP. Our results identify several structural differences between mGlu4 and mGlu5, as well as key differences responsible for MPEP-mediated positive and negative allosteric modulation, respectively. A novel mechanism of mGlu4 activation is revealed, which may apply to all mGluRs in general. This involves conformational changes in TM3, TM4 and TM5, separation of intracellular loop 2 (ICL2) from ICL1/ICL3, and destabilization of the ionic-lock. On the other hand, mGlu5 experiences little disturbance when MPEP binds, maintaining its inactive state with reduced conformational fluctuation. In addition, when MPEP is absent, a lipid molecule can enter the mGlu5 allosteric pocket.

Diverse activation pathways in class A GPCRs converge near the G-protein-coupling region

Class A G-protein-coupled receptors (GPCRs) are a large family of membrane proteins that mediate a wide variety of physiological functions, including vision, neurotransmission and immune responses. They are the targets of nearly one-third of all prescribed medicinal drugs such as beta blockers and antipsychotics. GPCR activation is facilitated by extracellular ligands and leads to the recruitment of intracellular G proteins. Structural rearrangements of residue contacts in the transmembrane domain serve as 'activation pathways' that connect the ligand-binding pocket to the G-protein-coupling region within the receptor. In order to investigate the similarities in activation pathways across class A GPCRs, we analysed 27 GPCRs from diverse subgroups for which structures of active, inactive or both states were available. Here we show that, despite the diversity in activation pathways between receptors, the pathways converge near the G-protein-coupling region. This convergence is mediated by a highly conserved structural rearrangement of residue contacts between transmembrane helices 3, 6 and 7 that releases G-protein-contacting residues. The convergence of activation pathways may explain how the activation steps initiated by diverse ligands enable GPCRs to bind a common repertoire of G proteins.