| 5 | 1/1 | 返回列表 |
| 查看: 1701 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
gytjyb新蟲 (小有名氣)
|
[交流]
【求助】寫了個MOE的動力學(xué)的過程,希望大家批評幾句! 已有3人參與
|
||
|
因為我是初學(xué)分子模擬,正好我我們學(xué)校有MOE的版權(quán)。 所以我試著用MOE為我的同源模型,跑了一個分子動力學(xué)的優(yōu)化。請各位高手看看一下過程,我也不知道該怎么用moe,我想這是通常一般情況下的過程,懇請各位師兄姐,評價一下我有哪些做得不對。可能我有很多不對的地方,也有許多不懂得地方,請?zhí)岢鰜怼Vx謝您了! 我是想構(gòu)建了120個氨基酸的蛋白模型,想對它做動力學(xué)優(yōu)化。我不知道對不對,有好多我都不太明白。 首先,我給蛋白用AMBER99力場加氫,修正電荷。然后又給蛋白加了一圓形層的水,最大原子到邊界的厚度為7,是在ph7的時候加的。(我不知道怎么加離子,也還不會判斷是否電中性) 然后我用500的力量固定蛋白,先優(yōu)化外面的水到系統(tǒng)為RMS=1,然后整個放開系統(tǒng)優(yōu)化到能量最小化,為RMS=0.05。 然后我再開始跑動力學(xué),因為我想到動力學(xué)先要吧外面的水跑好吧。所以,我先從50K升溫度到310K,用了50ps,選用NVT algorithrm用NPA(什么意思)。這時候在動力學(xué)選項里我選擇constrait all bond(這個選項默認(rèn)的意思除水以外都限制,有l(wèi)ights bond,all bond,none三種情況,我不知道這么做對不對,為什么這樣做),有關(guān)水的那項,我選擇用rigid(不知道什么意思),然后用50PS跑跑升溫的過程。 升溫完成以后,我繼續(xù)上次的模擬,我用1個NS,這時我想應(yīng)該是否選擇constrait none 什么都不限制了(為什么,這樣對嗎),維持300K, 其它都保持不變,一直到整個過程完,取最后200PS的平均構(gòu)像。 ( 然后我發(fā)現(xiàn)一個問題,不知道怎樣去看整個過程RMSD的變化。) 我就想用MOE做個模型優(yōu)化,我不知道可否這樣做!MOE的內(nèi)容很多,但是沒有動力學(xué)的列子! 謝謝大家,祝大家好運! |
木蟲 (正式寫手)
新蟲 (小有名氣)
|
為我是初學(xué)分子模擬,正好我我們學(xué)校有MOE的版權(quán)。 所以我試著用MOE為我的同源模型,跑了一個分子動力學(xué)的優(yōu)化。請各位高手看看一下過程,我也不知道該怎么用moe,我想這是通常一般情況下的過程,懇請各位師兄姐,評價一下我有哪些做得不對?赡芪矣泻芏嗖粚Φ牡胤剑灿性S多不懂得地方,請?zhí)岢鰜怼Vx謝您了! 我是想構(gòu)建了120個氨基酸的蛋白模型,想對它做動力學(xué)優(yōu)化。我不知道對不對,有好多我都不太明白。 首先,我給蛋白用AMBER99力場加氫,修正電荷。然后又給蛋白加了一圓形層的水,最大原子到邊界的厚度為7,是在ph7的時候加的。(我不知道怎么加離子,也還不會判斷是否電中性) 然后我用500的力量固定蛋白,先優(yōu)化外面的水到系統(tǒng)為RMS=1,然后整個放開系統(tǒng)優(yōu)化到能量最小化,為RMS=0.05。 然后我再開始跑動力學(xué),因為我想到動力學(xué)先要吧外面的水跑好吧。所以,我先從50K升溫度到310K,用了50ps,選用NVT algorithrm用NPA(什么意思)。這時候在動力學(xué)選項里我選擇constrait all bond(這個選項默認(rèn)的意思除水以外都限制,有l(wèi)ights bond,all bond,none三種情況,我不知道這么做對不對,為什么這樣做),有關(guān)水的那項,我選擇用rigid(不知道什么意思),然后用50PS跑跑升溫的過程。 升溫完成以后,我繼續(xù)上次的模擬,我用1個NS,這時我想應(yīng)該是否選擇constrait none 什么都不限制了(為什么,這樣對嗎),維持300K, 其它都保持不變,一直到整個過程完,取最后200PS的平均構(gòu)像。 ( 然后我發(fā)現(xiàn)一個問題,不知道怎樣去看整個過程RMSD的變化。) 我就想用MOE做個模型優(yōu)化,我不知道可否這樣做!MOE的內(nèi)容很多,但是沒有動力學(xué)的列子! 謝謝大家,祝大家好運! 本文來自: 小木蟲論壇 http://www.gaoyang168.com/bbs/viewthread.php?tid=1880455&fpage=1 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【求調(diào)劑】新能源材料本科,一志愿211,初試321 +4 | 求調(diào)劑學(xué)校, 2026-04-02 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿安徽大學(xué)計算機科學(xué)與技術(shù)學(xué)碩,331分求調(diào)劑 +3 | 蔣昌鵬qtj 2026-04-01 | 3/150 |
|
|
[考研] 318求調(diào)劑,計算材料方向 +7 | 吸喵有害笙命 2026-04-01 | 8/400 |
|
|
[考研] 085801 總分275 本科新能源 求調(diào)劑 +10 | bradoner 2026-04-01 | 11/550 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +19 | 1428151015 2026-03-27 | 19/950 |
|
|
[考研] 286分調(diào)劑 +18 | Faune 2026-03-30 | 20/1000 |
|
|
[考研] 353求調(diào)劑 +4 | 拉鉤不許變 2026-04-01 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-30 | 16/800 |
|
|
[考研]
|
廈大化工 2026-04-01 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 3/150 |
|
|
[考研] 一志愿華南師范361分,化學(xué)求調(diào)劑 +4 | Nicole88888 2026-04-01 | 4/200 |
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑,不限專業(yè),有學(xué)上就行 +25 | lllllos 2026-03-30 | 26/1300 |
|
|
[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +15 | Doleres 2026-03-31 | 16/800 |
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-29 | 7/350 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +9 | 手機用戶 2026-03-26 | 9/450 |
|
|
[考研] 086000調(diào)劑 +5 | 7901117076 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +6 | zzll03 2026-03-31 | 6/300 |
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|