| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

* Theoretical Biological Physics, Royal Institute of Technology, Stockholm, Sweden; and
Departments of Physics and Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania
Correspondence: Address reprint requests to Olle Edholm, Tel.: 46-8-553-7-8168; E-mail: oed{at}theophys.kth.se.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Even with all the detail provided by simulations, it is not obvious how to obtain the average area of each lipid in a mixed bilayer (2
,3
). As an example, we will consider mixtures of cholesterol and dipalmitoylphosphatidylcholine (DPPC) which were studied by simulations of simplified models more than 10 years ago (4
6
) and more recently in detailed molecular dynamics simulations that include explicit water (7
15
). We will focus on three articles that simulated several concentrations of cholesterol in DPPC (3
,16
,17
). One of the most easily accessible quantities to monitor in an NPT ensemble is the total area A of a bilayer with Nchol cholesterol molecules and NDPPC DPPC molecules in each monolayer. One can then define the area per total lipid a(x) as a function of the mole fraction of cholesterol, defined as x
Nchol/(NDPPC + Nchol),
![]() | (1) |
|
|
![]() | (2) |
| EXISTING METHODS |
|---|
|
|
|---|
Other methods use atomic details to divide the area using Voronoi tessellation. Voronoi tessellation methods come in a variety of types. One type projects the center of mass of each molecule onto a plane and then tessellates the area of the plane, but this has been explicitly criticized for overweighting the area of smaller molecules (2
). Instead they chose to tessellate using atoms near the hydrophobic/hydrophilic interface, three atoms for DPPC and one atom for cholesterol. In addition to the issue of how many atoms to use, this method essentially requires an arbitrary decision regarding which slice of area (i.e., at which z level from the bilayer center) is to be used for the tessellation. Although values for cholesterol/DPPC were not specifically given, the areas appear to be similar to those obtained by Chiu et al. (16
).
A recent article by Falck et al. (3
) shows quite nicely that there is indeed a z-dependence to the molecular areas when examined at an atomic level. Their method uses van der Waals surfaces around the lipid and cholesterol molecules and the leftover area is defined as "free area." This contrasts with the Voronoi tessellation, which assigns all the area to some molecule, and therefore has no free area. For z values near the cholesterol-ring structure, Fig. 11 from Falck et al. (3
) shows achol
0.3 nm2 and aDPPC
0.4 nm2 for x = 0 falling to aDPPC
0.3 nm2 for x = 0.5, but these are the "bare areas," none of which includes any free area.
We have reservations regarding the concept of free area (as well as with the corresponding free volume). One would not expect the free area to be largest in the water outside the bilayer, but this is the result shown in Fig. 12 of Falck et al. (3
). More generally, free area and free volume will be nonzero for two reasons. First, even a close-packed hard-sphere system will have a free volume corresponding to approximately one-third of the volume (smallest, 26% for a face-centered cubic lattice). Secondly, this lower free-volume limit is only reached at T = 0 for any condensed-matter system. Even in perfect crystals the anharmonicity of the interaction potentials generally lead to thermal expansion, but such expansion does not create any free regions into which other molecules could be inserted. That only occurs in crystals when vacancy defects occur, which is not an appropriate concept for fluids such as water or L
phase lipid bilayers. We therefore distrust definitions of free areas and volumes that rely only on the bare energetic description of molecules and that ignore the free energy aspect that pertains when T > 0. One motivation for a free-area description lies in its application to diffusion. After careful comparison, Falck et al. (3
) conclude that the free-area theory "tends to underestimate the changes in the values of the lateral diffusion coefficients" and that "it seems unnecessary to aim for a quantitative description with such a simple framework."
A different motivation for examining the variation of properties in the perpendicular direction is that such variations may affect the preference for forming different kinds of structures like bilayers, micelles, or hexagonal phases. Furthermore, it may modulate the free energies of conformational changes in functioning membrane proteins (18
). Such variations can, however, be better addressed by the lateral pressure profile perpendicular to the membrane (19
), which can be calculated in simulations (20
).
A third method for partitioning the area that was suggested by Hofsäß et al. (17
) could be clarified and phrased in the following way. The major assumption is that the bilayer can be characterized by a common thickness h(x) that can be expressed in terms of volume and area per lipid as
![]() | (3) |
![]() | (4) |
![]() | (5) |
3% (21
30%, from 0.47 nm2 in the gel phase (22
We now present an improvement of the volumetric part of the analysis of Hofsäß et al. (17
). This may be of some independent interest because volumes of the chemical components in mixtures are important, for example, in estimating electron density and neutron scattering lengths for diffraction experiments. We first define the volume per lipid v(x) by
![]() | (6) |
![]() | (7) |
|
| DEFINITION OF PARTIAL-SPECIFIC AREA |
|---|
|
|
|---|
![]() | (8) |
i. Because A is a homogeneous function of the Ni variables, it is guaranteed that
![]() | (9) |
Ni. (This is easily shown by taking the derivative with respect to
of both sides of A(
N1, ...,
Nm) =
A(N1, ..., Nm).) We apply this formalism to two-component mixtures of DPPC and cholesterol by dividing A(X) in Eq. 9 by NDPPC + Nchol, which then gives Eq. 2. We note that adding more water does not increase the total area of the membrane for fully hydrated membranes (excess water condition), so the specific area of water is identically zero, and therefore the water component can be ignored in Eq. 2.
Eq. 8 makes it clear that it is necessary to perform simulations for a set of values of x to determine the partial-specific areas ai(x). Given the total area per molecule a(x) from such a set of simulations, a convenient graphical way to obtain the ai(x) rewrites Eq. 2 as
![]() | (10) |
![]() | (11) |
| RESULTS FOR MIXTURES OF DPPC AND CHOLESTEROL |
|---|
|
|
|---|
These averages are shown in Fig. 1 compared to the individual simulated data sets from which they were derived.
Fig. 3 shows the plot suggested by Eq. 10 that can be used to obtain the specific area of cholesterol from the slope of the local tangent and the area of DPPC from the intercept of this tangent at x = 0. Fig. 3 shows two linear regimes that can be fitted with separate straight lines. For the smallest x the plot gives an area per DPPC of 0.64 nm2, consistent with the values for pure DPPC (1
) and an area per cholesterol of 0.81 nm2. One should, however, note that the three different simulations give quite different specific areas per cholesterol in this region. Using the data in Table 1 for pure DPPC and the lowest cholesterol concentration, achol(0) is 2.0 nm2 from Chiu et al. (16
), 0.71 nm2 from Hofsäß et al. (17
), and 0.11 nm2 from Falck et al. (3
). Nevertheless, even if the spread is quite large, due partly to difficulty in evaluating numerical derivatives, it is clear that all simulations give negative values of achol(0). For x > 0.3 the average curve gives achol = 0.26 nm2 and aDPPC = 0.51 nm2, and the difference between the individual simulations is <±0.04 nm2. In the intermediate region 0.07 < x < 0.30, the slope is small and positive and it might be tempting to consider a third region with the partial-specific areas achol = 0.095 nm2 and aDPPC = 0.58 nm2.
|
![]() | (12) |
|
| DISCUSSION AND CONCLUSIONS |
|---|
|
|
|---|
If the perturbation regions around each cholesterol molecule had specific cutoff lengths, then aDPPC(x) would remain constant at aDPPC(0) until these regions covered the entire surface. Of course, the perturbations decay smoothly instead of having rigid cutoff lengths, so aDPPC(x) gradually decreases with increasing x. Nevertheless, Fig. 4 suggests that there remain relatively unperturbed DPPC molecules up to x
0.1. However, in the high cholesterol regime as x approaches 0.5, aDPPC(x) approaches 0.50 nm2 and achol(x) approaches the constant value 0.27 nm2. These particular results, derived from the general partial-specific-area approach, suggest a more specific model for mixtures of cholesterol and DPPC that is presented in the Appendix.
The primary motivation for the partial-specific-area approach that we advocate in this article is that it involves none of the arbitrary assumptions that have had to be made in recent atomistic definitions. It is perfectly defined thermodynamically, with a well-founded and historical conceptual framework. It very elegantly displays the well-known condensing effect of cholesterol on lipid bilayers. However, we do not advocate its use exclusively. The other definitions have value for obtaining atomistic perspective, especially into the z variations along the bilayer normal. Nevertheless, we suggest that partial-specific area is the appropriate canonical quantity to report from simulations of mixtures that are performed as a function of mole fractions. It should also be considered for any experiments that are capable of obtaining area a(x).
| APPENDIX: A PARTICULAR MODEL FOR AREAS IN CHOLESTEROL/DPPC MIXTURES |
|---|
|
|
|---|
for all x; this parameter has been the goal of the atomistic approaches. The second parameter is the area
for DPPC when x = 0. The third parameter is the amount
a of area condensation of any DPPC molecule in contact with cholesterol. The fourth parameter is the maximum number n of DPPC molecules that can be condensed by a single cholesterol molecule. The partial-specific area results suggest that, for small x,
![]() | (13) |
at x = 0 gives the partial-specific area of cholesterol. For high cholesterol concentration, all DPPC are in contact with cholesterol and the results for the partial-specific areas suggest that the system behaves more like an ideal mixture with (
) being the area per DPPC and
the area per cholesterol; that is, for large x,
![]() | (14) |
The goal of the model is to provide a single function for all x that includes both these limits. To do this we assume that the cholesterol molecules are randomly distributed and if there is at least one cholesterol within distance d from a phospholipid, then the phospholipid becomes ordered and its area is reduced by
a. The maximum number of phospholipids, n, that become ordered by a single cholesterol is then given by
![]() | (15) |
![]() | (16) |
Let N and A go to infinity for fixed non-zero values of the other parameters, and using the fundamental definition of e (the base of the natural logarithm) gives
![]() | (17) |
![]() | (18) |
a(1 P0) gives the area reduction of each DPPC molecule that is within distance d of a cholesterol. Equations 17 and 18 then give the area per lipid as
![]() | (19) |
a = 0.14 nm2, and n = 7.5. The continuous curves corresponding to these values of the parameters are shown in Figs. 3 and 4, which indicate that this simple model is adequate to fit the simulation data. It is interesting to note that the value of n (which is the most uncertain parameter value with an error of the order ±2) is comparable to the number of nearest-neighbor phospholipids. The value of
is within the range of values that have been obtained from the atomistic approaches described in Existing Methods. Even though the model uses constant
and a mole-fraction-dependent decrease in the DPPC area, the cholesterol partial-specific-area achol(x) has a negative value of 0.78 nm2 at low mole fractions of cholesterol. Therefore, this model encompasses both the partial-specific-area perspective that emphasizes the dramatic condensation effect of cholesterol and the more atomistic perspective that seeks a fixed area for cholesterol. | ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Submitted on April 11, 2005; accepted for publication June 14, 2005.
| REFERENCES |
|---|
|
|
|---|
2. Pandit, S. A., S. Vasudevan, S. W. Chiu, R. J. Mashl, E. Jakobsson, and H. L. Scott. 2004. Sphingomyelin-cholesterol domains in phospholipid lipid membranes: atomistic simulation. Biophys. J. 87:10921100.
3. Falck, E., M. Patra, M. Karttunen, M. Hyvönen, and I. Vattulainen. 2004. Lessons of slicing membranes: interplay of packing, free area, and lateral diffusion in phospholipid/cholesterol bilayers. Biophys. J. 87:10761091.
4. Scott, H. L., and S. Kalaskar. 1989. Lipid chains and cholesterol in membranes: a Monte Carlo study. Biochemistry. 28:36873691.[CrossRef][Medline]
5. Scott, H. L. 1991. Lipid cholesterol interactions: Monte Carlo simulations and theory. Biophys. J. 59:445455.
6. Edholm, O., and A. Nyberg. 1992. Cholesterol in model membranes. Biophys. J. 63:10811089.
7. Robinson, A. J., G. Richards, P. J. Thomas, and M. M. Hann. 1995. Behavior of cholesterol and its effect on headgroup and chain conformations in lipid bilayers: a molecular dynamics study. Biophys. J. 68:164170.
8. Tu, K., M. L. Klein, and D. J. Tobias. 1998. Constant-pressure molecular dynamics investigation of cholesterol effects in a dipalmitoylphosphatidylcholine bilayer. Biophys. J. 75:21472156.
9. Smondyrev, A. M., and M. L. Berkowitz. 1999. Structure of dipalmitoylphosphatidylcholine/cholesterol bilayer at low and high cholesterol concentrations: molecular dynamics simulation. Biophys. J. 77:20752089.
10. Smondyrev, A. M., and M. L. Berkowitz. 2000. Molecular dynamics simulation of dipalmitoylphosphatidylcholine membrane with cholesterol sulfate. Biophys. J. 78:16721680.
11. Pasenkiewicz-Gierula, M., T. Róg, K. Kitamura, and A. Kusumi. 2000. Cholesterol effects on the phosphatidylcholine bilayer polar region: a molecular simulation study. Biophys. J. 78:13761389.
12. Smondyrev, A. M., and M. L. Berkowitz. 2001. Molecular dynamics simulation of the structure of dipalmitoylphosphatidylcholine bilayers with cholesterol, ergosterol and lanosterol. Biophys. J. 80:16491658.
13. Chiu, S. W., E. Jakobsson, and H. L. Scott. 2001. Combined Monte Carlo and molecular dynamics simulation of hydrated lipid-cholesterol lipid bilayers at low cholesterol concentration. Biophys. J. 80:11041114.
14. Chiu, S. W., E. Jakobsson, and H. L. Scott. 2001. Combined Monte Carlo and molecular dynamics simulation of hydrated dipalmitoyl-phosphatidylcholine-cholesterol lipid bilayers. J. Chem. Phys. 114:54355443.[CrossRef]
15. Róg, T., and M. Pasenkiewicz-Gierula. 2001. Cholesterol effects on the phosphatidylcholine bilayer nonpolar region: a molecular simulation study. Biophys. J. 81:21902202.
16. Chiu, S. W., E. Jakobsson, R. J. Mashl, and H. L. Scott. 2002. Cholesterol-induced modifications in lipid bilayers: a simulation study. Biophys. J. 83:18421853.
17. Hofsäß, C., E. Lindahl, and O. Edholm. 2003. Molecular dynamics simulations of phospholipid bilayers with cholesterol. Biophys. J. 84:21922206.
18. Gruner, S., and E. Shyamsunder. 1991. Is the mechanism of general anesthesia related to lipid membrane spontaneous curvature? Ann. N. Y. Acad. Sci. 625:685697.[Abstract]
19. Cantor, R. S. 1997. The lateral pressure profile in membranes: a physical mechanism of general anesthesia. Biochemistry. 36:23392344.[CrossRef][Medline]
20. Lindahl, E., and O. Edholm. 2000. Spatial and energetic-entropic decomposition of surface tension in lipid bilayers from molecular dynamics simulations. J. Chem. Phys. 113:38823893.[CrossRef]
21. Nagle, J. F., and D. A. Wilkinson. 1978. Lecithin bilayers: density measurements and molecular interactions. Biophys. J. 23:159175.
22. Tristram-Nagle, S., Y. Liu, J. Legleiter, and J. F. Nagle. 2002. Structure of gel phase DMPC determined by x-ray diffraction. Biophys. J. 83:33243335.
23. Kucerka, N., Y. Liu, N. Chu, H. I. Petrache, S. Tristram-Nagle, and J. F. Nagle. 2005. Structure of fully hydrated fluid phase DMPC and DLPC lipid bilayers using x-ray scattering from oriented multilamellar arrays and from unilamellar vesicles. Biophys. J. 88:26262637.
24. Moore, W. J. 1962. Physical Chemistry. Prentice-Hall, Englewood Cliffs, NJ. 118122.
25. Vist, M. R., and J. H. Davis. 1990. Phase equilibria of cholesterol/dipalmitoyl-phosphatidylcholine mixtures: 2H nuclear magnetic resonance and differential scanning calorimetry. Biochemistry. 29:451464.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
T. T. Mills, G. E. S. Toombes, S. Tristram-Nagle, D.-M. Smilgies, G. W. Feigenson, and J. F. Nagle Order Parameters and Areas in Fluid-Phase Oriented Lipid Membranes Using Wide Angle X-Ray Scattering Biophys. J., July 15, 2008; 95(2): 669 - 681. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Hodzic, M. Rappolt, H. Amenitsch, P. Laggner, and G. Pabst Differential Modulation of Membrane Structure and Fluctuations by Plant Sterols and Cholesterol Biophys. J., May 15, 2008; 94(10): 3935 - 3944. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Leekumjorn, Y. Wu, A. K. Sum, and C. Chan Experimental and Computational Studies Investigating Trehalose Protection of HepG2 Cells from Palmitate-Induced Toxicity Biophys. J., April 1, 2008; 94(7): 2869 - 2883. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Mathai, S. Tristram-Nagle, J. F. Nagle, and M. L. Zeidel Structural Determinants of Water Permeability through the Lipid Membrane J. Gen. Physiol., December 31, 2007; 131(1): 69 - 76. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. F. Nagle, J. C. Mathai, M. L. Zeidel, and S. Tristram-Nagle Theory of Passive Permeability through Lipid Bilayers J. Gen. Physiol., December 31, 2007; 131(1): 77 - 85. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Richardson, L. J. Cummings, H. J. Harris, and P. O'Shea Toward a Mathematical Model of the Assembly and Disassembly of Membrane Microdomains: Comparison with Experimental Models Biophys. J., June 15, 2007; 92(12): 4145 - 4156. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Zhao, T. Rog, A. A. Gurtovenko, I. Vattulainen, and M. Karttunen Atomic-Scale Structure and Electrostatics of Anionic Palmitoyloleoylphosphatidylglycerol Lipid Bilayers with Na+ Counterions Biophys. J., February 15, 2007; 92(4): 1114 - 1124. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Y. Bhide, Z. Zhang, and M. L. Berkowitz Molecular Dynamics Simulations of SOPS and Sphingomyelin Bilayers Containing Cholesterol Biophys. J., February 15, 2007; 92(4): 1284 - 1295. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Pandit, S.-W. Chiu, E. Jakobsson, A. Grama, and H. L. Scott Cholesterol Surrogates: A Comparison of Cholesterol and 16:0 Ceramide in POPC Bilayers Biophys. J., February 1, 2007; 92(3): 920 - 927. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Pandit, G. Khelashvili, E. Jakobsson, A. Grama, and H. L. Scott Lateral Organization in Lipid-Cholesterol Mixed Bilayers Biophys. J., January 15, 2007; 92(2): 440 - 447. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |