Dear Sir, <br><br>Yes it is not stated in your tutorial. I was mentioning <br>about the superposition of trajectory and visualization <br>in pymol. <br>Ok I will try that. Thank you for the suggestions sir. <br><br>Thanking you<br>
With Regards<br>M. Kavyashree<br><br><br><br><br><div class="gmail_quote">On Thu, Jun 23, 2011 at 2:43 PM, Tsjerk Wassenaar <span dir="ltr">&lt;<a href="mailto:tsjerkw@gmail.com">tsjerkw@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi Kavya,<br>
<br>
For certain I do not state in my tutorial that one should use g_nmtraj<br>
for this. It is intended for analysis of normal modes analysis, where<br>
the eigenvalues are related to the frequencies of the harmonic<br>
motions. You&#39;re using it on the modes with largest eigenvalues, which<br>
are interpreted as high frequencies, consequently giving small<br>
amplitudes. This explains that you&#39;re not seeing much.<br>
<br>
To get an idea of the extent of the motion, look at the extreme<br>
projections. To see how the trajectory moves along a certain mode, use<br>
g_anaeig -filt<br>
<br>
Cheers,<br>
<font color="#888888"><br>
Tsjerk<br>
</font><div><div></div><div class="h5"><br>
On Thu, Jun 23, 2011 at 10:56 AM, Kavyashree M &lt;<a href="mailto:hmkvsri@gmail.com">hmkvsri@gmail.com</a>&gt; wrote:<br>
&gt; Dear Sir,<br>
&gt;<br>
&gt; 1 simulation for 100ns -<br>
&gt; ED analysis proceeded as follows:<br>
&gt;<br>
&gt; g_covar -s &lt;input&gt;.tpr -f .&lt;input&gt;.xtc -o eigenvectors.xvg -v<br>
&gt; eigenvalues.trr -xpma covar.xpm<br>
&gt; g_anaeig -s &lt;input&gt;.tpr -f &lt;input&gt;.xtc -v eigenvectors.trr -eig<br>
&gt; eigenvalues.xvg -proj proj-evi.xvg -extr evi.pdb -rmsf rmsf-evi.xvg -first i<br>
&gt; -last i<br>
&gt; (here i refers to the eigen vector index)<br>
&gt; Thus generated 5 eigen vectors.<br>
&gt;<br>
&gt; Visualised the trajectory along these eigen vectors using<br>
&gt;<br>
&gt; g_nmtraj -s &lt;input&gt;tpr -v eigenvector.trr -o &lt;output&gt;pdb -eigen i  -phases x<br>
&gt; -temp x -amplitude x -nframes x<br>
&gt; (i referes to eigen vectors. In this I was not clear about the significance<br>
&gt; of phases, amplitude)<br>
&gt;<br>
&gt; I visualised the pbd output in pymol there was very minute movement of some<br>
&gt; loop regions.<br>
&gt;<br>
&gt; on the other side-<br>
&gt; When I visualised the .xt file output as pdb according to  Dr.Tjerk&#39;s<br>
&gt; tutorial or<br>
&gt; by superposing all the structures in trajectory over initial structure in<br>
&gt; pymol<br>
&gt; there was a significant degree of movement in those loops.<br>
&gt;<br>
&gt; so my question here is why these two does not correlate. I agree that<br>
&gt; trajectory<br>
&gt; along one eigenvector might not give movements of all the regions that is<br>
&gt; observed<br>
&gt; in superposed trajectory but whichever region it shows movements is<br>
&gt; extremely less<br>
&gt; compared to that when viewed from superposed structures in trajectory.<br>
&gt;<br>
&gt; I hope this is not confusing.<br>
&gt;<br>
&gt; Thanking you<br>
&gt; With regards<br>
&gt; M, Kavyashree<br>
&gt;<br>
&gt; On Thu, Jun 23, 2011 at 1:27 PM, Tsjerk Wassenaar &lt;<a href="mailto:tsjerkw@gmail.com">tsjerkw@gmail.com</a>&gt; wrote:<br>
&gt;&gt;<br>
&gt;&gt; Hi Kavya,<br>
&gt;&gt;<br>
&gt;&gt; For us to say anything sensible about it, we should at least know<br>
&gt;&gt; exactly what you&#39;ve tried. Copy-paste the commands exactly as you<br>
&gt;&gt; issued them, and provide the parts of the output that seem relevant.<br>
&gt;&gt;<br>
&gt;&gt; Cheers,<br>
&gt;&gt;<br>
&gt;&gt; Tsjerk<br>
&gt;&gt;<br>
&gt;&gt; On Thu, Jun 23, 2011 at 9:54 AM, Kavyashree M &lt;<a href="mailto:hmkvsri@gmail.com">hmkvsri@gmail.com</a>&gt; wrote:<br>
&gt;&gt; &gt; Dear Sir,<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; I tried fitting the proteins of the trajectory in pymol<br>
&gt;&gt; &gt; as mentioned in your (Dr. Tsjerk&#39;s) tutorial, but later<br>
&gt;&gt; &gt; I tried using the trjconv -fit rot+trans to fit the proteins<br>
&gt;&gt; &gt; in the trajectory as you had mentioned.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; I do not observe this degree of movement in protein<br>
&gt;&gt; &gt; when I view the trajectory along any 1-5 eigen vectors.<br>
&gt;&gt; &gt; This trajectory along any of these eigen vector is very minute<br>
&gt;&gt; &gt; and cannot be compared to that obtained after superposing<br>
&gt;&gt; &gt; the structure in the trajectory. Why is it so?<br>
&gt;&gt; &gt; Please let me know if I am understanding the concept wrong.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Thanking you<br>
&gt;&gt; &gt; With regards<br>
&gt;&gt; &gt; M. Kavyashree<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; On Thu, Jun 23, 2011 at 11:56 AM, Kavyashree M &lt;<a href="mailto:hmkvsri@gmail.com">hmkvsri@gmail.com</a>&gt;<br>
&gt;&gt; &gt; wrote:<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Thank you Sir!<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; With regards<br>
&gt;&gt; &gt;&gt; M. Kavyashree<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; On Thu, Jun 23, 2011 at 11:50 AM, Tsjerk Wassenaar &lt;<a href="mailto:tsjerkw@gmail.com">tsjerkw@gmail.com</a>&gt;<br>
&gt;&gt; &gt;&gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; trjconv -fit rot+trans<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; Cheers,<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; Tsjerk<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; On Jun 23, 2011 8:12 AM, &quot;Kavyashree M&quot; &lt;<a href="mailto:hmkvsri@gmail.com">hmkvsri@gmail.com</a>&gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; Dear Users,<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; Are there any tool for superposing the trajectory<br>
&gt;&gt; &gt;&gt;&gt; structures form MD. Please correct me if I am asking<br>
&gt;&gt; &gt;&gt;&gt; any illogical question.<br>
&gt;&gt; &gt;&gt;&gt; My previous question was regarding the trjconv output<br>
&gt;&gt; &gt;&gt;&gt; pdb trajectory, is there a way to superpose all these<br>
&gt;&gt; &gt;&gt;&gt; structures?<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; Thank you<br>
&gt;&gt; &gt;&gt;&gt; With regards<br>
&gt;&gt; &gt;&gt;&gt; M. Kavyashree<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; On Wed, Jun 22, 2011 at 10:52 PM, Kavyashree M &lt;<a href="mailto:hmkvsri@gmail.com">hmkvsri@gmail.com</a>&gt;<br>
&gt;&gt; &gt;&gt;&gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt; &gt; &gt; Dear users, &gt; &gt;    I...<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; --<br>
&gt;&gt; &gt;&gt;&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;&gt; &gt;&gt;&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt;&gt; &gt;&gt;&gt; Please search the archive at<br>
&gt;&gt; &gt;&gt;&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
&gt;&gt; &gt;&gt;&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt;&gt; &gt;&gt;&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt;&gt; &gt;&gt;&gt; Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; --<br>
&gt;&gt; &gt;&gt;&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;&gt; &gt;&gt;&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt;&gt; &gt;&gt;&gt; Please search the archive at<br>
&gt;&gt; &gt;&gt;&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
&gt;&gt; &gt;&gt;&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt;&gt; &gt;&gt;&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt;&gt; &gt;&gt;&gt; Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; --<br>
&gt;&gt; &gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;&gt; &gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt;&gt; &gt; Please search the archive at<br>
&gt;&gt; &gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
&gt;&gt; &gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt;&gt; &gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt;&gt; &gt; Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
&gt;&gt; &gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; --<br>
&gt;&gt; Tsjerk A. Wassenaar, Ph.D.<br>
&gt;&gt;<br>
&gt;&gt; post-doctoral researcher<br>
&gt;&gt; Molecular Dynamics Group<br>
&gt;&gt; * Groningen Institute for Biomolecular Research and Biotechnology<br>
&gt;&gt; * Zernike Institute for Advanced Materials<br>
&gt;&gt; University of Groningen<br>
&gt;&gt; The Netherlands<br>
&gt;&gt; --<br>
&gt;&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt;&gt; Please search the archive at<br>
&gt;&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
&gt;&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt;&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt;&gt; Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt; Please search the archive at<br>
&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt; Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
&gt;<br>
<br>
<br>
<br>
--<br>
Tsjerk A. Wassenaar, Ph.D.<br>
<br>
post-doctoral researcher<br>
Molecular Dynamics Group<br>
* Groningen Institute for Biomolecular Research and Biotechnology<br>
* Zernike Institute for Advanced Materials<br>
University of Groningen<br>
The Netherlands<br>
--<br>
gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>