| 7 | 1/1 | 返回列表 |
| 查看: 1686 | 回復(fù): 6 | |||
gytjyb新蟲 (小有名氣)
|
[交流]
【求助】寫了個MOE的動力學(xué)的過程,希望大家批評幾句! 已有3人參與
|
|
因?yàn)槲沂浅鯇W(xué)分子模擬,正好我我們學(xué)校有MOE的版權(quán)。 所以我試著用MOE為我的同源模型,跑了一個分子動力學(xué)的優(yōu)化。請各位高手看看一下過程,我也不知道該怎么用moe,我想這是通常一般情況下的過程,懇請各位師兄姐,評價(jià)一下我有哪些做得不對?赡芪矣泻芏嗖粚Φ牡胤,也有許多不懂得地方,請?zhí)岢鰜。謝謝您了! 我是想構(gòu)建了120個氨基酸的蛋白模型,想對它做動力學(xué)優(yōu)化。我不知道對不對,有好多我都不太明白。 首先,我給蛋白用AMBER99力場加氫,修正電荷。然后又給蛋白加了一圓形層的水,最大原子到邊界的厚度為7,是在ph7的時候加的。(我不知道怎么加離子,也還不會判斷是否電中性) 然后我用500的力量固定蛋白,先優(yōu)化外面的水到系統(tǒng)為RMS=1,然后整個放開系統(tǒng)優(yōu)化到能量最小化,為RMS=0.05。 然后我再開始跑動力學(xué),因?yàn)槲蚁氲絼恿W(xué)先要吧外面的水跑好吧。所以,我先從50K升溫度到310K,用了50ps,選用NVT algorithrm用NPA(什么意思)。這時候在動力學(xué)選項(xiàng)里我選擇constrait all bond(這個選項(xiàng)默認(rèn)的意思除水以外都限制,有l(wèi)ights bond,all bond,none三種情況,我不知道這么做對不對,為什么這樣做),有關(guān)水的那項(xiàng),我選擇用rigid(不知道什么意思),然后用50PS跑跑升溫的過程。 升溫完成以后,我繼續(xù)上次的模擬,我用1個NS,這時我想應(yīng)該是否選擇constrait none 什么都不限制了(為什么,這樣對嗎),維持300K, 其它都保持不變,一直到整個過程完,取最后200PS的平均構(gòu)像。 ( 然后我發(fā)現(xiàn)一個問題,不知道怎樣去看整個過程RMSD的變化。) 我就想用MOE做個模型優(yōu)化,我不知道可否這樣做!MOE的內(nèi)容很多,但是沒有動力學(xué)的列子! 謝謝大家,祝大家好運(yùn)! |
木蟲 (正式寫手)
新蟲 (小有名氣)
|
為我是初學(xué)分子模擬,正好我我們學(xué)校有MOE的版權(quán)。 所以我試著用MOE為我的同源模型,跑了一個分子動力學(xué)的優(yōu)化。請各位高手看看一下過程,我也不知道該怎么用moe,我想這是通常一般情況下的過程,懇請各位師兄姐,評價(jià)一下我有哪些做得不對。可能我有很多不對的地方,也有許多不懂得地方,請?zhí)岢鰜。謝謝您了! 我是想構(gòu)建了120個氨基酸的蛋白模型,想對它做動力學(xué)優(yōu)化。我不知道對不對,有好多我都不太明白。 首先,我給蛋白用AMBER99力場加氫,修正電荷。然后又給蛋白加了一圓形層的水,最大原子到邊界的厚度為7,是在ph7的時候加的。(我不知道怎么加離子,也還不會判斷是否電中性) 然后我用500的力量固定蛋白,先優(yōu)化外面的水到系統(tǒng)為RMS=1,然后整個放開系統(tǒng)優(yōu)化到能量最小化,為RMS=0.05。 然后我再開始跑動力學(xué),因?yàn)槲蚁氲絼恿W(xué)先要吧外面的水跑好吧。所以,我先從50K升溫度到310K,用了50ps,選用NVT algorithrm用NPA(什么意思)。這時候在動力學(xué)選項(xiàng)里我選擇constrait all bond(這個選項(xiàng)默認(rèn)的意思除水以外都限制,有l(wèi)ights bond,all bond,none三種情況,我不知道這么做對不對,為什么這樣做),有關(guān)水的那項(xiàng),我選擇用rigid(不知道什么意思),然后用50PS跑跑升溫的過程。 升溫完成以后,我繼續(xù)上次的模擬,我用1個NS,這時我想應(yīng)該是否選擇constrait none 什么都不限制了(為什么,這樣對嗎),維持300K, 其它都保持不變,一直到整個過程完,取最后200PS的平均構(gòu)像。 ( 然后我發(fā)現(xiàn)一個問題,不知道怎樣去看整個過程RMSD的變化。) 我就想用MOE做個模型優(yōu)化,我不知道可否這樣做!MOE的內(nèi)容很多,但是沒有動力學(xué)的列子! 謝謝大家,祝大家好運(yùn)! 本文來自: 小木蟲論壇 http://www.gaoyang168.com/bbs/viewthread.php?tid=1880455&fpage=1 |
新蟲 (小有名氣)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +11 | 妮妮ninicgb 2026-03-15 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 本人考085602 化學(xué)工程 專碩 +20 | 不知道叫什么! 2026-03-15 | 22/1100 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | yzysaa 2026-03-21 | 4/200 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +3 | zxts12 2026-03-21 | 3/150 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-17 | 8/400 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|