Positive allosteric modulators of type 5 metabotropic glutamate receptors (mGluR5) and their therapeutic potential for the treatment of CNS disorders

Studies utilizing selective pharmacological antagonists or targeted gene deletion have demonstrated that type 5 metabotropic glutamate receptors (mGluR5) are critical mediators and potential therapeutic targets for the treatment of numerous disorders of the central nervous system (CNS), including depression, anxiety, drug addiction, chronic pain, Fragile X syndrome, Parkinson's disease, and gastroesophageal reflux disease. However, in recent years, the development of positive allosteric modulators (PAMs) of the mGluR5 receptor have revealed that allosteric activation of this receptor may also be of potential therapeutic benefit for the treatment of other CNS disorders, including schizophrenia, cognitive deficits associated with chronic drug use, and deficits in extinction learning. Here we summarize the discovery and characterization of various mGluR5 PAMs, with an emphasis on those that are systemically active. We will also review animal studies showing that these molecules have potential efficacy as novel antipsychotic agents. Finally, we will summarize findings that suggest that mGluR5 PAMs have pro-cognitive effects such as the ability to enhance synaptic plasticity, improve performance in various learning and memory tasks, including extinction of drug-seeking behavior, and reverse cognitive deficits produced by chronic drug use.     

Two-track virtual screening approach to identify both competitive and allosteric inhibitors of human small C-terminal domain phosphatase 1

Despite a wealth of persuasive evidence for the involvement of human small C-terminal domain phosphatase 1 (Scp1) in the impairment of neuronal differentiation and in Huntington’s disease, small-molecule inhibitors of Scp1 have been rarely reported so far. This study aims to the discovery of both competitive and allosteric Scp1 inhibitors through the two-track virtual screening procedure. By virtue of the improvement of the scoring function by implementing a new molecular solvation energy term and by reoptimizing the atomic charges for the active-site Mg²⁺ ion cluster, we have been able to identify three allosteric and five competitive Scp1 inhibitors with low-micromolar inhibitory activity. Consistent with the results of kinetic studies on the inhibitory mechanisms, the allosteric inhibitors appear to be accommodated in the peripheral binding pocket through the hydrophobic interactions with the nonpolar residues whereas the competitive ones bind tightly in the active site with a direct coordination to the central Mg²⁺ ion. Some structural modifications to improve the biochemical potency of the newly identified inhibitors are proposed based on the binding modes estimated with docking simulations.     

In Silico Prediction and Validation of CB2 Allosteric Binding Sites to Aid the Design of Allosteric Modulators

Although the 3D structures of active and inactive cannabinoid receptors type 2 (CB2) are available, neither the X-ray crystal nor the cryo-EM structure of CB2-orthosteric ligand-modulator has been resolved, prohibiting the drug discovery and development of CB2 allosteric modulators (AMs). In the present work, we mainly focused on investigating the potential allosteric binding site(s) of CB2. We applied different algorithms or tools to predict the potential allosteric binding sites of CB2 with the existing agonists. Seven potential allosteric sites can be observed for either CB2-CP55940 or CB2-WIN 55,212-2 complex, among which sites B, C, G and K are supported by the reported 3D structures of Class A GPCRs coupled with AMs. Applying our novel algorithm toolset-MCCS, we docked three known AMs of CB2 including Ec2la (C-2), trans-β-caryophyllene (TBC) and cannabidiol (CBD) to each site for further comparisons and quantified the potential binding residues in each allosteric binding site. Sequentially, we selected the most promising binding pose of C-2 in five allosteric sites to conduct the molecular dynamics (MD) simulations. Based on the results of docking studies and MD simulations, we suggest that site H is the most promising allosteric binding site. We plan to conduct bio-assay validations in the future.     

Discovery of cryptic allosteric sites using reversed allosteric communication by a combined computational and experimental strategy

Allostery, which is one of the most direct and efficient methods to fine-tune protein functions, has gained increasing recognition in drug discovery. However, there are several challenges associated with the identification of allosteric sites, which is the fundamental cornerstone of drug design. Previous studies on allosteric site predictions have focused on communication signals propagating from the allosteric sites to the orthosteric sites. However, recent biochemical studies have revealed that allosteric coupling is bidirectional and that orthosteric perturbations can modulate allosteric sites through reversed allosteric communication. Here, we proposed a new framework for the prediction of allosteric sites based on reversed allosteric communication using a combination of computational and experimental strategies (molecular dynamics simulations, Markov state models, and site-directed mutagenesis). The desirable performance of our approach was demonstrated by predicting the known allosteric site of the small molecule MDL-801 in nicotinamide dinucleotide (NAD⁺)-dependent protein lysine deacetylase sirtuin 6 (Sirt6). A potential novel cryptic allosteric site located around the L116, R119, and S120 residues within the dynamic ensemble of Sirt6 was identified. The allosteric effect of the predicted site was further quantified and validated using both computational and experimental approaches. This study proposed a state-of-the-art computational pipeline for detecting allosteric sites based on reversed allosteric communication. This method enabled the identification of a previously uncharacterized potential cryptic allosteric site on Sirt6, which provides a starting point for allosteric drug design that can aid the identification of candidate pockets in other therapeutic targets.

Using reversed allosteric communication, we performed MD simulations, MSMs, and mutagenesis experiments, to discover allosteric sites. It reproduced the known allosteric site for MDL-801 on Sirt6 and uncovered a novel cryptic allosteric Pocket X.     

Computational Study of SARS-CoV-2 RNA Dependent RNA Polymerase Allosteric Site Inhibition

The COVID-19 pandemic has caused millions of fatalities since 2019. Despite the availability of vaccines for this disease, new strains are causing rapid ailment and are a continuous threat to vaccine efficacy. Here, molecular docking and simulations identify strong inhibitors of the allosteric site of the SARS-CoV-2 virus RNA dependent RNA polymerase (RdRp). More than one hundred different flavonoids were docked with the SARS-CoV-2 RdRp allosteric site through computational screening. The three top hits were Naringoside, Myricetin and Aureusidin 4,6-diglucoside. Simulation analyses confirmed that they are in constant contact during the simulation time course and have strong association with the enzyme’s allosteric site. Absorption, distribution, metabolism, excretion and toxicity (ADMET) data provided medicinal information of these top three hits. They had good human intestinal absorption (HIA) concentrations and were non-toxic. Due to high mutation rates in the active sites of the viral enzyme, these new allosteric site inhibitors offer opportunities to drug SARS-CoV-2 RdRp. These results provide new information for the design of novel allosteric inhibitors against SARS-CoV-2 RdRp.     

Design,Synthesis and Structural Analysis of Glucocerebrosidase Imaging Agents

Gaucher disease (GD) is a lysosomal storage disorder caused by inherited deficiencies in β-glucocerebrosidase (GBA). Current treatments require rapid disease diagnosis and a means of monitoring therapeutic efficacy, both of which may be supported by the use of GBA-targeting activity-based probes (ABPs). Here, we report the synthesis and structural analysis of a range of cyclophellitol epoxide and aziridine inhibitors and ABPs for GBA. We demonstrate their covalent mechanism-based mode of action and uncover binding of the new N-functionalised aziridines to the ligand binding cleft. These inhibitors became scaffolds for the development of ABPs; the O6-fluorescent tags of which bind in an allosteric site at the dimer interface. Considering GBA's preference for O6- and N-functionalised reagents, a bi-functional aziridine ABP was synthesized as a potentially more powerful imaging agent. Whilst this ABP binds to two unique active site clefts of GBA, no further benefit in potency was achieved over our first generation ABPs. Nevertheless, such ABPs should serve useful in the study of GBA in relation to GD and inform the design of future probes.     

Molecular Dynamics Approach to Probe the Allosteric Inhibition of PTP1B by Chlorogenic and Cichoric Acid

Protein tyrosine phosphatase 1B (PTP1B), a major negative regulator of the insulin and leptin signaling pathway, is a potential target for therapeutic intervention against diabetes and obesity. The recent discovery of an allosteric site in PTP1B has created an alternate strategy in the development of PTP1B targeted therapy. The current study investigates the molecular interactions between the allosteric site of PTP1B with two caffeoyl derivatives, chlorogenic acid (CGA) and cichoric acid (CHA), using computational strategies. Molecular docking analysis with CGA and CHA at the allosteric site of PTP1B were performed and the resulting protein-ligand complexes used for molecular dynamics simulation studies for a time scale of 10 ns. Results show stable binding of CGA and CHA at the allosteric site of PTP1B. The flexibility of the WPD loop was observed to be constrained by CGA and CHA in the open (inactive), providing molecular mechanism of allosteric inhibition. The allosteric inhibition of CGA and CHA of PTP1B was shown to be favorable due to no restriction by the α-7 helix in the binding of CGA and CHA at the allosteric binding site. In conclusion, our results exhibit an inhibitory pattern of CGA and CHA against PTP1B through potent binding at the allosteric site.     

Allostery Inhibition of BACE1 by Psychotic and Meroterpenoid Drugs in Alzheimer’s Disease Therapy

In over a century since its discovery, Alzheimer’s disease (AD) has continued to be a global health concern due to its incurable nature and overwhelming increase among older people. In this paper, we give an overview of the efforts of researchers towards identifying potent BACE1 exosite-binding antibodies and allosteric inhibitors. Herein, we apply computer-aided drug design (CADD) methods to unravel the interactions of some proposed psychotic and meroterpenoid BACE1 allosteric site inhibitors. This study is aimed at validating the allosteric potentials of these selected compounds targeted at BACE1 inhibition. Molecular docking, molecular dynamic (MD) simulations, and post-MD analyses are carried out on these selected compounds, which have been experimentally proven to exhibit allosteric inhibition on BACE1. The SwissDock software enabled us to identify more than five druggable pockets on the BACE1 structural surface using docking. Besides the active site region, a melatonin derivative (compound 1) previously proposed as a BACE1 allostery inhibitor showed appreciable stability at eight different subsites on BACE1. Refinement with molecular dynamic (MD) simulations shows that the identified non-catalytic sites are potential allostery sites for compound 1. The allostery and binding mechanism of the selected potent inhibitors show that the smaller the molecule, the easier the attachment to several enzyme regions. This finding hereby establishes that most of these selected compounds failed to exhibit strong allosteric binding with BACE1 except for compound 1. We hereby suggest that further studies and additional identification/validation of other BACE1 allosteric compounds be done. Furthermore, this additional allosteric site investigation will help in reducing the associated challenges with designing BACE1 inhibitors while exploring the opportunities in the design of allosteric BACE1 inhibitors.     

Designing Molecular Spanners to Throw in the Protein Networks

Proteins govern most aspects of cellular life and, through specific interfaces, are typically involved in intricate protein–protein interaction (PPI) networks and signaling pathways. Subtle up- or downregulation of key protein functions and PPIs results in disease; still, the preferred option to contrast the role of a protein in disease and healthy conditions alike remains its outright shutdown through orthosteric ligands that block its active site. Here, we explore subtler alternatives to modulate proteins and PPIs. Driven by a view of proteins as dynamic entities, we discuss ways to identify allosteric binding sites, which, when targeted by tailored ligands, can induce significant changes in the active site of a protein, and lead to agonistic or antagonistic effects. We also summarize the selective regulation of specific PPIs—either direct or allosteric—and show that effects can be stabilizing as well as destabilizing, depending on how the conformational equilibrium of a protein is shifted.     

Syntheses of mGluR5 PET radioligands through the radiofluorination of diaryliodonium tosylates

3-Fluoro-1-((thiazol-4-yl)ethynyl)benzenes constitute an important class of high-affinity metabotropic glutamate subtype 5 receptor (mGluR5) ligands, some of which have been labeled with fluorine-18 (t(1/2) = 109.7 min), to provide radioligands for molecular imaging of brain mGluR5 in living animal and human subjects with positron emission tomography (PET). Labeling in the 3-fluoro position of such ligands can be achieved through aromatic nucleophilic substitution of a halide leaving group with [(18)F]fluoride ion when a weakly activating m-nitrile group is present, but is generally very low yielding (<8%). Here we used a microfluidic reaction platform to show that greatly enhanced (up to 6-fold) radiochemical yields can be achieved from suitably synthesized diaryliodonium tosylate precursors. The presence of a m-nitrile or other activating group is not required. Similar conditions were adopted in a more conventional automated radiochemistry platform having a single-pot reactor, to produce mGluR5 radioligands with useful radioactivities for PET imaging.     

Metal ion induced allosteric transition in the catalytic activity of an artificial phosphodiesterase

An artificial phosphodiesterase () bearing two types of metal binding sites, a catalytic site and a regulatory bipyridine site showed a unique allosteric transition in the catalytic activity against the metal concentration. The rate constants for the hydrolysis reaction of 2-hydroxypropyl-p-nitrophenyl phosphate (HPNP) and RNA dimer (ApA) with and without an effector metal ion were evaluated; the k(obs) value of HPNP hydrolysis for .(Zn(2+))(3) (2.0 x 10(-4) s(-1)) is 3.3 times larger than that for .(Zn(2+))(2). In the case of and Cu(2+), a 19.4 times larger k(obs) value was obtained for .(Cu(2+))(3) (1.2 x 10(-3) s(-1)) against .(Cu(2+))(2). The increase in the catalytic activity is ascribed to the allosteric conformational transition of induced by the coordination of effector metal ion to the Bpy moiety. A detailed investigation revealed that a conformational change of induced by the third M(2+) complexation enhances the rate of hydrolysis rather than a change in the substrate affinity.     

Understanding cAMP-dependent allostery by NMR spectroscopy: comparative analysis of the EPAC1 cAMP-binding domain in its apo and cAMP-bound states

cAMP (adenosine 3',5'-cyclic monophosphate) is a ubiquitous second messenger that activates a multitude of essential cellular responses. Two key receptors for cAMP in eukaryotes are protein kinase A (PKA) and the exchange protein directly activated by cAMP (EPAC), which is a recently discovered guanine nucleotide exchange factor (GEF) for the small GTPases Rap1 and Rap2. Previous attempts to investigate the mechanism of allosteric activation of eukaryotic cAMP-binding domains (CBDs) at atomic or residue resolution have been hampered by the instability of the apo form, which requires the use of mixed apo/holo systems, that have provided only a partial picture of the CBD apo state and of the allosteric networks controlled by cAMP. Here, we show that, unlike other eukaryotic CBDs, both apo and cAMP-bound states of the EPAC1 CBD are stable under our experimental conditions, providing a unique opportunity to define at an unprecedented level of detail the allosteric interactions linking two critical functional sites of this CBD. These are the phosphate binding cassette (PBC), where cAMP binds, and the N-terminal helical bundle (NTHB), which is the site of the inhibitory interactions between the regulatory and catalytic regions of EPAC. Specifically, the combined analysis of the cAMP-dependent changes in chemical shifts, 2 degrees structure probabilities, hydrogen/hydrogen exchange (H/H) and hydrogen/deuterium exchange (H/D) protection factors reveals that the long-range communication between the PBC and the NTHB is implemented by two distinct intramolecular cAMP-signaling pathways, respectively, mediated by the beta2-beta3 loop and the alpha6 helix. Docking of cAMP into the PBC perturbs the NTHB inner core packing and the helical probabilities of selected NTHB residues. The proposed model is consistent with the allosteric role previously hypothesized for L273 and F300 based on site-directed mutagenesis; however, our data show that such a contact is part of a significantly more extended allosteric network that, unlike PKA, involves a tight coupling between the alpha- and beta-subdomains of the EPAC CBD. The proposed mechanism of allosteric activation will serve as a basis to understand agonism and antagonism in the EPAC system and provides also a general paradigm for how small ligands control protein-protein interfaces.     

Structure-activity analysis of the purine binding site of human liver glycogen phosphorylase

Human liver glycogen phosphorylase (HLGP) catalyzes the breakdown of glycogen to maintain serum glucose levels and is a therapeutic target for diabetes. HLGP is regulated by multiple interacting allosteric sites, each of which is a potential drug binding site. We used surface plasmon resonance (SPR) to screen for compounds that bind to the purine allosteric inhibitor site. We determined the affinities of a series of compounds and solved the crystal structures of three representative ligands with K(D) values from 17-550 microM. The crystal structures reveal that the affinities are partly determined by ligand-specific water-mediated hydrogen bonds and side chain movements. These effects could not be predicted; both crystallographic and SPR studies were required to understand the important features of binding and together provide a basis for the design of new allosteric inhibitors targeting this site.     

Mutant-Selective Allosteric EGFR Degraders are Effective Against a Broad Range of Drug-Resistant Mutations

Targeting epidermal growth factor receptor (EGFR) through an allosteric mechanism provides a potential therapeutic strategy to overcome drug-resistant EGFR mutations that emerge within the ATP binding site. Here, we develop an allosteric EGFR degrader, DDC-01-163, which can selectively inhibit the proliferation of L858R/T790M (L/T) mutant Ba/F3 cells while leaving wildtype EGFR Ba/F3 cells unaffected. DDC-01-163 is also effective against osimertinib-resistant cells with L/T/C797S and L/T/L718Q EGFR mutations. When combined with an ATP-site EGFR inhibitor, osimertinib, the anti-proliferative activity of DDC-01-163 against L858R/T790M EGFR-Ba/F3 cells is enhanced. Collectively, DDC-01-163 is a promising allosteric EGFR degrader with selective activity against various clinically relevant EGFR mutants as a single agent and when combined with an ATP-site inhibitor. Our data suggests that targeted protein degradation is a promising drug development approach for mutant EGFR.     

Parameterization of Org27569: An allosteric modulator of the cannabinoid CB1 G protein‐coupled receptor

The cannabinoid CB₁ receptor is a class A G protein‐coupled receptor (GPCR) that is the most widely expressed GPCR in the brain. Many GPCRs contain allosteric binding sites for endogenous and/or synthetic ligands, which are topographically distinct from the agonist‐binding site that is known as the orthosteric site. While both endogenous and synthetic ligands that act at the CB₁ orthosteric site have been known for some time, compounds that act at a CB₁ allosteric site have only recently been discovered. The most studied of these is 5‐chloro‐3‐ethyl‐1H‐indole‐2‐carboxylic acid [2‐(4‐piperidin‐1‐ylphenyl)ethyl]amide (Org27569). Because allosteric ligands are thought to act through conformational changes in the receptor that are transmitted from the allosteric to the orthosteric site, computational studies of the structural and dynamic interactions of Org27569 with the CB₁ receptor are crucial to achieve a molecular level understanding of the basis of action of this important new class of compounds. To date, such computational studies have not been possible due to the lack of a complete set of molecular mechanics force field parameters for Org27569. Here, we present the development of missing CHARMM force field parameters for Org27569 using previously published methods and the validation and application of these new parameters using normal mode analysis and molecular dynamics simulations combined with experimental infrared measurements. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Synthesis and evaluation of 1,2,4-methyltriazines as mGluR5 antagonists

In previous studies we showed that 3-(substituted phenylethynyl)-5-methyl[1,2,4]triazine analogues of MPEP were potent antagonists of glutamate-mediated mobilization of internal calcium in an mGluR5 in vitro efficacy assay. In the present study we report the synthesis and evaluation of six 3-(substituted biphenylethynyl)-5-methyl[1,2,4]triazines (5a-f), and five 3-(substituted phenoxyphenylethynyl)-5-methyltriazines (6a-e). Compound 2-(4-fluorophenyl-5-[2-(5-methyl[1,2,4]triazine-3-yl)ethynyl]benzonitrile (5f) with an IC(50) of 28.2 nM was the most potent analogue.     

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.

     

Preparation of a Series of 5-Methyl-3-(substituted)-[1,2,4]triazines

Procedures for the synthesis of thirty-six 5-methyl-3-(substituted)-[1,2,4]triazines have been described. These compounds were evaluated for antagonism at metabotropic glutamate receptor subtype 5 (mGluR5). Two compounds, 5b and 3c, were determined to be low micromolar inhibitors of mGluR5.     

Computational systematic selectivity of the Fasalog inhibitors between ROCK-I and ROCK-II kinase isoforms in Alzheimer’s disease

Human Rho-associated coiled-coil forming kinase (ROCK) is a class of essential neurokinases that consists of two structurally conserved isoforms ROCK-I and ROCK-II; they have been revealed to play distinct roles in the pathogenesis of Alzheimer’s disease (AD) and other neurological disorders. Selective targeting of the two kinase isoforms with small-molecule inhibitors is a great challenge due to the surprisingly high homology in kinase domain (92 %) and the full identity in kinase active site (100 %). Here, we describe a computational protocol to systematically profile the selectivity of Fasudil and its 25 analogs (termed as Fasalogs) between the two kinase isoforms. It is suggested that the substitution of Fasudil’s 1,4-diazepane moiety with rigid ring such as Ripasudil and Dimehtylfasudil would render the resulting inhibitors of ROCK-II over ROCK-I (II-o-I) selectivity, while the substitution with long, flexible group such as H-89 and BDBM92607 tends to have I-o-II selectivity. Structural analysis reveals that the inhibitor affinity is not only determined by the identical active site, but also contributed from the non-identical first and second shells of the site as well as other non-conserved kinase regions, which can indirectly influence the active site and inhibitor binding through allosteric effect. A further kinase assay basically confirms the computational findings, which also exhibits a good consistence with theoretical selectivity over 10 tested samples (R_p = 0.89). In particular, the Fasalog compounds Dimehtylfasudil and H-89 are identified as II-o-I and I-o-II selective inhibitors. They can be considered as promising lead molecular entities to develop new specific ROCK isoform-selective Fasalog inhibitors.