计算化学公社

 找回密码 Forget password
 注册 Register
Views: 765|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 结合自由能dG绝对值过大?

[复制链接 Copy URL]

3

帖子

0

威望

37

eV
积分
40

Level 2 能力者

各位前辈好,我用gromacs分别跑了10个蛋白-蛋白复合物结构的分子动力学模拟,并使用gmx_mmpbsa.bsh计算结合自由能。
得到的dG看起来有点离谱:最大值-158kj/mol,最小值-1047kj/mol。我看别人的教程貌似△G都是 负几十 的水平,请问我的数值正常吗?
另外,查看了RMSD发现2号结构在将近2ns时RMSD异常高,不知道是什么原因?VMD跟踪了一下2号结构的轨迹文件,有些帧结构里出现一根长长的线




1695014340397.png (83.69 KB, 下载次数 Times of downloads: 4)

其中一个自由能结果文件

其中一个自由能结果文件

rmsd.png (21.67 KB, 下载次数 Times of downloads: 5)

RMSD显示2号结构接近2ns处有异常值

RMSD显示2号结构接近2ns处有异常值

1695016189718.png (115.44 KB, 下载次数 Times of downloads: 3)

VMD查看2号结构运动轨迹

VMD查看2号结构运动轨迹

59

帖子

0

威望

703

eV
积分
762

Level 4 (黑子)

2#
发表于 Post on 2023-10-1 15:44:57 | 只看该作者 Only view this author
你的轨迹文件应该是没有做修正,所以出现了那种很长的线的情况;建议消除平动与转动;进行周期性矫正

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 19:21 , Processed in 0.190435 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list