Position of Transmembrane Helix 6 Determines Receptor G Protein Coupling Specificity

PepBinding: A Workflow for Predicting Peptide Binding Structures by Combining Peptide Docking and Peptide Gaussian Accelerated Molecular Dynamics Simulations

Predicting protein-peptide interactions is crucial for understanding peptide binding processes and designing peptide drugs. However, traditional computational modeling approaches face challenges in accurately predicting peptide-protein binding structures due to the slow dynamics and high flexibility of the peptides. Here, we introduce a new workflow termed "PepBinding" for predicting peptide binding structures, which combines peptide docking, all-atom enhanced sampling simulations using the Peptide Gaussian accelerated Molecular Dynamics (Pep-GaMD) method, and structural clustering. PepBinding has been demonstrated on seven distinct model peptides. In peptide docking using HPEPDOCK, the peptide backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures ranged from 3.8 to 16.0 Å, corresponding to the medium to inaccurate quality models according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. The Pep-GaMD simulations performed for only 200 ns significantly improved the docking models, resulting in five medium and two acceptable quality models. Therefore, PepBinding is an efficient workflow for predicting peptide binding structures and is publicly available at https://github.com/MiaoLab20/PepBinding.

Mechanistic insights into G-protein coupling with an agonist-bound G-protein-coupled receptor

G-protein-coupled receptors (GPCRs) activate heterotrimeric G proteins by promoting guanine nucleotide exchange. Here, we investigate the coupling of G proteins with GPCRs and describe the events that ultimately lead to the ejection of GDP from its binding pocket in the Gα subunit, the rate-limiting step during G-protein activation. Using molecular dynamics simulations, we investigate the temporal progression of structural rearrangements of GDP-bound Gs protein (Gs·GDP; hereafter GsGDP) upon coupling to the β2-adrenergic receptor (β2AR) in atomic detail. The binding of GsGDP to the β2AR is followed by long-range allosteric effects that significantly reduce the energy needed for GDP release: the opening of α1-αF helices, the displacement of the αG helix and the opening of the α-helical domain. Signal propagation to the Gs occurs through an extended receptor interface, including a lysine-rich motif at the intracellular end of a kinked transmembrane helix 6, which was confirmed by site-directed mutagenesis and functional assays. From this β2AR-GsGDP intermediate, Gs undergoes an in-plane rotation along the receptor axis to approach the β2AR-Gsempty state. The simulations shed light on how the structural elements at the receptor-G-protein interface may interact to transmit the signal over 30 Å to the nucleotide-binding site. Our analysis extends the current limited view of nucleotide-free snapshots to include additional states and structural features responsible for signaling and G-protein coupling specificity.

Determinants of Neutral Antagonism and Inverse Agonism in the β2-Adrenergic Receptor

Free-energy profiles for the activation/deactivation of the β2-adrenergic receptor (ADRB2) with neutral antagonist and inverse agonist ligands have been determined with well-tempered multiple-walker (MW) metadynamics simulations. The inverse agonists carazolol and ICI118551 clearly favor single inactive conformational minima in both the binary and ternary ligand-receptor-G-protein complexes, in accord with the inverse-agonist activity of the ligands. The behavior of neutral antagonists is more complex, as they seem also to affect the recruitment of the G-protein. The results are analyzed in terms of the conformational states of the well-known microswitches that have been proposed as indicators of receptor activity.

Molecular mechanism of muscarinic acetylcholine receptor M3 interaction with Gq

Muscarinic acetylcholine receptor M3 (M3) and its downstream effector Gq/11 are critical drug development targets due to their involvement in physiopathological processes. Although the structure of the M3-miniGq complex was recently published, the lack of information on the intracellular loop 3 (ICL3) of M3 and extensive modification of Gαq impedes the elucidation of the molecular mechanism of M3-Gq coupling under more physiological condition. Here, we describe the molecular mechanism underlying the dynamic interactions between full-length wild-type M3 and Gq using hydrogen-deuterium exchange mass spectrometry and NanoLuc Binary Technology-based cell systems. We propose a detailed analysis of M3-Gq coupling through examination of previously well-defined binding interfaces and neglected regions. Our findings suggest potential binding interfaces between M3 and Gq in pre-assembled and functionally active complexes. Furthermore, M3 ICL3 negatively affected M3-Gq coupling, and the Gαq AHD underwent unique conformational changes during M3-Gq coupling.

β-Arrestin-dependent and -independent endosomal G protein activation by the vasopressin type 2 receptor

The vasopressin type 2 receptor (V2R) is an essential G protein-coupled receptor (GPCR) in renal regulation of water homeostasis. Upon stimulation, the V2R activates Gαs and Gαq/11, which is followed by robust recruitment of β-arrestins and receptor internalization into endosomes. Unlike canonical GPCR signaling, the β-arrestin association with the V2R does not terminate Gαs activation, and thus, Gαs-mediated signaling is sustained while the receptor is internalized. Here, we demonstrate that this V2R ability to co-interact with G protein/β-arrestin and promote endosomal G protein signaling is not restricted to Gαs, but also involves Gαq/11. Furthermore, our data imply that β-arrestins potentiate Gαs/Gαq/11 activation at endosomes rather than terminating their signaling. Surprisingly, we found that the V2R internalizes and promote endosomal G protein activation independent of β-arrestins to a minor degree. These new observations challenge the current model of endosomal GPCR signaling and suggest that this event can occur in both β-arrestin-dependent and -independent manners.

β-arrestin-dependent and -independent endosomal G protein activation by the vasopressin type 2 receptor

The vasopressin type 2 receptor (V2R) is an essential GPCR in renal regulation of water homeostasis. Upon stimulation, the V2R activates Gαs and Gαq/11, which is followed by robust recruitment of β-arrestins and receptor internalization into endosomes. Unlike canonical GPCR signaling, the β-arrestin association with the V2R does not terminate Gαs activation, and thus, Gαs-mediated signaling is sustained while the receptor is internalized. Here, we demonstrate that this V2R ability to co-interact with G protein/β-arrestin and promote endosomal G protein signaling is not restricted to Gαs, but also involves Gαq/11. Furthermore, our data implies that β-arrestins potentiate Gαs/Gαq/11 activation at endosomes rather than terminating their signaling. Surprisingly, we found that the V2R internalizes and promote endosomal G protein activation independent of β-arrestins to a minor degree. These new observations challenge the current model of endosomal GPCR signaling and suggest that this event can occur in both β-arrestin-dependent and -independent manners.

7TM domain structures of adhesion GPCRs: what's new and what's missing?

Adhesion-type G protein-coupled receptors (aGPCRs) have long resisted approaches to resolve the structural details of their heptahelical transmembrane (7TM) domains. Single-particle cryogenic electron microscopy (cryo-EM) has recently produced aGPCR 7TM domain structures for ADGRD1, ADGRG1, ADGRG2, ADGRG3, ADGRG4, ADGRG5, ADGRF1, and ADGRL3. We review the unique properties, including the position and conformation of their activating tethered agonist (TA) and signaling motifs within the 7TM bundle, that the novel structures have helped to identify. We also discuss questions that the kaleidoscope of novel aGPCR 7TM domain structures have left unanswered. These concern the relative positions, orientations, and interactions of the 7TM and GPCR autoproteolysis-inducing (GAIN) domains with one another. Clarifying their interplay remains an important goal of future structural studies on aGPCRs.

Molecular Biophysics of Class A G Protein Coupled Receptors-Lipids Interactome at a Glance-Highlights from the A2A Adenosine Receptor

G protein-coupled receptors (GPCRs) are embedded in phospholipid membrane bilayers with cholesterol representing 34% of the total lipid content in mammalian plasma membranes. Membrane lipids interact with GPCRs structures and modulate their function and drug-stimulated signaling through conformational selection. It has been shown that anionic phospholipids form strong interactions between positively charged residues in the G protein and the TM5-TM6-TM 7 cytoplasmic interface of class A GPCRs stabilizing the signaling GPCR-G complex. Cholesterol with a high content in plasma membranes can be identified in more specific sites in the transmembrane region of GPCRs, such as the Cholesterol Consensus Motif (CCM) and Cholesterol Recognition Amino Acid Consensus (CRAC) motifs and other receptor dependent and receptor state dependent sites. Experimental biophysical methods, atomistic (AA) MD simulations and coarse-grained (CG) molecular dynamics simulations have been applied to investigate these interactions. We emphasized here the impact of phosphatidyl inositol-4,5-bisphosphate (PtdIns(4,5)P₂ or PIP₂), a minor phospholipid component and of cholesterol on the function-related conformational equilibria of the human A_2A adenosine receptor (A_2AR), a representative receptor in class A GPCR. Several GPCRs of class A interacted with PIP₂ and cholesterol and in many cases the mechanism of the modulation of their function remains unknown. This review provides a helpful comprehensive overview for biophysics that enter the field of GPCRs-lipid systems.

The role of G protein conformation in receptor-G protein selectivity

G protein-coupled receptors (GPCRs) selectively activate at least one of the four families of heterotrimeric G proteins, but the mechanism of coupling selectivity remains unclear. Structural studies emphasize structural complementarity of GPCRs and nucleotide-free G proteins, but selectivity is likely to be determined by transient intermediate-state complexes that exist before nucleotide release. Here we study coupling to nucleotide-decoupled G protein variants that can adopt conformations similar to receptor-bound G proteins without releasing nucleotide, and are therefore able to bypass intermediate-state complexes. We find that selectivity is degraded when nucleotide release is not required for GPCR-G protein complex formation, to the extent that most GPCRs interact with most nucleotide-decoupled G proteins. These findings demonstrate the absence of absolute structural incompatibility between noncognate receptor-G protein pairs, and are consistent with the hypothesis that transient intermediate states are partly responsible for coupling selectivity.

New insights into GPCR coupling and dimerisation from cryo-EM structures

Over the past three years (2020-2022) more structures of GPCRs have been determined than in the previous twenty years (2000-2019), primarily of GPCR complexes that are large enough for structure determination by single-particle cryo-EM. This review will present some structural highlights that have advanced our molecular understanding of promiscuous G protein coupling, how a G protein receptor kinase and β-arrestins couple to GPCRs, and GPCR dimerisation. We will also discuss advances in the use of gene fusions, nanobodies, and Fab fragments to facilitate the structure determination of GPCRs in the inactive state that, on their own, are too small for structure determination by single-particle cryo-EM.

Structures of Ric-8B in complex with Gα protein folding clients reveal isoform specificity mechanisms

Mammalian Ric-8 proteins act as chaperones to regulate the cellular abundance of heterotrimeric G protein α subunits. The Ric-8A isoform chaperones Gαi/o, Gα12/13, and Gαq/11 subunits, while Ric-8B acts on Gαs/olf subunits. Here, we determined cryoelectron microscopy (cryo-EM) structures of Ric-8B in complex with Gαs and Gαolf, revealing isoform differences in the relative positioning and contacts between the C-terminal α5 helix of Gα within the concave pocket formed by Ric-8 α-helical repeat elements. Despite the overall architectural similarity with our earlier structures of Ric-8A complexed to Gαq and Gαi1, Ric-8B distinctly accommodates an extended loop found only in Gαs/olf proteins. The structures, along with results from Ric-8 protein thermal stability assays and cell-based Gαolf folding assays, support a requirement for the Gα C-terminal region for binding specificity, and highlight that multiple structural elements impart specificity for Ric-8/G protein binding.

Ligand and G-protein selectivity in the κ-opioid receptor

The κ-opioid receptor (KOR) represents a highly desirable therapeutic target for treating not only pain but also addiction and affective disorders1. However, the development of KOR analgesics has been hindered by the associated hallucinogenic side effects2. The initiation of KOR signalling requires the Gi/o-family proteins including the conventional (Gi1, Gi2, Gi3, GoA and GoB) and nonconventional (Gz and Gg) subtypes. How hallucinogens exert their actions through KOR and how KOR determines G-protein subtype selectivity are not well understood. Here we determined the active-state structures of KOR in a complex with multiple G-protein heterotrimers-Gi1, GoA, Gz and Gg-using cryo-electron microscopy. The KOR-G-protein complexes are bound to hallucinogenic salvinorins or highly selective KOR agonists. Comparisons of these structures reveal molecular determinants critical for KOR-G-protein interactions as well as key elements governing Gi/o-family subtype selectivity and KOR ligand selectivity. Furthermore, the four G-protein subtypes display an intrinsically different binding affinity and allosteric activity on agonist binding at KOR. These results provide insights into the actions of opioids and G-protein-coupling specificity at KOR and establish a foundation to examine the therapeutic potential of pathway-selective agonists of KOR.

Rhodopsin, light-sensor of vision

The light sensor of vertebrate scotopic (low-light) vision, rhodopsin, is a G-protein-coupled receptor comprising a polypeptide chain with bound chromophore, 11-cis-retinal, that exhibits remarkable physicochemical properties. This photopigment is extremely stable in the dark, yet its chromophore isomerises upon photon absorption with 70% efficiency, enabling the activation of its G-protein, transducin, with high efficiency. Rhodopsin's photochemical and biochemical activities occur over very different time-scales: the energy of retinaldehyde's excited state is stored in <1 ps in retinal-protein interactions, but it takes milliseconds for the catalytically active state to form, and many tens of minutes for the resting state to be restored. In this review, we describe the properties of rhodopsin and its role in rod phototransduction. We first introduce rhodopsin's gross structural features, its evolution, and the basic mechanisms of its activation. We then discuss light absorption and spectral sensitivity, photoreceptor electrical responses that result from the activity of individual rhodopsin molecules, and recovery of rhodopsin and the visual system from intense bleaching exposures. We then provide a detailed examination of rhodopsin's molecular structure and function, first in its dark state, and then in the active Meta states that govern its interactions with transducin, rhodopsin kinase and arrestin. While it is clear that rhodopsin's molecular properties are exquisitely honed for phototransduction, from starlight to dawn/dusk intensity levels, our understanding of how its molecular interactions determine the properties of scotopic vision remains incomplete. We describe potential future directions of research, and outline several major problems that remain to be solved.

Gαi-derived peptide binds the µ-opioid receptor

BACKGROUND

G protein-coupled receptors (GPCRs) transduce external stimuli into the cell by G proteins via an allosteric mechanism. Agonist binding to the receptor stimulates GDP/GTP exchange within the heterotrimeric G protein complex, whereas recent structures of GPCR-G protein complexes revealed that the H5, S1 and S2 domains of Gα are involved in binding the active receptor, earlier studies showed that a short peptide analog derived from the C-terminus (H5) of the G protein transducin (G_t) is sufficient to stabilize rhodopsin in an active form.

METHODS

We have used Molecular Dynamics simulations along with biological evaluation by means of radio-ligand binding assay to study the interactions between Gα_i-derived peptide (G-peptide) and the µ-opioid receptor (µOR).

RESULTS

Here, we show that a Gα_i-derived peptide of 12 amino acids binds the µ-opioid receptor and acts as an allosteric modulator. The Gα_i-derived peptide increases µOR affinity for its agonist morphine in a dose-dependent way.

CONCLUSIONS

These results indicate that the GPCR-Gα peptide interaction observed so far for only rhodopsin can be extrapolated to µOR. In addition, we show that the C-terminal peptide of the Gα_i subunit is sufficient to stabilize the active conformation of the receptor. Our approach opens the possibility to investigate the GPCR-G protein interface with peptide modification.

Functional optimization of light-activatable Opto-GPCRs: Illuminating the importance of the proximal C-terminus in G-protein specificity

Introduction: G-protein coupled receptors (GPCRs) are the largest family of human receptors that transmit signals from natural ligands and pharmaceutical drugs into essentially every physiological process. One main characteristic of G-protein coupled receptors is their ability to specifically couple with different families of G-proteins, thereby triggering specific downstream signaling pathways. While an abundance of structural information is available on G-protein coupled receptorn interactions with G-proteins, little is known about the G-protein coupled receptor domains functionally mediating G-protein specificity, in particular the proximal C-terminus, the structure which cannot be predicted with high confidentiality due to its flexibility. Methods: In this study, we exploited OptoGPCR chimeras between lightgated G-protein coupled receptors (opsins) and ligand-gated G-protein coupled receptors to systematically investigate the involvement of the C-terminus steering G-protein specificity. We employed rhodopsin-beta2-adrenoceptor and melanopsin-mGluR6 chimeras in second messenger assays and developed structural models of the chimeras. Results: We discovered a dominant role of the proximal C-terminus, dictating G-protein selectivity in the melanopsin-mGluR6 chimera, whereas it is the intracellular loop 3, which steers G-protein tropism in the rhodopsin-beta2-adrenoceptor. From the functional results and structural predictions, melanopsin and mGluR6 use a different mechanism to bovine rhodopsin and b2AR to couple to a selective G-protein. Discussion: Collectively, this work adds knowledge to the G-protein coupled receptor domains mediating G-protein selectivity, ultimately paving the way to optogenetically elicited specific G-protein signaling on demand.

Evolutionary association of receptor-wide amino acids with G protein-coupling selectivity in aminergic GPCRs

Evolutionary analyses of aminergic G protein–coupled receptors reveal receptor-wide potential determinants of G protein–coupling selectivity.

G protein-coupled receptors (GPCRs) induce signal transduction pathways through coupling to four main subtypes of G proteins (G_s, G_i, G_q, and G_12/13), selectively. However, G protein selective activation mechanisms and residual determinants in GPCRs have remained obscure. Herein, we performed extensive phylogenetic analysis and identified specifically conserved residues for the aminergic receptors having similar coupling profiles. By integrating our methodology of differential evolutionary conservation of G protein–specific amino acids with structural analyses, we identified specific activation networks for G_s, G_i1, G_o, and G_q. To validate that these networks could determine coupling selectivity we further analyzed G_s-specific activation network and its association with G_s selectivity. Through molecular dynamics simulations, we showed that previously uncharacterized Glycine at position 7x41 plays an important role in receptor activation and it may determine G_s coupling selectivity by facilitating a larger TM6 movement. Finally, we gathered our results into a comprehensive model of G protein selectivity called “sequential switches of activation” describing three main molecular switches controlling GPCR activation: ligand binding, G protein selective activation mechanisms, and G protein contact.

The Conformational Dynamics of Heterotrimeric G Proteins During GPCR-Mediated Activation

Heterotrimeric G proteins (G proteins) are essential cellular signaling proteins that mediate extracellular signals to achieve various cellular functions. G-protein-coupled receptors (GPCRs) are the major guanine nucleotide exchange factors (GEFs) that induce G proteins to release guanosine diphosphate and rapidly bind to guanosine triphosphate, resulting in G protein activation. G proteins undergo dynamic conformational changes during the activation/inactivation process, and the precise structural mechanism of GPCR-mediated G protein activation is of great interest. Over the last decade, a number of GPCR-G protein complex structures have been identified, yet an understanding of the mechanisms underlying allosteric conformational changes during receptor-mediated G protein activation and GPCR-G protein coupling selectivity is only now emerging. This review discusses recent studies on the dynamic conformational changes of G proteins and provides insight into the structural mechanism of GPCR-mediated G protein activation.

DEER Analysis of GPCR Conformational Heterogeneity

G protein-coupled receptors (GPCRs) represent a large class of transmembrane helical proteins which are involved in numerous physiological signaling pathways and therefore represent crucial pharmacological targets. GPCR function and the action of therapeutic molecules are defined by only a few parameters, including receptor basal activity, ligand affinity, intrinsic efficacy and signal bias. These parameters are encoded in characteristic receptor conformations existing in equilibrium and their populations, which are thus of paramount interest for the understanding of receptor (mal-)functions and rational design of improved therapeutics. To this end, the combination of site-directed spin labeling and EPR spectroscopy, in particular double electron-electron resonance (DEER), is exceedingly valuable as it has access to sub-Angstrom spatial resolution and provides a detailed picture of the number and populations of conformations in equilibrium. This review gives an overview of existing DEER studies on GPCRs with a focus on the delineation of structure/function frameworks, highlighting recent developments in data analysis and visualization. We introduce "conformational efficacy" as a parameter to describe ligand-specific shifts in the conformational equilibrium, taking into account the loose coupling between receptor segments observed for different GPCRs using DEER.

The conformational transition during G protein-coupled receptor (GPCR) and G protein interaction

The precise structural mechanism of G protein-coupled receptor (GPCR)-G protein coupling has been of significant research interest because it provides fundamental knowledge on cellular signaling and valuable information for GPCR-targeted drug development. Over the last decade, several GPCR-G protein complex structures have been identified. However, these structures are mere snapshots of guanosine diphosphate (GDP)-released stable GPCR-G protein complexes, which have limited the understanding of the allosteric conformational transition during receptor binding to GDP release and the GPCR-G protein coupling selectivity. Recently, deeper insights into the mechanism underlying stepwise conformational changes during GPCR-G protein coupling were obtained using hydrogen/deuterium exchange mass spectrometry, hydroxyl radical footprinting-mass spectrometry, X-ray crystallography, cryoelectron microscopy, and molecular dynamics simulation techniques. This review summarizes these recent developments.

Functionally selective and biased agonists of muscarinic receptors

Disruption of cholinergic signalling via muscarinic receptors is associated with various pathologies, like Alzheimer's disease or schizophrenia. Selective muscarinic agonists possess therapeutic potential in the treatment of diabetes, pain or Sjögren's syndrome. The orthosteric binding site of all subtypes of the muscarinic receptor is structurally identical, making the development of affinity-based selective agonists virtually impossible. Some agonists, however, are functionally selective; they activate only a subset of receptors or signalling pathways. Others may stabilise specific conformations of the receptor leading to non-uniform modulation of individual signalling pathways (biased agonists). Functionally selective and biased agonists represent a promising approach for selective activation of individual subtypes of muscarinic receptors. In this work we review chemical structures, receptor binding and agonist-specific conformations of currently known functionally selective and biased muscarinic agonists in the context of their intricate intracellular signalling. Further, we take a perspective on the possible use of biased agonists for tissue and organ-specific activation of muscarinic receptors.

Receptor-Arrestin Interactions: The GPCR Perspective

Arrestins are a small family of four proteins in most vertebrates that bind hundreds of different G protein-coupled receptors (GPCRs). Arrestin binding to a GPCR has at least three functions: precluding further receptor coupling to G proteins, facilitating receptor internalization, and initiating distinct arrestin-mediated signaling. The molecular mechanism of arrestin–GPCR interactions has been extensively studied and discussed from the “arrestin perspective”, focusing on the roles of arrestin elements in receptor binding. Here, we discuss this phenomenon from the “receptor perspective”, focusing on the receptor elements involved in arrestin binding and emphasizing existing gaps in our knowledge that need to be filled. It is vitally important to understand the role of receptor elements in arrestin activation and how the interaction of each of these elements with arrestin contributes to the latter’s transition to the high-affinity binding state. A more precise knowledge of the molecular mechanisms of arrestin activation is needed to enable the construction of arrestin mutants with desired functional characteristics.

Ligand recognition and G-protein coupling selectivity of cholecystokinin A receptor

Cholecystokinin A receptor (CCK_AR) belongs to family A G-protein-coupled receptors and regulates nutrient homeostasis upon stimulation by cholecystokinin (CCK). It is an attractive drug target for gastrointestinal and metabolic diseases. One distinguishing feature of CCK_AR is its ability to interact with a sulfated ligand and to couple with divergent G-protein subtypes, including G_s, G_i and G_q. However, the basis for G-protein coupling promiscuity and ligand recognition by CCK_AR remains unknown. Here, we present three cryo-electron microscopy structures of sulfated CCK-8-activated CCK_AR in complex with G_s, G_i and G_q heterotrimers, respectively. CCK_AR presents a similar conformation in the three structures, whereas conformational differences in the 'wavy hook' of the Gα subunits and ICL3 of the receptor serve as determinants in G-protein coupling selectivity. Our findings provide a framework for understanding G-protein coupling promiscuity by CCK_AR and uncover the mechanism of receptor recognition by sulfated CCK-8.

Understanding the Mechanism of Activation/Deactivation of GLP-1R via Accelerated Molecular Dynamics Simulation

Glucagon-like peptide-1 receptor (GLP-1R), as a member of the class B G protein-coupled receptors (GPCRs), plays a crucial role in regulating blood glucose level signal recognition through its activation. The conformation changes during the activation pathway are of particular importance for its function. To investigate the activation mechanism of GLP-1R, the crystal structures of active and inactive forms are chosen to perform a total of 2 μs of accelerated molecular dynamics (aMD) simulations and 400ns of conventional molecular dynamics (cMD) simulations. With the aid of structural analysis and potential of mean force (PMF) calculations, we reveal the role of different helices in the activation and deactivation process and obtain the intermediate states during activation and deactivation that are difficult to capture in experiments. Protein structure network (PSN) was utilised to clarify the allosteric communication pathways of activation and deactivation and reveal the mechanisms of its activation and deactivation. The results could advance our understanding of the activation mechanism of GLP-1R and the related drug design.

G protein-coupled receptor-G protein interactions: a single-molecule perspective

G protein-coupled receptors (GPCRs) regulate many cellular and physiological processes, responding to a diverse range of extracellular stimuli including hormones, neurotransmitters, odorants, and light. Decades of biochemical and pharmacological studies have provided fundamental insights into the mechanisms of GPCR signaling. Thanks to recent advances in structural biology, we now possess an atomistic understanding of receptor activation and G protein coupling. However, how GPCRs and G proteins interact in living cells to confer signaling efficiency and specificity remains insufficiently understood. The development of advanced optical methods, including single-molecule microscopy, has provided the means to study receptors and G proteins in living cells with unprecedented spatio-temporal resolution. The results of these studies reveal an unexpected level of complexity, whereby GPCRs undergo transient interactions among themselves as well as with G proteins and structural elements of the plasma membrane to form short-lived signaling nanodomains that likely confer both rapidity and specificity to GPCR signaling. These findings may provide new strategies to pharmaceutically modulate GPCR function, which might eventually pave the way to innovative drugs for common diseases such as diabetes or heart failure.

The Dynamics of the Neuropeptide Y Receptor Type 1 Investigated by Solid-State NMR and Molecular Dynamics Simulation

We report data on the structural dynamics of the neuropeptide Y (NPY) G-protein-coupled receptor (GPCR) type 1 (Y1R), a typical representative of class A peptide ligand GPCRs, using a combination of solid-state NMR and molecular dynamics (MD) simulation. First, the equilibrium dynamics of Y1R were studied using 15N-NMR and quantitative determination of 1H-13C order parameters through the measurement of dipolar couplings in separated-local-field NMR experiments. Order parameters reporting the amplitudes of the molecular motions of the C-H bond vectors of Y1R in DMPC membranes are 0.57 for the Cα sites and lower in the side chains (0.37 for the CH2 and 0.18 for the CH3 groups). Different NMR excitation schemes identify relatively rigid and also dynamic segments of the molecule. In monounsaturated membranes composed of longer lipid chains, Y1R is more rigid, attributed to a higher hydrophobic thickness of the lipid membrane. The presence of an antagonist or NPY has little influence on the amplitude of motions, whereas the addition of agonist and arrestin led to a pronounced rigidization. To investigate Y1R dynamics with site resolution, we conducted extensive all-atom MD simulations of the apo and antagonist-bound state. In each state, three replicas with a length of 20 μs (with one exception, where the trajectory length was 10 μs) were conducted. In these simulations, order parameters of each residue were determined and showed high values in the transmembrane helices, whereas the loops and termini exhibit much lower order. The extracellular helix segments undergo larger amplitude motions than their intracellular counterparts, whereas the opposite is observed for the loops, Helix 8, and termini. Only minor differences in order were observed between the apo and antagonist-bound state, whereas the time scale of the motions is shorter for the apo state. Although these relatively fast motions occurring with correlation times of ns up to a few µs have no direct relevance for receptor activation, it is believed that they represent the prerequisite for larger conformational transitions in proteins.

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.

Agonist Binding and G Protein Coupling in Histamine H2 Receptor: A Molecular Dynamics Study

The histamine H2 receptor (H₂R) plays an important role in the regulation of gastric acid secretion. Therefore, it is a main drug target for the treatment of gastroesophageal reflux or peptic ulcer disease. However, there is as of yet no 3D-structural information available hampering a mechanistic understanding of H₂R. Therefore, we created a model of the histamine-H₂R-G_s complex based on the structure of the ternary complex of the β2-adrenoceptor and investigated the conformational stability of this active GPCR conformation. Since the physiologically relevant motions with respect to ligand binding and conformational changes of GPCRs can only partly be assessed on the timescale of conventional MD (cMD) simulations, we also applied metadynamics and Gaussian accelerated molecular dynamics (GaMD) simulations. A multiple walker metadynamics simulation in combination with cMD was applied for the determination of the histamine binding mode. The preferential binding pose detected is in good agreement with previous data from site directed mutagenesis and provides a basis for rational ligand design. Inspection of the H₂R-G_s interface reveals a network of polar interactions that may contribute to H₂R coupling selectivity. The cMD and GaMD simulations demonstrate that the active conformation is retained on a μs-timescale in the ternary histamine-H₂R-G_s complex and in a truncated complex that contains only G_s helix α5 instead of the entire G protein. In contrast, histamine alone is unable to stabilize the active conformation, which is in line with previous studies of other GPCRs.

Structural mechanism underlying primary and secondary coupling between GPCRs and the Gi/o family

Heterotrimeric G proteins are categorized into four main families based on their function and sequence, Gs, Gi/o, Gq/11, and G12/13. One receptor can couple to more than one G protein subtype, and the coupling efficiency varies depending on the GPCR-G protein pair. However, the precise mechanism underlying different coupling efficiencies is unknown. Here, we study the structural mechanism underlying primary and secondary Gi/o coupling, using the muscarinic acetylcholine receptor type 2 (M2R) as the primary Gi/o-coupling receptor and the β₂-adrenergic receptor (β₂AR, which primarily couples to Gs) as the secondary Gi/o-coupling receptor. Hydrogen/deuterium exchange mass spectrometry and mutagenesis studies reveal that the engagement of the distal C-terminus of Gαi/o with the receptor differentiates primary and secondary Gi/o couplings. This study suggests that the conserved hydrophobic residue within the intracellular loop 2 of the receptor (residue 34.51) is not critical for primary Gi/o-coupling; however, it might be important for secondary Gi/o-coupling.

G protein-coupled receptors (GPCRs) can couple to more than one G protein subtype, and the coupling efficiency varies depending on the GPCR-G protein pair. Here authors use hydrogen/deuterium exchange mass spectrometry and mutagenesis to study the structural mechanism underlying primary and secondary Gi/o coupling.

Structural basis of Gs and Gi recognition by the human glucagon receptor

Class B G protein-coupled receptors, an important class of therapeutic targets, signal mainly through the G_s class of heterotrimeric G proteins, although they do display some promiscuity in G protein binding. Using cryo-electron microscopy, we determined the structures of the human glucagon receptor (GCGR) bound to glucagon and distinct classes of heterotrimeric G proteins, G_s or G_i1 These two structures adopt a similar open binding cavity to accommodate G_s and G_i1 The G_s binding selectivity of GCGR is explained by a larger interaction interface, but there are specific interactions that affect G_i more than G_s binding. Conformational differences in the receptor intracellular loops were found to be key selectivity determinants. These distinctions in transducer engagement were supported by mutagenesis and functional studies.

Nanobody-enabled monitoring of kappa opioid receptor states

Recent studies show that GPCRs rapidly interconvert between multiple states although our ability to interrogate, monitor and visualize them is limited by a relative lack of suitable tools. We previously reported two nanobodies (Nb39 and Nb6) that stabilize distinct ligand- and efficacy-delimited conformations of the kappa opioid receptor. Here, we demonstrate via X-ray crystallography a nanobody-targeted allosteric binding site by which Nb6 stabilizes a ligand-dependent inactive state. As Nb39 stabilizes an active-like state, we show how these two state-dependent nanobodies can provide real-time reporting of ligand stabilized states in cells in situ. Significantly, we demonstrate that chimeric GPCRs can be created with engineered nanobody binding sites to report ligand-stabilized states. Our results provide both insights regarding potential mechanisms for allosterically modulating KOR with nanobodies and a tool for reporting the real-time, in situ dynamic range of GPCR activity.

Novel M2 -selective, Gi -biased agonists of muscarinic acetylcholine receptors

BACKGROUND AND PURPOSE

More than 30% of currently marketed medications act via GPCRs. Thus, GPCRs represent one of the most important pharmacotherapeutic targets. In contrast to traditional agonists activating multiple signalling pathways, agonists activating a single signalling pathway represent a new generation of drugs with increased specificity and fewer adverse effects.

EXPERIMENTAL APPROACH

We have synthesized novel agonists of muscarinic ACh receptors and tested their binding and function (on levels of cAMP and inositol phosphates) in CHO cells expressing individual subtypes of muscarinic receptors, primary cultures of rat aortic smooth muscle cells and suspensions of digested native tissues from rats. Binding of the novel compounds to M₂ receptors was modelled in silico.

KEY RESULTS

Two of the tested new compounds (1-(thiophen-2-ylmethyl)-3,6-dihydro-2H-pyridinium and 1-methyl-1-(thiophen-2-ylmethyl)-3,6-dihydro-2H-pyridinium) only inhibited cAMP synthesis in CHO cells, primary cultures, and native tissues, with selectivity for M₂ muscarinic receptors and displaying bias towards the G_i signalling pathway at all subtypes of muscarinic receptors. Molecular modelling revealed interactions with the orthosteric binding site in a way specific for a given agonist followed by agonist-specific changes in the conformation of the receptor.

CONCLUSIONS AND IMPLICATIONS

The identified compounds may serve as lead structures in the search for novel non-steroidal and non-opioid analgesics acting via M₂ and M₄ muscarinic receptors with reduced side effects associated with activation of the phospholipase C signalling pathway.

Structures of the M1 and M2 muscarinic acetylcholine receptor/G-protein complexes

Muscarinic acetylcholine receptors are G protein-coupled receptors that respond to acetylcholine and play important signaling roles in the nervous system. There are five muscarinic receptor subtypes (M1R to M5R), which, despite sharing a high degree of sequence identity in the transmembrane region, couple to different heterotrimeric GTP-binding proteins (G proteins) to transmit signals. M1R, M3R, and M5R couple to the G_{q/ 11} family, whereas M2R and M4R couple to the G_{i/ o} family. Here, we present and compare the cryo-electron microscopy structures of M1R in complex with G₁₁ and M2R in complex with G_oA The M1R-G₁₁ complex exhibits distinct features, including an extended transmembrane helix 5 and carboxyl-terminal receptor tail that interacts with G protein. Detailed analysis of these structures provides a framework for understanding the molecular determinants of G-protein coupling selectivity.

Improved Modeling of Peptide-Protein Binding Through Global Docking and Accelerated Molecular Dynamics Simulations

Peptides mediate up to 40% of known protein-protein interactions in higher eukaryotes and play a key role in cellular signaling, protein trafficking, immunology, and oncology. However, it is challenging to predict peptide-protein binding with conventional computational modeling approaches, due to slow dynamics and high peptide flexibility. Here, we present a prototype of the approach which combines global peptide docking using ClusPro PeptiDock and all-atom enhanced simulations using Gaussian accelerated molecular dynamics (GaMD). For three distinct model peptides, the lowest backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures obtained from PeptiDock were 3.3–4.8 Å, being medium quality predictions according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. GaMD simulations refined the peptide-protein complex structures with significantly reduced peptide backbone RMSDs of 0.6–2.7 Å, yielding two high quality (sub-angstrom) and one medium quality models. Furthermore, the GaMD simulations identified important low-energy conformational states and revealed the mechanism of peptide binding to the target proteins. Therefore, PeptiDock+GaMD is a promising approach for exploring peptide-protein interactions.

Non-Functional Trace Amine-Associated Receptor 1 Variants in Patients With Mental Disorders

Background: The G protein–coupled receptor (GPCR) trace amine-associated receptor 1 (TAAR1) is expressed across brain areas involved in emotions, reward and cognition, and modulates monoaminergic and glutamatergic neurotransmissions. TAAR1 is stimulated with nanomolar affinity by 3-iodothyronamine (T1AM), an endogenous messenger considered a novel branch of thyroid hormone signaling. The human gene for TAAR1 maps to locus 6q23, within a region associated with major mental disorders.

Materials and Methods: We screened a cohort of patients with major mental disorders (n = 104) and a group of healthy controls (n = 130) for TAAR1 variants. HEK293 cells were transiently transfected with: i) wild-type TAAR1 and ii) mutated TAAR1, either in homozygous or heterozygous state. Cell surface expression and Gs/adenylyl cyclase activation upon administration of β-phenylethylamine (PEA), T1AM, and RO5166017, were assessed.

Results: We detected 13 missense variants in TAAR1 coding region, with a significant enrichment in patients as compared to healthy controls (11 vs. 1, 1 variant in both groups, p < 0.01). In silico analysis identified four dysfunctional variants, all in patients. Three of these—R23C, Y131C, and C263R—were functionally characterized. In cells co-transfected with wild-type and mutated TAAR1, we observed a significant reduction of cell surface expression. In heterozygosity, the three TAAR1 variants substantially dampened Gs signaling in response to PEA, and, more robustly, to T1AM. Co-stimulation with PEA and RO5166017 did not yield any improvement in Gs signaling. R23C, Y131C, and C263R are rare in the general population and map in functionally important highly conserved positions across TAAR1 orthologous and paralogous genes.

Conclusions: Our findings suggest that disruptions of TAAR1 activity may be relevant to the pathophysiology of mental disorders, thereby providing a promising target for novel psychopharmacological interventions.

Mechanistic Insights into Specific G Protein Interactions with Adenosine Receptors

Coupling between G-protein-coupled receptors (GPCRs) and the G proteins is a key step in cellular signaling. Despite extensive experimental and computational studies, the mechanism of specific GPCR-G protein coupling remains poorly understood. This has greatly hindered effective drug design of GPCRs that are primary targets of ∼¹/₃ of currently marketed drugs. Here, we have employed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method to decipher the mechanism of the GPCR-G protein interactions. Adenosine receptors (ARs) were used as model systems based on very recently determined cryo-EM structures of the A₁AR and A_2AAR coupled with the G_i and G_s proteins, respectively. Changing the G_i protein to the G_s led to increased fluctuations in the A₁AR and agonist adenosine (ADO), while agonist 5'-N-ethylcarboxamidoadenosine (NECA) binding in the A_2AAR could be still stabilized upon changing the G_s protein to the G_i. Free energy calculations identified one stable low-energy conformation for each of the A₁AR-G_i and A_2AAR-G_s complexes as in the cryo-EM structures, similarly for the A_2AAR-G_i complex. In contrast, the ADO agonist and G_s protein sampled multiple conformations in the A₁AR-G_s system. GaMD simulations thus indicated that the A₁AR preferred to couple with the G_i protein to the G_s, while the A_2AAR could couple with both the G_s and G_i proteins, being highly consistent with experimental findings of the ARs. More importantly, detailed analysis of the atomic simulations showed that the specific AR-G protein coupling resulted from remarkably complementary residue interactions at the protein interface, involving mainly the receptor transmembrane 6 helix and the G_α α5 helix and α4-β6 loop. In summary, the GaMD simulations have provided unprecedented insights into the dynamic mechanism of specific GPCR-G protein interactions at an atomistic level.

Conformational plasticity of the intracellular cavity of GPCR-G-protein complexes leads to G-protein promiscuity and selectivity

While the dynamics of the intracellular surface in agonist-stimulated GPCRs is well studied, the impact of GPCR dynamics on G-protein selectivity remains unclear. Here, we combine molecular dynamics simulations with live-cell FRET and secondary messenger measurements, for 21 GPCR-G-protein combinations, to advance a dynamic model of the GPCR-G-protein interface. Our data show C terminus peptides of Gαs, Gαi, and Gαq proteins assume a small ensemble of unique orientations when coupled to their cognate GPCRs, similar to the variations observed in 3D structures of GPCR-G-protein complexes. The noncognate G proteins interface with latent intracellular GPCR cavities but dissociate due to weak and unstable interactions. Three predicted mutations in β2-adrenergic receptor stabilize binding of noncognate Gαq protein in its latent cavity, allowing promiscuous signaling through both Gαs and Gαq in a dose-dependent manner. This demonstrates that latent GPCR cavities can be evolved, by design or nature, to tune G-protein selectivity, giving insights to pluridimensional GPCR signaling.

Variable G protein determinants of GPCR coupling selectivity

G protein-coupled receptors (GPCRs) regulate a wide variety of important cellular processes and are targeted by a large fraction of approved drugs. GPCRs signal by activating heterotrimeric G proteins and must couple to a select subset of G proteins to produce appropriate intracellular responses. It is not known how GPCRs select G proteins, but it is generally accepted that the Gα subunit C terminus is the primary G protein determinant of coupling selectivity. We systematically studied coupling of GPCRs to four families of G proteins and chimeras with C terminal regions that were exchanged between families. We uncovered rules for coupling selectivity and found that different GPCRs can recognize different features of the same G protein for selective coupling.

G protein-coupled receptors (GPCRs) activate four families of heterotrimeric G proteins, and individual receptors must select a subset of G proteins to produce appropriate cellular responses. Although the precise mechanisms of coupling selectivity are uncertain, the Gα subunit C terminus is widely believed to be the primary determinant recognized by cognate receptors. Here, we directly assess coupling between 14 representative GPCRs and 16 Gα subunits, including one wild-type Gα subunit from each of the four families and 12 chimeras with exchanged C termini. We use a sensitive bioluminescence resonance energy transfer (BRET) assay that provides control over both ligand and nucleotide binding, and allows direct comparison across G protein families. We find that the G_s- and G_q-coupled receptors we studied are relatively promiscuous and always couple to some extent to G_i1 heterotrimers. In contrast, G_i-coupled receptors are more selective. Our results with Gα subunit chimeras show that the Gα C terminus is important for coupling selectivity, but no more so than the Gα subunit core. The relative importance of the Gα subunit core and C terminus is highly variable and, for some receptors, the Gα core is more important for selective coupling than the C terminus. Our results suggest general rules for GPCR-G protein coupling and demonstrate that the critical G protein determinants of selectivity vary widely, even for different receptors that couple to the same G protein.

Biased signaling of G protein coupled receptors (GPCRs): Molecular determinants of GPCR/transducer selectivity and therapeutic potential

G protein coupled receptors (GPCRs) convey signals across membranes via interaction with G proteins. Originally, an individual GPCR was thought to signal through one G protein family, comprising cognate G proteins that mediate canonical receptor signaling. However, several deviations from canonical signaling pathways for GPCRs have been described. It is now clear that GPCRs can engage with multiple G proteins and the line between cognate and non-cognate signaling is increasingly blurred. Furthermore, GPCRs couple to non-G protein transducers, including β-arrestins or other scaffold proteins, to initiate additional signaling cascades. Receptor/transducer selectivity is dictated by agonist-induced receptor conformations as well as by collateral factors. In particular, ligands stabilize distinct receptor conformations to preferentially activate certain pathways, designated 'biased signaling'. In this regard, receptor sequence alignment and mutagenesis have helped to identify key receptor domains for receptor/transducer specificity. Furthermore, molecular structures of GPCRs bound to different ligands or transducers have provided detailed insights into mechanisms of coupling selectivity. However, receptor dimerization, compartmentalization, and trafficking, receptor-transducer-effector stoichiometry, and ligand residence and exposure times can each affect GPCR coupling. Extrinsic factors including cell type or assay conditions can also influence receptor signaling. Understanding these factors may lead to the development of improved biased ligands with the potential to enhance therapeutic benefit, while minimizing adverse effects. In this review, evidence for ligand-specific GPCR signaling toward different transducers or pathways is elaborated. Furthermore, molecular determinants of biased signaling toward these pathways and relevant examples of the potential clinical benefits and pitfalls of biased ligands are discussed.

Computational modeling of the olfactory receptor Olfr73 suggests a molecular basis for low potency of olfactory receptor-activating compounds

The mammalian olfactory system uses hundreds of specialized G-protein-coupled olfactory receptors (ORs) to discriminate a nearly unlimited number of odorants. Cognate agonists of most ORs have not yet been identified and potential non-olfactory processes mediated by ORs are unknown. Here, we used molecular modeling, fingerprint interaction analysis and molecular dynamics simulations to show that the binding pocket of the prototypical olfactory receptor Olfr73 is smaller, but more flexible, than binding pockets of typical non-olfactory G-protein-coupled receptors. We extended our modeling to virtual screening of a library of 1.6 million compounds against Olfr73. Our screen predicted 25 Olfr73 agonists beyond traditional odorants, of which 17 compounds, some with therapeutic potential, were validated in cell-based assays. Our modeling suggests a molecular basis for reduced interaction contacts between an odorant and its OR and thus the typical low potency of OR-activating compounds. These results provide a proof-of-principle for identifying novel therapeutic OR agonists.

State-dependent Lipid Interactions with the A2a Receptor Revealed by MD Simulations Using In Vivo-Mimetic Membranes

Membranes are known to have modulatory effects on G protein-coupled receptors (GPCRs) via specific lipid interactions. However, the mechanisms of such modulations in physiological conditions and how they influence GPCR functions remain unclear. Here we report coarse-grained molecular dynamics simulations on the Adenosine A2a receptor in different conformational states embedded in an in vivo-mimetic membrane model. Nine lipid interaction sites were revealed. The strength of lipid interactions with these sites showed a degree of dependence on the conformational states of the receptor, suggesting that these lipids may regulate the conformational dynamics of the receptor. In particular, we revealed a dual role of PIP₂ on A2aR activation that involves both stabilization of the characteristic outward tilt of TM6 and enhancement of A2aR-mini-Gs association. Our results demonstrated that the bound lipids allosterically regulate the functional properties of GPCRs. These protein-lipid interactions provide a springboard for design of allosteric modulators of GPCRs.

Recent advances in computational studies of GPCR-G protein interactions

Protein-protein interactions are key in cellular signaling. G protein-coupled receptors (GPCRs), the largest superfamily of human membrane proteins, are able to transduce extracellular signals (e.g., hormones and neurotransmitters) to intracellular proteins, in particular the G proteins. Since GPCRs serve as primary targets of ~1/3 of currently marketed drugs, it is important to understand mechanisms of GPCR signaling in order to design selective and potent drug molecules. This chapter focuses on recent advances in computational studies of the GPCR-G protein interactions using bioinformatics, protein-protein docking and molecular dynamics simulation approaches.

Signal Transduction and Pathogenic Modifications at the Melanocortin-4 Receptor: A Structural Perspective

The melanocortin-4 receptor (MC4R) can be endogenously activated by binding of melanocyte-stimulating hormones (MSH), which mediates anorexigenic effects. In contrast, the agouti-related peptide (AgRP) acts as an endogenous inverse agonist and suppresses ligand-independent basal signaling activity (orexigenic effects). Binding of ligands to MC4R leads to the activation of different G-protein subtypes or arrestin and concomitant signaling pathways. This receptor is a key protein in the hypothalamic regulation of food intake and energy expenditure and naturally-occurring inactivating MC4R variants are the most frequent cause of monogenic obesity. In general, obesity is a growing problem on a global scale and is of social, medical, and economic relevance. A significant goal is to develop optimized pharmacological tools targeting MC4R without adverse effects. To date, this has not been achieved because of inter alia non-selective ligands across the five functionally different MCR subtypes (MC1-5R). This motivates further investigation of (i) the three-dimensional MC4R structure, (ii) binding mechanisms of various ligands, and (iii) the molecular transfer process of signal transduction, with the aim of understanding how structural features are linked with functional-physiological aspects. Unfortunately, experimentally elucidated structural information is not yet available for the MC receptors, a group of class A G-protein coupled receptors (GPCRs). We, therefore, generated MC4R homology models and complexes with interacting partners to describe approximate structural properties associated with signaling mechanisms. In addition, molecular insights from pathogenic mutations were incorporated to discriminate more precisely their individual malfunction of the signal transfer mechanism.

Role of membrane cholesterol in differential sensitivity of muscarinic receptor subtypes to persistently bound xanomeline

Xanomeline (3-(Hexyloxy)-4-(1-methyl-1,2,5,6-tetrahydropyridin-3-yl)-1,2,5-thiadiazole) is a muscarinic agonist that is considered to be functionally selective for the M₁/M₄ receptor subtypes. Part of xanomeline binding is resistant to washing. Wash-resistant xanomeline activates muscarinic receptors persistently, except for the M₅ subtype. Mutation of leucine 6.46 to isoleucine at M₁ or M₄ receptors abolished persistent activation by wash-resistant xanomeline. Reciprocal mutation of isoleucine 6.46 to leucine at the M₅ receptor made it sensitive to activation by wash-resistant xanomeline. Lowering of membrane cholesterol made M₁ and M₄ mutants and M₅ wild type receptors sensitive to activation by wash-resistant xanomeline. Molecular docking revealed a cholesterol binding site in the groove between transmembrane helices 6 and 7. Molecular dynamics showed that interaction of cholesterol with this binding site attenuates receptor activation. We hypothesize that differences in cholesterol binding to this site between muscarinic receptor subtypes may constitute the basis for xanomeline apparent functional selectivity and may have notable therapeutic implications. Differences in receptor-membrane interactions, rather than in agonist-receptor interactions, represent a novel possibility to achieve pharmacological selectivity. Our findings may be applicable to other G protein coupled receptors.

Rules of Engagement: GPCRs and G Proteins

G protein-coupled receptors (GPCRs) are a key drug target class. They account for over one-third of current pharmaceuticals, and both drugs that inhibit and promote receptor function are important therapeutically; in some cases, the same GPCR can be targeted with agonists and inhibitors, depending upon disease context. There have been major breakthroughs in understanding GPCR structure and drug binding through advances in X-ray crystallography, and membrane protein stabilization. Nonetheless, these structures have predominately been of inactive receptors bound to inhibitors. Efforts to capture structures of fully active GPCRs, in particular those in complex with the canonical, physiological transducer G protein, have been limited via this approach. Very recently, advances in cryo-electron microscopy have provided access to agonist:GPCR:G protein complex structures. These promise to revolutionize our understanding of GPCR:G protein engagement and provide insight into mechanisms of efficacy and coupling selectivity and how these might be controlled by biased agonists. Here we review what we have currently learned from the new GPCR:Gs and GPCR:Gi/o complex structures.

Development of an antibody fragment that stabilizes GPCR/G-protein complexes

Single-particle cryo-electron microscopy (cryo-EM) has recently enabled high-resolution structure determination of numerous biological macromolecular complexes. Despite this progress, the application of high-resolution cryo-EM to G protein coupled receptors (GPCRs) in complex with heterotrimeric G proteins remains challenging, owning to both the relative small size and the limited stability of these assemblies. Here we describe the development of antibody fragments that bind and stabilize GPCR-G protein complexes for the application of high-resolution cryo-EM. One antibody in particular, mAb16, stabilizes GPCR/G-protein complexes by recognizing an interface between Gα and Gβγ subunits in the heterotrimer, and confers resistance to GTPγS-triggered dissociation. The unique recognition mode of this antibody makes it possible to transfer its binding and stabilizing effect to other G-protein subtypes through minimal protein engineering. This antibody fragment is thus a broadly applicable tool for structural studies of GPCR/G-protein complexes.

The determination of high resolution structures of G protein coupled receptors (GPCRs) in complex with heterotrimeric G proteins is challenging. Here authors develop an antibody fragment, mAB16, which stabilizes GPCR/G-protein complexes and facilitates the application of high resolution cryo-EM.

Crystal structure of rhodopsin in complex with a mini-Go sheds light on the principles of G protein selectivity

Selective coupling of G protein (heterotrimeric guanine nucleotide-binding protein)-coupled receptors (GPCRs) to specific Gα-protein subtypes is critical to transform extracellular signals, carried by natural ligands and clinical drugs, into cellular responses. At the center of this transduction event lies the formation of a signaling complex between the receptor and G protein. We report the crystal structure of light-sensitive GPCR rhodopsin bound to an engineered mini-G_o protein. The conformation of the receptor is identical to all previous structures of active rhodopsin, including the complex with arrestin. Thus, rhodopsin seems to adopt predominantly one thermodynamically stable active conformation, effectively acting like a "structural switch," allowing for maximum efficiency in the visual system. Furthermore, our analysis of the well-defined GPCR-G protein interface suggests that the precise position of the carboxyl-terminal "hook-like" element of the G protein (its four last residues) relative to the TM7/helix 8 (H8) joint of the receptor is a significant determinant in selective G protein activation.

Exploring a new ligand binding site of G protein-coupled receptors

Identifying a target ligand binding site is an important step for structure-based rational drug design as shown here for G protein-coupled receptors (GPCRs), which are among the most popular drug targets. We applied long-time scale molecular dynamics simulations, coupled with mutagenesis studies, to two prototypical GPCRs, the M3 and M4 muscarinic acetylcholine receptors. Our results indicate that unlike synthetic antagonists, which bind to the classic orthosteric site, the endogenous agonist acetylcholine is able to diffuse into a much deeper binding pocket. We also discovered that the most recently resolved crystal structure of the LTB4 receptor comprised a bound inverse agonist, which extended its benzamidine moiety to the same binding pocket discovered in this work. Analysis on all resolved GPCR crystal structures indicated that this new pocket could exist in most receptors. Our findings provide new opportunities for GPCR drug discovery.