Justin,<br><br>It was a huge bilayer replacement on the magnitude
comparable with the protein size. In more detailes both of the lipid
leflets moved to the top and bottom side of the pb. The protein were
placed on the old place ( due to the posres) so it turned out isolated
from membrane :)<br><br><br>Also I've found that I could <font face="Arial"><font size="3"> use simulated annealing to slowly warm the system under <i>NPT</i> conditions</font></font><br>
<br>So in this case could I skip the nvt phase and do this mod NPT equilibration ?<br><br><br>also I've create posres for the lipid on Z direction only. I hopes this helps prevent such bilayer diffusion.<br>
<br>
James<br><br><div class="gmail_quote">2011/11/11 Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>></span><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div><br>
<br>
James Starlight wrote:<br>
<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>
What temperature for NVT phase should be in<br>
<br>
gen_temp = 323<br>
<br>
Does this parameere must be equal to the ref_t ?<br>
<br>
I've obtain strange results after NVT with<br>
ref_t = 297 297 297 <br>
for protein in DMPC bilayer.<br>
<br>
</blockquote>
<br></div>
You should generate velocities for the temperature you wish to use. Otherwise, equilibration will take longer than necessary. At worst, it may crash if your system is far from equilibrium as the temperature coupling algorithm may fail.<div>
<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
The bilayer with water parts diffused in up and down direction of the box relative protein<br>
Why this should occur ? Might I alpy posres for lipid in NVT phase ?<br>
<br>
</blockquote>
<br></div>
Given that I know nothing about the magnitude of such motion, I can't really say. If there are severe distortions or displacements, restraining headgroup atoms may be necessary for the initial phases of equilibration. Usually such restraints should not be necessary for most well-built systems.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>
2011/11/11 Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>><br>
<br>
<br>
<br>
James Starlight wrote:<br>
<br>
Thanks Justin<br>
<br>
I'll try to do that things.<br>
<br>
Some addition questions<br>
<br>
1) About nvt and apt equilibration<br>
<br>
As I understood ref_t of the system must be equal to<br>
temperature of the phase transition of the specific LIPID.<br>
But in the npt and nvt.mdp files there are severl enties for<br>
ref_t ( for some system groups in index.ndx )<br>
For my system composed of protein dmpc and water :<br>
<br>
ref_t = 323 297 323 ;<br>
<br>
Does all this values in mdp file must be equal ( to the<br>
temperature of the lipid phase transition)<br>
ref_t = 297 297 297<br>
or each value must correspond to the temperature of phase<br>
transition of the individual element ?<br>
ref_t = 323 297 323 ;<br>
<br>
Does thic correct for both NPT and NVT ensembles ?<br>
<br>
<br>
It is incorrect to set different groups at different temperatures.<br>
Whatever result you obtain would be completely unrealistic.<br>
<br>
-Justin<br>
<br>
<br>
James<br>
<br>
2011/11/11 Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a><br></div>
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a><div><div></div><div><br>
<mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>>><br>
<br>
<br>
<br>
<br>
James Starlight wrote:<br>
<br>
Justin,<br>
<br>
Could you tell me what difference in inflategro parameters (<br>
cut-off for lipid removing and scaling factors) should I make<br>
for insertion protein in another bilayes in comparison tu<br>
tutorial ?<br>
<br>
<br>
I've never made any changes; it's always worked just fine.<br>
<br>
<br>
Now I'm working with DMPC. I've inserted symmetrical alpha<br>
helices protein in this bilayer but inflategro deleate 2<br>
lipid<br>
from upper and just 1 from lower layer. It's strange<br>
because I<br>
suppose that equal bilayer molecules would be removed :o<br>
<br>
<br>
The number of lipids that are removed will depend on the starting<br>
configuration of the membrane and thus where lipids are<br>
relative to<br>
the protein. Theoretically, a symmetrical protein should<br>
delete an<br>
equal number of lipids from each leaflet, if the geometry of<br>
the top<br>
and bottom leaflets is comparable. Either try placing the<br>
protein<br>
at a different location to achieve equivalent removal or manually<br>
delete a lipid and be sure to equilibrate adequately.<br>
<br>
<br>
<br>
On Inflategro's sites I havenot found such specifity for<br>
above<br>
parameters for different bilayers :(<br>
<br>
<br>
The cutoffs specified are simple grid search parameters. They do<br>
not necessarily need to be altered depending upon the lipid type.<br>
<br>
<br>
-Justin<br>
<br></div></div>
-- ==============================<u></u>____==========<div><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></div>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540)<br>
231-9080<br>
<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>__<u></u>__<a href="http://vt.edu/Pages/Personal/justin" target="_blank">vt.edu/Pages/Personal/justin</a><br>
<<a href="http://vt.edu/Pages/Personal/justin" target="_blank">http://vt.edu/Pages/Personal/<u></u>justin</a>><br>
<<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>_<u></u>_<a href="http://vt.edu/Pages/Personal/justin" target="_blank">vt.edu/Pages/Personal/justin</a><br>
<<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a>>><br>
<br>
==============================<u></u>____==========<div><br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<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> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><u></u>><br>
<a href="http://lists.gromacs.org/____mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/____<u></u>mailman/listinfo/gmx-users</a><br>
<<a href="http://lists.gromacs.org/__mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/__<u></u>mailman/listinfo/gmx-users</a>><div><br>
<br>
<<a href="http://lists.gromacs.org/__mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/__<u></u>mailman/listinfo/gmx-users</a><br>
<<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a>>><br>
Please search the archive at<br></div>
<a href="http://www.gromacs.org/____Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/____<u></u>Support/Mailing_Lists/Search</a><br>
<<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a>><div><br>
<br>
<<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a>>> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>><br></div>
<mailto:<a href="mailto:gmx-users-request@" target="_blank">gmx-users-request@</a>__<a href="http://gromacs.org" target="_blank">gr<u></u>omacs.org</a><div><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>>>.<br>
Can't post? Read<br></div>
<a href="http://www.gromacs.org/____Support/Mailing_Lists" target="_blank">http://www.gromacs.org/____<u></u>Support/Mailing_Lists</a><br>
<<a href="http://www.gromacs.org/__Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists</a>><div><div></div><div><br>
<<a href="http://www.gromacs.org/__Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists</a><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a>>><br>
<br>
<br>
<br>
-- ==============================<u></u>__==========<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> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem." target="_blank">http://www.bevanlab.biochem.</a>__<a href="http://vt.edu/Pages/Personal/justin" target="_blank"><u></u>vt.edu/Pages/Personal/justin</a><br>
<<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a>><br>
<br>
==============================<u></u>__==========<br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<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/__<u></u>mailman/listinfo/gmx-users</a><br>
<<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a>><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/__Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists/Search</a><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a>> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>>.<br>
Can't post? Read <a href="http://www.gromacs.org/__Support/Mailing_Lists" target="_blank">http://www.gromacs.org/__<u></u>Support/Mailing_Lists</a><br>
<<a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a>><br>
<br>
<br>
</div></div></blockquote><div><div></div><div>
<br>
-- <br>
==============================<u></u>==========<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.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>
==============================<u></u>==========<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/<u></u>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/<u></u>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/<u></u>Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>