24小时热门版块排行榜

返回列表

sunnylfq

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 39
帖子: 3
在线: 4.4小时
虫号: 1961693
注册: 2012-08-29
专业: 化工系统工程

[求助] 关于GROMACS模拟蛋白质的重复性

对象是酶（不到700个氨基酸）与底物小分子的结合体
手动在二者之间添加共价键，通过MD观察小分子与活性位点间是否形成预测的氢键。

问题是：MD跑了2次，各20ns，一次出现氢键，一次没有。是GROMACS-MD的重复性不好么？

看文献里用AMBER跑蛋白质，模拟结构和晶体结构的RMSD才不到1A，但是GROMACS跑的话一般都是2~3A，是GROMACS自己的问题么？

回复此楼

» 收录本帖的淘帖专辑推荐

MDs-Gromacs

» 猜你喜欢

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

1楼 2013-07-17 16:24:49

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

dubo

金虫 (著名写手)

模拟EPI: 1
应助: 23 (小学生)
贵宾: 0.779
金币: 569.2
散金: 3220
红花: 31
沙发: 1
帖子: 1821
在线: 349.5小时
虫号: 559371
注册: 2008-05-17
性别: GG
专业: 高分子组装与超分子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
chaizhm: 金币+2, 谢谢~ 2013-07-18 16:16:58
sunnylfq: 金币+5, ★★★很有帮助 2013-07-19 08:59:39

你是用的sd跑的？每次的初始态都一样么？（初始速度是随机生成的？）

如果你用的积分算法是md,而且每次运行时初始速度都相同（都设置为0，或者输入一个初始速度），那么gromacs-md抛出的结果严格一样。

如果是做平行样本的话，建议多做几组，可以适当的延长模拟时间

赞一下

回复此楼

2楼2013-07-17 19:44:02

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

sunnylfq

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 39
帖子: 3
在线: 4.4小时
虫号: 1961693
注册: 2012-08-29
专业: 化工系统工程

引用回帖:

2楼: Originally posted by dubo at 2013-07-17 19:44:02
你是用的sd跑的？每次的初始态都一样么？（初始速度是随机生成的？）

如果你用的积分算法是md,而且每次运行时初始速度都相同（都设置为0，或者输入一个初始速度），那么gromacs-md抛出的结果严格一样。

如果是 ...

我仔细看了下mdp文件。
整个模拟过程包含EM，6次NVT模拟（从50到300K升温），一次NPT模拟，一次MD。

所有过程用的积分器都是md。但是NVT采用麦克斯韦分布随机生成初速度，NPT和MD用的是gen_vel=no.

这样的话，其实整个过程还是受到随机初速度的影响吧？
如果要得到精确重复性结果，NVT/NPT/MD都要设置gen_vel=no是么？

赞一下

回复此楼

3楼2013-07-18 11:35:13

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

相关版块跳转我要订阅楼主 sunnylfq 的主题更新

返回列表