计算化学公社

 找回密码 Forget password
 注册 Register
Views: 127|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] pdb2gmx生成聚合物拓扑报错,提示非标准残基中一个不存在的原子未找到

[复制链接 Copy URL]

4

帖子

0

威望

101

eV
积分
105

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师好,如题,报错段落如下:
  1. Fatal error:
  2. Atom 51 in residue PIL 1 was not found in rtp entry PIL with 51 atoms while sorting atoms.
复制代码
问题在于,我构建的非标残基(PIL)仅有50个原子,提交的.pdb文件中也不存在原子名称为“51”的原子,不知各位老师能否点拨一二,感谢!

Hint:我的体系为聚离子液体(可聚合单元为甲基丙烯酸酯),单体为50个原子,将单体用甲基饱和并结构优化后提交给ligpargen,手动构建了残基在oplsaa力场(文献绝大部分使用的都是oplsaa…)下的.rtp文件;将残基PIL定义为Protein,并定义了新的atomtypes;最后将端基未饱和的聚合物.pdb文件提交pdb2gmx,但是出现了上述的报错,然而经过检查,我并没有在我定义的这些文件中找到任何名称为“51”的原子,残基中也没有第51个原子。进一步尝试了在.rtp文件中随意添加一个名为“51”的原子,随后报错信息变为:
  1. Fatal error:
  2. Atom 64 in residue PIL 1 was not found in rtp entry PIL with 52 atoms while sorting atoms.
复制代码
实在是毫无头绪了,检查了.pdb格式好像也基本正确,各位老师了解可能是哪里出了问题吗?另外,有什么更简便的产生该体系聚合物oplsaa拓扑文件的方法吗?第一次发帖,尽请各位老师批评指正,恳请各位老师解惑,不胜感激!


PIL.rtp

3.32 KB, 下载次数 Times of downloads: 1

非标残基.rtp

PolyPIL.pdb

20.61 KB, 下载次数 Times of downloads: 2

聚合物.pdb

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-24 07:30 , Processed in 0.269064 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list