| 查看: 2112 | 回复: 2 | ||
[求助]
gromacs+GMXPBSAtool 计算结合自由能结果很奇怪啊。。。
|
|
各位老师好！ 楼主是做蛋白纳米粒子的小白半路出家接触了分子模拟，最近在做到结合自由能计算的时候遇到了问题，向大家请教。 我用的是gromacs和GMXPBSAtool，根据mmpbsa方法计算了蛋白蛋白之间的binding free energy（溶剂为水），结果如下 ΔGvwd -168.6(7) ΔGSA 0.1(0.8) ΔGelec 2698.9(53) ΔGPB -2571.0(35) ΔGnonpolar -168.5(3) ΔGpolar 127.9(5) ΔGbind -40.6(5) 对Gelec和Gpb两项，以上的结果似乎说明，我计算的这个蛋白系统中静电相互作用会非常严重地影响其稳定性，而该蛋白在溶剂（水）中的溶解又极有利于其稳定。但是我看到许多的参考文献中，Gelec这项通常是负值，而Gpb这项是正值。 因此我非常困惑，是否我的计算过程是有问题的，从而计算得到了这种本身就不合理的结果？还是说，计算过程和结果可能是正确的，而我计算的这个蛋白系统就如同计算所呈现的一般？ 请见谅问题可能比较粗浅……求大神们指点！！ |
回复此楼» 猜你喜欢
麻烦专家们看看评委们的意见（F口面上）
已经有4人回复
-
2026年叶企孙基金
已经有6人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然面上复盘~欢迎讨论
已经有13人回复
-
基金系统什么内容也没有
已经有9人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有11人回复
-
怎么查啊
已经有6人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2016-06-27 10:39:39
3楼2016-06-27 16:58:20
30