Molecular simulation of dioleoylphosphatidylcholine lipid bilayers at differing levels of hydration.

The structure and dynamics of the lipid and water components of dioleoylphosphatidylcholine bilayers at various levels of hydration were studied using molecular dynamics (MD) simulations. Equilibration of these systems proceeded by use of a hybrid MD and configurational-bias Monte Carlo technique using one atmosphere of pressure normal to the membrane and a set point for the lateral area derived from experimental Bragg spacings, combined with experimentally derived specific volumes for each of the system components. Membrane surface tensions were observed to be of the order of tens of dyn/cm. The transbilayer molecular fragment peak positions at low hydration were found to agree with experimental neutron and x-ray scattering profiles and previously published simulations. For hydration levels of 5.4, 11.4, and 16 waters/lipid, molecular fragment distributions and order parameters for the headgroup, lipid chains, and water were quantified. Spin-lattice relaxation rates and lateral self-diffusion coefficients of water agreed well with results from experimental nuclear magnetic resonance studies. Relaxation rates of the choline segments and chemical shift anisotropies for the phosphate and carbonyls were computed. Headgroup orientation, as measured by the P-N vector, showed enhanced alignment with the membrane surface at low hydration. The sign of the membrane dipole potential reversed at low hydration, with the membrane interior negative relative to the interlamellar region. Calculation of the number of water molecules in the headgroup hydration shell, as a function of hydration level, supports the hypothesis that the break point in the curve of Bragg spacing versus hydration level near 12 waters/lipid, observed experimentally by Hristova and White (1988. Biophys. J. 74:2419-2433), marks the completion of the first hydration shell.     

Structure and interactions of fully hydrated dioleoylphosphatidylcholine bilayers.

This study focuses on dioleoylphosphatidylcholine (DOPC) bilayers near full hydration. Volumetric data and high-resolution synchrotron x-ray data are used in a method that compares DOPC with well determined gel phase dipalmitoylphosphatidylcholine (DPPC). The key structural quantity obtained is fully hydrated area/lipid A0 = 72.2 +/- 1.1 A2 at 30 degrees C, from which other quantities such as thickness of the bilayer are obtained. Data for samples over osmotic pressures from 0 to 56 atmospheres give an estimate for the area compressibility of KA = 188 dyn/cm. Obtaining the continuous scattering transform and electron density profiles requires correction for liquid crystal fluctuations. Quantitation of these fluctuations opens an experimental window on the fluctuation pressure, the primary repulsive interaction near full hydration. The fluctuation pressure decays exponentially with water spacing, in agreement with analytical results for soft confinement. However, the ratio of decay length lambda(fl) = 5.8 A to hydration pressure decay length lambda = 2.2 A is significantly larger than the value of 2 predicted by analytical theory and close to the ratio obtained in recent simulations. We also obtain the traditional osmotic pressure versus water spacing data. Our analysis of these data shows that estimates of the Hamaker parameter H and the bending modulus Kc are strongly coupled.     

Lateral clustering of lipids in hydrated bilayers composed of dioleoylphosphatidylcholine and dipalmitoylphosphatidylcholine

Investigation of lateral heterogeneities (clusters) in cell membranes is an important step toward understanding the physical processes that lead to the formation of lipid domains and rafts. Computer modeling methods represent a powerful tool to solve the problem, since they can detect clusters containing only a few lipid molecules—the situation that still resists characterization with modern experimental techniques. Parameters of clustering depend on lipid composition of a membrane. In this work, we propose a computational method to detect and analyze parts of membrane with different packing densities. Series of one- and two-component fluid systems containing lipids with the same polar heads and different acyl chains, dioleoylphosphatidylcholine (18 : 1) and dipalmitoylphosphatidylcholine (16 : 0), were chosen as the objects under study. The developed algorithm is based on molecular dynamics simulation of hydrated lipid bilayers in all-atom mode. The method is universal and could be applied to any other membrane system with arbitrary lipid composition. Here, we demonstrated that the studied lipid bilayers reveal small lateral dynamic clusters composed of just several (most often, three) lipid molecules. This seems to be one of the most important reasons determining the “mosaic” nature of the membrane-water interface.     

Lateral interaction of cholesterol in diacylphosphatidylcholesterol bilayers

Thermotropic phase-transition properties of the aqueous dispersions of several diacylphosphatidylcholesterol (DRCh) analogs are examined. The aqueous dispersions of their calcium salts exhibit characteristic endothermic thermotropic transitions due to a change in the conformation of acyl chains. These dispersions consist of osmotically intact liposomes that trap ions, and at the transition temperature there is anomalous increase in the ion leakage. Wide-angle electron diffraction studies of DPCh · Ca monolayers also exhibit a transition from a sharp 4.25 Å band to a broad one centering at 4.7 Å, reflecting an order-disorder transition in the acyl chains. The long-range order in the organization of acyl chains of DRCh molecules could arise from intermolecular interactions between the cholesterol moleties to form a functional dimer, and such dimers are apparently cross-linked by Ca²⁺ to form a long-range interacting lattice of acyl chains. Evidence for this model is adduced from the fluorescence properties of the dispersions of dimyristoylphosphatidylcholesta-5,7,9-trienol. The phase-transition properties of DRCh are an ideal illustration of calcium-induced isothermal phase transition.     

Helix-helix interactions in lipid bilayers.

Using a continuum model, we calculated the electrostatic interaction free energy between two alpha-helices in three environments: the aqueous phase, a low dielectric alkane phase, and a simple representation of a lipid bilayer. As was found in previous work, helix-helix interactions in the aqueous phase are quite weak, because of solvent screening, and slightly repulsive, because of desolvation effects that accompany helix assembly. In contrast, the interactions can be quite strong in a hypothetical alkane phase because desolvation effects are essentially nonexistent and because helix-helix interactions are not well screened. In this type of environment, the antiparallel helix orientation is strongly favored over the parallel orientation. In previous work we found that the free energy penalty associated with burying helix termini in a bilayer is quite high, which is why the termini tend to protrude into the solvent. Under these conditions the electrostatic interaction is strongly screened by solvent; indeed, it is sufficient for the termini to protrude a few angstroms from the two surfaces of the bilayer for their interaction to diminish almost completely. The effect is consistent with the classical model of the helix dipole in which the dipole moment is represented by point charges located at either terminus. Our results suggest, in agreement with previous models, that there is no significant nonspecific driving force for helix aggregation and, hence, that membrane protein folding must be driven by specific interactions such as close packing and salt-bridge and hydrogen bond formation.     

Lateral diffusion in the liquid phases of dimyristoylphosphatidylcholine/cholesterol lipid bilayers: a free volume analysis.

The technique of fluorescence recovery after photobleaching is used to perform an extensive study of the lateral diffusion of a phospholipid probe in the binary mixture dimyristoylphosphatidylcholine/cholesterol, above the melting temperature of the phospholipid. In the regions of the phase diagram where a single liquid phase exists, diffusion can be quantitatively described by free volume theory, using a modified Macedo-Litovitz hybrid equation. In the liquid-liquid immiscibility region, the temperature dependence of the diffusion coefficient is in excellent agreement with current theories of generalized diffusivities in composite two-phase media. A consistent interpretation of the diffusion data can be provided based essentially on the idea that the primary effect of cholesterol addition to the bilayer is to occupy free volume. On this basis, a general interpretation of the phase behavior of this mixture is also proposed.     

Role of cholesterol flip-flop in oxidized lipid bilayers

We performed a series of molecular dynamics simulations of cholesterol (Chol) in nonoxidized 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) bilayer and in binary mixtures of PLPC-oxidized-lipid-bilayers with 0–50% Chol concentration and oxidized lipids with hydroperoxide and aldehyde oxidized functional groups. From the 60 unbiased molecular dynamics simulations (total of 161 μs), we found that Chol inhibited pore formation in the aldehyde-containing oxidized lipid bilayers at concentrations greater than 11%. For both pure PLPC bilayer and bilayers with hydroperoxide lipids, no pores were observed at any Chol concentration. Furthermore, increasing cholesterol concentration led to a change of phase state from the liquid-disordered to the liquid-ordered phase. This condensing effect of Chol was observed in all systems. Data analysis shows that the addition of Chol results in an increase in bilayer thickness. Interestingly, we observed Chol flip-flop only in the aldehyde-containing lipid bilayer but neither in the PLPC nor the hydroperoxide bilayers. Umbrella-sampling simulations were performed to calculate the translocation free energies and the Chol flip-flop rates. The results show that Chol’s flip-flop rate depends on the lipid bilayer type, and the highest rate are found in aldehyde bilayers. As the main finding, we shown that Chol stabilizes the oxidized lipid bilayer by confining the distribution of the oxidized functional groups.     

The condensing effect of cholesterol in lipid bilayers

The condensing effect of cholesterol on phospholipid bilayers was systematically investigated for saturated and unsaturated chains, as a function of cholesterol concentration. X-ray lamellar diffraction was used to measure the phosphate-to-phosphate distances, PtP, across the bilayers. The measured PtP increases nonlinearly with the cholesterol concentration until it reaches a maximum. With further increase of cholesterol concentration, the PtP remains at the maximum level until the cholesterol content reaches the solubility limit. The data in all cases can be quantitatively explained with a simple model that cholesterol forms complexes with phospholipids in the bilayers. The phospholipid molecules complexed with cholesterol are lengthened and this lengthening effect extends into the uncomplexed phospholipids surrounding the cholesterol complexes. This long-range thickening effect is similar to the effect of gramicidin on the thickness of lipid bilayers due to hydrophobic matching.     

Temperature dependence of structure, bending rigidity, and bilayer interactions of dioleoylphosphatidylcholine bilayers

X-ray diffuse scattering was measured from oriented stacks and unilamellar vesicles of dioleoylphosphatidylcholine lipid bilayers to obtain the temperature dependence of the structure and of the material properties. The area/molecule, A, was 75.5 Å² at 45°C, 72.4 Å² at 30°C, and 69.1 Å² at 15°C, which gives the area expansivity α_A = 0.0029/deg at 30°C, and we show that this value is in excellent agreement with the polymer brush theory. The bilayer becomes thinner with increasing temperature; the contractivity of the hydrocarbon portion was α_Dc = 0.0019/deg; the difference between α_A and α_Dc is consistent with the previously measured volume expansivity α_Vc = 0.0010/deg. The bending modulus K_C decreased as exp(455/T) with increasing T (K). Our area compressibility modulus K_A decreased with increasing temperature by 5%, the same as the surface tension of dodecane/water, in agreement again with the polymer brush theory. Regarding interactions between bilayers, the compression modulus B as a function of interbilayer water spacing D′_W was found to be nearly independent of temperature. The repulsive fluctuation pressure calculated from B and K_C increased with temperature, and the Hamaker parameter for the van der Waals interaction was nearly independent of temperature; this explains why the fully hydrated water spacing, D′_W, that we obtain from our structural results increases with temperature.     

Thermodynamics of lipid interactions in complex bilayers

The mutual interactions between lipids in bilayers are reviewed, including mixtures of phospholipids, and mixtures of phospholipids and cholesterol (Chol). Binary mixtures and ternary mixtures are considered, with special emphasis on membranes containing Chol, an ordered phospholipid, and a disordered phospholipid. Typically the ordered phospholipid is a sphingomyelin (SM) or a long-chain saturated phosphatidylcholine (PC), both of which have high phase transitions temperatures; the disordered phospholipid is 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) or dioleoylphosphatidylcholine (DOPC). The unlike nearest-neighbor interaction free energies (ω_AB) between lipids (including Chol), obtained by an variety of unrelated methods, are typically in the range of 0-400 cal/mol in absolute value. Most are positive, meaning that the interaction is unfavorable, but some are negative, meaning it is favorable. It is of special interest that favorable interactions occur mainly between ordered phospholipids and Chol. The interpretation of domain formation in complex mixtures of Chol and phospholipids in terms of phase separation or condensed complexes is discussed in the light of the values of lipid mutual interactions.     

Transbilayer distribution of bromine in fluid bilayers containing a specifically brominated analogue of dioleoylphosphatidylcholine

We describe in this paper the transbilayer distribution of the bromines of the specifically halogenated phospholipid 1-oleoyl-2-(9,10-dibromostearoyl)-sn-glycero-3- phosphocholine (OBPC). The distribution was determined by X-ray diffraction of oriented multilayers of mixtures of OBPC and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) at 66% relative humidity by the general approach of Franks et al. (1978) [Nature 276, 530-532]. The bromine distribution of OBPC in the fluid L alpha phase is described accurately by a pair of Gaussian functions located 7.97 +/- 0.27 A from the center of the bilayer with l/e half-widths of 4.96 +/- 0.62 A. We find that OBPC bilayers are accurately described as DOPC bilayers with an additional bromine distribution centered at the position of the double bond of DOPC and conclude that OBPC is an excellent structural isomorph for DOPC under the conditions of these experiments. The distribution obtained is the complete and fully resolved transbilayer image of the halogen label because the broad distribution of the bromines is due entirely to thermal disorder and not to experimental limitations [Wiener, M. C., & White, S. H. (1991a) Biophys. J. 59, 162-173]. The observed width of the bromine distribution indicates that virtually all of the hydrocarbon interior is accessible to the bromines. The distance between the bromine/double-bond position and the headgroup phosphate position was determined from one-dimensional Patterson maps and found to be approximately 12 A. The application of accurately determined bromine distributions to the quantitative interpretation of fluorescence quenching experiments is discussed. A method for the self-consistent global analysis of diffraction data from mixtures that permits the use of data sets with different instrumental scale factors is developed in an Appendix.     

Pressure effects on the lateral distribution of cholesterol in lipid bilayers: a time-resolved spectroscopy study.

The effects of hydrostatic pressure and temperature on the phase behavior and physical properties of the binary mixture palmitoyloleoylphosphatidylcholine/cholesterol, over the 0-40 molar % range of cholesterol compositions, were determined from the changes in the fluorescence lifetime distribution and anisotropy decay parameters of the natural lipid trans-parinaric acid (t-PnA). Pressurized samples were excited with a Ti-sapphire subpicosecond laser, and fluorescence decays were analyzed by the quantified maximum entropy method. Above the transition temperature (T(T) = -5 degrees C), at atmospheric pressure, two liquid-crystalline phases, alpha and beta, are formed in this system. At each temperature and cholesterol concentration below the transition pressure, the fluorescence lifetime distribution pattern of t-PnA was clearly modulated by the pressure changes. Pressure increased the fraction of the liquid-ordered beta-phase and its order parameter, but it decreased the amount of cholesterol in this phase. Palmitoyloleoylphosphatidylcholine/cholesterol phase diagrams were also determined as a function of temperature and hydrostatic pressure.     

Quantitation of cholesterol incorporation into extruded lipid bilayers

Cholesterol incorporation into lipid bilayers, in the form of multilamellar vesicles or extruded large unilamellar vesicles, has been quantitated. To this aim, the cholesterol contents of bilayers prepared from phospholipid:cholesterol mixtures 33-75 mol% cholesterol have been measured and compared with the original mixture before lipid hydration. There is a great diversity of cases, but under most conditions the actual cholesterol proportion present in the extruded bilayers is much lower than predicted. A quantitative analysis of the vesicles is thus required before any experimental study is undertaken.     

Effects of cholesterol and PIP2 on interactions between glycophorin A and Band 3 in lipid bilayers

In the erythrocyte membrane, the interactions between glycophorin A (GPA) and Band 3 are associated strongly with the biological function of the membrane and several blood disorders. In this work, using coarse-grained molecular-dynamics simulations, we systematically investigate the effects of cholesterol and phosphatidylinositol-4,5-bisphosphate (PIP2) on the interactions of GPA with Band 3 in the model erythrocyte membranes. We examine the dynamics of the interactions of GPA with Band 3 in different lipid bilayers on the microsecond time scale and calculate the binding free energy between GPA and Band 3. The results indicate that cholesterols thermodynamically favor the binding of GPA to Band 3 by increasing the thickness of the lipid bilayer and by producing an effective attraction between the proteins due to the depletion effect. Cholesterols also slow the kinetics of the binding of GPA to Band 3 by reducing the lateral mobility of the lipids and proteins and may influence the binding sites between the proteins. The anionic PIP2 lipids prefer binding to the surface of the proteins through electrostatic attraction between the PIP2 headgroup and the positively charged residues on the protein surface. Ions in the solvent facilitate PIP2 aggregation, which promotes the binding of GPA to Band 3.     

Phosphatidylcholine-cholesterol interactions: bilayers of heteroacid lipids containing linoleate lose calorimetric transitions at low cholesterol concentration

Model membranes composed of cholesterol plus one of two phosphatidylcholines (PC), each containing a saturated and a dienoic acyl chain, have been studied by differential scanning calorimetry. The gel to liquid-crystalline phase transition temperature of 1-palmitoyl-2-linoleoyl PC was -19.5 degrees C and that of 1-stearoyl-2-linoleoyl PC was -13.7 degrees C. The addition of cholesterol to the phosphatidylcholines in aqueous dispersion resulted in the progressive removal of the phase transition as observed by differential scanning calorimetry. Per mole of sterol in the membrane, cholesterol was more effective at reducing the enthalpy change of the phase transitions of these bilayers containing dienoic phosphatidylcholines than it is in eliminating the transition of membranes made with other phospholipids that contain more saturated chains. No transitions in membranes made with palmitoyl-linoleoyl PC or stearoyl-linoleoyl PC could be detected calorimetrically when 17 mol% cholesterol was present.     

Lateral diffusivity of lipid analogue excimeric probes in dimyristoylphosphatidylcholine bilayers.

The lateral mobility of pyrenyl phospholipid probes in dimyristoylphosphatidylcholine (DMPC) vesicles was determined from the dependence of the pyrene monomeric and excimeric fluorescence yields on the molar probe ratio. The analysis of the experimental data makes use of the milling crowd model for two-dimensional diffusivity and the computer simulated random walks of probes in an array of lipids. The fluorescence yields for 1-palmitoyl-2-(1'-pyrenedecanoyl)phosphatidylcholine (py10PC) in DMPC bilayers are well fitted by the model both below and above the fluid-gel phase transition temperature (Tc) and permit the evaluation of the probe diffusion rate (f), which is the frequency with which probes take random steps of length L, the host membrane lipid-lipid spacing. The lateral diffusion coefficient is then obtained from the relationship D = fL2/4. In passing through the fluid-gel phase transition of DMPC (Tc = 24 degrees C), the lateral mobility of py10PC determined in this way decrease only moderately, while D measured by fluorescence photobleaching recovery (FPR) experiments is lowered by two or more orders of magnitude in gel phase. This difference in gel phase diffusivities is discussed and considered to be related either to (a) the diffusion length in FPR experiments being about a micrometer or over 100 times greater than that of excimeric probes (approximately 1 nm), or (b) to nonrandomicity in the distribution of the pyrenyl probes in gel phase DMPC. At 35 degrees C, in fluid DMPC vesicles, the diffusion rate is f = 1.8 x 10(8) s-1, corresponding to D = 29 microns2 s-1, which is about three times larger than the value obtained in FPR experiments. The activation energy for lateral diffusion in fluid DMPC was determined to be 8.0 kcal/mol.     

Physicochemical interactions of amyloid beta-peptide with lipid bilayers

The aggregation and deposition onto neuronal cells of amyloid beta-peptide (Abeta) is central to the pathogenesis of Alzheimer's disease. Accumulating evidence suggests that membranes play a catalytic role in the aggregation of Abeta. This article summarizes the structures and properties of Abeta in solution and the physicochemical interaction of Abeta with lipid bilayers of various compositions. Reasons for discrepancies between results by different research groups are discussed. The importance of ganglioside clusters in the aggregation of Abeta is emphasized. Finally, a hypothetical physicochemical cascade in the pathogenesis of the disease is proposed.     

Hydration, structure, and molecular interactions in the headgroup region of dioleoylphosphatidylcholine bilayers: an electron spin resonance study

The relationship between bilayer hydration and the dynamic structure of headgroups and interbilayer water in multilamellar vesicles is investigated by electron spin resonance methods. Temperature variations of the order parameter of a headgroup spin label DPP-Tempo in DOPC in excess water and partially dehydrated (10 wt % water) show a cusp-like pattern around the main phase transition, Tc. This pattern is similar to those of temperature variations of the quadrupolar splitting of interbilayer D2O in PC and PE bilayers previously measured by 2H NMR, indicating that the ordering of the headgroup and the interbilayer water are correlated. The cusp-like pattern of these and other physical properties around Tc are suggestive of quasicritical fluctuations. Also, an increase (a decrease) in ordering of DPP-Tempo is correlated with water moving out of (into) interbilayer region into (from) the bulk water phase near the freezing point, Tf. Addition of cholesterol lowers Tf, which remains the point of increasing headgroup ordering. Using the small water-soluble spin probe 4-PT, it is shown that the ordering of interbilayer water increases with bilayer dehydration. It is suggested that increased ordering in the interbilayer region, implying a lowering of entropy, will itself lead to further dehydration of the interbilayer region until its lowered pressure resists further flow, i.e., an osmotic phenomenon.     

The effect of temperature on lipid-n-alkane interactions in lipid bilayers

The relationship between age-related alterations in the lipid composition of cultured rat-heart fibroblasts and several biochemical and biophysical parameters was investigated. Aged (14-15-day-old) cultures displayed higher mole ratios of sphingomyelin to phosphatidylcholine, as well as elevated cholesterol levels. A concomitant increase was observed in the total protein content of the cells and in the Vmax values of both membranal and cytoplasmic marker enzymes. Fluorescence photobleaching recovery was employed to study the lateral mobility of the lipid probe NBD-phosphatidylethanolamine and of membrane glycoproteins that bind succinylated concanavalin A. The mobile fractions of both probes were higher in aged cultures, while the lateral diffusion coefficients were lower. To further demonstrate the dependence of the above parameters on the cellular lipid composition, we have manipulated the lipid composition of old cultures by treatments with liposomes (small unilamellar vesicles) of specific compositions. Treatments which reversed the lipid composition towards that of young (5-6-day-old) cultures caused a concomitant reversal of the measured biochemical and biophysical parameters to the values observed in young cultures. These findings suggest that alterations in the organization and mobility of cell membrane constituents are involved in mediating changes in cellular functions. In view of our previous findings on cultures of rat-heart myocytes (Yechiel, E., Barenholz, Y. and Henis, Y.I. (1985) J. Biol. Chem. 260, 9132-9136), it appears that the modulation of cellular properties through the membrane lipid composition may be a general phenomenon in many cell types.     

NMR investigations of interactions between anesthetics and lipid bilayers

Interactions between anesthetics (lidocaine and short chain alcohols) and lipid membranes formed by dimyristoylphosphatidylcholine (DMPC) were studied using NMR spectroscopy. The orientational order of lidocaine was investigated using deuterium NMR on a selectively labelled compound whereas segmental ordering in the lipids was probed by two-dimensional ¹H-¹³C separated local field experiments under magic-angle spinning conditions. In addition, trajectories generated in molecular dynamics (MD) computer simulations were used for interpretation of the experimental results. Separate simulations were carried out with charged and uncharged lidocaine molecules. Reasonable agreement between experimental dipolar interactions and the calculated counterparts was observed. Our results clearly show that charged lidocaine affects significantly the lipid headgroup. In particular the ordering of the lipids is increased accompanied by drastic changes in the orientation of the P-N vector in the choline group.