Dear Amber users
I have been using PTRAJ to extract dihedral angles for a system using amber
generated pdb files and the corresponding prmtop file.
However, on comparing the ptraj calculated dihedral angles to those obtained
from VMD I noticed substantial differences. Indeed, it appears that the
ptraj derived dihedrals are wildly wrong.
The input file I used for ptraj was named ptraj.in:
*************************************************************
trajin ptraj.coords #file ptraj.coords is an "ambpdb" generated pdb file
#residue 11
dihedral phi_11 :11.C :12.N :12.CA :12.C out __phi_11
dihedral psi_11 :11.N :11.CA :11.C :12.N out __psi_11
#residue 12
dihedral phi_12 :12.C :13.N :13.CA :13.C out __phi_12
dihedral psi_12 :12.N :12.CA :12.C :13.N out __psi_12
#residue 13
dihedral phi_13 :13.C :14.N :14.CA :14.C out __phi_13
dihedral psi_13 :13.N :13.CA :13.C :14.N out __psi_13
**************************************************************
The above input was generated using a format akin to that specified in the
amber “PTRAJ tutorial” site:
http://amber.scripps.edu/tutorial/ptraj/index.html
To run the script I used: --> ptraj prm.top ptraj.in
The above procedure produced dihedrals that are wrong. Considering just the
phi angles I got:
Residue 11: 125.0
Residue 12:-168.6
Residue 13: 89.9
for the input lines shown in "ptraj.in" above.
I know they are wrong because I tested a structure that was already known to
have beta sheets, and the PTRAJ calculated torsion angles fell way outside
the range of allowed Ramachandran values. The values according to vmd are:
Residue 11: -129.8
Residue 12: -114.9
Residue 13: -93.7
which I presume are right since they fall into allowed Ramachandran phi
regions.
Has anyone else suffered problems of this nature, or is there something that
I've missed?
Many thanks in advance
hayden
______________________________________________
Hayden Eastwood
Perdita Barran Research Group
Joseph Black Building
Edinburgh University
West Mains Road
EH9 3JJ
Tel: 0131 650 4773
e-mail: s0237717.sms.ed.ac.uk
Research page:
http://homepages.ed.ac.uk/pbarran/PBRG/
"You have to be an academic to believe some things; no ordinary person would
be so stupid." -George Orwell
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jan 11 2006 - 06:10:05 PST
