Insight into the structural requirements of benzimidazole derivatives as interleukin‐2 inducible t‐cell kinase inhibitors by computational explorations

In the present work, a set of ligand‐ and receptor‐based 3D‐QSAR models were developed to explore the structure–activity relationship of 109 benzimidazole‐based interleukin‐2‐inducible T‐cell kinase (ITK) inhibitors. In order to reveal the requisite 3D structural features impacting the biological activities, a variety of in silico modeling approaches including the comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), docking, and molecular dynamics were applied. The results showed that the ligand‐based CoMFA model (Q² = 0.552, R²_ncv = 0.908, R²_pred = 0.787, SEE = 0.252, SEP = 0.558) and CoMSIA model (Q² = 0.579, R²_ncv = 0.914, R²_pred = 0.893, SEE = 0.240, SEP = 0.538) were superior to other models with greater predictive power. In addition, a combined analysis between the 3D contour maps and docking results showed that: (1) Compounds with bulky or hydrophobic substituents near ring D and electropositive or hydrogen acceptor groups around rings C and D could increase the activity. (2) The key amino acids impacting the receptor–ligand interactions in the binding pocket are Met438, Asp500, Lys391, and Glu439. The results obtained from this work may provide helpful guidelines in design of novel benzimidazole analogs as inhibitors of ITK. © 2013 Wiley Periodicals, Inc.     

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

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

楝酰胺(Rocaglamide)类化合物的三维定量构效关系

采用比较分子力场分析(CoMFA)和比较分子相似因子分析(CoMSIA)方法,对训练集中的26个楝酰胺(Rocaglamide)类化合物进行了三维定量构效关系(3D-QSAR)研究,最终建立的CoMFA模型和CoMSlA模型的q<'2>分别为0.593和0.656.并对测试集中的5个化合物的生物活性进行了预测,结果表明...     

芳基磺酰胺类丙酮酸激酶M2激动剂的基于多复合物的药效团模型和定量构效关系

丙酮酸激酶M2(PKM2)是肿瘤治疗中最具发展潜力的靶点之一. 本文以一系列丙酮酸激酶M2-激动剂复合物的晶体结构为基础, 采用基于多复合物的药效团(MCBP)方法产生了PKM2 的药效团模型. 并使用该药效团模型产生了62个芳基磺酰胺类PKM2激动剂的活性构象和分子叠合, 通过三维定量构效关系(3D-QSAR)方法研究了该类PKM2激动剂与PKM2蛋白的相互作用, 并建立了相关预测模型. 比较分子场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA)模型的交互验证相关系数q²分别为0.545 和0.653, 非交互验证相关系数r²分别为0.966和0.987. 本研究为进一步结构优化、设计和合成新型的PKM2激动剂提供了理论依据.     

3D‐QSAR studies and molecular design on a novel series of pyrimidine benzimidazoles as Lck inhibitors

The development of selective lymphocyte‐specific kinase (Lck) inhibitors has attracted much attention for the research of the treatment of T‐cell mediated autoimmune and inflammatory diseases. In the present work, three‐dimensional quantitative structure–activity relationship (3D‐QSAR) analyses are performed on a novel series of 4‐amino‐6‐benzimidazole‐pyrimidines acting as Lck inhibitors. The established 3D‐QSAR models show significant statistical quality and satisfactory predictive ability, with high q² and R² values: the comparative molecular field analysis (CoMFA) model (q² = 0.802, R² = 0.991), and the comparative molecular similarity indexes analysis (CoMSIA) model (q² = 0.731, R² = 0.982). The systemic external validation indicates that both CoMFA and CoMSIA models are quite robust and possess high predictive abilities with values of 0.881 and 0.877, values of 0.897 and 0.847, values of 0.897 and 0.850, and values of 0.897 and 0.854, respectively. Several key structural features accounting for the inhibitory activities of these compounds are discussed. Based on established models and design considerations, six new compounds with significantly improved activities are theoretically designed, which still await experimental confirmation and evaluation. These theoretical results may provide a useful reference for understanding the action mechanism and designing novel potential Lck inhibitors. © 2014 Wiley Periodicals, Inc.     

螺环吲哚类MDM2抑制剂的分子对接、定量构效关系和分子动力学模拟

采用分子对接、三维定量构效关系(3D-QSAR)和分子动力学方法研究了21个螺环吲哚类化合物与MDM2蛋白的相互作用, 并建立了相关预测模型. 比较分子场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA)模型的交互验证相关系数q²分别为0.573 和0.651, 非交互验证相关系数r²分别为0.948和0.980. 分子对接得到的结合模式与分子动力学模拟得到的结果一致, 结合模式表明该类螺环吲哚化合物主要通过疏水相互作用和氢键与MDM2结合. 基于上述相互作用模型设计并合成了6个新结构螺环吲哚化合物, 并在MDM2高表达的前列腺癌LNCaP细胞株上测定其活性, 结果表明化合物5, 6的半数抑制浓度均低于1μg·mL^-1, 可作为新的抗肿瘤药物先导化合物进一步深入研究. 本研究对以MDM2为靶点的新结构螺环吲哚类抑制剂的开发提供了理论和实验依据.     

五元杂环并嘧啶类胸苷酸合成酶抑制剂的构效关系和分子对接

用比较分子场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA)研究了38个五元杂环并嘧啶衍生物类胸苷酸合成酶抑制剂的三维定量构效关系(3D-QSAR), 建立了相关预测模型. CoMFA和CoMSIA模型的交互验证相关系数q²分别为0.662和0.672、非交互验证相关系数R²分别为0.921和0.884、外部交互验证相关系数Q_ext²分别为0.85和0.81. 分子对接得到的结合模式与三维定量构效关系得到的结果一致. 结果表明这两种模型都具有良好的预测能力, 可应用于指导化合物的设计和结构修饰, 为进一步设计新型胸苷酸合成酶抑制剂提供了理论依据.     

新磺酰脲类化合物除草活性的3D-QSAR分析

用比较分子力场分析 (CoMFA) 方法和比较分子相似性指数分析 (CoMSIA) 方法对所合成的新磺酰脲类化合物的除草活性进行了较为系统的3D-QSAR分析.两种方法所建立的模型对化合物的除草活性预测能力均较好,所得三维等值线图为合成高活性的化合物能提供指导作用     

苯并呋喃类N-肉豆蔻酰基转移酶抑制剂的三维定量构效关系研究

摘要采用比较分子力场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 系统地研究了40个苯并呋喃类N-肉豆蔻酰基转移酶(NMT)抑制剂的三维定量构效关系. 在CoMFA研究中, 考察了网格点步长对模型统计结果的影响. 在CoMSIA研究中, 研究了各种分子场组合、 网格点步长和衰减因子对模型统计结果的影响, 发现立体场、 静电场、 疏水场和氢键受体场的组合可得到最佳模型. 所建立的CoMFA和CoMSIA模型的交叉相关系数q2值分别为0.759和0.730, 均具有较强的预测能力. 利用CoMFA和CoMSIA模型的三维等值线图直观地解释了化合物的构效关系, 阐明了化合物结构中苯并呋喃环上各位置取代基对抑酶活性的影响, 为进一步结构优化提供了重要依据.     

3D-QSAR and molecular modeling of HIV-1 integrase inhibitors

Three-dimensional quantitative structure-activity relationship (3D QSAR) methods were applied on a series of inhibitors of HIV-1 integrase with respect to their inhibition of 3-processing and 3-end joining steps in vitro.The training set consisted of 27 compounds belonging to the class of thiazolothiazepines. The predictive ability of each model was evaluated using test set I consisting of four thiazolothiazepines and test set II comprised of seven compounds belonging to an entirely different structural class of coumarins. Maximum Common Substructure (MCS) based method was used to align the molecules and this was compared with other known methods of alignment. Two methods of 3D QSAR: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were analyzed in terms of their predictive abilities. CoMSIA produced significantly better results for all correlations. The results indicate a strong correlation between the inhibitory activity of these compounds and the steric and electrostatic fields around them. CoMSIA models with considerable internal as well as external predictive ability were obtained. A poor correlation obtained with hydrophobic field indicates that the binding of thiazolothiazepines to HIV-1 integrase is mainly enthalpic in nature. Further the most active compound of the series was docked into the active site using the crystal structure of integrase. The binding site was formed by the amino acid residues 64-67, 116, 148, 151-152, 155-156, and 159. The comparison of coefficient contour maps with the steric and electrostatic properties of the receptor shows high level of compatibility.     

咪唑啉酮类除草剂的三维构效关系研究

从三维角度出发,采取不同的构象搜索方法得到了咪唑酮类化合物分子的活性构象。利用比较分子场分析方法进一步证实了模板分子构象的正确性。并从静电场、立体场及活性关系等方面进行了三维定量构效关系研究,得到了具有较强预测能力的QSAR模型。     

Molecular electrostatic potentials as input for the alignment of HIV-1 integrase inhibitors in 3D QSAR

Comparative molecular similarity indices analysis (CoMSIA), a three-dimensional quantitative structure activity relationship (3D QSAR) paradigm, was used to examine the correlations between the calculated physicochemical properties and the in vitro activities (3'-processing and 3'-strand transfer inhibition) of a series of human immunodeficiency virus type 1 (HIV-1) integrase inhibitors. The training set consisted of 34 molecules from five structurally diverse classes: salicylpyrazolinones, dioxepinones, coumarins, quinones, and benzoic hydrazides. The data set was aligned using extrema of molecular electrostatic potentials (MEPs). The predictive ability of the resultant model was evaluated using a test set comprised of 7 molecules belonging to a different structural class of thiazepinediones. A CoMSIA model using an MEP-based alignment showed considerable internal as well external predictive ability (r2(cv) = 0.821, r2(pred) = 0.608 for 3'-processing; and r2(cv) = 0.759, r2(pred.) = 0.660 for 3'-strand transfer).     

Pharmacophore modelling,virtual screening and molecular docking studies on PLD1 inhibitors

Lipid metabolism plays a significant role in influenza virus replication and subsequent infection. The regulatory mechanism governing lipid metabolism and viral replication is not properly understood to date, but both Phospholipase D (PLD1 and PLD2) activities are stimulated in viral infection. In vitro studies indicate that chemical inhibition of PLD1 delays viral entry and reduction of viral loads. The current study reports a three-dimensional pharmacophore model based on 35 known PLD1 inhibitors. A sub-set of 25 compounds was selected as the training set and the remaining 10 compounds were kept in the test set. One hundred and twelve pharmacophore models were generated; a six-featured pharmacophore model (AADDHR.57) with survival score (2.69) produced a statistically significant three-dimensional quantitative structure–activity relationship model with r² = 0.97 (internal training set), r² = 0.71 (internal test set) and Q² = 0.64. The predictive power of the pharmacophore model was validated with an external test set (r² = 0.73) and a systematic virtual screening work-flow was employed showing an enrichment factor of 23.68 at the top 2% of the dataset (active and decoys). Finally, the model was used for screening of the filtered PubChem database to fetch molecules which can be proposed as potential PLD1 inhibitors for blocking influenza infection.     

Comparative Molecular Similarity Index Analysis (CoMSIA) to study hydrogen-bonding properties and to score combinatorial libraries

Comparative molecular field analysis has been applied to a data set of thermolysin inhibitors. Fields expressed in terms of molecular similarity indices (CoMSIA) have been used instead of the usually applied Lennard-Jones- and Coulomb-type potentials (CoMFA). Five different properties, assumed to cover the major contributions responsible for ligand binding, have been considered: steric, electrostatic, hydrophobic, and hydrogen-bond donor or acceptor properties. The statistical evaluation of the field properties by PLS analysis reveals a similar predictive potential to CoMFA. However, significantly improved and easily interpretable contour maps are obtained. The features in these maps intuitively suggest where to modify a molecular structure in terms of physicochemical properties and functional groups in order to improve its binding affinity. They can also be interpreted with respect to the known structural protein environment of thermolysin. Most of the highlighted regions in the maps are mirrored by features in the surrounding environment required for binding. Using the derived correlation model, different members of a combinatorial library designed for thermolysin inhibition have been scored for affinity. The results obtained demonstrate the prediction power of the CoMSIA method.     

基于分子对接的6-萘甲基取代HEPT类HIV-1逆转录酶抑制剂构效关系研究

采用分子对接方法得到了一系列6-萘甲基取代HEPT类逆转录酶抑制剂分子与HIV-1逆转录酶复合物模型,从中抽取出抑制剂分子的活性构象,进一步应用CoMFA和CoMSIA方法建立了具有较好预测能力的3D-QSAR模型,深入探讨了这些化合物的定量构效关系,为进一步的药物设计奠定了良好的基础.另外,以化合物13及其相应的β异构体24为代表,结合量子化学从头算分子轨道理论方法考察了它们的前线轨道,为阐明α和β系列化合物的活性差异提供了理论依据.     

Three-dimensional quantitative structure-activity relationship (3D QSAR) and pharmacophore elucidation of tetrahydropyran derivatives as serotonin and norepinephrine transporter inhibitors

Three-dimensional quantitative structure-activity relationship (3D QSAR) using comparative molecular field analysis (CoMFA) was performed on a series of substituted tetrahydropyran (THP) derivatives possessing serotonin (SERT) and norepinephrine (NET) transporter inhibitory activities. The study aimed to rationalize the potency of these inhibitors for SERT and NET as well as the observed selectivity differences for NET over SERT. The dataset consisted of 29 molecules, of which 23 molecules were used as the training set for deriving CoMFA models for SERT and NET uptake inhibitory activities. Superimpositions were performed using atom-based fitting and 3-point pharmacophore-based alignment. Two charge calculation methods, Gasteiger-Hückel and semiempirical PM3, were tried. Both alignment methods were analyzed in terms of their predictive abilities and produced comparable results with high internal and external predictivities. The models obtained using the 3-point pharmacophore-based alignment outperformed the models with atom-based fitting in terms of relevant statistics and interpretability of the generated contour maps. Steric fields dominated electrostatic fields in terms of contribution. The selectivity analysis (NET over SERT), though yielded models with good internal predictivity, showed very poor external test set predictions. The analysis was repeated with 24 molecules after systematically excluding so-called outliers (5 out of 29) from the model derivation process. The resulting CoMFA model using the atom-based fitting exhibited good statistics and was able to explain most of the selectivity (NET over SERT)-discriminating factors. The presence of −OH substituent on the THP ring was found to be one of the most important factors governing the NET selectivity over SERT. Thus, a 4-point NET-selective pharmacophore, after introducing this newly found H-bond donor/acceptor feature in addition to the initial 3-point pharmacophore, was proposed.     

含磷嘧啶类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作用过程中有重要作用.