Molecular dynamics simulations of phospholipid bilayers with cholesterol

To investigate the microscopic interactions between cholesterol and lipids in biological membranes, we have performed a series of molecular dynamics simulations of large membranes with different levels of cholesterol content. The simulations extend to 10 ns, and were performed with hydrated dipalmitoylphosphatidylcholine (DPPC) bilayers. The bilayers contain 1024 lipids of which 0-40% were cholesterol and the rest DPPC. The effects of cholesterol on the structure and mesoscopic dynamics of the bilayer were monitored as a function of cholesterol concentration. The main effects observed are a significant ordering of the DPPC chains (as monitored by NMR type order parameters), a reduced fraction of gauche bonds, a reduced surface area per lipid, less undulations--corresponding to an increased bending modulus for the membrane, smaller area fluctuations, and a reduced lateral diffusion of DPPC-lipids as well as cholesterols.     

Effect of cholesterol on the lactosylceramide domains in phospholipid bilayers

Lactosylceramide (LacCer) in the plasma membranes of immune cells is an important lipid for signaling in innate immunity through the formation of LacCer-rich domains together with cholesterol (Cho). However, the properties of the LacCer domains formed in multicomponent membranes remain unclear. In this study, we examined the properties of the LacCer domains formed in Cho-containing 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC) membranes by deuterium solid-state NMR and fluorescence lifetimes. The potent affinity of LacCer-LacCer (homophilic interaction) is known to induce a thermally stable gel phase in the unitary LacCer bilayer. In LacCer/Cho binary membranes, Cho gradually destabilized the LacCer gel phase to form the liquid-ordered phase by its potent order effect. In the LacCer/POPC binary systems without Cho, the ²H NMR spectra of 10′,10′-d₂-LacCer and 18′,18′,18′-d₃-LacCer probes revealed that LacCer was poorly miscible with POPC in the membranes and formed stable gel phases without being distributed in the liquid crystalline domain. The lamellar structure of the LacCer/POPC membrane was gradually disrupted at around 60°C, whereas the addition of Cho increased the thermal stability of the lamellarity. Furthermore, the area of the LacCer gel phase and its chain order were decreased in the LacCer/POPC/Cho ternary membranes, whereas the liquid-ordered domain, which was observed in the LacCer/Cho binary membrane, was not observed. Cho surrounding the LacCer gel domain liberated LacCer and facilitated forming the submicron to nano-scale small domains in the liquid crystalline domain of the LacCer/POPC/Cho membranes, as revealed by the fluorescence lifetimes of trans-parinaric acid and trans-parinaric acid-LacCer. Our findings on the membrane properties of the LacCer domains, particularly in the presence of Cho, would help elucidate the properties of the LacCer domains in biological membranes.     

The molecular packing and stability within highly curved phospholipid bilayers

It is shown that the area occupied per phospholipid molecule and the thickness of the bilayer are the same in vesicles as in a planar bilayer. From this it is concluded that the lower limit to the size of a vesicle depends on the packing of the head groups of the inner monolayer.     

Effect of free cholesterol on incorporation of triolein in phospholipid bilayers

Triacylglycerols are the major substrates for lipolytic enzymes that act at the surface of emulsion-like particles such as triglyceride-rich lipoproteins, chylomicrons, and intracellular lipid droplets. This study examines the effect of cholesterol on the solubility of a triacylglycerol, triolein, in phospholipid surfaces. Solubilities of [carbonyl-13C]triolein in phospholipid bilayer vesicles containing between 0 and 50 mol % free cholesterol, prepared by cosonication, were measured by 13C NMR. The carbonyl resonances from bilayer-incorporated triglyceride were shifted downfield in the 13C NMR spectra from those corresponding to excess, nonincorporated material. This enabled solubilities to be determined directly from carbonyl peak intensities at most cholesterol concentrations. The bilayer solubility of triolein was inversely proportional to the cholesterol/phospholipid mole ratio. In pure phospholipid vesicles the triolein solubility was 2.2 mol %. The triglyceride incorporation decreased to 1.1 mol % at a cholesterol/phospholipid mole ratio of 0.5, and at a mole ratio of 1.0 for the bilayer lipids, the triolein solubility was reduced to just 0.15 mol %. The effects of free cholesterol were more pronounced and progressive than observed previously on the bilayer solubility of cholesteryl oleate (Spooner, P. J. R., Hamilton, J. A., Gantz, D. L., & Small, D. M. (1986) Biochim. Biophys. Acta 860, 345-353]. As with cholesteryl oleate, we suggest that cholesterol also displaces solubilized triglyceride to deeper regions of the bilayer.     

Properties of unsaturated phospholipid bilayers: Effect of cholesterol

Properties of hydrated unsaturated phosphatidylcholine (PC) lipid bilayers containing 40 mol % cholesterol and of pure PC bilayers have been studied. Various methods were applied, including molecular dynamics simulations, self-consistent field calculations, and the pulsed field gradient nuclear magnetic resonance technique. Lipid bilayers were composed of 18:0/18:1(n-9)cis PC, 18:0/18:2(n-6)cis PC, 18:0/18:3(n-3)cis PC, 18:0/20:4(n-6)cis PC, and 18:0/22:6(n-3)cis PC molecules. Lateral self-diffusion coefficients of the lipids in all these bilayers, mass density distributions of atoms and atom groups with respect to the bilayer normal, the C-H and C-C bond order parameter profiles of each phospholipid hydrocarbon chain with respect to the bilayer normal were calculated. It was shown that the lateral self-diffusion coefficient of PC molecules of the lipid bilayer containing 40 mol % cholesterol is smaller than that for a corresponding pure PC bilayer; the diffusion coefficients increase with increasing the degree of unsaturation of one of the PC chains in bilayers of both types (i.e., in pure bilayers or in bilayers with cholesterol). The presence of cholesterol in a bilayer promoted the extension of saturated and polyunsaturated lipid chains. The condensing effect of cholesterol on the order parameters was more pronounced for the double C=C bonds of polyunsaturated chains than for single C-C bonds of saturated chains.     

EPR studies of phospholipid bilayers after lipoperoxidation. 1. Inner molecular order and fluidity gradient

The molecular order of fatty acyl chains in oriented lipid bilayers on solid support (SPB), made of either natural or synthetic phospholipids oxidized by Fenton reagent and probed with spin labeled lecithin (5-DSPC) was studied by means of EPR spectrometry. Phospholipids (ASPC, EYPC, mitochondrial extract) were oxidized as either aqueous buffer/methanol dispersions or reverse-phase evaporation vesicles (REV) suspensions. Oxidation was preliminarily revealed both by assaying MDA and by detecting conjugated dienes. Oxidized phospholipid species was quantified by preparative TLC. The degree of order in oriented lipid bilayers of samples containing oxidized phospholipids was estimated by the loss of EPR spectral anisotropy, and an empirical index of the related bilayer disorder was calculated from the second derivative spectra. Bilayers made of each non-oxidized phospholipid species from either ethanol solutions or REV suspensions showed the highest anisotropy, while the increasing presence of oxidized lipids in the samples resulted in progressive loss of EPR spectral anisotropy. In contrast, vesicles containing 40% of the oxidized species maintained an unaltered fluidity gradient, while REV formation was hindered by oxidized phospholipid percentages higher than 45% for ASPC and EYPC, and 35% for Mitochondrial lipids (MtL). It is concluded that the early stages of lipoperoxidation bring about disordering of the phospholipid bilayer interior rather than fluidity alterations, and that prolonged oxidation may result in loss of structural and chemical properties of the bilayer until the structure no longer holds.     

Interaction between dihydropyridines and phospholipid bilayers: a molecular dynamics simulation

Interaction of the calcium-channel antagonist dihydropyridines (DHPs), lacidipine and nifedipine, with a phospholipid bilayer was studied using 600 ps molecular dynamic simulations. We have constructed a double layer membrane model composed of 42 dimirystoyl-phosphatidylcholine molecules. The DHP molecules locate at about 7 Å from the centre of the membrane, inducing an asymmetry in the bilayer. While lacidipine did not induce significant local perturbations as judged by the gauche-trans isomerisation rate, nifedipine significantly decreased this rate, probably by producing a local rigidity of the membrane in the vicinity of the DHP.     

Effect of cholesterol on the structure and dynamic properties of unsaturated phospholipid bilayers

Molecular dynamics (MD) computer simulations of five different hydrated unsaturated phosphatidylcholine lipid bilayers built up by 18:0/18:1(n-9)cis PC, 18:0/18:2(n-6)cis PC, 18:0/18:3(n-3)cis PC, 18:0/20:4(n-6)cis PC, and 18:0/22:6(n-3)cis PC molecules with 40 mol% cholesterol, and the same five pure phosphatidylcholine bilayers have been performed at 303 K. The simulation box of a lipid bilayer contained 96 phosphatidylcholines, 64 cholesterols, and 3840 water molecules (48 phosphatidylcholine molecules and 32 cholesterols per layer and 24 water molecules per phospholipid or cholesterol in each case). The lateral self-diffusion coefficients of the lipids in these systems and mass density profiles with respect to the bilayer normal have been analyzed. It has been found that the lateral diffusion coefficients of phosphatidylcholine molecules increase with increasing number of double bonds in one of the lipid chains, both in pure bilayers and in bilayers with cholesterol. It has been found as well that the lateral diffusion coefficient of phosphatidylcholine molecules of a lipid bilayer with 40 mol% cholesterol is smaller than that for the corresponding pure phosphatidylcholine bilayer.     

Structure and interactive properties of highly fluorinated phospholipid bilayers.

Because liposomes containing fluoroalkylated phospholipids are being developed for in vivo drug delivery, the structure and interactive properties of several fluoroalkylated glycerophosphocholines (PCs) were investigated by x-ray diffraction/osmotic stress, dipole potential, and hydrophobic ion binding measurements. The lipids included PCs with highly fluorinated tails on both alkyl chains and PCs with one hydrocarbon chain and one fluoroalkylated chain. Electron density profiles showed high electron density peaks in the center of the bilayer corresponding to the fluorine atoms. The height and width of these high density peaks varied systematically, depending on the number of fluorines and their position on the alkyl chains, and on whether the bilayer was in the gel or liquid crystalline phase. Wide-angle diffraction showed that in both gel and liquid crystalline bilayers the distance between adjacent alkyl chains was greater in fluoroalkylated PCs than in analogous hydrocarbon PCs. For interbilayer separations of less than about 8 A, pressure-distance relations for fluoroalkylated PCs were similar to those previously obtained from PC bilayers with hydrocarbon chains. However, for bilayer separations greater than 8A, the total repulsive pressure depended on whether the fluoroalkylated PC was in a gel or liquid-crystalline phase. We argue that these pressure-distance relations contain contributions from both hydration and entropic repulsive pressures. Dipole potentials ranged from -680 mV for PCs with both chains fluoroalkylated to -180 mV for PCs with one chain fluoroalkylated, compared to +415 mV for egg PC. The change in dipole potential as a function of subphase concentration of tetraphenyl-boron was much larger for egg PC than for fluorinated PC monolayers, indicating that the fluorine atoms modified the binding of this hydrophobic anion. Thus, compared to conventional liposomes, liposomes made from fluoroalkylated PCs have different binding properties, which may be relevant to their use as drug carriers.     

Structural properties of a highly polyunsaturated lipid bilayer from molecular dynamics simulations.

The structure of a fully hydrated mixed (saturated/polyunsaturated) chain lipid bilayer in the biologically relevant liquid crystalline phase has been examined by performing a molecular dynamics study. The model membrane, a 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (SDPC, 18:0/22:6 PC) lipid bilayer, was investigated at constant (room) temperature and (ambient) pressure, and the results obtained in the nanosecond time scale reproduced quite well the available experimental data. Polyunsaturated fatty acids are found in high concentrations in neuronal and retinal tissues and are essential for the development of human brain function. The docosahexaenoic fatty acid, in particular, is fundamental for the proper function of the visual receptor rhodopsin. The lipid bilayer order has been investigated through the orientational order parameters. The water-lipid interface has been explored thoroughly in terms of its dimensions and the organization of the different components. Several types of interactions occurring in the system have been analyzed, specifically, the water-hydrocarbon chain, lipid-lipid and lipid-water interactions. The distribution of dihedral angles along the chains and the molecular conformations of the polyunsaturated chain of the lipids have also been studied. Special attention has been focused on the microscopic (molecular) origin of the effects of polyunsaturations on the different physical properties of membranes.     

Glycosphingolipid fatty acid arrangement in phospholipid bilayers: cholesterol effects.

Deuterium wide line NMR spectroscopy was used to study cholesterol effects on the ceramide portions of two glycosphingolipids (GSLs) distributed as minor components in fluid membranes. The common existence of very long fatty acids on GSLs was taken into account by including one glycolipid species with fatty acid chain length matching that of the host matrix, and one longer by 6 carbons. N-stearoyl and N-lignoceroyl galactosyl ceramide with perdeuterated fatty acid (18:0[d35] GalCer and 24:0[d47] GalCer) were prepared by partial synthesis. They were dispersed in bilayer membranes having the 18-carbon-fatty-acid phospholipid, 1-stearoyl-2-oleoyl-phosphatidylcholine (SOPC), as major component. Glycolipid fatty acid chain behavior and arrangement were analyzed using order profiles derived from their 2H-NMR spectra. Cholesterol effects on order parameter profiles for 18:0[d35] GalCer, with chain length equal to that of the host matrix, followed the pattern known for acyl chains of phospholipids. The presence of sterol led to restriction of trans/gauche isomerization along the length of the chain, with the largest absolute increase in order parameters being toward the surface, but somewhat greater relative effect just below the "plateau" region. In cholesterol-containing membranes, order parameter profiles for the long chain species, 24:0[d47] GalCer, showed a characteristic secondary "plateau" associated with carbon atoms C14 to C23, a feature also present in SOPC bilayers without cholesterol and in pure hydrated 24:0[d47] GalCer. Cholesterol-induced ordering effects on the long chain glycolipid were similar to those described for the shorter chain species, but were minimal at the methyl terminus.(ABSTRACT TRUNCATED AT 250 WORDS)     

Proton nuclear magnetic resonance studies on the molecular dynamics of plasmenylcholine/cholesterol and phosphatidylcholine/cholesterol bilayers

Physiologically relevant molecular species of plasmenylcholine and phosphatidylcholine were synthesized and their molecular dynamics and interactions with cholesterol were compared by determination of salient proton spin-lattice relaxation times and apparent activation energies for 1H-NMR observable motion. The molecular dynamics of PA PhosCho (1-hexadecanoyl-2-eicosatetra-5',8',11',14'-enoyl-sn-glycero-3-pho sphocholine) in multiple regions of the bilayer. Furthermore, the fluidity gradient of PA PhosCho was larger than that of PA PlasCho as ascertained by 1H spin-lattice relaxation time measurements. Introduction of cholesterol into each bilayer resulted in disparate effects on the dynamics of each subclass including: (1) increased motional freedom in the polar head group of PA PlasCho without substantial alterations in the dynamics of the polar head group of PA PhosCho; and (2) increased immobilization of the membrane interior in PA PlasCho in comparison to PA PhosCho. Analysis of Arrhenius plots of T1 relaxation times demonstrated that the apparent activation energies for vinyl and bisallylic methylene proton NMR observable motion in PA PhosCho were greater than that in PA PlasCho. Thus, comparisons of spin-lattice relaxation times and apparent activation energies demonstrate that vesicles comprised of PA PlasCho and PA PhosCho possess differential molecular dynamics and distinct interactions with cholesterol. Collectively, these results underscore the significance of the conjoint presence of the vinyl ether linkage and arachidonic acid as an important determinant of membrane dynamics in specialized mammalian membranes.     

Interactions of liquid crystal-forming molecules with phospholipid bilayers studied by molecular dynamics simulations

Recent experiments have shown that liquid crystals can be used to image mammalian cell membranes and to amplify structural reorganization in phospholipid-laden liquid crystal-aqueous interfaces. In this work, molecular dynamics simulations were employed to explore the interactions between commonly used liquid crystal-forming molecules and phospholipid bilayers. In particular, umbrella sampling was used to obtain the potential of mean force of 4-cyano-4'-pentylbiphenyl (5CB) and 4'-(3,4-difluor-phenyl)-4-pentyl-bicylohexyl (5CF) molecules partitioning into a dipalmitoylphosphatidylcholine bilayer. In addition, results of simulations are presented for systems consisting of a fully hydrated bilayer with 5CB or 5CF molecules at the lowest (4.5 mol %) and highest (20 mol %) concentrations used in recent laboratory experiments. It is found that mesogens preferentially partition from the aqueous phase into the membrane; the potential of mean force exhibits highly favorable free energy differences for partitioning (-18 k(B)T for 5CB and -26 k(B)T for 5CF). The location and orientation of mesogens associated with the most stable free energies in umbrella sampling simulations of dilute systems were found to be consistent with those observed in liquid-crystal-rich bilayers. It is found that the presence of mesogens in the bilayer enhances the order of lipid acyl tails, and changes the spatial and orientational arrangement of lipid headgroup atoms. These effects are more pronounced at higher liquid-crystal concentrations. In comparing the behavior of 5CB and 5CF, a stronger spatial correlation (i.e., possibly leading to aggregation) is observed between 5CB molecules within a bilayer than between 5CF molecules. Also, the range of molecular orientations and positions along the bilayer normal is larger for 5CB molecules. At the same time, 5CF molecules were found to bind more strongly to lipid headgroups, thereby slowing the lateral motion of lipid molecules.     

Hydration of phospholipid bilayers in the presence and absence of cholesterol

The number of water molecules bound (unfreezable) by a molecule of dipalmitoyl phosphatidylserine (DPPS) or by a molecule of dipalmitoyl phosphatidylcholine (DPPC) alone or in mixtures with cholesterol was determined by differential scanning calorimetry (DSC). When the phospholipids are in the gel state and in the absence of cholesterol, molecule of DPPS binds about 3.5 molecules of water and molecule of DPPC binds about 6 molecules of water. Number of water molecules bound increases when cholesterol crystallites are formed in the bilayer. For DPPS-cholesterol mixture at X(chol) -0.5, as well as for DPPC-cholesterol mixture at X(chol) -0.5 about 7 water molecules are bound.     

Enzymatic reduction of phospholipid and cholesterol hydroperoxides in artificial bilayers and lipoproteins

Lipid hydroperoxides (LOOHs) in various lipid assemblies are shown to be efficiently reduced and deactivated by phospholipid hydroperoxide glutathione peroxidase (PHGPX), the second selenoperoxidase to be identified and characterized. Coupled spectrophotometric analyses in the presence of NADPH, glutathione (GSH), glutathione reductase and Triton X-100 indicated that photochemically generated LOOHs in small unilamellar liposomes are substrates for PHGPX, but not for the classical glutathione peroxidase (GPX). PHGPX was found to be reactive with cholesterol hydroperoxides as well as phospholipid hydroperoxides. Kinetic iodometric analyses during GSH/PHGPX treatment of photoperoxidized liposomes indicated a rapid decay of total LOOH to a residual level of 35-40%; addition of Triton X-100 allowed the reaction to go to completion. The non-reactive LOOHs in intact liposomes were shown to be inaccessible groups on the inner membrane face. In the presence of iron and ascorbate, photoperoxidized liposomes underwent a burst of thiobarbituric acid-detectable lipid peroxidation which could be inhibited by prior GSH/PHGPX treatment, but not by GSH/GPX treatment. Additional experiments indicated that hydroperoxides of phosphatidylcholine, cholesterol and cholesteryl esters in low-density lipoprotein are also good substrates for PHGPX. An important role of PHGPX in cellular detoxification of a wide variety of LOOHs in membranes and internalized lipoproteins is suggested from these findings.     

Influence of cholesterol on phospholipid bilayers phase domains as detected by Laurdan fluorescence.

Coexisting gel and liquid-crystalline phospholipid phase domains can be observed in synthetic phospholipid vesicles during the transition from one phase to the other and, in vesicles of mixed phospholipids, at intermediate temperatures between the transitions of the different phospholipids. The presence of cholesterol perturbs the dynamic properties of both phases to such an extent as to prevent the detection of coexisting phases. 6-Lauroyl-2-dimethylaminopahthalene (Laurdan) fluorescence offers the unique advantage of well resolvable spectral parameters in the two phospholipid phases that can be used for the detection and quantitation of coexisting gel and liquid-crystalline domains. From Laurdan fluorescence excitation and emission spectra, the generalized polarization spectra and values were calculated. By the generalized polarization phospholipid phase domain coexistence can be detected, and each phase can be quantitated. In the same phospholipid vesicles where without cholesterol domain coexistence can be detected, above 15 mol% and, remarkably, at physiological cholesterol concentrations, > or = 30 mol%, no separate Laurdan fluorescence signals characteristic of distinct domains can be observed. Consequences of our results on the possible size and dynamics of phospholipid phase domains and their biological relevance are discussed.     

Short range order of hydrocarbon chains in fluid phospholipid bilayers studied by x-ray diffraction from highly oriented membranes

We present a study of the short range ordering of hydrocarbon chains in phospholipid bilayers. The x-ray peak associated with the hydrocarbon chains has been probed by means of reciprocal space mappings. Using 20 keV undulator radiation and samples of negligible mosaicity (orientational disorder), the intensity distribution is probed as a function of two coordinates, the momentum transfer parallel and perpendicular to the bilayer, over a wide range and at high resolution. Structural results are obtained concerning the distribution of tilted segments, the correlation length and the radial distribution function of the quasi two-dimensional liquid structure. A comparison is made with published molecular dynamics data (H. Heller, M. Schaefer, and K. Schulten. 1993. J. Phys. Chem. 97:8343-8360) by direct Fourier transformation of the atomic coordinates. The exact prefactor in the relationship between interchain distance and peak position is derived.     

Behavior of cholesterol and its effect on head group and chain conformations in lipid bilayers: a molecular dynamics study.

Cholesterol molecules were put into a computer-modeled hydrated bilayer of dimyristoyl phosphatidyl choline molecules, and molecular dynamics simulations were run to characterize the effect of this important molecule on membrane structure and dynamics. The effect was judged by observing differences in order parameters, tilt angles, and the fraction of gauche bonds along the hydrocarbon chains between lipids adjacent to cholesterol molecules and comparing them with those further away. It was observed that cholesterol causes an increase in the fraction of trans dihedrals and motional ordering of chains close to the rigid steroid ring system with a decrease in the kink population. The hydrogen-bonding interactions between cholesterol and lipid molecules were determined from radial distribution calculations and showed the cholesterol hydroxyl groups either solvated by water, or forming hydrogen bond contacts with the oxygens of lipid carbonyl and phosphate groups. The dynamics and conformation of the cholesterol molecules were investigated and it was seen that they had a smaller tilt with respect to the bilayer normal than the lipid chains and furthermore that the hydrocarbon tail of the cholesterol was conformationally flexible.