计算化学公社

标题: (已解决)gmx_MMPBSA结果异常(能到正的好几百)求助 [打印本页]

作者
Author:
wbqdssl    时间: 2025-7-12 23:32
标题: (已解决)gmx_MMPBSA结果异常(能到正的好几百)求助
本帖最后由 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文件设置有问题,但实在是没有相关经验,请各位老师指点,谢谢!
(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 0)


作者
Author:
wbqdssl    时间: 2025-7-13 15:49
问题找到了,是charmm-gui在构建配体的参数文件时出了问题,构建出来的配体不是完整的,是碎的;化学键不全。
如果其他人也遇到ΔGGAS这一项能量过大的情况可以先检查一下配体是否完整。




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3