| 查看: 3114 | 回复: 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吧？我试过可以得到但是不知道这种方法得到的结果是否正确。 |
回复此楼» 猜你喜欢
国自科送审了吗
已经有10人回复
-
博士招生
已经有5人回复
-
夜，静悄悄的
已经有8人回复
-
研究生做的很差，你们会让毕业吗？
已经有6人回复
-
26年博士申请自荐-电催化
已经有5人回复
-
2026博士或科研助理转27年博士
已经有5人回复
-
26年申博自荐-计算机视觉
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 呼叫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 (幼儿园)
- 金币: 17478.5
- 散金: 33
- 红花: 25
- 帖子: 815
- 在线: 350.8小时
- 虫号: 2441687
- 注册: 2013-04-29
- 性别: MM
- 专业: 理论和计算化学
8楼2014-12-06 16:56:00