<p>Hi Vijayan R.,</p>
<p>Eigenvectors are specific up to a sign. So it's not said the start will end up as the negative extreme. It seems your interpolation just shows what happens in reverse.<br>
As a sidenote, do check whether it is a biologically relevant motion, or merely relaxation.</p>
<p>Cheers,</p>
<p>Tsjerk </p>
<p><blockquote type="cite">On Dec 21, 2011 4:24 PM, "R.S.K.Vijayan" <<a href="mailto:biovijayan@gmail.com">biovijayan@gmail.com</a>> wrote:<br><br><br>Hi all<br><br>I carried out a distance analysis between two domains with respect to time based on the center of mass. Subsequently i also carried out a PCA based essential dynamics to identify the global/biologically relevant motions. The strange and inconsistent result i find here is <br>
<br>when i calculate the distance using the MD (progressive fitted) trajectory, I find that the two domains of interest move apart with respect to time. <br><br>on the contrary a ED followed by a smoothed interpolation between the two extreme motions revelas that the two domains move closer to each other in a concerted fashion.<br>
<br>does any one has an explanation for this discrepancy between ED motions and the MD motion ???<br><br>Thanking you in advance <br><br>Regards<br>Vijayan.R<br><br><br><br><div class="gmail_quote">On Tue, Dec 20, 2011 at 2:28 AM, bipin singh <span dir="ltr"><<a href="mailto:bipinelmat@gmail.com" target="_blank">bipinelmat@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Thanks a lot for your suggestions.<br><br><div class="gmail_quote">On Mon, Dec 19, 2011 at 21:00, David Mobley <span dir="ltr"><<a href="mailto:dmobley@gmail.com" target="_blank">dmobley@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Yes, changing the net charge of the system is something that is rather complicated in fact (one can plunge ahead and do it while ignoring the complications, but the results will typically be rather system-size dependent and essentially wrong). For more details refer to the Kastenholz and Hunenberger papers from J. Chem. Phys (2006 or 2007) and references therein. The good news is that for simple spherical ions K&H have worked out the relevant corrections. The bad news is that if you move away from simple spherical ions you've got problems.<div>
<br></div><div>David Mobley</div><div><br></div><div><div><div><br><br><div class="gmail_quote">On Fri, Dec 16, 2011 at 6:41 AM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><br>
<br>
bipin singh wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Hello,<br>
<br>
I am willing to study the free energy of binding of a cation (Ca++) to<br>
the protein and I am following the free energy tutorial<br>
provided by Justin<br>
(<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/free_energy" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin/<u></u>gmx-tutorials/free_energy</a><br>
).<br>
Please let me know whether the settings for this type of study would<br>
be same as given in the tutorial for ligand-Protein<br>
binding free energy calculation or it need some different approach:<br>
</blockquote>
<br></div>
By decoupling a Ca2+ ion, you are removing 2 charges from the system. I don't know how to properly treat such a case (perhaps someone else can comment), but likely you'll find such topics in the literature. A better approach may be umbrella sampling, but again the literature should point you to reasonable methodology. I'm sure others have dealt with such questions before.<br>
<br>
-Justin<div><br>
<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
The setting from the ligand-Protein binding free energy calculation<br>
are given as:<br>
<br>
<br>
van der Waals coupling:<br>
<br>
sc-alpha = 0.5 ; use soft-core for LJ (de)coupling<br>
sc-sigma = 0.3<br>
sc-power = 1<br>
couple-moltype = LIG<br>
couple-intramol = no<br>
couple-lambda0 = none ; non-interacting dummy in state A<br>
couple-lambda1 = vdw ; only vdW terms on in state B<br>
<br>
Coulombic coupling:<br>
<br>
sc-alpha = 0 ; soft-core during (dis)charging can<br>
be unstable!<br>
sc-sigma = 0<br>
couple-moltype = LIG<br>
couple-intramol = no<br>
couple-lambda0 = vdw ; only vdW terms in state A (the<br>
previous state B is now A)<br>
couple-lambda1 = vdw-q ; all nonbonded interactions are on in state B<br>
<br>
<br>
</blockquote>
<br>
-- <br></div>
==============================<u></u>==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | <a href="tel:%28540%29%20231-9080" value="+15402319080" target="_blank">(540) 231-9080</a><br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>
==============================<u></u>==========<div><div><br>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>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/<u></u>Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div></div></div><span><font color="#888888">-- <br>David Mobley<br><a href="mailto:dmobley@gmail.com" target="_blank">dmobley@gmail.com</a><br>
<a href="tel:504-383-3662" value="+15043833662" target="_blank">504-383-3662</a><br><br><br><span><font color="#888888">
</font></span></font></span></div><span><font color="#888888">
<br>--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">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'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" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br></font></span></blockquote></div><span><font color="#888888"><br>
<br clear="all"><font color="#888888"><br>-- <br><span style="color:rgb(51,51,255)">-----------------------</span><br style="color:rgb(51,51,255)">
<i><span style="color:rgb(51,51,255)">Regards,</span></i><br style="color:rgb(51,51,255)"><span style="color:rgb(51,51,255);font-family:garamond,serif">Bipin Singh</span><br><br>
</font></font></span><font color="#888888"><br>--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">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'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" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br></font></blockquote></div><br>
<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'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't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br></blockquote></p>