After adding water you can do energy minimization (EM) in two steps:<br><br>1. Constrain the protein backbone and do EM.<br>2. Now do EM on the full system.<br>3. Run a short MD simulation by constraining the protein backbone.<br>
The above three steps will help hydrate the protein molecule with minimal distortion of protein structure.<br><br>4. Now run a MD on full system.<br><br>for details looks here:<br><a href="http://www.google.com/url?sa=t&source=web&ct=res&cd=2&ved=0CBcQFjAB&url=http%3A%2F%2Feugen.leitl.org%2Fchem%2Fkerrigje%2Fpdf_files%2Ffwspidr_tutor.pdf&ei=jOPzS8a3Lab2MdX1_aAO&usg=AFQjCNGB_3mXSQRHuqehBSHXsRyXP1Gymg&sig2=bY3NqXHmruR7eSLVyAuCHQ">http://www.google.com/url?sa=t&source=web&ct=res&cd=2&ved=0CBcQFjAB&url=http%3A%2F%2Feugen.leitl.org%2Fchem%2Fkerrigje%2Fpdf_files%2Ffwspidr_tutor.pdf&ei=jOPzS8a3Lab2MdX1_aAO&usg=AFQjCNGB_3mXSQRHuqehBSHXsRyXP1Gymg&sig2=bY3NqXHmruR7eSLVyAuCHQ</a><br>
<br>-Gaurav<br><br><div class="gmail_quote">On Wed, May 19, 2010 at 8:18 AM, sonali dhindwal <span dir="ltr"><<a href="mailto:sonali11dhindwal@yahoo.co.in">sonali11dhindwal@yahoo.co.in</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;">
<table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font-family: inherit; font-style: inherit; font-variant: inherit; font-weight: inherit; font-size: inherit; line-height: inherit; font-size-adjust: inherit; font-stretch: inherit;" valign="top">
Sorry, but I couldnt get your question,<br>I have used this .mdp file for energy minimisation after addition of water and using <pre style="margin: 0em;">GROMOS96 43a1 force field :</pre><br>title = drg_trp<br>
cpp = /lib/cpp ; location of cpp on SGI<br>define = -DFLEX_SPC ; Use Ferguson’s Flexible water model [4]<br>constraints = none<br>integrator = steep<br>dt = 0.002 ; ps !<br>
nsteps =
2000<br>nstlist = 10<br>ns_type = grid<br>rlist = 0.9<br>coulombtype = PME ; Use particle-mesh ewald<br>rcoulomb = 0.9<br>rvdw = 1.0<br>fourierspacing = 0.12<br>fourier_nx = 0<br>
fourier_ny = 0<br>fourier_nz = 0<br>pme_order = 4<br>ewald_rtol = 1e-5<br>optimize_fft = yes<br>;<br>; Energy minimizing stuff<br>;<br>emtol =
1000.0<br>emstep = 0.01<br><br>I hope it will help you to guide me further<br>Thanks<br><div>--<br>Sonali Dhindwal</div><br><br>--- On <b>Wed, 19/5/10, Erik Marklund <i><<a href="mailto:erikm@xray.bmc.uu.se" target="_blank">erikm@xray.bmc.uu.se</a>></i></b> wrote:<br>
<blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Erik Marklund <<a href="mailto:erikm@xray.bmc.uu.se" target="_blank">erikm@xray.bmc.uu.se</a>><br>Subject: Re: [gmx-users] enegry minimisation<br>
To: "Discussion list for GROMACS users" <<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>Date: Wednesday, 19 May, 2010, 5:31 PM<div><div></div><div class="h5"><br><br><div>
sonali dhindwal skrev:<br>> Hello All<br>> This question may sound trivial to many, but as i am new to this field, please help.<br>> I want to ask a question regarding my previous query of distortion of protein strucutre after molecular dynamcs simulation.<br>
> I have noticed that after enegry minimisation using
steepest decent algorithm, using emtol of 1000 kJ mol^-1 nm^-1 , large amount of distortion occurs.<br>> So is it necessary to do enegry minimisation step before MD, because this is my modeled protein, and i have already done energy minimisation using different program and after that I have done refinement also.<br>
> Thanks and regards<br>> ^<br>> <br>> <br>> --<br>> Sonali Dhindwal<br>> <br>> <br>So how has your system setup changed since your previous EM? Addition of water? Cutoffs? PME?<br><br>-- -----------------------------------------------<br>
Erik Marklund, PhD student<br>Dept. of Cell and Molecular Biology, Uppsala University.<br>Husargatan 3, Box 596, 75124 Uppsala, Sweden<br>phone: +46 18 471 4537 fax: +46 18 511 755<br><a href="http://mc/compose?to=erikm@xray.bmc.uu.se" target="_blank">erikm@xray.bmc.uu.se</a> <a href="http://folding.bmc.uu.se/" target="_blank">http://folding.bmc.uu.se/</a><br>
<br>-- gmx-users mailing list <a href="http://mc/compose?to=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/search" target="_blank">http://www.gromacs.org/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="http://mc/compose?to=gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br></div></div></div></blockquote></td></tr></tbody></table><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/search" target="_blank">http://www.gromacs.org/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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br></blockquote></div><br>