[gmx-users] Combined essential dynamics

[Tsjerk Wassenaar]

Hi Ram,

Since you only use frames from 1ns - 5ns from the simulation, you'll 
miss the identical projections the starting structures would give. If 
you want to know how your individual trajectories run in the plane 
defined by these two eigenvectors, project the whole trajectories. Both 
projections should indeed start from the same point if the starting 
structures are identical.

Cheers,

Tsjerk


Sekar Ramachandran wrote:

>Hi,
>I am trying to carry out a COMB_ED analysis of two
>forms of the  protein. I have carried out 5ns
>simulations of the two proeins. I want to use 1 ns to
>5ns data from the two trajctories to run an essential
>dynamics analysis. I created a tpr file for both
>protiens containing only the CA atoms. I concatenated
>the two trajectories using the folloing command:
>trjconv -f protein1.trr protein2.trr -o comb.trr
>-settime
>I set the firsttime as zero (it was 1 ns int he
>original file).
>on the prompt for last time : I entered c. I obtained
>a trajectory file with 8 ns of data.
>I ran g_covar on this with the following command:
>g_covar -f comb.trr s protein1.tpr -o eigenval.xvg -v
>eigenvec.trr.
>I ran g_anaeig using the above eigenvectors and the
>indivdual trajectories to obtain projections of the
>indivdual trajectories of the two forms onto the first
>two eigenvectors.
>g_anaeig -v eigenvec.trr -f protein1.trr -s
>protein1.tpr -2d 2dproj.xvg -b 1000 -e 5000 -first 1
>-last 2.
>When  I look at the 2D projections of the two
>trajctories, they seem to be really far apart form
>each other. This seems to be odd because the starting
>structures for these forms are almost identical and  I
>would have expected that at the beginning of the
>simulations the two projections would be close to eact
>other. (the rmsd between the two starting structures
>is 0.0 angstroms). This makes me believe that there is
>something wrong with the way I ahve carried out my
>analysis. I would be really grateful to anyone who cna
>help me with this.
>Thanks a lot and am really sorry about this elaborate
>description.
>Ram
>
>
>=====
>Sekar Ramachandran
>Research Associate
>316 Baker Laboratories
>DEpartment of Chemistry and Chemical Biology
>Cornell University
>Ithaca 14853
>phone:607-254-6767
>
>
>		
>_______________________________
>Do you Yahoo!?
>Win 1 of 4,000 free domain names from Yahoo! Enter now.
>http://promotions.yahoo.com/goldrush
>_______________________________________________
>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.
>  
>


-- 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-- :)
-- :) 	Tsjerk A. Wassenaar, M.Sc.
-- :) 	Molecular Dynamics Group
-- :) 	Dept. of Biophysical Chemistry
-- :) 	University of Groningen
-- :) 	Nijenborgh 4
-- :) 	9747 AG Groningen
-- :) 	The Netherlands
-- :)
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-- :)
-- :) 	Hi! I'm a .signature virus!
-- :) 	Copy me into your ~/.signature to help me spread!
-- :)
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~