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.     

Effects of epicholesterol on the phosphatidylcholine bilayer: a molecular simulation study

Epicholesterol (Echol) is an epimeric form of cholesterol (Chol). A molecular dynamics simulation of the fully hydrated dimyristoylphosphatidylcholine-Echol (DMPC-Echol) bilayer membrane containing approximately 22 mol % of Echol was carried out for 5 ns. A 3-ns trajectory generated between 2 and 5 ns of molecular dynamics simulation was used for analyses to determine the effects of Echol on the membrane properties. As reference systems, pure DMPC and mixed DMPC-Chol bilayers were used. The study shows that Echol, like Chol, changes the organization of the bilayer/water interface and increases membrane order and condensation, but to a lesser degree. Effects of both sterols are based on the same atomic level mechanisms; their different strength arises from different vertical localizations of Echol and Chol hydroxyl groups in the membrane/water interface.     

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.     

Cholesterol effects on a mixed-chain phosphatidylcholine bilayer: a molecular dynamics simulation study

A molecular dynamics simulation of a mono-cis-unsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer containing approximately 22 mol% of cholesterol (POPC-Chol) was carried out for 15 ns. An 8-ns trajectory was analysed to determine the effects of Chol on the membrane properties and compare it with that on the fully saturated 1,2-dimyristoyl-phosphatidylcholine bilayer containing approximately 22 mol% of Chol (DMPC-Chol). The study suggests that the experimentally observed weaker effect of Chol on the POPC than DMPC bilayer might result from a different vertical localisation of the Chol hydroxyl group (OH-Chol) in both bilayers: in the POPC-Chol bilayer, OH-Chol is placed approximately 3 A higher in the bilayer interface than in the DMPC-Chol bilayer. Because of the rigid cis double bond in the beta-chain of POPC, Chol fits worse to the POPC-Chol membrane environment and is pushed up, in effect all Chol ring atoms are, on average, located above the double bond. Both in mono-cis-unsaturated and fully saturated PC bilayers, Chol induces stronger van der Waals interactions among the chains, whereas its interactions with the chains are weak. In contrast to DMPC, the smooth alpha-face of the Chol ring lowers the order of POPC chains, whereas the rough beta-face increases the order.     

Effects of phospholipid unsaturation on the membrane/water interface: a molecular simulation study.

Molecular dynamics (MD) simulations of fully hydrated bilayers in the liquid-crystalline state made of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) or 1-palmitoyl-2-elaidoyl-phosphatidylcholine (PEPC) were carried out to investigate the effect of the incorporation of a double bond in the phosphatidylcholine (PC) beta-chain (cis or trans) on the membrane/water interface. The bilayers reached thermal equilibrium after 3 and 1 ns of MD simulations, respectively, and productive runs were carried out for 3 ns for each bilayer. As reference systems, the 1,2-dimyristoyl-phosphatidylcholine (DMPC) bilayer (M. Pasenkiewicz-Gierula, Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi, 1999, Biophys. J. 76:1228-1240) and DMPC-cholesterol (Chol) bilayer containing 22 mol % Chol (M. Pasenkiewicz-Gierula, T. Róg, K. Kitamura, A. and Kusumi, 2000, Biophys. J. 78:1376-1389) were used. The study shows that at the interface of POPC, PEPC, and DMPC-Chol bilayers, average numbers of PC-water and PC-PC interactions are similar and, respectively, greater and smaller than in the DMPC bilayer. The average area/PC in mono-unsaturated bilayers is approximately 4 A(2) larger than in the DMPC bilayer; nevertheless, a strong correlation was found between a single molecular area (SMA) of a PC and the number of interactions this PC makes; i.e., PCs (either saturated or unsaturated) with the same SMA form similar numbers of intermolecular links. The numbers and corresponding SMAs are distributed about averages pertinent to each bilayer. No significant difference between cis and trans bonds was found.     

Nonpolar interactions between trans-membrane helical EGF peptide and phosphatidylcholines, sphingomyelins and cholesterol. Molecular dynamics simulation studies.

A molecular dynamics simulation study of four lipid bilayers with inserted trans-membrane helical fragment of epithelial growth factor (EGF) receptor (EGF peptide) was performed. The lipid bilayers differ in their lipid composition and consist of (i) unsaturated phosphatidylcholine (palmitoyloleoylphosphatidylcholine, POPC), (ii) POPC and 20 mol% of cholesterol (Chol), (iii) sphingomyelin (SM) and 20 mol% of Chol, and (iv) SM and 50 mol% of Chol. Only 1 out of 26 residues in the EGF-peptide sequence is polar (Thr). The hydrophobic thickness of each bilayer is different but shorter than the length of the peptide and so, due to hydrophobic mismatch, the inserted peptide is tilted in each bilayer. Additionally, in the POPC bilayer, which is the thinnest, the peptide loses its helical structure in a short three-amino acid fragment. This facilitates bending of the peptide and burying all hydrophobic amino acids inside the membrane core (Figure 1(b)). Bilayer lipid composition affects interactions between the peptide and lipids in the membrane core. Chol increases packing of atoms relative to the peptide side chains, and thus increases van der Waals interactions. On average, the packing around the peptide is higher in SM-based bilayers than POPC-based bilayers but for certain amino acids, packing depends on their position relative to the bilayer center. In the bilayer center, packing is higher in POPC-based bilayers, while in regions closer to the interface packing is higher in SM-based bilayers. In general, amino acids with larger side chains interact strongly with lipids, and thus the peptide sequence is important for the pattern of interactions at different membrane depths. This pattern closely resembles the shape of recently published lateral pressure profiles [Ollila et alJ. Struct. Biol. DOI:10.1016/j.jsb.2007.01.012].     

Cholesterol's inhibition effect on entering of chlorzoxazone into phosphatidylethanolamine bilayer: Relevance to cytochrome P450 drug metabolism at endoplasmic reticulum membranes

Many drugs are metabolized by cytochrome P450 (CYP) in the endoplasmic reticulum (ER) membrane. Recent studies have shown that CYP-substrate drugs reach the CYP active site after entering the lipid hydrophobic part of the ER membrane. To clarify the role of cholesterol (Chol) in the CYP-related drug metabolic process, we investigated the lipid bilayer entry of CYP-substrate drugs using a model membrane system as follows. The model membrane system comprised palmitoyl-oleoyl-phosphatidylethanolamine (POPE) and Chol. Phosphatidylethanolamine is the second major phospholipid component of ER membranes. Chlorzoxazone (CZX) was used as the CYP-substrate drug. Calorimetric measurements showed that the addition of CZX to POPE bilayers decreased the gel–liquid crystal phase transition temperature; X-ray diffraction indicated that CZX distributes into the liquid crystal phase bilayers but not practically the gel phase POPE bilayers. In the presence of Chol, dialysis and X-ray structural analyses showed that Chol inhibited CZX entry into the bilayer with an increase in Chol concentration. The Chol concentration in the ER membrane (5–10 mol%) is much lower than that in the plasma membrane (approximately 30 mol%). This fact may allow CYP-substrate drugs to enter the hydrophobic portion of the ER membrane more easily than other organelle membranes, yielding efficient drug metabolism.     

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.     

NMR study of the interaction of retinoids with phospholipid bilayers

The interaction of three vitamin A derivatives or retinoids: all-trans-retinoic acid, 13-cis-retinoic acid and retinol with multilamellar phospholipid bilayers was studied using a combination of 2H- and 31P-NMR measurements. The following model membrane systems were used: (1) dipalmitoylphosphatidylcholine (DPPC) bilayers; (2) bilayers composed of a mixture of DPPC and bovine heart phosphatidylcholine (PC); (3) mixed PC/phosphatidylethanolamine (PE) bilayers. Only a weak interaction was observed between 13-cis-retinoic acid and DPPC membranes. Addition of all-trans-retinoic acid at a molar ratio of 1:2 to the lipid causes a small decrease (5 C degrees) in the gel to liquid crystalline phase-transition temperature of DPPC, a small increase in the order parameters of the lipid side-chains of single component bilayers and no measurable effect in the other lipid systems studied. Considerably larger perturbation in the lipid bilayer structure is introduced by addition of retinol which, at a molar ratio of 1:2 to the lipid, lowered the gel to liquid crystalline phase-transition temperature of DPPC by 21 C degrees and caused a decrease of order parameters of the lipid side-chains in all three lipid bilayer systems. These effects are consistent with intercalation of retinol molecules into the bilayer interior. The results for the mixed PC/PE bilayers indicate that the presence of retinol caused lateral separation of PE- and retinol-enriched regions.     

Solid-state NMR study of membrane interactions of the pore-forming cytolysin, equinatoxin II

Equinatoxin II (EqtII) is a pore-forming protein from Actinia equina that lyses red blood cell and model membranes. Lysis is dependent on the presence of sphingomyelin (SM) and is greatest for vesicles composed of equimolar SM and phosphatidylcholine (PC). Since SM and cholesterol (Chol) interact strongly, forming domains or “rafts” in PC membranes, ³¹P and ²H solid-state NMR were used to investigate changes in the lipid order and bilayer morphology of multilamellar vesicles comprised of different ratios of dimyristoylphosphatidylcholine (DMPC), SM and Chol following addition of EqtII. The toxin affects the phase transition temperature of the lipid acyl chains, causes formation of small vesicle type structures with increasing temperature, and changes the T₂ relaxation time of the phospholipid headgroup, with a tendency to order the liquid disordered phases and disorder the more ordered lipid phases. The solid-state NMR results indicate that Chol stabilizes the DMPC bilayer in the presence of EqtII but leads to greater disruption when SM is in the bilayer. This supports the proposal that EqtII is more lytic when both SM and Chol are present as a consequence of the formation of domain boundaries between liquid ordered and disordered phases in lipid bilayers leading to membrane disruption.     

Molecular dynamics simulation of lipid packing and mobility in bilayer membranes

A model of lipid bilayer membrane in water has been developed. Parameters have been selected that allow molecular dynamics simulation of lipid bilayers in the all-atom approximation. The calculated indices of packing and mobility of lipid molecules for the liquid crystalline state of the bilayer agree well with the experimental data. Based on the model of the liquid crystalline state of the membrane, a system in the gel-like state has been constructed. The gel-state model reproduces well the packing of lipids in real bilayers, whereas the mobility of molecules proves to be overestimated.     

Coupling of cholesterol-rich lipid phases in asymmetric bilayers

We showed previously that cholesterol-rich liquid-ordered domains with lipid compositions typically found in the outer leaflet of plasma membranes induce liquid-ordered domains in adjacent regions of asymmetric lipid bilayers with apposed leaflets composed of typical inner leaflet lipid mixtures [Kiessling, V., Crane, J. M., and Tamm, L. K. (2006) Biophys. J. 91, 3313-26]. To further examine the nature of transbilayer couplings in asymmetric cholesterol-rich lipid bilayers, the effects on the lipid phase behavior in asymmetric bilayers of different lipid compositions were investigated. We established systems containing several combinations of natural extracted and synthetic lipids that exhibited coexisting liquid-ordered (lo) and liquid-disordered (ld) domains in a supported bilayer format. We find that lo phase domains are induced in all quaternary inner leaflet combinations composed of PCs, PEs, PSs, and cholesterol. Ternary mixtures of PCs/PEs/Chol, PCs/PSs/Chol also exhibit lo phases adjacent to outer leaflet lo phases. However, with the exception of brain PC extracts, binary PC/Chol mixtures are not induced to form lo phases by adjacent outer leaflet lo phases. Higher melting lipid ad-mixtures of PEs and PSs are needed for lo phase induction in the inner leaflet. It appears that the phase behavior of the inner leaflet mixtures is dominated by the intrinsic chain melting temperatures of the lipid components, rather than by their specific headgroup classes. In addition, similar studies with synthetic, completely saturated lipids and cholesterol show that lipid oxidation is not a factor in the observed phase behavior.     

Molecular dynamics simulations and 2H NMR study of the GalCer/DPPG lipid bilayer

Molecular dynamics simulations were performed on a two-component lipid bilayer system in the liquid crystalline phase at constant pressure and constant temperature. The lipid bilayers were composed of a mixture of neutral galactosylceramide (GalCer) and charged dipalmitoylphosphatidylglycerol (DPPG) lipid molecules. Two lipid bilayer systems were prepared with GalCer:DPPG ratio 9:1 (10%-DPPG system) and 3:1 (25%-DPPG system). The 10%-DPPG system represents a collapsed state lipid bilayer, with a narrow water space between the bilayers, and the 25%-DPPG system represents an expanded state with a fluid space of approximately 10 nm. The number of lipid molecules used in each simulation was 1024, and the length of the production run simulation was 10 ns. The simulations were validated by comparing the results with experimental data for several important aspects of the bilayer structure and dynamics. Deuterium order parameters obtained from (2)H NMR experiments for DPPG chains are in a very good agreement with those obtained from molecular dynamics simulations. The surface area per GalCer lipid molecule was estimated to be 0.608 +/- 0.011 nm(2). From the simulated electron density profiles, the bilayer thickness defined as the distance between the phosphorus peaks across the bilayer was calculated to be 4.21 nm. Both simulation systems revealed a tendency for cooperative bilayer undulations, as expected in the liquid crystalline phase. The interaction of water with the GalCer and DPPG oxygen atoms results in a strong water ordering in a spherical hydration shell and the formation of hydrogen bonds (H-bonds). Each GalCer lipid molecule makes 8.6 +/- 0.1 H-bonds with the surrounding water, whereas each DPPG lipid molecule makes 8.3 +/- 0.1 H-bonds. The number of water molecules per GalCer or DPPG in the hydration shell was estimated to be 10-11 from an analysis of the radial distribution functions. The formation of the intermolecular hydrogen bonds was observed between hydroxyl groups from the opposing GalCer sugar headgroups, giving an energy of adhesion in the range between -1.0 and -3.4 erg/cm(2). We suggest that this value is the contribution of the hydrogen-bond component to the net adhesion energy between GalCer bilayers in the liquid crystalline phase.     

Saturation with cholesterol increases vertical order and smoothes the surface of the phosphatidylcholine bilayer: a molecular simulation study

Molecular dynamics (MD) simulations of a mono-cis-unsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer and a POPC bilayer containing 50mol% cholesterol (POPC-Chol50) were carried out for 200ns to compare the spatial organizations of the pure POPC bilayer and the POPC bilayer saturated with Chol. The results presented here indicate that saturation with Chol significantly narrows the distribution of vertical positions of the center-of-mass of POPC molecules and POPC atoms in the bilayer. In the POPC-Chol50 bilayer, the same moieties of the lipid molecules are better aligned at a given bilayer depth, forming the following clearly separated membrane regions: the polar headgroup, the rigid core consisting of steroid rings and upper fragments of the acyl chains, and the fluid hydrocarbon core consisting of Chol chains and the lower fragments of POPC chains. The membrane surface of the POPC-Chol50 bilayer is smooth. The results have biological significance because the POPC-Chol50 bilayer models the bulk phospholipid portion of the fiber-cell membrane in the eye lens. It is hypothesized that in the eye lens cholesterol-induced smoothing of the membrane surface decreases light-scattering and helps to maintain lens transparency.     

Study of curcumin behavior in two different lipid bilayer models of liposomal curcumin using molecular dynamics simulation

Liposomal formulation of curcumin is an important therapeutic agent for the treatment of various cancers. Despite extensive studies on the biological effects of this formulation in cancer treatment, much remains unknown about curcumin–liposome interactions. Understanding how different lipid bilayers respond to curcumin molecule may help us to design more effective liposomal curcumin. Here, we used molecular dynamics simulation method to investigate the behavior of curcumin in two lipid bilayers commonly used in preparation of liposomal curcumin, namely dipalmitoylphosphatidylcholine (DPPC) and dimyristoylphosphatidylglycerol (DMPG). First, the free energy barriers for translocation of one curcumin molecule from water to the lipid bilayer were determined by using the potential of mean force (PMF). The computed free energy profile exhibits a global minimum at the solvent–headgroup interface (LH region) for both lipid membranes. We also evaluated the free energy difference between the equilibrium position of curcumin in the lipid bilayer and bulk water as the excess chemical potential. Our results show that curcumin has the higher affinity in DMPG compared to DPPC lipid bilayer (?8.39 vs. ?1.69 k_BT) and this is related to more hydrogen bond possibility for curcumin in DMPG lipid membrane. Next, using an unconstrained molecular dynamic simulation with curcumin initially positioned at the center of lipid bilayer, we studied various properties of each lipid bilayer system in the presence of curcumin molecule that was in full agreement with PMF and experimental data. The results of these simulation studies suggest that membrane composition could have a large effect on interaction of curcumin–lipid bilayer.     

Molecular dynamics simulations of charged and neutral lipid bilayers: treatment of electrostatic interactions

Molecular dynamics (MD) simulations complement experimental methods in studies of the structure and dynamics of lipid bilayers. The choice of algorithms employed in this computational method represents a trade-off between the accuracy and real calculation time. The largest portion of the simulation time is devoted to calculation of long-range electrostatic interactions. To speed-up evaluation of these interactions, various approximations have been used. The most common ones are the truncation of long-range interactions with the use of cut-offs, and the particle-mesh Ewald (PME) method. In this study, several multi-nanosecond cut-off and PME simulations were performed to establish the influence of the simulation protocol on the bilayer properties. Two bilayers were used. One consisted of neutral phosphatidylcholine molecules. The other was a mixed lipid bilayer consisting of neutral phosphatidylethanolamine and negatively charged phosphatidylglycerol molecules. The study shows that the cut-off simulation of a bilayer containing charge molecules generates artefacts; in particular the mobility and order of the charged molecules are vastly different from those determined experimentally. In the PME simulation, the bilayer properties are in general agreement with experimental data. The cut-off simulation of bilayers containing only uncharged molecules does not generate artefacts, nevertheless, the PME simulation gives generally better agreement with experimental data.     

Lipid bilayer deformation and the free energy of interaction of a Kv channel gating-modifier toxin

A number of membrane proteins act via binding at the water/lipid bilayer interface. An important example of such proteins is provided by the gating-modifier toxins that act on voltage-gated potassium (Kv) channels. They are thought to partition to the headgroup region of lipid bilayers, and so provide a good system for probing the nature of interactions of a protein with the water/bilayer interface. We used coarse-grained molecular dynamics simulations to compute the one-dimensional potential of mean force (i.e., free energy) profile that governs the interaction between a Kv channel gating-modifier toxin (VSTx1) and model phospholipid bilayers. The reaction coordinate sampled corresponds to the position of the toxin along the bilayer normal. The course-grained representation of the protein and lipids enabled us to explore extended time periods, revealing aspects of toxin/bilayer dynamics and energetics that would be difficult to observe on the timescales currently afforded by atomistic molecular dynamics simulations. In particular, we show for this model system that the bilayer deforms as it interacts with the toxin, and that such deformations perturb the free energy profile. Bilayer deformation therefore adds an additional layer of complexity to be addressed in investigations of membrane/protein systems. In particular, one should allow for local deformations that may arise due to the spatial array of charged and hydrophobic elements of an interfacially located membrane protein.     

Effects of a carane derivative local anesthetic on a phospholipid bilayer studied by molecular dynamics simulation

Molecular dynamics (MD) simulations of two hydrated palmitoyloleoylphosphatidylcholine (POPC) bilayers each containing eight carane derivative (KP-23) local anesthetic (LA) molecules in neutral (POPC-LA) or protonated (POPC-LAH) forms were carried out to investigate the effect of KP-23 and its protonation on the bilayer. 3-ns trajectories were used for analyses. A pure POPC bilayer was employed as a reference system. In both POPC-LA and POPC-LAH systems a few KP-23 molecules intercalated into the bilayer and moved near the bilayer/water interface. They were located on the hydrophobic core side of the interface in the POPC-LA bilayer, but on the water phase side in the POPC-LAH bilayer. The order of the POPC chains was higher in the POPC-LA bilayer than in the pure POPC bilayer and was lower in the POPC-LAH bilayer. Interactions between polar groups of KP-23 and POPC or water were responsible for a lower hydration of POPC headgroups in POPC bilayers containing KP-23 than in the pure POPC bilayer. KP-23 molecules were found to form aggregates both in POPC-LA and POPC-LAH bilayers. Due to higher amphiphilicity of LAH, the LAH aggregate was more micelle-like and larger than the LA one. The results demonstrate the rapid timescales of the initial processes that take place at and near the bilayer interface as well as details of the atomic level interactions between local anesthetic and the lipid matrix of a cell membrane.     

Cholesterol-sphingomyelin interactions: a molecular dynamics simulation study

Stearoylsphingomyelin (SSM) bilayers containing 0, 22, and 50 mol % cholesterol (Chol) and a pentadecanoyl-stearoylphosphatidylcholine (15SPC) bilayer containing 22 mol % Chol were molecular dynamics simulated at two temperatures (37 degrees C and 60 degrees C). 15SPC is the best PC equivalent of SSM. The Chol effect on the SSM bilayer differs significantly from that on the 15SPC bilayer. At the same temperature and Chol content, H-bonding of Chol with SSM is more extensive than with 15SPC. SSM-Chol H-bonding anchors the OH group of Chol in the lower regions of the SSM-Chol bilayer interface. Such a location strengthens the influence of Chol on the SSM chains. In effect, the phase of the SSM-Chol bilayer containing 22 mol % Chol at 37 degrees C is shifted from the gel to the liquid-ordered phase, and the bilayer displays similar properties below and above the main phase-transition temperature for a pure SSM bilayer of approximately 45 degrees C. In contrast, due to a higher location, Chol is not able to change the phase of the 15SPC-Chol bilayer, which at 37 degrees C remains in the gel phase. Chol affects both the core and interface of the SSM bilayer. With increasing Chol content, the order of SSM chains and hydration of SSM headgroups increase, whereas polar interactions between lipids decrease.     

Membrane hydration and structure on a subnanometer scale as seen by high resolution solid state nuclear magnetic resonance: POPC and POPC/C12EO4 model membranes.

The position on a subnanometer scale and the dynamics of structurally important water in model membranes was determined using a combination of proton magic-angle spinning NMR (MAS) with two-dimensional NOESY NMR techniques. Here, we report studies on phosphocholine lipid bilayers that were then modified by the addition of a nonionic surfactant that is shown to dehydrate the lipid. These studies are supplemented by 13C magic-angle spinning NMR investigations to get information on the dynamics of segmental motions of the membrane molecules. It can be shown that the hydrophilic chain of the surfactant is positioned at least partially within the hydrophobic core of the lipid bilayer. With the above NMR approach, we are able to establish molecular contacts between water and the lipid headgroup as well as with certain groups of the hydrocarbon chains and the glycerol backbone. This is possible because high resolution proton and 13C-NMR spectra of multilamellar bilayer membranes are obtained using MAS. A phase-sensitive NOESY must also be applied to distinguish positive and negative cross-peaks in the two-dimensional plot. These studies have high potential to investigate membrane proteins hydration and structural organization in a natural lipid bilayer surrounding.