 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 5586 | 回復(fù): 15 | ||||||
[交流]
NAMD自學(xué)筆記：psf文件的生成以及conf文件的參數(shù)設(shè)定 已有12人參與
|
||||||
|
NAMD自學(xué)筆記 一、NAMD中psf文件的生成 最近需要用到NAMD做一些分子模擬，主要是自己的蛋白結(jié)構(gòu)的動力學(xué)模擬。當(dāng)時看了NAMD官網(wǎng)上的tutorial，但還是有些東西不是很清楚，經(jīng)過多次嘗試和請教本站的大俠們終于有所認(rèn)識，在此將一些目前教程中沒有涉及的內(nèi)容整理一下，期望給廣大象我這樣的初學(xué)者以幫助。 這里我們先說psf文件的生成： 通常pdb文件中包含水分子，教程上只是將蛋白提出來，而如果有其他金屬離子則按教程的做法這些離子將不再包含在pdb文件中，所有第一步要仔細(xì)檢查自己的目標(biāo)蛋白文件。 我的做法： 1、先用別的軟件將pdb文件中的水去掉。如果存在金屬離子為單獨(dú)的一條鏈，將金屬離子合并到蛋白鏈，并修改其原子序號。 2、檢查pdb文件中是否有非天然氨基酸：我的pdb文件中有幾個MSE，先將其改為MET。 3、金屬離子的立場文件中命名和常用的pdb文件不同：比如Ca在NAMD自帶的tutorial的立場文件top_all27_prot_lipid.inp中為CAL，因此將vi修改pdb文件中的CA變?yōu)镃AL。只用空格，pdb文件中各列的值有空格分隔，其詳細(xì)格式請看http://blog.163.com/nm_myc_1013/ ... 455201331842724812/。 4、搞定pdb文件后，修改下教程中的pgn文件， package require psfgen topology top_all27_prot_lipid.inp pdbalias residue HIS HSE pdbalias atom ILE CD1 CD segment U {pdb ubqp.pdb} coordpdb ubqp.pdb U guesscoord writepdb ubq.pdb writepsf ubq.psf 這里的ubqp 及輸出的ubq改為自己的文件。 注意中文的教程里面： “topology top all27_prot_lipid.inp”這一行有錯誤，應(yīng)該是topology top_all27_prot_lipid.inp。立場文件名有錯誤，（小生一開始沒發(fā)現(xiàn)，整了好多次才注意到）。 Source name.pgn 生成所需的pdb及psf文件。 如果錯誤可能是工作目錄不對，請CD到對應(yīng)目錄再重新運(yùn)行。 這里順便提一下NAMD的安裝和編譯，在下載NAMD的時候有個安裝教程里面的東西很多，試著按照那個readme安裝，老出現(xiàn)錯誤，后來發(fā)現(xiàn)里面有個smart-install.pl(忘了可能是smart-config.pl)可以根據(jù)提示進(jìn)行charmm的安裝和編譯。其他按說明來就可以。 二、蛋白的溶劑化及相關(guān)conf文件的參數(shù)設(shè)定。 這里我們主要講water box 按照NAMD提供的教程生成蛋白的溶劑化模型。得到protein_wb.pdb以及protein_wb.psf。 然后確定次立方體系的邊長 a b c以及中心xyz坐標(biāo)。 根據(jù)教程來就可以 set everyone [atomselect top all] measure minmax $everyone measure center $everyone 邊長有xyz坐標(biāo)的最大值和最小值計算得到。 比如我的系統(tǒng)中 邊長分別為 73.1 69.1 72.2 這個邊長在conf文中用到，這里是我的理解： 在conf文件中下面部分的參數(shù)用到： # Periodic Boundary Conditions cellBasisVector1 74.0 0.0 0.0 cellBasisVector2 0.0 70.0 0.0 cellBasisVector3 0.0 0 73.0 cellOrigin 24.3 59.8 32.6 其中 cellBasisVector123的值應(yīng)該比3個邊長稍大一點(diǎn)就可以。不知道自己的理解是否正確，但可以運(yùn)行，且結(jié)果也可以。 蛋白溶劑化后可以添加對抗離子，采用VMD-Extention-modeling-add ions 添加0.145mol/L的 NaCl。 后面附conf文件一些參數(shù)的說明和修改情況： # Integrator Parameters timestep 2.0 ;# 2fs/step rigidBonds all ;# needed for 2fs steps nonbondedFreq 1 fullElectFrequency 2 stepspercycle 10 #這個后面3個參數(shù)改大一些可以節(jié)省計算時間，但是我嘗試了一下修改為100 200 1000是發(fā)生錯誤，后來一想干脆就安教程中的數(shù)據(jù)來。可能改為10 20 100可以，但后面最小化的步驟必學(xué)是stepspercycle 的整數(shù)倍。 # PME (for full-system periodic electrostatics) PME yes PMEGridSpacing 1.0 #manual grid definition #PMEGridSizeX 75 #PMEGridSizeY 72 #PMEGridSizeZ 75 這里的自定義PMEGrid的大小，教程中說這三個數(shù)的值最好是2,3,5的指數(shù)倍。不是很懂，干脆就不自己定義了。 # Output outputName $outputname restartfreq 500 ;# 500steps = every 1ps dcdfreq 250 xstFreq 250 outputEnergies 100 outputPressure 100 這里輸出頻率問題很多：按照這個數(shù)值輸出文件尤其是.dcd文件很大，想改一下跑5ns將這些值改為10000，但結(jié)果中蛋白各原子的位置出現(xiàn)錯誤，好像就是立場錯誤的樣子。將其值改為1000后結(jié)果正常。 ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization minimize 1000 #最小化步驟自己寫了3000但發(fā)現(xiàn)好只跑到1100多就停止了后面的動力學(xué)模擬就不跑了。改為1000后正常。 reinitvels $temperature run 2500000 ;# 5ns 三、中斷的NAMD如何續(xù)跑： 修改.conf ############################################################# ## ADJUSTABLE PARAMETERS ## ############################################################# structure ../MD/tse_ion.psf coordinates ../MD/tse_ion.pdb set previou tse_ion_eq #添加的內(nèi)容 set current tse_ion_rest #添加的內(nèi)容 這2個可能并不需要 set temperature 310 set outputname tse_ion_rest #重新計算后輸出到新的文件中 bincoordinates ../MD/tse_ion_eq.restart.coor binvelocities ../MD/tse_ion_eq.restart.vel extendedSystem ../MD/tse_ion_eq.restart.xsc #上面3條為需要添加參數(shù)。 firsttimestep 301000 #restart from 301000 step ############################################################# ## SIMULATION PARAMETERS ## ############################################################# # Input paraTypeCharmm on parameters ../MD/par_all27_prot_lipid.inp # temperature $temperature 溫度不需要這里需要注釋掉即前面加個#號。 ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization # minimize 1000 # reinitvels $temperature 最小化也不需要 這2個參數(shù)都加上# run 2500000 ;# 5ns 現(xiàn)在我用的和遇到的就這么多�，F(xiàn)在基本上可以進(jìn)行NAMD的動力學(xué)模擬了，期望這些能讓一些人少走一些彎路。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 科研之寶，所向披靡 | MD simulation | NAMD |
» 本帖已獲得的紅花（最新10朵）
yoghurt117» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有3人回復(fù)
-
311求調(diào)劑
已經(jīng)有10人回復(fù)
-
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- NAMD跨膜環(huán)肽管求助，Periodic cell has become too small for original patch grid!
已經(jīng)有7人回復(fù)
- NAMD安裝中的nanotubes.conf文件如何編寫
已經(jīng)有3人回復(fù)
- NAMD分子體系平衡和分子動力學(xué)模擬
已經(jīng)有10人回復(fù)
- topol文件不識別小分子
已經(jīng)有8人回復(fù)
- 【求助】在NAMD運(yùn)行中使用CHARMM格式的力場文件，但力場參數(shù)來自AMBER力場，會出問題?/a>
已經(jīng)有4人回復(fù)
- 【求助】namd 出錯Periodic cell has become too small for original patch grid!
已經(jīng)有10人回復(fù)
- 【求助】只是優(yōu)化后文件，沒有restart文件，能繼續(xù)跑嗎
已經(jīng)有22人回復(fù)
- 【求助】有人模擬過電泳嗎
已經(jīng)有9人回復(fù)
- 【求助】namd和vmd就只能識別.pdb.psf.coor.vel等文件嗎
已經(jīng)有6人回復(fù)
- 【求助】NAMD安裝問題
已經(jīng)有19人回復(fù)
- 【在線答疑】關(guān)于NAMD/VMD
已經(jīng)有210人回復(fù)
蘭S水晶
銀蟲 (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 570.7
- 散金: 5
- 帖子: 122
- 在線: 226.1小時
- 蟲號: 2584828
- 注冊: 2013-08-07
- 專業(yè): 理論和計算化學(xué)
3樓2014-07-28 07:14:10
?}果
銀蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 365.5
- 散金: 33
- 帖子: 70
- 在線: 44.7小時
- 蟲號: 3671519
- 注冊: 2015-01-30
- 性別: MM
- 專業(yè): 化學(xué)動力學(xué)
11樓2018-04-20 15:53:51
| 感謝分享 |
2樓2014-07-27 15:55:43
 |
4樓2014-08-20 10:19:27
5樓2015-03-25 20:01:43
yoghurt117
木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 8018.3
- 散金: 56
- 紅花: 8
- 帖子: 576
- 在線: 134.9小時
- 蟲號: 761174
- 注冊: 2009-05-01
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
6樓2016-07-30 09:16:48
zhangduo2016
新蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1759.8
- 散金: 50
- 紅花: 1
- 帖子: 328
- 在線: 65.9小時
- 蟲號: 4479337
- 注冊: 2016-03-08
- 性別: MM
- 專業(yè): 藥物設(shè)計與藥物信息
7樓2017-01-25 11:16:29
zhaojingzsx
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 808
- 帖子: 20
- 在線: 13.7小時
- 蟲號: 2052302
- 注冊: 2012-10-10
- 專業(yè): 馬克思主義其他學(xué)科
8樓2017-05-09 16:25:30
你好哈-
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2649.9
- 散金: 13
- 紅花: 3
- 帖子: 127
- 在線: 78.9小時
- 蟲號: 3023909
- 注冊: 2014-03-06
- 性別: MM
- 專業(yè): 藥物設(shè)計與藥物信息
9樓2017-12-29 23:49:52
zhouziyun
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1.5
- 帖子: 1
- 在線: 1.9小時
- 蟲號: 4639331
- 注冊: 2016-04-28
- 性別: MM
- 專業(yè): 生物信息學(xué)
10樓2018-01-30 00:42:19
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 supertype 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實(shí)用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
|
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 化工學(xué)碩306求調(diào)劑 +9 | 42838695 2026-03-12 | 9/450 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
|
[考博] 2026年博士申請 +3 | QwQwQW10 2026-03-11 | 3/150 |
|