| 3 | 1/1 | 返回列表 |
| 查看: 1662 | 回復: 2 | ||
little_lion木蟲 (小有名氣)
|
[求助]
蛋白質(zhì)中存在非標準殘基的話如何在Gromacs中生成力場 已有2人參與
|
| rt |

鐵蟲 (小有名氣)
|
Adding a new residue If you have the need to introduce a new residue into an existing force field so that you can use pdb2gmx, or modify an existing one, there are several files you will need to modify. You must consult the manual for description of the required format. Follow these steps: Add the residue to the .rtp file for your chosen force field. You might be able to copy an existing residue, rename it and modify it suitably, or you may need to use an external topology generation tool and adapt the results to the .rtp format. If you need hydrogens to be able to be added to your residue, create an entry in the relevant .hdb file. If you are introducing new atom types, add them to the atomtypes.atp and ffnonbonded.itp files. If you require any new bonded types, add them to ffbonded.itp. Add your residue to residuetypes.dat with the appropriate specification (Protein, DNA, Ion, etc). If the residue involves special connectivity to other residues, update specbond.dat. To avoid an issue with grompp, in your new residue do not use the atom names O1, O2, and probably H1, H2 and H3 also. These might otherwise be picked up as atoms from possible peptide termini, which is likely not what you'd want. 這個是Gromacs官網(wǎng)教程講的,我是按上面的步驟新添的非標準殘基。然后看你用的是什么力場,我用的是Gromos96力場,在該力場文件夾下找到.rtp文件,這個文件定義了所有標準殘基的信息,在這個文件的末尾手動追加一個你自己的新殘基信息。至于新殘基你得先手動畫一個你自己的殘基出來(有很多軟件可以畫,chem3D , gaussview5),畫完之后記得對這個殘基能量最小化一下。保存成PDB格式文件,丟到prodrg2網(wǎng)站上,它會直接生成出你要的拓撲信息,你要按照它生成出的拓撲信息為基準,手動修改成Gromos力場.rtp文件格式(你可以參照.rtp文件任何一個氨基酸來修改)。接著再去修改HDB文件,如果你有新的院子類型就得追加atomtypes.atp 文件,沒有就不用。最后,在residuetypes.dat上追加你新殘基的名字即可。我就是這樣手動修改的力場文件。 樓上說的acpype也可以用,支持Gromacs,我也用過。但是前提你得安裝Ambertools,裝了Ambertool就antechamber了,但是我裝起來有點麻煩,我是CENTOS系統(tǒng),差這個包差那個包,補了半天才裝上。還有一個問題,acpype好像是針對全原子力場的,比如opls。而Gromos96是聯(lián)合力場。。。這個……我就不知道了。。看你模擬的時候,選的哪個力場了。。。 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 食品專碩 一志愿雙一流 328 +3 | xiaom99 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 070300化學求調(diào)劑 +9 | 苑豆豆 2026-03-20 | 9/450 |
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +5 | 谷雨上岸 2026-03-23 | 6/300 |
|
|
[考研] 305分求調(diào)劑(食品工程) +5 | Sxy112 2026-03-21 | 7/350 |
|
|
[考研] 求材料,環(huán)境專業(yè)調(diào)劑 +3 | 18567500178 2026-03-18 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
[考研] 一志愿國科過程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 0703化學求調(diào)劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 一志愿北京化工大學070300 學碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 材料學碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 考研調(diào)劑求學校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 中南大學化學學碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|