求助：如何修改.pdb文件中原子名使其和.rtp中一致

日月明

报错：
Fatal error:Atom OCT1 in residue PRO 440 was not found in rtp entry PRO with 16 atoms
while sorting atoms.
请问可以直接把PDB文件中的1OCT和2OCT改成.rtp文件中的命名即O吗？
请各位大佬赐教，救救孩子
以下是对应部分的PDB和.rtp文件截图

Seyilaxa

删掉一个端基氧然后把OCT改成和rtp里一样的名称。
有SPDBV的话把.pdb导入SPDBV再另存为一下也可以

日月明

Seyilaxa 发表于 2024-4-16 13:27
删掉一个端基氧然后把OCT改成和rtp里一样的名称。
有SPDBV的话把.pdb导入SPDBV再另存为一下也可以

感谢大佬老师指导！！请问是删掉2OCT这个氧吗？可以再请教一下为什么要删掉吗？

Seyilaxa

日月明 发表于 2024-4-16 21:22
感谢大佬老师指导！！请问是删掉2OCT这个氧吗？可以再请教一下为什么要删掉吗？

有些力场可以自动识别C末端的羧基氧原子，但有的时候可能识别不出来。没有具体信息很难判断具体是什么原因
为了省事的话直接把残基改成和.rtp文件里一样的格式，就肯定不会出错

日月明

Seyilaxa 发表于 2024-4-16 22:48
有些力场可以自动识别C末端的羧基氧原子，但有的时候可能识别不出来。没有具体信息很难判断具体是什么原 ...

好的，了解，谢谢大佬老师指导！想再进一步请教下，这样对后面的步骤会有影响吗？比如如果我要设置质子化态，这个是不是会影响电荷什么的？

Seyilaxa

日月明 发表于 2024-4-17 09:13
好的，了解，谢谢大佬老师指导！想再进一步请教下，这样对后面的步骤会有影响吗？比如如果我要设置质子化 ...

不会。使用-ignh后氢原子全部由gromacs添加，和原始结构上的氢无关。gromacs也会根据需求补全末端（默认就是-COO-）

日月明

Seyilaxa 发表于 2024-4-17 16:21
不会。使用-ignh后氢原子全部由gromacs添加，和原始结构上的氢无关。gromacs也会根据需求补全末端（默认 ...

好的，明白了！再次感谢您的指导！！

Liota

Seyilaxa 发表于 2024-4-16 13:27
删掉一个端基氧然后把OCT改成和rtp里一样的名称。
有SPDBV的话把.pdb导入SPDBV再另存为一下也可以

如果这么修改的话，末端不就只相当于是一个羰基了吗？

Seyilaxa

Liota 发表于 2024-5-14 15:33
如果这么修改的话，末端不就只相当于是一个羰基了吗？

gmx pdb2gmx会自动补全末端羧基的，其实这里主要的问题是原子名不正确，删不删末端氧倒确实问题不大，可能我当时疏忽了

Liota

Seyilaxa 发表于 2024-5-14 21:13
gmx pdb2gmx会自动补全末端羧基的，其实这里主要的问题是原子名不正确，删不删末端氧倒确实问题不大，可 ...

那请问您的意思就是删除了末端的一个O，pdb2gmx也会修复成COO-嘛？

Seyilaxa

Liota 发表于 2024-5-14 21:39
那请问您的意思就是删除了末端的一个O，pdb2gmx也会修复成COO-嘛？

是的，还可以改成别的封端基团（除了Amber力场），可以随便找个pdb试一试

Liota

Seyilaxa 发表于 2024-5-14 23:41
是的，还可以改成别的封端基团（除了Amber力场），可以随便找个pdb试一试

好的好的，谢谢您啦，解答我的困惑啦