| 5 | 1/1 | 返回列表 |
| 查看: 1699 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
gytjyb新蟲 (小有名氣)
|
[交流]
【求助】寫了個(gè)MOE的動力學(xué)的過程,希望大家批評幾句! 已有3人參與
|
||
|
因?yàn)槲沂浅鯇W(xué)分子模擬,正好我我們學(xué)校有MOE的版權(quán)。 所以我試著用MOE為我的同源模型,跑了一個(gè)分子動力學(xué)的優(yōu)化。請各位高手看看一下過程,我也不知道該怎么用moe,我想這是通常一般情況下的過程,懇請各位師兄姐,評價(jià)一下我有哪些做得不對。可能我有很多不對的地方,也有許多不懂得地方,請?zhí)岢鰜。謝謝您了! 我是想構(gòu)建了120個(gè)氨基酸的蛋白模型,想對它做動力學(xué)優(yōu)化。我不知道對不對,有好多我都不太明白。 首先,我給蛋白用AMBER99力場加氫,修正電荷。然后又給蛋白加了一圓形層的水,最大原子到邊界的厚度為7,是在ph7的時(shí)候加的。(我不知道怎么加離子,也還不會判斷是否電中性) 然后我用500的力量固定蛋白,先優(yōu)化外面的水到系統(tǒng)為RMS=1,然后整個(gè)放開系統(tǒng)優(yōu)化到能量最小化,為RMS=0.05。 然后我再開始跑動力學(xué),因?yàn)槲蚁氲絼恿W(xué)先要吧外面的水跑好吧。所以,我先從50K升溫度到310K,用了50ps,選用NVT algorithrm用NPA(什么意思)。這時(shí)候在動力學(xué)選項(xiàng)里我選擇constrait all bond(這個(gè)選項(xiàng)默認(rèn)的意思除水以外都限制,有l(wèi)ights bond,all bond,none三種情況,我不知道這么做對不對,為什么這樣做),有關(guān)水的那項(xiàng),我選擇用rigid(不知道什么意思),然后用50PS跑跑升溫的過程。 升溫完成以后,我繼續(xù)上次的模擬,我用1個(gè)NS,這時(shí)我想應(yīng)該是否選擇constrait none 什么都不限制了(為什么,這樣對嗎),維持300K, 其它都保持不變,一直到整個(gè)過程完,取最后200PS的平均構(gòu)像。 ( 然后我發(fā)現(xiàn)一個(gè)問題,不知道怎樣去看整個(gè)過程RMSD的變化。) 我就想用MOE做個(gè)模型優(yōu)化,我不知道可否這樣做!MOE的內(nèi)容很多,但是沒有動力學(xué)的列子! 謝謝大家,祝大家好運(yùn)! |
新蟲 (小有名氣)
木蟲 (正式寫手)
新蟲 (小有名氣)
|
為我是初學(xué)分子模擬,正好我我們學(xué)校有MOE的版權(quán)。 所以我試著用MOE為我的同源模型,跑了一個(gè)分子動力學(xué)的優(yōu)化。請各位高手看看一下過程,我也不知道該怎么用moe,我想這是通常一般情況下的過程,懇請各位師兄姐,評價(jià)一下我有哪些做得不對?赡芪矣泻芏嗖粚Φ牡胤,也有許多不懂得地方,請?zhí)岢鰜。謝謝您了! 我是想構(gòu)建了120個(gè)氨基酸的蛋白模型,想對它做動力學(xué)優(yōu)化。我不知道對不對,有好多我都不太明白。 首先,我給蛋白用AMBER99力場加氫,修正電荷。然后又給蛋白加了一圓形層的水,最大原子到邊界的厚度為7,是在ph7的時(shí)候加的。(我不知道怎么加離子,也還不會判斷是否電中性) 然后我用500的力量固定蛋白,先優(yōu)化外面的水到系統(tǒng)為RMS=1,然后整個(gè)放開系統(tǒng)優(yōu)化到能量最小化,為RMS=0.05。 然后我再開始跑動力學(xué),因?yàn)槲蚁氲絼恿W(xué)先要吧外面的水跑好吧。所以,我先從50K升溫度到310K,用了50ps,選用NVT algorithrm用NPA(什么意思)。這時(shí)候在動力學(xué)選項(xiàng)里我選擇constrait all bond(這個(gè)選項(xiàng)默認(rèn)的意思除水以外都限制,有l(wèi)ights bond,all bond,none三種情況,我不知道這么做對不對,為什么這樣做),有關(guān)水的那項(xiàng),我選擇用rigid(不知道什么意思),然后用50PS跑跑升溫的過程。 升溫完成以后,我繼續(xù)上次的模擬,我用1個(gè)NS,這時(shí)我想應(yīng)該是否選擇constrait none 什么都不限制了(為什么,這樣對嗎),維持300K, 其它都保持不變,一直到整個(gè)過程完,取最后200PS的平均構(gòu)像。 ( 然后我發(fā)現(xiàn)一個(gè)問題,不知道怎樣去看整個(gè)過程RMSD的變化。) 我就想用MOE做個(gè)模型優(yōu)化,我不知道可否這樣做!MOE的內(nèi)容很多,但是沒有動力學(xué)的列子! 謝謝大家,祝大家好運(yùn)! 本文來自: 小木蟲論壇 http://www.gaoyang168.com/bbs/viewthread.php?tid=1880455&fpage=1 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +9 | 張zz111 2026-03-27 | 10/500 |
|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 6/300 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +8 | 舉個(gè)栗子oi 2026-03-24 | 9/450 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|