( S)-3,4-Dicarboxyphenylglycine (DCPG) was first reported in 2001 as a potent orthosteric agonist with high subtype selectivity for the mGlu8 receptor, but the structural basis for its high selectivity is not well understood. We have solved a cocrystal structure of recombinant human mGlu8 amino terminal domain (ATD) protein bound to ( S)-DCPG, which possesses the largest lobe opening angle observed to date among known agonist-bound mGlu ATD crystal structures. The binding conformation of ( S)-DCPG observed in the crystal structure is significantly different from that in the homology model built from an l-glutamate-bound rat mGlu1 ATD crystal structure, which has a smaller lobe opening angle. This highlights the importance of considering various lobe opening angles when modeling mGlu ATD-ligand complex. New homology models of other mGlu receptors based on the ( S)-DCPG-bound mGlu8 ATD crystal structure were explored to rationalize ( S)-DCPG's high mGlu8 receptor subtype selectivity.

Protein function prediction is an active area of research in computational biology. Function prediction can help biologists make hypotheses for characterization of genes and help interpret biological assays, and thus is a productive area for collaboration between experimental and computational biologists. Among various function prediction methods, predicting binding ligand molecules for a target protein is an important class because ligand binding events for a protein are usually closely intertwined with the proteins' biological function, and also because predicted binding ligands can often be directly tested by biochemical assays. Binding ligand prediction methods can be classified into two types: those which are based on protein-protein (or pocket-pocket) comparison, and those that compare a target pocket directly to ligands. Recently, our group proposed two computational binding ligand prediction methods, Patch-Surfer, which is a pocket-pocket comparison method, and PL-PatchSurfer, which compares a pocket to ligand molecules. The two programs apply surface patch-based descriptions to calculate similarity or complementarity between molecules. A surface patch is characterized by physicochemical properties such as shape, hydrophobicity, and electrostatic potentials. These properties on the surface are represented using three-dimensional Zernike descriptors (3DZD), which are based on a series expansion of a 3 dimensional function. Utilizing 3DZD for describing the physicochemical properties has two main advantages: (1) rotational invariance and (2) fast comparison. Here, we introduce Patch-Surfer and PL-PatchSurfer with an emphasis on PL-PatchSurfer, which is more recently developed. Illustrative examples of PL-PatchSurfer performance on binding ligand prediction as well as virtual drug screening are also provided.

As part of our ongoing efforts to identify novel ligands for the metabotropic glutamate 2 and 3 (mGlu2/3) receptors, we have incorporated substitution at the C3 and C4 positions of the (1S,2R,5R,6R)-2-amino-bicyclo[3.1.0]hexane-2,6-dicarboxylic acid scaffold to generate mGlu2/3 antagonists. Exploration of this structure-activity relationship (SAR) led to the identification of (1S,2R,3S,4S,5R,6R)-2-amino-3-[(3,4-difluorophenyl)sulfanylmethyl]-4-hydroxy-bicyclo[3.1.0]hexane-2,6-dicarboxylic acid hydrochloride (LY3020371·HCl, 19f), a potent, selective, and maximally efficacious mGlu2/3 antagonist. Further characterization of compound 19f binding to the human metabotropic 2 glutamate (hmGlu2) site was established by cocrystallization of this molecule with the amino terminal domain (ATD) of the hmGlu2 receptor protein. The resulting cocrystal structure revealed the specific ligand-protein interactions, which likely explain the high affinity of 19f for this site and support its functional mGlu2 antagonist pharmacology. Further characterization of 19f in vivo demonstrated an antidepressant-like signature in the mouse forced-swim test (mFST) assay when brain levels of this compound exceeded the cellular mGlu2 IC50 value.

Muscarinic M1-M5 acetylcholine receptors are G-protein-coupled receptors that regulate many vital functions of the central and peripheral nervous systems. In particular, the M1 and M4 receptor subtypes have emerged as attractive drug targets for treatments of neurological disorders, such as Alzheimer's disease and schizophrenia, but the high conservation of the acetylcholine-binding pocket has spurred current research into targeting allosteric sites on these receptors. Here we report the crystal structures of the M1 and M4 muscarinic receptors bound to the inverse agonist, tiotropium. Comparison of these structures with each other, as well as with the previously reported M2 and M3 receptor structures, reveals differences in the orthosteric and allosteric binding sites that contribute to a role in drug selectivity at this important receptor family. We also report identification of a cluster of residues that form a network linking the orthosteric and allosteric sites of the M4 receptor, which provides new insight into how allosteric modulation may be transmitted between the two spatially distinct domains.

As part of our ongoing research to identify novel agents acting at metabotropic glutamate 2 (mGlu2) and 3 (mGlu3) receptors, we have previously reported the identification of the C4α-methyl analog of mGlu2/3 receptor agonist 1 (LY354740). This molecule, 1S,2S,4R,5R,6S-2-amino-4-methylbicyclo[3.1.0]hexane-2,6-dicarboxylate 2 (LY541850), exhibited an unexpected mGlu2 agonist/mGlu3 antagonist pharmacological profile, whereas the C4β-methyl diastereomer (3) possessed dual mGlu2/3 receptor agonist activity. We have now further explored this structure-activity relationship through the preparation of cyclic and acyclic C4-disubstituted analogs of 1, leading to the identification of C4-spirocyclopropane 5 (LY2934747), a novel, potent, and systemically bioavailable mGlu2/3 receptor agonist which exhibits both antipsychotic and analgesic properties in vivo. In addition, through the combined use of protein-ligand X-ray crystallography employing recombinant human mGlu2/3 receptor amino terminal domains, molecular modeling, and site-directed mutagenesis, a molecular basis for the observed pharmacological profile of compound 2 is proposed.

The observation that cholinergic deafferentation of circuits projecting from forebrain basal nuclei to frontal and hippocampal circuits occurs in Alzheimer's disease has led to drug-targeting of muscarinic M1 receptors to alleviate cognitive symptoms. The high homology within the acetylcholine binding domain of this family however has made receptor-selective ligand development challenging. This work presents the synthesis scheme, pharmacokinetic and structure-activity-relationship study findings for M1-selective ligand, LY593093. Pharmacologically the compound acts as an orthosteric ligand. The homology modeling work presented however will illustrate that compound binding spans from the acetylcholine pocket to the extracellular loops of the receptor, a common allosteric vestibule for the muscarinic protein family. Altogether LY593093 represents a growing class of multi-topic ligands which interact with the receptors in both the ortho- and allosteric binding sites, but which exert their activation mechanism as an orthosteric ligand.

Virtual screening has been widely used in the drug discovery process. Ligand-based virtual screening (LBVS) methods compare a library of compounds with a known active ligand. Two notable advantages of LBVS methods are that they do not require structural information of a target receptor and that they are faster than structure-based methods. LBVS methods can be classified based on the complexity of ligand structure information utilized: one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D). Unlike 1D and 2D methods, 3D methods can have enhanced performance since they treat the conformational flexibility of compounds. In this paper, a number of 3D methods will be reviewed. In addition, four representative 3D methods were benchmarked to understand their performance in virtual screening. Specifically, we tested overall performance in key aspects including the ability to find dissimilar active compounds, and computational speed.

Identification of orthosteric mGlu(2/3) receptor agonists capable of discriminating between individual mGlu2 and mGlu3 subtypes has been highly challenging owing to the glutamate-site sequence homology between these proteins. Herein we detail the preparation and characterization of a series of molecules related to (1S,2S,5R,6S)-2-aminobicyclo[3.1.0]hexane-2,6-dicarboxylate 1 (LY354740) bearing C4-thiotriazole substituents. On the basis of second messenger responses in cells expressing other recombinant human mGlu2/3 subtypes, a number of high potency and efficacy mGlu2 receptor agonists exhibiting low potency mGlu3 partial agonist/antagonist activity were identified. From this, (1R,2S,4R,5R,6R)-2-amino-4-(1H-1,2,4-triazol-3-ylsulfanyl)bicyclo[3.1.0]hexane-2,6-dicarboxylic acid 14a (LY2812223) was further characterized. Cocrystallization of 14a with the amino terminal domains of hmGlu2 and hmGlu3 combined with site-directed mutation studies has clarified the underlying molecular basis of this unique pharmacology. Evaluation of 14a in a rat model responsive to mGlu2 receptor activation coupled with a measure of central drug disposition provides evidence that this molecule engages and activates central mGlu2 receptors in vivo.

Structure-based computational methods have been widely used in exploring protein-ligand interactions, including predicting the binding ligands of a given protein based on their structural complementarity. Compared to other protein and ligand representations, the advantages of a surface representation include reduced sensitivity to subtle changes in the pocket and ligand conformation and fast search speed. Here we developed a novel method named PL-PatchSurfer (Protein-Ligand PatchSurfer). PL-PatchSurfer represents the protein binding pocket and the ligand molecular surface as a combination of segmented surface patches. Each patch is characterized by its geometrical shape and the electrostatic potential, which are represented using the 3D Zernike descriptor (3DZD). We first tested PL-PatchSurfer on binding ligand prediction and found it outperformed the pocket-similarity based ligand prediction program. We then optimized the search algorithm of PL-PatchSurfer using the PDBbind dataset. Finally, we explored the utility of applying PL-PatchSurfer to a larger and more diverse dataset and showed that PL-PatchSurfer was able to provide a high early enrichment for most of the targets. To the best of our knowledge, PL-PatchSurfer is the first surface patch-based method that treats ligand complementarity at protein binding sites. We believe that using a surface patch approach to better understand protein-ligand interactions has the potential to significantly enhance the design of new ligands for a wide array of drug-targets.

As part of our ongoing interest in identifying novel agonists acting at metabotropic glutamate (mGlu) 2/3 receptors, we have explored the effect of structural modifications of 1S,2S,5R,6S-2-aminobicyclo[3.1.0]hexane-2,6-dicarboxylate (LY354740), a potent and pharmacologically balanced mGlu2/3 receptor agonist. Incorporation of relatively small substituents (e.g., F, O) at the C4 position of this molecule resulted in additional highly potent mGlu2/3 agonists that demonstrate excellent selectivity over the other mGlu receptor subtypes, while addition of larger C4-substituents (e.g., SPh) led to a loss of agonist potency and/or the appearance of weak mGlu2/3 receptor antagonist activity. Further characterization of the α-fluoro-substituted analogue (LY459477) in vivo revealed that this molecule possesses good oral bioavailability in rats and effectively suppresses phencyclidine-evoked locomotor activity at doses that do not impair neuromuscular coordination. This molecule therefore represents a valuable new addition to the arsenal of pharmacological tools competent to investigate mGlu2/3 receptor function both in vitro and in vivo.