论文部分内容阅读
蛋白质是生命活动的主要承载体,研究蛋白质的结构和功能有着非常重要的工业和医药价值。存在于蛋白质表面及其内部的静电场作为蛋白质内部主要相互作用的来源影响着蛋白质结构和功能的方方面面。对于蛋白质内部静电性质的研究吸引了大批实验和理论科学家的关注。振动斯塔克光谱在实验上被广泛的应用于探测蛋白质内部动力学行为的变化。通过向蛋白质内部插入含有某些特殊基团的探针的方法,通过线性斯塔克光谱可以在实验上比较精确的测量蛋白质内部静电场变化。理论上基于分子力场的动力学模拟方法为在微观尺度上揭示蛋白质的各种性质提供了非常有力的工具。基于分子力场的计算结果可以直接与实验上测得的数据进行比较,从而可以在分子层面上更加深入的分析和解释实验数据,甚至做出合理的预测。现在常用的分子力场模型比如AMBER,CHARMM,OPLS等在静电性质的计算上大多采用平均场近似的电荷模型的方法,缺少了由于蛋白质内部特殊的静电环境所带来的静电极化效应。虽然这种近似的电荷模型方法对于蛋白质动力学行为模拟的影响是有限的,但是越来越多的理论工作指出在蛋白质的某些性质的计算中已经达到了不可忽略的程度。在前人工作的基础上,我们把运用了量子化学计算方法从而考虑了蛋白质特殊的静电环境带来的极化效应的电荷polarized protein-specific charge (PPC)运用到人体醛糖还原酶(hALR2)线性斯塔克光谱的计算中,同时进行了基于传统的非极化AMBER力场的光谱计算。通过直接和实验值进行比较,我们发现基于极化电荷模型的力场比传统非极化的AMBER力场能给出更切合实验数据的计算结果,说明极化效应在蛋白斯塔克光谱的计算方面有着非常重要的作用,同时也揭示出现在常用的非极化的AMBER力场在与静电直接相关的性质的计算上需要进一步的提高。通过对蛋白质内部电场和stark光谱的研究,我们发现经典的力场模型在借助于高等级的量子化学计算后,能够更好重复试验的结果,而在现存的理论计算方法中,基于量子力学的计算方法是最为精确和可靠的。对于小的体系,在采用了适当的方法后,我们通常可以完成化学精度甚至于光谱精度的计算。理论化学家总是希望能够将量子化学的计算方法应用到尽可能大的分子体系中去。虽然近些年来计算机的发展对于含有几百甚至上千原子的分子,已经基本可以实现量子力学级别的计算,但是现在许多科学研究中要研究的分子所包含的原子的数目通常要远超出这个数值。比如像蛋白质分子,要进行传统意义上的量子化学的计算仍然是一项比较巨大的挑战。为了解决这样的问题,我们课题组自主开发和编写了一套基于分块的线性标度的量子化学的计算方法:Molecular Fractionation with Conjugate Caps (MFCC,分子碎片共轭帽),这种方法已经被应用到蛋白质性质计算的很多方面。为了提高这种方法对于蛋白质总能量的计算精度一种新的Generalized Molecular Fractionation with Conjugate Caps/Molecular Mechanics (GMFCC/MM,推广的分子碎片共轭帽)方法被提出来。基于这种方法我们进一步发展了静电嵌入的推广的分子碎片共轭帽的方法:Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Cap(EE-GMFCC,静电嵌入的推广的分子碎片共轭帽)以便更精确地计算蛋白质总能量。在这种方法中蛋白质的总能量由各个碎片计算得到的能量的线性组合得来。在所有碎片分子的计算中,通过简单的静电嵌入的方式来引入除去碎片分子以外的环境效应。通过对18个真实蛋白质体系在HF/6-31G*水平上的计算我们发现与传统意义上的全体系的量子化学的计算结果相比较,静电嵌入的方案相对于原来的推广的分子碎片共轭帽的方法把平均误差缩减到2.39kcal/mol。我们通过采用不同的两体截断半径的方法测试了几个含有几百原子具有三维结构的蛋白质体系,结果表明两体的相互作用矫正对于计算精度也有较大的贡献。EE-GMFCC的方法在MP2和DFT水平上的计算相对于全体系的量子计算的结果误差也都只有几个kcal/mol。而相对于全体系量子化学的计算对于含有几百个原子的蛋白质我们这种方法的计算效率有非常明显的提高,从而可以把量子化学的计算方法推广到任意大小尺寸的蛋白质分子中。蛋白质内部电场影响着蛋白质的结构和功能,而外电场对蛋白质空间结构的完整性同样有着重要的影响,有实验表明微波场能够在不引起热效应的前提下改变蛋白质的三维构象。外电场对蛋白质动力学的非热效应在实验上引起了广泛的关注。但是就现在的状况而言实验上的数据还不能对这些效应作出原子层面上的解释,借助于分子动力学的模拟我们能够在原子尺度上揭示其微观机理。应用分子动力学模拟的方法,我们在不同的外电场值下对牛胰岛素蛋白质进行了一系列长达1个微妙的动力学模拟,在模拟中我们发现0.15V/nm的外电场并没有对胰岛素蛋白的二级结构产生较大的影响,而0.25V/nm及以上强度的外电场可以破坏胰岛素的二级结构,当外加电场强度在0.15~0.60V/nm时,蛋白质二级结构开始出现破坏所需要的模拟时间与外加电场的大小没有十分明确的相关性。一系列的结果显示在这个外加电场的范围内(0.15~0.60V/nm)蛋白质长时间暴露在外电场中的时间累积效应是导致蛋白质二级结构破坏的主要因素。当外加电场的值大于0.60V/nm时,外电场的强度对蛋白质内部氢键的寿命有着非常重要的影响。蛋白质内部氢键在外加电场为1.0V/nnm时出现了快速演化显示出外加电场能够加速蛋白质的折叠或者破坏蛋白质的结构。