计算化学公社

标题: pdb2gmx 识别氨基酸名称错误问题 [打印本页]

作者
Author: cylqcjx 时间: 2022-5-16 13:08
标题: pdb2gmx 识别氨基酸名称错误问题
请教各位老师

运行pdb2gmx命令时报错（如图）

(, 下载次数 Times of downloads: 65)

但是检查输入的.pdb文件第76号是LYS而不是GLY（如下图）。于是在这里卡住了不知道该怎么修正。。

(, 下载次数 Times of downloads: 18)

这个.pdb通过pymol引入了几个突变位点，运行pdb2gmx的时候报过原子类型不能被amber99力场识别的问题，通过参照可被识别氨基酸的原子类型信息成功修正了。

但是76号没有突变过，核对了原始的.pdb文件（1jpz)是一样的。

于是就在这里卡住了。。想请教各位这个问题如何解决。

尝试用 -ignh命令继续，但是在最后开始长时间mdrun的时候频繁报blow up类型错误，于是才倒回来重新梳理。

谢谢大家！

作者
Author: Acee 时间: 2022-5-16 15:09
本帖最后由 Acee 于 2022-5-20 13:40 编辑

用pdbfixer修复一下试试

作者
Author: cylqcjx 时间: 2022-5-16 15:52

Acee 发表于 2022-5-16 15:09
把有C和O的那两行移动到NZ行的后面

请问 Atom序号需要同时调整吗？刚才试了一下还是报相同的错误

作者
Author: cylqcjx 时间: 2022-5-16 16:38
我按照力场里.rtp文件的类型，逐个修正了报错的残基信息，但是pdb2gmx会错误的把23号 ASP识别成HIS，一直报HIS incomplete ring.即使我把整个ASP23删掉了，还是在报错[attachimg]49268
[/attachimg]

(, 下载次数 Times of downloads: 13)

请教各位这个情况是bug吗

作者
Author: cylqcjx 时间: 2022-5-16 16:40
(, 下载次数 Times of downloads: 32)
附完整的.pdb文件

作者
Author: sobereva 时间: 2022-5-16 17:19
pdb里残基序号如果不是从1开始连续排的，pdb2gmx报错里的残基号就不是pdb里直接记录的残基号，搞清楚对应关系
残基里原子顺序随意，只看原子名
确保pdb里的残基中没缺重原子

作者
Author: cylqcjx 时间: 2022-5-20 12:51

sobereva 发表于 2022-5-16 17:19
pdb里残基序号如果不是从1开始连续排的，pdb2gmx报错里的残基号就不是pdb里直接记录的残基号，搞清楚对应关 ...

谢谢sob老师

仔细检查之后发现是一些格式问题，修正之后已经正常运行了

作者
Author: MOI_001 时间: 2023-12-20 10:26

cylqcjx 发表于 2022-5-20 12:51
谢谢sob老师

仔细检查之后发现是一些格式问题，修正之后已经正常运行了

请问你是用什么方法修正成功的呀？我也报了这个错，现在还没解决

作者
Author: ydp111~ 时间: 2024-9-30 18:49

MOI_001 发表于 2023-12-20 10:26
请问你是用什么方法修正成功的呀？我也报了这个错，现在还没解决

请问您修正成功了吗？我现在也出现了这个报错

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