| 7 | 1/1 | 返回列表 |
| 查看: 1690 | 回復(fù): 6 | |||
gytjyb新蟲 (小有名氣)
|
[交流]
【求助】寫了個(gè)MOE的動(dòng)力學(xué)的過程,希望大家批評幾句! 已有3人參與
|
|
因?yàn)槲沂浅鯇W(xué)分子模擬,正好我我們學(xué)校有MOE的版權(quán)。 所以我試著用MOE為我的同源模型,跑了一個(gè)分子動(dòng)力學(xué)的優(yōu)化。請各位高手看看一下過程,我也不知道該怎么用moe,我想這是通常一般情況下的過程,懇請各位師兄姐,評價(jià)一下我有哪些做得不對?赡芪矣泻芏嗖粚Φ牡胤剑灿性S多不懂得地方,請?zhí)岢鰜怼Vx謝您了! 我是想構(gòu)建了120個(gè)氨基酸的蛋白模型,想對它做動(dòng)力學(xué)優(yōu)化。我不知道對不對,有好多我都不太明白。 首先,我給蛋白用AMBER99力場加氫,修正電荷。然后又給蛋白加了一圓形層的水,最大原子到邊界的厚度為7,是在ph7的時(shí)候加的。(我不知道怎么加離子,也還不會判斷是否電中性) 然后我用500的力量固定蛋白,先優(yōu)化外面的水到系統(tǒng)為RMS=1,然后整個(gè)放開系統(tǒng)優(yōu)化到能量最小化,為RMS=0.05。 然后我再開始跑動(dòng)力學(xué),因?yàn)槲蚁氲絼?dòng)力學(xué)先要吧外面的水跑好吧。所以,我先從50K升溫度到310K,用了50ps,選用NVT algorithrm用NPA(什么意思)。這時(shí)候在動(dòng)力學(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)容很多,但是沒有動(dòng)力學(xué)的列子! 謝謝大家,祝大家好運(yùn)! |
木蟲 (正式寫手)
新蟲 (小有名氣)
|
為我是初學(xué)分子模擬,正好我我們學(xué)校有MOE的版權(quán)。 所以我試著用MOE為我的同源模型,跑了一個(gè)分子動(dòng)力學(xué)的優(yōu)化。請各位高手看看一下過程,我也不知道該怎么用moe,我想這是通常一般情況下的過程,懇請各位師兄姐,評價(jià)一下我有哪些做得不對?赡芪矣泻芏嗖粚Φ牡胤剑灿性S多不懂得地方,請?zhí)岢鰜。謝謝您了! 我是想構(gòu)建了120個(gè)氨基酸的蛋白模型,想對它做動(dòng)力學(xué)優(yōu)化。我不知道對不對,有好多我都不太明白。 首先,我給蛋白用AMBER99力場加氫,修正電荷。然后又給蛋白加了一圓形層的水,最大原子到邊界的厚度為7,是在ph7的時(shí)候加的。(我不知道怎么加離子,也還不會判斷是否電中性) 然后我用500的力量固定蛋白,先優(yōu)化外面的水到系統(tǒng)為RMS=1,然后整個(gè)放開系統(tǒng)優(yōu)化到能量最小化,為RMS=0.05。 然后我再開始跑動(dòng)力學(xué),因?yàn)槲蚁氲絼?dòng)力學(xué)先要吧外面的水跑好吧。所以,我先從50K升溫度到310K,用了50ps,選用NVT algorithrm用NPA(什么意思)。這時(shí)候在動(dòng)力學(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)容很多,但是沒有動(dòng)力學(xué)的列子! 謝謝大家,祝大家好運(yùn)! 本文來自: 小木蟲論壇 http://www.gaoyang168.com/bbs/viewthread.php?tid=1880455&fpage=1 |
新蟲 (小有名氣)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +6 | Tyoumou 2026-03-18 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +11 | 妮妮ninicgb 2026-03-15 | 15/750 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 生物學(xué)調(diào)劑 +3 | Surekei 2026-03-21 | 3/150 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|