[gmx-users] error in running pdb2gmx
hasthi
durgs7krish at gmail.com
Mon Nov 18 12:11:07 CET 2013
Dear Justin,
Thank you for your reply. Yeah I have modified my force
field to include the phosphorylated Serine residue.
Below is the pdb2gmx command employed
pdb2gmx -f alpha_ser129p.pdb -o alpha_sep129p.gro -ff gromos43a1 -ignh
error :
Fatal error:
Atom OXT in residue ALA 140 was not found in rtp entry ALA with 6 atoms
while sorting atoms.
snippet of my PDB file.
ATOM 2011 N ALA A 140 205.034 -53.935 -18.851 1.00
16.70 N
ATOM 2012 CA ALA A 140 204.971 -52.762 -19.777 1.00
17.17 C
ATOM 2013 C ALA A 140 203.730 -51.927 -19.454 1.00
17.37 C
ATOM 2014 O ALA A 140 202.844 -52.453 -18.799 1.00
17.40 O
ATOM 2015 CB ALA A 140 204.885 -53.253 -21.221 1.00
17.38 C
ATOM 2016 OXT ALA A 140 203.683 -50.781 -19.870 1.00
17.61 O1-
ATOM 2017 H01 ALA A 140 204.210 -54.107 -21.275 1.00
17.38 H
ATOM 2018 H02 ALA A 140 204.184 -54.311 -18.455 1.00
16.70 H
ATOM 2019 H03 ALA A 140 205.867 -52.154 -19.652 1.00
17.17 H
ATOM 2020 H04 ALA A 140 205.876 -53.550 -21.564 1.00
17.38 H
ATOM 2021 H05 ALA A 140 204.507 -52.451 -21.855 1.00
17.38 H
Besides this I also find the following lines in my output:
Identified residue MET1 as a starting terminus.
Warning: Residue SEP129 in chain has different type (Other) from starting
residue MET1 (Protein).
Warning: Residue GLU130 in chain has different type (Protein) from starting
residue MET1 (Protein).
Warning: Residue GLU131 in chain has different type (Protein) from starting
residue MET1 (Protein).
Warning: Residue GLY132 in chain has different type (Protein) from starting
residue MET1 (Protein).
Warning: Residue TYR133 in chain has different type (Protein) from starting
residue MET1 (Protein).
More than 5 unidentified residues at end of chain - disabling further
warnings.
Identified residue PRO128 as a ending terminus.
8 out of 8 lines of specbond.dat converted successfully
Special Atom Distance matrix:
MET1 MET5 HIS50 MET116
SD7 SD41 NE2360 SD817
MET5 SD41 1.034
HIS50 NE2360 5.849 5.344
MET116 SD817 7.199 7.369 10.906
MET127 SD905 8.542 8.838 12.786 2.065
Start terminus MET-1: NH3+
End terminus PRO-128: COO-
Why is the end termini is at Pro 128? Does it mean it does not consider the
rest of the residues?? Please explain me regarding this.
Regards,
Hasthi
More information about the gromacs.org_gmx-users
mailing list