CoMFA 3D-QSAR analysis of HIV-1 RT nonnucleoside inhibitors, TIBO derivatives based on docking conformation and alignment

HIV-1 RT is one of the key enzymes in the duplication of HIV-1. Inhibitors of HIV-1 RT are classified as nonnucleoside RT inhibitors (NNRTIs) and nucleoside analogues. NNRTIs bind in a region not associated with the active site of the enzyme. Within the NNRTI category, there is a set of inhibitors commonly referred to as TIBO inhibitors. Fifty TIBO inhibitors were used in the work to build 3-D QSAR models. The two known crystal structures of complexes are used to investigate and validate the docking protocol. The results show that the docking simulations reproduce the crystal complexes very well with RMSDs of approximately 1 A and approximately 0.6 A for 1REV and 1COU, respectively. The alignment of molecules and "active" conformation selection are the key to a successful 3D-QSAR model by CoMFA. The flexible docking (Autodock3) was used on determination of "active" conformation and molecular alignment, and CoMFA and CoMSIA were used to develop 3D-QSAR models of 50 TIBOs in the work. The 3D-QSAR models demonstrate a good ability to predict the activity of studied compounds (r2 = 0.972, 0.944, q2 = 0.704, 0.776). It is shown that the steric and electrostatic properties predicted by CoMFA contours can be related to the binding structure of the complex. The results demonstrate that the combination of ligand-based and receptor-based modeling is a powerful approach to build 3D-QSAR models.     

Three-dimensional Quantitative Structure-activity Relationship Models of HIV-1 Integrase Inhibitors of DKAs

As one of the three viral encoded enzymes of HIV-1 infection, HIV-1 integrase has become an attractive drug target for the treatment. Diketoacid compounds (DKAs) are one kind of potent and selective inhibitors of HIV-1 IN. In the present work, two three-dimensional QSAR techniques (CoMFA and CoMSIA) were employed to correlate the molecular structure with the activity of inhibiting the strand transfer for 147 DKAs. The all-oritation search (AOS) and all-placement search (APS) were used to optimize the CoMFA model. The diketo and keto-enol tautomers of DKAs were also used to establish the CoMFA models. The results indicated that the enol was the dominant conformation in the HIV-1 IN and DKAs complexes. It can provide a new method and reference to identify the bioactive conformation of drugs by using QSAR analysis. The best CoMSIA model, with five fields combined, implied that the hydrophobic field is very important as well as the steric and electrostatic fields. All models indicated favorable internal validation. A comparative analysis with the three models demonstrated that the CoMFA model seems to be more predictive. The contour maps could afford steric, electrostatic, hydrophobic and H-bond information about the interaction of ligand-receptor complex visually. The models would give some useful guidelines for designing novel and potent HIV-1 integrase inhibitors.     

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.     

二肽肽酶IV抑制剂的三维定量构效关系研究

二肽肽酶IV是一类用于治疗II型糖尿病具有潜在价值的关键酶, 很多此类酶的抑制剂用于处理此病具有相当好的有效性. 一系列N-取代的甘氨酰氰基吡咯烷衍生物对于二肽肽酶具有高的活性和选择性. 我们使用比较分子力场分析方法建立DPP-IV抑制剂——N-取代的甘氨酰氰基吡咯衍生物的三维定量构效关系, 该模型为设计用于治疗II型糖尿病的高效DPP-IV抑制剂提供结构信息. CoMFA模型的交叉验证相关系数q²＝0.575, 非交叉验证相关系数r²＝0.981, 绝对误差S＝0.184, F_9.68＝388.5. 使用七个预测集检验了模型的预测能力. 所得的模型解释了已有的构效关系, 并对同类化合物有较好的预测能力, 该模型可用于指导新型的DPP-IV抑制剂的设计与优化.     

HIV-1逆转录酶抑制剂的3D-QSAR研究及分子设计

采用Topomer Co MFA方法对24个二芳基苯胺衍生物进行三维定量构效关系研究,建立了3DQSAR模型,所得优化模型的非交叉相关系数、交互验证系数以及外部验证的复相关系数分别为0.928,0.654和0.940,结果表明该模型具有良好的稳定性和预测能力。采用分子对接技术对药物与受体的作用机制进行了研究,结果显示,药物与HIV-1逆转录酶的LYS172,GLU138,LYS101等位点作用明显。运用这些信息进行分子设计,在理论上获得了一些具有较高活性的新的二芳基苯胺类抗艾滋病药物,该QSAR的研究结果可为新药合成提供理论参考。     

Structural exploration of hydroxyethylamines as HIV-1 protease inhibitors: new features identified

The current study deals with chemometric modelling strategies (Naïve Bayes classification, hologram-based quantitative structure–activity relationship (HQSAR), comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA)) to explore the important features of hydroxylamine derivatives for exerting potent human immunodeficiency virus-1 (HIV-1) protease inhibition. Depending on the statistically validated reliable and robust quantitative structure–activity relationship (QSAR) models, important and crucial structural features have been identified that may be responsible for enhancing the activity profile of these hydroxylamine compounds. Arylsulfonamide function along with methoxy or fluoro substitution is important for enhancing activity. Bulky steric substitution at the sulfonamide nitrogen disfavours activity whereas smaller hydrophobic substitution at the same position is found to be favourable. Apart from the crucial oxazolidinone moiety, pyrrolidine, cyclic urea and methyl ester functions are also responsible for increasing the HIV-1 protease inhibitory profile. Observations derived from these modelling studies may be utilized further in designing promising HIV-1 protease inhibitors of this class.     

Three-dimensional quantitative structure activity relationship (QSAR) of cytotoxic active 3,5-diaryl-4,5-dihydropyrazole analogs: a comparative molecular field analysis (CoMFA) revisited study

ABSTRACT: In vitro antitumor evaluation of the synthesized 46 compounds of 3,5-diaryl-4,5-dihydropyrazoles against EAC cell lines and 3D QSAR study using pharmacophore and Comparative Molecular Field Analysis (CoMFA) methods were described. CoMFA derived QSAR model shows a good conventional squared correlation coefficient r2 and cross validated correlation coefficient r2 cv 0.896 and 0.568 respectively. In this analysis steric and electrostatic field contribute to the QSAR equation by 70% and 30% respectively, suggesting that variation in biological activity of the compounds is dominated by differences in steric (van der Waals) interactions. To visualize the CoMFA steric and electrostatic field from partial least squares (PLS) analysis, contour maps are plotted as percentage contribution to the QSAR equation and are associated with the differences in biological activity. BACKGROUND: Pyrazole derivatives exhibit a wide range of biological properties including promising antitumor activity. Furthermore, Aldol condensation assisted organic synthesis has delivered rapid routes to N-containing heterocycles, including pyrazoles. Combining these features, the use of chalconisation-assisted processes will provide rapid access to a targeted dihydropyrazoles library bearing a hydrazino 3D QSAR study using pharmacophore and Comparative Molecular Field Analysis (CoMFA) methods were described for evaluation of antioxidant properties. RESULTS: Chalcones promoted 1 of the 2 steps in a rapid, convergent synthesis of a small library of hydrazinyl pyrazole derivatives, all of which exhibited significant antitumor activity against Ehrlich Ascites Carcinoma (EAC) human tumor cell line comparable to that of the natural anticancer doxorubicin, as a reference standard during this study. In order to understand the observed pharmacological properties, quantitative structure-activity relationship (3D QSAR) study was initiated. CONCLUSIONS: Chalcones heating provides a rapid and expedient route to a series of pyrazoles to investigate their chracterization scavenging properties. Given their favorable properties, in comparison with known anticancer, these pyrazole derivatives are promising leads for further development and optimization.     

氟喹诺酮三唑酮衍生物抗增殖活性的CoMFA模型

应用比较分子力场分析法(CoMFA)研究了18种氟喹诺酮三唑啉酮衍生物对人白血病HL60细胞的体外抗增殖活性(pH)。建立的模型在高相关系数值(R2=0.969)和交叉验证系数值(R2cv=0.473)方面显示出良好的相关和预测能力。空间场和静电场对pH的贡献分别为58.5%和41.5%。基于CoMFA等高线图,揭示了该系列化合物抗增殖活性的一些关键结构因素,如疏水作用、空穴契合、库仑力及氢键等。这些结果为理解其作用机制、设计具有高抗增殖活性的新型氟喹诺酮三唑啉酮类化合物以及预测其活性提供了有益的理论参考。     

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.     

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

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

HEPT类逆转录酶抑制剂的三维定量构效关系

利用比较分子力场分析(CoMFA)方法对32个HEPT类HIV-1逆转录酶抑制剂(RTIs)的三维定量构效关系(3D-QSAR)进行了分析,建立了HIV-1逆转录酶抑制剂的3种3D-QSAR模型,发现影响其生物活性的主要因素为立体场因素,这与HIV-1RT的非底物结合部位(NNBS)的疏水性环境相吻合.进一步分析表明,适当长度的1-位侧链对保持化合物的抗病毒活性致关重要;增大5-位取代基的体积可增强生物活性;在1-位苄氧甲基的对位引入大体积基团有利于提高活性.同时考察立体场、静电场与生物活性的关系,表明,CoMFA模型为最佳预测模型,其交叉验证系数R_CV²=0.870,传统相关系数R²=0.986,标准偏差SE=0.146,F=294.546.用此模型预测了检验组3个HEPT类化合物的-lgEC₅₀,R_pred²=0.850,表明模型具有很好的预测能力,可为HEPT类HIV-1逆转录酶抑制剂的结构优化提供理论指导.     

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.     

苯并呋喃／噻吩联二苯类PTP1B抑制剂三维构效关系研究

主要采用比较分子力场分析方法（CoMFA）对苯并呋喃／噻吩联二苯类PTP1B (protein tyrosine phosphatase 1B)抑制剂进行了三维构效关系的研究，考察了 静电场、立体场和氢键场对构效关系的影响，交叉系数q^2的值达到0．58，表明 CoMFA得到的构效关系模型比较理想，同时test set中分子的预测活性也表明，模 型具有较好的预测能力，研究还表明，氢键场的加入不一定有利于模型的改善，通 过对分子场等值面图的分析，可以观察到叠合分子周围立体场和静电场对化合物活 性的影响，为改进原有化合物的结构，提高它们的活性提供了指导，还尝试采用比 较分子相似性指数分析方法（CoMFA）对这一系列化合物作了研究，结果表明虽然 CoMFA中加入了疏水场，但是对于研究的体系，CoMFA的模型质量并没有显著提高。