3D QSAR Modeling and Molecular Docking Studies on a Series of Triazole Analogues as Antibacterial Agents

The 3D QSAR analysis using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques is performed on novel nalidixic acid based 1,2,4-triazole derivatives suggested earlier as antibacterial agents. The CoMFA and CoMSIA models employed for a training set of 28 compounds gives reliable values of Q² (0.53 and 0.52, respectively) and R² (0.79 and 0.85, respectively). The contour maps produced by the CoMFA and CoMSIA models are used to determine a three-dimensional quantitative structure-activity relationship. Based on the 3D QSAR contours new molecules with high predicted activities are designed. In addition, surflex-docking is performed to confirm the stability of predicted molecules in the receptor.     

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.     

Combretastatins类微管蛋白抑制剂的定量构效关系与结合模式

以Combretastatins的B环改造化合物为研究对象, 采用遗传函数分析方法进行了二维定量构效关系研究. 研究结果表明, Apol, PMI-mag, Dipole-mag, Hbond donor和RadOfGyration等描述符对该系列抑制剂活性的贡献最大. 采用比较分子场分析方法(CoMFA)和比较分子相似因子分析方法(CoMSIA)进行了三维定量构效关系研究, 建立的CoMFA和CoMSIA模型的交叉验证相关系数q2分别为0.630和0.634, 具有较强的预测能力. 利用CoMFA和CoMSIA模型的三维等势图解析了Combretastatins类化合物的构效关系, 阐明了B环上各取代基对抑制微管蛋白聚合活性的影响, 同时应用分子对接方法分析并验证了定量构效关系模型.     

Modeling ligand-receptor interaction for some MHC class II HLA-DR4 peptide mimetic inhibitors using several molecular docking and 3D QSAR techniques

The ligand-receptor interaction between some peptidomimetic inhibitors and a class II MHC peptide presenting molecule, the HLA-DR4 receptor, was modeled using some three-dimensional (3D) quantitative structure-activity relationship (QSAR) methods such as the Comparative Molecular Field Analysis (CoMFA), Comparative Molecular Similarity Indices Analysis (CoMSIA), and a pharmacophore building method, the Catalyst program. The structures of these peptidomimetic inhibitors were generated theoretically, and the conformations used in the 3D QSAR studies were defined by docking them into the known structure of HLA-DR4 receptor through the GOLD, GLIDE Rigidly, GLIDE Flexible, and Xscore programs. Some of the parameters used in these docking programs were selected by docking an X-ray ligand into the receptor and comparing the root-means-square difference (RMSD) computed between the coordinates of the X-ray and docked structure. However, the goodness of a docking result for docking a series of peptidomimetic inhibitors into the HLA-DR4 receptor was judged by comparing the Spearman's rank correlation coefficient computed between each docking result and the activity data taken from the literature. The best CoMFA and CoMSIA models were constructed using the aligned structures of the best docking result. The CoMSIA was conducted in a stepwise manner to identify some important molecular features that were further employed in a pharmacophore building process by the Catalyst program. It was found that most inhibitors of the training set were accurately predicted by the best pharmacophore model, the Hypo1 hypothesis constructed. The deviation or conflict found between the actual and predicted activities of some inhibitors of both the training and the test sets were also investigated by mapping the Hypo1 hypothesis onto the corresponding structures of the inhibitors.     

3D QSAR selectivity analyses of carbonic anhydrase inhibitors: insights for the design of isozyme selective inhibitors

A 3D QSAR selectivity analysis of carbonic anhydrase (CA) inhibitors using a data set of 87 CA inhibitors is reported. After ligand minimization in the binding pockets of CA I, CA II, and CA IV isoforms, selectivity CoMFA and CoMSIA 3D QSAR models have been derived by taking the affinity differences (DeltapKi) with respect to two CA isozymes as independent variables. Evaluation of the developed 3D QSAR selectivity models allows us to determine amino acids in the respective CA isozymes that possibly play a crucial role for selective inhibition of these isozymes. We further combined the ligand-based 3D QSAR models with the docking program AUTODOCK in order to screen for novel CA inhibitors. Correct binding modes are predicted for various CA inhibitors with respect to known crystal structures. Furthermore, in combination with the developed 3D QSAR models we could successfully estimate the affinity of CA inhibitors even in cases where the applied scoring function failed. This novel strategy to combine AUTODOCK poses with CoMFA/CoMSIA 3D QSAR models can be used as a guideline to assess the relevance of generated binding modes and to accurately predict the binding affinity of newly designed CA inhibitors that could play a crucial role in the treatment of pathologies such as tumors, obesity, or glaucoma.     

3D QSAR studies on protein tyrosine phosphatase 1B inhibitors: comparison of the quality and predictivity among 3D QSAR models obtained from different conformer-based alignments

A set of 65 flexible peptidomimetic competitive inhibitors (52 in the training set and 13 in the test set) of protein tyrosine phosphatase 1B (PTP1B) has been used to compare the quality and predictive power of 3D quantitative structure-activity relationship (QSAR) comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models for the three most commonly used conformer-based alignments, namely, cocrystallized conformer-based alignment (CCBA), docked conformer-based alignment (DCBA), and global minima energy conformer-based alignment (GMCBA). These three conformers of 5-[(2S)-2-({(2S)-2-[(tert-butoxycarbonyl)amino]-3-phenylpropanoyl}amino)3-oxo-3-pentylamino)propyl]-2-(carboxymethoxy)benzoic acid (compound number 66) were obtained from the X-ray structure of its cocrystallized complex with PTP1B (PDB ID: 1JF7), its docking studies, and its global minima by simulated annealing. Among the 3D QSAR models developed using the above three alignments, the CCBA provided the optimal predictive CoMFA model for the training set with cross-validated r2 (q2)=0.708, non-cross-validated r2=0.902, standard error of estimate (s)=0.165, and F=202.553 and the optimal CoMSIA model with q2=0.440, r2=0.799, s=0.192, and F=117.782. These models also showed the best test set prediction for the 13 compounds with predictive r2 values of 0.706 and 0.683, respectively. Though the QSAR models derived using the other two alignments also produced statistically acceptable models in the order DCBA>GMCBA in terms of the values of q2, r2, and predictive r2, they were inferior to the corresponding models derived using CCBA. Thus, the order of preference for the alignment selection for 3D QSAR model development may be CCBA>DCBA>GMCBA, and the information obtained from the CoMFA and CoMSIA contour maps may be useful in designing specific PTP1B inhibitors.     

3D-QSAR studies of dihydropyrazole and dihydropyrrole derivatives as inhibitors of human mitotic kinesin Eg5 based on molecular docking

Human mitotic kinesin Eg5 plays an essential role in mitoses and is an interesting drug target against cancer. To find the correlation between Eg5 and its inhibitors, structure-based 3D-quantitative structure-activity relationship (QSAR) studies were performed on a series of dihydropyrazole and dihydropyrrole derivatives using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods. Based on the LigandFit docking results, predictive 3D-QSAR models were established, with cross-validated coefficient values (q2) up to 0.798 for CoMFA and 0.848 for CoMSIA, respectively. Furthermore, the CoMFA and CoMSIA models were mapped back to the binding sites of Eg5, which could provide a better understanding of vital interactions between the inhibitors and the kinase. Ligands binding in hydrophobic part of the inhibitor-binding pocket were found to be crucial for potent ligand binding and kinases selectivity. The analyses may be used to design more potent EG5 inhibitors and predict their activities prior to synthesis.     

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.     

Three-dimensional Quantitative Structure-Activity Relationship Analyses of a Series of Butenolide ETA Antagonists

Two kinds of Three-dimensional Quantitative Structure-activity Relationship(3D-QSAR) methods,comparative molecular filed analysis(CoMFA) and comparative molecular similarity indices analysis (CoMSIA) ,were applied to analyze the structure-activity relationship of a series of 63 butenolide ETA selective antagonists with respect to their inhibition against human ETA receptor,The CoMFA and CoMSIA models were developed for the conceivable alignment of the molecules based on a template structure from the crystallized data.The statistical results from the initial orientation of the aligned molecules show that the 3D-QSAR model from CoMFA(q^2=0.543) is obviously superior to that from the conventional CoMSIA(q^2=0.407).In order to refine the model,all-space search (ASS) was applied to minimize the field sampling process.By rotating and translating the molecular aggregate within the grid systematically,all the possible samplings of the molecular fields were tested and subsequently the one with the highest q^2 was picked out .The comparison of the sensitivity of CoMFA and CoMSIA to different space orientation shows that the CoMFA q^2 values are more sensitive to the translations and rotations of the aligned molecules with respect to the lattice than those of CoMSIA.The best CoMFA model from ASS was further refined by the region focused technique.The high quality of the best model is indicated by the high corss-validated correlation and the prediction on the external test set.The CoMFA coefficient contour plots identify several key features that explain the wide range of activities,which may help us to design new effective ETA selective antagonists.     

具有除草活性的大环内酯类衍生物的定量构效关系

研究了一系列结构新颖的具有除草活性的大环内酯衍生物的定量构效关系(QSAR). 构建的比较分子力场分析(CoMFA)、比较分子近似指数分析(CoMSIA)和全息定量构效关系(HQSAR)分子模型的交叉验证系数r2cv均大于0.5, 非交叉验证系数r2都超过0.8, 表明获取的QSAR模型具有可信的预测能力. 对CoMFA、CoMSIA模型的三维(3D)等势图分析, 发现除了立体场和静电场外, 疏水场和氢键受体场也是影响大环内酯类化合物除草活性的重要因素. 构建的HQSAR模型的原子贡献图提示的结构改造信息与三维QSAR的结果基本一致. 利用CoMFA、CoMSIA模型提供的信息,对目前已合成的活性最高化合物B1-3进行分子结构改造, 预测结果发现部分化合物可能具有更好的除草活性.     

Molecular modeling of triazine type MDR modulators using CoMFA and CoMSIA approaches

In the present study a series of 30 triazine derivatives was investigated by 3D QSAR methods with respect to their MDR reversing activity in vitro. Two approaches were applied and compared: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Molecular models with good predictive power were derived using steric, electrostatic and hydrophobic fields of the compounds. The results indicated the dominant role of the electrostatic and hydrophobic fields for MDR reversing activity of the investigated modulators. The obtained statistical parameters (Qcv2, Qpr2) showed that the CoMFA and CoMSIA models have similar predictivity. The CoMSIA models were slightly better than the CoMFA ones and obtained with lower number of principal components. The models were graphically interpreted using CoMFA and CoMSIA contour plots. The structural regions responsible for the differences in anti-MDR activity were analyzed in respect to their electrostatic and hydrophobic nature. An easier interpretation of the CoMSIA contour plots was noticed.     

褪黑激素受体拮抗剂的三维定量构效关系研究

采用两种分子场分析方法即比较分子场分析法（CoMFA）和比较分子相似因子分析法（CoMSIA）进行了37个褪黑激素受体拮抗剂的构效关系研究.计算结果表明，两种方法得到的构效关系模型都具有较好的预测能力.在计算中，还考察了不同格点距离和电荷计算方法对构效关系模型的影响.通过分析分子场等值面图在空间的分布，可以观察到叠合分子周围分子场特征对化合物活性的影响，为设计新的褪黑激素拮抗剂提供了一些理论依据.     

Molecular modeling of triazine type MDR modulators using CoMFA and CoMSIA approaches

In the present study a series of 30 triazine derivatives was investigated by 3D QSAR methods with respect to their MDR reversing activity in vitro . Two approaches were applied and compared: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Molecular models with good predictive power were derived using steric, electrostatic and hydrophobic fields of the compounds. The results indicated the dominant role of the electrostatic and hydrophobic fields for MDR reversing activity of the investigated modulators. The obtained statistical parameters ( Q cv 2 , Q pr 2 ) showed that the CoMFA and CoMSIA models have similar predictivity. The CoMSIA models were slightly better than the CoMFA ones and obtained with lower number of principal components. The models were graphically interpreted using CoMFA and CoMSIA contour plots. The structural regions responsible for the differences in anti-MDR activity were analyzed in respect to their electrostatic and hydrophobic nature. An easier interpretation of the CoMSIA contour plots was noticed.     

Use of 3D QSAR to investigate the mode of binding of pyrazinones to HIV-1 RT

Abstract  

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on the docked conformation were performed for 24 pyrazinone derivatives. All compounds were docked into the wild-type HIV-1 RT binding pocket and the lowest-energy docked configurations were used to construct the 3D QSAR models. The CoMFA and CoMSIA models enable good prediction of inhibition by the pyrazinones, with r_\textcv² r_{\text{cv}}^{2}  = 0.703 and 0.735. Results obtained from CoMFA and CoMSIA based on the docking conformation of the pyrazinones are, therefore, powerful means of elucidating the mode of binding of pyrazinones and suggesting the design of new potent NNRTIs.     

Comparing the quality and predictiveness between 3D QSAR models obtained from manual and automated alignment

A set of 113 flexible cyclic urea inhibitors of human immunodeficiency virus protease (HIV-1 PR) was used to compare the quality and predictive power of CoMFA and CoMSIA models for manually or automatically aligned inhibitor set. Inhibitors that were aligned automatically with molecular docking were in agreement with information obtained from existing X-ray structures. Both alignment methods produced statistically significant CoMFA and CoMSIA models, with the best q(2) value being 0.649 and the best predictive r(2) being 0.754. The manual alignment gave statistically higher values, whereas the automated alignment gave more robust models for predicting the activities of an external inhibitor set. Both models utilized similar amino acids in the HIV-1 PR active site, supporting the idea that hydrogen bonds form between an inhibitor and the backbone carbonyl oxygens of Gly48 and Gly48' and also the backbone NH group of Asp30, Gly48, Asp29', and Gly48' of the enzyme. These results suggest that an automated inhibitor alignment can yield predictive 3D QSAR models that are well comparable to manual methods. Thus, an automated alignment method in creating 3D QSAR models is encouragable when a well-characterized structure of the target protein is available.     

Combined 3D-QSAR, molecular docking, and molecular dynamics study on piperazinyl-glutamate-pyridines/pyrimidines as potent P2Y12 antagonists for inhibition of platelet aggregation

An unusually large data set of 397 piperazinyl-glutamate-pyridines/pyrimidines as potent orally bioavailable P2Y(12) antagonists for inhibition of platelet aggregation was studied for the first time based on the combination of three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, and molecular dynamics (MD) methods. The comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) studies have been performed with a training set of 317 compounds, estimating three superimposition methods. The best CoMFA and CoMSIA models, derived from superimposition I, shows leave-one-out cross-validation correlation coefficients (Q(2)) of 0.571 and 0.592 as well as the conventional correlation coefficients (R(2)(ncv)) of 0.814 and 0.834, respectively. In addition, the satisfactory results, based on the bootstrapping analysis and 10-fold cross-validation, further indicate the highly statistical significance of the optimal models. The external predictive abilities of these models were evaluated using a prediction set of 80 compounds, producing the predicted correlation coefficients (R(2)(pred)) of 0.664 and 0.668, respectively. The key amino acid residues were identified by molecular docking, and the stability and rationality of the derived molecular conformations were also validated by MD simulation. The good concordance between the docking results and CoMFA/CoMSIA contour maps provides helpful clues about the rational modification of molecules in order to design more potent P2Y(12) antagonists. We hope the developed models could provide some instructions for further synthesis of highly potent P2Y(12) antagonists.