24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1730  |  回复: 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的回帖

彭飞龙霸王

铁虫 (初入文坛)

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

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 335求调剂 +20 想上岸呀!! 2026-04-12 23/1150 2026-04-17 10:50 by cuisz
[考研] 收到复试调剂但是去不了 +6 小蜗牛* 2026-04-16 6/300 2026-04-17 10:05 by 涵竹刘
[考研] 22408 312求调剂 +23 门路摸摸 2026-04-14 25/1250 2026-04-16 21:21 by Art1977
[考研] 0831生医工第一轮调剂失败求助 +14 小熊睿睿_s 2026-04-11 18/900 2026-04-16 20:41 by lpl364211
[考研] 307中医考研调剂 +6 于以采蘩 2026-04-14 6/300 2026-04-16 16:20 by qingfeng258
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 复试调剂 +21 积极向上; 2026-04-10 23/1150 2026-04-15 12:50 by 西北望—风沙
[考研] 调剂 +12 月@163.com 2026-04-11 12/600 2026-04-14 15:37 by zs92450
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 农学0904 312求调剂 +4 Say Never 2026-04-11 4/200 2026-04-14 09:10 by zs92450
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 调剂结束 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 283求调剂 086004考英二数二 +17 那个噜子 2026-04-10 18/900 2026-04-11 16:27 by 明月此时有
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 085600材料与化工329分求调剂 +16 叶zilin 2026-04-10 16/800 2026-04-11 11:04 by may_新宇
信息提示
请填处理意见