Efficacy as an Intrinsic Property of the M2 Muscarinic Receptor in Its Tetrameric State

Preferential Gs protein coupling of the galanin Gal1 receptor in the µ-opioid-Gal1 receptor heterotetramer

Recent studies have proposed that heteromers of µ-opioid receptors (MORs) and galanin Gal₁ receptors (Gal₁Rs) localized in the mesencephalon mediate the dopaminergic effects of opioids. The present study reports converging evidence, using a peptide-interfering approach combined with biophysical and biochemical techniques, including total internal reflection fluorescence microscopy, for a predominant homodimeric structure of MOR and Gal₁R when expressed individually, and for their preference to form functional heterotetramers when co-expressed. Results show that a heteromerization-dependent change in the Gal₁R homodimeric interface leads to a switch in G-protein coupling from inhibitory Gi to stimulatory Gs proteins. The MOR-Gal₁R heterotetramer, which is thus bound to Gs via the Gal₁R homodimer and Gi via the MOR homodimer, provides the framework for a canonical Gs-Gi antagonist interaction at the adenylyl cyclase level. These novel results shed light on the intense debate about the oligomeric quaternary structure of G protein-coupled receptors, their predilection for heteromer formation, and the resulting functional significance.

Allostery of atypical modulators at oligomeric G protein-coupled receptors

Many G protein-coupled receptors (GPCRs) are therapeutic targets, with most drugs acting at the orthosteric site. Some GPCRs also possess allosteric sites, which have become a focus of drug discovery. In the M2 muscarinic receptor, allosteric modulators regulate the binding and functional effects of orthosteric ligands through a mix of conformational changes, steric hindrance and electrostatic repulsion transmitted within and between the constituent protomers of an oligomer. Tacrine has been called an atypical modulator because it exhibits positive cooperativity, as revealed by Hill coefficients greater than 1 in its negative allosteric effect on binding and response. Radioligand binding and molecular dynamics simulations were used to probe the mechanism of that modulation in monomers and oligomers of wild-type and mutant M2 receptors. Tacrine is not atypical at monomers, which indicates that its atypical effects are a property of the receptor in its oligomeric state. These results illustrate that oligomerization of the M2 receptor has functional consequences.

BRET- and fluorescence anisotropy-based assays for real-time monitoring of ligand binding to M2 muscarinic acetylcholine receptors

BRET and fluorescence anisotropy (FA) are two fluorescence-based techniques used for the characterization of ligand binding to G protein-coupled receptors (GPCRs) and both allow monitoring of ligand binding in real time. In this study, we present the first direct comparison of BRET-based and FA-based binding assays using the human M₂ muscarinic acetylcholine receptor (M₂R) and two TAMRA (5-carboxytetramethylrhodamine)-labeled fluorescent ligands as a model system. The determined fluorescent ligand affinities from both assays were in good agreement with results obtained from radioligand competition binding experiments. The assays yielded real-time kinetic binding data revealing differences in the mechanism of binding for the investigated fluorescent probes. Furthermore, the investigation of various unlabeled M₂R ligands yielded pharmacological profiles in accordance with earlier reported data. Taken together, this study showed that BRET- and FA-based binding assays represent valuable alternatives to radioactivity-based methods for screening purposes and for a precise characterization of binding kinetics supporting the exploration of binding mechanisms.

Shaky ground - The nature of metastable GPCR signalling complexes

How G protein-coupled receptors (GPCR) interact with one another remains an area of active investigation. Obligate dimers of class C GPCRs such as metabotropic GABA and glutamate receptors are well accepted, although whether this is a general feature of other GPCRs is still strongly debated. In this review, we focus on the idea that GPCR dimers and oligomers are better imagined as parts of larger metastable signalling complexes. We discuss the nature of functional oligomeric entities, their stabilities and kinetic features and how structural and functional asymmetries of such metastable entities might have implications for drug discovery. This article is part of the Special Issue entitled 'Receptor heteromers and their allosteric receptor-receptor interactions'.

Gaussian accelerated molecular dynamics for elucidation of drug pathways

INTRODUCTION

Understanding pathways and mechanisms of drug binding to receptors is important for rational drug design. Remarkable advances in supercomputing and methodological developments have opened a new era for application of computer simulations in predicting drug-receptor interactions at an atomistic level. Gaussian accelerated molecular dynamics (GaMD) is a computational enhanced sampling technique that works by adding a harmonic boost potential to reduce energy barriers. GaMD enables free energy calculations without the requirement of predefined collective variables. GaMD has proven useful in biomolecular simulations, in particular, the prediction of drug-receptor interactions. Areas covered: Herein, the authors review recent GaMD simulation studies that elucidated pathways of drug binding to proteins including the G-protein-coupled receptors and HIV protease. Expert opinion: GaMD is advantageous for enhanced simulations of, amongst many biological processes, drug binding to target receptors. Compared with conventional molecular dynamics, GaMD speeds up biomolecular simulations by orders of magnitude. GaMD enables routine drug binding simulations using personal computers with GPUs or common computing clusters. GaMD and, more broadly, enhanced sampling simulations are expected to dramatically increase our capabilities to determine the mechanisms of drug binding to a wide range of receptors in the near future. This will greatly facilitate computer-aided drug design.

Ligand-Induced Coupling between Oligomers of the M2 Receptor and the Gi1 Protein in Live Cells

Uncertainty over the mechanism of signaling via G protein-coupled receptors (GPCRs) relates in part to questions regarding their supramolecular structure. GPCRs and heterotrimeric G proteins are known to couple as monomers under various conditions. Many GPCRs form oligomers under many of the same conditions, however, and the biological role of those complexes is unclear. We have used dual-color fluorescence correlation spectroscopy to identify oligomers of the M₂ muscarinic receptor and of G_i1 in purified preparations and live Chinese hamster ovary cells. Measurements on differently tagged receptors (i.e., eGFP-M₂ and mCherry-M₂) and G proteins (i.e., eGFP-Gα_i1β₁γ₂ and mCherry-Gα_i1β₁γ₂) detected significant cross-correlations between the two fluorophores in each case, both in detergent micelles and in live cells, indicating that both the receptor and G_i1 can exist as homo-oligomers. Oligomerization of differently tagged G_i1 decreased upon the activation of co-expressed wild-type M₂ receptor by an agonist. Measurements on a tagged M₂ receptor (M₂-mCherry) and eGFP-Gα_i1β₁γ₂ co-expressed in live cells detected cross-correlations only in the presence of an agonist, which therefore promoted coupling of the receptor and the G protein. The effect of the agonist was retained when a fluorophore-tagged receptor lacking the orthosteric site (i.e., M₂(D103A)-mCherry) was co-expressed with the wild-type receptor and eGFP-Gα_i1β₁γ₂, indicating that the ligand acted via an oligomeric receptor. Our results point to a model in which an agonist promotes transient coupling of otherwise independent oligomers of the M₂ receptor on the one hand and of G_i1 on the other and that an activated complex leads to a reduction in the oligomeric size of the G protein. They suggest that GPCR-mediated signaling proceeds, at least in part, via oligomers.

Muscarinic receptor oligomerization

G protein-coupled receptors (GPCRs) have been classically described as monomeric entities that function by binding in a 1:1 stoichiometric ratio to both ligand and downstream signalling proteins. However, in recent years, a growing number of studies has supported the hypothesis that these receptors can interact to form dimers and higher order oligomers although the molecular basis for these interactions, the overall quaternary arrangements and the functional importance of GPCR oligomerization remain topics of intense speculation. Muscarinic acetylcholine receptors belong to class A of the GPCR family. Each muscarinic receptor subtype has its own particular distribution throughout the central and peripheral nervous systems. In the central nervous system, muscarinic receptors regulate several sensory, cognitive, and motor functions while, in the peripheral nervous system, they are involved in the regulation of heart rate, stimulation of glandular secretion and smooth muscle contraction. Muscarinic acetylcholine receptors have long been used as a model for the study of GPCR structure and function and to address aspects of GPCR dimerization using a broad range of approaches. In this review, the prevailing knowledge regarding the quaternary arrangement for the various muscarinic acetylcholine receptors has been summarized by discussing work ranging from initial results obtained using more traditional biochemical approaches to those generated with more modern biophysical techniques. This article is part of the Special Issue entitled 'Neuropharmacology on Muscarinic Receptors'.

The Dynamics of GPCR Oligomerization and Their Functional Consequences

The functional importance of G protein-coupled receptor (GPCR) oligomerization remains controversial. Although obligate dimers of class C GPCRs are well accepted, the generalizability of this phenomenon is still strongly debated with respect to other classes of GPCRs. In this review, we focus on understanding the organization and dynamics between receptor equivalents and their signaling partners in oligomeric receptor complexes, with a view toward integrating disparate viewpoints into a unified understanding. We discuss the nature of functional oligomeric entities, and how asymmetries in receptor structure and function created by oligomers might have implications for receptor function as allosteric machines and for future drug discovery.

Spatial Intensity Distribution Analysis: Studies of G Protein-Coupled Receptor Oligomerisation

Spatial intensity distribution analysis (SpIDA) is a recently developed approach for determining quaternary structure information on fluorophore-labelled proteins of interest in situ. It can be applied to live or fixed cells and native tissue. Using confocal images, SpIDA generates fluorescence intensity histograms that are analysed by super-Poissonian distribution functions to obtain density and quantal brightness values of the fluorophore-labelled protein of interest. This allows both expression level and oligomerisation state of the protein to be determined. We describe the application of SpIDA to investigate the oligomeric state of G protein-coupled receptors (GPCRs) at steady state and following cellular challenge, and consider how SpIDA may be used to explore GPCR quaternary organisation in pathophysiology and to stratify medicines.

GPCRs may exist and function as monomers: however, abundant evidence suggests they can form dimers/oligomers.

This concept has implications for drug discovery as it may offer opportunities to modulate the effects of known pharmaceuticals or identify new drug therapies.

A variety of approaches have been applied to this issue from traditional biochemical techniques, via resonance energy transfer approaches to recently developed image analysis-based techniques such as SpIDA. This uses mathematical analysis of confocal microscopy images to generate quantal brightness and density information for a fluorophore-tagged receptor.

SpIDA can be applied to live or fixed cells and native tissue.

SpIDA has been applied to GPCRs from each of the major subfamilies to explore their oligomerisation status at steady state and their regulation by receptor density and ligand binding.

Gaussian Accelerated Molecular Dynamics: Theory, Implementation, and Applications

A novel Gaussian Accelerated Molecular Dynamics (GaMD) method has been developed for simultaneous unconstrained enhanced sampling and free energy calculation of biomolecules. Without the need to set predefined reaction coordinates, GaMD enables unconstrained enhanced sampling of the biomolecules. Furthermore, by constructing a boost potential that follows a Gaussian distribution, accurate reweighting of GaMD simulations is achieved via cumulant expansion to the second order. The free energy profiles obtained from GaMD simulations allow us to identify distinct low energy states of the biomolecules and characterize biomolecular structural dynamics quantitatively. In this chapter, we present the theory of GaMD, its implementation in the widely used molecular dynamics software packages (AMBER and NAMD), and applications to the alanine dipeptide biomolecular model system, protein folding, biomolecular large-scale conformational transitions and biomolecular recognition.

Signalling assemblies: the odds of symmetry

The assembly of proteins into complexes is fundamental to nearly all biological signalling processes. Symmetry is a dominant feature of the structures of experimentally determined protein complexes, observed in the vast majority of homomers and many heteromers. However, some asymmetric structures exist, and asymmetry also often forms transiently, intractable to traditional structure determination methods. Here, we explore the role of protein complex symmetry and asymmetry in cellular signalling, focusing on receptors, transcription factors and transmembrane channels, among other signalling assemblies. We highlight a recurrent tendency for asymmetry to be crucial for signalling function, often being associated with activated states. We conclude with a discussion of how consideration of protein complex symmetry and asymmetry has significant potential implications and applications for pharmacology and human disease.

Graded activation and free energy landscapes of a muscarinic G-protein-coupled receptor

G-protein-coupled receptors (GPCRs) recognize ligands of widely different efficacies, from inverse to partial and full agonists, which transduce cellular signals at differentiated levels. However, the mechanism of such graded activation remains unclear. Using the Gaussian accelerated molecular dynamics (GaMD) method that enables both unconstrained enhanced sampling and free energy calculation, we have performed extensive GaMD simulations (∼19 μs in total) to investigate structural dynamics of the M₂ muscarinic GPCR that is bound by the full agonist iperoxo (IXO), the partial agonist arecoline (ARC), and the inverse agonist 3-quinuclidinyl-benzilate (QNB), in the presence or absence of the G-protein mimetic nanobody. In the receptor-nanobody complex, IXO binding leads to higher fluctuations in the protein-coupling interface than ARC, especially in the receptor transmembrane helix 5 (TM5), TM6, and TM7 intracellular domains that are essential elements for GPCR activation, but less flexibility in the receptor extracellular region due to stronger binding compared with ARC. Two different binding poses are revealed for ARC in the orthosteric pocket. Removal of the nanobody leads to GPCR deactivation that is characterized by inward movement of the TM6 intracellular end. Distinct low-energy intermediate conformational states are identified for the IXO- and ARC-bound M₂ receptor. Both dissociation and binding of an orthosteric ligand are observed in a single all-atom GPCR simulation in the case of partial agonist ARC binding to the M₂ receptor. This study demonstrates the applicability of GaMD for exploring free energy landscapes of large biomolecules and the simulations provide important insights into the GPCR functional mechanism.

Variable Dependence of Signaling Output on Agonist Occupancy of Ste2p, a G Protein-coupled Receptor in Yeast

We report here on the relationship between ligand binding and signaling responses in the yeast pheromone response pathway, a well characterized G protein-coupled receptor system. Responses to agonist (α-factor) by cells expressing widely varying numbers of receptors depend primarily on fractional occupancy, not the absolute number of agonist-bound receptors. Furthermore, the concentration of competitive antagonist required to inhibit α-factor-dependent signaling is more than 10-fold higher than predicted based on the known ligand affinities. Thus, responses to a particular number of agonist-bound receptors can vary greatly, depending on whether there are unoccupied or antagonist-bound receptors present on the same cell surface. This behavior does not appear to be due to pre-coupling of receptors to G protein or to the Sst2p regulator of G protein signaling. The results are consistent with a signaling response that is determined by the integration of positive signals from agonist-occupied receptors and inhibitory signals from unoccupied receptors, where the inhibitory signals can be diminished by antagonist binding.

Quaternary structures of opsin in live cells revealed by FRET spectrometry

Rhodopsin is a prototypical G-protein-coupled receptor (GPCR) that initiates phototransduction in the retina. The receptor consists of the apoprotein opsin covalently linked to the inverse agonist 11-cis retinal. Rhodopsin and opsin have been shown to form oligomers within the outer segment disc membranes of rod photoreceptor cells. However, the physiological relevance of the observed oligomers has been questioned since observations were made on samples prepared from the retina at low temperatures. To investigate the oligomeric status of opsin in live cells at body temperatures, we utilized a novel approach called Förster resonance energy transfer spectrometry, which previously has allowed the determination of the stoichiometry and geometry (i.e. quaternary structure) of various GPCRs. In the current study, we have extended the method to additionally determine whether or not a mixture of oligomeric forms of opsin exists and in what proportion. The application of this improved method revealed that opsin expressed in live Chinese hamster ovary (CHO) cells at 37°C exists as oligomers of various sizes. At lower concentrations, opsin existed in an equilibrium of dimers and tetramers. The tetramers were in the shape of a near-rhombus. At higher concentrations of the receptor, higher-order oligomers began to form. Thus, a mixture of different oligomeric forms of opsin is present in the membrane of live CHO cells and oligomerization occurs in a concentration-dependent manner. The general principles underlying the concentration-dependent oligomerization of opsin may be universal and apply to other GPCRs as well.

Allosteric modulation in monomers and oligomers of a G protein-coupled receptor

The M₂ muscarinic receptor is the prototypic model of allostery in GPCRs, yet the molecular and the supramolecular determinants of such effects are unknown. Monomers and oligomers of the M₂ muscarinic receptor therefore have been compared to identify those allosteric properties that are gained in oligomers. Allosteric interactions were monitored by means of a FRET-based sensor of conformation at the allosteric site and in pharmacological assays involving mutants engineered to preclude intramolecular effects. Electrostatic, steric, and conformational determinants of allostery at the atomic level were examined in molecular dynamics simulations. Allosteric effects in monomers were exclusively negative and derived primarily from intramolecular electrostatic repulsion between the allosteric and orthosteric ligands. Allosteric effects in oligomers could be positive or negative, depending upon the allosteric-orthosteric pair, and they arose from interactions within and between the constituent protomers. The complex behavior of oligomers is characteristic of muscarinic receptors in myocardial preparations.

DOI:http://dx.doi.org/10.7554/eLife.11685.001

Proteins called G protein-coupled receptors (GPCRs) are found on the surface of cells throughout the body. Hormones or other signal molecules – collectively known as ligands – from outside the cell can bind to the receptors to activate them. This causes a change in the structure of the receptor, which triggers a signal inside the cell to alter the cell’s behavior. GPCRs are known to form clusters of two or more receptor units, but it is not known if these clusters have unique properties or what role they play in cells.

Many drugs can bind to GPCRs and most of them block the activity of the receptors by taking the place of the natural ligand. Another way to alter the activity of a GPCR is with so-called 'allosteric' drugs. These bind to different sites on the receptor than the natural ligands do and can inhibit or enhance binding of the ligands by altering the shape of the receptor.

Shivnaraine et al. investigated how a type of GPCR called muscarinic cholinergic receptors interact within clusters. This involved developing a method to track the receptor in mammalian cells using a fluorescent sensor that detects changes in the allosteric site. The experiments show that two or more GPCRs need to interact for the receptors to respond to allosteric drugs in a manner that reflects the normal effect of the drugs on the body. This result is unexpected in light of the assumption that individual receptor molecules act independently. Shivnaraine et al.’s findings indicate that the clusters may play a role in the normal behavior of GPCRs in cells. A future challenge is to understand exactly how the GPCRs interact with each other.

DOI:http://dx.doi.org/10.7554/eLife.11685.002

Protease-activated receptors in hemostasis

Protease signaling in cells elicits multiple physiologically important responses via protease-activated receptors (PARs). There are 4 members of this family of G-protein-coupled receptors (PAR1-4). PARs are activated by proteolysis of the N terminus to reveal a tethered ligand. The rate-limiting step of PAR signaling is determined by the efficiency of proteolysis of the N terminus, which is regulated by allosteric binding sites, cofactors, membrane localization, and receptor dimerization. This ultimately controls the initiation of PAR signaling. In addition, these factors also control the cellular response by directing signaling toward G-protein or β-arrestin pathways. PAR1 signaling on endothelial cells is controlled by the activating protease and heterodimerization with PAR2 or PAR3. As a consequence, the genetic and epigenetic control of PARs and their cofactors in physiologic and pathophysiologic conditions have the potential to influence cellular behavior. Recent studies have uncovered polymorphisms that result in PAR4 sequence variants with altered reactivity that interact to influence platelet response. This further demonstrates how interactions within the plasma membrane can control the physiological output. Understanding the structural rearrangement following PAR activation and how PARs are allosterically controlled within the plasma membrane will determine how best to target this family of receptors therapeutically. The purpose of this article is to review how signaling from PARs is influenced by alternative cleavage sites and the physical interactions within the membrane. Going forward, it will be important to relate the altered signaling to the molecular arrangement of PARs in the cell membrane and to determine how these may be influenced genetically.

The GPCR heterotetramer: challenging classical pharmacology

Two concepts are gaining increasing acceptance in G protein-coupled receptor (GPCR) pharmacology: (i) pre-coupling of GPCRs with their preferred signaling molecules, and (ii) GPCR oligomerization. This is begging for the introduction of new models such as GPCR oligomer-containing signaling complexes with GPCR homodimers as functional building blocks. This model favors the formation of GPCR heterotetramers - heteromers of homodimers coupled to their cognate G protein. The GPCR heterotetramer offers an optimal framework for a canonical antagonistic interaction between activated Gs and Gi proteins, which can simultaneously bind to their respective preferred receptors and to adenylyl cyclase (AC) catalytic units. This review addresses the current evidence for pre-coupling of the various specific components that provide the very elaborate signaling machinery exemplified by the Gs-Gi-AC-coupled GPCR heterotetramer.

Oligomerization of Family B GPCRs: Exploration in Inter-Family Oligomer Formation

G-protein-coupled receptors (GPCRs) are classified into A to F subfamilies in which only families A, B, and C are present in mammals. Some of these GPCRs were found to form higher ordered structures such as oligomers with the discovery of interacting receptors in the form of homomers or heteromers. The importance of these oligomers on regulating receptor functions has recently been an intense research focus. It has been proposed that receptor oligomer formation has impact on its physiological importance on receptor trafficking, signaling, ligand-related regulation, and also is related to certain diseases. The present body of knowledge, however, comprises mainly intra-family oligomers formation and their consequences. Inter-family oligomers are recognized but there is limited information. This article aims to provide a current view regarding inter-family GPCR oligomerization in the subfamilies A, B, and C found in mammals.

Coupling of g proteins to reconstituted monomers and tetramers of the M2 muscarinic receptor

G protein-coupled receptors can be reconstituted as monomers in nanodiscs and as tetramers in liposomes. When reconstituted with G proteins, both forms enable an allosteric interaction between agonists and guanylyl nucleotides. Both forms, therefore, are candidates for the complex that controls signaling at the level of the receptor. To identify the biologically relevant form, reconstituted monomers and tetramers of the purified M2 muscarinic receptor were compared with muscarinic receptors in sarcolemmal membranes for the effect of guanosine 5'-[β,γ-imido]triphosphate (GMP-PNP) on the inhibition of N-[(3)H]methylscopolamine by the agonist oxotremorine-M. With monomers, a stepwise increase in the concentration of GMP-PNP effected a lateral, rightward shift in the semilogarithmic binding profile (i.e. a progressive decrease in the apparent affinity of oxotremorine-M). With tetramers and receptors in sarcolemmal membranes, GMP-PNP effected a vertical, upward shift (i.e. an apparent redistribution of sites from a state of high affinity to one of low affinity with no change in affinity per se). The data were analyzed in terms of a mechanistic scheme based on a ligand-regulated equilibrium between uncoupled and G protein-coupled receptors (the "ternary complex model"). The model predicts a rightward shift in the presence of GMP-PNP and could not account for the effects at tetramers in vesicles or receptors in sarcolemmal membranes. Monomers present a special case of the model in which agonists and guanylyl nucleotides interact within a complex that is both constitutive and stable. The results favor oligomers of the M2 receptor over monomers as the biologically relevant state for coupling to G proteins.