Thanks Justin. I will try again. But please refer to some protocol if you know and one last question that before doing umbrella sampling simulation how can one be sure that the pulling is good and one should go ahead with selecting window and doing umbrella sampling. In the end when you see histo.xvg and profile.xvg , it cost a lot computationally and if you don't get it right these computational resources are wasted. <div>
Shahid nayeem<br><br><div class="gmail_quote">On Mon, Feb 13, 2012 at 8:12 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im"><br>
<br>
shahid nayeem wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks for quick reply. I have created mutant of a complex by changing interface residue in VMD. These mutant are experimentally known to show less binding affinity. I want to reproduce these results with umbrella sampling. Now I am sending profile and histo file for wt and mutant. Please suggest where i am wrong.<br>
</blockquote>
<br></div>
The PMF curves look poorly converged. Your reaction coordinates are not the same for both the WT and mutant (you appear to have a far shorter reaction coordinate for the mutant). The energy minimum is also ill-defined for the mutant.<br>
<br>
As for the reason behind these phenomena, I cannot say, nor do I have time to sort through your data and try to work it out for you. Refer to the literature, find similar protocols, and proceed from there.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Shahid Nayeem<div><div></div><div class="h5"><br>
<br>
On Mon, Feb 13, 2012 at 7:15 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br>
<br>
<br>
shahid nayeem wrote:<br>
<br>
Dear Justin<br>
I am doing umbrella pulling simulation of a protein complex wt<br>
and mutant. I expect mutant to give lower deltG value. I am<br>
attaching a tif file of energy vs time curve of wt and mutant<br>
protein on pulling simulation. These energies are obtained by<br>
g_energy and selecting 11 option which is COM pulling energy. In<br>
this curve the first peak decreases and again rises at longer<br>
time. How many windows should be selected. As expected the peak<br>
of COM pulling energy is lower in mutants. Please explain why<br>
the energy again rises at higher time. Should I use the windows<br>
upto 160ps only because thereafter in both curve there is rise<br>
in energy value. the pull code used is as follows.<br>
<br>
<br>
What you have obtained is a path-dependent energy that may or may<br>
not signify anything useful - it almost certainly does not. I<br>
cannot offer an explanation of the sharp increase towards the end of<br>
the simulation other than to speculate that your box is of<br>
insufficient size and you're encountering PBC issues.<br>
<br>
As for how many windows are necessary, it's also impossible to tell.<br>
You need enough windows to adequately sample the reaction<br>
coordinate. Thus, it is decided based on how far you need to<br>
separate the two species, how strong the force constant is during<br>
the US simulations, the nature of the interactions in the system and<br>
how fast they converge, etc.<br>
<br>
-Justin<br>
<br>
<br>
pull_geometry = distance pull_dim = Y N Y<br>
pull_vec1 = 0.75 0 1<br>
pull_start = yes pull_ngroups = 1<br>
pull_group0 = Chain_B pull_group1 = Chain_A pull_rate1 = 0.01 pull_k1 = 1000 <br>
<br></div></div>
-- ==============================<u></u>__==========<div class="im"><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></div>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>__<a href="http://vt.edu/Pages/Personal/justin" target="_blank"><u></u>vt.edu/Pages/Personal/justin</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 class="im"><br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<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><div class="im"><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<br></div>
<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a><div class="im"><br>
<<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<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>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>
<<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a>><br>
<br>
<br>
</blockquote><div><div></div><div class="h5">
<br>
-- <br>
==============================<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> | (540) 231-9080<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>==========<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></div>