[gmx-users] problem pdb2gmx and acetylation

menica dibenedetto domenica.dibenedetto at gmail.com
Thu Sep 20 10:24:07 CEST 2012


Dear all,
I have a problem with pdb2gmx function.
I want to create the gro file of an acetylated at N-term protein.
I copy here the head of the file:
ATOM   2016  HC  ACE     0      24.770 -62.381 -11.080                       H
ATOM   2015  CT  ACE     0      25.644 -61.867 -10.739                       C
ATOM   2017  HC  ACE     0      26.283 -61.657 -11.571                       H
ATOM   2018  HC  ACE     0      26.167 -62.480 -10.034                       H
ATOM   2013   C  ACE     0      25.231 -60.546 -10.063                       C
ATOM   2014   O  ACE     0      24.058 -60.199  -9.966                       O
TER
ATOM   2019   N  MET     1      26.336 -59.728  -9.542                       N
ATOM   2020   H  MET     1      27.236 -59.762  -9.976                       H
ATOM      1  CA  MET     1      26.144 -59.704  -8.063                       C

The problem is that the output error told me that the atoms labels of
ACE not correspond to the labels in the rtp file,
but I checked it and them corresponds:
Program pdb2gmx_d, VERSION 4.5.5
Source code file: pdb2gmx.c, line: 655
Fatal error:
Atom HC in residue ACE 0 was not found in rtp entry ACE with 6 atoms
while sorting atoms.

I used the amber99sb-ildn.ff

What can I do?
I added the acetylation using gaussian

thanks for the attention
Regards
domenica



More information about the gromacs.org_gmx-users mailing list