Linear and nonlinear 3D-QSAR approaches in tandem with ligand-based homology modeling as a computational strategy to depict the pyrazolo-triazolo-pyrimidine antagonists binding site of the human adenosine A2A receptor

The integration of ligand- and structure-based strategies might sensitively increase the success of drug discovery process. We have recently described the application of Molecular Electrostatic Potential autocorrelated vectors (autoMEPs) in generating both linear (Partial Least-Square, PLS) and nonlinear (Response Surface Analysis, RSA) 3D-QSAR models to quantitatively predict the binding affinity of human adenosine A3 receptor antagonists. Moreover, we have also reported a novel GPCR modeling approach, called Ligand-Based Homology Modeling (LBHM), as a tool to simulate the conformational changes of the receptor induced by ligand binding. In the present study, the application of both linear and nonlinear 3D-QSAR methods and LBHM computational techniques has been used to depict the hypothetical antagonist binding site of the human adenosine A2A receptor. In particular, a collection of 127 known human A2A antagonists has been utilized to derive two 3D-QSAR models (autoMEPs/PLS&RSA). In parallel, using a rhodopsin-driven homology modeling approach, we have built a model of the human adenosine A2A receptor. Finally, 3D-QSAR and LBHM strategies have been utilized to predict the binding affinity of five new human A2A pyrazolo-triazolo-pyrimidine antagonists finding a good agreement between the theoretical and the experimental predictions.     

Using silico methods predicting ligands for orphan GPCRs

The G-protein coupled receptor (GPCR) superfamily is one of the most important drug target classes for the pharmaceutical industry. The completion of the human genome project has revealed that there are more than 300 potential GPCR targets of interest. The identification of their natural ligands can gain significant insights into regulatory mechanisms of cellular signaling networks and provide unprecedented opportunities for drug discovery. Much effort has been directed towards the GPCR ligand discovery study by both academic institutions and pharmaceutical industries. However, the endogenous ligands still remain unknown for about 150 GPCRs in the human genome. It is necessary to develop new strategies to predict candidate ligands for these so-called orphan receptors. Computational techniques are playing an increasingly important role in finding and validating novel ligands for orphan GPCRs (oGPCRs). In this paper, we focus on recent development in applying bioinformatics approaches for the discovery of GPCR ligands. In addition, some of the data resources for ligand identification are also provided.     

Development of biologically active compounds by combining 3D QSAR and structure-based design methods

One of the major challenges in computational approaches to drug design is the accurate prediction of the binding affinity of novel biomolecules. In the present study an automated procedure which combines docking and 3D-QSAR methods was applied to several drug targets. The developed receptor-based 3D-QSAR methodology was tested on several sets of ligands for which the three-dimensional structure of the target protein has been solved – namely estrogen receptor, acetylcholine esterase and protein-tyrosine-phosphatase 1B. The molecular alignments of the studied ligands were determined using the docking program AutoDock and were compared with the X-ray structures of the corresponding protein-ligand complexes. The automatically generated protein-based ligand alignment obtained was subsequently taken as basis for a comparative field analysis applying the GRID/GOLPE approach. Using GRID interaction fields and applying variable selection procedures, highly predictive models were obtained. It is expected that concepts from receptor-based 3D QSAR will be valuable tools for the analysis of high-throughput screening as well as virtual screening data     

Quantitative structure-based design: formalism and application of receptor-dependent RD-4D-QSAR analysis to a set of glucose analogue inhibitors of glycogen phosphorylase

A method for performing quantitative structure-based design has been developed by extending the current receptor-independent RI-4D-QSAR methodology to include receptor geometry. The resultant receptor-dependent RD-4D-QSAR approach employs a novel receptor-pruning technique to permit effective processing of ligands with the lining of the binding site wrapped about them. Data reduction, QSAR model construction, and identification of possible pharmacophore sites are achieved by a three-step statistical analysis consisting of genetic algorithm optimization followed by backward elimination multidimensional regression and ending with another genetic algorithm optimization. The RD-4D-QSAR method is applied to a series of glucose inhibitors of glycogen phosphorylase b, GPb. The statistical quality of the best RI- and RD-4D-QSAR models are about the same. However, the predictivity of the RD- model is quite superior to that of the RI-4D-QSAR model for a test set. The superior predictive performance of the RD- model is due to its dependence on receptor geometry. There is a unique induced-fit between each inhibitor and the GPb binding site. This induced-fit results in the side chain of Asn-284 serving as both a hydrogen bond acceptor and donor site depending upon inhibitor structure. The RD-4D-QSAR model strongly suggests that quantitative structure-based design cannot be successful unless the receptor is allowed to be completely flexible.     

Structure-based validation of the 3D-QSAR technique MaP

For three target proteins with different binding pocket characteristics (size and shape, hydrophobicity, hydrogen-bonding) a structure-based validation of the translationally and rotationally invariant 3D-QSAR technique MaP is performed (MaP: Mapping Property distributions of molecular surfaces). The structure-based validation procedure comprises two steps: first, QSAR models are derived without using the information of the target protein. Second, the models are back-projected into the crystal structure of the binding pockets and interpreted. It is demonstrated that MaP is able to identify characteristics important for ligand binding in the cases studied here. Moreover, it is demonstrated that MaP is a versatile 3D-QSAR technique since good, predictive models could be obtained for all three data sets showing distinct characteristics.     

G protein- and agonist-bound serotonin 5-HT2A receptor model activated by steered molecular dynamics simulations

A 5-HT(2A) receptor model was constructed by homology modeling based on the β(2)-adrenergic receptor and the G protein-bound opsin crystal structures. The 5-HT(2A) receptor model was transferred into an active conformation by an agonist ligand and a G(αq) peptide in four subsequent steered molecular dynamics (MD) simulations. The driving force for the transformation was the addition of several known intermolecular and receptor interhelical hydrogen bonds enforcing the necessary helical and rotameric movements. Subsquent MD simulations without constraints confirmed the stability of the activated receptor model as well as revealed new information about stabilizing residues and bonds. The active 5-HT(2A) receptor model was further validated by retrospective ligand screening of more than 9400 compounds, whereof 182 were known ligands. The results show that the model can be used in drug discovery for virtual screening and structure-based ligand design as well as in GPCR activation studies.     

Small Molecule Tools to Study Cellular Target Engagement for the Intracellular Allosteric Binding Site of GPCRs

A conserved intracellular allosteric binding site (IABS) has recently been identified at several G protein-coupled receptors (GPCRs). Ligands targeting the IABS, so-called intracellular allosteric antagonists, are highly promising compounds for pharmaceutical intervention and currently evaluated in several clinical trials. Beside co-crystal structures that laid the foundation for the structure-based development of intracellular allosteric GPCR antagonists, small molecule tools that enable an unambiguous identification and characterization of intracellular allosteric GPCR ligands are of utmost importance for drug discovery campaigns in this field. Herein, we discuss recent approaches that leverage cellular target engagement studies for the IABS and thus play a critical role in the evaluation of IABS-targeted ligands as potential therapeutic agents.     

两类促肾上腺皮质释放激素(CRH)拮抗剂的分子场分析

采用比较分子相似性指数分析方法(CoMSIA)及比较分子场分析方法(CoMSIA)研究了两组CRH拮抗剂结构与活性的关系。在两种方法中,都考虑了静电场、立体场以及氢键场对构效关系的影响,结果表明采用CoMSIA得到构效关系模型要明显优于采用CoMFA得到的构效关系模型,在CoMSIA计算中,当引入疏水场时,三维构效关系模型能得到明显的改善,通过这个三维构效关系模型,可以较为精确地预测化合物的活性。通过分析分子场等值面图在空间的分布,可以观察到叠合分子周围的立体、静电以及疏水特征对化合物活性的影响。     

Cheminformatics meets molecular mechanics: a combined application of knowledge-based pose scoring and physical force field-based hit scoring functions improves the accuracy of structure-based virtual screening

Poor performance of scoring functions is a well-known bottleneck in structure-based virtual screening (VS), which is most frequently manifested in the scoring functions' inability to discriminate between true ligands vs known nonbinders (therefore designated as binding decoys). This deficiency leads to a large number of false positive hits resulting from VS. We have hypothesized that filtering out or penalizing docking poses recognized as non-native (i.e., pose decoys) should improve the performance of VS in terms of improved identification of true binders. Using several concepts from the field of cheminformatics, we have developed a novel approach to identifying pose decoys from an ensemble of poses generated by computational docking procedures. We demonstrate that the use of target-specific pose (scoring) filter in combination with a physical force field-based scoring function (MedusaScore) leads to significant improvement of hit rates in VS studies for 12 of the 13 benchmark sets from the clustered version of the Database of Useful Decoys (DUD). This new hybrid scoring function outperforms several conventional structure-based scoring functions, including XSCORE::HMSCORE, ChemScore, PLP, and Chemgauss3, in 6 out of 13 data sets at early stage of VS (up 1% decoys of the screening database). We compare our hybrid method with several novel VS methods that were recently reported to have good performances on the same DUD data sets. We find that the retrieved ligands using our method are chemically more diverse in comparison with two ligand-based methods (FieldScreen and FLAP::LBX). We also compare our method with FLAP::RBLB, a high-performance VS method that also utilizes both the receptor and the cognate ligand structures. Interestingly, we find that the top ligands retrieved using our method are highly complementary to those retrieved using FLAP::RBLB, hinting effective directions for best VS applications. We suggest that this integrative VS approach combining cheminformatics and molecular mechanics methodologies may be applied to a broad variety of protein targets to improve the outcome of structure-based drug discovery studies.     

R-group template CoMFA combines benefits of “ad hoc” and topomer alignments using 3D-QSAR for lead optimization

Template CoMFA methodologies extend topomer CoMFA by allowing user-designated templates, for example the experimental receptor-bound conformation of a prototypical ligand, to help determine the alignment of training and test set structures for 3D-QSAR. The algorithms that generate its new structural modality, template-constrained topomers, are described. Template CoMFA's resolution of certain topomer CoMFA concerns, by providing user control of topological consistency and structural acceptability, is demonstrated for sixteen 3D-QSAR training sets, in particular the Selwood dataset.     

Design of Drug Efficacy Guided by Free Energy Simulations of the β2-Adrenoceptor

G-protein-coupled receptors (GPCRs) play important roles in physiological processes and are modulated by drugs that either activate or block signaling. Rational design of the pharmacological efficacy profiles of GPCR ligands could enable the development of more efficient drugs, but is challenging even if high-resolution receptor structures are available. We performed molecular dynamics simulations of the β₂ adrenergic receptor in active and inactive conformations to assess if binding free energy calculations can predict differences in ligand efficacy for closely related compounds. Previously identified ligands were successfully classified into groups with comparable efficacy profiles based on the calculated shift in ligand affinity upon activation. A series of ligands were then predicted and synthesized, leading to the discovery of partial agonists with nanomolar potencies and novel scaffolds. Our results demonstrate that free energy simulations enable design of ligand efficacy and the same approach can be applied to other GPCR drug targets.     

Homology model-based virtual screening for GPCR ligands using docking and target-biased scoring

The current study investigates the combination of two recently reported techniques for the improvement of homology model-based virtual screening for G-protein coupled receptor (GPCR) ligands. First, ligand-supported homology modeling was used to generate receptor models that were in agreement with mutagenesis data and structure-activity relationship information of the ligands. Second, interaction patterns from known ligands to the receptor were applied for scoring and rank ordering compounds from a virtual library using ligand-receptor interaction fingerprint-based similarity (IFS). Our approach was evaluated in retrospective virtual screening experiments for antagonists of the metabotropic glutamate receptor (mGluR) subtype 5. The results of our approach were compared to the results obtained by conventional scoring functions (Dock-Score, PMF-Score, Gold-Score, ChemScore, and FlexX-Score). The IFS lead to significantly higher enrichment rates, relative to the competing scoring functions. Though using a target-biased scoring approach, the results were not biased toward the chemical classes of the reference structures. Our results indicate that the presented approach has the potential to serve as a general setup for successful structure-based GPCR virtual screening.     

Use of drug discovery tools in rational organometallic catalyst design

A computational procedure is detailed where techniques common in the drug discovery process-2D- and 3D-quantitative structure-activity relationships (QSAR)-are applied to rationalize the catalytic activity of a synthetically flexible, Ti-N=P ethylene polymerization catalyst system. Once models relating molecular properties to catalyst activity are built with the two QSAR approaches, two database mining approaches are used to select a small number of ligands from a larger database that are likely to produce catalysts with high activity when grafted onto the Ti-N=P framework. The software employed throughout this work is freely available, is easy to use, and was applied in a "black box" approach to highlight areas where the drug discovery tools, designed to address organic molecules, have difficulty in addressing issues arising from the presence of a metal atom. In general, 3D-QSAR offers an efficient way to screen new potential ligands and separate those likely to lead to poor catalysts from those that are likely to contribute to highly active catalysts. The results for 2D-QSAR appear to be quantitatively unreliable, likely due to the presence of a metal atom; nonetheless, there is evidence that qualitative predictions from different models may be reliable. Pitfalls in the database mining techniques are identified, none of which are insurmountable. The lessons learned about the potential uses and drawbacks of the techniques described herein are readily applicable to other catalyst frameworks, thereby enabling a rational approach to catalyst improvement and design.     

Stalis: A Computational Method for Template-Based Ab Initio Ligand Design

Proteins interact with small molecules through specific molecular recognition, which is central to essential biological functions in living systems. Therefore, understanding such interactions is crucial for basic sciences and drug discovery. Here, we present S tructure t emplate-based a b initio li gand design s olution (Stalis), a knowledge-based approach that uses structure templates from the Protein Data Bank libraries of whole ligands and their fragments and generates a set of molecules (virtual ligands) whose structures represent the pocket shape and chemical features of a given target binding site. Our benchmark performance evaluation shows that ligand structure-based virtual screening using virtual ligands from Stalis outperforms a receptor structure-based virtual screening using AutoDock Vina, demonstrating reliable overall screening performance applicable to computational high-throughput screening. However, virtual ligands from Stalis are worse in recognizing active compounds at the small fraction of a rank-ordered list of screened library compounds than crystal ligands, due to the low resolution of the virtual ligand structures. In conclusion, Stalis can facilitate drug discovery research by designing virtual ligands that can be used for fast ligand structure-based virtual screening. Moreover, Stalis provides actual three-dimensional ligand structures that likely bind to a target protein, enabling to gain structural insight into potential ligands. Stalis can be an efficient computational platform for high-throughput ligand design for fundamental biological study and drug discovery research at the proteomic level. © 2019 Wiley Periodicals, Inc.     

Subtype selectivity of dopamine receptor ligands: insights from structure and ligand-based methods

Subtype selective dopamine receptor ligands have long been sought after as therapeutic and/or imaging agents for the treatment and monitoring of neurologic disorders. We report herein on a combined structure- and ligand-based approach to explore the molecular mechanism of the subtype selectivity for a large class of D?-like dopamine receptor ligands (163 ligands in total). Homology models were built for both human D(?L) and D? receptors in complex with haloperidol. Other ligands, which included multiple examples of substituted phenylpiperazines, were aligned against the binding conformations of haloperidol, and three-dimensional quantitative structure activity relationship (3D-QSAR) analyses were carried out. The receptor models show that although D? and D? share highly similar folds and 3D conformations, the slight sequence differences at their extracellular loop regions result in the binding cavity in D? being comparably shallower than in D?, which may explain why some larger ligands bind with greater affinity at D? compared to D? receptors. The QSAR models show excellent correlation and high predictive power even when evaluated by the most stringent criteria. They confirm that the origins of subtype selectivity for the ligands arise primarily due to differences in the contours of the two binding sites. The predictive models suggest that while both steric and electrostatic interactions contribute to the compounds' binding affinity, the major contribution arises from hydrophobic interactions, with hydrogen bonding conferring binding specificity. The current work provides clues for the development of more subtype selective dopamine receptor ligands. Furthermore, it demonstrates the possibility of being able to apply similar modeling methods to other subtypes or classes of receptors to study GPCR receptor-ligand interactions at a molecular level.     

Recent advances in QSAR and their applications in predicting the activities of chemical molecules, peptides and proteins for drug design

This review is to summarize three new QSAR (quantitative structure-activity relationship) methods recently developed in our group and their applications for drug design. Based on more solid theoretical models and advanced mathematical techniques, the conventional QSAR technique has been recast in the following three aspects. (1) In the fragment-based two dimensional QSAR, or abbreviated as FB-QSAR, the molecular structures in a family of drug candidates are divided into several fragments according to the substitutes being investigated. The bioactivities of drug candidates are correlated with physicochemical properties of the molecular fragments through two sets of coefficients: one is for the physicochemical properties and the other for the molecular fragments. (2) In the multiple field three dimensional QSAR, or MF-3D-QSAR, more molecular potential fields are integrated into the comparative molecular field analysis (CoMFA) through two sets of coefficients: one is for the potential fields and the other for the Cartesian three dimensional grid points. (3) In the AABPP (amino acid-based peptide prediction), the bioactivities of peptides or proteins are correlated with the physicochemical properties of all or partial residues of the sequence through two sets of coefficients: one is for the physicochemical properties of amino acids and the other for the weight factors of the residues. Meanwhile, an iterative double least square (IDLS) technique is developed for solving the two sets of coefficients in a training dataset alternately and iteratively. Using the two sets of coefficients, one can predict the bioactivity of a query peptide, protein, or drug candidate. Compared with the old methods, the new QSAR approaches as summarized in this review possess machine learning ability, can remarkably enhance the prediction power, and provide more structural information. Meanwhile, the future challenge and possible development in this area have been briefly addressed as well.