Molecular details of dimerization kinetics reveal negligible populations of transient µ-opioid receptor homodimers at physiological concentrations

Characterization of binding kinetics and intracellular signaling of new psychoactive substances targeting cannabinoid receptor using transition-based reweighting method

New psychoactive substances (NPS) targeting cannabinoid receptor 1 pose a significant threat to society as recreational abusive drugs that have pronounced physiological side effects. These greater adverse effects compared to classical cannabinoids have been linked to the higher downstream β-arrestin signaling. Thus, understanding the mechanism of differential signaling will reveal important structure-activity relationship essential for identifying and potentially regulating NPS molecules. In this study, we simulate the slow (un)binding process of NPS MDMB-Fubinaca and classical cannabinoid HU-210 from CB1 using multi-ensemble simulation to decipher the effects of ligand binding dynamics on downstream signaling. The transition-based reweighing method is used for the estimation of transition rates and underlying thermodynamics of (un)binding processes of ligands with nanomolar affinities. Our analyses reveal major interaction differences with transmembrane TM7 between NPS and classical cannabinoids. A variational autoencoder-based approach, neural relational inference (NRI), is applied to assess the allosteric effects on intracellular regions attributable to variations in binding pocket interactions. NRI analysis indicate a heightened level of allosteric control of NPxxY motif for NPS-bound receptors, which contributes to the higher probability of formation of a crucial triad interaction (Y7.53-Y5.58-T3.46) necessary for stronger β-arrestin signaling. Hence, in this work, MD simulation, data-driven statistical methods, and deep learning point out the structural basis for the heightened physiological side effects associated with NPS, contributing to efforts aimed at mitigating their public health impact.

Computational drug development for membrane protein targets

The application of computational biology in drug development for membrane protein targets has experienced a boost from recent developments in deep learning-driven structure prediction, increased speed and resolution of structure elucidation, machine learning structure-based design and the evaluation of big data. Recent protein structure predictions based on machine learning tools have delivered surprisingly reliable results for water-soluble and membrane proteins but have limitations for development of drugs that target membrane proteins. Structural transitions of membrane proteins have a central role during transmembrane signaling and are often influenced by therapeutic compounds. Resolving the structural and functional basis of dynamic transmembrane signaling networks, especially within the native membrane or cellular environment, remains a central challenge for drug development. Tackling this challenge will require an interplay between experimental and computational tools, such as super-resolution optical microscopy for quantification of the molecular interactions of cellular signaling networks and their modulation by potential drugs, cryo-electron microscopy for determination of the structural transitions of proteins in native cell membranes and entire cells, and computational tools for data analysis and prediction of the structure and function of cellular signaling networks, as well as generation of promising drug candidates.

Dimerization of melanocortin 4 receptor controls puberty onset and body size polymorphism

Xiphophorus fish exhibit a clear phenotypic polymorphism in puberty onset and reproductive strategies of males. In X. nigrensis and X. multilineatus, puberty onset is genetically determined and linked to a melanocortin 4 receptor (Mc4r) polymorphism of wild-type and mutant alleles on the sex chromosomes. We hypothesized that Mc4r mutant alleles act on wild-type alleles by a dominant negative effect through receptor dimerization, leading to differential intracellular signaling and effector gene activation. Depending on signaling strength, the onset of puberty either occurs early or is delayed. Here, we show by Förster Resonance Energy Transfer (FRET) that wild-type Xiphophorus Mc4r monomers can form homodimers, but also heterodimers with mutant receptors resulting in compromised signaling which explains the reduced Mc4r signaling in large males. Thus, hetero- vs. homo- dimerization seems to be the key molecular mechanism for the polymorphism in puberty onset and body size in male fish.

Computationally designed GPCR quaternary structures bias signaling pathway activation

Communication across membranes controls critical cellular processes and is achieved by receptors translating extracellular signals into selective cytoplasmic responses. While receptor tertiary structures can be readily characterized, receptor associations into quaternary structures are challenging to study and their implications in signal transduction remain poorly understood. Here, we report a computational approach for predicting receptor self-associations, and designing receptor oligomers with various quaternary structures and signaling properties. Using this approach, we designed chemokine receptor CXCR4 dimers with reprogrammed binding interactions, conformations, and abilities to activate distinct intracellular signaling proteins. In agreement with our predictions, the designed CXCR4s dimerized through distinct conformations and displayed different quaternary structural changes upon activation. Consistent with the active state models, all engineered CXCR4 oligomers activated the G protein Gi, but only specific dimer structures also recruited β-arrestins. Overall, we demonstrate that quaternary structures represent an important unforeseen mechanism of receptor biased signaling and reveal the existence of a bias switch at the dimer interface of several G protein-coupled receptors including CXCR4, mu-Opioid and type-2 Vasopressin receptors that selectively control the activation of G proteins vs β-arrestin-mediated pathways. The approach should prove useful for predicting and designing receptor associations to uncover and reprogram selective cellular signaling functions.

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.

The M1 muscarinic receptor is present in situ as a ligand-regulated mixture of monomers and oligomeric complexes

Although it is appreciated that members of the large family of rhodopsin-like cell surface receptors can form dimeric or larger protein complexes when expressed at high levels in cultured cells, their organizational state within native cells and tissues of the body is largely unknown. We assessed this in neurons of the central nervous system by replacing the M₁ muscarinic acetylcholine receptor in mice with a form of this receptor with an added fluorescent protein. Receptor function was unaltered by this change, and the biophysical approach we used demonstrated that the receptor exists as a mixture of monomers and dimers or oligomers. Drug treatments that target this receptor promote its monomerization, which may have significance for receptor function.

The quaternary organization of rhodopsin-like G protein-coupled receptors in native tissues is unknown. To address this we generated mice in which the M₁ muscarinic acetylcholine receptor was replaced with a C-terminally monomeric enhanced green fluorescent protein (mEGFP)–linked variant. Fluorescence imaging of brain slices demonstrated appropriate regional distribution, and using both anti-M₁ and anti–green fluorescent protein antisera the expressed transgene was detected in both cortex and hippocampus only as the full-length polypeptide. M₁-mEGFP was expressed at levels equal to the M₁ receptor in wild-type mice and was expressed throughout cell bodies and projections in cultured neurons from these animals. Signaling and behavioral studies demonstrated M₁-mEGFP was fully active. Application of fluorescence intensity fluctuation spectrometry to regions of interest within M₁-mEGFP–expressing neurons quantified local levels of expression and showed the receptor was present as a mixture of monomers, dimers, and higher-order oligomeric complexes. Treatment with both an agonist and an antagonist ligand promoted monomerization of the M₁-mEGFP receptor. The quaternary organization of a class A G protein-coupled receptor in situ was directly quantified in neurons in this study, which answers the much-debated question of the extent and potential ligand-induced regulation of basal quaternary organization of such a receptor in native tissue when present at endogenous expression levels.

Estimation of binding rates and affinities from multiensemble Markov models and ligand decoupling

Accurate and efficient simulation of the thermodynamics and kinetics of protein-ligand interactions is crucial for computational drug discovery. Multiensemble Markov Model (MEMM) estimators can provide estimates of both binding rates and affinities from collections of short trajectories but have not been systematically explored for situations when a ligand is decoupled through scaling of non-bonded interactions. In this work, we compare the performance of two MEMM approaches for estimating ligand binding affinities and rates: (1) the transition-based reweighting analysis method (TRAM) and (2) a Maximum Caliber (MaxCal) based method. As a test system, we construct a small host-guest system where the ligand is a single uncharged Lennard-Jones (LJ) particle, and the receptor is an 11-particle icosahedral pocket made from the same atom type. To realistically mimic a protein-ligand binding system, the LJ ϵ parameter was tuned, and the system was placed in a periodic box with 860 TIP3P water molecules. A benchmark was performed using over 80 µs of unbiased simulation, and an 18-state Markov state model was used to estimate reference binding affinities and rates. We then tested the performance of TRAM and MaxCal when challenged with limited data. Both TRAM and MaxCal approaches perform better than conventional Markov state models, with TRAM showing better convergence and accuracy. We find that subsampling of trajectories to remove time correlation improves the accuracy of both TRAM and MaxCal and that in most cases, only a single biased ensemble to enhance sampled transitions is required to make accurate estimates.

Modulation of adenosine A2a receptor oligomerization by receptor activation and PIP2 interactions

GPCRs have been shown to form oligomers, which generate distinctive signaling outcomes. However, the structural nature of the oligomerization process remains uncertain. We have characterized oligomeric configurations of the adenosine A2a receptor (A2aR) by combining large-scale molecular dynamics simulations with Markov state models. These oligomeric structures may also serve as templates for studying oligomerization of other class A GPCRs. Our simulation data revealed that receptor activation results in enhanced oligomerization, more diverse oligomer populations, and a more connected oligomerization network. The active state conformation of the A2aR shifts protein-protein association interfaces to those involving intracellular loop ICL3 and transmembrane helix TM6. Binding of PIP₂ to A2aR stabilizes protein-protein interactions via PIP₂-mediated association interfaces. These results indicate that A2aR oligomerization is responsive to the local membrane lipid environment. This, in turn, suggests a modulatory effect on A2aR whereby a given oligomerization profile favors the dynamic formation of specific supramolecular signaling complexes.

Structure-Based Design and Development of Chemical Probes Targeting Putative MOR-CCR5 Heterodimers to Inhibit Opioid Exacerbated HIV-1 Infectivity

Crystal structures of ligand-bound G-protein-coupled receptors provide tangible templates for rationally designing molecular probes. Herein, we report the structure-based design, chemical synthesis, and biological investigations of bivalent ligands targeting putative mu opioid receptor C-C motif chemokine ligand 5 (MOR-CCR5) heterodimers. The bivalent ligand VZMC013 possessed nanomolar level binding affinities for both the MOR and CCR5, inhibited CCL5-stimulated calcium mobilization, and remarkably improved anti-HIV-1BaL activity over previously reported bivalent ligands. VZMC013 inhibited viral infection in TZM-bl cells coexpressing CCR5 and MOR to a greater degree than cells expressing CCR5 alone. Furthermore, VZMC013 blocked human immunodeficiency virus (HIV)-1 entry in peripheral blood mononuclear cells (PBMC) cells in a concentration-dependent manner and inhibited opioid-accelerated HIV-1 entry more effectively in phytohemagglutinin-stimulated PBMC cells than in the absence of opioids. A three-dimensional molecular model of VZMC013 binding to the MOR-CCR5 heterodimer complex is constructed to elucidate its mechanism of action. VZMC013 is a potent chemical probe targeting MOR-CCR5 heterodimers and may serve as a pharmacological agent to inhibit opioid-exacerbated HIV-1 entry.

Structural dynamics bridge the gap between the genetic and functional levels of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly all physiological processes in the human body and represent an important drug targeting class. The genes encoding the different GPCR (sub)types determine their specific functionality, which can be altered by natural genetic variants and isoforms. Deciphering the molecular link between sequence diversity and its functional consequences is a current challenge and critical for the comprehension of the physiological response of GPCRs. It requires a global understanding of how protein sequence translates into protein structure, how this impacts the structural motions of the protein, and, finally, how all these factors determine the receptor functionality. Here, we discuss available resources and state-of-the-art computational approaches to address this question.

Differences in interactions between transmembrane domains tune the activation of metabotropic glutamate receptors

The metabotropic glutamate receptors (mGluRs) form a family of neuromodulatory G-protein-coupled receptors that contain both a seven-helix transmembrane domain (TMD) and a large extracellular ligand-binding domain (LBD) which enables stable dimerization. Although numerous studies have revealed variability across subtypes in the initial activation steps at the level of LBD dimers, an understanding of inter-TMD interaction and rearrangement remains limited. Here, we use a combination of single molecule fluorescence, molecular dynamics, functional assays, and conformational sensors to reveal that distinct TMD assembly properties drive differences between mGluR subtypes. We uncover a variable region within transmembrane helix 4 (TM4) that contributes to homo- and heterodimerization in a subtype-specific manner and tunes orthosteric, allosteric, and basal activation. We also confirm a critical role for a conserved inter-TM6 interface in stabilizing the active state during orthosteric or allosteric activation. Together this study shows that inter-TMD assembly and dynamic rearrangement drive mGluR function with distinct properties between subtypes.

Single-molecule FRET imaging of GPCR dimers in living cells

Class C G protein-coupled receptors (GPCRs) are known to form stable homodimers or heterodimers critical for function, but the oligomeric status of class A and B receptors, which constitute >90% of all GPCRs, remains hotly debated. Single-molecule fluorescence resonance energy transfer (smFRET) is a powerful approach with the potential to reveal valuable insights into GPCR organization but has rarely been used in living cells to study protein systems. Here, we report generally applicable methods for using smFRET to detect and track transmembrane proteins diffusing within the plasma membrane of mammalian cells. We leverage this in-cell smFRET approach to show agonist-induced structural dynamics within individual metabotropic glutamate receptor dimers. We apply these methods to representative class A, B and C receptors, finding evidence for receptor monomers, density-dependent dimers and constitutive dimers, respectively.

A critical perspective on Markov state model treatments of protein-protein association using coarse-grained simulations

Atomic-level information is essential to explain the specific interactions governing protein–protein recognition in terms of structure and dynamics. Of particular interest is a characterization of the time-dependent kinetic aspects of protein–protein association and dissociation. A powerful framework to characterize the dynamics of complex molecular systems is provided by Markov State Models (MSMs). The central idea is to construct a reduced stochastic model of the full system by defining a set of conformational featured microstates and determining the matrix of transition probabilities between them. While a MSM framework can sometimes be very effective, different combinations of input featurization and simulation methods can significantly affect the robustness and the quality of the information generated from MSMs in the context of protein association. Here, a systematic examination of a variety of MSMs methodologies is undertaken to clarify these issues. To circumvent the uncertainties caused by sampling issues, we use a simplified coarse-grained model of the barnase–barstar protein complex. A sensitivity analysis is proposed to identify the microstates of an MSM that contribute most to the error in conjunction with the transition-based reweighting analysis method for a more efficient and accurate MSM construction.

Allosteric Modulation of GPCRs of Class A by Cholesterol

G-protein coupled receptors (GPCRs) are membrane proteins that convey extracellular signals to the cellular milieu. They represent a target for more than 30% of currently marketed drugs. Here we review the effects of membrane cholesterol on the function of GPCRs of Class A. We review both the specific effects of cholesterol mediated via its direct high-affinity binding to the receptor and non-specific effects mediated by cholesterol-induced changes in the properties of the membrane. Cholesterol binds to many GPCRs at both canonical and non-canonical binding sites. It allosterically affects ligand binding to and activation of GPCRs. Additionally, it changes the oligomerization state of GPCRs. In this review, we consider a perspective of the potential for the development of new therapies that are targeted at manipulating the level of membrane cholesterol or modulating cholesterol binding sites on to GPCRs.

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.

Interclass GPCR heteromerization affects localization and trafficking

Membrane trafficking processes regulate G protein-coupled receptor (GPCR) activity. Although class A GPCRs are capable of activating G proteins in a monomeric form, they can also potentially assemble into functional GPCR heteromers. Here, we showed that the class A serotonin 5-HT2A receptors (5-HT2ARs) affected the localization and trafficking of class C metabotropic glutamate receptor 2 (mGluR2) through a mechanism that required their assembly as heteromers in mammalian cells. In the absence of agonists, 5-HT2AR was primarily localized within intracellular compartments, and coexpression of 5-HT2AR with mGluR2 increased the intracellular distribution of the otherwise plasma membrane-localized mGluR2. Agonists for either 5-HT2AR or mGluR2 differentially affected trafficking through Rab5-positive endosomes in cells expressing each component of the 5-HT2AR-mGluR2 heterocomplex alone, or together. In addition, overnight pharmacological 5-HT2AR blockade with clozapine, but not with M100907, decreased mGluR2 density through a mechanism that involved heteromerization between 5-HT2AR and mGluR2. Using TAT-tagged peptides and chimeric constructs that are unable to form the interclass 5-HT2AR-mGluR2 complex, we demonstrated that heteromerization was necessary for the 5-HT2AR-dependent effects on mGluR2 subcellular distribution. The expression of 5-HT2AR also augmented intracellular localization of mGluR2 in mouse frontal cortex pyramidal neurons. Together, our data suggest that GPCR heteromerization may itself represent a mechanism of receptor trafficking and sorting.

Endogenous Opiates and Behavior: 2018

This paper is the forty-first consecutive installment of the annual anthological review of research concerning the endogenous opioid system, summarizing articles published during 2018 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides and receptors as well as effects of opioid/opiate agonists and antagonists. The review is subdivided into the following specific topics: molecular-biochemical effects and neurochemical localization studies of endogenous opioids and their receptors (2), the roles of these opioid peptides and receptors in pain and analgesia in animals (3) and humans (4), opioid-sensitive and opioid-insensitive effects of nonopioid analgesics (5), opioid peptide and receptor involvement in tolerance and dependence (6), stress and social status (7), learning and memory (8), eating and drinking (9), drug abuse and alcohol (10), sexual activity and hormones, pregnancy, development and endocrinology (11), mental illness and mood (12), seizures and neurologic disorders (13), electrical-related activity and neurophysiology (14), general activity and locomotion (15), gastrointestinal, renal and hepatic functions (16), cardiovascular responses (17), respiration and thermoregulation (18), and immunological responses (19).

Single-molecule analysis reveals agonist-specific dimer formation of µ-opioid receptors

G-protein-coupled receptors (GPCRs) are key signaling proteins that mostly function as monomers, but for several receptors constitutive dimer formation has been described and in some cases is essential for function. Using single-molecule microscopy combined with super-resolution techniques on intact cells, we describe here a dynamic monomer-dimer equilibrium of µ-opioid receptors (µORs), where dimer formation is driven by specific agonists. The agonist DAMGO, but not morphine, induces dimer formation in a process that correlates both temporally and in its agonist- and phosphorylation-dependence with β-arrestin2 binding to the receptors. This dimerization is independent from, but may precede, µOR internalization. These data suggest a new level of GPCR regulation that links dimer formation to specific agonists and their downstream signals.

Selective and Wash-Resistant Fluorescent Dihydrocodeinone Derivatives Allow Single-Molecule Imaging of μ-Opioid Receptor Dimerization

μ-Opioid receptors (μ-ORs) play a critical role in the modulation of pain and mediate the effects of the most powerful analgesic drugs. Despite extensive efforts, it remains insufficiently understood how μ-ORs produce specific effects in living cells. We developed new fluorescent ligands based on the μ-OR antagonist E-p-nitrocinnamoylamino-dihydrocodeinone (CACO), that display high affinity, long residence time and pronounced selectivity. Using these ligands, we achieved single-molecule imaging of μ-ORs on the surface of living cells at physiological expression levels. Our results reveal a high heterogeneity in the diffusion of μ-ORs, with a relevant immobile fraction. Using a pair of fluorescent ligands of different color, we provide evidence that μ-ORs interact with each other to form short-lived homodimers on the plasma membrane. This approach provides a new strategy to investigate μ-OR pharmacology at single-molecule level.

Prediction and targeting of GPCR oligomer interfaces

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.

Molecular Basis of Opioid Action: From Structures to New Leads

Since the isolation of morphine from the opium poppy over 200 years ago, the molecular basis of opioid action has remained the subject of intense inquiry. The identification of specific receptors responsible for opioid function and the discovery of many chemically diverse molecules with unique opioid-like efficacies have provided glimpses into the molecular logic of opioid action. Recent revolutions in the structural biology of transmembrane proteins have, for the first time, yielded high-resolution views into the 3-dimensional shapes of all 4 opioid receptors. These studies have begun to decode the chemical logic that enables opioids to specifically bind and activate their receptor targets. A combination of spectroscopic experiments and computational simulations has provided a view into the molecular movements of the opioid receptors, which itself gives rise to the complex opioid pharmacology observed at the cellular and behavioral levels. Further diversity in opioid receptor structure is driven by both genetic variation and receptor oligomerization. These insights have enabled computational drug discovery efforts, with some evidence of success in the design of completely novel opioids with unique efficacies. The combined progress over the past few years provides hope for new, efficacious opioids devoid of the side effects that have made them the scourge of humanity for millennia.

Quantitative Single-Residue Bioorthogonal Labeling of G Protein-Coupled Receptors in Live Cells

High-end microscopy studies of G protein-coupled receptors (GPCRs) require installing onto the receptors bright and photostable dyes. Labeling must occur in quantitative yields, to allow stoichiometric data analysis, and in a minimally invasive fashion, to avoid perturbing GPCR function. We demonstrate here that the genetic incorporation of trans-cyclooct-2-ene lysine (TCO*) allows achieving quantitative single-residue labeling of the extracellular loops of the β₂-adrenergic and the muscarinic M₂ class A GPCRs, as well as of the corticotropin releasing factor class B GPCR. Labeling occurs within a few minutes by reaction with dye-tetrazine conjugates on the surface of live cells and preserves the functionality of the receptors. To precisely quantify the labeling yields, we devise a method based on fluorescence fluctuation microscopy that extracts the number of labeling sites at the single-cell level. Further, we show that single-residue labeling is better suited for studies of GPCR diffusion than fluorescent-protein tags, since the latter can affect the mobility of the receptor. Finally, by performing dual-color competitive labeling on a single TCO* site, we devise a method to estimate the oligomerization state of a GPCR without the need for a biological monomeric reference, which facilitates the application of fluorescence methods to oligomerization studies. As TCO* and the dye-tetrazines used in this study are commercially available and the described microscopy techniques can be performed on a commercial microscope, we expect our approach to be widely applicable to fluorescence microscopy studies of membrane proteins in general.

Allostery in G protein-coupled receptors investigated by molecular dynamics simulations

G-protein-coupled receptors (GPCRs) are allosteric signaling machines that trigger distinct functional responses depending on the particular conformational state they adopt upon binding. This so-called GPCR functional selectivity is prompted by ligands of different efficacy binding at orthosteric or allosteric sites on the receptor, as well as by interactions with intracellular protein partners or other receptor types. Molecular dynamics (MD) simulations can provide important mechanistic, thermodynamic, and kinetic insights into these interactions at a level of molecular detail that is necessary to rightly inform modern drug discovery. Here, we review the most recent MD contributions to understanding GPCR allostery, with an emphasis on their strengths and limitations.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

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'.

Allosteric Modulation of Class A GPCRs: Targets, Agents, and Emerging Concepts

G-protein-coupled receptors (GPCRs) have been tractable drug targets for decades with over one-third of currently marketed drugs targeting GPCRs. Of these, the class A GPCR superfamily is highly represented, and continued drug discovery for this family of receptors may provide novel therapeutics for a vast range of diseases. GPCR allosteric modulation is an innovative targeting approach that broadens the available small molecule toolbox and is proving to be a viable drug discovery strategy, as evidenced by recent FDA approvals and clinical trials. Numerous class A GPCR allosteric modulators have been discovered recently, and emerging trends such as the availability of GPCR crystal structures, diverse functional assays, and structure-based computational approaches are improving optimization and development. This Perspective provides an update on allosterically targeted class A GPCRs and their disease indications and the medicinal chemistry approaches toward novel allosteric modulators and highlights emerging trends and opportunities in the field.

GPCR homo-oligomerization

G protein-coupled receptors (GPCRs) are an extensive class of trans-plasma membrane proteins that function to regulate a wide range of physiological functions. Despite a general perception that GPCRs exist as monomers an extensive literature has examined whether GPCRs can also form dimers and even higher-order oligomers, and if such organization influences various aspects of GPCR function, including cellular trafficking, ligand binding, G protein coupling and signalling. Here we focus on recent studies that employ approaches ranging from computational methods to single molecule tracking and both quantal brightness and fluorescence fluctuation measurements to assess the organization, stability and potential functional significance of dimers and oligomers within the class A, rhodopsin-like GPCR family.