amber中非標準氨基酸殘基的引數生成

2021-07-25 13:55:27 字數 1112 閱讀 1531

在amber教程中給出了相關的教程,具體流程都是根據這個教程來製作的。

但是他的電荷使用的是bcc電荷,但是並沒有介紹如何把電荷轉換為resp。

流程:2.antechamber -fi ccif -i cro.cif -bk cro -fo ac -o cro_step1.ac -at amber #和教程上一致,目的為了生成參照的ac檔案,-at 必須制定,否則原子型別預設gaff;

3.antechamber -fi gout -i cro.log -rn cro -fo ac -o cro_step2.ac -c resp -at amber #賦予resp電荷到新的ac檔案中;

4.文字操作將cro_step2.ac中的座標以及電荷資訊列覆蓋到cro_step1.ac中,同時修改主鏈n原子的原子型別從nt到n,儲存為cro_step3.ac;

5.製作mc檔案,格式如下圖,參照cro_step3.ac檔案製作;

6.prepgen -i cro_step3.ac -o cro.prepin -m cro.mc -rn cro #得到主鏈資訊檔案;

7.parmchk2 -i cro.prepin -f prepi -p frcmod.cro -a y -p $amberhome/dat/leap/parm/parm10.dat #檢查生成的引數資訊,檢視frcmod檔案,若其中出現某一行引數為0且有 "attn,need revision" 標識則說明在parm10.dat中引數不完整需要用gaff進行補充;

8.grep -v "attn" frcmod.cro > frcmod1.cro #剔除不包含的引數,其餘存成引數檔案;

9.parmchk2 -i cro.prepin -f prepi -o frcmod2.cro #剩餘的引數利用gaff生成引數,由於命令中沒有指定 "-a y" 則只生成了原力場中未包含的引數檔案;

10.在leap中匯入檔案即可,loadamberprep cro.prepin ,loadamberparams frcmod1.cro ,loadamberparams frcmod2.cro ,對蛋白無需額外處理,(需要去除h,pdb4amber)。

非標準程式設計能力

一切的一切都基於我的個人經驗,歡迎拍磚和補充。想從coder稱為programmer嗎?想讓工資成幾何數的增長嗎?想。以下題目可以不完全地測試您系統程式設計功力,如果給了您丁點啟發或者幫助都是我的榮幸。一。windows和 nix作業系統的區別在 tips 答案不包括圖形使用者介面。二。什麼是記憶體...

非標準路徑的 CMAKE ROOT

由於系統本身的 cmake 版本太低,需要更新乙個。沒有管理員全限,只能在 home下安裝軟體。修改 bashrc,增加 export path home local bin path export ld library path home local lib ld library path 然後編...

BeautifulSoup解析非標準HTML的問題

beautifulsoup版本 4.3.2 在用beautifulsoup.find all 搜尋html時,遇到下面的 a href shipin donghuapian 2012 07 25 23404.html title 謙謙君子 target blank 溫潤如玉 a 可以看出 中a標籤的...