On Sat, November 21, 2009 at 2:36:08 AM, leila karami <karami.lei...@gmail.com>
wrote:
> I used amber03 force field in gromacs to study pr-dna interaction.
> I compare my pdb file and ffamber03.rtp file.
> I should correct my pdb file to tally with rtp file. N-terminal (NUMBER 1)
> and C-terminal
> (NUMBER 70) of my protein is GLY.
>
pdb file and rtp file relating to gly1 and gly70 are placed in follow.
> pdb file of gly1:
> 1 N NGLY 1
> 2 CA NGLY 1
> 3 C NGLY 1
> 4 O NGLY 1
> 5 1H NGLY 1
> 6 2HA NGLY 1
> 7 2H NGLY 1
> 8 3H NGLY 1
> 9 1HA NGLY 1
> rtp file of gly1:
> [ NGLY ]
> [ atoms ]
> N amber99_39 0.29430 1
> H1 amber99_17 0.16420 2
> H2 amber99_17 0.16420 3
> H3 amber99_17 0.16420 4
> CA amber99_11 -0.01000 5
>
HA1 amber99_28 0.08950 6
> HA2 amber99_28 0.08950 7
> C amber99_2 0.61630 8
> O amber99_41 -0.57220 9
I could be wrong, but I believe that the atom names in the pdb file are
supposed to match
the names in rtp entry. There may be some formatting issues as well in your
actual pdb file
since the error message references a residue GLY1 instead of residue NGLY1, but
in the
absence of the actual pdb entry line, I cannot comment on that.
Sincerely
--
Bruce D. Ray, Ph.D.
Associate Scientist
IUPUI
Physics Dept.
402 N. Blackford St.
Indianapolis, IN 46202-3273
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
