| 查看: 2070 | 回复: 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这项是正值。 因此我非常困惑，是否我的计算过程是有问题的，从而计算得到了这种本身就不合理的结果？还是说，计算过程和结果可能是正确的，而我计算的这个蛋白系统就如同计算所呈现的一般？ 请见谅问题可能比较粗浅……求大神们指点！！ |
回复此楼» 猜你喜欢
水蒸气蒸馏法乳化现象
已经有3人回复
-
时间戳他又来了
已经有8人回复
-
小木虫看见有人已经知道结果了
已经有8人回复
-
J-J-W不为人知的一面
已经有11人回复
-
自然基金已经彻底沦为学术资源的交换工具
已经有36人回复
-
全人类对于AI的整体认知，尚且停留在幼儿园启蒙阶段
已经有14人回复
-
没消息就是被刷了呗
已经有10人回复
-
周末思考-婚姻长久底色
已经有12人回复
-
27考博求助
已经有4人回复
-
深圳大学材料学院“智能分子材料团队”招聘博士后
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2016-06-27 10:39:39
3楼2016-06-27 16:58:20