| 查看: 3389 | 回复: 7 | ||
[求助]
gromacs中MM-GBSA计算 已有1人参与
|
|
大家好。我现在装的是gromacs4.5.4，这个版本的gromacs已经可以做GBSA了。所以我现在想用MM-GBSA的方法计算protein-ligand之间的结合自由能。 MM-GBSA方法计算结合自由能的公式如下 deltaG=EMM-TdeltaSMM+deltaGsolv 其中EMM就是protein与ligand之间的静电能和范德华能之和，我用PME方法处理长程静电力。轨迹跑完后定义energygrps=protein ligand，用rerun重新跑轨迹得到protein和ligand之间的LJ-SR,LJ-LR,COUL-SR这三项，还有系统的coul-recip项。我觉得EMM=LJ-SR+LJ-LR+COUL-SR+COUL-LR.可是coul-recip是系统总的长程静电能，而得不到protein和ligand之间的长程静电能COUL-LR。 请问大家这个问题应该怎么解决呢？ 我用显示溶剂跑的轨迹，之后应该可以用rerun来跑GBSA得到Gsolv吧？我试过可以得到但是不知道这种方法得到的结果是否正确。 |
回复此楼» 猜你喜欢
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 呼叫gromacs大神
已经有9人回复
- 请教GMX中做QM/MM
已经有1人回复
2楼2012-11-13 14:56:55
piaoluodeyu
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1242.6
- 红花: 2
- 帖子: 288
- 在线: 213.1小时
- 虫号: 1051343
- 注册: 2010-07-03
- 专业: 合成药物化学
3楼2012-11-13 15:45:50
4楼2012-11-13 20:28:31
piaoluodeyu
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1242.6
- 红花: 2
- 帖子: 288
- 在线: 213.1小时
- 虫号: 1051343
- 注册: 2010-07-03
- 专业: 合成药物化学
5楼2012-11-14 08:38:10
6楼2012-11-14 14:11:00
piaoluodeyu
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1242.6
- 红花: 2
- 帖子: 288
- 在线: 213.1小时
- 虫号: 1051343
- 注册: 2010-07-03
- 专业: 合成药物化学
7楼2012-11-14 20:26:02
JurryPeng
至尊木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 17853.5
- 散金: 33
- 红花: 25
- 帖子: 815
- 在线: 354.5小时
- 虫号: 2441687
- 注册: 2013-04-29
- 性别: MM
- 专业: 理论和计算化学
8楼2014-12-06 16:56:00
