[gmx-users] PCoorA and Energies

[Tsjerk Wassenaar]

Hi Andrea,

You have a different energy landscape for each temperature. For a
given point in conformational space, the energy of your system will be
higher with higher temperatures. However, the energies will not only
be shifted, but the shift will depend on the position in
conformational space (so wells will become shallower and barriers will
be crossed more easily). In short, no, you can't just plot the
energies over the "principal coordinates" and construct an energy
landscape, you'll have to do it for the different temperatures
seperately (which would be nice for an illustration though). Note that
distance in PCoorA will be "conformational distance", i.e. the RMSD.

Ciao,

Tsjerk

On 3/20/07, andrea spitaleri <spitaleri.andrea at hsr.it> wrote:
> Hi guys,
> principal coordinate analysis (PCoorA) is a variant of the classical PCA. It is the covariance
> matrix of the distances. Becker et al. (JCP 1998,19,1255 and 2001,114,993) sentences that the PCoorA
> reflects the relationship between the conformations rather than the coordinates as in the classical
> PCA study. In his study he reported that plotting PCoorA-1 and/or PCoorA-2 (1st and 2nd principal
> coordinate axis) versus the Energies one can extract a potential energy surface in order to study
> the different weight of the various conformers.
> Now, I have performed a REMD (2ns x 15 replicas) of my small peptide and I did the usual analysis on
> the 30ns trr trajectory.
> I am wondering now whether plotting the PCoorA-1/2 versus the Energies from the REMD makes
> physically sense, since the total energy.edr is coming from different temperatures.
> Any hints? any suggestions?
>
> Thanks in advance
>
> Regards,
>
> andrea
>
>
> --
> -------------------------------
> Andrea Spitaleri PhD
> Dulbecco Telethon Institute
> c/o DIBIT Scientific Institute
> Biomolecular NMR, 1B4
> Via Olgettina 58
> 20132 Milano (Italy)
> http://biomolecularnmr.ihsr.dom/
> -------------------------------
>
> ********************************************************************
> Sostieni la ricerca del San Raffaele con il 5permille!
> E' SEMPLICE E NON COSTA NULLA.
> Basta indicare nell'apposito riquadro della dichiarazione dei
> redditi ("Ricerca sanitaria")
> il codice fiscale della
> Fondazione Centro S. Raffaele del Monte Tabor:
> 03 06 42 80 153 e ricordarsi di firmare.
> Se vuoi saperne di piu' scrivi a 5permille at hsr.it o vai sul sito
> www.5xmille.org
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>


-- 
Tsjerk A. Wassenaar, Ph.D.
Junior UD (post-doc)
Biomolecular NMR, Bijvoet Center
Utrecht University
Padualaan 8
3584 CH Utrecht
The Netherlands
P: +31-30-2539931
F: +31-30-2537623