[gmx-users] TWO AMINOACIDS WITH SAME NUMBER IN .GRO FILE

Michetti Davide dmichetti at gmail.com
Thu Oct 16 14:13:13 CEST 2014


Dear gmx users,
I have just realized, at the end of a molecular dynamics simulation, that I had in my .gro file to aminoacids with the same number. This happened cause in my pdb file, what was supposed to be residue 53, was instead 52A (due to an insertion of Lys in that position, compared to other homologs). Those are the lines in the .pdb I am talking about.

ATOM    482  HA  ASN A  52      11.594  -3.380  37.114  1.00  0.00           H
ATOM    483  HB3 ASN A  52      10.188  -5.588  35.526  1.00  0.00           H
ATOM    484  HB2 ASN A  52      11.571  -5.800  36.555  1.00  0.00           H
ATOM    485 HD21 ASN A  52       7.924  -5.965  38.224  1.00  0.00           H
ATOM    486 HD22 ASN A  52       8.207  -5.978  36.494  1.00  0.00           H
ATOM    487  N   LYS A  52A     11.898  -2.478  34.227  1.00 26.03           N
ATOM    488  CA  LYS A  52A     12.656  -2.134  33.020  1.00 26.03           C
ATOM    489  C   LYS A  52A     12.428  -3.160  31.890  1.00 25.58           C
ATOM    490  O   LYS A  52A     13.333  -3.399  31.090  1.00 25.54           O
ATOM    491  CB  LYS A  52A     14.147  -1.837  33.318  1.00 26.46           C

So basically in the .gro file the residue 52 is composed by the ASN 52 and LYS 52. The question is if this can constitute a problem for the simulation? It seems that everything went smooth and I got no error, but I would like to double check..

Thanks for your time,

Davide


More information about the gromacs.org_gmx-users mailing list