Molecular Modeling and Design of Arylthioindole Derivatives as Tubulin Inhibitors

Three-dimensional quantitative structure activity relationship （3D-QSAR） and docking studies of a series of arylthioindole derivatives as tubulin inhibitors against human breast cancer cell line MCF-7 have been carried out. An optimal 3D-QSAR model from the comparative molecular field analysis （CoMFA） for training set with significant statistical quality （R2=0.898） and predictive ability （q2=0.654） was established. The same model was further applied to predict pIC50 values of the compounds in test set, and the resulting predictive correlation coefficient R2（pred） reaches 0.816, further showing that this CoMFA model has high predictive ability. Moreover, the appropriate binding orientations and conformations of these compounds interacting with tubulin are located by docking study, and it is very interesting to find the consistency between the CoMFA field distribution and the 3D topology structure of active site of tubulin. Based on CoMFA along with docking results, some important factors improving the activities of these compounds were discussed in detail and were summarized as follows： the substituents R3-R5 （on the phenyl ring） with higher electronegativity, the substituent R6 with higher eleetropositivity and bigger bulk, the substituent R7 with smaller bulk, and so on. In addition, five new compounds with higher activities have been designed. Such results can offer useful theoretical references for experimental works.     

3D QSAR and molecular docking studies of benzimidazole derivatives as hepatitis C virus NS5B polymerase inhibitors

The urgent need for novel HCV antiviral agents has provided an impetus for understanding the structural requisites of NS5B polymerase inhibitors at the molecular level. Toward this objective, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of 67 HCV NS5B polymerase inhibitors were performed using two methods. First, ligand-based 3D QSAR studies were performed based on the lowest energy conformations employing the atom fit alignment method. Second, receptor-based 3D QSAR models were derived from the predicted binding conformations obtained by docking all NS5B inhibitors at the allosteric binding site of NS5B (PDB ID: 2dxs). Results generated from the ligand-based model were found superior (r2cv values of 0.630 for CoMFA and 0.668 for CoMSIA) to those obtained by the receptor-based model (r2cv values of 0.536 and 0.561 for CoMFA and CoMSIA, respectively). The predictive ability of the models was validated using a structurally diversified test set of 22 compounds that had not been included in a preliminary training set of 45 compounds. The predictive r2 values for the ligand-based CoMFA and CoMSIA models were 0.734 and 0.800, respectively, while the corresponding predictive r2 values for the receptor-based CoMFA and CoMSIA models were 0.538 and 0.639, respectively. The greater potency of the tryptophan derivatives over that of the tyrosine derivatives was interpreted based on CoMFA steric and electrostatic contour maps. The CoMSIA results revealed that for a NS5B inhibitor to have appreciable inhibitory activity it requires hydrogen bond donor and acceptor groups at the 5-position of the indole ring and an R substituent at the chiral carbon, respectively. Interpretation of the CoMFA and CoMSIA contour maps in context of the topology of the allosteric binding site of NS5B provided insight into NS5B-inhibitor interactions. Taken together, the present 3D QSAR models were found to accurately predict the HCV NS5B polymerase inhibitory activity of structurally diverse test set compounds and to yield reliable clues for further optimization of the benzimidazole derivatives in the data set.     

Molecular modeling studies of human immunodeficiency virus type 1 protease inhibitors using three‐dimensional quantitative structure‐activity relationship,virtual screening,and docking simulations

In this paper, two 3‐dimensional quantitative structure‐activity relationship models for 60 human immunodeficiency virus (HIV)‐1 protease inhibitors were established using random sampling analysis on molecular surface and translocation comparative molecular field vector analysis (Topomer CoMFA). The non–cross‐validation (r²), cross‐validation (q²), correlation coefficient of external validation (Q²_ext), and F of 2 models were 0.94, 0.80, 0.79, and 198.84 and 0.94, 0.72, 0.75, and 208.53, respectively. The results indicated that 2 models were reasonable and had good prediction ability. Topomer Search was used to search R groups in the ZINC database, 20 new compounds were designed, and the Topomer CoMFA model was used to predicate the biological activity. The results showed that 18 new compounds were more active than the template molecule. So the Topomer Search is effective in screening and can guide the design of new HIV/AIDS drugs. The mechanism of action was studied by molecular docking, and it showed that the protease inhibitors and Ile50, Asp25, and Arg8 sites of HIV‐1 protease have interactions. These results have provided an insight for the design of new potent inhibitors of HIV‐1 protease.     

Mapping the binding site of a large set of quinazoline type EGF-R inhibitors using molecular field analyses and molecular docking studies

In the current work, three-dimensional QSAR studies for one large set of quinazoline type epidermal growth factor receptor (EGF-R) inhibitors were conducted using two types of molecular field analysis techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). These compounds belonging to six different structural classes were randomly divided into a training set of 122 compounds and a test set of 13 compounds. The statistical results showed that the 3D-QSAR models derived from CoMFA were superior to those generated from CoMSIA. The most optimal CoMFA model after region focusing bears significant cross-validated r(2)(cv) of 0.60 and conventional r(2) of 0.92. The predictive power of the best CoMFA model was further validated by the accurate estimation to these compounds in the external test set, and the mean agreement of experimental and predicted log(IC(50)) values of the inhibitors is 0.6 log unit. Separate CoMFA models were conducted to evaluate the influence of different partial charges (Gasteiger-Marsili, Gasteiger-Hückel, MMFF94, ESP-AM1, and MPA-AM1) on the statistical quality of the models. The resulting CoMFA field map provides information on the geometry of the binding site cavity and the relative weights of various properties in different site pockets for each of the substrates considered. Moreover, in the current work, we applied MD simulations combined with MM/PBSA (Molecular mechanics/Possion-Boltzmann Surface Area) to determine the correct binding mode of the best inhibitor for which no ligand-protein crystal structure was present. To proceed, we define the following procedure: three hundred picosecond molecular dynamics simulations were first performed for the four binding modes suggested by DOCK 4.0 and manual docking, and then MM/PBSA was carried out for the collected snapshots. The most favorable binding mode identified by MM/PBSA has a binding free energy about 10 kcal/mol more favorable than the second best one. The most favorable binding mode identified by MM/PBSA can give satisfactory explanation of the SAR data of the studied molecules and is in good agreement with the contour maps of CoMFA. The most favorable binding mode suggests that with the quinazoline-based inhibitor, the N3 atom is hydrogen-bonded to a water molecule which, in turn, interacts with Thr 766, not Thr 830 as proposed by Wissner et al. (J. Med. Chem. 2000, 43, 3244). The predicted complex structure of quinazoline type inhibitor with EGF-R as well as the pharmacophore mapping from CoMFA can interpret the structure activities of the inhibitors well and afford us important information for structure-based drug design.     

Homology model-guided 3D-QSAR studies of HIV-1 integrase inhibitors

In the present study, we report the exploration of binding modes of potent HIV-1 integrase (IN) inhibitors MK-0518 (raltegravir) and GS-9137 (elvitegravir) as well as chalcone and related amide IN inhibitors we recently synthesized and the development of 3D-QSAR models for integrase inhibition. Homology models of DNA-bound HIV-1 IN were constructed on the basis of the X-ray crystal structure of the foamy virus IN-DNA complex (PDB ID: 3L2T ) and used for docking. The binding modes of raltegravir and elvitegravir in our homology models are in accordance with those in the foamy virus structure revealing interactions important for inhibitor-IN binding. To gain further insights into the structural requirements for IN inhibition, three-dimensional quantitative structure activity relationship (3D-QSAR) studies were conducted using raltegravir, elvitegravir, and their analogs; our synthesized 3-keto salicylic acid IN inhibitor series; as well as other structurally related HIV-1 IN inhibitors. In the first part of the study with 103 compounds, atom-fit alignments, I and II, and docking-based alignment, III, were used to develop 3D-QSAR models 1, 2, and 3, respectively, each comprising comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) 3D-QSARs. This initial analysis indicated that the docking-based (structure-based) model 3 performed better than the atom-fit (ligand-based) models 1 and 2, in terms of statistical significance and robustness. Thus, the docking-based alignment was then subsequently used with an expanded data set of 296 compounds for building a more comprehensive 3D-QSAR, model 4. Model 4 afforded good q2 values of 0.70 and 0.75 for CoMFA and CoMSIA 3D-QSARs, respectively, and showed good predictive performance on an external validation test set of 59 compounds with predictive r2 values up to 0.71. The HIV IN-DNA homology model of biological relevance and the comprehensive 3D-QSAR models developed in the present study provide insights and new predictive tools for structure-based design and optimization of IN inhibitors.     

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of thiazolone derivatives as hepatitis C virus NS5B polymerase allosteric inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models for a series of thiazolone derivatives as novel inhibitors bound to the allosteric site of hepatitis C virus (HCV) NS5B polymerase were developed based on CoMFA and CoMSIA analyses. Two different conformations of the template molecule and the combinations of different CoMSIA field/fields were considered to build predictive CoMFA and CoMSIA models. The CoMFA and CoMSIA models with best predictive ability were obtained by the use of the template conformation from X-ray crystal structures. The best CoMFA and CoMSIA models gave q (2) values of 0.621 and 0.685, and r (2) values of 0.950 and 0.940, respectively for the 51 compounds in the training set. The predictive ability of the two models was also validated by using a test set of 16 compounds which gave r (pred) (2) values of 0.685 and 0.822, respectively. The information obtained from the CoMFA and CoMSIA 3D contour maps enables the interpretation of their structure-activity relationship and was also used to the design of several new inhibitors with improved activity.     

Receptor-independent 4D-QSAR analysis of a set of norstatine derived inhibitors of HIV-1 protease

A training set of 27 norstatine derived inhibitors of HIV-1 protease, based on the 3(S)-amino-2(S)-hydroxyl-4-phenylbutanoic acid core (AHPBA), for which the -log IC(50) values were measured, was used to construct 4D-QSAR models. Five unique RI-4D-QSAR models, from two different alignments, were identified (q(2) = 0.86-0.95). These five models were used to map the atom type morphology of the lining of the inhibitor binding site at the HIV protease receptor site as well as predict the inhibition potencies of seven test set compounds for model validation. The five models, overall, predict the -log IC(50) activity values for the test set compounds in a manner consistent with their q(2) values. The models also correctly identify the hydrophobic nature of the HIV protease receptor site, and inferences are made as to further structural modifications to improve the potency of the AHPBA inhibitors of HIV protease. The finding of five unique, and nearly statistically equivalent, RI-4D-QSAR models for the training set demonstrates that there can be more than one way to fit structure-activity data even within a given QSAR methodology. This set of unique, equally good individual models is referred to as the manifold model.     

3D-QSAR and molecular docking study of LRRK2 kinase inhibitors by CoMFA and CoMSIA methods

Three-dimensional quantitative structure–activity relationship (3D-QSAR) modelling was conducted on a series of leucine-rich repeat kinase 2 (LRRK2) antagonists using CoMFA and CoMSIA methods. The data set, which consisted of 37 molecules, was divided into training and test subsets by using a hierarchical clustering method. Both CoMFA and CoMSIA models were derived using a training set on the basis of the common substructure-based alignment. The optimum PLS model built by CoMFA and CoMSIA provided satisfactory statistical results (q² = 0.589 and r² = 0.927 and q² = 0.473 and r² = 0.802, respectively). The external predictive ability of the models was evaluated by using seven compounds. Moreover, an external evaluation set with known experimental data was used to evaluate the external predictive ability of the porposed models. The statistical parameters indicated that CoMFA (after region focusing) has high predictive ability in comparison with standard CoMFA and CoMSIA models. Molecular docking was also performed on the most active compound to investigate the existence of interactions between the most active inhibitor and the LRRK2 receptor. Based on the obtained results and CoMFA contour maps, some features were introduced to provide useful insights for designing novel and potent LRRK2 inhibitors.     

HIV-1整合酶链转移抑制剂的3D-QSAR、分子对接和分子动力学模拟研究

为了获得高活性、结构新颖的整合酶链转移(INST)抑制剂,本文采用Co MFA和Co MSIA两种方法对32个萘啶类INST抑制剂进行了三维定量构效关系研究,并建立了相关模型,其交叉验证系数分别为q~2=0. 809和q~2=0. 816,拟合验证系数分别为r~2=0. 998和r~2=0. 981,表明所建立的模型是可靠的且具有一定的预测能力。利用分子对接探讨小分子化合物与INST蛋白的相互作用模式,结果表明,萘啶类化合物主要通过疏水作用和氢键作用与INSTIs蛋白结合。最后通过分子动力学模拟进一步验证对接结果发现,对接的结合模式与分子动力学模拟得到的结果是一致的。本研究获得的综合模型和推论可以为开发有效的HIV INSTIs提供重要的理论信息。     

Ab initio calculation of the interaction energy in the P2 binding pocket of HIV‐1 protease

Interactions at the P2 binding pocket of human immunodeficiency virus type 1 (HIV‐1) protease have been studied using calculated interaction energies for model systems that mimic this binding pocket. Models were built for the P2 pocket of HIV‐1 protease in complex with TMC114, nelfinavir, and amprenavir. A two‐step procedure was applied. In the first step, the size of the model system was confined to ～40 atoms, and the interaction energy was calculated at different computational levels. In the second step, the size of the system was increased to 138 atoms, and the calculations were only performed at the HF/6‐31G** level. The interaction energy of the HIV‐1 protease/TMC114 complex was found to be more favorable than the interaction energies of the other complexes because of the additional hydrogen bond interaction this inhibitor is able to make with the HIV‐1 protease backbone. The results of the calculations are supported by stockholder charges and electrostatic potential maps. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

3C-like蛋白酶抑制剂的构效关系、分子对接和分子动力学

3C-like蛋白酶是中东呼吸综合征冠状病毒（MERS-CoV）等其它冠状病毒的繁殖过程中极为重要的蛋白酶。它已成为人类在抗冠状病毒领域中的研究热点。本文基于计算生物学方法对与MERS-CoV同属的蝙蝠冠状病毒HKU4（HKU4-CoV）的43个肽类3C-like蛋白酶抑制剂分子,建立三维定量构效关系（3D-QSAR）模型。在基于配体叠合的基础上,发现比较分子相似性指数分析法（CoMSIA）中的四个场组合（位阻场、静电场、氢键供体场与氢键受体场）为最优的模型（Q²=0.522,R_ncv²=0.996,R_pre²=0.904;Q²：交叉验证相关系数,R_ncv²：非交叉验证相关系数,R_pre²：验证集分子的预测值相关系数）,并借助该模型通过分子对接（docking）与分子动力学（MD）方法阐明了配受体结合作用。实验结果表明：（1）基于最优的CoMSIA模型基础上的三维等势图形象地说明了分子基团的位阻作用、静电作用、氢键供体与氢键受体作用对分子生物活性的影响;（2）分子对接研究结果显示了疏水性以及结晶水、氨基酸His166和Glu169在配体和受体结合过程中产生重要作用;（3）分子动力学模拟进一步验证了分子对接模型的可靠性,并发现了两个新的关键氨基酸Ser24与Gln192,它们与配体产生了两个较强的氢键。此外,根据这些结果,一些新的具有潜在抑制活性的肽类化合物作为3C-like蛋白酶抑制剂被获得。以上结果能够帮助深入了解3C-like蛋白酶与肽类抑制剂的作用机理,并且能够为今后的抗MERS-CoV药物设计提供有价值的参考。     

Structure-based CoMFA and CoMSIA study of indolinone inhibitors of PDK1

Phosphoinositide-dependent protein kinase-1 (PDK1) is a Ser/Thr kinase which phosphorylates and activates members of the AGC kinase group known to control processes such as tumor cell growth, protection from apoptosis, and tumor angiogenesis. In this paper, CoMFA and CoMSIA studies were carried out on a training set of 56 conformationally rigid indolinone inhibitors of PDK1. Predictive 3D QSAR models, established using atom fit alignment rule based on crystallographic-bound conformation, had cross-validated (r _cv²) values of 0.738 and 0.816 and non-cross-validated (r _ncv²) values of 0.912 and 0.949 for CoMFA and CoMSIA models, respectively. The predictive ability of the CoMFA and CoMSIA models was determined using a test set of 14 compounds, which gave predictive correlation coefficients (r _pred²) of 0.865 and 0.837, respectively. Structure-based interpretation of the CoMFA and CoMSIA field properties provided further insights for the rational design of new PDK1 inhibitors. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A ligand-based comparative molecular field analysis (CoMFA) and homology model based molecular docking studies on 3′, 4′-dihydroxyflavones as rat 5-lipoxygenase inhibitors: Design of new inhibitors

In this study, ligand based comparative molecular field analysis (CoMFA) with five principal components was performed on class of 3′, 4′-dihydroxyflavone derivatives for potent rat 5-LOX inhibitors. The percentage contributions in building of CoMFA model were 91.36% for steric field and 8.6% for electrostatic field. R² values for training and test sets were found to be 0.9320 and 0.8259, respectively. In case of LOO, LTO and LMO cross validation test, q² values were 0.6587, 0.6479 and 0.5547, respectively. These results indicate that the model has high statistical reliability and good predictive power. The extracted contour maps were used to identify the important regions where the modification was necessary to design a new molecule with improved activity. The study has developed a homology model for rat 5-LOX and recognized the key residues at the binding site. Docking of most active molecule to the binding site of 5-LOX confirmed the stability and rationality of CoMFA model. Based on molecular docking results and CoMFA contour plots, new inhibitors with higher activity with respect to the most active compound in data set were designed.     

Insight into the structural requirements of narlaprevir-type inhibitors of NS3/NS4A protease based on HQSAR and molecular field analyses

In the life cycle of hepatitis C virus (HCV), NS3/NS4A protease has been proved to play a vital role in the replication of the HCV virus. Narlaprevir and its derivatives, the inhibitors of NS3/NS4A, would be potentially developed as important anti-HCV drugs in the future. In this study, quantitative structure-activity relationship (QSAR) analyses for 190 narlaprevir derivatives were conducted using comparative molecular field analysis (CoMFA), comparative molecular indices analysis (CoMSIA) and hologram quantitative structure-activity relationship (HQSAR) techniques. Both of the best CoMFA and HQSAR models showed statistical significance for the training set and good predictive accuracy for the test set, which strongly manifested the robustness of the CoMFA and HQSAR models. The CoMFA contour maps and the HQSAR contribution maps were both presented. Furthermore, based on the essential factors for ligand binding derived from the QSAR models, sixteen new derivatives were designed and some of them showed higher inhibitory activities confirmed by our models and molecular docking studies. General speaking, this study provides useful suggestions for the design of potential anti-HCV drugs.