[gmx-users] Re: dPCA
Mu Yuguang (Dr)
YGMu at ntu.edu.sg
Fri Oct 14 03:41:49 CEST 2005
Dear Gia, David,
I think in some cases dPCA definitely better then normal PCA, as our
paper (Proteins) described. One shortcoming of it is that just like the
NMR structure determination, only the short range parameters (dihedrals)
involved. So for big proteins some important global motions may be mixed
or equally weighted with local motions.
I do not believe the many minima shown by dPCA are artifacts.
Gia, sorry to reply your requests late.
Here also send you the fortran program I used for dPCA analysis.
Similarly as David said, it read a trajectory file with the format
time cos(a1) sin(a1) cos(a2) sin(a2) ....
In case you want compare.
Dr. Yuguang Mu
Division of Structural and Computational Biology
School of Biological Sciences
Nanyang Technological University
60 Nanyang Drive
From: gmx-users-bounces at gromacs.org
[mailto:gmx-users-bounces at gromacs.org] On Behalf Of David
Sent: Friday, October 14, 2005 4:45 AM
To: Gia Maisuradze
Cc: gmx-users at gromacs.org
Subject: [gmx-users] Re: dPCA
On Thu, 2005-10-13 at 13:36 -0700, Gia Maisuradze wrote:
> Thank you very much for your reply. As you said, dPCA shows more
> free energy landscape is more rugged than in "normal" PCA. According
> paper that I mentioned in my previous email, dPCA is more realistic,
> Cartesian PCA represents an artifacts due to the mixing of internal
> overall motion. On the other hand, there are other opinions, for
> Erik Lindahl mentioned that the resulting eigenvectors are no longer
> orthogonal in the cartesian norm and dPCA is problematic.
> I am going to test dPCA and compare the results to normal PCA. Maybe
> clarify some problems.
We are currently considering some other possibilities as well. A recent
paper by the Levy group used PCA on a set of internal distances (e.g.
Ca-Ca) which has the advantage that you don't need to do a LSQ
superposition, but also that long distances weigh differently than short
distances. A problem in any method is that even if you obtain orthogonal
eigenvectors they may still be coupled. If you find a general solution
for these problems it may be very useful indeed.
> With best wishes,
> ----- Original Message -----
> From: "David" <spoel at xray.bmc.uu.se>
> To: "Gia Maisuradze" <gia at chem.unr.edu>
> Cc: <gmx-users at gromacs.org>
> Sent: Thursday, October 13, 2005 1:03 PM
> Subject: Re: dPCA
> > On Thu, 2005-10-13 at 11:48 -0700, Gia Maisuradze wrote:
> >> Dear David,
> >> I have checked the changes in new Gromacs 3.3. It says that
> >> PCA is supported in g_angle (you are the author). Does it mean that
> >> can do dihedral principal component analysis for proteins? Is it
> >> same type of analysis that has been done by Gerhard Stock' group
> >> (PROTEINS, 58, 45, 2005)? If yes, I don't clearly understand how I
> >> do that. Normal PCA (in cartesian coordinates) can be done by using
> >> g_covar, where we build the covariance matrix and diagonalize it,
> >> get eigenvectors and eigenvalues. Then analyze them.
> >> Could you, please, explain how I can do dPCA on Gromacs, if it is
> >> possible.
> > Yes, you use g_angle -oc to produce a new trajectory file
> > cos(phi) sin(phi) cos(psi) sin(psi) etc.
> > This you run through g_covar.
> > According to some people this information is not really useful for
> > protein analysis, as it can show many false minima and it is hard to
> > back to real space. Barriers you see may be trivial due to
> > atoms.
> >> The same question I have asked Erik Lindahl, but he did not know
> >> dPCA in new Gromacs version.
> >> Thanks,
> >> Gia
> > --
> > David.
> > David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
> > Dept. of Cell and Molecular Biology, Uppsala University.
> > Husargatan 3, Box 596, 75124 Uppsala, Sweden
> > phone: 46 18 471 4205 fax: 46 18 511 755
> > spoel at xray.bmc.uu.se spoel at gromacs.org
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596, 75124 Uppsala, Sweden
phone: 46 18 471 4205 fax: 46 18 511 755
spoel at xray.bmc.uu.se spoel at gromacs.org http://xray.bmc.uu.se/~spoel
gmx-users mailing list
gmx-users at gromacs.org
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-request at gromacs.org.
-------------- next part --------------
A non-text attachment was scrubbed...
Size: 18444 bytes
More information about the gromacs.org_gmx-users