计算化学公社

标题: 求助在自行替换核苷酸中部分原子后,pdb2gmx运行报错缺少原子 [打印本页]

作者
Author:
HYUAN    时间: 2023-5-18 17:37
标题: 求助在自行替换核苷酸中部分原子后,pdb2gmx运行报错缺少原子
老师好,因为研究需要,我需要在一个DNA中的首端G核苷酸上五号碳原子上的羟基-OH替换为-S-Au,力场选择为amber的OL15力场,对应的即O5'替换为S,H5T替换为Au,替换是在高斯中替换,然后导出为.pdb文件的,之后在GROMACS中运行pdb2gmx命令时,出现以下报错:Fatal error:Residue 1 named DG of a molecule in the input file was mapped
to an entry in the topology database, but the atom O5' used in
that entry is not found in the input file. Perhaps your atom
and/or residue naming needs to be fixed.
我认为应该其中的O5'被替换所以导致缺失,这种情况应该怎么处理呢?
以下文件中1PB_1T.pdb为替换前的文件,1PB_1T_SAU.pdb为替换后的文件




作者
Author:
sobereva    时间: 2023-5-18 19:52
pdb里的残基里的原子名必须和rtp文件里相应残基对应。自己改了标准残基就相当于非标准残基,需要自己在rtp里创建新的残基信息与之对应。GROMACS基本常识性知识




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3