ABSTRACT Experimental evidences have indicated that cholesterol may adapt highly regular lateral distributions (i.e., superlattices) in a phospholipid bilayer. We investigated the formations of superlattices at cholesterol mole fraction of 0.154, 0.25, 0.40, and 0.5 using Monte Carlo simulation. We found that in general, conventional pairwise-additive interactions cannot produce superlattices. Instead, a multibody (nonpairwise) interaction is required. Cholesterol superlattice formation reveals that although the overall interaction between cholesterol and phospholipids is favorable, it contains two large opposing components: an interaction favoring cholesterol-phospholipid mixing and an unfavorable acyl chain multibody interaction that increases nonlinearly with the number of cholesterol contacts. The magnitudes of interactions are in the order of kT. The physical origins of these interactions can be explained by our umbrella model. They most likely come from the requirement for polar phospholipid headgroups to cover the nonpolar cholesterol to avoid the exposure of cholesterol to water and from the sharp decreasing of acyl chain conformation entropy due to cholesterol contact. This study together with our previous work demonstrate that the driving force of cholesterol-phospholipid mixing is a hydrophobic interaction, and multibody interactions dominate others over a wide range of cholesterol concentration.
INTRODUCTION
Cholesterol is a major constituent of the mammalian plasma membranes. The molecular interactions between cholesterol and other lipid molecules have been the subjects of many studies (Finegold, 1993). We seek a general picture of cholesterol-phospholipids interaction, which can capture the key molecular interactions, and would allow us to understand a wide range of experimental results or even to predict new phenomena. In this study, we focus on the molecular interactions, which produce an interesting phenomenon in cholesterol containing lipid membranes: cholesterol superlattices.
Regular distribution of certain molecules in a lipid bilayer was initially proposed based on the observation of a series of "kinks" or "dips" in the ratio of excimer-to-monomer fluorescence of pyrene-phosphatidylcholine at some particular mole fractions (Somerharju et al., 1985; Tang and Chong, 1992; Chong et al., 1994). The bulky pyrene moieties were thought to form hexagonal superlattices to maximize separation from each other. Later, fluorescence data on cholesterol/phospholipid mixtures indicated that cholesterol molecules might also form superlattices in lipid bilayers (Chong, 1994; Virtanen et al., 1995; Liu et al., 1997). Recently, it has been suggested that lipid headgroups may also adopt superlattice distribution in phosphatidylethanolamtine/phosphatidylcholine (PE/PC) bilayers (Cheng et al., 1997, 1999).
Numerous superlattice patterns have been suggested based on geometrical symmetry arguments, such as at cholesterol mole fraction (chi^sub chol^) of 0.118, 0.154, 0.20, 0.25, 0.33, 0.40, and 0.5. Sugar et al. (1994) explored the formation of superlattice patterns using Monte Carlo simulations and introduced a long-range pairwise-additive repulsive interaction. They found that the long-range pairwise-additive repulsion can generate a superlattice pattern at chi^sub chol^ = 0.5 but cannot produce large-scale superlattices of any other compositions. Thus, some key issues about the superlattices remain unsolved: How could such crystal-like structures exist in a bilayer without rigid chemical bonds between molecules? What kinds of molecular interactions are generally required to produce cholesterol superlattices? What are the magnitudes of the interaction energies? What are the physical origins of these interactions?
Recently, using x-ray diffraction and novel sample preparation procedures, we have measured the solubility limits of cholesterol in several different phospholipid bilayers (Huang et al., 1999). Interestingly, these solubility limits occur at cholesterol concentrations that correspond to welldefined cholesterol/phospholipid mole ratios: 1/1 for PE bilayers and 2/1 for PC bilayers.
We have developed a model of cholesterol-phospholipid interaction, which explains these discrete solubility limits. Our Monte Carlo simulations showed that pairwise-additive interaction was inadequate. Instead, the data can only be explained if cholesterol molecules take part in certain type of multibody interactions. The unfavorable cholesterol-cholesterol multibody interaction can be explained by our "umbrella model": in a bilayer environment, nonpolar cholesterols must be covered by neighboring polar phospholipid headgroups to avoid the unfavorable free energy of exposing cholesterol to water. In a lattice model, this requirement can only be expressed in terms of multibody (or nonpair-- wise) interactions. At high cholesterol concentrations, this multibody interaction dominates all others. Thus, only those phospholipid/cholesterol lateral distributions, which meet this coverage requirement, would be allowed. As the concentration of cholesterol increases, fewer and fewer lateral distributions become possible. Near the solubility limit, cholesterol and phospholipid molecules can only adapt some highly regular lateral distributions (i.e., superlattices). The solubility limit is reached when surrounding phospholipid headgroups can no longer completely cover any more cholesterol: the chemical potential of cholesterol jumps steeply, which leads to cholesterol crystal precipitation. Our model predicted that depending on the ability of phospholipid headgroups covering the neighboring cholesterol, cholesterol precipitation is most likely to occur near three discrete values of cholesterol mole fraction, 0.50, 0.57, and 0.67, which correspond to cholesterol/phospholipid mole ratios of 1/1, 4/3, and 2/1, respectively. Thus, the hydrophobic interaction has been implicated as the key driving force in the lateral organization of cholesterol in biomembranes (Huang and Feigenson, 1999).
The author thanks Dr. Gerald W. Feigenson for many productive discussions. This work was supported by the National Science Foundation Grant MCB-9722818 and the Petroleum Research Fund Grant PRF-34872-AC7.
REFERENCES
Almeida, P. F. F., W. L. C. Vaz, and T. E. Thompson. 1993. Percolation and diffusion in three-component lipid bilayers: effect of cholesterol on an equimolar mixture of two phosphatidylcholines. Biophys. J. 64: 399-412.
Cheng, K. H., M. Ruonala, J. Virtanen, and P. Somerharju. 1997. Evidence for superlattice arrangements in fluid phosphatidylcholine/ phosphatidylethanolamine bilayers. Biophys. J. 73:1967-1976.
Cheng, K. H., J. Virtanen, and P. Somerharju. 1999. Fluorescence studies of dehydroergosterol in phosphatidylethanolamine phosphatidylcholine bilayers. Biophys. J 77:3108-3119.
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:1104-1114.
Chong, P. L.-G. 1994. Evidence for regular distribution of sterols in liquid crystalline phosphatidylcholine bilayers. Proc. Natl. Acad. Sci. U.S.A. 91:10069-10073.
Chong, P. L.-G., D. Tang, and I. P. Sugar. 1994. Exploration of physical principles underlying lipid regular distribution: effects of pressure, temperature, and radius of curvature on E/M dips in pyrene-labeled PC/ DMPC binary mixtures. Biophys. J. 66:2029-2038.
Demel, R. A., L. L. M. van Deenen, and B. A. Pethica. 1967. Monolayer interactions of phospholipids and cholesterol. Biochim. Biophys. Acta. 135:11-19.
Feigenson, G. W, and J. T. Buboltz. 2001. Ternary phase diagram of dipalmitoyl-PC/dilauroyl-PC/cholesterol: nanoscopic domain formation driven by cholesterol. Biophys. J. 80:2775-2788.
Finegold, L., editor, 1993. Cholesterol in membrane models. CRC Press, Boca Raton, FL.
Huang, J., J. T. Buboltz, and G. W. Feigenson. 1999. Maximum solubility of cholesterol in phosphatidylcholine and phosphatidylethanolamine bilayers. Biochim. Biophys. Acta. 1417:89-100.
Huang, J., and G. W. Feigenson. 1999. A microscopic interaction model of maximum solubility of cholesterol in lipid bilayers. Biophys. J. 76: 2142-2157.
Kinsky, S. C., S. A. Lose, D. Zopf, L. L. M. van Deepen, and J. Haxby. 1967. Interaction of filipin and derivatives with erythrocyte membranes and lipid dispersions: electron microscopic observations. Biochim. Biophys. Acta. 135:844-861.
Leathes, J. B. 1925. Role of fats in vital phenomena. Lancet. 208:853-856. Liu, F., I. P. Sugar, and P. L.-G. Chong. 1997. Cholesterol and ergosterol superlattices in three-component liquid crystalline lipid bilayers as revealed by dehydroergosterol fluorescence. Biophys. J. 72:2243-2254.
Metropolis, N., A. W. Rosenbluth, M. N. Rosenbluth, and A. Teller. 1953. Equation of state calculations by fast computing machines. J. Chem. Phys. 21:1087-1092.
Moore, W. J. 1962. Physical Chemistry, 3rd Ed. Prentice-Hall, Englewood Cliffs, NJ. 128.
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: 2075-2089.
Somerharju, P., J. A. Virtanen, and K. H. Cheng. 1999. Lateral organization of membrane lipids: the superlattice view. Biochim. Biophys. Acta. 1440:32-48.
Somerharju, P., J. A. Virtanen, K. K. Eklund, P. Vainio, and P. K. Kinnunen. 1985. 1-Palmitoyl-2-pyrenedecanoyl glycerophospholipids as
membrane probes: evidence for regular distribution in liquid-crystalline phosphatidylcholine bilayers. Biochemistry. 24:2773-2781.
Spink, C. H., S. Manley, and M. Breed. 1999. Thermodynamics of transfer of cholesterol from gel to fluid phases of phospholipid bilayers. Biochim. Biophys. Acta. 1279:190-196.
Stockton, B. W., and I. C. P. Smith. 1976. A deuterium NMR study of the condensing effect of cholesterol on egg phosphatidylcholine bilayer membranes. Chem. Phys. Lipids. 17:251-261.
Sugar, I. P., D. Tang, and P. L.-G. Chong. 1994. Monte Carlo simulation of lateral distribution of molecules in a two-component lipid membrane: effect of long-range repulsive interactions. J. Phys. Chem. 98: 7201-7210.
Tang, D., and P. L.-G. Chong. 1992. E/M dips: evidence for lipids regularly distributed into hexagonal super-lattices in pyrene-PC/DMPC binary mixtures at specific concentrations. Biophys. J. 63:903-910.
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:2147-2156.
Virtanen, J. A., M. Ruonala, M. Vauhkonen, and P. Somerharju. 1995. Lateral organization of liquid-crystalline cholesterol-dimyristoylphosphatidylcholine bilayers: evidence for domains with hexagonal and centered rectangular cholesterol superlattices. Biochemistry. 34: 11568-11581.
Vist, M., and J. H. Davis. 1990. Phase equilibria of cholesterol/ dipalmitoylphosphatidylcholine mixtures: ZH nuclear magnetic resonance and differential scanning calorimetry. Biochemistry. 29:451-464.
Wang, M. M., I. P. Sugar, and P. L.-G. Chong. 2002. Effect of double bond position on dehydroergosterol fluorescence intensity dips in phosphatidylcholine bilayers with saturated sn-I and monenoic sn-2 acyl chains. J. Phys. Chem., In press.
Zuckermann, M. J., J. H. Ipsen, and 0. G. Mouritsen. 1993. Theoretical studies of the phase behavior of lipid bilayers containing cholesterol. In Cholesterol in Membrane Models. L. Finegold, editor. CRC Press, Boca Raton, FL. 223-257.
Juyang Huang
Department of Physics, Texas Tech University, Lubbock, Texas 79409 USA
Submitted January 22, 2002, and accepted for publication April 3, 2002.
Address reprint requests to Juyang Huang, Department of Physics, Texas Tech University, Box 41051 Lubbock, TX 79409. Tel: 806-742-4780; Fax: 806-742-1182; E-mail: juhuang@ttacs.ttu.edu.
Copyright Biophysical Society Aug 2002
Provided by ProQuest Information and Learning Company. All rights Reserved