顺式新烟碱类衍生物的三维定量构效关系研究

采用比较分子力场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA)对34个顺式新烟碱类衍生物的杀虫活性进行三维定量构效关系(3D-QSAR)研究.构建的CoMFA和CoMSIA模型的交叉验证系数rc2v分别为0.877和0.862,非交叉验证系数r2分别为0.970和0.961,表明建立的3D-QSAR模型具有较好的统计相关性和预测能力.一系列的研究结果指出:立体场、静电场和氢键受体场是描述顺式新烟碱类衍生物的化学结构与杀虫活性关系的重要参数;在咪唑啉环的3,4位不宜引入较大的取代基,提高咪唑啉环的电负性或增强硝基一个端氧的氢键受体特征有利于提高顺式新烟碱类衍生物的杀虫活性.     

毒蕈碱受体激动剂的三维定量构效关系研究

采用比较分子场分析法(CoMFA)研究了55个四氢吡啶类毒蕈碱受体激动剂的三维定量构效关系(3D-QSAR), 建立了具有较强预测能力的3D-QSAR模型. 所得模型的交叉验证相关系数(q²)为0.507, 常规相关系数(R²)为0.982 , 标准方差为0.218, 说明系列化合物分子周围立体场和静电场的分布与生物活性间存在良好的相关性. 模型不仅很好地预测了训练集和测试集化合物的活性, 而且为设计活性更高的受体激动剂提供了理论依据.     

含呋喃环双酰脲类衍生物的三维定量构效关系研究

采用比较分子力场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 对27个新型双酰基脲类化合物的杀蚊幼虫(Aedes aegypti L.)活性进行三维定量构效关系(3D-QSAR)研究. 在CoMFA研究中, 考察了网格点步长对统计结果的影响. 在CoMSIA研究中, 系统考察了各种分子场组合、网格点步长和衰减因子对模型统计结果的影响, 发现立体场和氢键供体场的组合得到最佳模型. 所建立的CoMFA和CoMSIA模型的非交叉验证相关系数r2值分别为0.828和0.841, 并都具有较强的预测能力. CoMFA和CoMSIA模型的三维等值图不仅直观地解释了结构与活性的关系, 而且为后续优化该系列化合物提供了理论依据.     

2,4-二氨基嘧啶类黏着斑激酶抑制剂的3D-QSAR研究

本文选取42个2,4-二氨基嘧啶类黏着斑激酶(FAK)小分子抑制剂,分别以比较分子场分析法(CoMFA)与相似性指数分析法(CoMSIA)构建3D-QSAR模型,评价模板分子、公共骨架点、最小能量优化参数、分子构象等因素对模型优化的影响。分析最优模型中立体场、静电场以及氢键等因素对抑制剂活性的影响,并应用分子对接分析该类抑制剂与FAK蛋白的相互作用。结果表明,选择化合物16作为模板分子,骨架A作为公共骨架点,最小能量优化参数中电荷、最大迭代系数、最低能量限定值分别为MMFF94、1000、0.01kcal/mol时所构建的模型最优。以CoMFA和CoMSIA构建的3D-QSAR模型的交叉验证系数(q~2)分别为0.666和0.736,非交叉验证系数(R~2)分别为0.990和0.989,表明此模型具有良好的预测能力。分子对接分析显示,其与FAK的氨基酸残基CYS502、ASP564形成了重要的氢键作用,并与周围残基形成了较强的疏水作用。通过3D-QSAR的构建与分子对接分析,可指导2,4-二氨基嘧啶类FAK小分子抑制剂的进一步结构优化设计。     

广谱沙粒病毒抑制剂的3D-QSAR分析

沙粒病毒(Arenaviruses)遍布全球,其中的拉沙热病毒可引起致命的拉沙热.通过应用比较分子场分析(CoMFA)和比较相似性指数分析法(CoMSIA)对47个广谱沙粒病毒抑制剂进行了三维定量构效关系(3D-QSAR)分析.使用立体场、静电场、疏水场和氢键受体场组合获得最优模型CoMSIA的统计结果为Q2=0.518,R2ncv=0.972,R2pre=0.911,说明该模型的可靠性和较好预测能力.此外,模型等势线图直观地解释了分子结构与其活性的关系,为进一步设计新型高效的沙粒病毒抑制剂提供了理论依据.     

STAT3环烷烃并噻吩类衍生物抑制剂的三维定量构效关系研究

本文对STAT3抑制剂的化学结构与生物活性之间的关系进行研究。采用三维定量构效关系(3D-QSAR)中的比较分子力场分析(CoMFA)和比较分子相似性指数分析(CoMSIA)方法针对52个STAT3抑制剂建立3D-QSAR模型,阐明了抑制剂化学结构与其生物活性之间的关系。所构建的CoMFA模型交叉验证系数为0.548,非交叉验证系数为0.754,标准偏差为0.278,显著系数为58.297;所构建的CoMSIA模型交叉验证系数为0.892,非交叉验证系数为0.597,标准偏差为0.192,显著系数为57.794。结果显示CoMFA和CoMSIA模型具有良好的稳定性和预测能力。3D-QSAR模型等势图提供的相关场信息对新型STAT3抑制剂的设计具有指导意义。     

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

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

新型三唑类抗真菌化合物的三维定量构效关系研究

采用比较分子力场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA), 系统研究了40个新型三唑类化合物抗真菌活性的三维定量构效关系. 在CoMFA研究中, 研究了两种药效构象对模型的影响, 并考察了网格点步长对统计结果的影响. 在CoMSIA研究中, 系统考察了各种分子场组合、网格点步长和衰减因子对模型统计结果的影响, 发现立体场、静电场、疏水场和氢键受体场的组合得到最佳模型. 所建立CoMFA和CoMSIA模型的交叉相关系数q²值分别为0.718和0.655, 并都具有较强的预测能力. CoMFA和CoMSIA模型的三维等值线图直观地解释了化合物的构效关系, 阐明了化合物结构中苯环上各位置取代基对抗真菌活性的影响, 为进一步结构优化提供了重要依据.     

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

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

2,4-二氨基嘧啶类FAK抑制剂的3D-QSAR研究

摘要 本文选取42个2,4-二氨基嘧啶类FAK小分子抑制剂,分别以比较分子场分析法（CoMFA）与相似性指数分析法（CoMSIA）构建3D-QSAR模型,评价模板分子、公共骨架点、最小能量优化参数、分子构象等因素对模型优化的影响。分析最优模型中立体场、静电场以及氢键等因素对抑制剂活性的影响,并应用分子对接分析该类抑制剂与FAK蛋白的相互作用。结果表明选择16号化合物作为模板分子,骨架A作为公共骨架点,最小能量优化参数中电荷、最大迭代系数、最低能量限定值分别为MMFF94、1000、0.01 Kcal/mol时所构建的模型最优。以CoMFA和CoMSIA构建的3D-QSAR模型的交叉验证系数（q2）分别为0.666和0.736,非交叉验证系数（R2）分别为0.990和0.989,表明此模型具有良好的预测能力。分子对接分析显示,其与FAK的氨基酸残基CYS502、ASP564形成了重要的氢键作用,并与周围残基形成了较强的疏水作用。通过3D-QSAR的构建与分子对接分析,可指导2,4-二氨基嘧啶类FAK小分子抑制剂的进一步结构优化设计。     

Molecular modeling studies on 11H-dibenz[b,e]azepine and dibenz[b,f][1,4]oxazepine derivatives as potent agonists of the human TRPA1 receptor

A computational strategy based on comparative molecular fields analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) was performed on a series of the 11H-dibenz[b,e]azepine and dibenz[b,f][1,4]oxazepine derivatives as potent agonists of the human TRPA1 receptor. The CoMFA and CoMSIA models resulting from a 21 molecule training set gave r2(cv) values of 0.631 and 0.542 and r2 values of 0.986 and 0.981, respectively. The statistically significant models were validated by a test set of five compounds with predictive r2(pred). values of 0.967 and 0.981 for CoMFA and CoMSIA, respectively. A systemic external validation was also performed on the established models. The information obtained from 3D counter maps could facilitate the design of more potent human TRPA1 receptor agonists.     

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.     

CBP/P300溴结构域联芳基类抑制剂三维定量构效关系研究

CREB结合蛋白(CBP)和与其高度同源的P300蛋白是组蛋白乙酰化酶的两个亚型,两者通过它们的溴结构域(bromodomain,BRD)与染色质结合,目前,CBP/P300已经成为人类在肿瘤靶点领域中的研究热点。本研究基于CBP/P300溴结构域联芳基类抑制剂建立三维定量构效关系,采用比较分子力场分析法(Co MFA)和比较分子相似性指数分析法(Co MSIA)分别建立35个已知活性抑制剂的3D-QSAR模型,以确定CBP/P300溴结构域联芳基类抑制剂分子结构与生物活性之间的定量关系。Co MFA和Co MSIA模型活性数据p IC50的预测值与实验值基本一致,说明这两个模型具有较高的预测能力和统计学意义。根据Co MFA和Co MSIA模型所提供的立体场、静电场、疏水场、氢键给体场、氢键供体场等信息提出了改善此类抑制剂活性的药物设计思路,为指导设计具有更高活性的新分子和预测更加有效的CBP/P300溴结构域抑制剂提供理论依据。     

Molecular Modeling Studies of Vascular Endothelial Growth Factor Receptor Tyrosine Kinase Inhibitors Combining Molecular Docking and 3D-QSAR Methods

The vascular endothelial growth factor (VEGF) and its receptor tyrosine kinases VEGFR-2 or kinase insertdomain receptor (KDR) have emerged as attractive targets for the design of novel anticancer agents. In the present work, molecular docking method combined with three dimensional quantitative structure-activity relationships (comparative molecular field analysis (CoMFA) and comparative molecular similarity indice analysis (CoMSIA)) to analyze the possible interactions between KDR and those derivatives which acted as selective inhibitors. The CoMFA and CoMSIA models gave a cross-validated coefficient Q2 of 0.713 and 0.549, non-cross-validated R2 values of 0.974 and 0.878, and predicted R2 values of 0.966 and 0.823, respectively. The 3D contour maps generated by the CoMFA and CoMSIA models were used to identify the key structural requirements responsible for the biological activity. The information obtained from 3D-QSAR and docking studies were very helpful to design novel selective inhibitors of KDR with desired activity and good chemical property.     

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.     

Interpretation of scoring functions using 3D molecular fields. Mapping the diacyl-hydrazine-binding pocket of an insect ecdysone receptor

A method is presented for the interpretation of receptor docking score values (rough measures of binding affinities) of ligands in terms of 3D molecular field interaction contributions. The FlexX and FlexX-Pharm methods were used to dock the structures of designed sets of ligands into the ligand-binding pocket of a selected receptor. In the next step the relationship was investigated between the FlexX and CScore scores and 3D molecular fields obtained for the docked conformations of the ligands, using the CoMFA (Comparative Molecular Field Analysis) and CoMSIA (Comparative Molecular Similarity Indices Analysis) methods. The approach yielded highly significant CoMFA and CoMSIA models demonstrating that a high portion of the variance in the docking score values of the ligands can be explained by steric, electrostatic, hydrophobic, and hydrogen bond donor and acceptor molecular field interaction contributions. The approach was exemplified by using the crystal structure of the ligand-binding domain of the ecdysone receptor (EcR) of the moth Heliotis virescens as well as virtual molecule libraries of analogues of known diacyl-hydrazine (DAH) type ecdysteroid agonists. By docking appropriately designed virtual compound libraries into the DAH binding pocket of EcR followed by CoMFA and CoMSIA of the docked conformations, hitherto unexplored regions of the receptor cavity could be mapped. By mapping the significant molecular field interaction contributions onto the model of the receptor-ligand complex, important receptor-ligand interactions could be highlighted that may help the design of novel highly scored receptor ligands. An advantage of the method is that no experimental biological activity data are required to exhaustively map the receptor-binding site.     

3D-QSAR Study on the Inhibitory Activity of Flavonoids on PIM-1 Kinase

20 Typical flavonoids were selected for study on the interaction between them and PIM-1 kinase with the comparative molecular field analysis method(CoMFA) as well as the comparative molecular similarity index analysis method(CoMSIA) based on molecule docking.3D-QSAR models between these flavonoids and receptor PIM-1 kinase were established.The obtained optimal cross-validation correlation coefficient Q2 for CoMFA model was 0.582,and the non-cross-validation correlation coefficient R2 was 0.955;the corresponding values for CoMSIA model were 0.790 and 0.974,respectively.These two models showed fairly fine stability and predictive ability.In addition,molecule docking results revealed the key residues in the receptor cavity and their specific action ways with flavonoids.     

两类促肾上腺皮质释放激素(CRH)拮抗剂的分子场分析

采用比较分子相似性指数分析方法(CoMSIA)及比较分子场分析方法(CoMSIA)研究了两组CRH拮抗剂结构与活性的关系。在两种方法中,都考虑了静电场、立体场以及氢键场对构效关系的影响,结果表明采用CoMSIA得到构效关系模型要明显优于采用CoMFA得到的构效关系模型,在CoMSIA计算中,当引入疏水场时,三维构效关系模型能得到明显的改善,通过这个三维构效关系模型,可以较为精确地预测化合物的活性。通过分析分子场等值面图在空间的分布,可以观察到叠合分子周围的立体、静电以及疏水特征对化合物活性的影响。     

Three-dimensional quantitative structure-activity relationship of nucleosides acting at the A3 adenosine receptor: analysis of binding and relative efficacy

The binding affinity and relative maximal efficacy of human A3 adenosine receptor (AR) agonists were each subjected to ligand-based three-dimensional quantitative structure-activity relationship analysis. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) used as training sets a series of 91 structurally diverse adenosine analogues with modifications at the N6 and C2 positions of the adenine ring and at the 3', 4', and 5' positions of the ribose moiety. The CoMFA and CoMSIA models yielded significant cross-validated q2 values of 0.53 (r2 = 0.92) and 0.59 (r2 = 0.92), respectively, and were further validated by an external test set (25 adenosine derivatives), resulting in the best predictive r2 values of 0.84 and 0.70 in each model. Both the CoMFA and the CoMSIA maps for steric or hydrophobic, electrostatic, and hydrogen-bonding interactions well reflected the nature of the putative binding site previously obtained by molecular docking. A conformationally restricted bulky group at the N6 or C2 position of the adenine ring and a hydrophilic and/or H-bonding group at the 5' position were predicted to increase A3AR binding affinity. A small hydrophobic group at N6 promotes receptor activation. A hydrophilic and hydrogen-bonding moiety at the 5' position appears to contribute to the receptor activation process, associated with the conformational change of transmembrane domains 5, 6, and 7. The 3D-CoMFA/CoMSIA model correlates well with previous receptor-docking results, current data of A3AR agonists, and the successful conversion of the A3AR agonist into antagonists by substitution (at N6) or conformational constraint (at 5'-N-methyluronamide).