[gmx-users] error in pdb2gmx

Sadaf Rani sadafrani6 at gmail.com
Sun Apr 19 15:58:37 CEST 2020


Dear Gromacs users
I am facing a strange problem while doing pdb2gmx:-
When I use the command:-
 gmx pdb2gmx -f 74A-G.pdb -o 74A-G.gro
It shows the following:-
going to rename ./amber99sb-ildn.ff/aminoacids.r2b
Opening force field file ./amber99sb-ildn.ff/aminoacids.r2b
going to rename ./amber99sb-ildn.ff/dna.r2b
Opening force field file ./amber99sb-ildn.ff/dna.r2b
going to rename ./amber99sb-ildn.ff/rna.r2b
Opening force field file ./amber99sb-ildn.ff/rna.r2b
Reading 74A-G.pdb...
Read '', 4082 atoms
Analyzing pdb file
Splitting chemical chains based on TER records or chain id changing.
There are 2 chains and 0 blocks of water and 492 residues with 4082 atoms

  chain  #res #atoms
  1 'A'   489   3970
  2 'B'     3    112

All occupancies are one
Opening force field file ./amber99sb-ildn.ff/atomtypes.atp

-------------------------------------------------------
Program:     gmx pdb2gmx, version 2020-UNCHECKED
Source file: src/gromacs/gmxpreprocess/resall.cpp (line 98)

Fatal error:
Invalid atomtype format: ''

However, when I use without gmx, as below:-
pdb2gmx -f 74A-G.pdb -o 74A-G.gro
It generates topology file successfully. What's wrong with the format in
atomtypes.atp I am unable to understand. Could you please help me to find
out?
I have added a new residue to my forcefield file as mentioned in gromacs
manual.
Please correct me where I am wrong, I have added following in
atomtypes.atp file:-
;[ atomtypes ]
; name   bond_type       mass      charge    ptype      sigma      eps
      nh 14.01000     0.000       A     3.25000e-1   7.11280e-1
      hn  1.00800     0.000       A     1.06908e-1   6.56888e-2
      ca 12.01000     0.000       A     3.39967e-1   3.59824e-1
      nb 14.01000     0.000       A     3.25000e-1   7.11280e-1
      h5  1.00800     0.000       A     2.42146e-1   6.27600e-2
      nc 14.01000     0.000       A     3.25000e-1   7.11280e-1
      cd 12.01000     0.000       A     3.39967e-1   3.59824e-1
      na 14.01000     0.000       A     3.25000e-1   7.11280e-1
      c3 12.01000     0.000       A     3.39967e-1   4.57730e-1
      h2  1.00800     0.000       A     2.29317e-1   6.56888e-2
      os 16.00000     0.000       A     3.00001e-1   7.11280e-1
      h1  1.00800     0.000       A     2.47135e-1   6.56888e-2
      p5 30.97000     0.000       A     3.74177e-1   8.36800e-1
       o 16.00000     0.000       A     2.95992e-1   8.78640e-1
      oh 16.00000     0.000       A     3.06647e-1   8.80314e-1
      ho  1.00800     0.000       A     0.00000e+0   0.00000e+0
      h4  1.00800     0.000       A     2.51055e-1   6.27600e-2
       c 12.01000     0.000       A     3.39967e-1   3.59824e-1
       n 14.01000     0.000       A     3.25000e-1   7.11280e-1
      ha  1.00800     0.000       A     2.59964e-1   6.27600e-2

Thanks.

Sadaf


More information about the gromacs.org_gmx-users mailing list