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 ...

science.uva.nl
from science.uva.nl More from this publisher
30.08.2013 Views

72 Phase behavior of coarse-grained lipid bilayers As we increase the tail length the gel phases are stabilized and the transition shifts to higher temperatures. The effect of increasing the head-head repulsion on the gel to liquid crystalline transition temperature is much more pronounced for the Lβ → Lα compared to LβI → Lα. This can be understood from the fact that in the interdigitated phase the average distance between the heads is already much larger compared to the non-interdigitated phase, and a further increase in this distance does not have a dramatic effect on the stability of the gel phase. For lipids ht8 and ht9 the LβI phase occurs at slightly lower repulsion parameters than for lipids ht6 and ht7. This is consistent with experimental results [91]. Since the interdigitated phase is more closely packed than the non-interdigitated phase, the van der Waals energy is greater. This energy gain is proportional to the number of carbon atoms in the phospholipid chain and thus interdigitation becomes energetically more favorable for longer chains. Also in our simulations we observe that the interdigitated phase is more compact and hence a∗ hh decreases slightly with increasing tail length. It is interesting to compare these results with the experimental data. Misquitta and Caffrey in [108] systematically investigate the phase diagrams of monoacylglycerols, a single-tail lipid, and show a similar tail length dependence for the Lβ →Lα transition. Interestingly, as we will show later, for a similar model of a double-tail lipid we do not observe the spontaneous formation of an interdigitated phase. This corresponds to the experimental observation that for the most common double-tail lipids the interdigitated phase does not form spontaneously, but should be induced by the addition of, for example, alcohol [98]. The effect of adding salt on the gel to liquid crystalline transition has been studied for double-tail lipids [109] and recently for single-tail lipids [110]. These studies show that adding so-called kosmotropic salts increases the Lβ →Lα transition temperature, while chaotropic salts decrease this transition temperature. Similar effects have been observed for nonionic single-tail lipids [111]. Takahashi et al. [110] explain these observations by assuming that kosmotropes tend to be excluded from the interfacial region and hence reduce the amount of interfacial water, while chaotropic salts have the inverse effects, i.e. are adsorbed at the interfacial region and increase the amount of interfacial water. In our model a similar effect can be achieved by changing the head-head interactions; increasing or decreasing a hh corresponds to adding chaotropes or kosmotropes, respectively. Our simulations show that decreasing the head-head repulsion stabilizes the Lβ phase, which corresponds to the case that water is excluded from the interface. Adding chaotropic salts has the reverse effect: it increases the head-head repulsion and stabilizes the Lα phase. Our simulations show that at sufficiently high head-head repulsion the interdigitated phase (LβI) is stable. This suggests that it might be possible in experiments to induce the Lβ →LβI phase transition by adding chaotropic salts to the systems.

5.3 Double-tail lipid bilayers 73 5.3 Double-tail lipid bilayers In the previous section we have discussed the phase behavior of single-tail lipid bilayers. In this section we extend the model and investigate the phase behavior of double-tail lipids. In particular, we consider a coarse-grained lipid with three headbeads and two tails of five beads each. This lipid will be denoted as h3(t5)2. As discussed in section 2.2 of Chapter 2, if a DPD bead is taken to represent the volume of three water molecules, i.e. 90 ˚A 3 , then the coarse-grained lipid h3(t5)2 can be considered as a model for dimyristoylphosphatidylcholine (DMPC) (see figure 2.2 in Chapter 2). 5.3.1 Computational details For this study we consider bilayers of 900 h3(t5)2 lipids, and 25 water particles per lipid, corresponding to fully hydrated conditions (i.e. no interaction of the bilayer with its periodic images). The values of the interaction parameters are aww = att = 25, a hh = 35, awt = a ht = 80. Two consecutive beads in the lipid chain are connected by a harmonic spring, with equilibrium distance ro = 0.7, and elastic constant Kr = 100. Consecutive bonds in the lipid chain are subjected to a harmonic bond-bending potential. The values of the parameters related to this bond-bending potential were derived from the comparison of the CG model with MD simulations on all-atom model for a DMPC lipid bilayer [112]. The resulting values for the bending constant and the equilibrium angle in the lipid tails are Kθ=6 and θo=180 o , respectively. About the bond-bending potential between the head-bead connected to the lipid tails and the first beads in the tails values of Kθ=3 and θo=90 o were found to reproduce the correct configurational distribution and structure of the atomistic detailed phospholipid. All the simulations were carried at zero surface tension conditions. To explore the phase diagram of the bilayer, the temperature of the system was gradually decreased from T ∗ = 1.0 to T ∗ = 0.2 . At each temperature, a total of 100,000 DPD-MC cycles was performed, of which the first 20,000 cycles were needed for equilibration. Statistical averages were then collected over the next 80,000 DPD-MC cycles. The phase boundaries were detected, as described in the previous sections for the single-tail lipids, by monitoring the temperature behavior of the area per lipid, the bilayer thickness and the order of the tails. 5.3.2 Results and Discussion The area per lipid, AL, bilayer hydrophobic thickness, Dc, and tail order parameter S tail, as function of reduced temperature, T ∗ , are shown in Figure 5.15.

72 Phase behavior <strong>of</strong> coarse-grained <strong>lipid</strong> <strong>bilayers</strong><br />

As we increase the tail length the gel phases are stabilized <strong>and</strong> the transition shifts<br />

to higher temperatures. The effect <strong>of</strong> increasing the head-head repulsion on the gel to<br />

liquid crystalline transition temperature is much more pronounced for the Lβ → Lα<br />

compared to LβI → Lα. This can be understood from the fact that in the interdigitated<br />

phase the average distance between the heads is already much larger compared<br />

to the non-interdigitated phase, <strong>and</strong> a further increase in this distance does not have<br />

a dramatic effect on the stability <strong>of</strong> the gel phase.<br />

For <strong>lipid</strong>s ht8 <strong>and</strong> ht9 the LβI phase occurs at slightly lower repulsion parameters<br />

than for <strong>lipid</strong>s ht6 <strong>and</strong> ht7. This is consistent <strong>with</strong> experimental results [91]. Since<br />

the interdigitated phase is more closely packed than the non-interdigitated phase,<br />

the van der Waals energy is greater. This energy gain is proportional to the number<br />

<strong>of</strong> carbon atoms in the phospho<strong>lipid</strong> chain <strong>and</strong> thus interdigitation becomes energetically<br />

more favorable for longer chains. Also in our simulations we observe that<br />

the interdigitated phase is more compact <strong>and</strong> hence a∗ hh decreases slightly <strong>with</strong> increasing<br />

tail length.<br />

It is interesting to compare these results <strong>with</strong> the experimental data. Misquitta<br />

<strong>and</strong> Caffrey in [108] systematically investigate the phase diagrams <strong>of</strong> monoacylglycerols,<br />

a single-tail <strong>lipid</strong>, <strong>and</strong> show a similar tail length dependence for the Lβ →Lα<br />

transition. Interestingly, as we will show later, for a similar model <strong>of</strong> a double-tail<br />

<strong>lipid</strong> we do not observe the spontaneous formation <strong>of</strong> an interdigitated phase. This<br />

corresponds to the experimental observation that for the most common double-tail<br />

<strong>lipid</strong>s the interdigitated phase does not form spontaneously, but should be induced<br />

by the addition <strong>of</strong>, for example, alcohol [98].<br />

The effect <strong>of</strong> adding salt on the gel to liquid crystalline transition has been studied<br />

for double-tail <strong>lipid</strong>s [109] <strong>and</strong> recently for single-tail <strong>lipid</strong>s [110]. These studies<br />

show that adding so-called kosmotropic salts increases the Lβ →Lα transition temperature,<br />

while chaotropic salts decrease this transition temperature. Similar effects<br />

have been observed for nonionic single-tail <strong>lipid</strong>s [111]. Takahashi et al. [110] explain<br />

these observations by assuming that kosmotropes tend to be excluded from the interfacial<br />

region <strong>and</strong> hence reduce the amount <strong>of</strong> interfacial water, while chaotropic<br />

salts have the inverse effects, i.e. are adsorbed at the interfacial region <strong>and</strong> increase<br />

the amount <strong>of</strong> interfacial water. In our model a similar effect can be achieved by<br />

changing the head-head interactions; increasing or decreasing a hh corresponds to<br />

adding chaotropes or kosmotropes, respectively. Our simulations show that decreasing<br />

the head-head repulsion stabilizes the Lβ phase, which corresponds to the case<br />

that water is excluded from the interface. Adding chaotropic salts has the reverse<br />

effect: it increases the head-head repulsion <strong>and</strong> stabilizes the Lα phase. Our simulations<br />

show that at sufficiently high head-head repulsion the interdigitated phase<br />

(LβI) is stable. This suggests that it might be possible in experiments to induce the<br />

Lβ →LβI phase transition by adding chaotropic salts to the systems.

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

Saved successfully!

Ooh no, something went wrong!