计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › Gromacs中CHARMM力场无法正确识别钙离子

Gromacs中CHARMM力场无法正确识别钙离子

查看数: 279 | 评论数: 5 | 收藏 Add to favorites 0

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

帖子模式

fastsleep

发布时间: 2026-6-10 15:30

正文摘要:

本帖最后由 fastsleep 于 2026-6-10 15:30 编辑 Gromacs里使用CHARMM36力场生成拓扑文件时报错，Amber力场无这个问题，看报错应该是pdb文件中的钙离子无法被CHARMM力场正确识别，但是我的蛋白含有钙离子结合位点 ...

回复 Reply

dayu8278 发表于 Post on 2026-6-13 10:50:23

fastsleep 发表于 2026-6-10 21:31
好的，之前修改之后，在pymol里打开位置变了，不知道有没有影响

如果CA位置变了，显然就是错了

fastsleep 发表于 Post on 2026-6-10 21:31:39

student0618 发表于 2026-6-10 17:14
检查ions.itp的原子名，修改对应的钙离子原子名/残基名

也可以生成蛋白拓朴后手动处理拓朴include对应离 ...

好的，之前修改之后，在pymol里打开位置变了，不知道有没有影响

fastsleep 发表于 Post on 2026-6-10 21:29:40

sobereva 发表于 2026-6-10 17:13
“识别成另一条链” 不意味着有任何错误，只要拓扑信息正确就完事，要搞清楚内在逻辑。北京科音分子动力学 ...

好的，谢谢老师，我去试试看

student0618 发表于 Post on 2026-6-10 17:14:22

检查ions.itp的原子名，修改对应的钙离子原子名/残基名

也可以生成蛋白拓朴后手动处理拓朴include对应离子

sobereva 发表于 Post on 2026-6-10 17:13:33

“识别成另一条链” 不意味着有任何错误，只要拓扑信息正确就完事，要搞清楚内在逻辑。北京科音分子动力学与GROMACS培训班（http://www.keinsci.com/KGMX）里讲钙调素模拟的例子里，钙调素结合了四个钙，把pdb里钙离子的残基名改成与要用的力场的目录下的rtp里钙离子对应的名字后，四个钙都被pdb2gmx当做“一条链”单独产生了一个itp文件，模拟完全正常。
若不想被这么处理，就把pdb里钙离子先去掉后用pdb2gmx产生拓扑文件，然后自己把钙离子加入gro里，主拓扑文件里自己include ions.itp

某某GUI那货根本没必要试，完全多余

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