计算化学公社

标题: 已解决：amber14处理完非标准氨基酸残基后tleap出现报错 [打印本页]

作者
Author: fuchg 时间: 2022-4-25 23:25
标题: 已解决：amber14处理完非标准氨基酸残基后tleap出现报错
本帖最后由 fuchg 于 2022-4-28 16:08 编辑

各位老师，我的体系中有一个非标准氨基酸残基SEP，现在按照教程的方法进行处理，顺利生成sep.prepin、frcmod1.sep、frcmod2.sep，但是在tleap中却出现如下错误：
(, 下载次数 Times of downloads: 22)
下面是我的sep.ac、sep.mc、sep.prepin、frcmod1.sep、frcmod2.sep、leap.in、leap.log文件，请问各位老师这种问题该如何解决呢？、
(, 下载次数 Times of downloads: 3) (, 下载次数 Times of downloads: 9) (, 下载次数 Times of downloads: 6) (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 5) (, 下载次数 Times of downloads: 3)

作者
Author: Frozen-Penguin 时间: 2022-4-26 10:29
本帖最后由 Frozen-Penguin 于 2022-4-26 10:40 编辑

.mc文件可能写错了，主链只指定CA原子，尾部应该是C而不是CB，CB是侧链，需要忽略的原子只指定NH2中的H和COOH中的O和H，如果要去掉磷酸里的H，应该在计算电荷之前去掉，否则可能会影响电荷的计算结果。请重新看一下教程中关于主链文件的说明

作者
Author: fuchg 时间: 2022-4-26 15:37

Frozen-Penguin 发表于 2022-4-26 10:29
.mc文件可能写错了，主链只指定CA原子，尾部应该是C而不是CB，CB是侧链，需要忽略的原子只指定NH2中的H和CO ...

好哒，谢谢，我在重新尝试一下

作者
Author: Frozen-Penguin 时间: 2022-4-26 18:35

fuchg 发表于 2022-4-26 15:37
好哒，谢谢，我在重新尝试一下

另外还要注意电荷量，如果去掉了磷酸上的H，电荷量应该是-2

作者
Author: fuchg 时间: 2022-4-28 15:00

Frozen-Penguin 发表于 2022-4-26 18:35
另外还要注意电荷量，如果去掉了磷酸上的H，电荷量应该是-2

好哒，谢谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)