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

54 Structural characterization of lipid bilayers gion where the lipid headgroup are hydrated by the water, and is due to the fact that double-tail lipids have a larger headgroup (three beads) respect to the one of singletail lipids (one bead). Also, at the HT interface the total density shows more marked minima in the case of double-tail lipids, and a small difference between the two bilayers is observed in the bilayer hydrophobic core where the total density of the tail beads is slightly larger for the single tail lipids. A remarkable fact is that, despite the π(z) 0.1 0.0 −0.1 −0.2 −0.3 −6 −4 −2 0 2 4 6 z (a) (b) Figure 4.11: Figure (a) shows the lateral pressure profile distribution in a bilayer of singletail lipids (gray line) and a bilayer of double-tail lipids (black) line. In figure (b) the pressure profile from atomistic MD simulations of a DPPC bilayer is shown. Figure (b) is reproduced from reference [77] with kind permission of the authors. differences in the density profiles of the two bilayers being not too large, the distribution of lateral pressure in the double-tail lipid bilayer is very different compared with the single-tail lipid bilayer, as shown in figure 4.11(a). In the case of double-tail lipids, the peaks at the WH and HT interface become more pronounced, indicating that the lipids are more densely packed. Also, because of the larger headgroup size of the double-tail lipids, the peak at the WH interface broadens. The central peak (at the bilayer midplane) decreases in size, due to the lower density in the bilayer center. Most noticeably, we observe the appearance of two secondary peaks at the edge of the hydrophobic core, which are not present in the case of the single-tail lipids. These peaks correspond to the region of maximum density for the tail beads (excluded the terminal one). This shape of the pressure profile is remarkably similar to the one calculated by Lindahl and Edholm in atomistic MD simulations of a dipalmitoylphosphatidylcholine bilayer [77], and shown in figure 4.11(b). In Chapter 2 we have shown that, if a DPD bead is chosen to represent three methyl groups, the lipid h3(t5)2 can be mapped onto the DMPC phospholipid. The structure of DPPC is very similar to the one of DMPC, both having a phosphocholine

4.4 Results and discussion 55 headgroup and two saturated hydrocarbon chains, with DPPC having 16 carbons in each tail, and DMPC 14. This result indicates that the coarse-grained model we have developed for double-tail lipids well reproduces the distribution of pressure in a bilayer of saturated phosphocholines. However, care must be used in such a conclusion. We want to point out that, given the differences between the atomistic and the CG representations, and the different way in which the inter-and intra-molecular interactions are implemented in the two models, the remarkable similarity between the pressure profiles in the two model might be fortuitous. To further address this issue, it would be interesting to compare in a more systematic way pressure profiles calculated from MD and CG simulations. Unfortunately, the published pressure profiles in lipid bilayers from MD simulations are still very few.

54 Structural characterization <strong>of</strong> <strong>lipid</strong> <strong>bilayers</strong><br />

gion where the <strong>lipid</strong> headgroup are hydrated by the water, <strong>and</strong> is due to the fact that<br />

double-tail <strong>lipid</strong>s have a larger headgroup (three beads) respect to the one <strong>of</strong> singletail<br />

<strong>lipid</strong>s (one bead). Also, at the HT interface the total density shows more marked<br />

minima in the case <strong>of</strong> double-tail <strong>lipid</strong>s, <strong>and</strong> a small difference between the two <strong>bilayers</strong><br />

is observed in the bilayer hydrophobic core where the total density <strong>of</strong> the tail<br />

beads is slightly larger for the single tail <strong>lipid</strong>s. A remarkable fact is that, despite the<br />

π(z)<br />

0.1<br />

0.0<br />

−0.1<br />

−0.2<br />

−0.3<br />

−6 −4 −2 0 2 4 6<br />

z<br />

(a) (b)<br />

Figure 4.11: Figure (a) shows the lateral pressure pr<strong>of</strong>ile distribution in a bilayer <strong>of</strong> singletail<br />

<strong>lipid</strong>s (gray line) <strong>and</strong> a bilayer <strong>of</strong> double-tail <strong>lipid</strong>s (black) line. In figure (b) the pressure<br />

pr<strong>of</strong>ile from atomistic MD simulations <strong>of</strong> a DPPC bilayer is shown. Figure (b) is reproduced<br />

from reference [77] <strong>with</strong> kind permission <strong>of</strong> the authors.<br />

differences in the density pr<strong>of</strong>iles <strong>of</strong> the two <strong>bilayers</strong> being not too large, the distribution<br />

<strong>of</strong> lateral pressure in the double-tail <strong>lipid</strong> bilayer is very different compared<br />

<strong>with</strong> the single-tail <strong>lipid</strong> bilayer, as shown in figure 4.11(a). In the case <strong>of</strong> double-tail<br />

<strong>lipid</strong>s, the peaks at the WH <strong>and</strong> HT interface become more pronounced, indicating<br />

that the <strong>lipid</strong>s are more densely packed. Also, because <strong>of</strong> the larger headgroup size<br />

<strong>of</strong> the double-tail <strong>lipid</strong>s, the peak at the WH interface broadens. The central peak<br />

(at the bilayer midplane) decreases in size, due to the lower density in the bilayer<br />

center. Most noticeably, we observe the appearance <strong>of</strong> two secondary peaks at the<br />

edge <strong>of</strong> the hydrophobic core, which are not present in the case <strong>of</strong> the single-tail<br />

<strong>lipid</strong>s. These peaks correspond to the region <strong>of</strong> maximum density for the tail beads<br />

(excluded the terminal one). This shape <strong>of</strong> the pressure pr<strong>of</strong>ile is remarkably similar<br />

to the one calculated by Lindahl <strong>and</strong> Edholm in atomistic MD simulations <strong>of</strong> a<br />

dipalmitoylphosphatidylcholine bilayer [77], <strong>and</strong> shown in figure 4.11(b).<br />

In Chapter 2 we have shown that, if a DPD bead is chosen to represent three<br />

methyl groups, the <strong>lipid</strong> h3(t5)2 can be mapped onto the DMPC phospho<strong>lipid</strong>. The<br />

structure <strong>of</strong> DPPC is very similar to the one <strong>of</strong> DMPC, both having a phosphocholine

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

Saved successfully!

Ooh no, something went wrong!