5-羧基苯并咪唑类HCV NS5B聚合酶抑制剂的HQSAR研究及分子设计

利用分子全息技术研究了129个5-羧基苯并咪唑类HCV NS5B聚合酶抑制剂的结构与活性之间的关系.讨论了分子碎片大小、碎片区分参数及全息长度对模型质量的影响.利用偏最小二乘法(partial least square,PLS)建立了一组以99个化合物为训练集的最优模型,该模型的交叉验证相关系数q~2=0.820,非交叉验证相关系数r~2=0.963,标准偏差SEE=0.213;用最优模型对由30个化合物组成的测试集进行预测,得到其相关系数r_(pred)~2=0.98,表明了该模型具有良好的预测能力及拟合能力.利用色码图对模型中不同原子及不同结构的贡献进行了解释,在此基础上根据最优HQSAR模型设计了几种具有良好抗HCV活性的苯并咪唑类HCV NS5B聚合酶抑制剂分子,为新型HCV NS5B聚合酶抑制剂的设计和优化提供了参考.     

3-吡啶基醚类化合物的分子全息QSAR研究

采用分子全息定量构效关系(HQSAR, hologram quantitative structure-activity relationship)方法, 研究了28个3-吡啶基醚类化合物对乙酰胆碱α4β2受体的亲和性与它们的分子结构之间的关系, 讨论了分子碎片大小、分子碎片亚结构类型以及分子全息长度对QSAR的影响, 得到了较好的HQSAR模型, 模型的交叉验证系数平方q²＝0.670, 非交叉相关系数平方r²＝0.965, 偏差S＝0.093. 利用HQSAR的颜色编码, 对化合物中不同基团对亲和活性的影响进行了讨论, 对新配体的合成具有一定的指导作用.     

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

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

苯并咪唑类缓蚀剂的HQSAR研究及分子设计

采用分子全息定量构效关系(HQSAR)方法, 构建苯并咪唑衍生物在酸性环境中的缓蚀性能与结构之间的定量构效关系模型, 研究不同碎片区分参数及碎片大小对模型质量的影响, 寻找最优HQSAR模型, 并对其稳定性及预测能力进行评价. 结果显示: 选取碎片区分参数为原子类型(A)、化学键类型(B)、连接性(C)、氢原子(H)、手性(Ch)、氢键给体和受体(D&;A), 碎片大小为1-3 建模时, 得到的HQSAR 模型(r²(非交叉验证系数)=0.996, q²(交叉验证系数)=0.960, SE_cv(交叉验证标准误差)=3.709)具有良好的统计学稳定性及预测能力. 根据最优HQSAR模型图设计出的38种苯并咪唑类化合物理论上均具有较好的缓蚀性能. 本研究为油气田新型高效缓蚀剂研发提供可靠的理论依据.     

三苯胺类染料敏化太阳能电池能量转化效率的全息定量构效关系研究

研究构建了三苯胺类化合物分子结构与相应染料敏化太阳能电池能量转化效率(PCE)之间的全息定量构效关系(HQSAR)模型。当fragment distinction、fragment size、hologram length和principal components分别为“C、DA”、“4-7”、“199”和“6”时,可以获得最优HQSAR模型。采用外部测试集验证和留一交叉验证对所建立模型进行检验,外部测试集验证中CCC和 Q2F3分别为0.933和0.892,留一交叉验证中其q2cv和r2分别为0.791和0.902,表明所建立模型具有较好的拟合效果和预测能力。通过所建立HQSAR模型的分子贡献图可知,环戊并二噻吩基团的存在有利于提高PCE值,长烷基链的存在可能降低PCE值。     

基于5-芳基乙内酰脲类化合物的CoMFA和HQSAR研究

用CoMFA和HQSAR两种QSAR方法研究了50个乙内酰脲类分子的定量构效关系.本研究从构象搜索所得的低能结构出发构建化合物分子的构象, 建立CoMFA模型,并进行了全空间搜索. HQSAR本质上是一种二维的QSAR方法,与CoMFA方法相比,该方法在数据处理方面,比CoMFA方法快捷,并且可重复性好.两种方法均得到了较好分析结果, CoMFA的交叉验证相关系数q2 值为0.815, HQSAR的q2值为0.893.这些方程有力地说明了该类分子在(R,R)-N-3,5-dinitrobenzoyl-1,2-diamine型手性固定相上拆分过程中的影响因素,对今后类似拆分的实验研究提供了理论支持.     

6-异位-5,8-O-二甲基乙酰紫草素衍生物的2D/3D-QSAR研究

本文综合运用密度泛函理论(DFT)、分子力学(MM2)和统计学等方法对22个具有抗人体乳腺癌(MDA-MB-231)活性的6-异位-5,8-O-二甲基乙酰紫草素类衍生物进行二维(2D)定量构效关系(QSAR)研究,同时运用比较分子相似性指数分析(Co MSIA)方法进行三维(3D)QSAR研究。所建最优2D-QSAR方程的留一法交叉验证系数(q2)和拟合相关系数(R2)分别为0.833和0.900;Co MSIA(SEHA)模型的q2和非交叉验证系数(r2)分别为0.624和0.999,预测相关系数R2pred为0.838,表明所建立的2D/3D-QSAR模型都具有良好的统计学意义及合理、可信的预报能力,可以预测未知化合物的活性。研究结果为理解该类化合物的作用机理和设计合成更高活性的新化合物提供理论参考。     

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

取代喹啉类化合物抗菌活性的3D-QSAR研究及分子设计

通过比较分子力场分析方法(Co MFA)研究取代喹啉类化合物对金黄色葡萄球菌抑菌活性(p M)的三维定量结构-活性相关(3D-QSAR)。12个化合物建立了预测模型,7个化合物作为验证集(含模板分子)。训练集的Co MFA模型显示立体场、静电场对生物活性贡献依次为49.8%、50.2%。该模型的交叉验证相关系数R2cv=0.650,非交叉验证相关系数R2=0.918,对测试集中的7个化合物的生物活性进行了预测,显示出较强的稳定性和良好的预测能力。通过分析Co MFA三维等势图发现,在取代喹啉类化合物抑菌机理中,R4取代基的强吸电性起主要作用,其次是其他取代基的疏水性作用。应用上述规律进行分子设计,获得了3个在理论上具有较高抑菌活性的新的取代喹啉衍生物,期待实验的验证。     

1,2-萘醌类化合物抑制PTP1B的三维定量构效关系研究

蛋白酪氨酸磷酸酶1B (protein tyrosine phosphatase 1B, PTP-1B)是近年来发现的治疗II型糖尿病的新靶点, 1,2-萘醌类化合物对PTP-1B有较好的抑制活性, 具有良好的药用前景. 为了设计出本类化合物抑制效果更好的分子构型, 用比较分子力场分析(CoMFA)和比较分子相似性指数分析(CoMSIA)对该类化合物进行了三维定量构效关系(3D-QSAR)的研究, 并建立了相关的预测模型. 其中, CoMFA模型的交叉验证相关系数(q²)为0.555, 非交叉验证相关系数(r²)为0.991, 标准偏差(SEE)为0.049, F值为564.910. CoMSIA模型的q²为0.558, r²为0.991, SEE为0.050, F值为542.773. 计算结果表明, 获得的CoMFA和CoMSIA模型具有良好的预测能力, 可以应用于指导该类化合物的设计.     

QSAR and Pharmacophore Studies of Thiazolidine-4-carboxylic Acid Derivatives as Novel Influenza Neuraminidase Inhibitors Using HQSAR, Topomer CoMFA and CoMSIA

In order to understand the chemical-biological interactions governing their activities toward neuraminidase (NA), QSAR models of 28 thiazolidine-4-carboxylic acid derivatives with inhibitory influenza A virus were developed. The obtained HQSAR (hologram quantitative structure activity relationship), Topomer CoMFA and CoMSIA (comparative molecular similarity indices analysis) models were robust and had good exterior predictive capabilities. Moreover, QSAR modeling results elucidated that hydrogen bonds highly contributed to the inhibitory activity, then electrostatic and hydrophobic factors. Squared multiple correlation coefficients (R2) of HQSAR, Topomer CoMFA and CoMSIA models were 0.994, 0.978 and 0.996, respectively. Squared cross-validated correlation coefficients (Q2) of HQSAR, Topomer CoMFA and CoMSIA models were in turn 0.951, 919 and 0.820. Furthermore, squared multiple correlation coefficients for the test set (R2test) of HQSAR, CoMFA and CoMSIA models were 0.879, 0.912 and 0.953, respectively. Squared cross-validated correlation coefficients for the test set (Q2ext) of HQSAR, Topomer CoMFA and CoMSIA models were 0.867, 0.884 and 0.899, correspondingly.     

Topomer CoMFA,HQSAR Studies and Molecular Docking of 2,5-Diketopiperazine Derivatives as Oxytocin Inhibitors

Topomer comparative molecular field analysis(Topomer Co MFA) and holographic quantitative structure-activity relationship(HQSAR) for 130 2,5-diketopiperazine derivatives were used to build a three-dimensional quantitative structure-activity relationship(3D-QSAR) model. The results show that the models have high predictive ability. For Topomer CoMFA, the cross-validated q~2 value is 0.710 and the non-cross-validated r~2 value is 0.834. The most effective HQSAR model shows that the cross-validation q~2 value is 0.700, the non-cross-validated r~2 value is 0.815, and the best hologram length value is 353 using connections and bonds as fragment distinctions. 50 highly active 2,5-diketopiperazine derivatives were designed based on the three-dimensional equipotential map and HQSAR color code map. Finally, the molecular docking method was also used to study the interactions of these new molecules by docking the ligands into the diketopiperazine active site, which revealed the likely bioactive conformations. This study showed that there are extensive interactions between the new molecule and Arg156, Arg122 residues in the active site of diketopiperazine. These results provide useful insights for the design of potent of the new 2,5-diketopiperazine derivatives.     

Predicting toxicity of benzene derivatives by molecular hologram derived quantitative structure-activity relationships (QSARS)

Holographic quantitative structure-activity relationship (HQSAR) is an emerging QSAR technique with the combined application of molecular hologram, which encodes the frequency of occurrence of various molecular fragment types, and the subsequent partial least squares (PLS) regression analysis. Based on molecular hologram, alignment-free QSAR models could be rapidly and easily developed with highly statistical significance and predictive ability. In this paper, the toxicity data for a series of 83 benzene derivatives to the autotrophic Chlorella vulgaris (IGC50, negative logarithmic form of 6-h 50% population growth inhibition concentration in mmol/l) were subjected to HQSAR analysis and this resulted in a model with a high predictive ability. The robustness and predictive ability of the model were validated by "leave-one-out" (LOO) cross-validation procedure and an external testing set. The influence of fragment distinction parameters and fragment size on the quality of the HQSAR model have been also discussed.     

A comparative study of quantitative structure-activity relationship methods based on gallic acid derivatives

By using hologram quantitative structure-activity relationship (HQSAR) and comparative molecular field analysis (CoMFA) methods, the relationships between the structures of 49 gallic acid derivatives and their analgesic activity have been investigated to yield statistically reliable models with considerable predictive power. The best HQSAR model was generated using atoms, bond and connectivity as fragment distinction parameters and fragment size 5-7 from a hologram length of 307 with 3 components. High conventional r2 (r2 = 0.825) and cross-validation r2 (r2(cv) = 0.726) values were obtained. CoMFA analyses varying lattice size and location, grid spacing, probe charges and using, Tripos standard and Indicator force field were performed. The best model was developed with 4 components using sp3-hybridized carbon atom with +1.0 charge as probe, grid spacing (2 A), lattice offset (1.0, 3.0, -2.5). The CoMFA model showed a conventional correlation coefficient r2 of 0.889 and across-validation r2(cv) equals to 0.633. The robustness and predictive ability of the HQSAR and CoMFA models have been validated by means of an external test set. The results indicate that both models possess high statistical quality in the prediction of analgesic potency of novel gallic acid analogs.     

Predicting toxicity of benzene derivatives by molecular hologram derived quantitative structure–activity relationships (QSARS)

Holographic quantitative structure–activity relationship (HQSAR) is an emerging QSAR technique with the combined application of molecular hologram, which encodes the frequency of occurrence of various molecular fragment types, and the subsequent partial least squares (PLS) regression analysis. Based on molecular hologram, alignment-free QSAR models could be rapidly and easily developed with highly statistical significance and predictive ability. In this paper, the toxicity data for a series of 83 benzene derivatives to the autotrophic Chlorella vulgaris (IGC⁵⁰, negative logarithmic form of 6-h 50% population growth inhibition concentration in mmol/l) were subjected to HQSAR analysis and this resulted in a model with a high predictive ability. The robustness and predictive ability of the model were validated by “leave-one-out” (LOO) cross-validation procedure and an external testing set. The influence of fragment distinction parameters and fragment size on the quality of the HQSAR model have been also discussed.