本发明涉及蛋白质工程和分子动力学模拟技术领域，具体涉及一种通过构建charmmrotamers力场预测突变氨基酸侧链结构的方法。
背景技术
利用传统的蛋白质工程实验方法改造蛋白质，耗时费力，可行性不高。计算机运算速度的快速发展使基于力场的分子动力学模拟变成现实。
charmm(chemistryatharvardmacromolecularmechanics)是分子动力学场集合的名称，以及与他们相关的分子动力学模拟和分析计算机软件包的名称。它作为引领软件，对蛋白质分子动力学模拟的发展起到了至关重要的作用，但是在预测突变氨基酸的侧链结构时，charmm力场仍有不足之处。
蛋白质的结构复杂多样，基本组成单位是氨基酸，各个氨基酸具有不同rotamers(rotamer是仅由绕分子內单键的内旋转构象不同而产生的旋转异构体，在实际的晶体结构里面有许多rotamers)的三维结构，主链是重合的，侧链在不同的位置。当把某一氨基酸突变为目标氨基酸的时候，因为charmm力场本身只有氨基酸的一种rotamer，不包括其他rotamers，所以charmm力场不可以准确快速的确定目标氨基酸侧链的位置。
技术实现要素：
本发明为了克服以上技术的不足，提供了一种使新charmm力场能够更好的预测突变氨基酸的侧链结构，并通过能量最小化和分子动力学模拟，得出结合能最小的突变体结构，与已知突变体的晶体结构进行对比，验证新charmm力场的正确性和优越性的预测突变氨基酸侧链结构的方法。
本发明克服其技术问题所采用的技术方案是：
一种通过构建charmmrotamers力场预测突变氨基酸侧链结构的方法，包括如下步骤：
a)下载并打开rotamers库中含有某氨基酸全部rotamers的pdb格式文件，将其中的n种旋转异构体独立保存并分别命名为氨基酸英文简称1-n；
b)使用vmd视图软件打开pdb格式文件，选用vpk模型，在vmd视图软件中测量氨基酸的键长bonds、测量键角angles、测量二面角dihedrals，构建rotamercharmm力场并得到n种氨基酸稳定的旋转异构体结构；
c)将测量完成的n种氨基酸旋转异构体的数据添加到旧charmm力场top模板文件中所对应氨基酸的位置；
d)在winscp软件中打开蛋白质文件并找到氨基酸所在位置，删除氨基酸侧链，将氨基酸名字末位字母改为1-n，在putty软件中得到氨基酸1-n的n个输出文件fullprot.pdb；
e)使用vmd视图软件打开步骤d)中的fullprot.pdb文件，将得到的晶体结构与旧的charmm力场中对应的氨基酸结构和n种氨基酸旋转异构体结构分别进行对比；
f)通过公式
g)取步骤f)中数值最小的rmsd(v,w)为最接近蛋白质晶体结构的旋转异构体；
h)将步骤g)中最接近蛋白质晶体结构的旋转异构体、蛋白质晶体结构以及旧charmm力场下结构文件上传到winscp软件，通过putty得到输出结果e127_min.pdb文件并分别重新命名；
i)在putty软件中运行./charmm<equ.inp指令进行分子动力学模拟，分别让三种结构运行输出，得到输出文件e127_dyn0并分别重命名。
进一步的，步骤d)中在putty软件中输入/charmm<step1_pdbreader.inp>step1_pdbreader.out指令运行charmm，得到氨基酸1-n的n个输出文件fullprot.pdb。
进一步的，步骤h)中在putty软件中输入visolv1.inp命令进入insert模式，修改溶剂化中心后退出insert模式，依次输入命令./charmm<solv1.inp和./charmm<solv2.inp进行溶剂化，得到已加上水球的solv3文件，再依次输入【./charmm<patch.inp】、【./charmm<sb1.inp】、【./charmm<sb2.inp】、【./charmm<min.inp】四个指令得到输出结果e127_min.pdb文件。
本发明的有益效果是：通过运用vmd图形显示软件及putty与winscp运行软件的结合，进行结构对比，我们找到了氨基酸不同旋转异构体中与晶体结构最相近且与旧charmm力场区分明显的稳定结构。通过能量最小化和分子动力学模拟研究，进行进一步测量与计算得出rmsd值，通过比较，进而证实了所得结构的准确性。我们的研究结果在数据准确丰富的基础上，快速的找到了氨基酸中能量最低、结构稳定的旋转异构体，节约了实验时间和成本，为蛋白质药物设计提供了可靠依据。
具体实施方式
下面对本发明做进一步说明。
一种通过构建charmmrotamers力场预测突变氨基酸侧链结构的方法，包括如下步骤：
a)下载并打开rotamers库中含有某氨基酸全部rotamers的pdb格式文件，将其中的n种旋转异构体独立保存并分别命名为氨基酸英文简称1-n，如酪氨酸tyr即为ty1-4。
b)使用vmd视图软件打开pdb格式文件，选用vpk模型，此时可以观察到氨基酸的结构是由主链和侧链两部分组成，在vmd视图软件中测量氨基酸的键长bonds、测量键角angles、测量二面角dihedrals，构建rotamercharmm力场并得到n种氨基酸稳定的旋转异构体结构。在vmd视图软件中，atoms表示原子，可以知道每个原子的命名；bonds表示键长，按顺序单击两个原子可测出其键长；angles表示键角，此选项下依次单击三个原子即可测出其键角；dihedrals表示二面角，依次单击选中四个原子即可测量出其二面角的角度。用vmd测量数据时每个ic需要四个连接的原子，它们有两种连接方式，正常模型需要测出i-j与k-l的键长、i-j-k与j-k-l的键角、i-j-k-l的二面角；非正常模型需要测出i-k*与k*-l的键长、i-k*-j与l-k*-j的键角、i-j-k*-l的二面角。
c)通过运用vmd软件对氨基酸n个旋转异构体进行显示，测量各有效原子之间的键长、键角、二面角，将这些数据记录分析，可以快速找到氨基酸侧链的变化，进而方便找到氨基酸稳定的旋转异构体结构。将测量完成的n种氨基酸旋转异构体的数据添加到旧charmm力场top模板文件中所对应氨基酸的位置。
d)在winscp软件中打开蛋白质文件并找到氨基酸所在位置，删除氨基酸侧链，将氨基酸名字末位字母改为1-n，在putty软件中得到氨基酸1-n的n个输出文件fullprot.pdb；
e)使用vmd视图软件打开步骤d)中的fullprot.pdb文件，将得到的晶体结构与旧的charmm力场中对应的氨基酸结构和n种氨基酸旋转异构体结构分别进行对比。
f)通过公式
g)取步骤f)中数值最小的rmsd(v,w)为最接近蛋白质晶体结构的旋转异构体。
h)将步骤g)中最接近蛋白质晶体结构的旋转异构体、蛋白质晶体结构以及旧charmm力场下结构文件上传到winscp软件，通过putty得到输出结果e127_min.pdb文件并分别重新命名。
i)在putty软件中运行./charmm<equ.inp指令进行分子动力学模拟，分别让三种结构运行输出，得到输出文件e127_dyn0并分别重命名。
补充说明上述步骤结合取得的效果。
通过运用vmd图形显示软件及putty与winscp运行软件的结合，进行结构对比，我们找到了氨基酸不同旋转异构体中与晶体结构最相近且与旧charmm力场区分明显的稳定结构。通过能量最小化和分子动力学模拟研究，进行进一步测量与计算得出rmsd值，通过比较，进而证实了所得结构的准确性。我们的研究结果在数据准确丰富的基础上，快速的找到了氨基酸中能量最低、结构稳定的旋转异构体，节约了实验时间和成本，为蛋白质药物设计提供了可靠依据。具体的，步骤d)中在putty软件中输入/charmm<step1_pdbreader.inp>step1_pdbreader.out指令运行charmm，得到氨基酸1-n的n个输出文件fullprot.pdb。
具体的，步骤h)中在putty软件中输入visolv1.inp命令进入insert模式，修改溶剂化中心后退出insert模式，依次输入命令./charmm<solv1.inp和./charmm<solv2.inp进行溶剂化，得到已加上水球的solv3文件，再依次输入【./charmm<patch.inp】、【./charmm<sb1.inp】、【./charmm<sb2.inp】、【./charmm<min.inp】四个指令得到输出结果e127_min.pdb文件。