(已解决)gmx_MMPBSA结果异常（能到正的好几百）求助

wbqdssl · 发表于 Post on 2025-7-12 23:32:58

本帖最后由 wbqdssl 于 2025-7-13 15:47 编辑

各位老师好，我在用gmx_MMPBSA计算蛋白-有机配体的MMPBSA时结果异常（能到正的好几百），但我初次接触这个程序，还不了解具体是哪方面的问题。输入文件和结果文件都在下面，恳请各位老师指点！
用的charmm-gui网站建体系，使用charmm36m力场；2开头的文件的体系是一个含LPH（LPH is a positively charged virtual particle attached to halogen atoms.）的体系，照官网的说明进行了处理，也更改了in文件里面的参数；1开头的体系是正常蛋白-有机配体体系。观察轨迹（周期性校正和平动转动都消除了）发现体系已经平衡，配体蛋白都极其稳定，配体在结合位点；输入给gmx_MMPBSA的轨迹只有蛋白和配体，没有水和离子；在动力学模拟时盐浓度是0.15M。
猜测是in文件设置有问题，但实在是没有相关经验，请各位老师指点，谢谢！

2mmpbsa_all_halo.in (8.4 KB, 下载次数 Times of downloads: 1)

2FINAL_RESULTS_MMPBSA.dat (9.61 KB, 下载次数 Times of downloads: 0)

1mmpbsa_all.in (8.4 KB, 下载次数 Times of downloads: 2)

1FINAL_RESULTS_MMPBSA.dat (9.61 KB, 下载次数 Times of downloads: 2)

体系pdb文件.zip (256.39 KB, 下载次数 Times of downloads: 0)

wbqdssl · 发表于 Post on 2025-7-13 15:49:51

问题找到了，是charmm-gui在构建配体的参数文件时出了问题，构建出来的配体不是完整的，是碎的；化学键不全。
如果其他人也遇到ΔGGAS这一项能量过大的情况可以先检查一下配体是否完整。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GROMACS] (已解决)gmx_MMPBSA结果异常（能到正的好几百）求助

评分 Rate

浏览过的版块