[發明專利]通過構建CHARMM ROTAMERS力場預測突變氨基酸側鏈結構的方法

有效

• 申請號: 201810591726.5
• 申請日: 2018-06-08
• 公開（公告）號: CN109002690B
• 公開（公告）日: 2021-05-18
• 發明（設計）人: 姚建莊;李春雨;朱莉
• 申請（專利權）人: 濟南大學
• 主分類號: G16B20/50
• 分類號: G16B20/50;G16B5/00;G16B50/00
• 代理公司: 濟南泉城專利商標事務所 37218
• 代理人: 支文彬

1.一種通過構建CHARMM ROTAMERS力場預測突變氨基酸側鏈結構的方法，其特征在于，包括如下步驟：

a)下載并打開ROTAMERS庫中含有某氨基酸全部rotamers的pdb格式文件，將其中的N種旋轉異構體獨立保存并分別命名為氨基酸英文簡稱1-N；

b)使用VMD視圖軟件打開pdb格式文件，選用VPK模型，在VMD視圖軟件中測量氨基酸的鍵長Bonds、測量鍵角Angles、測量二面角Dihedrals，構建CHARMM ROTAMERS力場并得到N種旋轉異構體的穩定結構；

c)將測量完成的N種氨基酸旋轉異構體的數據添加到舊CHARMM力場top模板文件中所對應氨基酸的位置；

d)在WinSCP軟件中打開蛋白質文件并找到氨基酸所在位置，刪除氨基酸側鏈，將氨基酸名字末位字母改為1-N，在PuTTY軟件中得到氨基酸1-N的N個輸出文件fullprot.pdb；

e)使用VMD視圖軟件打開步驟d)中的fullprot.pdb文件，將得到的晶體結構與舊的CHARMM力場中對應的氨基酸結構和N種氨基酸旋轉異構體結構分別進行對比；

f)通過公式分別計算N種旋轉異構體的旋轉異構體與晶體偏轉數值RMSD(v,w)，式中v表示蛋白質晶體結構，w為步驟b)中氨基酸1-N或舊的CHARMM力場下該氨基酸的結構，x,y,z為蛋白質中每個重原子的坐標，n表示蛋白質側鏈結構的原子數，i為蛋白質中第i個重原子；

g)取步驟f)中數值最小的RMSD(v,w)為最接近蛋白質晶體結構的旋轉異構體；

h)將步驟g)中最接近蛋白質晶體結構的旋轉異構體、蛋白質晶體結構以及舊CHARMM力場下結構文件上傳到WinSCP軟件，通過PuTTY得到輸出結果e127_min.pdb文件并分別重新命名；

i)在PuTTY軟件中運行./charmmequ.inp指令進行分子動力學模擬，分別讓三種結構運行輸出，得到輸出文件e127_dyn0并分別重命名。

2.根據權利要求1所述的通過構建CHARMM ROTAMERS力場預測突變氨基酸側鏈結構的方法，其特征在于：步驟d)中在Putty軟件中輸入/charmmstep1_pdbreader.inpstep1_pdbreader.out指令運行charmm，得到氨基酸1-N的N個輸出文件fullprot.pdb。

3.根據權利要求1所述的通過構建CHARMM ROTAMERS力場預測突變氨基酸側鏈結構的方法，其特征在于：步驟h)中在PuTTY軟件中輸入vi solv1.inp命令進入insert模式，修改溶劑化中心后退出insert模式，依次輸入命令./charmmsolv1.inp和./charmmsolv2.inp進行溶劑化，得到已加上水球的solv3文件，再依次輸入【./charmmpatch.inp】、【./charmmsb1.inp】、【./charmmsb2.inp】、【./charmmmin.inp】四個指令得到輸出結果e127_min.pdb文件。