计算化学公社

标题: gromacs输出的pdb文件加氢数目不对怎么办 [打印本页]

作者
Author:
Yoghurt    时间: 2019-12-6 20:17
标题: gromacs输出的pdb文件加氢数目不对怎么办
用GMX跑完蛋白质的动力学,trjconv出pdb文件,发现每个残基HIS上的芳香环上的一个N原子都多加一个H(真实应该是NC双键,但是输出的默认为单键,因此自动补氢的时候N上多出一个H),请问如何能避免这种问题?谢谢

作者
Author:
puzhongji    时间: 2019-12-6 20:45
AMBER Histidine residues
Histidine (HIS in normal pdb files) is really one of three possible residues:
HID: Histidine with hydrogen on the delta nitrogen
HIE: Histidine with hydrogen on the epsilon nitrogen
HIP: Histidine with hydrogens on both nitrogens; this is positively charged.
It is up to the user to inspect the environment of each histidine and identify the type that is appropriate.
作者
Author:
sobereva    时间: 2019-12-7 20:24
trjconv过程本身不会修改结构。弄清楚导致氢不对是哪步所致。pdb文件里本身不记录形式键




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