Respected Sir,<br><br>Greetings from Pawan.<br>I did the Inflategro procedure for the POPC bilayer generated using genconf.<br>It took around 26 compressions for coming near the initial area (just above it).<br>The minimization were all converged to Fmax < 1350.<br>
If I decrease the Fmax less than this I am getting machine precision.<br>But when I proceeded with the final compressed structure for pr mdrun it gave lincs warnings and ended with segmentation fault.<br>As an alternative to this bilayer I used the DPPC bilayer (pre-equilibrated) which is given in GMX-Benchmark distribution.<br>
I carried out the same steps of Inflategro. I used the cutoff of 14 A in the inflation and compression step also.<br>In this case the the Steepest Descents converged to Fmax < 1000 in all the steps. The maximum force war never above 850.<br>
But when I did position restraint mdrun with the last compressed file I got Lincs warnings and Segementation fault after few steps (30 - 40 steps).<br>Can you please help how to proceed ?<br><br>Thanking you,<br><br>Yours sincerely,<br>
Pawan<br><br><br><div class="gmail_quote">On Mon, Mar 23, 2009 at 6:26 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;">
<br>
<br>
Pawan Kumar wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Respected Sir,<br>
<br><div class="im">
Greetings from Pawan.<br>
I read in the manual about the freeze-grps but dint find a good explanation about the usage.<br>
Can you please tell me how to use freeze-grps in all dimensions during the minimizations ?<br>
<br>
</div></blockquote>
<br>
If you search for "freezegrps" within the list archive (the "Search" feature of the Gromacs homepage), you will find several hundred results, many of which provide examples of what users are trying and what may or may not work. This should be instructive.<br>
<br>
In a general sense, freezegrps are just like anything else under the group theory in Gromacs. Make an index group for the group of atoms you want to freeze (i.e., headgroups) and the dimensions in which you want them frozen (x, y, z - yes(Y) or no(N) for each).<br>
<br>
If you are choosing the route of POPE, it is very difficult to get these systems built using InflateGRO. At least, that has been my experience. Prepare to spend a lot of time with trial and error and very careful equilibration.<br>
<br>
-Justin<br>
<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">
Thanking you,<br>
Pawan<br>
<br>
On Mon, Mar 23, 2009 at 5:58 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a><br></div><div><div></div><div class="h5">
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br>
<br>
Pawan Kumar wrote:<br>
<br>
Respected Sir,<br>
<br>
Thanks for all your help, suggestions and guidance.<br>
I have few more queries.<br>
Which type of box is appropriate in editconf - cubic or triclinic ?<br>
<br>
<br>
That depends entirely upon the dimensions of your system and what is<br>
adequate to accommodate the size of your embedded protein.<br>
<br>
<br>
I have read in one of your mails in the archives that POPE<br>
bilayer contains some bad interactions. As it contains more no.<br>
of lipids (340) it covers the full protein ( trimer which I am<br>
working on ) other than the few loop regions. How far that<br>
bilayer can be used or is it fine to use the biger bilayer of<br>
popc128a whcih I obtained by genconf ?<br>
<br>
<br>
The choice of lipid should not be based on a size convenience. It<br>
should be based on a valid model system, and experimental evidence<br>
to which you can correlate your results, if possible.<br>
<br>
POPE presented a terrible challenge, since during the in vacuo<br>
InflateGRO minimizations, the headgroups folded in on themselves to<br>
form intra-molecule hydrogen bonds that resulted in a collapse of<br>
the molecule. The solution was to use freeze-grps in all dimensions<br>
during these minimizations, and then equilibrate very carefully.<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></div></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<div class="im"><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></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div class="im"><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><br>
before posting!<br>
Please don't post (un)subscribe requests to the list. Use<br>
the www interface or send it to<br>
<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@gromacs.org</a>>.<div class="im"><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>
</div></blockquote>
<br>
-- <br>
========================================<div><div></div><div class="h5"><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>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">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 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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br>