Novel Inhibitory Site Revealed by Mode of Action on the Metabotropic Glutamate 7 Receptor Venus Flytrap Domain

Metabotropic glutamate (mGlu) receptors play a key role in modulating most synapses in the brain. The mGlu7 receptors inhibit presynaptic neurotransmitter release and offer therapeutic possibilities for post-traumatic stress disorders or epilepsy. Screening campaigns provided mGlu7-specific allosteric modulators as the inhibitor XAP044 (Gee et al. J. Biol. Chem. 2014). In contrast to other mGlu receptor allosteric modulators, XAP044 does not bind in the transmembrane domain but to the extracellular domain of the mGlu7 receptor and not at the orthosteric site. Here, we identified the mode of action of XAP044, combining synthesis of derivatives, modeling and docking experiments, and mutagenesis. We propose a unique mode of action of these inhibitors, preventing the closure of the Venus flytrap agonist binding domain. While acting as a noncompetitive antagonist of L-AP4, XAP044 and derivatives act as apparent competitive antagonists of LSP4-2022. These data revealed more potent XAP044 analogues and new possibilities to target mGluRs.

Amino Acids Bearing Aromatic or Heteroaromatic Substituents as a New Class of Ligands for the Lysosomal Sialic Acid Transporter Sialin

Sialin, encoded by the SLC17A5 gene, is a lysosomal sialic acid transporter defective in Salla disease, a rare inherited leukodystrophy. It also enables metabolic incorporation of exogenous sialic acids, leading to autoantibodies against N-glycolylneuraminic acid in humans. Here, we identified a novel class of human sialin ligands by virtual screening and structure-activity relationship studies. The ligand scaffold is characterized by an amino acid backbone with a free carboxylate, an N-linked aromatic or heteroaromatic substituent, and a hydrophobic side chain. The most potent compound, 45 (LSP12-3129), inhibited N-acetylneuraminic acid 1 (Neu5Ac) transport in a non-competitive manner with IC₅₀ ≈ 2.5 μM, a value 400-fold lower than the K_M for Neu5Ac. In vitro and molecular docking studies attributed the non-competitive character to selective inhibitor binding to the Neu5Ac site in a cytosol-facing conformation. Moreover, compound 45 rescued the trafficking defect of the pathogenic mutant (R39C) causing Salla disease. This new class of cell-permeant inhibitors provides tools to investigate the physiological roles of sialin and help develop pharmacological chaperones for Salla disease.

LSP5-2157 a new inhibitor of vesicular glutamate transporters

Vesicular glutamate transporters (VGLUT1-3) mediate the uptake of glutamate into synaptic vesicles. VGLUTs are pivotal actors of excitatory transmission and of almost all brain functions. Their implication in various pathologies has been clearly documented. Despite their functional importance, the pharmacology of VGLUTs is limited to a few dyes such as Trypan Blue, Rose Bengal or Brilliant Yellow type. Here, we report the design and evaluation of new potent analogs based on Trypan Blue scaffold. Our best compound, named LSP5-2157, has an EC50 of 50 nM on glutamate vesicular uptake. Using a 3D homology model of VGLUT1 and docking experiments, we determined its putative binding subdomains within vesicular glutamate transporters and validated the structural requirement for VGLUT inhibition. To better estimate the specificity and potency of LSP5-2157, we also investigated its ability to block glutamatergic transmission in autaptic hippocampal cells. Neither glutamate receptors nor GABAergic transmission or transmission machinery were affected by LSP5-2157. Low doses of compound reversibly reduce glutamatergic neurotransmission in hippocampal autpases. LSP5-2157 had a low and depressing effect on synaptic efficacy in hippocampal slice. Furthermore, LSP5-2157 had no effect on NMDA-R- mediated fEPSP but reduce synaptic plasticity induced by 3 trains of 100 Hz. Finally, LSP5-2157 had the capacity to inhibit VGLUT3-dependent auditory synaptic transmission in the guinea pig cochlea. In this model, it abolished the compound action potential of auditory nerve at high concentration showing the limited permeation of LSP5-2157 in an in-vivo model. In summary, the new ligand LSP5-2157, has a high affinity and specificity for VGLUTs and shows some permeability in isolated neuron, tissue preparations or in vivo in the auditory system. These findings broaden the field of VGLUTs inhibitors and open the way to their use to assess glutamatergic functions in vitro and in vivo.

Successful prediction of substrate-binding pocket in SLC17 transporter sialin

Secondary active transporters from the SLC17 protein family are required for excitatory and purinergic synaptic transmission, sialic acid metabolism, and renal function, and several members are associated with inherited neurological or metabolic diseases. However, molecular tools to investigate their function or correct their genetic defects are limited or absent. Using structure-activity, homology modeling, molecular docking, and mutagenesis studies, we have located the substrate-binding site of sialin (SLC17A5), a lysosomal sialic acid exporter also recently implicated in exocytotic release of aspartate. Human sialin is defective in two inherited sialic acid storage diseases and is responsible for metabolic incorporation of the dietary nonhuman sialic acid N-glycolylneuraminic acid. We built cytosol-open and lumen-open three-dimensional models of sialin based on weak, but significant, sequence similarity with the glycerol-3-phosphate and fucose permeases from Escherichia coli, respectively. Molecular docking of 31 synthetic sialic acid analogues to both models was consistent with inhibition studies. Narrowing the sialic acid-binding site in the cytosol-open state by two phenylalanine to tyrosine mutations abrogated recognition of the most active analogue without impairing neuraminic acid transport. Moreover, a pilot virtual high-throughput screening of the cytosol-open model could identify a pseudopeptide competitive inhibitor showing >100-fold higher affinity than the natural substrate. This validated model of human sialin and sialin-guided models of other SLC17 transporters should pave the way for the identification of inhibitors, glycoengineering tools, pharmacological chaperones, and fluorescent false neurotransmitters targeted to these proteins.

Cinnabarinic acid, an endogenous metabolite of the kynurenine pathway, activates type 4 metabotropic glutamate receptors

Cinnabarinic acid is an endogenous metabolite of the kynurenine pathway that meets the structural requirements to interact with glutamate receptors. We found that cinnabarinic acid acts as a partial agonist of type 4 metabotropic glutamate (mGlu4) receptors, with no activity at other mGlu receptor subtypes. We also tested the activity of cinnabarinic acid on native mGlu4 receptors by examining 1) the inhibition of cAMP formation in cultured cerebellar granule cells; 2) protection against excitotoxic neuronal death in mixed cultures of cortical cells; and 3) protection against 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine toxicity in mice after local infusion into the external globus pallidus. In all these models, cinnabarinic acid behaved similarly to conventional mGlu4 receptor agonists, and, at least in cultured neurons, the action of low concentrations of cinnabarinic acid was largely attenuated by genetic deletion of mGlu4 receptors. However, high concentrations of cinnabarinic acid were still active in the absence of mGlu4 receptors, suggesting that the compound may have off-target effects. Mutagenesis and molecular modeling experiments showed that cinnabarinic acid acts as an orthosteric agonist interacting with residues of the glutamate binding pocket of mGlu4. Accordingly, cinnabarinic acid did not activate truncated mGlu4 receptors lacking the N-terminal Venus-flytrap domain, as opposed to the mGlu4 receptor enhancer, N-phenyl-7-(hydroxyimino)cyclopropa[b]chromen-1a-carboxamide (PHCCC). Finally, we could detect endogenous cinnabarinic acid in brain tissue and peripheral organs by high-performance liquid chromatography-tandem mass spectrometry analysis. Levels increased substantially during inflammation induced by lipopolysaccharide. We conclude that cinnabarinic acid is a novel endogenous orthosteric agonist of mGlu4 receptors endowed with neuroprotective activity.

Predicting fragment binding poses using a combined MCSS MM-GBSA approach

Improved methods are required to predict the position and orientation (pose) of binding to the target protein of low molecular weight compounds identified in fragment screening campaigns. This is particularly important to guide initial chemistry to generate structure-activity relationships for the cases where a high resolution structure cannot be obtained. We have assessed the benefit of an implicit solvent method for assessment of fragment binding poses generated by the Multiple Copy Simultaneous Search (MCSS) method in CHARMm. Additionally, the effect of using multiple receptor structures for a flexible receptor is investigated. The original MCSS performance -50% of fragment positions accurately predicted and scored - was increased up to 67% by scoring MCSS energy minima with a Molecular Mechanics Generalized Born approach with molecular volume integration and Surface Area model (MM-GBSA). The same increase in performance (but occasionally for different targets) was observed when using the docking program GOLD followed by MM-GBSA rescoring. The combined results from both methods resulted in a higher success rate emphasizing that a comparison of different docking methods can increase the correct identification of binding poses. For a receptor where multiple structures are available, Hsp90, the average performance on randomly adding receptor structures was also investigated. The results suggest that predictions using these docking methods can be used with some confidence to guide chemical optimization, if the structure of the target either remains relatively fixed on ligand binding, or if a number of crystal structures are available with diverse ligands bound and there is information on the positions of key water molecules in the binding site.

A virtual screening hit reveals new possibilities for developing group III metabotropic glutamate receptor agonists

(R)-PCEP (3-amino-3-carboxypropyl-2'-carboxyethyl phosphinic acid, 1), a new metabotropic glutamate receptor 4 (mGlu4R) agonist, was discovered in a previously reported virtual screening. The (S)-enantiomer and a series of derivatives were synthesized and tested on recombinant mGlu4 receptors. A large number of derivatives activated this receptor but was not able to discriminate between mGlu4 and mGlu8 receptors. The most potent ones 6 and 12 displayed an EC(50) of 1.0 +/- 0.2 microM at mGlu4R. Interestingly these agonists with longer alkyl chains revealed a new binding pocket adjacent to the glutamate binding site, which is lined with residues that differ among the mGluR subtypes and that will allow the design of more selective compounds. Additionally 6 was able to activate mGlu7 receptor with an EC(50) of 43 +/- 16 microM and is thus significantly more potent than L-AP4 (EC(50) of 249 +/- 106 microM).

High-potency olfactory receptor agonists discovered by virtual high-throughput screening: molecular probes for receptor structure and olfactory function

The detection of diverse chemical structures by the vertebrate olfactory system is accomplished by the recognition of odorous ligands by their cognate receptors. In the present study, we used computational screening to discover novel high-affinity agonists of an olfactory G protein-coupled receptor that recognizes amino acid ligands. Functional testing of the top candidates validated several agonists with potencies higher than any of the receptor's known natural ligands. Computational modeling revealed molecular interactions involved in ligand binding and further highlighted interactions that have been conserved in evolutionarily divergent amino acid receptors. Significantly, the top compounds display robust activities as odorants in vivo and include a natural product that may be used to signal the presence of bacteria in the environment. Our virtual screening approach should be applicable to the identification of new bioactive molecules for probing the structure of chemosensory receptors and the function of chemosensory systems in vivo.

Synthesis and biological evaluation of 1-amino-2-phosphonomethylcyclopropanecarboxylic acids, new group III metabotropic glutamate receptor agonists

Stereoisomers of 1-amino-2-phosphonomethylcyclopropanecarboxylic acid (APCPr), conformationally restricted analogues of L-AP4 (2-amino-4-phosphonobutyric acid), have been prepared and evaluated at recombinant group III metabotropic glutamate receptors. They activate these receptors over a broad range of potencies. The most potent isomer (1S,2R)-APCPr displays a similar pharmacological profile as that of L-AP4 (EC50 0.72, 1.95, >500, 0.34 microM at mGlu4, 6, 7, 8 receptors, respectively, and no effect at group I/II mGluRs). It was characterized on native receptors located in the basal ganglia (BG) where it induced a robust and reversible inhibition of synaptic transmission. It was tested in vivo in haloperidol-induced catalepsy, a model of Parkinsonian akinesia, by direct infusion in the globus pallidus of the BG. At a dose of 0.5 nmol/microL, catalepsy was significantly antagonized. This study reveals that (1S,2R)-APCPr is a potent group III mGluR agonist and confirms that these receptors may be considered as a therapeutic target in the Parkinson's disease.

Amino-pyrrolidine tricarboxylic acids give new insight into group III metabotropic glutamate receptor activation mechanism

Like most class C G-protein-coupled receptors, metabotropic glutamate (mGlu) receptors possess a large extracellular domain where orthosteric ligands bind. Crystal structures revealed that this domain, called Venus FlyTrap (VFT), adopts a closed or open conformation upon agonist or antagonist binding, respectively. We have described amino-pyrrolidine tricarboxylic acids (APTCs), including (2S,4S)-4-amino-1-[(E)-3-carboxyacryloyl]pyrrolidine-2,4-dicarboxylic acid (FP0429), as new selective group III mGlu agonists. Whereas FP0429 is an almost full mGlu4 agonist, it is a weak and partial agonist of the closely related mGlu8 subtype. To get more insight into the activation mechanism of mGlu receptors, we aimed to elucidate why FP0429 behaves differently at these two highly homologous receptors by focusing on two residues within the binding site that differ between mGlu4 and mGlu8. Site-directed mutagenesis of Ser157 and Gly158 of mGlu4 into their mGlu8 homologs (Ala) turned FP0429 into a weak partial agonist. Conversely, introduction of Ser and Gly residues into mGlu8 increased FP0429 efficacy. Docking of FP0429 in mGlu4 VFT 3D model helped us characterize the role of each residue. Indeed, mGlu4 Ser157 seems to have an important role in FP0429 binding, whereas Gly158 may allow a deeper positioning of this agonist in the cavity of lobe I, thereby ensuring optimal interactions with lobe II residues in the fully closed state of the VFT. In contrast, the presence of a methyl group in mGlu8 (Ala instead of Gly) weakens the interactions with the lobe II residues. This probably results in a less stable or a partially closed form of the mGlu8 VFT, leading to partial receptor activation.

Design and synthesis of APTCs (aminopyrrolidinetricarboxylic acids): identification of a new group III metabotropic glutamate receptor selective agonist

A new family of mGlu receptor orthosteric ligands called APTCs was designed and synthesized using a parallel chemistry approach. Amongst 65 molecules tested on mGlu4, mGlu6 and mGlu8 subtypes, (2S,4S)-4-amino-1-[(E)-3-carboxyacryloyl]pyrrolidine-2,4-dicarboxylic acid (8a06-FP0429) has been shown to be a full mGlu4 agonist and a partial mGlu8 agonist. In addition, 8a06 was shown to be selective versus group I and II mGlu subtypes. A possible explanation for this efficacy difference is proposed by docking experiment performed with molecular model of the receptor.

Amino acid recognition by Venus flytrap domains is encoded in an 8-residue motif

A motif foramino acid recognition by proteins or domains of the periplasmic binding protein-like I superfamily has been identified. An initial pattern of 5 residues was based on a multiple sequence alignment of selected proteins of that fold family and on common structural features observed in the crystal structure of some members of the family [leucine isoleucine valine binding protein (LIVBP), leucine binding protein (LBP), and metabotropic glutamate receptor type 1 (mGlu1R) amino terminal domain)]. This pattern was used against the PIR-NREF sequence database and further refined to retrieve all sequences of proteins that belong to the family and eliminate those that do not belong to it. A motif of 8 residues was finally selected to build up the general signature. A total of 232 sequences were retrieved. They were found to belong to only three families of proteins: bacterial periplasmic binding proteins (PBP, 71 sequences), family 3 (or C) of G-protein coupled receptor (GPCR) (146 sequences), and plant putative ionotropic glutamate receptors (iGluR, 15 sequences). PBPs are known to adopt a bilobate structure also named Venus flytrap domain, or LIVBP domain in the present case. Family 3/C GPCRs are also known to hold such a domain. However, for plant iGluRs, it was previously detected by classical similarity searches but not specifically described. Thus plant iGluRs carry two Venus flytrap domains, one that binds glutamate and an additional one that would be a modulatory LIVBP domain. In some cases, the modulator binding to that domain would be an amino acid.

Molecular determinants of ligand selectivity in a vertebrate odorant receptor

The identification of the chemical structure of an odorant by the vertebrate olfactory system is thought to occur through the combinatorial activity from multiple receptors, each tuned to recognize different chemical features. What are the molecular determinants underlying the selectivity of individual odorant receptors for their cognate ligands? To address this question, we performed molecular modeling and site-directed mutagenesis on the ligand-binding region of two orthologous amino acid odorant receptors belonging to the "C family" of G-protein-coupled receptors in goldfish and zebrafish. We identified the critical ligand-receptor interactions that afford ligand binding as well as selectivity for different amino acids. Moreover, predictions regarding binding pocket structure allowed us to alter, in a predictable manner, the receptor preferences for different ligands. These results reveal how this class of odorant receptor has evolved to accommodate ligands of varying chemical structure and further illuminate the molecular principles underlying ligand recognition and selectivity in this family of chemosensory receptors.

Virtual Screening Workflow Development Guided by the “Receiver Operating Characteristic” Curve Approach. Application to High-Throughput Docking on Metabotropic Glutamate Receptor Subtype 4

The "receiver operating characteristic" (ROC) curve method is a well-recognized metric used as an objective way to evaluate the ability of a given test to discriminate between two populations. This facilitates decision-making in a plethora of fields in which a wrong judgment may have serious consequences including clinical diagnosis, public safety, travel security, and economic strategies. When virtual screening is used to speed-up the drug discovery process in pharmaceutical research, taking the right decision upon selecting or discarding a molecule prior to in vitro evaluation is of paramount importance. Characterizing both the ability of a virtual screening workflow to select active molecules and the ability to discard inactive ones, the ROC curve approach is well suited for this critical decision gate. As a case study, the first virtual screening workflow focused on metabotropic glutamate receptor subtype 4 (mGlu4R) agonists is reported here. Six compounds out of 38 selected and tested in vitro were shown to have agonist activity on this target of therapeutic interest.

Structure-Based Virtual Screening: An Application to Human Topoisomerase II α

The eukaryotic topoisomerase II is involved in several vital processes, such as replication, transcription, and recombination. Many compounds interfering with the catalytic action of this enzyme are efficient in human cancer chemotherapy. We applied a methodology combining molecular modeling and virtual screening techniques to identify human topoisomerase II alphainhibitors. Data from structural biology and enzymatic assays together with a good background on the enzyme mechanism of action were helpful in the approach. A human topoisomerase II alpha model provided an insight into the structural features responsible for the activity of the enzyme. A protocol comprising several substructural and protein structure-based three-dimensional pharmacophore filters enabled the successful retrieving of inhibitors of the enzyme from large databases of compounds, thus validating the approach. A subset of protein structural features required for the enzyme inhibition at the protein-DNA interface were identified and incorporated into the pharmacophore models. Compounds sharing a DNA-intercalating chromophore and a moiety interfering with the protein active site emerged as good inhibitors.

Common and selective molecular determinants involved in metabotopic glutamate receptor agonist activity

Several potent and group selective agonists of metabotropic glutamate receptors (mGluRs) have been docked at mGlu1,2,4R binding sites in the closed conformation of the bilobate extracellular domain. Quisqualic acid and (S)-3,5-dihydroxyphenylglycine (3,5-DHPG) were selected for mGlu1R, dicarboxycyclopropylglycine (DCG-IV), LY354740, (S)-4-carboxyphenylglycine (4CPG) for mGlu2R, and (S)-2-amino-4-phosphonobutyric acid (AP4), 1-aminocyclopentane-1,3,4-tricarboxylic acid (ACPT-I), (S)-4-phosphonophenylglycine (PPG) for mGlu4R. The models show a conserved binding pattern for the glycine moiety (alpha-amino and alpha-acidic functions) and group specific bindings for the distal acidic function. The best agonists allow optimized interaction with both lobes of the binding domain. Interlobe connections around the ligand are also described and participate in stabilizing the closed form of the amino-terminal domain. Altogether, the docking models support the proposal that the stabilization of a closed state represents a key step in agonist activation of mGluRs.

Mapping the agonist-binding site of GABAB type 1 subunit sheds light on the activation process of GABAB receptors

The gamma-amino-n-butyric acid type B (GABA(B)) receptor is composed of two subunits, GABA(B)1 and GABA(B)2, belonging to the family 3 heptahelix receptors. These proteins possess two domains, a seven transmembrane core and an extracellular domain containing the agonist binding site. This binding domain is likely to fold like bacterial periplasmic binding proteins that are constituted of two lobes that close upon ligand binding. Here, using molecular modeling and site-directed mutagenesis, we have identified residues in the GABA(B)1 subunit that are critical for agonist binding and activation of the heteromeric receptor. Our data suggest that two residues (Ser(246) and Asp(471)) located within lobe I form H bonds and a salt bridge with carboxylic and amino groups of GABA, respectively, demonstrating the pivotal role of lobe I in agonist binding. Interestingly, our data also suggest that a residue within lobe II (Tyr(366)) interacts with the agonists in a closed form model of the binding domain, and its mutation into Ala converts the agonist baclofen into an antagonist. These data demonstrate the pivotal role played by the GABA(B)1 subunit in the activation of the heteromeric GABA(B) receptor and are consistent with the idea that a closed state of the binding domain of family 3 receptors is required for their activation.

Three-dimensional model of the extracellular domain of the type 4a metabotropic glutamate receptor: new insights into the activation process

Metabotropic glutamate receptors (mGluRs) belong to the family 3 of G-protein-coupled receptors. On these proteins, agonist binding on the extracellular domain leads to conformational changes in the 7-transmembrane domains required for G-protein activation. To elucidate the structural features that might be responsible for such an activation mechanism, we have generated models of the amino terminal domain (ATD) of type 4 mGluR (mGlu4R). The fold recognition search allowed the identification of three hits with a low sequence identity, but with high secondary structure conservation: leucine isoleucine valine-binding protein (LIVBP) and leucine-binding protein (LBP) as already known, and acetamide-binding protein (AmiC). These proteins are characterized by a bilobate structure in an open state for LIVBP/LBP and a closed state for AmiC, with ligand binding in the cleft. Models for both open and closed forms of mGlu4R ATD have been generated. ACPT-I (1-aminocyclopentane 1,3,4-tricarboxylic acid), a selective agonist, has been docked in the two models. In the open form, ACPT-I is only bound to lobe I through interactions with Lys74, Arg78, Ser159, and Thr182. In the closed form, ACPT-I is trapped between both lobes with additional binding to Tyr230, Asp312, Ser313, and Lys317 from lobe II. These results support the hypothesis that mGluR agonists bind a closed form of the ATDs, suggesting that such a conformation of the binding domain corresponds to the active conformation.

Structure of the human glycogen-associated protein phosphatase 1 regulatory subunit hGM: homology modeling revealed an (alpha/beta)8-barrel-like fold in the multidomain protein

Protein phosphatase 1 (PP1) is widely distributed among tissues and species and acts as a regulator of many important cellular processes. By targeting the catalytic part of PP1 (PP1C) toward particular loci and substrates, regulatory subunits constitute key elements conferring specificity to the holoenzyme. Here, we report the identification of an (alpha/beta)8-barrel-like structure within the N-ter stretch of the human PP1 regulatory subunit hGM, which is part of the family of diverse proteins associated with glycogen metabolism. Protein homology modeling gave rise to a three-dimensional (3D) model for the 381 N-ter residue stretch of hGM, based on sequence similarity with Streptomyces olivochromogenes xylose isomerase, identified by using FASTA. The alignment was subsequently extended by using hydrophobic cluster analysis. The homology-derived model includes the putative glycogen binding area located within the 142-230 domain of hGM as well as a structural characterization of the PP1C interacting domain (segment 51-67). Refinement of the latter by molecular dynamics afforded a topology that is in agreement with previous X-ray studies (Egloff et al., 1997). Finite difference Poisson-Boltzmann calculations performed on the interacting domains of PP1C and hGM confirm the complementarity of the local electrostatic potentials of the two partners. This work highlights the presence of a conserved fold among distant species (mammalian, Caenorhabditis elegans, yeast) and, thus, emphasizes the involvement of PP1 in crucial basic cellular functions.