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 ...
60 Phase behavior of coarse-grained lipid bilayers ht 6 ht7 ht8 ht9 Figure 5.2: Schematic drawing of the model lipids used in this study with their nomenclature. The black particles represent the head beads and the white particles the tail beads. The repulsion parameters used are aww = att = 25, a wh = 15, and awt = 80 (see equation 2.2, Chapter 2). In addition, we vary the head-head interaction parameter (a hh) to study the effect of changing the interactions between the hydrophilic segments of a lipid. In a real system the head-head interactions can be changed by, for example, adding salt to the system. All our simulations are performed on tensionless bilayers of 200 lipids. The total number of particles was 3500. The overall density of the system is ρ = 3. We initialize our system by distributing lipids randomly in water and we observe the self-assembly of a bilayer using DPD simulation only. After the bilayer is formed, we perform, in addition to the DPD moves, Monte Carlo moves in which we change the area as well. A typical simulation required 100,000 cycles of which 20,000 cycles were needed for equilibration. All the results are expressed in the usual reduced units, i.e. using Rc as the unit of length and kBTo = 1, with To room temperature, as unit of the energy. In the following we will denote as T ∗ the temperature expressed in this unit. 5.2.2 Results and Discussion In this section we first describe in detail the different phases of a bilayer formed by single tail lipids consisting of one head bead and nine tail beads (ht9), which we study at different reduced temperatures, from T ∗ = 0.8 to T ∗ = 1.5. At a fixed head-head repulsion of a hh = 35 an interdigitated gel phase is formed at low temperatures, while at a hh = 15 the non-interdigitated Lβ phase is formed. We then investigate the influence of changing the interactions between the headgroups and we show that we obtain the non-interdigitated Lβ phase or the interdigitated LβI phase dependent on the head-head repulsion parameter (see figure 5.3). Finally, we investigate the influence of tail length and we compare our results with experimental data on single-tail lipids.
5.2 Single-tail lipid bilayers 61 ahh 40 30 20 L EI L E 10 0.7 0.8 0.9 1.0 1.1 Figure 5.3: Computed phase diagram of the lipid ht9 as function of head-head repulsion parameter ahh and reduced temperature T ∗ . At high values of the head-head repulsion parameters the interdigitated LβI phase is formed, while at low values the non-interdigitated Lβ phase is formed. Increasing temperature causes the melting of the bilayer to the Lα phase. The lipid ht9 Head-head repulsion a hh = 35 In figure 5.4 the average area per lipid AL and the bilayer thickness, Dc, are plotted as function of temperature. The error bars have been calculated with the block averages method [76, 105]. In all the other plots of the area per lipid or bilayer thickness, we will not include error bars, which, however, have been estimated as ≤ 5%. A L 0.95 0.90 0.85 0.80 0.8 1.0 1.2 T* 1.4 1.6 (a) D c T* 8.20 7.80 7.40 L D 7.00 0.8 1.0 1.2 T* 1.4 1.6 Figure 5.4: Area per lipid AL and (b) bilayer thickness Dc as function of reduced temperature T ∗ for lipid type ht9. (b)
- Page 16 and 17: 10 Simulation method for coarse-gra
- Page 18 and 19: 12 Simulation method for coarse-gra
- Page 20 and 21: 14 Simulation method for coarse-gra
- Page 22 and 23: 16 Simulation method for coarse-gra
- Page 24 and 25: 18 Simulation method for coarse-gra
- Page 27 and 28: III Surface tension in lipid bilaye
- Page 29 and 30: 3.2 Method of calculation of surfac
- Page 31 and 32: 3.2 Method of calculation of surfac
- Page 33 and 34: 3.3 Constant surface tension ensemb
- Page 35 and 36: 3.3 Constant surface tension ensemb
- Page 37 and 38: 3.4 Surface tension in lipid bilaye
- Page 39 and 40: 3.4 Surface tension in lipid bilaye
- Page 41 and 42: 3.4 Surface tension in lipid bilaye
- Page 43 and 44: IV Structural characterization of l
- Page 45 and 46: 4.2 Structural quantities 39 been r
- Page 47 and 48: 4.3 Computational details 41 lipid
- Page 49 and 50: 4.4 Results and discussion 43 ρ(z)
- Page 51 and 52: 4.4 Results and discussion 45 one l
- Page 53 and 54: 4.4 Results and discussion 47 WH HT
- Page 55 and 56: 4.4 Results and discussion 49 Shill
- Page 57 and 58: 4.4 Results and discussion 51 chain
- Page 59 and 60: 4.4 Results and discussion 53 4.4.4
- Page 61: 4.4 Results and discussion 55 headg
- Page 64 and 65: 58 Phase behavior of coarse-grained
- Page 68 and 69: 62 Phase behavior of coarse-grained
- Page 70 and 71: 64 Phase behavior of coarse-grained
- Page 72 and 73: 66 Phase behavior of coarse-grained
- Page 74 and 75: 68 Phase behavior of coarse-grained
- Page 76 and 77: 70 Phase behavior of coarse-grained
- Page 78 and 79: 72 Phase behavior of coarse-grained
- Page 80 and 81: 74 Phase behavior of coarse-grained
- Page 82 and 83: 76 Phase behavior of coarse-grained
- Page 85 and 86: VI Interaction of small molecules w
- Page 87 and 88: 6.2 Computational details 81 For re
- Page 89 and 90: 6.3 Results and Discussion 83 withi
- Page 91 and 92: 6.3 Results and Discussion 85 ρ(z)
- Page 93 and 94: 6.3 Results and Discussion 87 ρ(z)
- Page 95 and 96: 6.3 Results and Discussion 89 S m 0
- Page 97 and 98: 6.3 Results and Discussion 91 π(z)
- Page 99: 6.3 Results and Discussion 93 the l
- Page 102 and 103: 96 Mesoscopic model for lipid bilay
- Page 104 and 105: 98 Mesoscopic model for lipid bilay
- Page 106 and 107: 100 Mesoscopic model for lipid bila
- Page 108 and 109: 102 Mesoscopic model for lipid bila
- Page 110 and 111: 104 Mesoscopic model for lipid bila
- Page 112 and 113: 106 Mesoscopic model for lipid bila
- Page 114 and 115: 108 Mesoscopic model for lipid bila
5.2 Single-tail <strong>lipid</strong> <strong>bilayers</strong> 61<br />
ahh<br />
40<br />
30<br />
20<br />
L EI<br />
L E<br />
10<br />
0.7 0.8 0.9 1.0 1.1<br />
Figure 5.3: Computed phase diagram <strong>of</strong> the <strong>lipid</strong> ht9 as function <strong>of</strong> head-head repulsion parameter<br />
ahh <strong>and</strong> reduced temperature T ∗ . At high values <strong>of</strong> the head-head repulsion parameters<br />
the interdigitated LβI phase is formed, while at low values the non-interdigitated Lβ phase<br />
is formed. Increasing temperature causes the melting <strong>of</strong> the bilayer to the Lα phase.<br />
The <strong>lipid</strong> ht9<br />
Head-head repulsion a hh = 35<br />
In figure 5.4 the average area per <strong>lipid</strong> AL <strong>and</strong> the bilayer thickness, Dc, are plotted as<br />
function <strong>of</strong> temperature. The error bars have been calculated <strong>with</strong> the block averages<br />
method [76, 105]. In all the other plots <strong>of</strong> the area per <strong>lipid</strong> or bilayer thickness, we<br />
will not include error bars, which, however, have been estimated as ≤ 5%.<br />
A L<br />
0.95<br />
0.90<br />
0.85<br />
0.80<br />
0.8 1.0 1.2<br />
T*<br />
1.4 1.6<br />
(a)<br />
D c<br />
T*<br />
8.20<br />
7.80<br />
7.40<br />
L D<br />
7.00<br />
0.8 1.0 1.2<br />
T*<br />
1.4 1.6<br />
Figure 5.4: Area per <strong>lipid</strong> AL <strong>and</strong> (b) bilayer thickness Dc as function <strong>of</strong> reduced temperature<br />
T ∗ for <strong>lipid</strong> type ht9.<br />
(b)