24小时热门版块排行榜

>论坛更新日志 (3940)
>文献求助 (426)
>导师招生 (365)
>虫友互识 (363)
>休闲灌水 (239)
>考博 (166)
>硕博家园 (120)
>招聘信息布告栏 (116)
>博后之家 (95)
>论文投稿 (89)
>基金申请 (67)
>考研 (61)
>公派出国 (60)
>论文道贺祈福 (53)
>找工作 (42)
>教师之家 (36)

返回列表

飞飞是胖纸

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 61.5
红花: 1
帖子: 16
在线: 8.5小时
虫号: 2468481
注册: 2013-05-17
专业: 生物化学

[求助] 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这项是正值。
因此我非常困惑，是否我的计算过程是有问题的，从而计算得到了这种本身就不合理的结果？还是说，计算过程和结果可能是正确的，而我计算的这个蛋白系统就如同计算所呈现的一般？
请见谅问题可能比较粗浅……求大神们指点！！

回复此楼