> 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.
As you might have imagined, we tested the dihedral routine prior to
releasing it which suggests that something else is funny, either with
respect to the coordinates, atom selection or something else.
What version of ptraj are you running?
Just to check, I created a small blocked peptide (ACE-ALA-ALA-NME), ran
short GB dynamics and then calculated phi/psi and matched the results to
those calculated with Midas-Plus. The results agreed. So, the question
is why are you getting funny values?
(1) Note that (unless you are offset by one intentionally) unless you have
a blocking residue (like an acetyl or ACE) there is no phi for the first
residue as the definition is C-1, N, CA, C (in contrast to what the
ptraj tutorial says) so that for residue 11 you would want:
dihedral phi_11 :10.C :11.N :11.CA :11.C out __phi_11
(2) Are you sure that the coordinates being read in and the prmtop
correspond?
You can check this by:
ptraj prm.top << EOF
trajin parm.coords
trajout test_pdb pdb
EOF
and then look at the resulting pdb file to see if it makes sense. Also,
if there is a mismatch between the number of coordinates read in from the
PDB and what ptraj expects, a warning is printed to the ptraj output.
(3) Are you sure the atom selections were recognized? You can see this by
looking at the output:
PTRAJ: dihedral phi2 :1.C :2.N :2.CA :2.C out phi2.dat
Mask [:1.C] represents 1 atoms
Mask [:2.N] represents 1 atoms
Mask [:2.CA] represents 1 atoms
Mask [:2.C] represents 1 atoms
...
ACTIONS
1> DIHEDRAL: saved to array named phi2
Atom selection 1 is :1.C
Atom selection 2 is :2.N
Atom selection 3 is :2.CA
Atom selection 4 is :2.C
Data will be dumped to a file named phi2.dat
(4) Another way to test is to simply build a default peptide with LEaP.
If coordinates are not specified, LEaP builds an extended (fully beta)
peptide with angles very close to 180.0 for phi/psi (which you can compare
with ptraj and vmd)
tleap -f $AMBERHOME/dat/leap/cmd/leaprc.ff94
aa = sequence { ACE ALA ALA NME}
saveamberparm aa prmtop inpcrd
quit
ptraj prmtop << EOF
trajin inpcrd
dihedral psi1 :1.C :2.N :2.CA :2.C out phi1.dat
dihedral phi1 :2.N :2.CA :2.C :3.N out psi1.dat
dihedral psi2 :2.C :3.N :3.CA :3.C out phi2.dat
dihedral phi2 :3.N :3.CA :3.C :4.N out psi2.dat
EOF
head *.dat
==> phi1.dat <==
1.00 -179.999776
==> psi1.dat <==
1.00 179.999921
==> phi2.dat <==
1.00 -179.999769
==> psi2.dat <==
1.00 179.999923
If you try this simple test and it doesn't work, please let me know what
version of the ptraj program you are using (i.e. off the WWW version 6.?
or integrated with what version of AMBER) and I can try to track it down.
--tom
\-/ Thomas E. Cheatham, III (Assistant Professor) College of Pharmacy
-/- Departments of Med. Chem. and of Pharmaceutics and Pharm. Chem.
/-\ Adjunct Asst Prof of Bioeng.; Center for High Performance Computing
\-/ University of Utah, 30 S 2000 E, SH 201, Salt Lake City, UT 84112
-/-
/-\ tec3.utah.edu (801) 587-9652; FAX: (801) 585-9119
\-/ BPRP295A
http://www.chpc.utah.edu/~cheatham
-----------------------------------------------------------------------
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:06 PST
