[AMBER] Too large numbers in PDB when using CPPtraj to generate PDB

From: Myron <myron.qmandmd.gmail.com>
Date: Tue, 12 Dec 2017 06:26:16 -0600

I am using Amber to simulate a long protein link. Although the simulation
goes well, but the simulation result looks every puzzling. I tried to use
Amber to generate a series of PDB files but the generated PDB file looks
like this:

ATOM 6 HA LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   H
ATOM 7 CB LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   C
ATOM 8 HB2 LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   H
ATOM 9 HB3 LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   H
ATOM 10 CG LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   C
ATOM 11 HG2 LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   H
ATOM 12 HG3 LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   H
ATOM 13 CD LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   C
ATOM 14 HD2 LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   H
ATOM 15 HD3 LYS 1
-234821288583424.000-610838750948.00023925014113984.000 1.00 0.00
   H

So I think this is due the drifting of protein, so I try to add the* center
origin *command in CPPtraj but the problem remains.

I am very confused why this error happens.
And here is my CPPtraj scripts:

parm ../mys.top
trajin ../mys001.coord 1 50000 10
#trajin ../mys002.coord 1 50000 100
#trajin ../mys003.coord 1 50000 100
#trajin ../mys004.coord 1 50000 1000

center origin

trajout tmp.nc
trajout pdb/mdp.pdb pdb multi


It would be great if you can tell me where I use the wrong options in the
CPPtraj scripts.

Best,
Myron
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Dec 12 2017 - 04:30:05 PST
Custom Search