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.     

Microscopic binding of butyrylcholinesterase with quinazolinimine derivatives and the structure–activity correlation

Butyrylcholinesterase (BChE) is not only an important protein for development of anti-cocaine medication but also an established drug target to develop new treatment for Alzheimer’s disease (AD). The molecular basis of interaction of a new series of quinazolinimine derivatives as BChE inhibitors has been studied by molecular docking and molecular dynamics (MD) simulations. The molecular docking and MD simulations revealed that all of these inhibitors bind with BChE in similar binding mode. Based on the similar binding mode, we have carried out three-dimensional quantitative structure–activity relationship (3D-QSAR) studies on these inhibitors using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), to understand the structure–activity correlation of this series of inhibitors and to develop predictive models that could be used in the design of new inhibitors of BChE. The study has resulted in satisfactory 3D-QSAR models. We have also developed ligand-based 3D-QSAR models. The contour maps obtained from the 3D-QSAR models in combination with the simulated binding structures help to better interpret the structure–activity relationship and is consistent with available experimental activity data. The satisfactory 3D-QSAR models strongly suggest that the determined BChE-inhibitor binding modes are reasonable. The identified binding modes and developed 3D-QSAR models for these BChE inhibitors are expected to be valuable for rational design of new BChE inhibitors that may be valuable in the treatment of Alzheimer’s disease.     

Combined 3D-QSAR modeling and molecular docking study on indolinone derivatives as inhibitors of 3-phosphoinositide-dependent protein kinase-1

3-Phosphoinositide-dependent protein kinase-1 (PDK1) is a promising target for developing novel anticancer drugs. In order to understand the structure-activity correlation of indolinone-based PDK1 inhibitors, we have carried out a combined molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) modeling study. The study has resulted in two types of satisfactory 3D-QSAR models, including the CoMFA model (r(2)=0.907; q(2)=0.737) and CoMSIA model (r(2)=0.991; q(2)=0.824), for predicting the biological activity of new compounds. The detailed microscopic structures of PDK1 binding with inhibitors have been studied by molecular docking. We have also developed docking-based 3D-QSAR models (CoMFA with q(2)=0.729; CoMSIA with q(2)=0.79). The contour maps obtained from the 3D-QSAR models in combination with the docked binding structures help to better interpret the structure-activity relationship. All of the structural insights obtained from both the 3D-QSAR contour maps and molecular docking are consistent with the available experimental activity data. This is the first report on 3D-QSAR modeling of PDK1 inhibitors. The satisfactory results strongly suggest that the developed 3D-QSAR models and the obtained PDK1-inhibitor binding structures are reasonable for the prediction of the activity of new inhibitors and in future drug design.     

Virtual Screening and Structure Generation Applied to Drug Design

The methods of computer-aided drug design can be divided into two categories according to whether or not the structures of receptors are known1, corresponding to two principal strategies: (1) searching the bio-active ligands against virtual combinatorial libraries and calculating the affinity energy between ligand and receptor by docking ; (2) QSAR and 3D-structure data-mining. 3D-QSAR method is now applied widely to drug discovery, but this method is generally limited to refine the structu…     

Integration of Ligand-Based and Structure-Based Methods for the Design of Small-Molecule TLR7 Antagonists

Toll-like receptor 7 (TLR7) is activated in response to the binding of single-stranded RNA. Its over-activation has been implicated in several autoimmune disorders, and thus, it is an established therapeutic target in such circumstances. TLR7 small-molecule antagonists are not yet available for therapeutic use. We conducted a ligand-based drug design of new TLR7 antagonists through a concerted effort encompassing 2D-QSAR, 3D-QSAR, and pharmacophore modelling of 54 reported TLR7 antagonists. The developed 2D-QSAR model depicted an excellent correlation coefficient (R²_training: 0.86 and R²_test: 0.78) between the experimental and estimated activities. The ligand-based drug design approach utilizing the 3D-QSAR model (R²_training: 0.95 and R²_test: 0.84) demonstrated a significant contribution of electrostatic potential and steric fields towards the TLR7 antagonism. This consolidated approach, along with a pharmacophore model with high correlation (R_training: 0.94 and R_test: 0.92), was used to design quinazoline-core-based hTLR7 antagonists. Subsequently, the newly designed molecules were subjected to molecular docking onto the previously proposed binding model and a molecular dynamics study for a better understanding of their binding pattern. The toxicity profiles and drug-likeness characteristics of the designed compounds were evaluated with in silico ADMET predictions. This ligand-based study contributes towards a better understanding of lead optimization and the future development of potent TLR7 antagonists.     

Fragment-based strategy for structural optimization in combination with 3D-QSAR

Fragment-based drug design has emerged as an important methodology for lead discovery and drug design. Different with other studies focused on fragment library design and active fragment identification, a fragment-based strategy was developed in combination with three-dimensional quantitative structure–activity relationship (3D-QSAR) for structural optimization in this study. Based on a validated scaffold or fragment hit, a series of structural optimization was conducted to convert it to lead compounds, including 3D-QSAR modelling, active site analysis, fragment-based structural optimization and evaluation of new molecules. 3D-QSAR models and active site analysis provided sufficient information for confirming the SAR and pharmacophoric features for fragments. This strategy was evaluated through the structural optimization on a c-Met inhibitor scaffold 5H-benzo[4,5]cyclohepta[1,2-b]pyridin-5-one, which resulted in an c-Met inhibitor with high inhibitory activity. Our study suggested the effectiveness of this fragment-based strategy and the druggability of our newly explored active region. The reliability of this strategy indicated it could also be applied to facilitate lead optimization of other targets.     

In vitro antitumor activity,ADME-Tox and 3D-QSAR of synthesized and selected natural styryl lactones

Prostate cancer is a common cause of death in men and a novel treating methods should be developed. In order to find a new drug for prostate cancer, a series of novel conformationally constrained analogues of (+)-goniofufurone and 7-epi-(+)-goniofufurone, as well as the newly synthesized styryl lactones containing the cinnamic acid ester groups were evaluated for in vitro cytotoxicity against prostate cancer cell (PC-3). Furthermore, prediction of physicochemical characteristics and drugability as well as in silico ADME-Tox tests of investigated compounds were performed. The 3D-QSAR model was established using the comparative molecular field analysis method. According to obtained results, the tricyclic compounds 9 and 10 had the highest potency with IC₅₀ < 20 μM. Evaluation of structural features through 3D-QSAR model identified steric field feature on the cinnamic acid ester groups at C-7 as a crucial for the cytotoxic activity. This research suggests that most of the analysed compounds have desirable properties for drug candidates and high potential in drug development, which recommend them for further research in treatment of prostate cancer. Furthermore, obtained 3D-QSAR model is able to successfully identify styryl lactones that have significant cytotoxic activity and provide information for screening and design of novel inhibitors against PC-3 cell line that could be used as drugs in treatment of the prostate cancer.     

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.     

3D-QSAR and receptor modeling of tyrosine kinase inhibitors with flexible atom receptor model (FLARM)

A set of epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors was investigated with the aim of developing 3D-QSAR models using the Flexible Atom Receptor Model (FLARM) method. Some 3D-QSAR models were built with high correlation coefficients, and the FLARM method predicted the biological activities of compounds in test set well. The FLARM method also gave the pseudoreceptor model, which indicates the possible interactions between the receptor and the ligand. The possible interactions include two hydrogen bonds, one hydrophobic interaction, and one sulfur-aromatic interaction, which are in accord with those in the pharmacophore model given by the scientists at Novartis. This shows that the FLARM method can bridge 3D-QSAR and receptor modeling in computer-aided drug design. Pharmacophore can be obtained according to these results, and 3D searching can then be done with databases to find the lead compound of EGFR tyrosine kinase inhibitors.     

A 3D-QSAR-driven approach to binding mode and affinity prediction

A method for predicting the binding mode of a series of ligands is proposed. The procedure relies on three-dimensional quantitative structure-activity relationships (3D-QSAR) and does not require structural knowledge of the binding site. Candidate alignments are automatically built and ranked according to a consensus scoring function. 3D-QSAR analysis based on the selected binding mode enables affinity prediction of new drug candidates having less than 10 rotatable bonds.     

A selectivity study of benzenesulfonamide derivatives on human carbonic anhydrase II/IX by 3D-QSAR,Molecular Docking and Molecular Dynamics Simulation

Nowadays, different approaches have been pursued with the intent to develop sulfonamide-like carbonic anhydrase inhibitors that possess better selectivity profiles toward the different human isoforms of the enzyme. Here, we used conventional 3D-QSAR methods, including comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), and Topomer CoMFA, to construct three-dimensional quantitative structure-activity relationship (3D-QSAR) models for benzenesulfonamide derivatives as human carbonic anhydrase (hCA) II/IX inhibitors. The theoretical models had good reliability (R²>0.75) and predictability (Q²>0.55), and the contour maps could graphically present the contributions of the force fields for activity and identify the structural divergence between human carbonic anhydrase II inhibitors and human carbonic anhydrase IX inhibitors. Consequently, we explored the selectivity of inhibitor for human carbonic anhydrase II and IX through molecular docking, and the difference of activity coincides with the potential binding mode well. According to the results of the predicted values and the molecule docking, we found that the inhibitors published in the literature had stronger inhibition on the hCA IX; based on the theoretical models, we designed seven new compounds with good potential activity and reasonably good ADMET profile, which could selectively inhibit hCA IX. Molecular Dynamics Simulation showed that newly-designed compound D7 had good selectivity on hCA IX. The findings from 3D-QSAR and docking studies maybe helpful in the rational drug design of isoform-selective inhibitors.     

STAT3环烷烃并噻吩类衍生物抑制剂的三维定量构效关系研究

本文对STAT3抑制剂的化学结构与生物活性之间的关系进行研究。采用三维定量构效关系(3D-QSAR)中的比较分子力场分析(CoMFA)和比较分子相似性指数分析(CoMSIA)方法针对52个STAT3抑制剂建立3D-QSAR模型,阐明了抑制剂化学结构与其生物活性之间的关系。所构建的CoMFA模型交叉验证系数为0.548,非交叉验证系数为0.754,标准偏差为0.278,显著系数为58.297;所构建的CoMSIA模型交叉验证系数为0.892,非交叉验证系数为0.597,标准偏差为0.192,显著系数为57.794。结果显示CoMFA和CoMSIA模型具有良好的稳定性和预测能力。3D-QSAR模型等势图提供的相关场信息对新型STAT3抑制剂的设计具有指导意义。     

A three-dimensional chemical phase pharmacophore mapping,QSAR modelling and electronic feature analysis of benzofuran salicylic acid derivatives as LYP inhibitors

Lymphoid tyrosine phosphatase (LYP), encoded by the PTPN22 gene, has a critical negative regulatory role in T-cell antigen receptor (TCR) and emerged as a promising drug target for human autoimmune diseases. A five-point pharmacophore with two hydrogen bond acceptors, one hydrogen bond donor and two aromatic ring features was generated for a series of benzofuran salicylic acid derivatives as LYP inhibitors in order to elucidate their anti-autoimmune activity. The generated pharmacophore yielded a significant 3D-QSAR model with r² of 0.9146 for a training set of 27 compounds. The model also showed excellent predictive power with Q² of 0.7068 for a test set of eight compounds. The investigation of the 3D-QSAR model has revealed the structural insights which could lead to more potent analogues. The most active and inactive compounds were further subjected to electronic structure analysis using density functional theory (DFT) at B3LYP/3?21^?G level to support the 3D-QSAR predictions. The results obtained from this study are expected to be useful in the proficient design and development of benzofuran salicylic acid derivatives as inhibitors of LYP.     

GPR40 受体苯丙酸类激动剂三维定量构效关系研究

苯丙酸类化合物是G蛋白偶联受体40(GPR40)潜在的生物活性药物。本文基于比较分子力场分析法(Co MFA)和比较分子相似性指数分析法(CoMSIA),分别建立了40个已知活性的GPR40受体苯丙酸类激动剂的三维定量构效关系(3D-QSAR)模型,研究该类激动剂与生物活性之间的关系。CoMFA和CoMSIA模型的交叉验证系数(q~2)分别为0. 527和0. 500,拟合验证系数(r~2)分别为0. 901和0. 860,两个3D-QSAR模型预测值与实验值基本一致,表明模型具有良好的可信度和预测能力。根据两个3D-QSAR模型提供的立体场、静电场、疏水场、氢键供体场和氢键受体场所提供的信息提出优化该类抑制剂结构的药物设计思路,为指导设计更高活性的GPR40激动剂以及GRR40新分子激动活性的预测提供理论依据。     

Virtual screening and rational drug design method using structure generation system based on 3D-QSAR and docking

An efficient virtual and rational drug design method is presented. It combines virtual bioactive compound generation with 3D-QSAR model and docking. Using this method, it is possible to generate a lot of highly diverse molecules and find virtual active lead compounds. The method was validated by the study of a set of anti-tumor drugs. With the constraints of pharmacophore obtained by DISCO implemented in SYBYL 6.8, 97 virtual bioactive compounds were generated, and their anti-tumor activities were predicted by CoMFA. Eight structures with high activity were selected and screened by the 3D-QSAR model. The most active generated structure was further investigated by modifying its structure in order to increase the activity. A comparative docking study with telomeric receptor was carried out, and the results showed that the generated structures could form more stable complexes with receptor than the reference compound selected from experimental data. This investigation showed that the proposed method was a feasible way for rational drug design with high screening efficiency.     

Virtual screening and rational drug design method using structure generation system based on 3D-QSAR and docking

An efficient virtual and rational drug design method is presented. It combines virtual bioactive compound generation with 3D-QSAR model and docking. Using this method, it is possible to generate a lot of highly diverse molecules and find virtual active lead compounds. The method was validated by the study of a set of anti-tumor drugs. With the constraints of pharmacophore obtained by DISCO implemented in SYBYL 6.8, 97 virtual bioactive compounds were generated, and their anti-tumor activities were predicted by CoMFA. Eight structures with high activity were selected and screened by the 3D-QSAR model. The most active generated structure was further investigated by modifying its structure in order to increase the activity. A comparative docking study with telomeric receptor was carried out, and the results showed that the generated structures could form more stable complexes with receptor than the reference compound selected from experimental data. This investigation showed that the proposed method was a feasible way for rational drug design with high screening efficiency.     

Exploration of a binding mode of benzothiazol-2-yl acetonitrile pyrimidine core based derivatives as potent c-Jun N-terminal kinase-3 inhibitors and 3D-QSAR analyses

C-Jun N-terminal kinase (JNK) is a therapeutic target for inhibitors which may provide clinical benefit in the pathogenesis of rheumatoid arthritis (RA) as well as in various apoptosis-related disorders. The benzothiazol-2-yl acetonitrile derivatives, recently reported by Pascale et al. (J. Med. Chem. 2005, 48, 4596-4607), are the first generation JNK inhibitors of this class. To understand inhibitory mechanisms and elucidate pharmacophoric properties of these derivatives molecular docking and 3D-QSAR studies were performed on a set of 44 compounds. Ligand Fit module of Cerius2 (4.9) was employed to locate the binding orientations of all the compounds within the JNK-3 ATP binding site. A good correlation (r2=0.810) between the calculated binding free energies (-PMF score) and the experimental inhibitory activities suggests that the identified binding conformations of these potential inhibitors are reliable. Based on the binding conformations, robust and highly predictive 3D-QSAR models were developed with conventional r2 0.886 and 0.802, full cross-validation r2 0.980 and 0.788, and predictive r2 0.965 and 0.968 for MFA and MSA, respectively. The interaction mode was demonstrated taking into consideration inhibitor conformation, hydrogen bonding, and electrostatic interaction. The 3D-QSAR model built in this study will provide clear guidelines for a novel inhibitor design based on the benzothiazole derivatives against JNK-3 for the treatment of inflammatory disorders.     

Receptor-based modeling and 3D-QSAR for a quantitative production of the butyrylcholinesterase inhibitors based on genetic algorithm

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models have been constructed using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) for a series of structurally related steroidal alkaloids as butyrylcholinesterase (BuChE) inhibitors. Docking studies were employed to position the inhibitors into the BuChE active site to determine the most probable binding mode. The strategy was to explore multiple inhibitor conformations in producing a more reliable 3D-QSAR model. These multiple conformations were derived using the FlexS program. The conformation selection step for CoMFA was done by genetic algorithm. The genetic algorithm based CoMFA approach was found to be the best. Both CoMFA and CoMSIA yielded significant cross-validated q(2) values of 0.701 and 0.627 and the r(2) values of 0.979 and 0.982, respectively. These statistically significant models were validated by a test set of five compounds. Comparison of CoMFA and CoMSIA contour maps helped to identify structural requirements for the inhibitors and serves as a basis for the design of the next generation of the inhibitor analogues. The results demonstrate that the combination of ligand-based and receptor-based modeling with use of a genetic algorithm is a powerful approach to build 3D-QSAR models. These data can be used for the lead optimization process with respect to inhibition enhancement which is important for the drug discovery and development for Alzheimer's disease.     

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