AMBER Archive (2006)

Subject: AMBER: ptraj PCA

From: Vlad Cojocaru (Vlad.Cojocaru_at_eml-r.villa-bosch.de)
Date: Tue Jan 24 2006 - 10:02:20 CST

Dear ambers,
   When using the inputs below to perform PCA on a protein trajectory
(file.mdcrd containing a full protein but no ions or water), ptraj
always fails to create the projection file (${fname}.ppj -- second
step), the error being listed below. Obviously, the error has something
to do with the mask because if I replace @~H* with @CA, the projection
file is properly created. However I would like to be able to use all the
heavy atoms so that I can visualize the results in IED using a pdb file
with all the heavy atoms in the protein (I believe, IED doesn't allow me
to use a pdb of the entire protein but eigenvectors calculated only for
the CA atoms).
   If somebody has experienced such a problem, could you give me an
advice on how to deal with this?
   I have to mention also that the modesfile ${fname}.pev created in the
first step looks OK with 25 eigenvectors output and in the ptraj log
file of the 2nd step the mask seems to be properly recognized (PTRAJ:
projection modes ${fname}.pev out ${fname}.ppj beg 1 end 25 @~H* \
Mask [~@ H*] represents 3736 atoms), with 3736 being the appropriate
number of heavy atoms in my protein. The error comes on STDOUT output.

Thanks in advance
vlad

/input1:/

trajin file.mdcrd
rms first @~H*
matrix covar name mc @~H*
analyze matrix mc out ${fname}.pev vecs 25

/input2:/

trajin file.mdcrd
projection modes ${fname}.pev out ${fname}.ppj beg 1 end 25 @~H*

/ptraj error (file ${fname}.ppj is blank -- only the header present)/

FYI: Number of read evecs is 1, number of requested evecs is 25
FYI: Number of read evecs is 1, number of requested evecs is 25
WARNING in ptraj(), transformProjection: no. of atom coords does not
agree with no. of vect elementsignoring command 

-- 
Dr. Vlad Cojocaru
EML Research gGmbH
Molecular and Cellular Modeling Group
Schloss-Wolfsbrunnenweg 33
69118 Heidelberg, Germany
Phone: +49-6221-533266
Fax: +49-6221-533298
e-mail: Vlad.Cojocaru_at_eml-r.villa-bosch.de
http://projects.villa-bosch.de/mcm/people/cojocaru/

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu