Hello Justin Sir,<br><br>Greetings from Pawan<br>Thanks for your valuable suggestion and reply.<br>After inserting the protein in the bilayer using genbox I have minimized the whole system without using any position restraints (i.e. define = -DFLEXIBLE in em.mdp file). I used vanderwaal's distance parameter ( -vdwd of 0.6 ) in the genbox step.<br>
After running mdrun for energy minimization I got the output as :<br>Steepest Descents converged to Fmax<2250 in 14 steps.<br>Potential energy = - 6.9484700e+05<br>Maximum force = 2.2114819e+03 on atom 34277<br>Norm. of force = 5.0103039e+04<br>
<br>I tried decreasing the emtol value in the em.mdp file but it ended with machine precision.<br>I have read in literature that 5000 steps of Steetest Descents run is required after inserting the protein in the bilayer which should be followed by atleast 1000 steps of conjugate gradients. How can I accomplish this ? Is there any parameter to be given in the em.mdp file ? I use steep as the integrator in the mdp file for energy minimization. Please help with some suggestions.<br>
Thanks in advance.<br><br>Thanking you,<br>Pawan<br><br><br><br><div class="gmail_quote">On Thu, Mar 19, 2009 at 5:41 PM, <span dir="ltr"></span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Pawan Kumar wrote:<br>
> Hello Justin Sir,<br>
><br>
> Thanks for your reply.<br>
> I changed the dt value as per the suggestion.<br>
> But after 200 steps the same kind of warnings came like "pressure<br>
> scaling more than 1%" and "1-4 interactions" and then it stopped after<br>
> writing few pdb files. Then I decreased the dt value still lesser<br>
> (0.00001) and this time it continued with all the given steps but the<br>
> lipids' structure were distorted into small fragments.<br>
<br>
If your system is not properly energy-minimized, there is no change in dt that<br>
will help you. Ensure that your EM converged to reasonable values of Epot and<br>
Fmax before continuing.<br>
<br>
> The other query I have is what should the optimal value for these<br>
> parameters for membrane proteins : rlist, rcoulomb, and rvdw. I have<br>
> used the value of 1 before.<br>
<br>
These parameters depend upon the force field you are using. Read about the<br>
parameter set, as well as modifications that may be made to it. Studies have<br>
suggested that longer cut-off's should be used with membrane proteins. There is<br>
no substitute for a thorough literature search before starting a project!<br>
<br>
> One last question is how to overcome the lincs warnings in position<br>
> restraint mdrun.<br>
<br>
That will depend upon whether or not your system is properly constructed and<br>
energy-minimized.<br>
<br>
-Justin<br>
<br>
> Thanks in advance.<br>
><br>
> Thanking you,<br>
> Pawan<br>
><br>
> On Thu, Mar 19, 2009 at 5:01 PM,<br>
><br>
> Pawan Kumar wrote:<br>
> > Hi Xavier sir,<br>
> ><br>
> > Thanks for your valuable reply.<br>
> > How can I refine the non-bonded set-up ?<br>
><br>
> rlist, rcoulomb, rvdw as well as coulombtype. *Never* use cut-off<br>
> electrostatics for a membrane system (or really any other, for that<br>
> matter).<br>
> Your results will be far less accurate than with PME (this is in the<br>
> literature).<br>
><br>
> > And where can I specify the time step to 0.0001 ?<br>
><br>
> dt (read the manual).<br>
><br>
> -Justin<br>
><br>
> > I am new to gromacs.<br>
> > Sorry to ask such questions.<br>
> ><br>
> > Thanking you,<br>
> > Pawan<br>
> ><br>
> ><br>
><br>
><br>
><br>
> ------------------------------------------------------------------------<br>
><br>
> _______________________________________________<br>
> gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
> <a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.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>
<br>
--<br>
========================================<br>
<br>
Justin A. Lemkul<br>
Graduate Research Assistant<br>
ICTAS Doctoral Scholar<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>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
gmx-users mailing list<br>
<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.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>
<br>
End of gmx-users Digest, Vol 59, Issue 137<br>
******************************************<br>
</blockquote></div><br>