 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 2923 | 回復(fù): 6 | ||||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ，點(diǎn)擊這里進(jìn)行查看 | ||||
[交流]
gromacs自己添加residue遇到的問題 已有3人參與
|
||||
|
大家好，最近一直在研究gromacs的residue這塊，我根據(jù)prodrg給出的top file自己在.rtp中添加了新的residue，具體步驟如下，在選擇的force field文件夾中的.rtp中添加residue，在ffbonded.itp里增加新的鍵的類型，最后更新residuetypes.dat，我認(rèn)為添加的內(nèi)容都挺對(duì)的，但用pdb2gmx生成top文件時(shí)出現(xiàn)這樣的錯(cuò)誤，我想請(qǐng)問我少了那步呢，不知道還要在什么地方添加?xùn)|西呢？ Fatal error: atom N not found in buiding block 2MOL while combining tdb and rtp For more information and tips for troubleshooting, please check the GROMACS |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 高分子學(xué)習(xí)--qlh |
» 猜你喜歡
312求調(diào)劑
已經(jīng)有7人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有16人回復(fù)
-
280求調(diào)劑
已經(jīng)有3人回復(fù)
-
266求調(diào)劑
已經(jīng)有7人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有5人回復(fù)
-
環(huán)境工程調(diào)劑
已經(jīng)有8人回復(fù)
-
302求調(diào)劑
已經(jīng)有9人回復(fù)
-
297求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有4人回復(fù)
-
焦慮
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- gromacs，用pdb文件生成top文件時(shí)遇到的問題
已經(jīng)有4人回復(fù)
- topol文件不識(shí)別小分子
已經(jīng)有8人回復(fù)
- 【專家會(huì)診】bay 討論namd,gmx處理水溶液,蛋白等MD 模擬;討論水及水通道
已經(jīng)有50人回復(fù)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
chaizhm: 金幣+5, 謝謝~ 2013-01-20 10:03:31
zh1987hs: 金幣+15, 模擬EPI+1, jiao版真是好娃啊~哈哈 2013-01-20 10:38:01
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
chaizhm: 金幣+5, 謝謝~ 2013-01-20 10:03:31
zh1987hs: 金幣+15, 模擬EPI+1, jiao版真是好娃啊~哈哈 2013-01-20 10:38:01
|
你還是應(yīng)該先非常仔細(xì)的研究一下手冊(cè)，真的！ 你的意思就是定義了一個(gè)新的 分子但是gromacs里面沒有這個(gè)residue 吧： 下面我以charmm力場下添加 三氧化二鋁的模型為例子，說明如何添加： 1、使用MS建立滿足自己體系大小的三氧化二鋁模型的PDB文件 2、導(dǎo)入到VMD中建立PSF文件 3、通過三氧化二鋁的psf和pdb文件，轉(zhuǎn)化成gromacs的 itp 文件，保存成alo_pore_charged.itp，并且放到charmmm27.ff文件夾內(nèi)（部分內(nèi)容如下）： [moleculetype] ; molname nrexcl al2o3Nanopore 1 [atoms] ; id type resnr residu atom cgnr charge mass 1 OAL 1 ALO O8 1 -1.05 15.9994 2 OAL 1 ALO O9 2 -1.05 15.9994 3 OAL 1 ALO O10 3 -1.05 15.9994 4 AL 1 ALO AL7 4 1.601296 26.9815 5 OAL 1 ALO O2 5 -1.05 15.9994 6 OAL 1 ALO O13 6 -1.05 15.9994 7 OAL 1 ALO O14 7 -1.05 15.9994 8 OAL 1 ALO O15 8 -1.05 15.9994 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 20062 AL 1 ALO AL4 20062 1.601296 26.9815 20063 AL 1 ALO AL5 20063 1.601296 26.9815 20064 AL 1 ALO AL6 20064 1.601296 26.9815 [bonds] ; i j funct length force.c. [angles] ; i j k funct angle force.c. 這里具體的bond angle 信息根據(jù)psf文件中轉(zhuǎn)化，鍵長，鍵能自己在這里設(shè)定。 4、由于我這里的使用的原子名字，不是charmm力場的標(biāo)準(zhǔn)原子名，因此需要把這些添加到 charmm力場的 atomtypes.atp （charmm27.ff目錄內(nèi)）下面， ; atom types below are under my definition ALS 26.9815 ; the atoms on the AL2O3 pore's surface AL 26.9815 ; AL atoms in the AL2O3 OAL 15.9994 ; O atoms in the AL2O3 5、把建立好的pdb文件導(dǎo)入到gromacs中，生成。top文件就可以往下面繼續(xù)了： genbox -cp sys_box.pdb -cs spc216.gro -o sys_water.gro -p sys.top （.top文件可以手動(dòng)修改） 最后的sys.top如下： ; Include forcefield parameters #include "charmm27.ff/forcefield.itp" ; Include your topologies #include "charmm27.ff/alo_pore_charged.itp" ; Include chain topologies #include "dsdna45_DNA_chain_A.itp" #include "dsdna45_DNA_chain_B.itp" ; Include water topology #include "charmm27.ff/spc.itp" #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcx fcy fcz 1 1 1000 1000 1000 #endif ; Include topology for ions #include "charmm27.ff/ions.itp" [ system ] ; Name dsDNA in nanopore in water [ molecules ] ; Compound #mols DNA_chain_A 1 DNA_chain_B 1 al2o3Nanopore 1 SOL 88342 K 1803 CL 1847 注意注意注意： .top 中的 [ molecules ] 中的順序，一定要和坐標(biāo)文件（pdb 或者 gro）文件中的順序一致 ok思路就是這樣，具體的多看手冊(cè)，“atom N not found ”，我想主要就是第四步出現(xiàn)的問題吧 [ Last edited by jiaoyixiong on 2013-1-20 at 09:27 ] |
2樓2013-01-20 09:25:24
3樓2013-01-20 16:52:54
4樓2013-01-20 16:58:17
5樓2013-01-21 20:29:08
youngfi
木蟲 (小有名氣)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 2419.9
- 散金: 15
- 紅花: 1
- 帖子: 144
- 在線: 100.9小時(shí)
- 蟲號(hào): 1850824
- 注冊(cè): 2012-06-07
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
6樓2013-06-19 22:55:41
SunnyRain0.0
金蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1564.1
- 散金: 1
- 紅花: 7
- 帖子: 624
- 在線: 146.8小時(shí)
- 蟲號(hào): 1711069
- 注冊(cè): 2012-03-23
- 性別: MM
- 專業(yè): 中藥學(xué)其他科學(xué)問題
7樓2013-10-14 20:12:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ilxmc 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +3 | 陽陽哇塞 2026-03-14 | 7/350 |
|
|---|---|---|---|---|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 10/500 |
|
|
|
[考研] 281求調(diào)劑（0805） +3 | 煙汐憶海 2026-03-16 | 9/450 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
|
[考研] 工科，求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|