[gmx-users] order of atoms in pdb-file
Mr.Sridhar
sridhar at www.cdfd.org.in
Sun Jan 11 06:04:01 CET 2004
Dear gromacs users,
I faced this problem while running g_chi.
I wanted to calculate the the number of violations of ramachandran angles
of the structures in trajectory file.
The program outputs residue wise files of violations, that was fine. but
there were error messages saying
"No chi1 & chi2 angle for ALA356 "
and so on for almost all residues.
When I observed one of the pdb file of one of the snapshots(I worte it
with trjconv), I found that the order of atoms (backbone and sidechain
atoms) in residues does not match the standard pdb format. For example,
#############order in gromacs pdb file################################
ATOM 4754 N VAL 469 52.310 49.270 35.790 1.00 0.00
ATOM 4755 H VAL 469 51.440 49.270 35.300 1.00 0.00
ATOM 4756 CA VAL 469 53.460 48.600 35.170 1.00 0.00
ATOM 4757 CB VAL 469 53.170 47.120 34.910 1.00 0.00
ATOM 4758 CG1 VAL 469 54.410 46.380 34.390 1.00 0.00
ATOM 4759 CG2 VAL 469 52.690 46.340 36.140 1.00 0.00
ATOM 4760 C VAL 469 53.740 49.280 33.840 1.00 0.00
ATOM 4761 O VAL 469 52.960 49.690 32.980 1.00 0.00
###############order in a standard pdb file##############################
ATOM 219 N VAL A 60 19.146 69.476 49.828 1.00 43.93
ATOM 220 CA VAL A 60 18.468 69.172 51.091 1.00 39.04
ATOM 221 C VAL A 60 17.398 68.085 51.001 1.00 40.95
ATOM 222 O VAL A 60 17.209 67.353 51.969 1.00 36.94
ATOM 223 CB VAL A 60 17.888 70.420 51.791 1.00 39.98
ATOM 224 CG1 VAL A 60 18.437 71.732 51.199 1.00 38.84
ATOM 225 CG2 VAL A 60 16.378 70.406 51.791 1.00 39.59
Is this the reason for error message in g_chi ?
When I checked the starting pdb file I used for running dynamics it
contained the same order as that given by gromacs (above).
Does this affect the analysis of other properties or the dynamics run
itself ?
Waiting for your suggessions.
sridhar.
More information about the gromacs.org_gmx-users
mailing list