世界上几乎近三分之一的蛋白质都含有金属离子,金属离子对于蛋白质的结构和功能甚至大多数生命过程都是极其关键的。金属蛋白在细胞中发挥着重要作用,如酶催化,蛋白的传输和储存,信号传递等。例如,钙结合蛋白与钙离子结合是细胞信号传递的重要步骤。铜蛋白与体内的氧化还原反应密切相关。含铁蛋白如肌红蛋白和血红蛋白具有携氧和储氧的能力。锌指蛋白与多种生物功能相关,如DNA识别、RNA聚合、转录,激酶活性调控等等。目前对于蛋白质等生物大分子的研究主要局限在经典力学的范畴。常用的分子力场,比如AMBER, CHARMM, GROMOS, OPLS等,在蛋白、核酸、糖类以及脂质等方面的研究取得了巨大的成功。然而目前对于金属蛋白的理论研究尚处于初级阶段。现在较成熟的处理金属离子的方法是键连模型,即金属离子显式地同与之配位的原子相连。但是在具有催化功能的金属蛋白的活性区域,金属离子价态会随所处状态变化,这就需要应用非键连模型。但是目前流行的分子力场中的非键连模型处理配位数可变体系均不够精确。因此,开展对金属离子以及金属蛋白的定量化的理论计算研究对于我们深入了解相关生物物理性质,预测和设计新的蛋白及功能等有着非常重要的意义。为精确描述蛋白质内金属离子的行为,本文基于最近开发的一种新型电荷,蛋白质专一的极化电荷(PPC),着力发展了可应用于金属蛋白的电荷方案。我们利用PPC的思想,用全量子力学的方法计算金属离子的配位区域,利用包含极化效应的电子密度拟合出金属以及其第一配位层的电荷。本文中,我们首先详细研究了拟使用的量化计算方法的有效性,和分子动力学模拟中的抽样方法,最后基于量化计算和动力学模拟中获得的经验,将金属体系的PPC电荷应用于双锌蛋白DFsc以及DF1这两个体系。受以上思路启发,第三章中,我们考察了一个具有催化氧化活性的人工设计的蛋白质DFsc。DFsc是一个人工设计的小型蛋白,包含两个金属离子。这两个金属离子可以是锌、铁、锰、铬等,但是只有铁和锌配位的DFsc蛋白才具有催化氧化活性。在该蛋白中,两个相邻的锌离子位于溶剂可及通道底部,不仅具有催化活性,而且配位数可变。这种情况下只能应用非键模型来对其模拟。在该章中,我们应用了可极化电荷PPC模拟DFsc,并使电荷随时间更新以模拟电荷的变化。然后我们用四种AMBER力场(AMBER94, AMBER99, AMBER99SB, AMBER03)进行模拟作为对比。MD模拟的结果表明,首先对RMSD的分析发现,DFsc在我们设计的模拟过程中结构稳定。和核磁共振(NMR)耦合数据的比对发现,我们模拟得到的结构同NMR信号具有很高的一致性。这说明我们的方案对此蛋白质的骨架的模拟是成功的。其次,对金属配位集团的结构分析(配位键的键长平均值,以及它们随时间的变化)表明,两个锌离子在模拟过程中的配位数发生了变化,平均之后发现一个为四配位,另一个为五配位。这说明两个锌离子的不对称性得到了保持。而且,我们的模拟正确反映出,两个锌离子即使在蛋白质骨架达到平衡之后,仍具有很大的涨落。根据前人的结果,我们推测该涨落是由于锌离子具有反应活性引起的。但是在对照的AMBER力场的结果中,锌离子极其稳定。最后,我们通过对锌离子间距和附近水分子运动的关系的交叉相关分析,发现水分子对锌离子的扰动是造成该体系中锌离子位置涨落的原因。但是,DFsc中的锌离子的位置的涨落是否是APPC电荷本身带来的效果,而不是体系内在的特征,尚不完全明确。于是在第六章中,我们选择DFsc的一个类似物,DF1,进行了进一步的模拟。DFl与DFsc的不同之处在于,在DFsc中存在的溶剂可及通道在DFl中被变异的残基阻隔,使得锌离子被包裹在蛋白质内部,不能跟水接触,这样就没有催化氧化活性。我们经过模拟发现,我们模拟得到的锌离子的涨落,在DF1中并不存在。DF1中的两个锌离子之间的距离稳定在3.77A左右(实验值为3.92A)。这进一步证明了我们的电荷方案可以正确反映金属离子在蛋白质中的动力学行为。而AMBER力场高估了锌离子在蛋白质的催化部位的稳定性。目前密度泛函方法是最流行的量子化学计算方法。但是对于含金属体系的性质的计算,比如光谱,并没有完全普适的泛函与基组的组合方案。要决定采用何种泛函与基组的组合,需要预先在小体系中进行测试。第五章中,我们详细地研究了目前流行的密度泛函对含金属的体系二茂铁(ferrocene)的几何结构和激发光谱的影响。计算结果表明,流行的泛函和基组的组合几乎均可给出正确的几何结构。但是不同的组合对紫外可见光谱的计算值影响则非常大。在该章中我们发现PBEO/6-31++G**配合可极化导体模型的组合可以给出非常精确地紫外-可见吸收谱,并成功应用于二茂铁的类似化合物二苯铬,以及几个实验合成的化合物(acetylferrocene(AcFc), hydroxyethylferrocene(HyFc), vinylferrocene(ViFc), and ethynylferrocene(EtFc))的光谱模拟。抽样是分子动力学模拟的一个重要环节,正确而高效的抽样方案可以得到非常精确的自由能。这对蛋白质-配体相互作用能的计算非常重要。为了测试在已有力场下应用高精度抽样技术,是否可以得到足够精确的自由能,第六章中,我们使用了目前高精度的热力学积分(TI)方法,在传统的AMBER力场下,研究了蛋白质Avidin与带电配体biotin形成的复合物Avidin/biotin的配位作用能,和该蛋白质与另一个带电配体形成的复合物Avidin/BTN2的配位作用能的差值,△△G。我们发现,即使模拟的轨迹非常稳定,模拟时间足够长,计算得到的△△G(约为lkcal/mol)仍然比实验值(6kcal/mol)小很多,反而类似于此前另一研究组(B. Kuhn and P. A. Kullman)使用AMBER力场,用另外的方法(MM-PB/SA)计算得到的数值。在模拟过程中,体系中的四个氢键发生了变化。其中两个是配体作为氢键受体,而另外两个是配体作为氢键给体。这使得主要由氢键的变化导致的自由能变正负抵销,最终使A△G过小。这证明了传统不包含极化效果的力场处理带电体系会有根本性的偏差。