30.08.2013 Views

Mesoscopic models of lipid bilayers and bilayers with embedded ...

Mesoscopic models of lipid bilayers and bilayers with embedded ...

Mesoscopic models of lipid bilayers and bilayers with embedded ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

5.3 Double-tail <strong>lipid</strong> <strong>bilayers</strong> 73<br />

5.3 Double-tail <strong>lipid</strong> <strong>bilayers</strong><br />

In the previous section we have discussed the phase behavior <strong>of</strong> single-tail <strong>lipid</strong> <strong>bilayers</strong>.<br />

In this section we extend the model <strong>and</strong> investigate the phase behavior <strong>of</strong><br />

double-tail <strong>lipid</strong>s. In particular, we consider a coarse-grained <strong>lipid</strong> <strong>with</strong> three headbeads<br />

<strong>and</strong> two tails <strong>of</strong> five beads each. This <strong>lipid</strong> will be denoted as h3(t5)2. As<br />

discussed in section 2.2 <strong>of</strong> Chapter 2, if a DPD bead is taken to represent the volume<br />

<strong>of</strong> three water molecules, i.e. 90 ˚A 3 , then the coarse-grained <strong>lipid</strong> h3(t5)2 can be<br />

considered as a model for dimyristoylphosphatidylcholine (DMPC) (see figure 2.2 in<br />

Chapter 2).<br />

5.3.1 Computational details<br />

For this study we consider <strong>bilayers</strong> <strong>of</strong> 900 h3(t5)2 <strong>lipid</strong>s, <strong>and</strong> 25 water particles per<br />

<strong>lipid</strong>, corresponding to fully hydrated conditions (i.e. no interaction <strong>of</strong> the bilayer<br />

<strong>with</strong> its periodic images).<br />

The values <strong>of</strong> the interaction parameters are aww = att = 25, a hh = 35, awt = a ht =<br />

80. Two consecutive beads in the <strong>lipid</strong> chain are connected by a harmonic spring,<br />

<strong>with</strong> equilibrium distance ro = 0.7, <strong>and</strong> elastic constant Kr = 100.<br />

Consecutive bonds in the <strong>lipid</strong> chain are subjected to a harmonic bond-bending<br />

potential. The values <strong>of</strong> the parameters related to this bond-bending potential were<br />

derived from the comparison <strong>of</strong> the CG model <strong>with</strong> MD simulations on all-atom<br />

model for a DMPC <strong>lipid</strong> bilayer [112]. The resulting values for the bending constant<br />

<strong>and</strong> the equilibrium angle in the <strong>lipid</strong> tails are Kθ=6 <strong>and</strong> θo=180 o , respectively. About<br />

the bond-bending potential between the head-bead connected to the <strong>lipid</strong> tails <strong>and</strong><br />

the first beads in the tails values <strong>of</strong> Kθ=3 <strong>and</strong> θo=90 o were found to reproduce the<br />

correct configurational distribution <strong>and</strong> structure <strong>of</strong> the atomistic detailed phospho<strong>lipid</strong>.<br />

All the simulations were carried at zero surface tension conditions. To explore the<br />

phase diagram <strong>of</strong> the bilayer, the temperature <strong>of</strong> the system was gradually decreased<br />

from T ∗ = 1.0 to T ∗ = 0.2 . At each temperature, a total <strong>of</strong> 100,000 DPD-MC cycles was<br />

performed, <strong>of</strong> which the first 20,000 cycles were needed for equilibration. Statistical<br />

averages were then collected over the next 80,000 DPD-MC cycles.<br />

The phase boundaries were detected, as described in the previous sections for<br />

the single-tail <strong>lipid</strong>s, by monitoring the temperature behavior <strong>of</strong> the area per <strong>lipid</strong>,<br />

the bilayer thickness <strong>and</strong> the order <strong>of</strong> the tails.<br />

5.3.2 Results <strong>and</strong> Discussion<br />

The area per <strong>lipid</strong>, AL, bilayer hydrophobic thickness, Dc, <strong>and</strong> tail order parameter<br />

S tail, as function <strong>of</strong> reduced temperature, T ∗ , are shown in Figure 5.15.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!