 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 3471 | 回復(fù): 5 | |||||||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)4次，作者cndenis增加金幣 3 個(gè) | |||||||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ，點(diǎn)擊這里進(jìn)行查看 | |||||||
[資源]
使用Amber進(jìn)行QM/MM動(dòng)力學(xué)模擬
|
|||||||
|
使用Amber進(jìn)行QM/MM動(dòng)力學(xué)模擬 cnDenis，2012年5月11日 分子動(dòng)力學(xué)軟件Amber可以進(jìn)行QM/MM混合模擬，所得到的軌跡也可以用QM/MM混合的方法計(jì)算自由能。 由于分子力學(xué)（MM）方法是將原子視為一個(gè)整體，不考慮核外的電子分布的影響，量子力學(xué)（QM）方法考慮了核外電子，看起來給人感覺要高級(jí)和準(zhǔn)確一些。但并不總是如此，在構(gòu)像分析的一些性質(zhì)中，用合適的力子力場(chǎng)所得到的結(jié)果可能會(huì)比用半經(jīng)驗(yàn)量化方法做出來的更好。而且QM比MM昂貴得多，一般而言體系中QM部分原子不要超過80-100個(gè)，否則運(yùn)行起來會(huì)非常慢。 Amber中內(nèi)置的QM模塊源于Mopac的代碼，支持 PM3，AM1，RM1，MNDO，PDDG/PM3，PDDG/MNDO， PM3CARB1和DFTB/SCC-DFTB等半經(jīng)驗(yàn)量化方法。具體支持的元素類型請(qǐng)參考AmberTool手冊(cè)。默認(rèn)的方法為PM3。 Amber手冊(cè)使用QM/MM模擬的sander輸入文件例子如下(括號(hào)中的只是注釋，實(shí)際應(yīng)用時(shí)要去掉）: &cntrl imin=0, nstlim=10000, (perform MD for 10,000 steps) dt=0.002, (2 fs time step) ntt=1, tempi=0.1, temp0=300.0 (Berendsen temperature control) ntb=1, (Constant volume periodic boundaries) ntf=2, ntc=2, (Shake hydrogen atoms) cut=8.0, (8 angstrom classical non-bond cut off) ifqnt=1 (Switch on QM/MM coupled potential) / &qmmm qmmask=':753' (Residue 753 should be treated using QM) qmcharge=-2, (Charge on QM region is -2) qm_theory='PM3', (Use the PM3 semi-empirical Hamiltonian) qmcut=8.0 (Use 8 angstrom cut off for QM region) / 在&cntrl中設(shè)置ifqnt=1表示使用QM/MM模擬，然后添加&qmmm段，設(shè)置QM的各種參數(shù)。 qmmask 是設(shè)置使用QM的原子，冒號(hào)后面表示殘基序號(hào)。還有其它的表示方法，語法請(qǐng)參考Amber手冊(cè)ambmask部分。注意如果需要從一個(gè)分子中選取部分原子的話，QM原子和MM原子之間最好是非極性鍵，如C-C鍵。不可以把相連的C-H原子分別放在兩個(gè)區(qū)。 qmcharge 是QM部分所帶的電荷。 qm_theory 是QM計(jì)算使用的方法，默認(rèn)值是PM3，如上所述。 qmcut 是QM原子受MM原子影響的范圍，注意這個(gè)值是對(duì)于QM原子整體而言的，如果一個(gè)MM原子與任何一個(gè)QM原子的距離小于qmcut，所有的QM原子都會(huì)受到它影響。對(duì)于GB隱式水體系，設(shè)15~20埃，對(duì)于周期性顯式水體系，一般8~9埃就可了。 目前只有sander支持QM/MM，PMEMD還不支持，所以CUDA加速是沒法用的。QM/MM可以使用MPI并行加速，但由于并不是所有代碼都支持并行化，手冊(cè)上說超過8個(gè)CPU加速效果就不大了。 動(dòng)力學(xué)軌跡也可以使用QM/MM處理。以下是我用MMPBSA.py計(jì)算自由能的一個(gè)輸入文件。 calculate dG with QM/MM &general verbose=2, / &gb ifqnt=1, qm_residues=1, qm_theory='PM3', / 與普通純MM的計(jì)算相比，輸入文件在&gb段多了ifqnt、qm_residue和qm_theory，要注意的是qm_residues的表示方法與sander輸入文件中不同，在這里不能使用Mask，而是直接寫上殘基序號(hào)，如果有多個(gè)殘基，要用分號(hào)分開。 目前Amber中的QM/MM不支持PB模型，所以MMPBSA.py中的輸入只能有g(shù)b段，不能有pb段，也就是只能算MM/GBSA，不能算MM/PBSA。 一般來說，模擬時(shí)和自由能計(jì)算的QM部分設(shè)置要一致，否則可能結(jié)果會(huì)很不合理。 [ Last edited by cndenis on 2012-5-11 at 21:59 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 收藏 | amber | MD/QM/MM | 吸附 |
| QM/MM |
» 猜你喜歡
291求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑
已經(jīng)有8人回復(fù)
-
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷
已經(jīng)有12人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有6人回復(fù)
-
考研化學(xué)308分求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑
已經(jīng)有5人回復(fù)
-
工科0856求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國石油大學(xué)（華東） 本科齊魯工業(yè)大學(xué)
已經(jīng)有4人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有16人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有5人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 關(guān)于amber動(dòng)力學(xué)模擬中MM/GBSA的問題
已經(jīng)有7人回復(fù)
- 分子模擬時(shí)間越長(zhǎng)越好嗎？
已經(jīng)有44人回復(fù)
- AMBER中qm/mm的問題
已經(jīng)有5人回復(fù)
- [求助]QM/MM聯(lián)用
已經(jīng)有11人回復(fù)
3樓2012-05-13 16:49:09
2樓2012-05-12 15:35:36
4樓2013-12-11 09:53:26
6樓2017-10-27 06:30:10
簡(jiǎn)單回復(fù)
2013-12-18 16:12
回復(fù)
五星好評(píng) 
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cndenis 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研化學(xué)308分求調(diào)劑 +6 | 你好明天你好 2026-03-23 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 316求調(diào)劑 +7 | 梁茜雯 2026-03-19 | 7/350 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +6 | 墨墨漠 2026-03-23 | 7/350 |
|
|
|
[考研] 306求調(diào)劑 +9 | chuanzhu川燭 2026-03-18 | 9/450 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛生物生物愛?/a> 2026-03-17 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|