24小时热门版块排行榜    

查看: 1907  |  回复: 6

gytjyb

新虫 (小有名气)

[交流] 【求助】写了个MOE的动力学的过程,希望大家批评几句! 已有3人参与

因为我是初学分子模拟,正好我我们学校有MOE的版权。
       所以我试着用MOE为我的同源模型,跑了一个分子动力学的优化。请各位高手看看一下过程,我也不知道该怎么用moe,我想这是通常一般情况下的过程,恳请各位师兄姐,评价一下我有哪些做得不对。可能我有很多不对的地方,也有许多不懂得地方,请提出来。谢谢您了!
      我是想构建了120个氨基酸的蛋白模型,想对它做动力学优化。我不知道对不对,有好多我都不太明白。
      首先,我给蛋白用AMBER99力场加氢,修正电荷。然后又给蛋白加了一圆形层的水,最大原子到边界的厚度为7,是在ph7的时候加的。(我不知道怎么加离子,也还不会判断是否电中性)
       然后我用500的力量固定蛋白,先优化外面的水到系统为RMS=1,然后整个放开系统优化到能量最小化,为RMS=0.05。
       然后我再开始跑动力学,因为我想到动力学先要吧外面的水跑好吧。所以,我先从50K升温度到310K,用了50ps,选用NVT algorithrm用NPA(什么意思)。这时候在动力学选项里我选择constrait all bond(这个选项默认的意思除水以外都限制,有lights bond,all bond,none三种情况,我不知道这么做对不对,为什么这样做),有关水的那项,我选择用rigid(不知道什么意思),然后用50PS跑跑升温的过程。
       升温完成以后,我继续上次的模拟,我用1个NS,这时我想应该是否选择constrait none 什么都不限制了(为什么,这样对吗),维持300K, 其它都保持不变,一直到整个过程完,取最后200PS的平均构像。 ( 然后我发现一个问题,不知道怎样去看整个过程RMSD的变化。)   
       我就想用MOE做个模型优化,我不知道可否这样做!MOE的内容很多,但是没有动力学的列子!
       谢谢大家,祝大家好运!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dragontj

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
请问你是怎么加水和平衡电荷的。
2楼2010-07-02 15:01:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gytjyb

新虫 (小有名气)

solvate 里面似乎可以加水,至于平衡电荷,我没有看到,倒是参数设置里面可以指定PH
3楼2010-08-31 14:36:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

彭飞龙霸王

铁虫 (初入文坛)

ghcacj:请不要复制别人的回帖和LZ发帖内容。 2010-08-31 15:49:59
为我是初学分子模拟,正好我我们学校有MOE的版权。
       所以我试着用MOE为我的同源模型,跑了一个分子动力学的优化。请各位高手看看一下过程,我也不知道该怎么用moe,我想这是通常一般情况下的过程,恳请各位师兄姐,评价一下我有哪些做得不对。可能我有很多不对的地方,也有许多不懂得地方,请提出来。谢谢您了!
      我是想构建了120个氨基酸的蛋白模型,想对它做动力学优化。我不知道对不对,有好多我都不太明白。
      首先,我给蛋白用AMBER99力场加氢,修正电荷。然后又给蛋白加了一圆形层的水,最大原子到边界的厚度为7,是在ph7的时候加的。(我不知道怎么加离子,也还不会判断是否电中性)
       然后我用500的力量固定蛋白,先优化外面的水到系统为RMS=1,然后整个放开系统优化到能量最小化,为RMS=0.05。
       然后我再开始跑动力学,因为我想到动力学先要吧外面的水跑好吧。所以,我先从50K升温度到310K,用了50ps,选用NVT algorithrm用NPA(什么意思)。这时候在动力学选项里我选择constrait all bond(这个选项默认的意思除水以外都限制,有lights bond,all bond,none三种情况,我不知道这么做对不对,为什么这样做),有关水的那项,我选择用rigid(不知道什么意思),然后用50PS跑跑升温的过程。
       升温完成以后,我继续上次的模拟,我用1个NS,这时我想应该是否选择constrait none 什么都不限制了(为什么,这样对吗),维持300K, 其它都保持不变,一直到整个过程完,取最后200PS的平均构像。 ( 然后我发现一个问题,不知道怎样去看整个过程RMSD的变化。)   
       我就想用MOE做个模型优化,我不知道可否这样做!MOE的内容很多,但是没有动力学的列子!
       谢谢大家,祝大家好运!
本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=1880455&fpage=1
4楼2010-08-31 15:14:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

彭飞龙霸王

铁虫 (初入文坛)

请问你是怎么加水和平衡电荷。
5楼2010-08-31 15:15:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

彭飞龙霸王

铁虫 (初入文坛)

你做了试验了吗????
6楼2010-08-31 15:17:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gytjyb

新虫 (小有名气)

★ ★
lei0736(金币+2):欢迎发起讨论 2010-09-01 22:26:47
edit---build--water soak  构建溶液。
protonate 可以设定PH,T,盐等信息。
我就知道这些,还可能不对,你们可以参考下!
7楼2010-08-31 16:43:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gytjyb 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 7s8du2bt8y 2026-06-26 7/350 2026-06-27 18:48 by ztgu5ulw9z
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 7s8du2bt8y 2026-06-26 8/400 2026-06-27 18:47 by ztgu5ulw9z
[考博] 27年博士招生信息 +5 rvnc 2026-06-26 8/400 2026-06-27 18:32 by rvnc
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 7s8du2bt8y 2026-06-26 5/250 2026-06-27 17:07 by ztgu5ulw9z
[基金申请] 无聊看看filecode +7 流流伤 2026-06-23 10/500 2026-06-27 16:19 by 北京中袖冉老师
[基金申请] 无聊看看时间戳打发时间 +8 虎鹤双形 2026-06-23 8/400 2026-06-27 15:53 by miaochunhui
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 7s8du2bt8y 2026-06-26 4/200 2026-06-27 12:47 by 9g0rmhtq5w
[论文投稿] 求推荐期刊,重谢 +4 girlbaby 2026-06-23 4/200 2026-06-26 16:52 by 不打工牛马
[有机交流] 求助!! 5+3 我啥都没看见 2026-06-24 4/200 2026-06-26 09:35 by 951037019
[微米和纳米] kh550接枝Sio2失败求助 70+3 娃儿方便面 2026-06-20 3/150 2026-06-26 09:23 by 猫鱼不是鱼
[论文投稿] 职称论文投稿 170+4 guoj5292 2026-06-22 12/600 2026-06-25 22:27 by Yanyanoo
[有机交流] 反应求助 10+3 slz_1986 2026-06-24 6/300 2026-06-25 21:38 by nBu锂
[论文投稿] 基于自然哲学类比的风化壳型稀土矿 +4 太一新韵 2026-06-22 16/800 2026-06-25 18:53 by 太一新韵
[文学芳草园] 看《给阿ma的情书》有感 +6 myrtle 2026-06-21 10/500 2026-06-25 17:54 by myrtle
[基金申请] 2026年WR青拔进展 +5 chs564851482 2026-06-24 7/350 2026-06-24 18:15 by chs564851482
[基金申请] 中!中!中! +10 zhse276 2026-06-22 10/500 2026-06-24 16:48 by zjhzf5201018
[基金申请] 会评什么时候开始? +3 Vivilian 2026-06-24 4/200 2026-06-24 16:30 by jurkat.1640
[基金申请] 国自然申请五篇代表作大比拼,感觉这个是最重要的 +7 naalan7001 2026-06-22 12/600 2026-06-24 14:02 by naalan7001
[基金申请] 咨询 +3 _xyan818 2026-06-24 3/150 2026-06-24 08:12 by Equinoxhua
[基金申请] 评委有多少概率知道其他专家手中有哪些人的本子? +6 huitong441 2026-06-22 6/300 2026-06-23 15:45 by 新城子曾
信息提示
请填处理意见