金属蛋白力场的开发与应用

《金属蛋白力场的开发与应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、Thedevelopment11··r’anC【aDDnCatlonSOtmetalloproteinforcefieldbyYongleLiDirectedbyProfessorDaiqianXieandJohnZHZhangSchoolofChemistryandChemicalEngineeringNanjingUniversityAug2011SubmittedinpartialfulfilmentoftherequirementsforthedegreeofPhDinPhysicalChemistry毕业论文题目：金属蛋白力场的开发与应用物理化学专业200

2、8级：李永乐指导教师(姓名、职称)：谢代前教授博士摘要世界上几乎近三分之一的蛋白质都含有金属离子，金属离子对于蛋白质的结构和功能甚至大多数生命过程都是极其关键的。金属蛋白在细胞中发挥着重要作用，如酶催化，蛋白的传输和储存，信号传递等。例如，钙结合蛋白与钙离子结合是细胞信号传递的重要步骤。铜蛋白与体内的氧化还原反应密切相关。含铁蛋白如肌红蛋白和血红蛋白具有携氧和储氧的能力。锌指蛋白与多种生物功能相关，如DNA识别、RNA聚合、转录，激酶活性调控等等。目前对于蛋白质等生物大分子的研究主要局限在经典力学的范畴。常用的分子力场，比如AMBER，CHARMM，GROMOS，O

3、PLS等，在蛋白、核酸、糖类以及脂质等方面的研究取得了巨大的成功。然而目前对于金属蛋白的理论研究尚处于初级阶段。现在较成熟的处理金属离子的方法是键连模型，即金属离子显式地同与之配位的原子相连。但是在具有催化功能的金属蛋白的活性区域，金属离子价态会随所处状态变化，这就需要应用非键连模型。但是目前流行的分子力场中的非键连模型处理配位数可变体系均不够精确。因此，开展对金属离子以及金属蛋白的定量化的理论计算研究对于我们深入了解相关生物物理性质，预测和设计新的蛋白及功能等有着非常重要的意义。为精确描述蛋白质内金属离子的行为，本文基于最近开发的一种新型电荷，蛋白质专一的极化电荷

4、(PPC)，着力发展了可应用于金属蛋白的电荷方案。我们利用PPC的思想，用全量子力学的方法计算金属离子的配位区域，利用包含极化效应的电子密度拟合出金属以及其第一配位层的电荷。本文中，我们首先详细研究了拟使用的量化计算方法的有效性，和分子动力学模拟中的抽样方法，最后基于量化计算和动力学模拟中获得的经验，将金属体系的PPC电荷应用于双锌蛋白DFsc以及DFl这两个体系。受以上思路启发，第三章中，我们考察了一个具有催化氧化活性的人工设计的蛋白质DFsc。DFsc是一个人工设计的小型蛋白，包含两个金属离子。这两个金属离子可以是锌、铁、锰、铬等，但是只有铁和锌配位的DFsc蛋

5、白才具有催化氧化活性。在该蛋白中，两个相邻的锌离子位于溶剂可及通道底部，不仅具有催化活性，而且配位数可变。这种情况下只能应用非键模型来对其模拟。在该章中，我们应用了可极化电荷PPC模拟DFsc，并使电荷随时间更新以模拟电荷的变化。然后我们用四种AMBER力场(AMBER94，AMBER99，AMBER99SB，AMBER03)进行模拟作为对比。MD模拟的结果表明，首先对RMSD的分析发现，DFsc在我们设计的模拟过程中结构稳定。和核磁共振(NMR)耦合数据的比对发现，我们模拟得到的结构同NMR信号具有很高的一致性。这说明我们的方案对此蛋白质的骨架的模拟是成功的。其次

6、，对金属配位集团的结构分析(配位键的键长平均值，以及它们随时间的变化)表明，两个锌离子在模拟过程中的配位数发生了变化，平均之后发现一个为四配位，另一个为五配位。这说明两个锌离子的不对称性得到了保持。而且，我们的模拟正确反映出，两个锌离子即使在蛋白质骨架达到平衡之后，仍具有很大的涨落。根据前人的结果，我们推测该涨落是由于锌离子具有反应活性引起的。但是在对照的AMBER力场的结果中，锌离子极其稳定。最后，我们通过对锌离子间距和附近水分子运动的关系的交叉相关分析，发现水分子对锌离子的扰动是造成该体系中锌离子位置涨落的原因。但是，DFsc中的锌离子的位置的涨落是否是APPC

7、电荷本身带来的效果，而不是体系内在的特征，尚不完全明确。于是在第六章中，我们选择DFsc的一个类似物，DFl，进行了进一步的模拟。DFl与DFsc的不同之处在于，在DFsc中存在的溶剂可及通道在DFl中被变异的残基阻隔，使得锌离子被包裹在蛋白质内部，不能跟水接触，这样就没有催化氧化活性。我们经过模拟发现，我们模拟得到的锌离子的涨落，在DFl中并不存在。DFl中的两个锌离子之间的距离稳定在3．77A左右(实验值为3．92／鼙)。这进一步证明了我们的电荷方案可以正确反映金属离子在蛋白质中的动力学行为。而AMBER力场高估了锌离子在蛋白质的催化部位的稳定性。目前密度泛