mdrun输出.gro文件ligand散开

查看数: 1033 | 评论数: 7 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2023-3-19 09:28

尝试做一个配体蛋白水溶液模型，但是在能量最小化后发现ligand散开（如下图）。，.trr文件因为体积比较大不便上传，请谅解，整个过程用3080跑大概15s，并不会您耽搁太久 ...

回复 Reply

MGET127 发表于 Post on 2023-3-21 08:05:07

sobereva 发表于 2023-3-21 02:55
太缺乏力场常识
GROMOS和AMERR/GAFF完全不兼容，要include也是include AMBER的forcefield.itp

好的老师，我之后会注意这个问题。

sobereva 发表于 Post on 2023-3-21 02:55:57

MGET127 发表于 2023-3-20 22:51
目前得出的结论：我引入的gmx的"gromos54a7.ff/forcefield.itp"中的[default]和sobtop得出的.top中的[defau ...

太缺乏力场常识
GROMOS和AMERR/GAFF完全不兼容，要include也是include AMBER的forcefield.itp

MGET127 发表于 Post on 2023-3-20 22:51:46

目前得出的结论：我引入的gmx的"gromos54a7.ff/forcefield.itp"中的[default]和sobtop得出的.top中的[default]值不同，我之前选用的为系统自带参数导致.top文件错误。

MGET127 发表于 Post on 2023-3-20 08:46:57

sobereva 发表于 2023-3-20 04:56
从你的第三张图看，结构已经严重扭曲了，诸如甲基碳都没有sp3杂化应有的结构。不管散不散开，当前结果都已 ...

谢谢老师回答，我是通过sobtop，124顺序得到的.top文件。我再尝试下其他方法。

sobereva 发表于 Post on 2023-3-20 04:56:50

从你的第三张图看，结构已经严重扭曲了，诸如甲基碳都没有sp3杂化应有的结构。不管散不散开，当前结果都已经没法用了，分子的拓扑文件一定存在硬伤。

MGET127 发表于 Post on 2023-3-19 12:54:21

牧生发表于 2023-3-19 12:24
大概率是.top文件里面，各个物质出现的顺序不对。

谢谢老师回答，但是我检查了我.top文件中的[molecules]栏，顺序是蛋白-配体-SOL-NA，和.gro文件中的顺序是一样的，数目也能对上。不太明白是什么顺序不对。

牧生 发表于 Post on 2023-3-19 12:24:16

大概率是.top文件里面，各个物质出现的顺序不对。

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