Study on structure-activity relationship of mutation-dependent herbicide resistance acetohydroxyacid synthase through 3D-QSAR and mutation

Seventy-four sulfonylureas were synthesized and tested for their inhibitory activity against the whole enzyme of E. coli acetohydroxyacid synthase (AHAS, EC 2.2.1.6) isoenzyme II, and 3D-QSAR analyses were performed based on these inhibitory activities. The binding conformation of chlorimuron-ethyl, a commercial herbicide of AHAS, in the crystal structure of AHAS complex was extracted and used as template to build the initial three-dimensional structure of other sulfonylureas, and then all structures were fully geometry optimized. After systematic optimization of the alignment rule, molecular orienta- tion, grid space and attenuation factor, two satisfactory models with excellent performances (CoMFA: q2 = 0.735, r2 = 0.954, n = 7, r 2 pred = 0.832; CoMSIA: q2 = 0.721, r2 = 0.913, n = 8, r 2pred = 0.844) were estab- lished. By mapping the 3D contour maps of CoMFA and CoMSIA models into the possible inhibitory active site in the crystal structure of catalytic subunit of yeast AHAS, a plausible binding model for AHAS, with best fit QSAR in the literature so far, was proposed. Moreover, the results of 3D-QSAR were further utilized to interpret resistance of site-directed mutants. A relative activity index (RAI) for AHAS enzyme mutant was defined for the first time to relate the 3D-QSAR and resistance of mutants. This study, for the first time, demonstrated that combination of 3D-QSAR and enzyme mutation can be used to decipher the molecular basis of ligand-receptor interaction mechanism. This study refined our understanding of the ligand-receptor interaction and resistance mechanism in AHAS-sulfonylurea system, and provided basis for designing new potent herbicides to combat the herbicide resistance.     

HLA-A*0201限制性CTL表位肽的分子力场分析

运用比较分子力场(CoMFA)和比较相似性指数分析(CoMSIA)方法研究了50个HLA-A* 0201限制性CTL表位九肽结构与亲和性间的关系,另外15个表位九肽作为预测集用于检验模型的预测能力。结果表明,采用CoMSIA得到的构效关系模型(q2 =0.608,r2 =0.987,F=440.4,SD=0. 111)要明显优于采用CoMFA得到的构效关系模型。在CoMSIA计算中,当引入疏水场时,三维构效关系模型得到明显改善,通过该三维构效关系模型,可较精确地估算预测集中15个CTL表位肽与HLA-A* 0201间的亲和力(r2pred=0. 703,SD=0. 368)。通过分析分子场等势面图在空间的分布,可以观察到表位肽分子周围的立体及疏水特征对表位肽与HLA-A* 0201间结合亲和力的影响,从而为进一步对CTL表位肽进行结构改造并基于此进行治疗性疫苗分子设计提供理论基础。     

2-苯基4-吡喃酮类COX-2抑制剂的三维构效关系研究

 环氧化酶（cyclooxygenase,COX）家族是非甾体类抗炎药物的作用靶酶,其中COX-2是 NSAIDs药物发挥作用与治疗炎症有关的靶点。使用三维定量构效关系研究方法研究2-苯基4-吡喃酮类COX-2抑制剂的分子结构与生物活性之间的定量关系,并构建三维定量构效关系模型。CoMSIA模型证实该类化合物的结构参数和其抑制活性存在明显的相关性（R²为0.929）,且有较好预测能力（q²为0.746）;同时使用“留五法”证实模型的稳定和可靠。三维构效关系模型系数等势图显示Subset A和B的R取代基4位上官能团变化对活性影响最大,在结构修饰时需要优先考虑;Subset A的R取代基2位上官能团的影响次之。     

黄酮类化合物与雌激素受体作用的三维定量构效关系

黄酮类化合物是一种普遍存在于食物和植物中的内分泌干扰物. 选取20个典型的黄酮类化合物, 采用基于分子对接(docking)的比较分子场分析法(CoMFA)和比较分子相似性指数(CoMSIA)方法, 研究它们与雌激素受体(ER?)的相互作用, 建立了这些黄酮类化合物和雌激素受体之间的三维定量构效关系(QSAR)模型, 得到CoMFA模型的最佳交互验证相关系数Q²_LOO为0.845, 非交互验证相关系数R²为0.988, 相应的CoMSIA模型的交互验证相关系数Q²_LOO为0.670, 非交互验证相关系数R²为0.990, 说明模型具有很强的稳定性和预测能力. 此外, 分子对接结果揭示了受体腔中的关键残基及其与黄酮类化合物作用的具体方式.     

一类细胞凋亡诱导剂的结构特征研究

许多抗癌药物通过引起癌细胞的凋亡来达到治疗目的,因此开发能诱导癌细胞凋亡的药物一直是癌症研究的热点.以153个结构各异的具有caspase-3激活作用的细胞凋亡诱导剂为研究对象,通过三维定量构效关系(3D-QSAR)分析,得到了一个基于位阻场、静电场、疏水场、氢键供体场和受体场参数的最优比较分子相似性指数法(CoMSIA)模型.该模型(Q2=0.51、R2ncv=0.89和R2pred=0.82)具有良好的内部一致性和预测能力,通过对模型的等值线图分析发现:(1)R2取代基提高诱导剂活性的因素包括位阻大、负电性大、有亲水性和具有氢键供体作用;(2)位阻适中和/或具有氢键供体作用的R4取代基,以及位阻小和/或亲水的R1取代基对提高分子的活性是有利的;(3)N-甲基-N-苯基-1-萘胺类细胞凋亡诱导剂具有正电性的A环3号位或C环7号位取代基,以及具有负电性的B环1号位取代基有利于提高诱导剂的生物活性.对该模型的分析更有利于认识细胞凋亡诱导剂的分子特征,并为设计和开发新型细胞凋亡诱导剂提供理论指导.     

部分苯酚衍生物生物降解性的2D和3D-QSBR研究

应用B3LYP/6-311G**水平计算得到的量子化学参数,在实验数据的基础上,建立了部分苯酚衍生物生物降解性能(二氧化碳生成量(PCD))的定量结构-生物降解性相关模型(2D-QSBR),该模型包括分子平均极化率(α)和熵(SΘ)两个参数,其中α对PCD的影响更为显著,模型的相关系数R2=0.933,交叉验证相关系数q2=0.894.该结果优于midix,6-31G*计算水平下得到的2D模型,此外,应用基于分子模拟技术的比较分子相似性指数(CoMSIA)方法,建立了3D-QSBR模型,模型的R2=0.964,q2=0.716,化合物的氢键供受体特性对降解能力的影响较为显著,其次是化合物的疏水性.结果表明,2D和3D-QSBR模型都具有良好的稳定性和预测能力,可以预测同类化合物的PCD值.     

5-HT 6受体拮抗剂三维定量构效关系研究

为指导合成高效的5-HT6受体拮抗剂,采用比较分子场分析(CoMFA)和比较相似性指数分析(CoMSIA)对143个5-HT6受体拮抗剂数据进行了三维定量构效关系(3D-QSAR)研究,分别得到了具有良好可靠性和预测能力的CoMFA(Q2=0.513,R2ncv=0.864,R2pre=0.731)和CoMSIA模型(Q2=0.515,R2ncv=0.844,R2pre=0.777).由模型的等势线图分析,可得如下结论:大体积及/或电负性较大的R2取代基、疏水性R3和疏水性苯环R1位取代基、可作氢键受体的R2取代基、可作氢键供体的R3取代基有助于增大活性.这些结论能够更好地帮助理解5-HT6受体拮抗剂的抑制机理,并为今后的药物设计与合成提供新思路.     

6-(1-萘甲基)取代S-DABO类逆转录酶抑制剂的3D-QSAR研究

以活性化合物15为例,采用分子对接的方法模拟了化合物与HIV-1逆转录酶的作用,从而探讨了系列6-萘甲基取代S-DABO类化合物的作用机制.并基于分子对接后的活性构象,应用比较分子力场分析(CoMFA)和比较分子相似因子分析(CoMSIA)法对该类化合物的三维定量构效关系进行了研究,建立了有较好预测能力3D-QSAR模型,为该类化合物进一步的结构优化提供理论指导.     

地棘蛙素类似物的CoMFA和CoMSIA的3D-QSAR研究

采用比较分子力场 (CoMFA)分析法和比较分子相似性指标场(CoMSIA)分析法,进行偏最小二乘法(PLS)分析,建立了地棘蛙素类化合物的三维构效模型.通过验证,说明该系列化合物分子立体场、静电场、氢键场和疏水场的分布与生物活性之间有良好的相关性.用模型预测同系列测试集分子,结果与实验值偏差较小.预测结果表明,该力场模型有一定的预测能力.并得出了新的药效团定义模型,可用来指导设计新的地棘蛙素类化合物烟碱型乙酰胆碱配体.     

一类鬼臼毒素衍生物抗癌(KB细胞)活性的三维定量构效关系

采用比较分子场分析(CoMFA)和比较分子相似因子分析(CoMSIA)方法,对99个鬼臼毒素衍生物的抗癌活性系统地作了三维定量构效关系研究,均得到了较为理想的模型.coMFA所得模型的交叉验证系数q2=0.644,非交叉验证相关系数r2=0.930,标准偏差S=0.134,F=161.325,立体场与静电场贡献分别为64.6%和35.4%.CoMSIA所得模型的q2=0.764,r2=0.964,S=0.096,F=324.112,立体场、静电场、疏水场、氢键给体场、氖键受体场贡献分别为9.5%,22.0%,21.9%,23.9%和22.7%.两模型的等值面图为改造此类化合物的结构、提高药物活性提供了理论依据和研究方向.