[gmx-users] Resizing heterogeneous membrane (POPG/POPE mixture)
shayamra at post.tau.ac.il
Wed Jun 3 10:53:39 CEST 2009
By some fluke one of the layers consists of the lipid ratio I require (while
the other does not).
I will try your suggestion regarding to horizontal flip of that layer.
From: gmx-users-bounces at gromacs.org [mailto:gmx-users-bounces at gromacs.org]
On Behalf Of Justin A. Lemkul
Sent: Tuesday, June 02, 2009 23:49 PM
To: Discussion list for GROMACS users
Subject: Re: [gmx-users] Resizing heterogeneous membrane (POPG/POPE mixture)
Shay Amram wrote:
> Dear GROMACS users,
> I have been trying to expand a membrane by a non-integer multiplier (for
> example X1.5). This can be done using
> *genbox -cs Membrane.gro -o Expanded_membrane.gro -box 1.5X 1.5Y Z***
> This has worked very good with _homogenous_ membranes, and only very
> short equilibrium times were required to re-equilibrate the membrane
> Now I am trying to deal the same way with _heterogeneous_ membrane
> (POPG/POPE mixture, Mikko Karttunen et. al).
> When trying to invoke the above command to the POPG/POPE membrane I get
> a membrane with different compositions of lipids in the upper and lower
> This happens (so I suspect) because the arrangement of lipids in each
> leaflet is somewhat different so that the molecules that "get' to be
> multiplied do not conserve the
> same ratio of different lipids as in the original structure (which has
> ratio of 1:3 POPG/POPE).
If you have an excess of POPE, you can (either manually, or by using a
randomly choose some POPE to be converted to (or replaced with) POPG by some
clever means. This is not a trivial task.
> Is there a way to multiply the membrane by a _non-integer_ number AND
> somehow conserve the ratio of different lipid species too?
The output is dependent upon the configuration of the input and the box
no, there is no straightforward, automatic, way of doing so given the
constraints of your box dimensions.
> Or at the very least, to have both upper and lower leaflets identical
> after multiplication? (This, too, would help immensely).
If you can construct a satisfactory monolayer, you can manipulate its
coordinates (i.e. horizontal flip and translation) to create a bilayer. You
would have to do plenty of equilibration after doing so, of course.
> gmx-users mailing list gmx-users at gromacs.org
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Justin A. Lemkul
ICTAS Doctoral Scholar
Department of Biochemistry
jalemkul[at]vt.edu | (540) 231-9080
gmx-users mailing list gmx-users at gromacs.org
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-request at gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gromacs.org_gmx-users