论文部分内容阅读
基于物理学原理的计算化学,其方法和技术可广泛应用于生物分子体系的模拟研究。在蛋白质功能和结构的研究中,可用于蛋白质三维结构的同源建模,蛋白质-蛋白质复合物的分子对接,和分子动力学模拟蛋白的柔性及功能分析。在计算机辅助药物设计中,基于结构的药物虚拟筛选已经被证明是非常有效的方法。在本论文的研究工作中,我们针对两个很有意义的生物学问题开发了相应的计算方法并开展了研究,一是开发了一种基于物理学原理的预测激酶抑制剂选择性的算法;二是设计了一套基于结构信息的从头预测蛋白质-蛋白质复合物的对接策略。如何合理设计具有选择性的激酶小分子抑制剂是激酶药物开发中一个具有挑战性但是很重要的问题。我们开发的算法主要包括以下步骤:1.选择具有代表性的处于活性构象的激酶晶体结构;2.用比较对接的方法构建配体受体复合物;3.在配体存在下搜索结合位点的氨基酸构象;4.用MM-GB/SA的打分函数评估最终的配体和受体之间的结合自由能。我们首先回顾性地评价了这种算法在预测17个常见的激酶小分子抑制剂和143个激酶选择性的表现。两个相关性分析参数,排名前20%的富集因子(EF20)平均值2.81和预测指数(PI)平均值0.35表明我们的算法能够预测一个激酶小分子抑制剂的选择性,尤其是对有选择性好的抑制剂具有更好的预测能力。然后我们将这种算法应用到化疗药物米托蒽醌(MX)可能结合的激酶靶点中。根据预测结果,我们取排名前五位的激酶进行体外活性检测。结果表明米托蒽醌对PIM1激酶具有低纳摩尔级的抑制活性。同时,我们发现米托蒽醌也能在肿瘤细胞中抑制PIM1底物的磷酸化。我们进一步解析了米托蒽醌和PIM1的复合物晶体结构,从而揭示了米托蒽醌抗癌机制的一种新的结构基础和作用机制。虽然米托蒽醌最初被认为是通过嵌入DNA和抑制Ⅱ型拓扑异构酶起作用的,我们的研究结果第一次表明米托蒽醌还可能通过抑制PIM1激酶发挥其抗癌疗效和特异性。蛋白激酶在细胞内行使功能时是被严格调控的,包括激酶自身二聚化或寡聚化,其它结构域或蛋白的相互作用。例如JAK2激酶结构域(JH1结构域)就受到相邻的假激酶结构域(JH2结构域)的调控,但是具体两个结构域之间的作用机理依然不清楚。我们首先采用MuInf,一种能够鉴定蛋白扭转角相关性位点的算法来预测JAK2假激酶结构域和激酶结构域可能的变构调节位点,然后将这些变构调节位点作为蛋白质与蛋白质相互作用的界面来构建起始模型,随后我们用蛋白质对接和分子动力学模拟来更合理地修正我们的结构模型。我们构建的JAK2复合物模型揭示了JAK2假激酶结构域自我抑制的可能的一种作用机制,它通过重要的疏水性作用和广泛的静电相互作用将激酶结构域中的活化茎环限制在非活性状态,同时也阻碍αC螺旋的转动。进一步的分子动力学模拟表明在大多数骨髓增生性肿瘤患者中发现的JAK2的一个获得性突变(V617F),可能通过减弱假激酶结构域的这种抑制作用来激动JAK2。重要的是我们用的假激酶结构域在结合界面上的氨基酸的突变实验证实,这些重要氨基酸的突变确实能够激活JAK2激酶活性,并能促进用这些突变体转化的BaF3/EpoR细胞增殖。尽管还可能存在其它的假激酶结构域介导的JAK2激酶的调控机制,我们构建的JAK2假激酶和激酶复合物模型提供了一种合理的假说,这有助于进一步设计实验来阐明假激酶结构域在正常和病理中的JAK2的调控机理。同时我们这种逐级蛋白对接的策略能够更广泛地应用到研究其它蛋白质-蛋白质相互作用的体系。