Modeling the Binding and Function of Metabotropic Glutamate Receptors

Allosteric ligands control the activation of a class C GPCR heterodimer by acting at the transmembrane interface

G protein-coupled receptors (GPCRs) are among the most promising drug targets. They often form homo- and heterodimers with allosteric cross-talk between receptor entities, which contributes to fine-tuning of transmembrane signaling. Specifically controlling the activity of GPCR dimers with ligands is a good approach to clarify their physiological roles and validate them as drug targets. Here, we examined the mode of action of positive allosteric modulators (PAMs) that bind at the interface of the transmembrane domains of the heterodimeric GABA_B receptor. Our site-directed mutagenesis results show that mutations of this interface impact the function of the three PAMs tested. The data support the inference that they act at the active interface between both transmembrane domains, the binding site involving residues of the TM6s of the GABA_B1 and the GABA_B2 subunit. Importantly, the agonist activity of these PAMs involves a key region in the central core of the GABA_B2 transmembrane domain, which also controls the constitutive activity of the GABA_B receptor. This region corresponds to the sodium ion binding site in class A GPCRs that controls the basal state of the receptors. Overall, these data reveal the possibility of developing allosteric compounds able to specifically modulate the activity of GPCR homo- and heterodimers by acting at their transmembrane interface.

Allosteric interactions via the orthosteric ligand binding sites in a constitutive G-protein-coupled receptor homodimer

I interpret some recent data to indicate that co-operative effects take place between the (identical) orthosteric binding sites in a G-protein-coupled receptor dimer. In the current study, the reasonability of this concept was tested by creating a mathematical model. The model is composed of a symmetrical constitutive receptor dimer in which the protomers are able to affect each other allosterically, and it includes binding, receptor activation and signal amplification steps. The model was utilized for analyses of previous data as well as simulations of predicted behaviour. The model demonstrates the behaviour stated in the hypotheses, i.e. even an apparently neutral receptor ligand can allosterically affect agonist binding or receptor activation by binding to the normal orthosteric ligand binding site. Therewith the speculated allosteric action originating from the orthosteric binding site of the dimeric receptor is a realistic possibility. The results of the simulations and curve fitting constitute a reasonable starting point for further studies, and the model can be utilized to design meaningful experiments to investigate these questions.

Exploring the Activation Mechanism of the mGlu5 Transmembrane Domain

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

Reinterpreting anomalous competitive binding experiments within G protein-coupled receptor homodimers using a dimer receptor model

An increasing number of G protein-coupled receptors (GPCRs) have been reported to be expressed in the plasma membrane as dimers. Since most ligand binding data are currently fitted by classical equations developed only for monomeric receptors, the interpretation of data could be misleading in the presence of GPCR dimers. On the other hand, the equations developed from dimer receptor models assuming the existence of two orthosteric binding sites within the dimeric molecule offer the possibility to directly calculate macroscopic equilibrium dissociation constants for the two sites, an index of cooperativity (D_C) that reflects the molecular communication within the dimer and, importantly, a constant of radioligand-competitor allosteric interaction (K_DAB) in competitive assays. Here, we provide a practical way to fit competitive binding data that allows the interpretation of apparently anomalous results, such as competition curves that could be either bell-shaped, monophasic or biphasic depending on the assay conditions. The consideration of a radioligand-competitor allosteric interaction allows fitting these curve patterns both under simulation conditions and in real radioligand binding experiments, obtaining competitor affinity parameters closer to the actual values. Our approach is the first that, assuming the formation of receptor homodimers, is able to explain several experimental results previously considered erroneous due to their impossibility to be fitted. We also deduce the radioligand concentration responsible for the conversion of biphasic to monophasic or to bell-shaped curves in competitive radioligand binding assays. In conclusion, bell-shaped curves in competitive binding experiments constitute evidence for GPCR homodimerization.

Quantifying the allosteric interactions within a G-protein-coupled receptor heterodimer

G-protein-coupled receptors are central to signal transduction and cell communication. The possibility that cells use receptor heteromerization to modulate individual receptor pathways is a surmise that cannot be precluded. Given the complexity of these processes, mathematical models contribute to understanding how receptors and their respective ligands regulate signaling. Here, a mathematical model is presented that quantifies the allosteric interactions within a receptor heterodimer. The model is based on the operational model of allosterism including constitutive receptor activity, which provides the pharmacological analysis of heteromerization with well-established and widely used modeling and fitting procedures.

Pharmacological evidence for a metabotropic glutamate receptor heterodimer in neuronal cells

Metabotropic glutamate receptors (mGluRs) are mandatory dimers playing important roles in regulating CNS function. Although assumed to form exclusive homodimers, 16 possible heterodimeric mGluRs have been proposed but their existence in native cells remains elusive. Here, we set up two assays to specifically identify the pharmacological properties of rat mGlu heterodimers composed of mGlu2 and 4 subunits. We used either a heterodimer-specific conformational LRET-based biosensor or a system that guarantees the cell surface targeting of the heterodimer only. We identified mGlu2-4 specific pharmacological fingerprints that were also observed in a neuronal cell line and in lateral perforant path terminals naturally expressing mGlu2 and mGlu4. These results bring strong evidence for the existence of mGlu2-4 heterodimers in native cells. In addition to reporting a general approach to characterize heterodimeric mGluRs, our study opens new avenues to understanding the pathophysiological roles of mGlu heterodimers.

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

Metabotropic glutamate receptor subtype 5: molecular pharmacology, allosteric modulation and stimulus bias

The metabotropic glutamate receptor subtype 5 (mGlu5 ) is a family C GPCR that has been implicated in various neuronal processes and, consequently, in several CNS disorders. Over the past few decades, GPCR-based drug discovery, including that for mGlu5 receptors, has turned considerable attention to targeting allosteric binding sites. Modulation of endogenous agonists by allosteric ligands offers the advantages of spatial and temporal fine-tuning of receptor activity, increased selectivity and reduced adverse effects with the potential to elicit improved clinical outcomes. Further, with greater appreciation of the multifaceted nature of the transduction of mGlu5 receptor signalling, it is increasingly apparent that drug discovery must take into consideration unique receptor conformations and the potential for stimulus-bias. This novel paradigm proposes that different ligands may differentially modulate distinct signalling pathways arising from the same receptor. We review our current understanding of the complexities of mGlu5 receptor signalling and regulation, and how these relate to allosteric ligands. Ultimately, a deeper appreciation of these relationships will provide the foundation for targeted drug design of compounds with increased selectivity, not only for the desired receptor but also for the desired signalling outcome from the receptor. Linked Articles This article is part of a themed section on Molecular Pharmacology of G Protein-Coupled Receptors. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v173.20/issuetoc.

Computational analysis of negative and positive allosteric modulator binding and function in metabotropic glutamate receptor 5 (in)activation

Metabotropic glutamate receptors (mGluRs) are high-profile G-protein coupled receptors drug targets because of their involvement in several neurological disease states, and mGluR5 in particular is a subtype whose controlled allosteric modulation, both positive and negative, can potentially be useful for the treatment of schizophrenia and relief of chronic pain, respectively. Here we model mGluR5 with a collection of positive and negative allosteric modulators (PAMs and NAMs) in both active and inactive receptor states, in a manner that is consistent with experimental information, using a specialized protocol that includes homology to increase docking accuracy, and receptor relaxation to generate an individual induced fit with each allosteric modulator. Results implicate two residues in particular for NAM and PAM function: NAM interaction with W785 for receptor inactivation, and NAM/PAM H-bonding with S809 for receptor (in)activation. Models suggest the orientation of the H-bond between allosteric modulator and S809, controlled by PAM/NAM chemistry, influences the position of TM7, which in turn influences the shape of the allosteric site, and potentially the receptor state. NAM-bound and PAM-bound mGluR5 models also reveal that although PAMs and NAMs bind in the same pocket and share similar binding modes, they have distinct effects on the conformation of the receptor. Our models, together with the identification of a possible activation mechanism, may be useful in the rational design of new allosteric modulators for mGluR5.

Collective behaviours: from biochemical kinetics to electronic circuits

In this work we aim to highlight a close analogy between cooperative behaviors in chemical kinetics and cybernetics; this is realized by using a common language for their description, that is mean-field statistical mechanics. First, we perform a one-to-one mapping between paradigmatic behaviors in chemical kinetics (i.e., non-cooperative, cooperative, ultra-sensitive, anti-cooperative) and in mean-field statistical mechanics (i.e., paramagnetic, high and low temperature ferromagnetic, anti-ferromagnetic). Interestingly, the statistical mechanics approach allows a unified, broad theory for all scenarios and, in particular, Michaelis-Menten, Hill and Adair equations are consistently recovered. This framework is then tested against experimental biological data with an overall excellent agreement. One step forward, we consistently read the whole mapping from a cybernetic perspective, highlighting deep structural analogies between the above-mentioned kinetics and fundamental bricks in electronics (i.e. operational amplifiers, flashes, flip-flops), so to build a clear bridge linking biochemical kinetics and cybernetics.

Allosteric transition: a comparison of two models

Introduction

Two recent models are in use for analysis of allosteric drug action at receptor sites remote from orthosteric binding sites. One is an allosteric two-state mechanical model derived in 2000 by David Hall. The other is an extended operational model developed in 2007 by Arthur Christopoulos’s group. The models are valid in pharmacology, enzymology, transportology as well as several other fields of biology involving allosteric concentration effects.

Results

I show here that Hall’s model for interactions between an orthoster, an alloster, and a receptive unit is the best choice of model both for simulation and analysis of allosteric concentration-responses at equilibrium or steady-state.

Conclusions

As detailed knowledge of receptors systems becomes available, systems with several pathways and states and/ or more than two binding sites should be analysed by extended forms of the Hall model rather than for instance a Hill type exponentiation of terms as introduced in non-mechanistic (operational) model approaches; yielding semi-quantitative estimates of actual system parameters based on Hill’s unlikely simultaneity model for G protein-coupled receptors.

Modeling cooperativity effects in dimeric G protein-coupled receptors

G protein-coupled receptors organize into oligomeric arrangements to exert their function. In this chapter, three models of dimeric receptors, the two-state dimer receptor model, the metabotropic glutamate receptor model, and the asymmetric/symmetric three-state dimer receptor model are revisited focusing on the cooperative effects between their binding sites and the subunits they are composed of. The mathematical analysis reveals the complexity of the intra-receptor interactions providing insights on the mechanistic aspects of receptor function.

Group III and subtype 4 metabotropic glutamate receptor agonists: discovery and pathophysiological applications in Parkinson's disease

Restoring the balance between excitatory and inhibitory circuits in the basal ganglia, following the loss of dopaminergic (DA) neurons of the substantia nigra pars compacta, represents a major challenge to treat patients affected by Parkinson's disease (PD). The imbalanced situation in favor of excitation in the disease state may also accelerate excitotoxic processes, thereby representing a potential target for neuroprotective therapies. Reducing the excitatory action of glutamate, the major excitatory neurotransmitter in the basal ganglia, should lead to symptomatic improvement for PD patients and may promote the survival of DA neurons. Recent studies have focused on the modulatory action of metabotropic glutamate (mGlu) receptors on neurodegenerative diseases including PD. Group III mGlu receptors, including subtypes 4, 7 and 8, are largely expressed in the basal ganglia. Recent studies highlight the use of selective mGlu4 receptor positive allosteric modulators (PAMs) for the treatment of PD. Here we review the effects of newly-designed group-III orthosteric agonists on neuroprotection, neurorestoration and reduction of l-DOPA induced dyskinesia in animal models of PD. The combination of orthosteric mGlu4 receptor selective agonists with PAMs may open new avenues for the symptomatic treatment of PD. This article is part of a Special Issue entitled 'Metabotropic Glutamate Receptors'.

Possible new targets for GPCR modulation: allosteric interactions, plasma membrane domains, intercellular transfer and epigenetic mechanisms

It has been estimated that at least 50% of the drugs available on the market act on G-protein coupled receptors (GPCRs) and most of these are basically or agonists or antagonists of this type of receptors. Herein, we propose new putative targets for drug development based on recent data on GPCR allosterism and on the existence of receptor mosaics (RMs). The main target for drug development is still GPCRs, but the focus is not the orthosteric binding pocket. According to the mosaic model of the plasma membrane, we mainly discuss the possibility of indirect modulatory pharmacological actions on expression/function of GPCRs. In particular, the following two new targets will be analyzed: a) The possibility of pharmacological interventions on the roamer-type of volume transmission (VT), which allow the intercellular transfer of set of signal molecules such as GPCRs, tetraspanins and ribonucleic acids. Thus, there is the possibility of pharmacological interventions on the decoding capabilities of neurons and/or glial cells by means of an action on composition and release of micro-vesicles. b) The possibility of pharmacological interventions on epigenetic mechanisms by taking into account their inter-relationships with GPCRs. As a matter of fact, there are epigenetic changes that are characteristic of periods of developmental plasticity that could provide a target for therapeutic intervention in the event of brain damage. We believe that almost all the biochemical knowledge presently available on GPCRs can be used in the development of these new pharmacological approaches.

The promiscuous mGlu5 receptor--a range of partners for therapeutic possibilities?

The issue of non-specific effects for potential therapeutics is particularly salient in neurological/psychiatric disorders, where adverse drug reactions could impair critical brain functions. The issue of specificity is not limited to candidate molecules, as receptor targets themselves often influence physiological as well as pathological outcomes. Metabotropic glutamate receptor 5 (mGlu5) is an example of a "promiscuous" receptor target that has been implicated in addiction, but also many other processes. However, if receptor modulation could be restricted to specific pathways/brain regions, mGlu5 may still prove to be a viable therapeutic target for various indications. Using this premise, a number of possible methods to refine drug development strategy are discussed, including exploiting specific interactions of mGlu5 with other receptors to narrow the influence of pharmacological agents, and also the use of RNA interference targeted to specific cells/regions of the brain.

Modelling the interdependence between the stoichiometry of receptor oligomerization and ligand binding for a coexisting dimer/tetramer receptor system

Many G protein-coupled receptors have been shown to exist as oligomers, but the oligomerization state and the effects of this on receptor function are unclear. For some G protein-coupled receptors, in ligand binding assays, different radioligands provide different maximal binding capacities. Here we have developed mathematical models for co-expressed dimeric and tetrameric species of receptors. We have considered models where the dimers and tetramers are in equilibrium and where they do not interconvert and we have also considered the potential influence of the ligands on the degree of oligomerization. By analogy with agonist efficacy, we have considered ligands that promote, inhibit or have no effect on oligomerization. Cell surface receptor expression and the intrinsic capacity of receptors to oligomerize are quantitative parameters of the equations. The models can account for differences in the maximal binding capacities of radioligands in different preparations of receptors and provide a conceptual framework for simulation and data fitting in complex oligomeric receptor situations.

On the fitting of binding data when receptor dimerization is suspected

Mechanistic and empirical modelling are compared in context of dimeric receptors. In particular, the supposed advantages of the two-state dimer model for fitting of binding data with respect to classical approaches such as the two-independent sites model are investigated. The two models are revisited from both the mechanistic and empirical point of views. The problem of overparameterized models and the benefits of the concurrent use of mechanistic and empirical models for mechanism analysis are discussed. The pros and cons of mathematical models are examined with special emphasis given to the interpretation of the connection between the shapes of the curves and receptor cooperativity. It is shown that a given pharmacological phenotype (curve shape) can be obtained from different receptor genotypes (as, for instance, non-interconvertible monomeric receptor species, receptor-G protein interactions and dimeric receptors), though values of the Hill coefficient greater than one are indicative of receptor oligomerization. The existence of a relationship between the recently defined dimer cooperativity index and the more familiar Hill coefficient is proven.