Molecular docking and 3D-QSAR study of pyranmycin derivatives against 16S rRNA A site

To understand pharmacophore properties of pyranmycin derivatives and to design novel inhibitors of 16S rRNA A site, comparative molecular field analysis (CoMFA) approach was applied to analyze three-dimensional quantitative structure–activity relationship (3D-QSAR) of 17 compounds. AutoDock 3.0.5 program was employed to locate the orientations and conformations of the inhibitors interacting with 16S rRNA A site. The interaction mode was demonstrated in the aspects of inhibitor conformation, hydrogen bonding and electrostatic interaction. Similar binding conformations of these inhibitors and good correlations between the calculated binding free energies and experimental biological activities suggest that the binding conformations of these inhibitors derived from docking procedure were reasonable. Robust and predictive 3D-QSAR model was obtained by CoMFA with q² values of 0.723 and 0.993 for cross-validated and non-cross-validated, respectively. The 3D-QSAR model built here will provide clear guidelines for novel inhibitors design based on the Pyranmycin derivatives against 16S rRNA A site.     

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.     

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.     

CoMFA and CoMSIA 3D-quantitative structure-activity relationship model on benzodiazepine derivatives, inhibitors of phosphodiesterase IV

Recently, we reported structurally novel PDE4 inhibitors based on 1,4-benzodiazepine derivatives. The main interest in developing bezodiazepine-based PDE4 inhibitors is in their lack of adverse effects of emesis with respect to rolipram-like compounds. A large effort has thus been made toward the structural optimization of this series. In the absence of structural information on the inhibitor binding mode into the PDE4 active site, 2D-QSAR (H-QSAR) and two 3D-QSAR (CoMFA and CoMSIA) methods were applied to improve our understanding of the molecular mechanism controlling the PDE4 affinity of the benzodiazepine derivatives. As expected, the CoMSIA 3D contour maps have provided more information on the benzodiazepine interaction mode with the PDE4 active site whereas CoMFA has built the best tool for activity prediction. The 2D pharmacophoric model derived from CoMSIA fields is consistent with the crystal structure of the PDE4 active site reported recently. The combination of the 2D and 3D-QSAR models was used not only to predict new compounds from the structural optimization process, but also to screen a large library of bezodiazepine derivatives.     

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.     

Molecular dynamics simulation, free energy calculation and structure-based 3D-QSAR studies of B-RAF kinase inhibitors

(V600E)B-RAF kinase is the most frequent onco-genic protein kinase mutation in melanoma and is a promising target to treat malignant melanoma. In this work, a molecular modeling study combining QM-polarized ligand docking, molecular dynamics, free energy calculation, and three-dimensional quantitative structure-activity relationships (3D-QSAR) was performed on a series of pyridoimidazolone compounds as the inhibitors of (V600E)B-RAF kinase to understand the binding mode between the inhibitors and (V600E)B-RAF kinase and the structural requirement for the inhibiting activity. 3D-QSAR models, including CoMFA and CoMSIA, were developed from the conformations obtained by QM-polarized ligand docking strategy. The obtained models have a good predictive ability in both internal and external validation. Furthermore, molecular dynamics simulation and free energy calculations were employed to determine the detailed binding process and to compare the binding mode of the inhibitors with different activities. The binding free energies calculated by MM/PBSA gave a good correlation with the experimental biological activity. The decomposition of free energies by MM/GBSA indicates the van der Waals interaction is the major driving force for the interaction between the inhibitors and (V600E)B-RAF kinase. The hydrogen bond interactions between the inhibitors with Glu501 and Asp594 of the (V600E)B-RAF kinase help to stabilize the DFG-out conformation. The results from this study can provide some insights into the development of novel potent (V600E)B-RAF kinase inhibitors.     

3D-QSAR and Docking Studies of Pyrido[2,3-d]pyrimidine Derivatives as Wee1 Inhibitors

In order to investigate the inhibiting mechanism and obtain some helpful information for de-signing functional inhibitors against Wee1, three-dimensional quantitative structure-activity relationship (3D-QSAR) and docking studies have been performed on 45 pyrido[2,3-d] pyrim-idine derivatives acting as Wee1 inhibitors. Two optimal 3D-QSAR models with significant statistical quality and satisfactory predictive ability were established, including the CoMFA model (q²=0.707, R²=0.964) and CoMSIA model (q²=0.645, R²=0.972). The external val-idation indicated that both CoMFA and CoMSIA models were quite robust and had high predictive power with the predictive correlation coefficient values of 0.707 and 0.794, essen-tial parameter r²_m values of 0.792 and 0.826, the leave-one-out r²_m(LOO) values of 0.781 and 0.809, r²_m(overall) values of 0.787 and 0.810, respectively. Moreover, the appropriate binding orientations and conformations of these compounds interacting with Wee1 were revealed by the docking studies. Based on the CoMFA and CoMSIA contour maps and docking analyses, several key structural requirements of these compounds responsible for inhibitory activity were identified as follows: simultaneously introducing high electropositive groups to the sub-stituents R1 and R5 may increase the activity, the substituent R2 should be smaller bulky and higher electronegative, moderate-size and strong electron-withdrawing groups for the substituent R3 is advantageous to the activity, but the substituent X should be medium-size and hydrophilic. These theoretical results help to understand the action mechanism and design novel potential Wee1 inhibitors.     

基于计算机模拟技术对JAK2抑制剂的定量构效关系、分子设计和分子动力学研究

本文选取了52个对Janus激酶2(JAK2)有抑制作用的小分子化合物,分别使用3D-QSAR中的CoMFA和CoMSIA方法构建了两个可靠的、具有预测能力的模型,并利用分子对接分析数据集化合物与JAK2蛋白的相互作用,表明化合物主要通过氢键和范德华作用与JAK2靶蛋白结合。根据3D-QSAR模型的分析结果,设计了40个化合物,利用构建的模型预测其抑制活性;使用软件预测了化合物的药代动力学(ADME)参数,开展分子对接模拟,最终选择化合物D01和D22与JAK2靶蛋白进行了分子动力学模拟研究,结果显示两个复合物结合构象稳定,与分子对接结果趋势一致。本研究的结果可以为JAK2抑制剂的研发提供一些新的思路,为临床开发此类药物提供理论支撑。     

含磷嘧啶类CDK9抑制剂的分子对接、3D-QSAR和分子动力学模拟

选取64个具有潜力的含磷嘧啶类细胞周期依赖性蛋白激酶(CDK9)小分子抑制剂,采用分子对接方法研究了该类小分子与CDK9的结合作用,结果表明,分子构象、氢键形成、疏水性和氨基酸残基Cys106在此类抑制剂与CDK9的结合过程中具有重要作用.在配体叠合的基础上,运用比较分子力场分析(Co MFA)、比较分子相似性指数分析(Co MSIA)和Topomer Co MFA(T-COMFA)研究了分子结构与抑制活性的关系,发现由训练集立体场、静电场和疏水场组合的Co MSIA模型为最优模型,其内部交叉验证相关系数(Q2=0.557)、非交叉验证相关系数(R2=0.959)和外部预测相关系数(r2=0.863)具有统计学意义,该模型的三维等值线图直观显示了化合物的活性与其三维结构的关系.根据这些结果设计了10个具有新结构的含磷嘧啶类化合物,分子对接和分子动力学模拟结果表明,新化合物和CDK9的结合模式与原化合物64相同,自由能分析从理论上证明了新化合物64d的CDK9抑制活性优于化合物64,并且显示含磷基团与残基Asp109的静电场能在化合物与CDK9作用过程中有重要作用.     

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.     

DCN1-UBE2M相互作用抑制剂的3D-QSAR、分子对接和分子动力学研究

类泛素化是一种蛋白质翻译后修饰,其异常会导致神经退行性疾病和多种肿瘤的发生,因此它被视为有希望的抗肿瘤靶标。研究表明,抑制DCN1-UBE2M相互作用可选择性阻遏类泛素化。本文基于哌啶基脲类DCN1-UBE2M相互作用抑制剂进行3D-QSAR、分子对接和分子动力学模拟研究。利用3D-QSAR中的CoMFA和CoMSIA方法构建了相关模型,其交叉验证系数q~2分别为0.686、0.682,拟合验证系数r~2分别为0.966、0.931,表明模型是可靠的且预测能力较好。接着运用分子对接分析哌啶基脲类化合物与DCN1的相互作用,结果表明,它们主要通过氢键和疏水作用与靶蛋白结合。通过分子动力学模拟研究进一步了解结合模型和验证对接结果。本文所得的研究结果可为今后此类化合物的结构优化提供有效信息。     

Docking and 3D-QSAR Analysis on a Series of Pyridone-based EZH2 Inhibitors

Enhancer of Zeste homolog 2(EZH2) is closely correlated with malignant tumor and regarded as a promising target to treat B-cell lymphoma. In our research, the molecular docking and three-dimensional quantitative structure-activity relationships(3D-QSAR) studies were performed on a series of pyridone-based EZH2 compounds. Molecular docking allowed us to study the critical interactions at the binding site of EZH2 protein with inhibitors and identify the practical conformations of ligands in binding pocket. Moreover, the docking-based alignment was applied to derive the reliable 3D-QSAR models. Comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA) provided available ability of visualization. All the derived 3D-QSAR models were considered to be statistically significant with respect to the internal and external validation parameters. For the CoMFA model, q~2 = 0.649, r~2 = 0.961 and r~2 pred = 0.877. For the CoMSIA model, q~2 = 0.733, r~2 = 0.980 and r~2 pred = 0.848. With the above arguments, we extracted the correlation between the biological activity and structure. Based on the binding interaction and 3D contour maps, several new potential inhibitors with higher biological activity predicted were designed, which still awaited experimental validation. These theoretical conclusions could be helpful for further research and exploring potential EZH2 inhibitors.     

苯并咪唑类缓蚀剂的3D-QSAR研究及分子设计

采用比较分子场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 对苯并咪唑衍生物抗盐酸腐蚀的缓蚀性能进行了三维定量构效关系研究, 并使用留一法交叉验证手段对3D-QSAR模型的稳定性及预测能力进行了分析. 结果表明, 立体场、静电场和氢键供体场(电子给体)是影响苯并咪唑缓蚀剂缓蚀性能的主要因素; 所构建的CoMFA模型(q2=0.541, R2=0.996)和CoMSIA模型(q2=0.581, R2=0.987)均具有较好的统计学稳定性和预测能力. 基于3D-QSAR等势图设计出了几种具有较好缓蚀性能的苯并咪唑化合物, 为油气田新型缓蚀剂的研发提供了一种新思路.     

芳香噻嗪类衍生物作为选择性醛糖还原酶抑制剂的分子对接和基于受体的三维定量构效关系研究

芳香噻嗪类衍生物被证明是一类选择性较好的高活性醛糖还原酶抑制剂(ARIs).本文对44个芳香噻嗪类化合物进行了分子对接(docking)和三维定量构效关系(3D-QSAR)研究,并探索了此类化合物与醛糖还原酶(ALr²)的作用机理.醛糖还原酶与醛还原酶(ALR1)活性位点的叠加结果显示, ALr²中残基Leu 300和Cys298的存在是化合物1m具有高选择性的原因.分别建立了比较分子场分析方法(CoMFA, q² = 0.649, r² =0.934; q²:交叉验证相关系数, r²:非交叉验证相关系数)和比较分子相似性指数分析方法(CoMSIA, q² = 0.746, r² = 0.971)模型,并对影响此类化合物生物活性的结构进行了鉴定.结果显示,两个模型均具有较高预测能力,并通过测试集中的7个化合物进行了验证,其中CoMFA模型和CoMSIA模型的预测相关系数(r_Pred²)分别为0.748和0.828. 3D-QSAR模型中的三维等值线图表明,在化合物1m的苄基环上C3和C4位置以及苯并噻嗪母核上C5和C7位置进行改进可能对生物活性的提高有利,此预测与我们前期报道的苯并噻嗪母核C7位改进结果一致.本文所建3D-QSAR模型能够在理性设计具有更高生物活性的新型ARIs中发挥重要作用.     

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.     

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.     

6-氮杂吲唑类Pim-1激酶抑制剂的3D-QSAR、药效团模型研究和分子设计

Pim-1 激酶通过作用于多种信号通路或靶点影响肿瘤的发生发展,近年来被认为是肿瘤治疗的良好靶标。本文采用SYBYL-X2. 1. 1软件中的TopomerCoMFA、GALAHAD模块建立计算机模型,研究39个基于6-氮杂吲唑环的Pim-1激酶抑制剂的三维定量构效关系及药效团特征元素。结果显示,TopomerCoMFA建模所得交叉验证系数（q2）和相关系数（r2）分别为0. 756和0. 951,结合外部验证表明此3D-QSAR模型具有较高预测能力及较好的统计学稳定性,同时,用等势图描述了R1、R2基团处立体场、静电场对活性的具体影响。药效团研究结果表明,含氢键受体的芳香杂环母核结构,以及侧链取代基中含有芳香杂环结构对化合物的活性贡献较大。最后根据上述模型信息新设计了15个Pim-1激酶抑制剂分子并完成活性预测及分子对接模式研究,其中4个分子的预测pIC50高于建模分子中活性最好的化合物17,Surflex-Dock分析显示新设计分子均与Pim-1激酶形成较强氢键相互作用。基于6-氮杂吲唑环的Pim-1激酶抑制剂的3D-QSAR模型以及药效团模型可用于指导新型抑制剂的结构优化,为设计和开发具有较高活性的新型Pim-1激酶抑制剂提供有效帮助。