论文摘要
肿瘤目前已是严重威胁人类的生命健康的常见病、多发病。全球每年新发病人数已超过1000万,而死亡人数也上升到了700万以上。而且近年来由于环境污染以及不良的生活习惯,使得肿瘤尤其是恶性肿瘤的发病率不仅逐年上升,同时还逐渐呈现年轻化的趋势。因此,对肿瘤的防治工作已经非常迫切。抗肿瘤药物是目前恶性肿瘤的主要治疗手段。但由于目前临床应用的抗肿瘤药物存在诸多问题（水溶性差、毒副作用严重以及易产生获得性耐药性）,不能满足临床治疗的需要。因此,新型的抗肿瘤药物的研究开发成为药物研究的重点领域之一,其中微管蛋白抑制剂又是目前抗肿瘤药物研究的热点方向。微管(microtubulel是真核细胞的重要组分,也是重要的抗肿瘤药物作用靶标。现有研究表明,微管蛋白（tubulin）中存在三个主要的药物结合位点:紫杉醇结合位点（Taxol site）、长春碱结合位点（Vinblastine site）以及秋水仙碱结合位点（Colchicinesite）。但在这三个位点中,秋水仙碱位点由于自身结合腔的体积较小,有利于小分子抗肿瘤抑制剂的研究。本课题以微管蛋白中的秋水仙碱结合位点为靶标,分析研究结合位点的结构性质,并模拟其与抑制剂的结合,确定各亚区域的性质、关键作用残基以及潜在的结合部位;从抑制剂结构中的共有基团入手,逐步扩充共有基团的性质,从而假设了抑制剂的结构模板。再从活性构象出发并结合已有的构效关系研究,提出抑制剂的三维药效团模型以及影响活性的部分结构因素。在此基础上,开展新型微管抑制剂的设计、合成、结构改造以及体外活性研究。一.微管蛋白中秋水仙碱药物结合位点的研究1.结合位点的性质采用多拷贝同时搜寻（MCSS）方法,以不同性质的官能团为探针对结合位点进行分析。根据官能团的空间分布位置,将结合位点分为三个疏水口袋（Ⅰ,Ⅱ和Ⅲ）以及三个极性区域（Ⅳ,Ⅴ和Ⅵ）。疏水口袋Ⅰ和Ⅱ为主要的疏水基团分布区,但是在Ⅱ中芳香性基团的分布在能量上占优且优势构象为垂直的插入其中。Ⅲ位于结合位点的底部,受其空间所限,只能接受体积较小的脂肪族基团且能量较高。极性基团主要分布在极性区域Ⅳ和Ⅴ中,前者中氢键受体、供体基团无明显差别,但后者中氢键受体基团分布的能量与数量更有优势。Ⅵ则位于Ⅰ和Ⅲ之间,极性基团有分布但能量较高。此外,β-Cys241的-SH因为带有微弱的极性,因而附近可以分布两种不同的基团,如芳香性和氢键受体基团。属于极性区域Ⅳ中α-Ser178的骨架O、α-Thr181的骨架-NH,以及位于Ⅴ中的β-Ala250～Leu255间的骨架-NH的周围均富集了不少极性基团的分布,提示了这些残基的重要性。对结合位点的静电势研究表明,整个结合位点主要表现为电中性,但是在Ⅰ的底部与Ⅵ具有负电特性,而在极性区域Ⅳ和Ⅴ却为正电特性。这一分析,与MCSS的研究结果相符,且也同晶体结构中Colchicine的构象较为匹配。而呈负电的和正电的区域,却没有被已有抑制剂占据,因而可以作为潜在的部位参与结合。2.蛋白与抑制剂的结合模式研究采用AutoDock软件,模拟了蛋白与抑制剂的结合。按照结合能量以及出现的次数,确定抑制剂的结合构象。结果表明,抑制剂主要通过与疏水口袋Ⅰ和Ⅱ间的疏水作用以及与极性区域Ⅳ间的氢键作用而与蛋白结合。抑制剂结构中体积较大且非平面的疏水部分（如三甲氧基苯环）位于疏水口袋Ⅰ中,其中的氢键受体原子（如O）或者苯环与β-Cys241的-SH存在极性作用（氢键或者阳离子-π作用）。平面的疏水部分则垂直的姿态插入Ⅱ中,其中的疏水侧链伸向Ⅱ的内部,而氢键供体原子或者受体原子,分别与α-Ser178或者α-Thr181形成氢键。此外,抑制剂中位于Ⅰ和Ⅱ二者之间的结构部分上的O原子（如=O）的空间位置比较相似,均指向loop区的β-Ala250的骨架-NH。对AutoDock对接得到的抑制剂结合构象又进行了Affinity的柔性对接方法考察结合位点中残基的构象改变情况。结果表明:主要的结合残基β-Cys241、β-Leu248、α-Thr179、α-Val181以及β-Asn352侧链的构象改变较大,可能是因为随着周围配体中的基团位置的改变而产生相应的变化。但是结合位点中的其余残基的构象均未发生明显改变。3.确定关键作用残基以及潜在药物结合部位根据MCSS中官能团的分布位置与能量,提示了极性区域Ⅰ中的β-Cys241、Ⅳ中的α-Scr178和α-Thr181,以及Ⅴ中的β-Ala250～Leu255等残基可能具有一定的重要性。在抑制剂的结合构象中,β-Cys241、Ⅳ中α-Ser178和α-Thr181均以氢键的形式参与了结合。同时,参考已有的构效关系研究,失去与β-Cys241或者α-Ser178以及α-Thr181间的极性作用,对活性具有一定的负面影响,因而明确的指出它们为关键的结合残基。由于Ⅴ中的β-Ala250～Leu255之间的骨架-NH位置集中且又有开阔的空间,而且MCSS计算中氢键受体官能团的能量较低,因而是一个比较理想的氢键供体区域。然而,对接研究表明,现有的抑制剂均未与该区域结合。此外,现有抑制剂一个共同问题是较低的脂水分配系数,因此,该区域可作为潜在的药物部位,通过引入新的极性基团改善化合物的水溶性及增加氢键作用而增强与蛋白的结合。疏水口袋Ⅲ和极性区域Ⅵ也并没有参与结合,但是因为它们直接位于重要的疏水口袋Ⅰ的下方,因而也可以作为新的结合部位为改造Ⅰ中的基团提供指导。二.秋水仙碱位点抑制剂的结构特征研究1.二维结构特征与结构模板从经典的抑制剂结构中的共有基团受到启发,假设了其结构特征。随后按照性质特点,对共有基团进行扩展,进而此特征逐步扩大,即:由1～4个原子作为桥连连接起来的2个疏水环或基团（分别以A和B表示）,其一与TMP类似,另一与邻甲氧基酚环或者1,2-亚甲二氧基苯环类似。最终提出了二维的结构模板,并结合已有的结构,更进一步指出了具体的结构类型。2.抑制剂的三维药效团模型模型的构建按照对接模拟得到的抑制剂结合构象中官能团的空间位置,并参考关键的结合残基的性质,提出了由以下六个元素组成的药效团模型,即:2个疏水中心（H1,H2）、疏水侧链H3、2个氢键受体（A1,A2）和极性原子（P）。H1和A1位于疏水口袋Ⅰ中,A1与β-Cys241存在氢键作用;H2和H3位于疏水口袋Ⅱ中;P位于极性区域Ⅳ中,与α-Ser178或α-Thr181形成氢键;A2靠近loop区的β-Ala250的骨架-NH。同时由于极性区域Ⅴ可以作为潜在的结合部位,因而在A2和loop之间,还假设存在另一个氢键受体点A3以同loop区中的骨架-NH形成氢键。其中H1和A1组成了平面P1,H2、H3以及P构成平面P2,但是平面P1-P2之间存在42～60的夹角。而结合已提出了结构特点,药效团元素中的H1和A1组成了模板的A部分位于疏水口袋Ⅰ中,H2、H3以及P构成B部分位于疏水口袋Ⅱ和极性区域Ⅳ中,A2或A3则属于桥连部分以连接A和B部分。3.影响活性的其它结构因素在模拟了不同结构改造的现有抑制剂与蛋白的结合时发现,改变H1部分疏水基团的体积,改变H2部分的平面性,以及改变结构模板中的AB部分位于桥连的空间位置时,对活性具有重大影响。因而参考已有的结构改造研究,指出了影响活性其它的结构因素为:除了必须的疏水性外,H1部分的体积大小、H2部分的平面性以及桥连部分是否能保证将AB位于其同侧。三.新型小分子抑制剂的计算机辅助设计以上述研究为基础开展先导化合物的设计工作。首先根据三维药效团模型中各个药效元素的性质与特点选择不同的官能团;随后按照结合位点各个区域的性质,采用LUDI挑选不同的连接片段以手工方法将上述官能团连接起来,得到了多种不同的先导模板。最后通过接模拟它们与蛋白的结合模式,并综合考虑其与蛋白的结合能量、分子间相互作用能、LUDI打分结果、形成的氢键个数以及酯水分配系数,最终确定了咪唑酮环为桥连的结构作为先导化合物。在设计的时候,还考虑了以下几个因素:遵循提出的药效团模型,保证结构中存在六个药效点;符合提出了结构特征与结构模板,结构中仍然存在A、B以及桥连这三个部分,且连接片段必须保证A、B部分位于桥连的同侧;匹配结合位点的性质,按照疏水口袋Ⅰ和Ⅱ,以及Ⅳ的特征选择合适的官能团;利用新发现的潜在结合区域,引入新的亲水基团;合成的难度以及结构改造的空间。四.合成与结构修饰研究参考相关文献,提出了一条反应条件温和,后处理简便且易于结构改造的合成路线。即:以三甲氧基苯甲酸为原料,经酰氯化后与乙酰乙酸乙酯发生酰基化反应,随后脱除乙酰基后再与NBS反应得到的溴代产物,在碱性条件下和脲直接环合得到关键中间体（Ⅵ）。此后,中间体（Ⅵ）按照桥连部分改造的需要,又分别与不同的芳胺、酚等反应;或者还原为醇或醛后,得到反酯、烯胺或者乙烯结构。结构改造中,还对B部分进行结构修饰,考察了B部分中取代基的体积大小、供电子/吸电子、取代基的位置以及个数等因素对活性的影响。最终得到了20个化合物,经数据库检索所有结构均为首次报道。同时对中间体咪唑酮的合成方法进行了改进,提出了先溴代再与脲环合的2步反应路线,以乙酸乙酯:石油醚（v:v=1:3）重结晶即可得到产物,总收率达到55.4%。与已有文献中的先亚硝化、经还原后再与异氰酸反应的合成方法相比,新路线的反应步骤简短、处理方便,而且提高了收率。五.抑制剂的活性测试研究肿瘤细胞生长抑制实验实验表明,所有化合物对六种不同的肿瘤细胞都具有一定的生长抑制作用,体现了广谱的抗肿瘤活性。其中对LOVO（人结肠癌细胞）、CEM（人白血病细胞）、MDA-MB-435（人乳腺癌细胞）抑制作用明显。化合物CON-2、CON-3、CON-8、CON-10、CON-11以及CON-13对这三个细胞的活性低于5μM,其中CON-7、CON-14最高活性分别达到了176nM和330nM。此外,CON-7、CON-14对除SKOV-3（人卵巢癌细胞）外的其它五种瘤细胞具有明显的生长抑制作用,活性范围主要集中在0.17nM～5μM之间。抗血管实验表明,大部分的化合物都能较好的抑制人脐静脉内皮细胞（ECV-304）的生长,9个化合物的活性低于16μM。同时,此实验结果与抑瘤实验活性存在一定的相关性。初步的构效关系分析表明:B部分的芳环中有供电子基团或者氢键供体基团的存在对活性有增强作用,而间位取代的活性又强于对位和邻位取代。生物活性数据不仅验证了最初的设计思想,更为进一步的研究设计新结构类型的微管蛋白抑制剂提供了重要参考。