Small-scale lipid-membrane structure: simulation versus experiment.

Recently, it has become obvious that the conventional picture of the fluid lipid-bilayer component of biological membranes being a fairly structureless 'fluid mosaic' solvent is far from correct. The lipid bilayer displays distinct static and dynamic structural organization on a small scale, for example in terms of differentiated lipid domains, and evidence is accumulating that these structures are of importance for the functioning of biological membranes, including the activity of membrane-bound enzymes and receptors and morphological changes at the cell surface. Insight into the relationship between this small-scale structure and biological functioning holds promise for a more rational approach to modulate function via manipulation of the lipid microenvironment and the lipid/protein interface in particular. Computer simulation has proved to be a useful tool in investigating membrane structure on a small scale-specifically the nanometer scale (1-100 nm), which is in between the molecular scale accessible by various spectroscopic techniques and molecular dynamics calculations, and the micrometer scale accessible by scattering and microscopy techniques.

Studies on the lack of cooperativity in the melting of lipid bilayers.

The gel-to-fluid first-order melting transition of lipid bilayers is simulated by the use of a microscopic interaction model which includes a variable number of lipid-chain conformational states. The results suggest that the experimental observation of 'continuous melting' in pure wet lipid bilayers, rather than being ascribed to the presence of impurities, may be explained as a result of kinetically caused metastability of intermediate lipid-chain conformations.

Modelling the phase equilibria in two-component membranes of phospholipids with different acyl-chain lengths.

A phenomenological model is proposed to describe the membrane phase equilibria in binary mixtures of saturated phospholipids with different acyl-chain lengths. The model is formulated in terms of thermodynamic and thermomechanic properties of the pure lipid bilayers, specifically the chain-melting transition temperature and enthalpy, the hydrophobic bilayer thickness, and the lateral area compressibility modulus. The model is studied using a regular solution theory made up of a set of interaction parameters which directly identify that part of the lipid-lipid interaction which is due to hydrophobic mismatch of saturated chains of different lengths. It is then found that there is effectively a single universal interaction parameter which, in the full composition range, describes the phase equilibria in mixtures of DMPC/DPPC, DPPC/DSPC, DMPC/DSPC, and DLPC/DSPC, in excellent agreement with experimental measurements. The model is used to predict the variation with temperature and composition of the specific heat, as well as of the average membrane thickness and area in each of the phases. Given the value of the universal interaction parameter, the model is then used to predict the phase diagrams of binary mixtures of phospholipids with different polar head groups, e.g., DPPC/DPPE, DMPC/DPPE and DMPE/DSPC. By comparison with experimental results for these mixtures, it is shown that difference in acyl-chain lengths gives the major contribution to deviation from ideal mixing. Application of the model to mixtures with non-saturated lipids is also discussed.

Passive ion permeability of lipid membranes modelled via lipid-domain interfacial area.

A microscopic interaction model of the gel-to-fluid chain-melting phase transition of fully hydrated lipid bilayer membranes is used as a basis for modelling the temperature dependence of passive transmembrane permeability of small ions, e.g. Na+. Computer simulation of the model shows that the phase transition is accompanied by strong lateral density fluctuations which manifest themselves in the formation of inhomogeneous equilibrium structures of coexisting gel and fluid domains. The interfaces of these domains are found to be dominated by intermediate lipid-chain conformations. The interfacial area is shown to have a pronounced peak at the phase transition. By imposing a simple model for ion diffusion through membranes which assigns a high relative permeation rate to the domain interfaces, the interfacial area is then identified as a membrane property which has the proper temperature variation to account for the peculiar experimental observation of a strongly enhanced passive ion permeability at the phase transition. The excellent agreement with the experimental data for Na+-permeation, taken together with recent experimental results for the phase transition kinetics, provides new insight into the microphysical mechanisms of reversible electric breakdown. This insight indicates that there is no need for aqueous pore-formation to explain the experimental observation of a dramatic increase in ion conductance subsequent to electric pulses.

Lindane suppresses the lipid-bilayer permeability in the main transition region.

The effects of a small molecule, the insecticide lindane, on unilamellar DMPC bilayers in the phase transition region, have been studied by means of differential scanning calorimetry and fluorescence spectroscopy. The calorimetric data show that increasing concentrations of lindane broaden the transition and lower the transition temperature, without changing the transition enthalpy significantly. Lindane therefore enhances the thermal fluctuations of the bilayer. The calorimetric data furthermore suggest that the bilayer structure is intact and not disrupted by even high concentrations (32 mol%) of lindane. Fluorescence spectroscopy was used to measure the passive permeability of unilamellar DMPC bilayers to Co2+ ions. The data show that lindane seals the bilayer for Co2+ penetration and that this effect increases with increasing lindane concentration. The results are discussed in relation to the effects on the permeability of other small molecules, e.g., anesthetics.

Calorimetric and theoretical studies of the effects of lindane on lipid bilayers of different acyl chain length.

The effects of the insecticide lindane on the phase transition in multilamellar bilayers of saturated diacylphosphatidylcholines of different acyl chain length (DC14PC, DC16PC, and DC18PC) have been studied by means of differential scanning calorimetry (DSC), as well as computer-simulation calculations on a molecular interaction model. The calorimetric data show that increasing concentrations of lindane lower the transition temperature and lead to a broadening of the specific heat in a systematic way depending on the lipid acyl chain length. Kinetic effects in the observed calorimetric traces indicate that the incorporation of lindane into multilamellar lipid bilayers is slow, but faster for the shorter lipid species. Large unilamellar vesicles do not show such kinetic effects. The transition enthalpy is for all three lipid species found to be independent of the lindane concentration which implies that the entropy of mixing is vanishingly small. This lends support to a microscopic molecular interaction model which assigns the absorbed lindane molecules to interstitial sites in the bilayer. Computer-simulation calculations on this model, which assumes a specific interaction between lindane and certain excited acyl chain configurations, lead to predictions of the lipid-water partition coefficient in qualitative agreement with experimental measurements (Antunes-Madeira and Madeira (1985) Biochim. Biophys. Acta 820, 165-172). The partition coefficient has a peak near the phase transition which is a consequence of enhanced interfacial adsorption of lindane at lipid-domain interfaces.

Phospholipase A(2) activity towards vesicles of DPPC and DMPC-DSPC containing small amounts of SMPC.

Phospholipase A(2) (PLA(2)) is an interfacially active enzyme whose hydrolytic activity is known to be enhanced in one-component phospholipid bilayer substrates exhibiting dynamic micro-heterogeneity. In this study the activity of PLA(2) towards large unilamellar vesicles composed of DPPC:SMPC and DMPC:DSPC:SMPC is investigated using fluorescence and HPLC techniques. Phase diagrams of the mixtures are established by differential scanning calorimetry and the PLA(2) activity, monitored by the lag time, is correlated with the phase behavior of the mixtures. In addition, the degree of lipid hydrolysis in the DMPC:DSPC:SMPC lipid mixtures is detected by HPLC. The PLA(2) activity is found to be significantly increased in the temperature range of the coexistence region where the lipid mixtures exhibit lateral gel-fluid phase separation. Furthermore, in the entire temperature range it is demonstrated that PLA(2) preferentially hydrolyzes the short chain DMPC lipid. This discriminative effect becomes less pronounced when the asymmetric lipid SMPC is present in the lipid substrate. Inclusion of SMPC into either DPPC or DMPC:DSPC vesicles prolongs the lag time. The results clearly show that the PLA(2) activity is significantly enhanced by lipid bilayer micro-heterogeneity in both one-component and multi-component lipid bilayer substrates. The PLA(2) activity measurements are discussed in terms of dynamic gel-fluid lipid domain formation due to density fluctuations and static lipid domain formation due to gel-fluid phase separation.

Phase behavior and permeability properties of phospholipid bilayers containing a short-chain phospholipid permeability enhancer.

The thermodynamic phase behavior and trans-bilayer permeability properties of multilamellar phospholipid vesicles containing a short-chain DC10PC phospholipid permeability enhancer have been studied by means of differential scanning calorimetry and fluorescence spectroscopy. The calorimetric scans of DC14PC lipid bilayer vesicles incorporated with high concentrations of DC10PC demonstrate a distinct influence on the lipid bilayer thermodynamics manifested as a pronounced freezing-point depression and a narrow phase coexistence region. Increasing amounts of DC10PC lead to a progressive lowering of the melting enthalpy, implying a mixing behavior of the DC10PC in the bilayer matrix similar to that of a substitutional impurity. The phase behavior of the DC10PC-DC14PC mixture is supported by fluorescence polarization measurements which, furthermore, in the low-temperature gel phase reveal a non-monotonic concentration-dependent influence on the structural bilayer properties; small concentrations of DC10PC induce a disordering of the acyl chains, whereas higher concentrations lead to an ordering. Irreversible fluorescence quench measurements demonstrate a substantial increase in the trans-bilayer permeability over broad temperature and composition ranges. At temperatures corresponding to the peak positions of the heat capacity, a maximum in the trans-bilayer permeability is observed. The influence of DC10PC on the lipid bilayer thermodynamics and the associated permeability properties is discussed in terms of microscopic effects on the lateral lipid organization and heterogeneity of the bilayer.

The permeability and the effect of acyl-chain length for phospholipid bilayers containing cholesterol: theory and experiment.

The model of Cruzeiro-Hansson et al. (Biochim. Biophys. Acta (1989) 979, 166-1176) for lipid-cholesterol bilayers at low cholesterol concentrations is used to predict the thermodynamic properties and the passive ion permeability of lipid bilayers as a function of acyl-chain length and cholesterol concentration. Numerical simulations based on the Monte Carlo method are used to determine the equilibrium state of the system near the main gel-fluid phase transition. The permeability is calculated using an ansatz which relates the passive permeability to the amount of interfaces formed in the bilayer when cholesterol is present. The model predicts at low cholesterol contents an increase in the membrane permeability in the transition region both for increasing cholesterol concentration and for decreasing chain length at a given value of the reduced temperature. This is in contrast to the case of lipid bilayers containing high cholesterol concentrations where the cholesterol strongly suppresses the permeability. Experimental results for the Na+ permeability of C15PC and DPPC (C16PC) bilayers containing cholesterol are presented which confirm the theoretical predictions at low cholesterol concentrations.

Intrinsic molecules in lipid membranes change the lipid-domain interfacial area: cholesterol at domain interfaces.

A theoretical analysis of the effects of intrinsic molecules on the lateral density fluctuations in lipid bilayer membranes is carried out by means of computer simulations on a microscopic interaction model of the gel-to-fluid chain-melting phase transition. The inhomogeneous equilibrium structures of gel and fluid domains, which in previous work (Cruzeiro-Hansson, L. and Mouritsen, O.G. (1988) Biochim. Biophys. Acta 944, 63-72) were shown to characterize the transition region of pure lipid membranes, are here shown to be enhanced by intrinsic molecules such as cholesterol. Cholesterol is found to increase the interfacial area and to accumulate in the interfaces. The interfacial area, the average cluster size, the lateral compressibility, and the membrane area are calculated as functions of temperature and cholesterol concentration. It is shown that the enhancement by cholesterol of the lateral density fluctuations and the lipid-domain interfacial area is most pronounced away from the transition temperature. The implications of the results are discussed in relation to passive ion permeability and function of interfacially active enzymes such as phospholipase.

Association of ethanol with lipid membranes containing cholesterol, sphingomyelin and ganglioside: a titration calorimetry study.

The association of ethanol at physiologically relevant concentrations with lipid bilayers of different lipid composition has been investigated by use of isothermal titration calorimetry (ITC). The liposomes examined were composed of combinations of lipids commonly found in neural cell membranes: dimyristoyl phosphatidylcholine (DMPC), ganglioside (GM(1)), sphingomyelin and cholesterol. The calorimetric results show that the interaction of ethanol with fluid lipid bilayers is endothermic and strongly dependent on the lipid composition of the liposomes. The data have been used to estimate partitioning coefficients for ethanol into the fluid lipid bilayer phase and the results are discussed in terms of the thermodynamics of partitioning. The presence of 10 mol% sphingomyelin or ganglioside in DMPC liposomes enhances the partitioning coefficient by a factor of 3. Correspondingly, cholesterol (30 mol%) reduces the partitioning coefficient by a factor of 3. This connection between lipid composition and partitioning coefficient correlates with in vivo observations. Comparison of the data with the molecular structure of the lipid molecules suggests that ethanol partitioning is highly sensitive to changes in the lipid backbone (glycerol or ceramide) while it appears much less sensitive to the nature of the head group.

A microscopic model for lipid/protein bilayers with critical mixing.

A statistical mechanical lattice model is proposed to describe the phase diagram of phospholipid bilayers with small transmembrane proteins or polypeptides. The model is based on the extended Pink-Green-Chapman model (Zhang et al. (1992) Phys. Rev. A 45, 7560-7567) for pure lipid bilayers which undergo a first-order gel-fluid phase transition. The interaction between the lipid bilayer and the protein or polypeptide is modelled using the concept of hydrophobic matching. The phase diagram has been derived by computer-simulation techniques which fully account for thermal density fluctuations and which operate on the level of the free-energy thereby permitting an accurate identification of the phase boundaries. The calculations predict a closed loop of gel-fluid coexistence with a lower critical mixing point. Specific-heat traces across the phase diagram are also presented. The theoretical results for the phase diagram, the specific-heat function, and the transition enthalpy are related to recent experimental measurements on phospholipid bilayers mixed with synthetic transmembrane amphiphilic polypeptides or with gramicidin A.

Phase equilibria in the phosphatidylcholine-cholesterol system.

A thermodynamic and a microscopic interaction model are proposed to describe the phase equilibria in the phosphatidylcholine-cholesterol system. The model calculations allow for a solid phase with conformationally ordered acyl chains and liquid phases with conformationally ordered as well as disordered chains. The resulting phase diagram is in excellent agreement with the experimental phase diagram for dipalmitoylphosphatidylcholine bilayers with cholesterol as determined by a recent NMR and calorimetry study. It is thus demonstrated that the phase behaviour of phosphatidylcholine-cholesterol mixtures can be rationalized using only a few basic assumptions: (i) Cholesterol interacts favourably with phosphatidylcholine chains in an extended conformation, (ii) the main transition of pure phosphatidylcholine bilayers takes place in terms of translational variables as well acyl-chain conformational variables, and (iii) cholesterol disturbs the translational order in the crystalline (gel) state of phosphatidylcholine. These results suggest that the occurrence of specific phosphatidylcholine-cholesterol complexes is not implied by the experimental thermodynamic data.

The effects of density fluctuations on the partitioning of foreign molecules into lipid bilayers: application to anaesthetics and insecticides.

An extensive computer-simulation study is performed on a simple but general molecular model recently proposed (Jørgensen et al. (1991) Biochem. Biophys. Acta 1062, 277-238) to describe foreign molecules interacting with lipid bilayers. The model is a multi-state lattice model of the main bilayer transition in which the foreign molecules are assumed to intercalate at interstitial lattice positions. Specific as well as non-specific interactions between the foreign molecules and the lipid acyl chains are considered. Particular attention is paid to the fluctuating properties of the membrane and how the presence of the foreign molecules modulates these fluctuations in the transition region. By means of computer-stimulation techniques, a detailed account is given of the macroscopic as well as microscopic consequences of the fluctuations. The macroscopic consequences of the fluctuations are seen in the thermal anomalies of the specific heat and the passive trans-membrane permeability. Microscopically, the fluctuations manifest themselves in lipid-domain formation in the transition region which implies an effective dynamic membrane heterogeneity. Within the model it is found that certain anaesthetics and insecticides which are characterised by specific interactions with the lipids have a strong effect on the heterogeneity of the membrane inducing regions of locally very high concentration of the foreign molecules. This leads to a broadening of the specific heat peak and a maximum in the membrane/water partition coefficient. These results are in accordance with available experimental data for volatile general anaesthetics like halothane, local anaesthetics like cocain derivatives, and insecticides like lindane.

A general model for the interaction of foreign molecules with lipid membranes: drugs and anaesthetics.

A general microscopic interaction model is proposed to describe the changes in the physical properties of phospholipid bilayer membranes due to foreign molecules which, to different degrees, partition between the membrane phases and the aqueous environment. The model is a multi-state lattice model for the main phase transition of lipid bilayers and the foreign molecules are assumed to intercalate as interstitials in the lattice. By varying the model parameters, the diversity in the thermodynamic properties of the model is explored using computer-simulation techniques which faithfully take account of the thermal fluctuations. The calculations are performed in both the canonical and the grand canonical ensembles corresponding to the cases where the concentration of foreign molecules in the membrane is either fixed or varies as the external conditions are changed. A classification of the diverse thermal behaviour, specifically with regard to the phase diagram, the specific heat, the density fluctuations, and the partition coefficient, is suggested with a view to rationalizing a large body of experimental measurements of the effects of different foreign molecules on membrane properties. The range of foreign molecules considered includes compounds as diverse as volatile general anaesthetics like halothane, cocaine-derived local anaesthetics like procaine, calcium-channel blocking drugs like verapamil, antidepressants like chlorpromazine, and anti-cancer agents like adriamycin.

Theoretical analysis of protein organization in lipid membranes.

The fundamental physical principles of the lateral organization of trans-membrane proteins and peptides as well as peripheral membrane proteins and enzymes are considered from the point of view of the lipid-bilayer membrane, its structure, dynamics, and cooperative phenomena. Based on a variety of theoretical considerations and model calculations, the nature of lipid-protein interactions is considered both for a single protein and an assembly of proteins that can lead to aggregation and protein crystallization in the plane of the membrane. Phenomena discussed include lipid sorting and selectivity at protein surfaces, protein-lipid phase equilibria, lipid-mediated protein-protein interactions, wetting and capillary condensation as means of protein organization, mechanisms of two-dimensional protein crystallization, as well as non-equilibrium organization of active proteins in membranes. The theoretical findings are compared with a variety of experimental data.

Model of a sub-main transition in phospholipid bilayers.

A recently discovered submain phase transition in multi-lamellar bilayers of long-chain saturated diacyl phosphatidylcholines (Jørgensen, K. (1995) Biochim. Biophys. Acta 1240, 111-114) is discussed in terms of a theoretical molecular interaction model using computer simulation techniques. The model interprets the transition to be due to a decoupling of the acyl-chain melting from the melting of the pseudo-two-dimensional crystalline lattice of the P beta' phase. A two-stage melting process is predicted by the calculations suggesting that the sub-main transition involves a lattice melting whereas the acyl-chain melting takes place at a higher temperature at the main transition. The calculated heat contents of the two transitions as well as the chain-length dependence compare favorably with experimental data for multi-lamellar phosphatidylcholine lipid bilayers.

Phase equilibria and local structure in binary lipid bilayers.

A molecular interaction model is used to describe the phase diagram of two-component phospholipid bilayer membranes of saturated phospholipids, DCnPC, with different acyl-chain lengths, n = 12,14,18,20. The interaction between acyl chains of different length is formulated in terms of a hydrophobic mismatch which permits the series of binary phase diagrams to be calculated in terms of a single 'universal' interaction parameter. The properties of the model are calculated by computer-simulation techniques which not only permit determination of the specific-heat function and the phase diagram but also reveal the local structure of the mixture in the different parts of the phase diagram. The local structure is described pictorially and characterized quantitatively in terms of a correlation function. It is shown that the non-ideal mixing of lipid species due to mismatch in the hydrophobic lengths leads to a progressively increasing local ordering as the chain-length difference is increased. A pronounced local structure is found to persist deep inside the fluid phase of the mixture. The local structure is discussed in relation to the features observed in the specific-heat function, for which theoretical data, as well as experimental data obtained from differential-scanning calorimetry are presented.

Acyl-coenzyme A organizes laterally in membranes and is recognized specifically by acyl-coenzyme A binding protein.

Long chain acyl-coenzyme A (acyl-CoA) is a biochemically important amphiphilic molecule that is known to partition strongly into membranes by insertion of the acyl chain. At present, microscopically resolved evidence is lacking on how acyl-CoA influences and organizes laterally in membranes. By atomic force microscopy (AFM) imaging of membranes exposed to acyl-CoA in microM concentrations, it is shown that aggregate formation takes place within the membrane upon long-time exposure. It is known that acyl-CoA is bound by acyl-CoA binding protein (ACBP) with high affinity and specificity and that ACBP may bind and desorb membrane-bound acyl-CoA via a partly unknown mechanism. Following incubation with acyl-CoA, it is shown that ACBP is able to reverse the formation of acyl-CoA aggregates and to associate peripherally with acyl-CoA on the membrane surface. Our microscopic results point to the role of ACBP as an intermembrane transporter of acyl-CoA and demonstrate the ability of AFM to reveal the remodelling of membranes by surfactants and proteins.

Anesthetics affect membrane heterogeneity.

A computer-simulation study is performed on a simple but general molecular model of the interaction of general and local anesthetics with lipid membranes. In the neighborhood of the gel-to-fluid membrane phase transition it is found that anesthetics have a strong effect on the heterogeneity of the membrane-inducing regions of locally high anesthetics concentration. This leads to a broadening of the specific heat peak and a maximum in the membrane/water partition coefficient. These results are in accordance with available experimental data for volatile general anesthetics like halothane and local anesthetics like cocaine derivatives.