论文部分内容阅读
蛋白质分子动态模拟最常用的方法是分子动力学模拟。在过去二十年中,分子动力学模拟的算法及其在生物大分子中的应用都有很大的发展。同时,随着模拟的发展和经验的积累,人们越来越认识到模拟技术在生物大分子动态模拟中的重要性。对于大系统的分子动力学模拟,常用的方法有拉伸分子动力学(steered molecular dynamics,SMD),靶向分子动力学(targeted molecular dynamics,TMD),主成分动力学(essential dynamics,ED),构象涨跃动力学(conformational flooding,CF),随机排斥分子动力学(randomexpulsion molecular dynamics,REMD),探针分子动力学(force probe molecular dynamics,FPMD),交互式动力学(Interactive molecular dynamics,IMD)等方法。通过在模拟中引入外力、简谐势和减少蛋白质构象空间自由度等技术使得模拟的时间范围比常规动力学宽,因而可以模拟配体和受体结合、蛋白质的异构效应等生物现象。例如,拉伸分子动力学作为分子动力学的一个分支,20世纪90年代以来在研究配体的结合/解离过程中起到了重要的作用;靶向分子动力学是研究蛋白质分子不同构象之间的转变过程的重要工具。本文旨在对分子动力学模拟中的一些改进技术的研究,论文由五个部分组成。第一章,概述了分子动力学理论和自由能计算方法的原理和概况,在此基础上,简要地介绍了本文的主要工作和发现。第二章,作为预备知识,介绍了分子动力学和自由能的理论和应用。回顾了分子动力学的基本原理和发展历史,分子动力学模拟中的一些常用技术和方法,重点阐述了拉伸分子动力学和靶向分子动力学两种模拟技术。最后,从原理和应用上详细介绍了自由能计算的五种方法。第三章提出了一种方向自适应调整的拉伸分子动力学方法。首先,建立了拉伸方向调整的多目标优化模型,然后,采用基于信息熵的多种群遗传算法优化拉伸方向。这种对拉伸分子动力学的技术改进不仅显著提高了拉伸分子动力学的模拟效率,而且使模拟的成功率大为增强。新方法已经应用在细胞色素P450 3A4-甲吡酮的复合物的解离上,获得了更合理的解离通道,与传统的拉伸分子动力学方法比较,新方法具有模拟时间少,需要的解离力小,解离路径的能垒低等优点。对于传统方法不能解离的拉伸方向,新方法也可以通过对拉伸方向的优化和调整,找到一条成功的解离路径。第四章提出一种基于路径优化的非平衡自由能预测方法。首先,基于非平衡自由能预测原理建立一个分子解离的多目标优化模型,使用基于信息熵的多种群自适应遗传算法优化解离路径。然后,沿着优化的解离路径对拉伸分子动力学轨迹计算平均力势。采用基于Jarzynski方程的外推方法得到蛋白质复合物的自由能。结果表明,本方法不仅有较高的预测精度和效率,而且可以模拟解离全过程,这也为药物设计提供了重要的结构信息。第五章用靶向分子动力学研究了钙调蛋白的构象变化路径。钙调蛋白(CaM)是生物体内功能广泛的一种重要信号蛋白。根据已发现的钙调蛋白(CaM)打开和闭合两种不同的构象,研究了拮抗剂对构象变化的影响以及CaM全局构象变化路径。首先,进行了含拮抗剂和不含拮抗剂的两种常规分子动力学模拟,结果表明:CaM在独立存在时具有从闭合状态开启的趋势,它的构象动态变化推动了CaM变构功能的实现;拮抗剂具有将CaM的构象变化“锁住”在闭合状态的功能,有利于CaM控制一些激酶和磷酸酶的活性。在此基础上,进一步用靶向分子动力学模拟了CaM从闭合到打开的构象变化过程,得到一条稳定的变化路径和4个可能的过渡态构象。论文最后部分给出了本文的工作总结和对后续研究的展望。