New Allosteric Modulators of Metabotropic Glutamate Receptor 5 (mGluR5) Found by Ligand-Based Virtual Screening

Computational Chemistry for the Identification of Lead Compounds for Radiotracer Development

The use of computer-aided drug design (CADD) for the identification of lead compounds in radiotracer development is steadily increasing. Traditional CADD methods, such as structure-based and ligand-based virtual screening and optimization, have been successfully utilized in many drug discovery programs and are highlighted throughout this review. First, we discuss the use of virtual screening for hit identification at the beginning of drug discovery programs. This is followed by an analysis of how the hits derived from virtual screening can be filtered and culled to highly probable candidates to test in in vitro assays. We then illustrate how CADD can be used to optimize the potency of experimentally validated hit compounds from virtual screening for use in positron emission tomography (PET). Finally, we conclude with a survey of the newest techniques in CADD employing machine learning (ML).

Visualization of Topological Pharmacophore Space with Graph Edit Distance

A topological pharmacophore (TP) is a chemical graph-based pharmacophore representation, where nodes are pharmacophoric features (PF) and edges are topological distances between PFs. Previously proposed sparse pharmacophore graphs (SPhGs) for TPs were shown to be effective in identifying structurally different active compounds while maintaining the interpretability of the graphs. However, one limitation of using SPhGs as queries is that many structurally similar SPhGs can be identified from a set of active compounds, requiring the classification and visualization of SPhGs, followed by an understanding of the pharmacophore hypotheses. In this study, we propose a scheme for SPhG analysis based on dimensionality reduction techniques with the graph edit distance (GED) metric. This metric enables measuring similarities among SPhGs in a quantitative manner. The visualization of SPhGs, which themselves are the graphs shared by active compounds, can help us understand the pharmacophore hypotheses as well as the data set. As a proof-of-concept study, we generated two-dimensional SPhG-maps using three dimensionality reduction techniques for six biological targets. A comparison with other pharmacophore representations was also conducted. We demonstrated knowledge extraction (interpretation of the data set) from the generated maps. Our findings include a suitable mapping algorithm as well as a pharmacophore hypothesis analysis procedure using an SPhG-map.

Sparse Topological Pharmacophore Graphs for Interpretable Scaffold Hopping

The aim of scaffold hopping (SH) is to find compounds consisting of different scaffolds from those in already known active compounds, giving an opportunity for unexplored regions of chemical space. We previously demonstrated the usefulness of pharmacophore graphs (PhGs) for this purpose through proof-of-concept virtual screening experiments. PhGs consist of nodes and edges corresponding to pharmacophoric features (PFs) and their topological distances. Although PhGs were effective in SH, they are hard to interpret as they are complete graphs. Herein, we introduce an intuitive representation of a molecule, termed as sparse pharmacophore graphs (SPhG) by keeping the topological distances among PFs as much as possible while reducing the number of edges in the graphs. Several benchmark calculations quantitatively confirmed the sparseness of the graphs and the preservation of topological distances among pharmacophoric points. As proof-of-concept applications, virtual screening (VS) trials for SH were conducted using active and inactive compounds from ChEMBL and PubChem databases for three biological targets: thrombin, tyrosine kinase ABL1, and κ-opioid receptor. The performances of VS were comparable with using fully connected PhGs. Furthermore, highly ranked SPhGs were interpretable for the three biological targets, in particular for thrombin, for which selected SPhGs were in agreement with the structure-based interpretation.

Discovery of potential negative allosteric modulators of mGluR5 from natural products using pharmacophore modeling, molecular docking, and molecular dynamics simulation studies

mGluR5, which belongs to the G-protein-coupled receptor superfamily, is believed to be associated with many human diseases, such as a wide range of neurological disorders, gastroesophageal reflux disease, and cancer. Comparing with compounds that target on the orthosteric binding site, significant roles have been established for mGluR5 negative allosteric modulators (NAMs) due to their higher subtype selectivity and more suitable pharmacokinetic profiles. Nevertheless, to date, none of them have come to market for various reasons. In this study, a 3D quantitative pharmacophore model was generated by using the HypoGen module in Discovery Studio 4.0. With several validation methods ultilized, the optimal pharmacophore model Hypo2 was selected to discover potential mGluR5 NAMs from natural products. Two hundred and seventeen potential NAMs were obtained after being filtered by Lipinski’s rule (≥4). Then, molecular docking was used to refine the pharmacophore-based screening results and analyze the binding mode of NAMs and mGluR5. Three compounds, aglaiduline, 5-O-ethyl-hirsutanonol, and yakuchinone A, with good ADMET properties, acceptable Fit value and estimated value, and high docking score, were reserved for a molecular dynamics simulation study. All of them have stability of ligand binding. From our computational results, there might exhibit drug-like negative allosteric moderating effects on mGluR5 in these natural products. This work provides a reliable method for discovering mGluR5 NAMs from natural products.

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.

Bioisosteric Replacement and Scaffold Hopping in Lead Generation and Optimization

Bioisosteric replacement and scaffold hopping are twin methods used in drug design to improve the synthetic accessibility, potency and drug like properties of a compound and to move into novel chemical space. Bioisosteric replacement involves swapping functional groups of a molecule with other functional groups that have similar biological properties. Scaffold hopping is the replacement of the core framework of a molecule with another scaffold that will improve the properties of the molecule or to find similar potent compounds that exist in novel chemical space. This review outlines the key concepts, importance and challenges of both methods using examples and comparisons of techniques available for finding bioisosteric replacements and scaffold hops. There are many methods available for bioisosteric replacement and scaffold hopping, all with their own advantages and disadvantages. Drug design projects would benefit from a combination of these methods to retrieve diverse and complimentary results. Continuing progress in these fields will allow further validation of both methods as well as the accumulation of knowledge on bioisosteres and possible scaffold replacements.

Synergism of virtual screening and medicinal chemistry: identification and optimization of allosteric antagonists of metabotropic glutamate receptor 1

We report the identification of novel potent and selective metabotropic glutamate receptor 1 (mGluR1) antagonists by virtual screening and subsequent hit optimization. For ligand-based virtual screening, molecules were represented by a topological pharmacophore descriptor (CATS-2D) and clustered by a self-organizing map (SOM). The most promising compounds were tested in mGluR1 functional and binding assays. We identified a potent chemotype exhibiting selective antagonistic activity at mGluR1 (functional IC(50)=0.74+/-0.29 microM). Hit optimization yielded lead structure 16 with an affinity of K(i)=0.024+/-0.001 microM and greater than 1000-fold selectivity for mGluR1 versus mGluR5. Homology-based receptor modelling suggests a binding site compatible with previously reported mutation studies. Our study demonstrates the usefulness of ligand-based virtual screening for scaffold-hopping and rapid lead structure identification in early drug discovery projects.

Scaffold hopping in drug discovery using inductive logic programming

In chemoinformatics, searching for compounds which are structurally diverse and share a biological activity is called scaffold hopping. Scaffold hopping is important since it can be used to obtain alternative structures when the compound under development has unexpected side-effects. Pharmaceutical companies use scaffold hopping when they wish to circumvent prior patents for targets of interest. We propose a new method for scaffold hopping using inductive logic programming (ILP). ILP uses the observed spatial relationships between pharmacophore types in pretested active and inactive compounds and learns human-readable rules describing the diverse structures of active compounds. The ILP-based scaffold hopping method is compared to two previous algorithms (chemically advanced template search, CATS, and CATS3D) on 10 data sets with diverse scaffolds. The comparison shows that the ILP-based method is significantly better than random selection while the other two algorithms are not. In addition, the ILP-based method retrieves new active scaffolds which were not found by CATS and CATS3D. The results show that the ILP-based method is at least as good as the other methods in this study. ILP produces human-readable rules, which makes it possible to identify the three-dimensional features that lead to scaffold hopping. A minor variant of a rule learnt by ILP for scaffold hopping was subsequently found to cover an inhibitor identified by an independent study. This provides a successful result in a blind trial of the effectiveness of ILP to generate rules for scaffold hopping. We conclude that ILP provides a valuable new approach for scaffold hopping.

A new series of pyridinyl-alkynes as antagonists of the metabotropic glutamate receptor 5 (mGluR5)

Synthesis and some structure-activity relationships for a new series of propargyl ethers as mGluR5 antagonists are reported.