| 24小時(shí)熱門版塊排行榜 |
| 3 | 1/1 | 返回列表 |
| 查看: 1663 | 回復(fù): 2 | ||
little_lion木蟲(chóng) (小有名氣)
|
[求助]
蛋白質(zhì)中存在非標(biāo)準(zhǔn)殘基的話如何在Gromacs中生成力場(chǎng) 已有2人參與
|
| rt |

鐵蟲(chóng) (小有名氣)
|
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. 這個(gè)是Gromacs官網(wǎng)教程講的,我是按上面的步驟新添的非標(biāo)準(zhǔn)殘基。然后看你用的是什么力場(chǎng),我用的是Gromos96力場(chǎng),在該力場(chǎng)文件夾下找到.rtp文件,這個(gè)文件定義了所有標(biāo)準(zhǔn)殘基的信息,在這個(gè)文件的末尾手動(dòng)追加一個(gè)你自己的新殘基信息。至于新殘基你得先手動(dòng)畫一個(gè)你自己的殘基出來(lái)(有很多軟件可以畫,chem3D , gaussview5),畫完之后記得對(duì)這個(gè)殘基能量最小化一下。保存成PDB格式文件,丟到prodrg2網(wǎng)站上,它會(huì)直接生成出你要的拓?fù)湫畔,你要按照它生成出的拓(fù)湫畔榛鶞?zhǔn),手動(dòng)修改成Gromos力場(chǎng).rtp文件格式(你可以參照.rtp文件任何一個(gè)氨基酸來(lái)修改)。接著再去修改HDB文件,如果你有新的院子類型就得追加atomtypes.atp 文件,沒(méi)有就不用。最后,在residuetypes.dat上追加你新殘基的名字即可。我就是這樣手動(dòng)修改的力場(chǎng)文件。 樓上說(shuō)的acpype也可以用,支持Gromacs,我也用過(guò)。但是前提你得安裝Ambertools,裝了Ambertool就antechamber了,但是我裝起來(lái)有點(diǎn)麻煩,我是CENTOS系統(tǒng),差這個(gè)包差那個(gè)包,補(bǔ)了半天才裝上。還有一個(gè)問(wèn)題,acpype好像是針對(duì)全原子力場(chǎng)的,比如opls。而Gromos96是聯(lián)合力場(chǎng)。。。這個(gè)……我就不知道了。?茨隳M的時(shí)候,選的哪個(gè)力場(chǎng)了。。。 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 +6 | α____ 2026-03-22 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚(yú)ly 2026-03-21 | 11/550 |
|
|
[考研] 求調(diào)劑 +5 | 林之夕 2026-03-24 | 5/250 |
|
|
[考研]
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-24 | 5/250 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +4 | 曼111 2026-03-24 | 5/250 |
|
|
[考研] 材料與化工328分調(diào)劑 +4 | 。,。,。,。i 2026-03-23 | 4/200 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生(材料、化學(xué)、電化學(xué),環(huán)境) +3 | 我愛(ài)學(xué)電池 2026-03-23 | 3/150 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語(yǔ) 2026-03-20 | 3/150 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|