[gmx-users] missing of methyl at N-terminal
Mark Abraham
Mark.Abraham at anu.edu.au
Fri Nov 7 14:03:52 CET 2008
Bhawana Gupta wrote:
> thankyou for ur reply.
> exact error which i m getting, when running pdb2gmx command with the pdb
> which i had taken from insightII is
> Program pdb2gmx, VERSION 4.0
> Source code file: pdb2gmx.c, line: 429
>
> Fatal error:
> Atom C in residue NH2 6 not found in rtp entry with 3 atoms
> while sorting atoms
So from this you can work out that pdb2gmx is guessing that residue 6
(named "NHM" in your pdb file) is intended to be an NH2 residue (which
has 3 atoms), but then gets confused trying to match one of them to the
C atom in the residue named "NHM". You need to name residues according
to the names in the .rtp database. Here, "NHM" presumably isn't right
for N-methylamine your force field. Go and look up what it should be and
edit your input to help pdb2gmx along.
Computers can't be programmed to guess correctly in all cases what the
user meant. You need to take some responsibility for telling it useful
things :-) A good working knowledge of the contents of chapter 5 of the
manual is invaluable here.
Mark
More information about the gromacs.org_gmx-users
mailing list