[gmx-users] Topology for hetroatom (PRODRG)

Alok alokjain at iitk.ac.in
Tue Feb 22 14:47:31 CET 2005


hello gmx users,
               i want to simulate a peptide-linker-peptide system. The
system is XXX-NH-CH2-CH2-NH-XXX where X represents the
standard amino acids.
 To simulate this system i want to first generate a topology file for the
linker section of the system.I did so by using  PRODRG server by suppling
the linker section as the input..
                  The input coordinate file to the server and the output
from the server are as follows:

Input:
HETATM   14  N     2     2      -0.564   0.106   1.392
HETATM   15  C     2     2       0.779   0.380   0.903
HETATM   16  C     2     2       1.220   1.824   1.201
HETATM   17  N     2     2       2.587   2.069   0.782

Output:

[ moleculetype ]
; Name nrexcl
 _2      3

[ atoms ]
;   nr      type  resnr resid  atom  cgnr   charge     mass
     1        NL     1   _2     NAA     1    0.704  14.0067
     2         H     1   _2     HAA     1    0.059   1.0080
     3         H     1   _2     HAB     1    0.059   1.0080
     4         H     1   _2     HAC     1    0.059   1.0080
     5       CH2     1   _2     CAB     1    0.119  14.0270
     6       CH2     1   _2     CAC     2    0.119  14.0270
     7        NL     1   _2     NAD     2    0.704  14.0067
     8         H     1   _2     HAE     2    0.059   1.0080
     9         H     1   _2     HAF     2    0.059   1.0080
    10         H     1   _2     HAD     2    0.059   1.0080

The queries to the output are :

1)   The linker now has three hydrogens attached to nitrogen whereas the 
linker with the peptide will have only one hydrogen attached to
nitrogen.So can we restrict the number of hydrogens to one or should i now
delete two hydrogens from the topology file before using it for my
system????

2)   Should I use the same parameters including the charges directly as
given in the output in the topology file for my simulation or are there
other issues to look into before using this topology file?????

3)  These parameters were developed by the server for an isolated linker
(N-C-C-N).So can i use the same parameters for the linker when it is bound
to my peptide (XXX-Linker-XXX)????

Please shed some light into this issue.

Thanking all in advance.
                        Alok Jain






More information about the gromacs.org_gmx-users mailing list