<br><br><div class="gmail_quote">On Sat, Mar 5, 2011 at 12:22 AM, 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: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div class="im">
<br></div>
It is difficult to say for certain how OPLS should be used with MD, since it was originally designed for different software doing MC simulations. I commonly see rlist=rcoulomb=rvdw=1.0 with vdwtype = cutoff and coulombtype = PME in the literature.<div class="im">
<br>
<br></div>
During EM there are no temperatures. Please consult some textbook material if this concept is foreign to you.<div class="im"><br></div></blockquote><div> </div><div> I am sorry sir I meant during production MD run? sorry for the mistake.<br>
<br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
<br></div>
There is no single metric that, in isolation, will tell you if a simulation is equilibrated. You have to look at many variables and judge for yourself. First, is the ensemble stable (temperature, pressure, energy, whatever else)? Then, is the protein structure stable (RMSD, radius of gyration, secondary structure, hydrogen bonds, many more)? From your RMSD plots, I would not be convinced that either are fully equilibrated, as one is still trending upwards. More generally, as I said, you can't simply judge the simulation based on one criterion.<div>
<br></div></blockquote><div>The Temperature , kinetic, potential and total energy had equilibrated (I mean they were oscillating about a mean value). My <br>doubt in that was this equilibration had reached right from the beginning which is surprising as it has to take some time for the <br>
eqilibration to occur. But radius of gyration was quite bad there was no such sign of such equilibration.<br> <br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div> </div></blockquote><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
4. If MD is run on a protein which had intrinsic domain motion, how will that equilibrate during MD? i.e., how will the rmsd plot appear when it is equilibrated?<br>
</blockquote>
<br></div>
Depends entirely upon the types of motions that are involved. Domain motions are usually very slow, and you may not see such movement on the nanosecond time scale with conventional MD. Normal modes calculations are generally better-suited for these types of studies.<div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Will it equlibrate at all? will not the domain movement or movement of a large loop in the structure prevent it from equilibration?<br>
<br>
</blockquote>
<br></div>
Ideally, in the limit of infinite sampling, you would see your protein's structure interconvert between whatever conformations are expected at some predictable time interval. I doubt you'll ever be able to run a conventional simulation long enough to see such behavior for even the smallest protein domains, unless you spend several years collected many microseconds of data. You may be able to see such motion faster with implicit solvent, but then the kinetics are meaningless in an absolute sense.<br>
<br></blockquote><div>Yes I got this point. Thanks. So the question here is, does the cut offs that i have chosen had given rise to such anomalies or any other major factor is there?<br> </div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
-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><font color="#888888">
-- <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></font></blockquote><div><br>Thanking you <br>With Regards<br>M. Kavyashree <br>
</div></div><br>