Stepwise activation of a metabotropic glutamate receptor

Metabotropic glutamate receptors belong to a family of G protein-coupled receptors that are obligate dimers and possess a large extracellular ligand-binding domain that is linked via a cysteine-rich domain to their 7-transmembrane domain1. Upon activation, these receptors undergo a large conformational change to transmit the ligand binding signal from the extracellular ligand-binding domain to the G protein-coupling 7-transmembrane domain2. In this manuscript, we propose a model for a sequential, multistep activation mechanism of metabotropic glutamate receptor subtype 5. We present a series of structures in lipid nanodiscs, from inactive to fully active, including agonist-bound intermediate states. Further, using bulk and single-molecule fluorescence imaging, we reveal distinct receptor conformations upon allosteric modulator and G protein binding.

Domain coupling in activation of a family C GPCR

The G protein-coupled metabotropic glutamate receptors form homodimers and heterodimers with highly diverse responses to glutamate and varying physiological function. The molecular basis for this diversity remains poorly delineated. We employ molecular dynamics, single-molecule spectroscopy, and hydrogen-deuterium exchange to dissect the pathway of activation triggered by glutamate. We find that activation entails multiple loosely coupled steps and identify a novel pre-active intermediate whose transition to the active state forms dimer interactions that set signaling efficacy. Such subunit interactions generate functional diversity that differs across homodimers and heterodimers. The agonist-bound receptor is remarkably dynamic, with low occupancy of G protein-coupling conformations, providing considerable headroom for modulation of the landscape by allosteric ligands. Sites of sequence diversity within the dimerization interface and diverse coupling between activation rearrangements may contribute to precise decoding of glutamate signals and transients over broad spatial and temporal scales.

Step-wise activation of a Family C GPCR

Metabotropic glutamate receptors belong to a family of G protein-coupled receptors that are obligate dimers and possess a large extracellular ligand-binding domain (ECD) that is linked via a cysteine-rich domain (CRDs) to their 7-transmembrane (TM) domain. Upon activation, these receptors undergo a large conformational change to transmit the ligand binding signal from the ECD to the G protein-coupling TM. In this manuscript, we propose a model for a sequential, multistep activation mechanism of metabotropic glutamate receptor subtype 5. We present a series of structures in lipid nanodiscs, from inactive to fully active, including agonist-bound intermediate states. Further, using bulk and single-molecule fluorescence imaging we reveal distinct receptor conformations upon allosteric modulator and G protein binding.

A non-canonical mechanism of GPCR activation

The goal of designing safer, more effective drugs has led to tremendous interest in molecular mechanisms through which ligands can precisely manipulate signaling of G-protein-coupled receptors (GPCRs), the largest class of drug targets. Decades of research have led to the widely accepted view that all agonists-ligands that trigger GPCR activation-function by causing rearrangement of the GPCR's transmembrane helices, opening an intracellular pocket for binding of transducer proteins. Here we demonstrate that certain agonists instead trigger activation of free fatty acid receptor 1 by directly rearranging an intracellular loop that interacts with transducers. We validate the predictions of our atomic-level simulations by targeted mutagenesis; specific mutations which disrupt interactions with the intracellular loop convert these agonists into inverse agonists. Further analysis suggests that allosteric ligands could regulate signaling of many other GPCRs via a similar mechanism, offering rich possibilities for precise control of pharmaceutically important targets.

Negative allosteric modulation of the glucagon receptor by RAMP2

Receptor activity-modifying proteins (RAMPs) modulate the activity of many Family B GPCRs. We show that RAMP2 directly interacts with the glucagon receptor (GCGR), a Family B GPCR responsible for blood sugar homeostasis, and broadly inhibits receptor-induced downstream signaling. HDX-MS experiments demonstrate that RAMP2 enhances local flexibility in select locations in and near the receptor extracellular domain (ECD) and in the 6th transmembrane helix, whereas smFRET experiments show that this ECD disorder results in the inhibition of active and intermediate states of the intracellular surface. We determined the cryo-EM structure of the GCGR-Gs complex at 2.9 Å resolution in the presence of RAMP2. RAMP2 apparently does not interact with GCGR in an ordered manner; however, the receptor ECD is indeed largely disordered along with rearrangements of several intracellular hallmarks of activation. Our studies suggest that RAMP2 acts as a negative allosteric modulator of GCGR by enhancing conformational sampling of the ECD.

A saturation-mutagenesis analysis of the interplay between stability and activation in Ras

Cancer mutations in Ras occur predominantly at three hotspots: Gly 12, Gly 13, and Gln 61. Previously, we reported that deep mutagenesis of H-Ras using a bacterial assay identified many other activating mutations (Bandaru et al., 2017). We now show that the results of saturation mutagenesis of H-Ras in mammalian Ba/F3 cells correlate well with the results of bacterial experiments in which H-Ras or K-Ras are co-expressed with a GTPase-activating protein (GAP). The prominent cancer hotspots are not dominant in the Ba/F3 data. We used the bacterial system to mutagenize Ras constructs of different stabilities and discovered a feature that distinguishes the cancer hotspots. While mutations at the cancer hotspots activate Ras regardless of construct stability, mutations at lower-frequency sites (e.g. at Val 14 or Asp 119) can be activating or deleterious, depending on the stability of the Ras construct. We characterized the dynamics of three non-hotspot activating Ras mutants by using NMR to monitor hydrogen-deuterium exchange (HDX). These mutations result in global increases in HDX rates, consistent with destabilization of Ras. An explanation for these observations is that mutations that destabilize Ras increase nucleotide dissociation rates, enabling activation by spontaneous nucleotide exchange. A further stability decrease can lead to insufficient levels of folded Ras - and subsequent loss of function. In contrast, the cancer hotspot mutations are mechanism-based activators of Ras that interfere directly with the action of GAPs. Our results demonstrate the importance of GAP surveillance and protein stability in determining the sensitivity of Ras to mutational activation.

Corrole–protein interactions in H-NOX and HasA

Mutagenesis was utilised to reveal corrole–protein interactions in H-NOX and HasA. The key interaction is a hydrogen bond between the PO unit of the corrole and a protonated histidine residue.

How GPCR Phosphorylation Patterns Orchestrate Arrestin-Mediated Signaling

Binding of arrestin to phosphorylated G-protein-coupled receptors (GPCRs) controls many aspects of cell signaling. The number and arrangement of phosphates may vary substantially for a given GPCR, and different phosphorylation patterns trigger different arrestin-mediated effects. Here, we determine how GPCR phosphorylation influences arrestin behavior by using atomic-level simulations and site-directed spectroscopy to reveal the effects of phosphorylation patterns on arrestin binding and conformation. We find that patterns favoring binding differ from those favoring activation-associated conformational change. Both binding and conformation depend more on arrangement of phosphates than on their total number, with phosphorylation at different positions sometimes exerting opposite effects. Phosphorylation patterns selectively favor a wide variety of arrestin conformations, differently affecting arrestin sites implicated in scaffolding distinct signaling proteins. We also reveal molecular mechanisms of these phenomena. Our work reveals the structural basis for the long-standing "barcode" hypothesis and has important implications for design of functionally selective GPCR-targeted drugs.

Structural and functional characterization of G protein-coupled receptors with deep mutational scanning

The >800 human G protein–coupled receptors (GPCRs) are responsible for transducing diverse chemical stimuli to alter cell state- and are the largest class of drug targets. Their myriad structural conformations and various modes of signaling make it challenging to understand their structure and function. Here, we developed a platform to characterize large libraries of GPCR variants in human cell lines with a barcoded transcriptional reporter of G protein signal transduction. We tested 7800 of 7828 possible single amino acid substitutions to the beta-2 adrenergic receptor (β₂AR) at four concentrations of the agonist isoproterenol. We identified residues specifically important for β₂AR signaling, mutations in the human population that are potentially loss of function, and residues that modulate basal activity. Using unsupervised learning, we identify residues critical for signaling, including all major structural motifs and molecular interfaces. We also find a previously uncharacterized structural latch spanning the first two extracellular loops that is highly conserved across Class A GPCRs and is conformationally rigid in both the inactive and active states of the receptor. More broadly, by linking deep mutational scanning with engineered transcriptional reporters, we establish a generalizable method for exploring pharmacogenomics, structure and function across broad classes of drug receptors.

Molecular mechanism of biased signaling in a prototypical G protein-coupled receptor

Biased signaling, in which different ligands that bind to the same G protein-coupled receptor preferentially trigger distinct signaling pathways, holds great promise for the design of safer and more effective drugs. Its structural mechanism remains unclear, however, hampering efforts to design drugs with desired signaling profiles. Here, we use extensive atomic-level molecular dynamics simulations to determine how arrestin bias and G protein bias arise at the angiotensin II type 1 receptor. The receptor adopts two major signaling conformations, one of which couples almost exclusively to arrestin, whereas the other also couples effectively to a G protein. A long-range allosteric network allows ligands in the extracellular binding pocket to favor either of the two intracellular conformations. Guided by this computationally determined mechanism, we designed ligands with desired signaling profiles.

Angiotensin and biased analogs induce structurally distinct active conformations within a GPCR

Biased agonists of G protein-coupled receptors (GPCRs) preferentially activate a subset of downstream signaling pathways. In this work, we present crystal structures of angiotensin II type 1 receptor (AT1R) (2.7 to 2.9 angstroms) bound to three ligands with divergent bias profiles: the balanced endogenous agonist angiotensin II (AngII) and two strongly β-arrestin-biased analogs. Compared with other ligands, AngII promotes more-substantial rearrangements not only at the bottom of the ligand-binding pocket but also in a key polar network in the receptor core, which forms a sodium-binding site in most GPCRs. Divergences from the family consensus in this region, which appears to act as a biased signaling switch, may predispose the AT1R and certain other GPCRs (such as chemokine receptors) to adopt conformations that are capable of activating β-arrestin but not heterotrimeric G_q protein signaling.

Structure of the M2 muscarinic receptor-β-arrestin complex in a lipid nanodisc

Following agonist activation, G protein-coupled receptors (GPCRs) recruit β-arrestin, which desensitizes heterotrimeric G protein signaling and promotes receptor endocytosis¹. Additionally, β-arrestin directly regulates many cell signaling pathways that can induce cellular responses distinct from that of G proteins². Here we present a cryo-electron microscopy (cryoEM) structure of β-arrestin1 (βarr1) in complex with muscarinic acetylcholine-2-receptor (M2R) reconstituted in lipid nanodiscs. The M2R-βarr1 structure shows a multimodal network of flexible interactions, including binding of the βarr1 N-domain to phosphorylated receptor residues and βarr1 finger loop insertion into the M2R seven-transmembrane bundle, which adopts a conformation similar to that in the M2R-heterotrimeric G_o protein structure³. Moreover, the cryoEM map reveals that the βarr1 C-domain edge engages the lipid bilayer. Through atomistic simulations, biophysical, biochemical, and cellular assays, we show that the C-edge is critical for stable complex formation, βarr1 recruitment, receptor internalization, and desensitization of G protein activation. Taken together, these data suggest the cooperative interactions of β-arrestin with both the receptor and phospholipid bilayer contribute to its functional versatility.

Conformational transitions of a neurotensin receptor 1-Gi1 complex

Neurotensin receptor 1 (NTSR1) is a G-protein-coupled receptor (GPCR) that engages multiple subtypes of G protein, and is involved in the regulation of blood pressure, body temperature, weight and the response to pain. Here we present structures of human NTSR1 in complex with the agonist JMV449 and the heterotrimeric Gi1 protein, at a resolution of 3 Å. We identify two conformations: a canonical-state complex that is similar to recently reported GPCR-Gi/o complexes (in which the nucleotide-binding pocket adopts more flexible conformations that may facilitate nucleotide exchange), and a non-canonical state in which the G protein is rotated by about 45 degrees relative to the receptor and exhibits a more rigid nucleotide-binding pocket. In the non-canonical state, NTSR1 exhibits features of both active and inactive conformations, which suggests that the structure may represent an intermediate form along the activation pathway of G proteins. This structural information, complemented by molecular dynamics simulations and functional studies, provides insights into the complex process of G-protein activation.

Quantitative mapping of protein-peptide affinity landscapes using spectrally encoded beads

Transient, regulated binding of globular protein domains to Short Linear Motifs (SLiMs) in disordered regions of other proteins drives cellular signaling. Mapping the energy landscapes of these interactions is essential for deciphering and perturbing signaling networks but is challenging due to their weak affinities. We present a powerful technology (MRBLE-pep) that simultaneously quantifies protein binding to a library of peptides directly synthesized on beads containing unique spectral codes. Using MRBLE-pep, we systematically probe binding of calcineurin (CN), a conserved protein phosphatase essential for the immune response and target of immunosuppressants, to the PxIxIT SLiM. We discover that flanking residues and post-translational modifications critically contribute to PxIxIT-CN affinity and identify CN-binding peptides based on multiple scaffolds with a wide range of affinities. The quantitative biophysical data provided by this approach will improve computational modeling efforts, elucidate a broad range of weak protein-SLiM interactions, and revolutionize our understanding of signaling networks.

Smoothened stimulation by membrane sterols drives Hedgehog pathway activity

Hedgehog signalling is fundamental to embryonic development and postnatal tissue regeneration1. Aberrant postnatal Hedgehog signalling leads to several malignancies, including basal cell carcinoma and paediatric medulloblastoma2. Hedgehog proteins bind to and inhibit the transmembrane cholesterol transporter Patched-1 (PTCH1), which permits activation of the seven-transmembrane transducer Smoothened (SMO) via a mechanism that is poorly understood. Here we report the crystal structure of active mouse SMO bound to both the agonist SAG21k and to an intracellular binding nanobody that stabilizes a physiologically relevant active state. Analogous to other G protein-coupled receptors, the activation of SMO is associated with subtle motions in the extracellular domain, and larger intracellular changes. In contrast to recent models3-5, a cholesterol molecule that is critical for SMO activation is bound deep within the seven-transmembrane pocket. We propose that the inactivation of PTCH1 by Hedgehog allows a transmembrane sterol to access this seven-transmembrane site (potentially through a hydrophobic tunnel), which drives the activation of SMO. These results-combined with signalling studies and molecular dynamics simulations-delineate the structural basis for PTCH1-SMO regulation, and suggest a strategy for overcoming clinical resistance to SMO inhibitors.

Structure of a Signaling Cannabinoid Receptor 1-G Protein Complex

Cannabis elicits its mood-enhancing and analgesic effects through the cannabinoid receptor 1 (CB1), a G protein-coupled receptor (GPCR) that signals primarily through the adenylyl cyclase-inhibiting heterotrimeric G protein G_i. Activation of CB1-G_i signaling pathways holds potential for treating a number of neurological disorders and is thus crucial to understand the mechanism of G_i activation by CB1. Here, we present the structure of the CB1-G_i signaling complex bound to the highly potent agonist MDMB-Fubinaca (FUB), a recently emerged illicit synthetic cannabinoid infused in street drugs that have been associated with numerous overdoses and fatalities. The structure illustrates how FUB stabilizes the receptor in an active state to facilitate nucleotide exchange in G_i. The results compose the structural framework to explain CB1 activation by different classes of ligands and provide insights into the G protein coupling and selectivity mechanisms adopted by the receptor.

Structure of the µ-opioid receptor-Gi protein complex

The μ-opioid receptor (μOR) is a G-protein-coupled receptor (GPCR) and the target of most clinically and recreationally used opioids. The induced positive effects of analgesia and euphoria are mediated by μOR signalling through the adenylyl cyclase-inhibiting heterotrimeric G protein Gi. Here we present the 3.5 Å resolution cryo-electron microscopy structure of the μOR bound to the agonist peptide DAMGO and nucleotide-free Gi. DAMGO occupies the morphinan ligand pocket, with its N terminus interacting with conserved receptor residues and its C terminus engaging regions important for opioid-ligand selectivity. Comparison of the μOR-Gi complex to previously determined structures of other GPCRs bound to the stimulatory G protein Gs reveals differences in the position of transmembrane receptor helix 6 and in the interactions between the G protein α-subunit and the receptor core. Together, these results shed light on the structural features that contribute to the Gi protein-coupling specificity of the µOR.

Catalytic activation of β-arrestin by GPCRs

β-arrestins are critical regulator and transducer proteins for G-protein-coupled receptors (GPCRs). β-arrestin is widely believed to be activated by forming a stable and stoichiometric GPCR-β-arrestin scaffold complex, which requires and is driven by the phosphorylated tail of the GPCR. Here we demonstrate a distinct and additional mechanism of β-arrestin activation that does not require stable GPCR-β-arrestin scaffolding or the GPCR tail. Instead, it occurs through transient engagement of the GPCR core, which destabilizes a conserved inter-domain charge network in β-arrestin. This promotes capture of β-arrestin at the plasma membrane and its accumulation in clathrin-coated endocytic structures (CCSs) after dissociation from the GPCR, requiring a series of interactions with membrane phosphoinositides and CCS-lattice proteins. β-arrestin clustering in CCSs in the absence of the upstream activating GPCR is associated with a β-arrestin-dependent component of the cellular ERK (extracellular signal-regulated kinase) response. These results delineate a discrete mechanism of cellular β-arrestin function that is activated catalytically by GPCRs.

Mechanism of intracellular allosteric β2AR antagonist revealed by X-ray crystal structure

G-protein-coupled receptors (GPCRs) pose challenges for drug discovery efforts because of the high degree of structural homology in the orthosteric pocket, particularly for GPCRs within a single subfamily, such as the nine adrenergic receptors. Allosteric ligands may bind to less-conserved regions of these receptors and therefore are more likely to be selective. Unlike orthosteric ligands, which tonically activate or inhibit signalling, allosteric ligands modulate physiologic responses to hormones and neurotransmitters, and may therefore have fewer adverse effects. The majority of GPCR crystal structures published to date were obtained with receptors bound to orthosteric antagonists, and only a few structures bound to allosteric ligands have been reported. Compound 15 (Cmpd-15) is an allosteric modulator of the β₂ adrenergic receptor (β₂AR) that was recently isolated from a DNA-encoded small-molecule library. Orthosteric β-adrenergic receptor antagonists, known as beta-blockers, are amongst the most prescribed drugs in the world and Cmpd-15 is the first allosteric beta-blocker. Cmpd-15 exhibits negative cooperativity with agonists and positive cooperativity with inverse agonists. Here we present the structure of the β₂AR bound to a polyethylene glycol-carboxylic acid derivative (Cmpd-15PA) of this modulator. Cmpd-15PA binds to a pocket formed primarily by the cytoplasmic ends of transmembrane segments 1, 2, 6 and 7 as well as intracellular loop 1 and helix 8. A comparison of this structure with inactive- and active-state structures of the β₂AR reveals the mechanism by which Cmpd-15 modulates agonist binding affinity and signalling.

Identification of Phosphorylation Codes for Arrestin Recruitment by G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) mediate diverse signaling in part through interaction with arrestins, whose binding promotes receptor internalization and signaling through G protein-independent pathways. High-affinity arrestin binding requires receptor phosphorylation, often at the receptor's C-terminal tail. Here, we report an X-ray free electron laser (XFEL) crystal structure of the rhodopsin-arrestin complex, in which the phosphorylated C terminus of rhodopsin forms an extended intermolecular β sheet with the N-terminal β strands of arrestin. Phosphorylation was detected at rhodopsin C-terminal tail residues T336 and S338. These two phospho-residues, together with E341, form an extensive network of electrostatic interactions with three positively charged pockets in arrestin in a mode that resembles binding of the phosphorylated vasopressin-2 receptor tail to β-arrestin-1. Based on these observations, we derived and validated a set of phosphorylation codes that serve as a common mechanism for phosphorylation-dependent recruitment of arrestins by GPCRs.

Continuum approaches to understanding ion and peptide interactions with the membrane

Experimental and computational studies have shown that cellular membranes deform to stabilize the inclusion of transmembrane (TM) proteins harboring charge. Recent analysis suggests that membrane bending helps to expose charged and polar residues to the aqueous environment and polar head groups. We previously used elasticity theory to identify membrane distortions that minimize the insertion of charged TM peptides into the membrane. Here, we extend our work by showing that it also provides a novel, computationally efficient method for exploring the energetics of ion and small peptide penetration into membranes. First, we show that the continuum method accurately reproduces energy profiles and membrane shapes generated from molecular simulations of bare ion permeation at a fraction of the computational cost. Next, we demonstrate that the dependence of the ion insertion energy on the membrane thickness arises primarily from the elastic properties of the membrane. Moreover, the continuum model readily provides a free energy decomposition into components not easily determined from molecular dynamics. Finally, we show that the energetics of membrane deformation strongly depend on membrane patch size both for ions and peptides. This dependence is particularly strong for peptides based on simulations of a known amphipathic, membrane binding peptide from the human pathogen Toxoplasma gondii. In total, we address shortcomings and advantages that arise from using a variety of computational methods in distinct biological contexts.

Membrane bending is critical for the stability of voltage sensor segments in the membrane

The interaction between membrane proteins and the surrounding membrane is becoming increasingly appreciated for its role in regulating protein function, protein localization, and membrane morphology. In particular, recent studies have suggested that membrane deformation is needed to stably accommodate proteins harboring charged amino acids in their transmembrane (TM) region, as it is energetically prohibitive to bury charge in the hydrophobic core of the bilayer. Unfortunately, current computational methods are poorly equipped for describing such deformations, as atomistic simulations are often too short to observe large-scale membrane reorganization and most continuum approaches assume a flat membrane. Previously, we developed a method that overcomes these shortcomings by using elasticity theory to characterize equilibrium membrane distortions in the presence of a TM protein, while using traditional continuum electrostatic and nonpolar energy models to determine the energy of the protein in the membrane. Here, we linked the elastostatics, electrostatics, and nonpolar numeric solvers to permit the calculation of energies for nontrivial membrane deformations. We then coupled this procedure to a robust search algorithm that identifies optimal membrane shapes for a TM protein of arbitrary chemical composition. This advance now permits us to explore a host of biological phenomena that were beyond the scope of our original method. We show that the energy required to embed charged residues in the membrane can be highly nonadditive, and our model provides a simple mechanical explanation for this nonadditivity. Our results also predict that isolated voltage sensor segments do not insert into rigid membranes, but membrane bending dramatically stabilizes these proteins in the bilayer despite their high charge content. Additionally, we use the model to explore hydrophobic mismatch with regard to nonpolar peptides and mechanosensitive channels. Our method is in quantitative agreement with molecular dynamics simulations at a tiny fraction of the computational cost.