Lateral organization in lipid-cholesterol mixed bilayers

Interactions between lipid and cholesterol molecules in membranes play an important role in the structural and functional properties of cell membranes. Although structural properties of lipid-cholesterol mixtures have been extensively studied, an understanding of the role of cholesterol in the lateral organization of bilayers has been elusive. In this article, we propose a simple yet powerful model, based on self-consistent mean-field theory and molecular dynamics simulations, for lipid bilayers containing cholesterol. Properties predicted by our model are shown to be in excellent agreement with experimental data. Our model predicts that cholesterol induces structural changes in the bilayer through the formation of regions of ordered lipids surrounding each cholesterol molecule. We find that the "smooth" and "rough" sides of cholesterol play crucial roles in formation and distribution of the ordered regions. Our model is predictive in that input parameters are obtained from independent atomistic molecular dynamics simulations. The model and method are general enough to describe other heterogeneous lipid bilayers, including lipid rafts.     

Adhesion and hemifusion of cytoplasmic myelin lipid membranes are highly dependent on the lipid composition

We report the effects of calcium ions on the adhesion and hemifusion mechanisms of model supported myelin lipid bilayer membranes of differing lipid composition. As in our previous studies Min et al. [1,2], the lipid compositions used mimic "healthy" and "diseased-like" (experimental autoimmune encephalomyelitis, EAE) membranes. Our results show that the interaction forces as a function of membrane separation distance are well described by a generic model that also (and in particular) includes the hydrophobic interaction arising from the hydrophobically exposed (interior) parts of the bilayers. The model is able to capture the mechanical instability that triggers the onset of the hemifusion event, and highlights the primary role of the hydrophobic interaction in membrane fusion. The effects of lipid composition on the fusion mechanism, and the adhesion forces between myelin lipid bilayers, can be summarized as follows: in calcium-free buffer, healthy membranes do not present any signs of adhesion or hemifusion, while diseased membranes hemifuse easily. Addition of 2mM calcium favors adhesion and hemifusion of the membranes independently of their composition, but the mechanisms involved in the two processes were different: healthy bilayers systematically presented stronger adhesion forces and lower energy barriers to fusion compared to diseased bilayers. These results are of particular relevance for understanding lesion development (demyelination, swelling, vacuolization and/or vesiculation) in myelin associated diseases such as multiple sclerosis and its relationship to lipid domain formation in myelin membranes.     

Quantification of phase transitions of lipid mixtures from bilayer to non-bilayer structures: Model,experimental validation and implication on membrane fusion

Lipid bilayers provide a solute-proof barrier that is widely used in living systems. It has long been recognized that the structural changes of lipids during the phase transition from bilayer to non-bilayer have striking similarities with those accompanying membrane fusion processes. In spite of this resemblance, the numerous quantitative studies on pure lipid bilayers are difficult to apply to real membranes. One reason is that in living matter, instead of pure lipids, lipid mixtures are involved and there is currently no model that establishes the connection between pure lipids and lipid mixtures. Here, we make this connection by showing how to obtain (i) the short-range repulsion between bilayers made of lipid mixtures and, (ii) the pressure at which transition from bilayer phase to non-bilayer phases occur. We validated our models by fitting the experimental data of several lipid mixtures to the theoretical data calculated based on our model. These results provide a useful tool to quantitatively predict the behavior of complex membranes at low hydration.     

Studies of short-wavelength collective molecular motions in lipid bilayers using high resolution inelastic X-ray scattering

We summarize a series of experimental results made with the newly developed high resolution X-ray scattering (IXS) instrument on two pure lipid bilayers, including dimyristoylphosphatidylcholine (DMPC) and dilauroylphosphatidylcholine (DLPC) in both gel and liquid crystal phases, and lipid bilayers containing cholesterol. By analyzing the IXS data based on the generalized three effective eigenmode model (GTEE), we obtain dispersion relations of the high frequency density oscillations (phonons) of lipid molecules in these bilayers. We then compare the dispersion relations of pure lipid bilayers of different chain lengths among themselves and the dispersion relations of pure lipid bilayers with those of the cholesterol containing bilayers. We also compare our experimental results with collective dynamics data generated by computer molecular dynamics (MD) simulations for dipalmitoylphosphatidylcholine (DPPC) in gel phase and DMPC in liquid crystal phase.     

Lipid-cholesterol interactions in the P beta' phase. Application of a statistical mechanical model.

We describe a statistical mechanical model for lipid-cholesterol mixtures in the P beta' (ripple) phase of lipid bilayers. The model is a simple extension of an earlier model for the ripple phase in pure lipid bilayers. The extension consists of adding a degree of freedom to allow for the occupation of underlying lattice sites by cholesterol molecules, and adding a lipid-cholesterol interaction term to the model Hamiltonian. The interaction term was constructed based on numerical calculations of lipid-cholesterol energies for several different packing juxtapositions of the two molecules. Other than the lipid-cholesterol interactions, the extended model uses the same parameter set as the earlier model, so that comparison of the properties of the extended model with experimental data serves as a test of the validity of the original model. Properties of the model were calculated using the Monte Carlo method. Results are displayed as snapshots of the ripple configurations at different cholesterol concentrations. The spacing of the ripples increases with increasing cholesterol concentration and the rate of increase compares very well with experimental data. The success of this model supports the conclusion drawn earlier that frustration arising from anisotropic packing interactions is responsible for the ripple phase in lipid bilayers. In the extended model these packing interactions are responsible for the selective partitioning of cholesterol in the regions between the ripples.     

Combined influence of cholesterol and synthetic amphiphillic peptides upon bilayer thickness in model membranes.

Deuterium (2H) NMR was used to study bilayer hydrophobic thickness and mechanical properties when cholesterol and/or synthetic amphiphillic polypeptides were added to deuterated POPC lipid bilayer membranes in the liquid-crystalline (fluid) phase. Smoothed acyl chain orientational order profiles were used to calculate bilayer hydrophobic thickness. Addition of 30 mol% cholesterol to POPC at 25 degrees C increased the bilayer thickness from 2.58 to 2.99 nm. The peptides were chosen to span the bilayers with more or less mismatch between the hydrophobic peptide length and membrane hydrophobic thickness. The average thickness of the pure lipid bilayers was significantly perturbed upon addition of peptide only in cases of large mismatch, being increased (decreased) when the peptide hydrophobic length was greater (less) than that of the pure bilayer, consistent with the "mattress" model of protein lipid interactions (Mouritsen, O.G., and M. Bloom. 1984. Biophys. J. 46:141-153). The experimental results were also used to examine the combined influence of the polypeptides and cholesterol on the orientational order profile and thickness expansivity of the membranes. A detailed model for the spatial distribution of POPC and cholesterol molecules in the bilayers was proposed to reconcile the general features of these measurements with micromechanical measurements of area expansivity in closely related systems. Experiments to test the model were proposed.     

Insertion and pore formation driven by adsorption of proteins onto lipid bilayer membrane-water interfaces

We describe the binding of proteins to lipid bilayers in the case for which binding can occur either by adsorption to the lipid bilayer membrane-water interface or by direct insertion into the bilayer itself. We examine in particular the case when the insertion and pore formation are driven by the adsorption process using scaled particle theory. The adsorbed proteins form a two-dimensional "surface gas" at the lipid bilayer membrane-water interface that exerts a lateral pressure on the lipid bilayer membrane. Under conditions of strong intrinsic binding and a high degree of interfacial converge, this pressure can become high enough to overcome the energy barrier for protein insertion. Under these conditions, a subtle equilibrium exists between the adsorbed and inserted proteins. We propose that this provides a control mechanism for reversible insertion and pore formation of proteins such as melittin and magainin. Next, we discuss experimental data for the binding isotherms of cytochrome c to charged lipid membranes in the light of our theory and predict that cytochrome c inserts into charged lipid bilayers at low ionic strength. This prediction is supported by titration calorimetry results that are reported here. We were furthermore able to describe the observed binding isotherms of the pore-forming peptides endotoxin (alpha 5-helix) and of pardaxin to zwitterionic vesicles from our theory by assuming adsorption/insertion equilibrium.     

Electron spin resonance and steady-state fluorescence polarization studies of lipid bilayers containing integral proteins

We derive equations that describe changes in the steady-state fluorescence polarization of the probe 1,6-diphenyl-1,3,5-hexatriene (DPH) or in the spectrum of electron spin resonance (ESR) nitroxide spin-labeled lipid probes as a function of the intrinsic molecule concentration in lipid bilayer membranes. We make use of an assumption used by us in an earlier paper. The equations are independent of any membrane model. They are valid when a DPH probe or a spin-labeled chain is equivalent to an unlabeled lipid hydrocarbon chain only as far as their general space-filling properties are concerned. We consider cases where the bilayer is either in a single homogeneous phase or in a two-phase region. We apply our equations to analyze ESR data from delipidated sarcoplasmic reticulum membranes and from egg yolk phosphatidylcholine bilayers containing Ca2+-ATPase, and DPH data from dipalmitoylphosphatidylcholine (DPPC) bilayers containing Ca2+-ATPase, both for T greater than Tc. The following conclusions were derived: (i) Ca2+-ATPase oligomers are "randomly" distributed, for the concentrations studied, in the fluid phase. (ii) There is no fixed stoichiometric ratio of "boundary" lipids and oligomers. (iii) Between 24k and 28k lipid molecules are able to surround each isolated oligomer composed of k Ca2+-ATPase monomers. Finally, we apply our equations to analyze DPH studies on DPPC bilayers containing Ca2+-ATPase for T less than Tc. We find that the results reported are in accord with the predictions of the model. In the Appendix, we show that an analytical expression for probabilities used by us is in very good agreement with the results of computer simulation.     

Galactosylceramide domain microstructure: impact of cholesterol and nucleation/growth conditions

Galactosylceramide (GalCer), a glycosphingolipid, is believed to exist in the extracellular leaflet of cell membranes in nanometer-sized domains or rafts. The local clustering of GalCer within rafts is thought to facilitate the initial adhesion of certain viruses, including HIV-1, and bacteria to cells through multivalent interactions between receptor proteins (gp120 for HIV-1) and GalCer. Here we use atomic force microscopy (AFM) to study the effects of cholesterol on solid-phase GalCer domain microstructure and miscibility with a fluid lipid 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) in supported lipid bilayers. Using "slow-cooled vesicle fusion" to prepare the supported lipid bilayers, we were able to overcome the nonequilibrium effects of the substrate (verified by comparison to results for giant unilamellar vesicles) and accurately quantify the dramatic effect of cholesterol on the GalCer domain surface area/perimeter ratio (A(D)/P) and DLPC-GalCer miscibility. We compare these results to a supported lipid bilayer system in which the bilayer is rapidly cooled (nonequilibrium conditions), "quenched vesicle fusion", and find that the microstructures are remarkably similar above a cholesterol mol fraction of approximately 0.06. We determined that GalCer domains were contained in one leaflet distal to the mica substrate through qualitative binding experiments with Trichosanthes kirilowii agglutinin (TKA), a galactose-specific lectin, and AFM of Langmuir-Blodgett deposited GalCer/DLPC supported lipid bilayers. In addition, GalCer domains in bilayers containing cholesterol rearranged upon tip-sample contact. Our results further serve to clarify why discrepancies exist between different model membrane systems and between model membranes and cell membranes. In addition, these results offer new insight into the effect of cholesterol and surrounding lipid on domain microstructure and behavior. Finally, our observations may be pertinent to cell membrane structure, dynamics, and HIV infection.     

Theoretical conformational analysis of phospholipids bilayers

We present a computational approach describing the conformation of lipid molecules (1-2-dipalmitoyl-sn-glycero-3 phosphocholine (DPPC)) organized in bilayers. The classical semi-empirical method used in peptide conformational analysis has been extended successfully to lipids. The excellent agreement between our theoretical predictions and recent experimental data on the molecular organization of lipid bilayers suggests that the method could be a valuable tool in the lipid conformational analysis but also in the prediction of orientation and mode of insertion of amphiphilic molecules into the lipid bilayer.     

Interaction of LL-37 with model membrane systems of different complexity: influence of the lipid matrix

As the main difference between bacterial and mammalian cell membranes is their net charge, the focal point of consideration in many model membrane experiments with antimicrobial peptides is lipid headgroup charge. We studied the interaction of the human multifunctional peptide LL-37 with single phospholipid monolayers, bilayers, and bilayers composed of binary mixtures of the four phospholipid species predominantly used in model membrane experiments (phosphatidylcholine, phosphatidylethanolamine, phosphatidylglycerol, and phosphatidylserine). We found that 1), the effects on single lipid monolayers are not comparable to those on the corresponding bilayers; 2), there are four different effects of LL-37 on bilayers of the four lipids; 3), the preference of LL-37 for the specific lipids is roughly inversely related to chain packing density; and 4), in the binary lipid mixtures, one lipid—and not necessarily the charged one—generally governs the mode of lipid/peptide interaction. Thus, our results show that lipid net charge is not the decisive factor determining the membrane-perturbing mechanism of LL-37, but only one of several parameters, among them packing density, the ability to form intermolecular H-bonds, and lipid molecular shape, which emphasizes how profoundly the choice of the model system can influence the outcome of a study of lipid/peptide interaction.     

Structure of lipid bilayers

The quantitative experimental uncertainty in the structure of fully hydrated, biologically relevant, fluid (L(alpha)) phase lipid bilayers has been too large to provide a firm base for applications or for comparison with simulations. Many structural methods are reviewed including modern liquid crystallography of lipid bilayers that deals with the fully developed undulation fluctuations that occur in the L(alpha) phase. These fluctuations degrade the higher order diffraction data in a way that, if unrecognized, leads to erroneous conclusions regarding bilayer structure. Diffraction measurements at high instrumental resolution provide a measure of these fluctuations. In addition to providing better structural determination, this opens a new window on interactions between bilayers, so the experimental determination of interbilayer interaction parameters is reviewed briefly. We introduce a new structural correction based on fluctuations that has not been included in any previous studies. Updated measurements, such as for the area compressibility modulus, are used to provide adjustments to many of the literature values of structural quantities. Since the gel (L(beta)') phase is valuable as a stepping stone for obtaining fluid phase results, a brief review is given of the lower temperature phases. The uncertainty in structural results for lipid bilayers is being reduced and best current values are provided for bilayers of five lipids.     

Electrophysiological recordings of single ion channels in planar lipid bilayers using a polymethyl methacrylate microfluidic chip

Planar lipid bilayers are used for functional studies of ion channel proteins using electrophysiological techniques. We have been developing a plastic micro-fluidic device for the reconstitution of planar lipid bilayers and electrophysiological recordings toward a "membrane protein chip" for high-throughput screening. In the previous report [Suzuki, H., Tabata, K.V., Noji, H., Takeuchi, S., 2006. Highly reproducible method of planar lipid bilayer reconstitution in polymethyl methacrylate microfluidic chip. Langmuir 22 (4), 1937-1942], we presented the method and device in which the reproducibility of planar lipid bilayers reached 90%, and multiple bilayers were formed simultaneously. In this communication, we show that our device has excellent electric properties suitable for ion channel analysis down to single molecular level. Additional aspects on the optical accessibility and controllability on lipid bilayer formation are also presented.     

Effect of hydrodynamic interactions on the diffusion of integral membrane proteins: tracer diffusion in organelle and reconstituted membranes.

A persistent discrepancy exists between theoretical predictions and experimental observations for the diffusion coefficients of integral membrane proteins in lipid bilayers free of immobilized proteins. Current thermodynamic theories overestimate tracer diffusion coefficients at high area fractions. We explore the hypothesis that the combined effect of hydrodynamic and thermodynamic interactions reconciles theory with experiment. We have determined previously the short- and long-time tracer diffusivities, Ds and Dl, respectively, of integral membrane proteins in lipid bilayers as a function of their area fraction, phi. The results are based on two-particle hydrodynamic and thermodynamic interactions and are precise to O(phi). Here we extend the results for Dl to high phi by combining the hydrodynamic results for Ds into theories for Dl based on many-particle thermodynamic interactions. The results compare favorably with the experimental measurements of Dl as a function of protein area fraction for bacteriorhodopsin in reconstituted membranes and for complex III of the mitochondrial inner membrane. The agreement suggests that both hydrodynamic and thermodynamic interactions are important determinants of diffusion coefficients of proteins in lipid bilayers. Additional experiments are required to verify the role of hydrodynamic interactions in protein diffusion in reconstituted systems.     

The effect of cholesterol on short- and long-chain monounsaturated lipid bilayers as determined by molecular dynamics simulations and X-ray scattering

We investigate the structure of cholesterol-containing membranes composed of either short-chain (diC14:1PC) or long-chain (diC22:1PC) monounsaturated phospholipids. Bilayer structural information is derived from all-atom molecular dynamics simulations, which are validated via direct comparison to x-ray scattering experiments. We show that the addition of 40 mol % cholesterol results in a nearly identical increase in the thickness of the two different bilayers. In both cases, the chain ordering dominates over the hydrophobic matching between the length of the cholesterol molecule and the hydrocarbon thickness of the bilayer, which one would expect to cause a thinning of the diC22:1PC bilayer. For both bilayers there is substantial headgroup rearrangement for lipids directly in contact with cholesterol, supporting the so-called umbrella model. Importantly, in diC14:1PC bilayers, a dynamic network of hydrogen bonds stabilizes long-lived reorientations of some cholesterol molecules, during which they are found to lie perpendicular to the bilayer normal, deep within the bilayer’s hydrophobic core. Additionally, the simulations show that the diC14:1PC bilayer is significantly more permeable to water. These differences may be correlated with faster cholesterol flip-flop between the leaflets of short-chain lipid bilayers, resulting in an asymmetric distribution of cholesterol molecules. This asymmetry was observed experimentally in a case of unilamellar vesicles (ULVs), and reproduced through a set of novel asymmetric simulations. In contrast to ULVs, experimental data for oriented multilamellar stacks does not show the asymmetry, suggesting that it results from the curvature of the ULV bilayers.     

Model answers to lipid membrane questions

Ever since it was discovered that biological membranes have a core of a bimolecular sheet of lipid molecules, lipid bilayers have been a model laboratory for investigating physicochemical and functional properties of biological membranes. Experimental and theoretical models help the experimental scientist to plan experiments and interpret data. Theoretical models are the theoretical scientist's preferred toys to make contact between membrane theory and experiments. Most importantly, models serve to shape our intuition about which membrane questions are the more fundamental and relevant ones to pursue. Here we review some membrane models for lipid self-assembly, monolayers, bilayers, liposomes, and lipid-protein interactions and illustrate how such models can help answering questions in modern lipid cell biology.     

Effects of monovalent anions of the hofmeister series on DPPC lipid bilayers Part I: swelling and in-plane equations of state

Aiming to improve understanding of the mechanisms behind specific anion effects in biological systems we have studied the effects of sodium salts of simple monovalent anions belonging to the Hofmeister series on the bilayers of the zwitterionic lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine using small-angle x-ray scattering and the osmotic stress technique. NaCl, NaBr, NaNO(3), NaI, and NaSCN were used in this investigation. The electrolytes were found to swell the bilayers and to increase the area per lipid headgroup at each value of the osmotic pressure, suggesting the association of anions with the bilayer-lipid interfaces. The effects follow the Hofmeister series with SCN(-) inducing the most pronounced changes. "Ion competition" experiments with mixed NaI/NaCl solutions at total salinity 0.1 and 0.5 M revealed that the effect of ions on the lipid equation-of-state is roughly linear at low concentrations, but strongly nonlinear at high concentrations. The experimental results are fitted in a companion article to provide "binding" or "partitioning" constants of anions in the lipid bilayers.     

Effects of Natural and Enantiomeric Cholesterol on the Thermotropic Phase Behavior and Structure of Egg Sphingomyelin Bilayer Membranes

Phospholipids, sphingolipids, and sterols are the major lipid components of the plasma membranes of eukaryotic cells. Because these three lipid classes occur naturally as enantiomerically pure compounds, enantiospecific lipid-lipid and lipid-sterol interactions could in principle occur in the lipid bilayers of eukaryotic plasma membranes. Although previous biophysical studies of phospholipid and phospholipid-sterol model membrane systems have consistently failed to observe such enantiomerically selective interactions, a recent monolayer study of the interactions of natural and enantiomeric cholesterol with egg sphingomyelin has apparently revealed the existence of enantiospecific sterol-sphingolipid interactions. To determine whether enantiospecific sterol-sphingolipid interactions also occur in more biologically relevant lipid-bilayer systems, differential scanning calorimetric, x-ray diffraction, and neutral buoyant-density measurements were utilized to study the effects of natural and enantiomeric cholesterol on the thermotropic phase behavior and structure of egg sphingomyelin bilayers. The calorimetry experiments show that the natural and enantiomeric cholesterol have essentially identical effects on the temperature, enthalpy, and cooperativity of the gel/liquid-crystalline phase transition of egg sphingomyelin bilayers within the limits of experimental error. As well, the x-ray diffraction and neutral buoyancy experiments indicate that bilayers formed from mixtures of natural or enantiomeric cholesterol and egg sphingomyelin have, within experimental uncertainty, the same structure and mass density. We thus conclude that significant enantioselective cholesterol-sphingolipid interactions do not occur in this lipid-bilayer model membrane system.     

Characterization of the solubilization of lipid bilayers by surfactants

This communication addresses the state of aggregation of lipid-detergent mixed dispersions. Analysis of recently published data suggest that for any given detergent-lipid mixture the most important factor in determining the type of aggregates (mixed vesicles or mixed micelles) and the size of the aggregate is the detergent to lipid molar ratio in these aggregates, herein denoted the effective ratio, Re. For mixed bilayers this effective ratio has been previously shown to be a function of the lipid and detergent concentrations and of an equilibrium partition coefficient, K, which describes the distribution of the detergent between the bilayers and the aqueous phase. We show that, similar to mixed bilayers, the size of mixed micelles is also a function of the effective ratio, but for these dispersions the distribution of detergent between the mixed micelles and the aqueous medium obeys a much higher partition coefficient. In practical terms, the detergent concentration in the mixed micelles is equal to the difference between the total detergent concentration and the critical micelle concentration (cmc). Thus, the effective ratio is equal to this difference divided by the lipid concentration. Transformation of mixed bilayers to mixed micelles, commonly denoted solubilization, occurs when the surfactant to lipid effective ratio reaches a critical value. Experimental evaluation of this critical ratio can be based on the linear dependence of detergent concentration, required for solubilization, on the lipid concentration. According to the 'equilibrium partition model', the dependence of the 'solubilizing detergent concentration' on the lipid concentration intersects with the lipid axis at -1/K, while the slope of this dependence is the critical effective ratio. On the other hand, assuming that when solubilization occurs the detergent concentration in the aqueous phase is approximately equal to the critical micelle concentration, implies that the above dependence intersects with the detergent axis at the critical micelle concentration, while its slope, again, is equal to the critical effective ratio. Analysis of existing data suggests that within experimental error both these distinctively different approaches are valid, indicating that the critical effective ratio at which solubilization occurs is approximately equal to the product of the critical micelle concentration and the distribution coefficient K. Since the nature of detergent affects K and the critical micelle concentration in opposite directions, the critical ('solubilizing') effective ratio depends upon the nature of detergent less than any of these two factors.     

A model for phase transitions in lipid bilayers and biological membranes

We wish to present an order-disorder model for the observed phase transitions in lipid bilayers and biological membranes. We show that the model may, under certain circumstances, exhibit two phase transitions, one corresponding to positional disordering of entire lipid molecules, and the other corresponding to orientational disordering in the hydrocarbon chains. We then give results of our numerical analysis of the model and compare them with experimental data. Shortcomings of the model and future directions for analyses of this type are also discussed.