Loading0760 发表于 2024-9-29 16:26
跟src/gromacs/gmxpreprocess/pdb2gmx.cpp (line 376)没关系,读取的时候就是/opt/gromacs-2023.4/share/gro ...

用你的方法成功了，谢谢！！

Loading0760 发表于 2024-10-1 15:09
是B链段的ARG,  这个蛋白的结构你是怎么得到的, 我检查了一下, B链段有氨基酸缺失, 如图所示, 氨基酸直接从 ...

琢磨了几天，终于成功生成拓扑了，非常感谢大神

是B链段的ARG,  这个蛋白的结构你是怎么得到的, 我检查了一下, B链段有氨基酸缺失, 如图所示, 氨基酸直接从74号跳到了78号, 与此同时, ARG也少了很多原子.
这个结构应该是有PDB ID的吧,可以用同源建模把缺失的补齐.
我是直接把B链删除了,然后测试了一下,能跑. 后续三个报错都是关于力场参数的, 这是amber14的老问题的.

hdb文件里,第三行加氢的文本处
1        1        H1        C        O        CA
修改成H1

Loading0760 发表于 2024-9-30 09:16
把PCA.hdb的内容加到aminoacid.hdb里,不要用专门的PCA.hdb文件.

我加进去了，还是一样的错误呢

ydp111~ 发表于 2024-9-29 20:48
谢谢大神，按照您的建议在aminoacids.r2b文件加上了，之前的报错没有啦，但是又出现了新的报错

把PCA.hdb的内容加到aminoacid.hdb里,不要用专门的PCA.hdb文件.

跟src/gromacs/gmxpreprocess/pdb2gmx.cpp (line 376)没关系,读取的时候就是/opt/gromacs-2023.4/share/gromacs/top/amber99sb-ildn这读取的.
我建议你不要动gromacs-2023.4/share/gromacs/top/下的拓扑文件, 可以自己拷贝到工作目录下,然后在工作的目录修改,简言之,就是别动系统给的默认力场.

'PCA' as a starting terminus
尝试一下这样,在aminoacids.r2b的文件结尾加一行这个.
PS: 我自己的非标准残基定义的名称是PCB,有点巧

ydp111~ 发表于 2024-9-28 09:01
老师，我是加入Amber99sb-ildn力场里，用的也是这个力场

看清楚pdb2gmx运行时候的提示，读的是哪个目录的
安装了不止一个gromacs的时候老有人弄错

8成是没加对地方，诸如加到了A力场的rtp里，但选择的是B力场

Loading0760 发表于 2024-9-26 17:59
你把这个加进Amber99sb-ildn的力场文件了吗

我之前用的是amber14sb_parmbsc1力场,当时是把对应的文件加 ...

加进去啦