[gmx-users] pdb2gmx + non-protein

Alberto Imparato Alberto.Imparato at na.infn.it
Thu Jan 30 20:27:18 CET 2003


Hi All
I'm trying to obtain a topology file starting from a very simple .gro
file, using pdb2gmx .
My 2_dppc.gro file contains 2 DPPC molecules, each containing 50 atoms:

DPPC
  100
    1DPPC  C1  1   1.305  0.021  4.907   -0.4017   -0.1211   -0.18
    1DPPC  C2  2   1.345  0.856  4.733   -0.171   -0.6351   -0.2036
    1DPPC  C3  3   1.191  0.812  4.904   0.1336   0.3719   0.2801
    1DPPC  N4  4   1.244  0.919  4.819   -0.0492   -0.1049   -0.4439
    1DPPC  C5  5   1.131  0.985  4.751   -0.7172   -0.308   0.4587
    ....
    ....
    2DPPC  C1  51   1.305  0.021  0.093   -0.4017   -0.1211   -0.18
    2DPPC  C2  52   1.345  0.856  0.267   -0.171   -0.6351   -0.2036
    2DPPC  C3  53   1.191  0.812  0.096   0.1336   0.3719   0.2801
    2DPPC  N4  54   1.244  0.919  0.181   -0.0492   -0.1049   -0.4439
    2DPPC  C5  55   1.131  0.985  0.249   -0.7172   -0.308   0.4587
    ....
    ....

then I type "pdb2gmx -f 2_dppc.gro -p topol -o out.gro"

the topology file topol.top I obtain, looks like

#######################
#include "ffgmx.itp"

[ moleculetype ]
; Name            nrexcl
Protein             3

[ atoms ]
;   nr       type  resnr residue  atom   cgnr     charge       mass  typeB
chargeB      massB
     1        CH3      1   DPPC     C1      1      0.124     15.035   ;
qtot 0.124
     2        CH3      1   DPPC     C2      1      0.124     15.035   ;
qtot 0.2
....
....
[ system ]
; Name
DPPC

[ molecules ]
; Compound        #mols
Protein             1
########################
WHY the 2 molecules are not recognized as DPPC molecules, but are treated
as 1 single Protein ???
The same happens when I first transform 2_dppc.gro to  2_dppc.pdb using
editconf and the process the latter file with pdb2gmx.
Thak you
Alberto








More information about the gromacs.org_gmx-users mailing list