J. Chem. Phys. 135, 244701 (2011); http://dx.doi.org/10.1063/1.3660673 (22 pages)
Thermal fluctuations in shape, thickness, and molecular orientation in lipid bilayers
(Received 19 July 2011; accepted 26 October 2011; published online 22 December 2011)
© 2011 American Institute of Physics
Article Outline
- INTRODUCTION
- BILAYER ENERGETICS
- GEOMETRIC DESCRIPTION
- Macroscopic bilayer free energy
- Protrusions
- FLUCTUATION SPECTRA OF HOMOGENEOUS MEMBRANES
- Decoupled protrusion/bending approximation
- γλ = 0 approximation
- SIMULATION DETAILS
- Implicit solvent model
- MARTINI force field model (DPPC)
- COMPARISON TO THEORY
- CONCLUSION
RELATED DATABASES
KEYWORDS and PACS
ARTICLE DATA
- E. Browicz, Zbl. Med. Wiss. 28, 625 (1890).
- F. Brochard and J. F. Lennon, J. Phys. (France) 36, 1035 (1975). [Inspec] [ISI]
- J. Israelachvili and H. Wennerström, J. Phys. Chem. 96, 520 (1992).
- R. Lipowsky, Nature (London) 349, 475 (1991). [MEDLINE]
- C. S. Peskin, G. M. Odell, and G. F. Oster, Biophys. J. 65, 316 (1993). [MEDLINE] [CAS]
- R. Bruinsma and P. Pincus, Curr. Opin. Solid State Mater. Sci. 1, 401 (1996). [Inspec] [ISI]
- R. Golestanian, M. Goulian, and M. Kardar, Phys. Rev. E 54, 6725 (1996). [MEDLINE] [CAS]
- W. Helfrich, Z. Naturforsch. C 28, 693 (1973). [MEDLINE]
- P. B. Canham, J. Theor. Biol. 26, 61 (1970). [MEDLINE]
- F. C. Frank, Discuss. Faraday Soc. 25, 19 (1958).
- R. Lipowsky and S. Grotehans, Biophys. Chem. 49, 27 (1994). [ISI]
- J. N. Israelachvili, Intermolecular and Surface Forces (Academic, New York, 1985).
- F. C. MacKintosh and T. C. Lubensky, Phys. Rev. Lett. 67, 1169 (1991). [MEDLINE]
- M. Hamm and M. M. Kozlov, Eur. J. Phys. E. 3, 323 (2000). [ISI]
- Y. Kozlovsky and M. M. Kozlov, Biophys. J. 82, 882 (2002). [ISI] [MEDLINE] [CAS]
- P. I. Kuzmin, S. A. Akimov, Y. A. Chizmadzhev, J. Zimmerberg, and F. S. Cohen, Biophys. J. 88, 1120 (2005). [Inspec] [MEDLINE]
- A. Kessel, N. Ben-Tal, and S. May, Biophys. J. 81, 643 (2001). [Inspec] [ISI] [MEDLINE] [CAS]
- G. Khelashvili, G. Pabst, and D. Harries, J. Phys. Chem. B 114, 7524 (2010).
- J. Aittoniemi, T. Róg, P. Niemelä, M. Pasenkiewicz-Gierula, M. Karttunen, and I. Vattulainen, J. Phys. Chem. B 110, 25562 (2006).
- S. May and A. Ben-Shaul, Phys. Chem. Chem. Phys. 2, 4494 (2000).
- S. May and A. Ben-Shaul, Biophys. J. 76, 751 (1999). [Inspec] [ISI] [MEDLINE] [CAS]
- S. May, Eur. Biophys. J. 29, 17 (2000). [Inspec] [MEDLINE]
- S. May, Langmuir 18, 6356 (2002). [ISI] [CAS]
- M. Fosňarič, A. Iglič, and S. May, Phys. Rev. E 74, 051503 (2006).
- K. Bohinc, V. Kralj-Iglič, and S. May, J. Chem. Phys. 119, 7433 (2000)JCPSA6000119000014007435000001.
- Y. Kozlovsky, J. Zimmerberg, and M. M. Kozlov, Biophys. J. 87, 999 (2004). [MEDLINE]
- J.-B. Fournier, Eur. J. Phys. E 11, 261 (1999).
- J.-B. Fournier, Europhys. Lett. 43, 725 (1998). [Inspec] [ISI]
- B. West, F. L. H. Brown, and F. Schmid, Biophys. J. 96, 101 (2009). [Inspec] [MEDLINE] [CAS]
- M. Hamm and M. M. Kozlov, Eur. J. Phys. E 6, 519 (1998). [ISI]
- S. Tristram-Nagle, R. Zhang, R. M. Suter, C. R. Worthington, W. J. Sun, and J. F. Nagle, Biophys. J. 64, 1097 (1993).
- T. C. Lubensky and F. C. MacKintosh, Phys. Rev. Lett. 71, 1565 (1993). [MEDLINE]
- C.-M. Chen, T. C. Lubensky, and F. C. MacKintosh, Phys. Rev. E 51, 504 (1995). [MEDLINE]
- K. Sengupta, V. A. Raghunathan, and Y. Hatwalne, Phys. Rev. Lett. 87, 055705 (2001). [ISI] [MEDLINE] [CAS]
- U. Seifert, J. Shillcock, and P. Nelson, Phys. Rev. Lett. 77, 5237 (1996). [MEDLINE]
- E. A. G. Aniansson, S. N. Wall, M. Almgren, H. Hoffmann, I. Kielmann, W. Ulbricht, R. Zana, J. Lang, and C. Tondre, J. Phys. Chem. 80, 905 (1976). [ISI] [CAS]
- T. J. McIntosh, S. Advani, R. E. Burton, D. V. Zhelev, D. Needham, and S. A. Simon, Biochemistry 34, 8520 (1995).
- A. Halperin and O. Mouritsen, Eur. Biophys. J. 34, 967 (2005).
- R. Goetz, G. Gompper, and R. Lipowsky, Phys. Rev. Lett. 82, 221 (1999).
- E. Lindahl and O. Edholm, J. Chem. Phys. 113, 3882 (2000)JCPSA6000113000009003882000001.
- E. Lindahl and O. Edholm, Biophys. J. 79, 426 (2000). [Inspec] [ISI] [MEDLINE] [CAS]
- S. J. Marrink and A. E. Mark, J. Phys. Chem. 105, 6122 (2001). [CAS]
- S. W. Chiu, S. Vasudevan, E. Jakobsson, R. J. Mashl, and H. L. Scott, Biophys. J 85, 3624 (2003).
- G. Brannigan and F. Brown, Biophys. J. 90, 1501 (2006). [Inspec] [MEDLINE] [CAS]
- G. Brannigan, P. F. Philips, and F. L. H. Brown, Phys. Rev. E 72, 011915 (2005). [ISI]
- G. Brannigan, A. C. Tamboli, and F. L. H. Brown, J. Chem. Phys. 121 (2004)JCPSA6000121000007003259000001. [ISI] [MEDLINE]
- M. J. Stevens, J. Chem. Phys. 121, 11942 (2004)JCPSA6000121000023011942000001. [ISI] [MEDLINE]
- I. R. Cooke, K. Kremer, and M. Deserno, Phys. Rev. E 72, 011506 (2005).
- W. K. den Otter, J. Chem. Phys. 123, 214906 (2005)JCPSA6000123000021214906000001. [MEDLINE]
- J. Neder, B. West, P. Nielaba, and F. Schmid, J. Chem. Phys. 132, 115101 (2009)JCPSA6000132000011115101000001. [MEDLINE]
- B. West and F. Schmid, Soft Matter 6, 1275 (2010).
- E. R. May, A. Narang, and D. I. Kopelevich, Phys. Rev. E 76, 021913 (2007). [ISI]
- E. R. May, A. Narang, and D. I. Kopelevich, Mol. Simul. 33, 787 (2007).
- F. M. Thakkar, P. K. Maiti, V. Kumaran, and K. G. Ayappa, Soft Matter 7, 3963 (2011).
- S. J. Marrink, J. H. Risselada, S. Yefimov, P. D. Tieleman, and A. H. de Vries, J. Phys. Chem. B 111, 7812 (2007). [Inspec] [MEDLINE]
- S. J. Marrink, A. H. de Vries, and A. E. Mark, J. Phys. Chem. B 108, 750 (2004).
- I. W. Stewart, The Static and Dynamic Continuum Theory of Liquid Crystals: A Mathematical Introduction (Taylor & Francis, London, 2004).
- A. Imparato, J. Chem. Phys. 124, 154714 (2006)JCPSA6000124000015154714000001. [MEDLINE]
- J.-B. Fournier and C. Barbetta, Phys. Rev. Lett. 100, 078103 (2008).
- W. Cai, T. C. Lubensky, P. Nelson, and T. Powers, J. Phys. II 4, 931 (1994). [ISI] [CAS]
- O. Farago and P. Pincus, Eur. Phys. J. E 11, 399 (2003). [MEDLINE]
- O. Farago and P. Pincus, J. Chem. Phys. 120, 2934 (2004)JCPSA6000120000006002934000001. [MEDLINE]
- J. Stecki, J. Chem. Phys. 125, 154902 (2006)JCPSA6000125000015154902000001. [MEDLINE]
- F. Schmid, Eur. Phys. Lett. 95, 28008 (2011).
- S. A. Safran, Statistical Thermodynamics of Surfaces, Interfaces and Membranes (Westview Press, Boulder, CO, 1994).
- N. G. V. Kampen, Stochastic Processes in Physics and Chemistry (Amsterdam, North Holland, 2007).
- S. J. Plimpton, J. Comput. Phys. 117, 1 (1995).
- R. Goetz and R. Lipowsky, J. Chem. Phys. 108, 7397 (1998).
- P. J. Steinbach and B. B. Brooks, J. Comput. Chem. 15, 667 (1994).
- O. Farago, J. Chem. Phys. 119, 596 (2003).
- D. Marsh, Biochim. Biophys. Acta. 1286, 183 (1996).
- E. Lindahl and O. Edholm, J. Chem. Phys. 113, 3882 (2000)JCPSA6000113000009003882000001.
- M. Bhandarkar et al., NAMD User's Guide. Version 2.7b2, University of Illinois and Beckman Institute, Urbana, IL (2009), sec. 12.2.
- A. Harasima, Adv. Chem. Phys. 1, 203 (1958).
- D. van der Spoel, E. Lindahl, B. Hess, G. Groenhof, A. E. Mark, and H. J. C. Berendsen, J. Comput. Chem. 26, 1701 (2005). [Inspec] [MEDLINE] [CAS]
- H. J. C. Berendsen, D. van der Spoel, and R. van Drunen, Comput. Phys. Commun. 91, 43 (1995).
- B. Hess, C. Kutzner, D. van der Spoel, and E. Lindahl, J. Chem. Theory Comput. 4, 435 (2008).
- See http://md.chem.rug.nl/cgmartini to download pre-assembled bilayer coordinates and force field files.
- S. O. H. Ollila, J. H. Risselada, M. Louhivuori, E. Lindahl, I. Vattulainen, and S. J. Marrink, Phys. Rev. Lett. 102, 078101 (2009). [MEDLINE]
- J. H. Irving and J. G. Kirkwood, J. Chem. Phys. 18, 817 (1950)JCPSA6000018000006000817000001.
- J. Sonne, F. Y. Hansen, and G. H. Peters, J. Chem. Phys. 122, 124903 (2005)JCPSA6000122000012124903000001. [ISI] [MEDLINE]
- W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes: The Art of Scientific Computing, 3rd ed. (Cambridge University Press, Cambridge, England, 2007).
- S. E. Feller and R. W. Pastor, J. Chem. Phys. 111, 1281 (1999)JCPSA6000111000003001281000001.
- S. Shkulipa, Computer Simulations of Lipid Bilayer Dynamics, Ph.D. dissertation, University of Twente, The Netherlands, 2006).
- Q. Waheed and O. Edholm, Biophys. J. 97, 2754 (2009). [Inspec] [MEDLINE] [CAS]
- J. S. Rowlinson and B. Widom, Molecular Theory of Capillarity (Dover, New York, 2003).
- G. Brannigan and F. L. H. Brown, Biophys. J. 92, 864 (2007). [MEDLINE]
- M. Watson and F. L. H. Brown, Biophys. J 98, L9 (2010). [MEDLINE] [CAS]
- E. G. Brandt, A. R. Braun, J. N. Sachs, J. F. Nagle, and O. Edholm, Biophys. J. 100, 2104 (2011).
- S. May, Y. Kozlovsky, A. Ben-Shaul, and M. M. Kozlov, Eur. Phys. J. E 14, 299 (2004). [MEDLINE]
- M. Frigo and S. G. Johnson, Proccedings of IEEE 93, 216 (2005).
Figures (click on thumbnails to view enlargements)
(α) are normal to z(α) and point toward the interior of the bilayer. The unit vectors field
(α) points along the hydrocarbon chains. b(α)
(α) extend from z(α) to the surface separating the two leaflets, z(m). In other words, the top monolayer is bounded by z(1) and z(m), while the bottom monolayer is bounded by z(m) and z(2). The mean height z+ is the average of z(1) and z(2). Left: a bilayer in its minimal energy configuration, in which λ(α) = 0, z(m) = z+, N(α) = n(α) = 0, the thickness is 2b0 and the area per molecule is Σ0 (dotted red). The volume per lipid v satisfies v = b0Σ0. Since there are no protrusions, h(α) = z(α). Right: an arbitrarily deformed bilayer. On short length scales, the polar-nonpolar interfaces h(α) are not smooth (dashed curves). The protrusion fields λ(α) displace the interface in the normal direction, so that (−1)αλ(α)
(α) (dashed vectors) extend from z(α) to h(α). The fields z(m) (black) and z+ (blue) differ in general. The thickness deformations are exaggerated for illustrative purposes. We will assume throughout the paper that the absolute values of the following quantities are much less than one: Nj(α), nj(α), ∇N(α), ∇n(α), ∇z(m), Σ(α)/Σ0−1, b(α)/b0−1, |z(α)−z(m)|/b0−1, (z+−z(m))/b0,∇λ(α),λ(α)/b0.
FIG.1 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
(1) points from the displaced surface toward the end of the hydrocarbon chain (along the red dashed vector). The dashed vector only shows the direction of
(1) but not its magnitude, since
(1) is a unit vector. The monolayer thickness is large compared to λ(1), so
(1) ≈
(1). Bottom:
(1) ≠
(1) due to a gradient in the protrusion field.
FIG.2 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
(α),21 , 22 which is a function of the hydrocarbon chain length b(α), the cross sectional area at the polar-nonpolar interface Σ(α) (green), and the cross sectional area at the center of the head group Σh(α) (yellow). The two surfaces which contain Σ(α) and Σh(α) are separated by a fixed distance ℓh. The volume v of the hydrocarbon chain region (finely dashed box) is constant to enforce chain incompressibility.
FIG.3 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
FIG.4 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
FIG.5 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
FIG.6 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
(α) is measured within each model. Note that for visual purposes
(α) is not normalized. Left: Our implicit solvent model (CG). Each lipid consists of a hydrophilic head bead (green), an interfacial bead (orange), and three hydrophobic tail beads (gray). Though the potentials were designed to mimic generic intermolecular forces, no explicit electrostatic interactions are present.
(α) points from the interfacial bead toward the last tail bead. Right: For the MARTINI force field (DPPC), each lipid consists of a positively charged bead representing the choline group (red), a negatively charged bead representing the phosphate group (purple), two beads of intermediate hydrophobicity representing the glycerol ester linkage (green), and two chains of four hydrophobic beads each, representing the hydrocarbons (brown). The charged beads interact via a shifted Coulombic potential energy function.55
(α) points from the midpoint of the interfacial beads toward the midpoint of the last tail beads.
FIG.7 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
FIG.8 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
q⊥|2〉, 〈|
q⊥|2〉, 〈|
q∥|2〉, 〈|
q∥|2〉 fluctuations for the CG model. Simulation data are displayed as circles. The data points used for the fits (black) are evenly spaced along the qx and qy axes. Data points which do not lie on the axes are shown in orange. The solid curves represent best fits of the data to Eqs. ( 39 , 40 , 41 , 42 , 43 ) with γλ = 0. Fit parameters are listed in Table 3. 〈|hq|2〉 is also plotted on a semi-log scale in the inset. Dashed curves correspond to the protrusion-free limit by using Eqs. ( 39 , 40 , 41 , 42 ) and the same values in Table 3 but with kλ → ∞. While only a small contribution to the height and thickness modes, protrusions are essential for explaining the qualitative behavior of 〈|
q∥|2〉 and 〈|
q∥|2〉.
FIG.9 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
FIG.10 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
−∥q〉, 〈tq
−∥q〉 and 〈|
q(α)|2〉 agree with the theoretical predictions very closely. Both the theory and data show that the fluctuations in the cross-terms are purely imaginary. This may be understood mathematically since
q∥* = −
−q∥ and hq* = h−q and similarly for the peristaltic term (see Sec. 3). The monolayer tilt averaged over all directions 〈|
q(α)|2〉 is also plotted in MNK.
FIG.11 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
1,
2,
3, and
4.
FIG.12 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
5,
6,
7,
8) are located at the end of the lipid tails whose polar-nonpolar interfaces lie at (
1,
2,
3,
4), respectively. For graphical clarity we write x′ ≡ x + dx and y′ ≡ y + dy. z(1) is shown in green.
FIG.13 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
(1) extends from z(1) (green) to the head group surface (yellow). The coloring corresponds to Fig. 3.
FIG.14 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
FIG.15 Download High Resolution Image (.zip file) |
Export Figure to PowerPoint
Tables
,
}.
View Table
View Table
View Table















This Publication
Scitation
SPIN
Google Scholar
PubMed