Thank you.<br><br>no,I think the problem is not zero point.<br>Actually I maean some thing like this one.please have a look at new attached file.<br>I just added 3 new windows in 0.3 - 0.8<br>But the whoe of my curve shifted upward.<br>
<br>There is another question:<br>Actually I pulled my drug from 0.18 - 3.18<br>then,extracted some windows.umbrella0 for 0.18 and umbrella20 for 3.18<br>Besides as we know drug is oscilating in it's windows.<br>so ,it must cover places as like as 0.12 and 3.22 too.<br>
Am I right?<br><br>But looking at curve it seems my drug couln't cover the region around 0.18<br>besides it ha covered completely(untill 0.03) the region around 0.18<br><br>What is the problem?<br><br>Since I started from 0.18 (which energy is negative there),I am in doubt about starting point.<br>
Because it can be eighter 0.03 OR 0.18<br><br>Sorry,my english is not very good<br>Thanks in advance for your reply<br><br><div class="gmail_quote">On Mon, May 30, 2011 at 7:14 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="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im"><br>
<br>
mohsen ramezanpour wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Thank you for your reply<br>
<br>
Actually my system is not the same as tutorial.<br>
</blockquote>
<br></div>
Then you should not equate the results, or be concerned when things look different. Every system is different.<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
I tried to obtain protein-drug delta G.<br>
I obtained that.<br>
Yes,you are right,since I don't have good computaional systems I couldn't run more than 1 ns for each windows.<br>
<br>
when there are some minimum and maximum in our curve ,How do I must evaluate Delta G?<br>
the difference between minimum and maximum values OR difference between values corresponding to<br>
0 and last windows(where they must be converged(be flatted))?<br>
</blockquote>
<br></div>
Revisit fundamental thermodynamics. DeltaG is the difference in free energy between the final and initial states. It doesn't matter how you get there - DeltaG is a state function. Like I said, there may be local maxima or minima along the reaction coordinate, but the DeltaG for the process is still the highest point (where the PMF has leveled off, presumably in solution) minus the energy value at your starting point (your global minimum, it would seem).<div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Besides when I add a new windows in places which there are not enough sampling ,a vertical shift is resultedin the whole of my curve.<br>
Is not it strong(wonderfull?)?<br>
<br>
</blockquote>
<br></div>
I don't completely understand this comment. You can control the position of the zero value of the PMF using g_wham -zprof0. Since DeltaG is relative, you can choose a location for the zero point however you like. I generally set the zero point somewhere out where the PMF curve has leveled off (i.e., no interaction between the pulled group and the reference). If your curve is shifting dramatically when adding windows, you probably aren't anywhere close to achieving a stable result and you will need substantially more simulation time in each window.<div>
<div></div><div class="h5"><br>
<br>
-Justin<br>
<br>
-- <br>
========================================<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.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<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 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>
</div></div></blockquote></div><br>