Re: Lipid bilayer tears apart during simulation

From: Aron Broom (broomsday_at_gmail.com)
Date: Fri Jun 14 2013 - 15:34:31 CDT

I haven't done much with lipids, but I when I did they were pretty stable.
Two things to consider:

1) Is there a particular reason to use NVT versus NPT? Your bilayer may be
tearing apart because of incorrect pressures. On that note, which
forcefield are you using? Some lipid forcefields need you to use a
constant surface tension adjustment in order to be stable. I think the
latest CHARMM ones are OK, but certainly the older ones, and AMBER ones
still need that.

2) Have you neutralized the charge in the system and did you give the ions
time to equilibrate to their preferred positions near the headgroups?

~Aron

On Fri, Jun 14, 2013 at 4:14 PM, Peterson J <peterson.jjjj_at_gmail.com> wrote:

> Hi all,
>
> I am currently simulating a membrane protein homology model embedded in
> lipid bilayer and performed simulated annealing and lipid melting steps
> successfully. But when the system is subjected under equilibration
> simulation (NVT) with constraints only on protein Ca atoms, I see the lipid
> bilayer tearing apart.
> What are the possible solution to solve this?
>
> Thanks,
>
> Peterson
>

-- 
Aron Broom M.Sc
PhD Student
Department of Chemistry
University of Waterloo

This archive was generated by hypermail 2.1.6 : Wed Dec 31 2014 - 23:21:18 CST