[gmx-users] how to arbitrarily increase the XY size of an existing bilayer system (adding more lipids)
Christopher Neale
chris.neale at alum.utoronto.ca
Wed Feb 1 18:30:37 CET 2017
Another great idea. Thanks Xavier! I'll report back on what I eventually do in a few weeks with info on issues and timescales necessary for relaxation so that it's available if anyone ever wants to do it in the future.
________________________________________
From: gromacs.org_gmx-users-bounces at maillist.sys.kth.se <gromacs.org_gmx-users-bounces at maillist.sys.kth.se> on behalf of XAvier Periole <x.periole at rug.nl>
Sent: 01 February 2017 11:18:36
To: gmx-users at gromacs.org
Subject: Re: [gmx-users] how to arbitrarily increase the XY size of an existing bilayer system (adding more lipids)
What about tilling 2x2 and then removing lipids that have x&y larger than a cutoff … then minimisation+equilibration of the interface even atomistic would be quite trivial.
X-
> On 01 Feb 2017, at 17:14, Christopher Neale <chris.neale at alum.utoronto.ca> wrote:
>
> Dear Users:
>
> I have some atomistic systems of a membrane protein embedded in a lipid bilayer. I currently have N lipids and I would like to increase that to 1.5N or 2N lipids (distributed equally in the bilayer plane) without disturbing the existing structure. Increasing to 4N lipids would be relatively easy by tiling the exisintg lipids 2x2 in XY. I am writing to ask if anybody has ever done this and if they have any suggestions.
>
> The ideas that I have at present are:
> a) add lipids and equilibrate (using restaints on existing lipids and protein). Increase in size could be by simply increasing XY dimensions and then using gmx genbox, knowing that I will end up with a bilayer defect at the edges that I will have to heal by MD.
> b) convert to a CG model, add lipids and equilibrate, then back to an AA model (using restaints on existing lipids and protein)
>
> I realize it sounds like an absurd requirement, but for scientific reasons in this particular case I can not simply start over by embedding the protein in a new, larger bilayer (i.e., it’s not just about saving compute time when generating the larger systems).
>
> Thank you for any advice,
> Chris.
> --
> Gromacs Users mailing list
>
> * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.
--
Gromacs Users mailing list
* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.
More information about the gromacs.org_gmx-users
mailing list