黄酮类多药耐药相关蛋白抑制活性的构效关系研究.docVIP

黄酮类多药耐药相关蛋白抑制活性的构效关系研究 　　[摘要] 该文应用比较分子相似性指数法（comparative molecular similarity index analysis， CoMSIA）对29个黄酮类成分的结构与其多药耐药相关蛋白（multidrug-resistance associated protein，MRP）MRP1和MRP2的抑制活性进行了定量构效关系（quantitative strucuture-activity relationship， QSAR）研究。通过考查不同分子力场组合对CoMSIA建模的影响，获得了对MRP1和MRP2抑制活性起重要作用的分子力场，构建了最优QSAR模型，并探讨了黄酮类多药耐药相关蛋白抑制剂的结构修饰方法。研究结果不仅可为新药研发提供指导，也有助于部分探讨MRP1和MRP2的受体结构和含黄酮类中药的协同作用机制。 　　[关键词] 多药耐药相关蛋白；比较分子相似性指数法；黄酮类化合物；定量构效关系 　　[收稿日期] 2013-07-30 　　[基金项目] 国家自然科学基金项目 　　[通信作者] *张燕玲，副研究员，Tel：（010E-mail：collean_zhang@163.com 　　[作者简介] 乔连生，E-mail：2284607267@ 　　多药耐药（multidrug resistance，MDR）是导致肿瘤治疗失败的主要原因之一，引起MDR的原因很多，其中多药耐药相关蛋白（multidrug-resistance associated protein，MRP）介导的药物外排是MDR的重要原因之一。MRP主要包括9个成员，分别为MRP1～9。其中MRP1和MRP2能转运与葡萄糖醛酸、谷胱甘肽、硫酸盐以及其他有机阴离子相结合形成的共轭化合物而引起MDR[1]。黄酮类化合物存在于多种植物中，研究表明其对MRP1和MRP2具有抑制作用，可能开发成为MDR调节剂，但目前对于黄酮类化合物抑制剂的作用机制及构效关系尚不完全明确[2-3]。 　　文献[4]测得了黄酮类化合物对MRP1和MRP2的抑制指数，利用2D-QSAR的方法建立了黄酮类化合物活性与结构之间的关系，但其并未充分考虑化合物的三维结构特征以及与受体结合的影响。3D-QSAR不仅有助于获得药效的基团模型，还有助于推演“假想”的受体模型，指导合理药物设计[5]。比较分子相似性指数分析方法（comparative molecular similarity index analysis， CoMSIA）是较为成功且广为应用的3D-QSAR分析方法，其可以分析药物分子的立体场、静电场、疏水场、氢键给体场和氢键受体场与生物活性之间的关系[6]。因此，本研究采用CoMSIA方法分析黄酮类化合物的结构特点，探讨MRP1与MRP2受体的结构特征及黄酮类化合物抑制其活性的作用机制，为黄酮类多药耐药调节剂的研发提供理论指导。 　　1 材料与方法 　　1.1 化合物活性数据及构象优化 黄酮类化合物的活性数据来源于文献[4]，其结构母核见图1，其二维结构及抑制活性数值见表1，其中实测值为25 μmol·L-1下黄酮类化合物的抑制百分数。采用SYBYL软件构建分子的三维结构，使用Tripos标准力场进行初步能量优化，负载Gasteiger-Huckel电荷，最大迭代次数设置为1 000次，并利用系统搜寻法获得分子的最低能量构象，其他参数均采用系统的默认值。 　　图1 化合物的基本结构 　　Fig.1 The basic structure of the compounds 　　1.2 分子叠合 将文献[4]中的29个黄酮类化合物随机挑选20个化合物作为训练集，其余9个化合物为测试集（表1中标记为星号的化合物），按照其公共骨架进行叠合，模板分子选择结构最简单的黄酮分子17，分子叠合结果见图2，3。 　　1.3 CoMSIA分析 用偏最小二乘法（partial least square， PLS）对黄酮类化合物的抑制指数和5种力 　　表1 MRP1和MRP2抑制剂的活性值和预测结果 　　Table 1 The activity values and prediction of MRP1 and MRP2 inhibitor 　　续表1 　　续表1 　　注：*化合物为测试集化合物。 　　图2 公共骨架 　　Fig.2 Public skeleton 　　图3 分子叠合结果 　　Fig.3 Molecular alignment results 　　场指标进行回归。为了发现最佳的力场组合方式，推测受体的结构及其与配体的结合方式，本文同时采用不同力场组合进