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

Exploring the binding mechanism of positive allosteric modulators in human metabotropic glutamate receptor 2 using molecular dynamics simulations

Positive allosteric modulators (PAMs) of human metabotropic glutamate receptor 2 (hmGlu₂) are well-known in the treatment of psychiatric disorders for their higher selectivity and lower tolerance risk. A variety of PAMs have been reported over the last decade and two compounds were in Phase II clinical trials for schizophrenia and anxiety. These trials were discontinued on account of the unsatisfactory therapeutic efficacy, but PAMs were explored as novel treatments for addiction and epilepsy. Thus, it is still important to explore novel hmGlu₂ PAMs in the near future. Nowadays, the challenges in optimizing drug potency and improving scaffold diversity for PAMs are the noncomprehensive character analyses of multiple scaffolds; the exploration of the binding modes of PAMs in the allosteric binding site have been proposed to reduce this difficulty. However, there has been no comprehensive research about the binding profiles of PAMs in the hmGlu₂ receptor. To address this issue, this work explores the binding characters of eight PAMs representing five chemical series by multiple computational methods. As a result, the shared binding modes of the eight studied PAMs interacting with 15 residues in the allosteric binding site were defined. In addition, the reduced hydrophobicity with low electronegativity of R₁, increased hydrophobicity with low negative electron density of R₂ and the electronegativity of the linker were identified as indicators that regulate the affinity of PAMs. This finding agrees well with the physicochemical properties of reported multiple series PAMs. This comprehensive work sheds additional light on the binding mechanism and physicochemical regularity underlining PAMs affinity and could be further utilized as a structural and energetic blueprint for discovering and assessing novel PAMs for hmGlu₂.

Allosteric Molecular Switches in Metabotropic Glutamate Receptors

Metabotropic glutamate receptors (mGlu) are class C G protein-coupled receptors of eight subtypes that are omnipresently expressed in the central nervous system. mGlus have relevance in several psychiatric and neurological disorders, therefore they raise considerable interest as drug targets. Allosteric modulators of mGlus offer advantages over orthosteric ligands owing to their increased potential to achieve subtype selectivity, and this has prompted discovery programs that have produced a large number of reported allosteric mGlu ligands. However, the optimization of allosteric ligands into drug candidates has proved to be challenging owing to induced-fit effects, flat or steep structure-activity relationships and unexpected changes in theirpharmacology. Subtle structural changes identified as molecular switches might modulate the functional activity of allosteric ligands. Here we review these switches discovered in the metabotropic glutamate receptor family..

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.

Identification of Potential Dual Negative Allosteric Modulators of Group I mGluR Family: A Shape Based Screening, ADME Prediction, Induced Fit Docking and Molecular Dynamics Approach Against Neurodegenerative Diseases

BACKGROUND

Glutamate is the principal neurotransmitter in the human brain that exerts its effects through ionotropic glutamate receptors (iGluRs) and metabotropic glutamate receptors (mGluRs). The mGluRs are a class of C GPCRs that play a vital role in various neurobiological functions, mGluR1 and mGluR5 are the two receptors distributed throughout the brain involved in cognition, learning, memory, and other important neurological processes. Dysfunction of these receptors can cause neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, X-fragile syndrome, anxiety, depression, etc., hence these receptors are high profile targets for pharmaceutical industries.

OBJECTIVE

The objective of our study is to find the novel dual negative allosteric modulators to regulate both mGluR1 and mGluR5.

METHODS

In this study, shape screening protocol was used to find the dual negative allosteric modulators for both mGluR1 and mGluR5 followed by ADME prediction, induced-fit docking (IFD) and molecular dynamics simulations. Further, DFT analysis and MESP studies were carried out for the selected compounds.

RESULTS

Around 247 compounds were obtained from the eMolecules database and clustered through the CANVAS module and filtered with ADME properties. Furthermore, IFD revealed that the top four compounds (16059796, 25004252, 4667236 and 11670690) having good protein-ligand interactions and binding free energies. The molecular electrostatic potential of the top compounds shows interactions in the amine group and the oxygen atom in the negative potential regions. Finally, molecular dynamics simulations were performed with all the selected as well as the reported compound 29 indicates that the screened hits have better stability of protein ligand complex.

CONCLUSION

Hence, from the results, it is evident that top hits 16059796, 25004252, 4667236 and 11670690 could be a novel and potent dual negative allosteric modulators for mGluR1 and mGluR5.

Recent Advances in the Drug Discovery and Development of Dualsteric/ Bitopic Activators of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) represent the largest family of proteins targeted by drug design and discovery efforts. Of these efforts, the development of GPCR agonists is highly desirable, due to their therapeutic robust utility in treating diseases caused by deficient receptor signaling. One of the challenges in designing potent and selective GPCR agonists lies in the inability to achieve combined high binding affinity and subtype selectivity, due to the high homology between orthosteric sites among GPCR subtypes. To combat this difficulty, researchers have begun to explore the utility of targeting topographically distinct and less conserved binding sites, namely "allosteric" sites. Pursuing these sites offers the benefit of achieving high subtype selectivity, however, it also can result in a decreased binding affinity and potency as compared to orthosteric agonists. Therefore, bitopic ligands comprised of an orthosteric agonist and an allosteric modulator connected by a spacer and allowing binding with both the orthosteric and allosteric sites within one receptor, have been developed. It may combine the high subtype selectivity of an allosteric modulator with the high binding affinity of an orthosteric agonist and provides desired advantages over orthosteric agonists or allosteric modulators alone. Herein, we review the recent advances in the development of bitopic agonists/activators for various GPCR targets and their novel therapeutic potentials.

Fragment-Based Approaches for Allosteric Metabotropic Glutamate Receptor (mGluR) Modulators

Metabotropic glutamate receptors (mGluR) are members of the class C G-Protein Coupled Receptors (GPCR-s) and have eight subtypes. These receptors are responsible for a variety of functions in the central and peripheral nervous systems and their modulation has therapeutic utility in neurological and psychiatric disorders. It was previously established that selective orthosteric modulation of these receptors is challenging, and this stimulated the search for allosteric modulators. Fragment-Based Drug Discovery (FBDD) is a viable approach to find ligands binding at allosteric sites owing to their limited size and interactions. However, it was also observed that the structure-activity relationship of allosteric modulators is often sharp and inconsistent. This can be attributed to the characteristics of the allosteric binding site of mGluRs that is a water channel where ligand binding is accompanied with induced fit and interference with the water network, both playing a role in receptor activation. In this review, we summarize fragment-based drug discovery programs on mGluR allosteric modulators and their contribution identifying of new mGluR ligands with better activity and selectivity.

The role of water and protein flexibility in the structure-based virtual screening of allosteric GPCR modulators: an mGlu5 receptor case study

Stabilizing unique receptor conformations, allosteric modulators of G-protein coupled receptors (GPCRs) might open novel treatment options due to their new pharmacological action, their enhanced specificity and selectivity in both binding and signaling. Ligand binding occurs at intrahelical allosteric sites and involves significant induced fit effects that include conformational changes in the local protein environment and water networks. Based on the analysis of available crystal structures of metabotropic glutamate receptor 5 (mGlu₅) we investigated these effects in the binding of mGlu₅ receptor negative allosteric modulators. A large set of retrospective virtual screens revealed that the use of multiple protein structures and the inclusion of selected water molecules improves virtual screening performance compared to conventional docking strategies. The role of water molecules and protein flexibility in ligand binding can be taken into account efficiently by the proposed docking protocol that provided reasonable enrichment of true positives. This protocol is expected to be useful also for identifying intrahelical allosteric modulators for other GPCR targets.

Computational Drug Design Applied to the Study of Metabotropic Glutamate Receptors

Metabotropic glutamate (mGlu) receptors are a family of eight GPCRs that are attractive drug discovery targets to modulate glutamate action and response. Here we review the application of computational methods to the study of this family of receptors. X-ray structures of the extracellular and 7-transmembrane domains have played an important role to enable structure-based modeling approaches, whilst we also discuss the successful application of ligand-based methods. We summarize the literature and highlight the areas where modeling and experiment have delivered important understanding for mGlu receptor drug discovery. Finally, we offer suggestions of future areas of opportunity for computational work.

Menthol Stereoisomers Exhibit Different Effects on α4β2 nAChR Upregulation and Dopamine Neuron Spontaneous Firing

Menthol contributes to poor cessation rates among smokers, in part because menthol enhances nicotine reward and reinforcement. Mentholated tobacco products contain (-)-menthol and (+)-menthol, in varying proportions. We examined these two menthol stereoisomers for their ability to upregulate α4β2 nAChRs and to alter dopamine neuron firing frequency using long-term, low-dose (≤500 nm) exposure that is pharmacologically relevant to smoking. We found that (-)-menthol upregulates α4β2 nAChRs while (+)-menthol does not. We also found that (-)-menthol decreases dopamine neuron baseline firing and dopamine neuron excitability, while (+)-menthol exhibits no effect. We then examined both stereoisomers for their ability to inhibit α4β2 nAChR function at higher concentrations (>10 µm) using the Xenopus oocyte expression system. To probe for the potential binding site of menthol, we conducted flooding simulations and site-directed mutagenesis. We found that menthol likely binds to the 9´ position on the TM2 (transmembrane M2) helix. We found that menthol inhibition is dependent on the end-to-end distance of the side chain at the 9´ residue. Additionally, we have found that (-)-menthol is only modestly (∼25%) more potent than (+)-menthol at inhibiting wild-type α4β2 nAChRs and a series of L9´ mutant nAChRs. These data reveal that menthol exhibits a stereoselective effect on nAChRs and that the stereochemical effect is much greater for long-term, submicromolar exposure in mice than for acute, higher-level exposure. We hypothesize that of the two menthol stereoisomers, only (-)-menthol plays a role in enhancing nicotine reward through nAChRs on dopamine neurons.

Exploring the Binding Mechanism of Metabotropic Glutamate Receptor 5 Negative Allosteric Modulators in Clinical Trials by Molecular Dynamics Simulations

Metabotropic glutamate receptor 5 (mGlu₅) plays a key role in synaptic information storage and memory, which is a well-known target for a variety of psychiatric and neurodegenerative disorders. In recent years, the increasing efforts have been focused on the design of allosteric modulators, and the negative allosteric modulators (NAMs) are the front-runners. Recently, the architecture of the transmembrane (TM) domain of mGlu₅ receptor has been determined by crystallographic experiment. However, it has been not well understood how the pharmacophores of NAMs accommodated into the allosteric binding site. In this study, molecular dynamics (MD) simulations were performed on mGlu₅ receptor bound with NAMs in preclinical or clinical development to shed light on this issue. In order to identify the key residues, the binding free energies as well as per-residue contributions for NAMs binding to mGlu₅ receptor were calculated. Subsequently, the in silico site-directed mutagenesis of the key residues was performed to verify the accuracy of simulation models. As a result, the shared common features of the studied 5 clinically important NAMs (mavoglurant, dipraglurant, basimglurant, STX107, and fenobam) interacting with 11 residues in allosteric site were obtained. This comprehensive study presented a better understanding of mGlu₅ receptor NAMs binding mechanism, which would be further used as a useful framework to assess and discover novel lead scaffolds for NAMs.

Atom-based 3D-QSAR, induced fit docking, and molecular dynamics simulations study of thieno[2,3-b]pyridines negative allosteric modulators of mGluR5

Atom-based three dimensional-quantitative structure-activity relationship (3D-QSAR) model was developed on the basis of 5-point pharmacophore hypothesis (AARRR) with two hydrogen bond acceptors (A) and three aromatic rings for the derivatives of thieno[2,3-b]pyridine, which modulates the activity to inhibit the mGluR5 receptor. Generation of a highly predictive 3D-QSAR model was performed using the alignment of predicted pharmacophore hypothesis for the training set (R² = 0.84, SD = 0.26, F = 45.8, N = 29) and test set (Q² = 0.74, RMSE = 0.235, Pearson-R = 0.94, N = 9). The best pharmacophore hypothesis AARRR was selected, and developed three dimensional-quantitative structure activity relationship (3D-QSAR) model also supported the outcome of this study by means of favorable and unfavorable electron withdrawing group and hydrophobic regions of most active compound 42d and least active compound 18b. Following, induced fit docking and binding free energy calculations reveals the reliable binding orientation of the compounds. Finally, molecular dynamics simulations for 100 ns were performed to depict the protein-ligand stability. We anticipate that the resulted outcome could be supportive to discover potent negative allosteric modulators for metabotropic glutamate receptor 5 (mGluR5).

Metabotropic glutamate receptors as a new therapeutic target for malignant gliomas

Metabotropic glutamate receptors (mGluR) are predominantly involved in maintenance of cellular homeostasis of central nervous system. However, evidences have suggested other roles of mGluR in human tumors. Aberrant mGluR signaling has been shown to participate in transformation and maintenance of various cancer types, including malignant brain tumors. This review intends to summarize recent findings regarding the involvement of mGluR-mediated intracellular signaling pathways in progression, aggressiveness, and recurrence of malignant gliomas, mainly glioblastomas (GBM), highlighting the potential therapeutic applications of mGluR ligands. In addition to the growing number of studies reporting mGluR gene or protein expression in glioma samples (resections, lineages, and primary cultures), pharmacological blockade in vitro of mGluR1 and mGluR3 by selective ligands has been shown to be anti-proliferative and anti-migratory, decreasing activation of MAPK and PI3K pathways. In addition, mGluR3 antagonists promoted astroglial differentiation of GBM cells and also enabled cytotoxic action of temozolomide (TMZ). mGluR3-dependent TMZ toxicity was supported by increasing levels of MGMT transcripts through an intracellular signaling pathway that sequentially involves PI3K and NF-κB. Further, continuous pharmacological blockade of mGluR1 and mGluR3 have been shown to reduced growth of GBM tumor in two independent in vivo xenograft models. In parallel, low levels of mGluR3 mRNA in GBM resections may be a predictor for long survival rate of patients. Since several Phase I, II and III clinical trials are being performed using group I and II mGluR modulators, there is a strong scientifically-based rationale for testing mGluR antagonists as an adjuvant therapy for malignant brain tumors.

Shining Light on an mGlu5 Photoswitchable NAM: A Theoretical Perspective

Metabotropic glutamate receptors (mGluRs) are important drug targets because of their involvement in several neurological diseases. Among mGluRs, mGlu5 is a particularly high-profile target because its positive or negative allosteric modulation can potentially treat schizophrenia or anxiety and chronic pain, respectively. Here, we computationally and experimentally probe the functional binding of a novel photoswitchable mGlu5 NAM, termed alloswitch-1, which loses its NAM functionality under violet light. We show alloswitch-1 binds deep in the allosteric pocket in a similar fashion to mavoglurant, the co-crystallized NAM in the mGlu5 transmembrane domain crystal structure. Alloswitch-1, like NAM 2-Methyl-6-(phenylethynyl)pyridine (MPEP), is significantly affected by P655M mutation deep in the allosteric pocket, eradicating its functionality. In MD simulations, we show alloswitch-1 and MPEP stabilize the co-crystallized water molecule located at the bottom of the allosteric site that is seemingly characteristic of the inactive receptor state. Furthermore, both NAMs form H-bonds with S809 on helix 7, which may constitute an important stabilizing interaction for NAM-induced mGlu5 inactivation. Alloswitch-1, through isomerization of its amide group from trans to cis is able to form an additional interaction with N747 on helix 5. This may be an important interaction for amide-containing mGlu5 NAMs, helping to stabilize their binding in a potentially unusual cis-amide state. Simulated conformational switching of alloswitch-1 in silico suggests photoisomerization of its azo group from trans to cis may be possible within the allosteric pocket. However, photoexcited alloswitch-1 binds in an unstable fashion, breaking H-bonds with the protein and destabilizing the co-crystallized water molecule. This suggests photoswitching may have destabilizing effects on mGlu5 binding and functionality.

Selective Negative Allosteric Modulation Of Metabotropic Glutamate Receptors – A Structural Perspective of Ligands and Mutants

The metabotropic glutamate receptors have a wide range of modulatory functions in the central nervous system. They are among the most highly pursued drug targets, with relevance for several neurological diseases, and a number of allosteric modulators have entered clinical trials. However, so far this has not led to a marketed drug, largely because of the difficulties in achieving subtype-selective compounds with desired properties. Very recently the first crystal structures were published for the transmembrane domain of two metabotropic glutamate receptors in complex with negative allosteric modulators. In this analysis, we make the first comprehensive structural comparison of all metabotropic glutamate receptors, placing selective negative allosteric modulators and critical mutants into the detailed context of the receptor binding sites. A better understanding of how the different mGlu allosteric modulator binding modes relates to selective pharmacological actions will be very valuable for rational design of safer drugs.

Computational Advances for the Development of Allosteric Modulators and Bitopic Ligands in G Protein-Coupled Receptors

Allosteric modulators of G protein-coupled receptors (GPCRs), which target at allosteric sites, have significant advantages against the corresponding orthosteric compounds including higher selectivity, improved chemical tractability or physicochemical properties, and reduced risk of receptor oversensitization. Bitopic ligands of GPCRs target both orthosteric and allosteric sites. Bitopic ligands can improve binding affinity, enhance subtype selectivity, stabilize receptors, and reduce side effects. Discovering allosteric modulators or bitopic ligands for GPCRs has become an emerging research area, in which the design of allosteric modulators is a key step in the detection of bitopic ligands. Radioligand binding and functional assays ([(35)S]GTPγS and ERK1/2 phosphorylation) are used to test the effects for potential modulators or bitopic ligands. High-throughput screening (HTS) in combination with disulfide trapping and fragment-based screening are used to aid the discovery of the allosteric modulators or bitopic ligands of GPCRs. When used alone, these methods are costly and can often result in too many potential drug targets, including false positives. Alternatively, low-cost and efficient computational approaches are useful in drug discovery of novel allosteric modulators and bitopic ligands to help refine the number of targets and reduce the false-positive rates. This review summarizes the state-of-the-art computational methods for the discovery of modulators and bitopic ligands. The challenges and opportunities for future drug discovery are also discussed.

Allosteric Binding Site and Activation Mechanism of Class C G-Protein Coupled Receptors: Metabotropic Glutamate Receptor Family

Metabotropic glutamate receptors (mGluR) are mainly expressed in the central nervous system (CNS) and contain eight receptor subtypes, named mGluR1 to mGluR8. The crystal structures of mGluR1 and mGluR5 that are bound with the negative allosteric modulator (NAM) were reported recently. These structures provide a basic model for all class C of G-protein coupled receptors (GPCRs) and may aid in the design of new allosteric modulators for the treatment of CNS disorders. However, these structures are only combined with NAMs in the previous reports. The conformations that are bound with positive allosteric modulator (PAM) or agonist of mGluR1/5 remain unknown. Moreover, the structural information of the other six mGluRs and the comparisons of the mGluRs family have not been explored in terms of their binding pockets, the binding modes of different compounds, and important binding residues. With these crystal structures as the starting point, we built 3D structural models for six mGluRs by using homology modeling and molecular dynamics (MD) simulations. We systematically compared their allosteric binding sites/pockets, the important residues, and the selective residues by using a series of comparable dockings with both the NAM and the PAM. Our results show that several residues played important roles for the receptors' selectivity. The observations of detailed interactions between compounds and their correspondent receptors are congruent with the specificity and potency of derivatives or compounds bioassayed in vitro. We then carried out 100 ns MD simulations of mGluR5 (residue 26-832, formed by Venus Flytrap domain, a so-called cysteine-rich domain, and 7 trans-membrane domains) bound with antagonist/NAM and with agonist/PAM. Our results show that both the NAM and the PAM seemed stable in class C GPCRs during the MD. However, the movements of "ionic lock," of trans-membrane domains, and of some activation-related residues in 7 trans-membrane domains of mGluR5 were congruent with the findings in class A GPCRs. Finally, we selected nine representative bound structures to perform 30 ns MD simulations for validating the stabilities of interactions, respectively. All these bound structures kept stable during the MD simulations, indicating that the binding poses in this present work are reasonable. We provided new insight into better understanding of the structural and functional roles of the mGluRs family and facilitated the future structure-based design of novel ligands of mGluRs family with therapeutic potential.

Overlapping binding sites drive allosteric agonism and positive cooperativity in type 4 metabotropic glutamate receptors

Type 4 metabotropic glutamate (mGlu4) receptors are emerging targets for the treatment of various disorders. Accordingly, numerous mGlu4-positive allosteric modulators (PAMs) have been identified, some of which also display agonist activity. To identify the structural bases for their allosteric action, we explored the relationship between the binding pockets of mGlu4 PAMs with different chemical scaffolds and their functional properties. By use of innovative mGlu4 biosensors and second-messenger assays, we show that all PAMs enhance agonist action on the receptor through different degrees of allosteric agonism and positive cooperativity. For example, whereas VU0155041 and VU0415374 display equivalent efficacies [log(τ(B)) = 1.15 ± 0.38 and 1.25 ± 0.44, respectively], they increase the ability of L-AP4 to stabilize the active conformation of the receptor by 4 and 39 times, respectively. Modeling and docking studies identify 2 overlapping binding pockets as follows: a first site homologous to the pocket of natural agonists of class A GPCRs linked to allosteric agonism and a second one pointing toward a site topographically homologous to the Na(+) binding pocket of class A GPCRs, occupied by PAMs exhibiting the strongest cooperativity. These results reveal that intrinsic efficacy and cooperativity of mGlu4 PAMs are correlated with their binding mode, and vice versa, integrating structural and functional knowledge from different GPCR classes.