Influence of cholesterol on dynamics of dimyristoylphosphatidylcholine bilayers as studied by deuterium NMR relaxation

Polarizable atomic multipole-based force field for cholesterol

Cholesterol is one of the essential component of lipid in membrane. We present a polarizable atomic multipole force field (FF) for the molecular dynamic simulation of cholesterol. The FF building process follows the computational framework as the atomic multipole optimized energetics for biomolecular applications (AMOEBA) model. In this framework, the electronics parameters, including atomic monopole moments, dipole moments, and quadrupole moments calculated from ab initio calculations in the gas phase, are applied to represent the charge distribution. Furthermore, the many-body polarization is modeled by following the same pattern of distributed atomic polarizabilities. Then, the bilayers composed of two typical phospholipid molecules, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), in a range of different cholesterol concentrations are built and implemented by molecular dynamics (MD) simulations based on the proposed polarizable FF. The simulation results are statistically analyzed to validate the feasibility of the proposed FF. The structural properties of the bilayers are calculated to compare with the related experimental values. The MD values show the same trend of experimental values changes.Communicated by Ramaswamy H. Sarma.

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

Biophysics of Membrane Stiffening by Cholesterol and Phosphatidylinositol 4,5-bisphosphate (PIP2)

Cell membranes regulate a wide range of phenomena that are implicated in key cellular functions. Cholesterol, a critical component of eukaryotic cell membranes, is responsible for cellular organization, membrane elasticity, and other critical physicochemical parameters. Besides cholesterol, other lipid components such as phosphatidylinositol 4,5-bisphosphate (PIP2) are found in minor concentrations in cell membranes yet can also play a major regulatory role in various cell functions. In this chapter, we describe how solid-state deuterium nuclear magnetic resonance (²H NMR) spectroscopy together with neutron spin-echo (NSE) spectroscopy can inform synergetic changes to lipid molecular packing due to cholesterol and PIP2 that modulate the bending rigidity of lipid membranes. Fundamental structure-property relations of molecular self-assembly are illuminated and point toward a length and time-scale dependence of cell membrane mechanics, with significant implications for biological activity and membrane lipid-protein interactions.

Interdomain exchange and the flip-flop of cholesterol in ternary component lipid membranes and their effects on heterogeneous cholesterol diffusion

Cell membranes are heterogeneous with a variety of lipids, cholesterol, and proteins and are composed of domains of different compositions. Such heterogeneous environments make the transport of cholesterol complicated: cholesterol not only diffuses within a particular domain but also travels between domains. Cholesterol also flip-flops between upper and lower leaflets such that cholesterol may reside both within leaflets and in the central region between two leaflets. How the presence of multiple domains and the interdomain exchange of cholesterol would affect the cholesterol transport, however, remains elusive. In this study, therefore, we perform molecular dynamics simulations up to 100μs for ternary component lipid membranes, which consist of saturated lipids (dipalmitoylphosphatidylcholine, DPPC), unsaturated lipids (dilinoleylphosphatidylcholine, DIPC), and cholesterol. The ternary component membranes in our simulations form two domains readily: DPPC and DIPC domains. We find that the diffusion of cholesterol molecules is much more heterogeneous and non-Gaussian than expected for binary component lipid membranes of lipids and cholesterol. The non-Gaussian parameter of the cholesterol molecules is about four times larger in the ternary component lipid membranes than in the binary component lipid membranes. Such non-Gaussian and heterogeneous transport of cholesterol arises from the interplay among the interdomain kinetics, the different diffusivity of cholesterol in different domains, and the flip-flop of cholesterol. This suggests that in cell membranes that consist of various domains and proteins, the cholesterol transport can be very heterogeneous. We also find that the mechanism of the interdomain exchange differs for different domains: cholesterol tends to exit the DIPC domain along the central region of the membrane for the DIPC-to-DPPC transition, while the cholesterol is likely to exit the DPPC domain within the membrane leaflet for the DPPC-to-DIPC transition. Also, the interdomain exchange kinetics of cholesterol for the DPPC-to-DIPC transition is up to 7.9 times slower than the DIPC-to-DPPC transition.

Binding and dynamics of melatonin at the interface of phosphatidylcholine-cholesterol membranes

The characterization of interactions between melatonin, one main ingredient of medicines regulating sleeping rhythms, and basic components of cellular plasma membranes (phospholipids, cholesterol, metal ions and water) is very important to elucidate the main mechanisms for the introduction of melatonin into cells and also to identify its local structure and microscopic dynamics. Molecular dynamics simulations of melatonin inside mixtures of dimyristoylphosphatidylcholine and cholesterol in NaCl solution at physiological concentration have been performed at 303.15 K to systematically explore melatonin-cholesterol, melatonin-lipid and melatonin-water interactions. Properties such as the area per lipid and thickness of the membrane as well as selected radial distribution functions, binding free energies, angular distributions, atomic spectral densities and translational diffusion of melatonin are reported. The presence of cholesterol significantly affects the behavior of melatonin, which is mainly buried into the interfaces of membranes. Introducing cholesterol into the system helps melatonin change from folded to extended configurations more easily. Our results suggest that there exists a competition between the binding of melatonin to phospholipids and to cholesterol by means of hydrogen-bonds. Spectral densities of melatonin reported in this work, in overall good agreement with experimental data, revealed the participation of each atom of melatonin to its complete spectrum. Melatonin self-diffusion coefficients are of the order of 10-7 cm2/s and they significantly increase when cholesterol is addeed to the membrane.

The Phase Behavior and Organization of Sphingomyelin/Cholesterol Membranes: A Deuterium NMR Study

We have studied the dependence of the phase and domain characteristics of sphingomyelin (SM)/cholesterol model membranes on sterol content and temperature using deuterium nuclear magnetic resonance. NMR spectra of N-palmitoyl(D31)-D-erythro-sphingosylphosphorylcholine (PSM-d31) were taken for temperatures from 25 to 70°C and cholesterol concentrations of 0-40%. Analogous experiments were performed using 1-palmitoyl,2-palmitoyl(D31)-sn-glycero-3-phosphocholine (DPPC-d31)/cholesterol membranes to carefully compare the data obtained using palmitoyl chains that have similar "kinked" conformations. The constructed phase diagrams exhibit both solid-ordered (so) + liquid-ordered (lo) and liquid-disordered (ld) + lo phase-coexistence regions with a clear three-phase line. Macroscopic (micron-sized) coexistence of ld and lo phases was not observed; instead, line-broadening in the ld+lo region was characterized by intermediate exchange of lipids between the two types of domains. The length scales associated with the domains were estimated to be 75-150 nm for PSM-d31/cholesterol and DPPC-d31/cholesterol model membranes.

Flexible lipid nanomaterials studied by NMR spectroscopy

Advances in solid-state nuclear magnetic resonance spectroscopy inform the emergence of material properties from atomistic-level interactions in membrane lipid nanostructures.

Cholesterol Effects on the Physical Properties of Lipid Membranes Viewed by Solid-state NMR Spectroscopy

In this chapter, we review the physical properties of lipid/cholesterol mixtures involving studies of model membranes using solid-state NMR spectroscopy. The approach allows one to quantify the average membrane structure, fluctuations, and elastic deformation upon cholesterol interaction. Emphasis is placed on understanding the membrane structural deformation and emergent fluctuations at an atomistic level. Lineshape measurements using solid-state NMR spectroscopy give equilibrium structural properties, while relaxation time measurements study the molecular dynamics over a wide timescale range. The equilibrium properties of glycerophospholipids, sphingolipids, and their binary and tertiary mixtures with cholesterol are accessible. Nonideal mixing of cholesterol with other lipids explains the occurrence of liquid-ordered domains. The entropic loss upon addition of cholesterol to sphingolipids is less than for glycerophospholipids, and may drive formation of lipid rafts. The functional dependence of ²H NMR spin-lattice relaxation (R _1Z) rates on segmental order parameters (S _CD) for lipid membranes is indicative of emergent viscoelastic properties. Addition of cholesterol shows stiffening of the bilayer relative to the pure lipids and this effect is diminished for lanosterol. Opposite influences of cholesterol and detergents on collective dynamics and elasticity at an atomistic scale can potentially affect lipid raft formation in cellular membranes.

Effects of cholesterol on the binding of the precursor neurotransmitter tryptophan to zwitterionic membranes

The characterization of the microscopical forces between the essential α-amino-acid tryptophan, precursor of the neurotransmitter serotonin and of the hormone melatonin, and the basic components of cell membranes and their environments (phospholipids, cholesterol, ionic species, and water) is of central importance to elucidate their local structure and dynamics as well as the mechanisms responsible for the access of tryptophan to the interior of the cell. We have performed nanosecond molecular dynamics simulations of tryptophan embedded in model zwitterionic bilayer membranes made by di-palmitoyl-phosphatidyl-choline and cholesterol inside aqueous sodium-chloride solution in order to systematically examine tryptophan-lipid, tryptophan-cholesterol, and tryptophan-water interactions under liquid-crystalline phase conditions. Microscopic properties such as the area per lipid, lipid thickness, radial distribution functions, hydrogen-bonding lengths, atomic spectral densities, and self-diffusion coefficients have been evaluated. Our results show that the presence of tryptophan significantly affects the structure and dynamics of the membrane. Tryptophan spends long periods of time at the water-membrane interface, and it plays a central role by bridging a few lipids and cholesterol chains by means of hydrogen-bonds. The computed spectral densities, in excellent agreement with experimental infrared and Raman data, revealed the participation of each atomic site of tryptophan to the complete spectrum of the molecule. Tryptophan self-diffusion coefficients have been found to be in between 10^-7 and 10^-6 cm²/s and strongly depending of the concentration of cholesterol in the system.

Facilitated and Non-Gaussian Diffusion of Cholesterol in Liquid Ordered Phase Bilayers Depends on the Flip-Flop and Spatial Arrangement of Cholesterol

The diffusion of cholesterol in biological membranes is critical to cellular processes such as the formation of cholesterol-enriched domains. The cholesterol diffusion may be complicated especially when cholesterol flip-flops and/or stays at the membrane center. Understanding the diffusion mechanism of cholesterol at a molecular level should be, therefore, a topic of interest. We perform molecular dynamics simulations up to 100 μs for lipid bilayers with various concentrations of cholesterol. We find that cholesterol diffusion in the liquid ordered phase depends on whether it is within leaflets or at the bilayer center, is non-Gaussian for several microseconds, and is enhanced significantly compared to that of lipids. Cholesterol at the bilayer center diffuses fast, while cholesterol in the hydrocarbon region with upright orientation diffuses relatively slowly. Such position-dependent dynamics of cholesterol leads to facilitated and non-Gaussian diffusion.

Concepts and Methods of Solid-State NMR Spectroscopy Applied to Biomembranes

Concepts of solid-state NMR spectroscopy and applications to fluid membranes are reviewed in this paper. Membrane lipids with ²H-labeled acyl chains or polar head groups are studied using ²H NMR to yield knowledge of their atomistic structures in relation to equilibrium properties. This review demonstrates the principles and applications of solid-state NMR by unifying dipolar and quadrupolar interactions and highlights the unique features offered by solid-state ²H NMR with experimental illustrations. For randomly oriented multilamellar lipids or aligned membranes, solid-state ²H NMR enables direct measurement of residual quadrupolar couplings (RQCs) due to individual C-²H-labeled segments. The distribution of RQC values gives nearly complete profiles of the segmental order parameters S_CD⁽ⁱ⁾ as a function of acyl segment position (i). Alternatively, one can measure residual dipolar couplings (RDCs) for natural abundance lipid samples to obtain segmental S_CH order parameters. A theoretical mean-torque model provides acyl-packing profiles representing the cumulative chain extension along the normal to the aqueous interface. Equilibrium structural properties of fluid bilayers and various thermodynamic quantities can then be calculated, which describe the interactions with cholesterol, detergents, peptides, and integral membrane proteins and formation of lipid rafts. One can also obtain direct information for membrane-bound peptides or proteins by measuring RDCs using magic-angle spinning (MAS) in combination with dipolar recoupling methods. Solid-state NMR methods have been extensively applied to characterize model membranes and membrane-bound peptides and proteins, giving unique information on their conformations, orientations, and interactions in the natural liquid-crystalline state.

Probing Hydronium Ion Histidine NH Exchange Rate Constants in the M2 Channel via Indirect Observation of Dipolar-Dephased 15N Signals in Magic-Angle-Spinning NMR

Water-protein chemical exchange in membrane-bound proteins is an important parameter for understanding how proteins interact with their aqueous environment, but has been difficult to observe in membrane-bound biological systems. Here, we demonstrate the feasibility of probing specific water-protein chemical exchange in membrane-bound proteins in solid-state MAS NMR. By spin-locking the 1H magnetization along the magic angle, the 1H spin diffusion is suppressed such that a water-protein chemical exchange process can be monitored indirectly by dipolar-dephased 15N signals through polarization transfer from 1H. In the example of the Influenza A full length M2 protein, the buildup of dipolar-dephased 15N signals from the tetrad of His37 side chains have been observed as a function of spin-lock time. This confirms that hydronium ions are in exchange with protons in the His37 NH bonds at the heart of the M2 proton conduction mechanism, with an exchange rate constant of ∼1750 s-1 for pH 6.2 at -10 °C.

Cholesterol-induced suppression of membrane elastic fluctuations at the atomistic level

Applications of solid-state NMR spectroscopy for investigating the influences of lipid-cholesterol interactions on membrane fluctuations are reviewed in this paper. Emphasis is placed on understanding the energy landscapes and fluctuations at an emergent atomistic level. Solid-state (2)H NMR spectroscopy directly measures residual quadrupolar couplings (RQCs) due to individual C-(2)H labeled segments of the lipid molecules. Moreover, residual dipolar couplings (RDCs) of (13)C-(1)H bonds are obtained in separated local-field NMR spectroscopy. The distributions of RQC or RDC values give nearly complete profiles of the order parameters as a function of acyl segment position. Measured equilibrium properties of glycerophospholipids and sphingolipids including their binary and tertiary mixtures with cholesterol show unequal mixing associated with liquid-ordered domains. The entropic loss upon addition of cholesterol to sphingolipids is less than for glycerophospholipids and may drive the formation of lipid rafts. In addition relaxation time measurements enable one to study the molecular dynamics over a wide time-scale range. For (2)H NMR the experimental spin-lattice (R1Z) relaxation rates follow a theoretical square-law dependence on segmental order parameters (SCD) due to collective slow dynamics over mesoscopic length scales. The functional dependence for the liquid-crystalline lipid membranes is indicative of viscoelastic properties as they emerge from atomistic-level interactions. A striking decrease in square-law slope upon addition of cholesterol denotes stiffening relative to the pure lipid bilayers that is diminished in the case of lanosterol. Measured equilibrium properties and relaxation rates infer opposite influences of cholesterol and detergents on collective dynamics and elasticity at an atomistic scale that potentially affects lipid raft formation in cellular membranes.

Molecular, dynamic, and structural origin of inhomogeneous magnetization transfer in lipid membranes

PURPOSE

To elucidate the dynamic, structural, and molecular properties that create inhomogeneous magnetization transfer (ihMT) contrast.

METHODS

Amphiphilic lipids, lamellar phospholipids with cholesterol, and bovine spinal cord (BSC) specimens were examined along with nonlipid systems. Magnetization transfer (MT), enhanced MT (eMT, obtained with double-sided radiofrequency saturation), ihMT (MT - eMT), and dipolar relaxation, T_1D , were measured at 2.0 and 11.7 T.

RESULTS

The amplitude of ihMT ratio (ihMTR) is positively correlated with T_1D values. Both ihMTR and T_1D increase with increasing temperature in BSC white matter and in phospholipids and decrease with temperature in other lipids. Changes in ihMTR with temperature arise primarily from alterations in MT rather than eMT. Spectral width of MT, eMT, and ihMT increases with increasing carbon chain length.

CONCLUSIONS

Concerted motions of phospholipids in white matter decrease proton spin diffusion leading to increased proton T_1D times and increased ihMT amplitudes, consistent with decoupling of Zeeman and dipolar spin reservoirs. Molecular specificity and dynamic sensitivity of ihMT contrast make it a suitable candidate for probing myelin membrane disorders. Magn Reson Med 77:1318-1328, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Formation of Gel-like Nanodomains in Cholesterol-Containing Sphingomyelin or Phosphatidylcholine Binary Membrane As Examined by Fluorescence Lifetimes and (2)H NMR Spectra

In this study, we measured the time-resolved fluorescence of trans-parinaric acid (tPA), steady-state fluorescence anisotropy of diphenylhexatriene (DPH), and (2)H NMR of 10,10-d2-stearoyl lipids in stearoyl sphingomyelin with cholesterol (SSM/Chol) and l-palmitoyl-2-stearoyl-sn-glycero-3-phosphocholine with Chol (PSPC/Chol) binary membranes. The results suggest that the membrane order obtained from the fluorescence experiments shows a similar temperature dependency as those of the (2)H NMR data. More importantly, the time-resolved fluorescence data implied the presence of at least two types of domains, cholesterol-poor gel-like domains (CPGLD) and cholesterol-enriched liquid-ordered (Lo) domains. These domains appear on a nano-to-micro second time scale for both SSM-Chol and PSPC-Chol membranes. The relative size of the gel-like domain was also estimated from the temperature-dependent lifetime measurements and (2)H NMR spectral changes. The results imply that the size of the gel-like domains is very small, probably on the nanometer scale, and smaller in SSM-Chol membrane than those in PSPC-Chol bilayers, which could account for the higher thermal stability of SM-Chol membranes. The present study demonstrates that gel-like nanodomains occur in SM-Chol binary membrane even with Chol content of over 33 mol %, which has been thought to consist exclusively of Lo phase, implying that not only Lo domains but also gel-like nanodomains are important for formation of lipid-ordered phase in SM-Chol and PC-Chol membranes.

A Parameterization of Cholesterol for Mixed Lipid Bilayer Simulation within the Amber Lipid14 Force Field

The Amber Lipid14 force field is expanded to include cholesterol parameters for all-atom cholesterol and lipid bilayer molecular dynamics simulations. The General Amber and Lipid14 force fields are used as a basis for assigning atom types and basic parameters. A new RESP charge derivation for cholesterol is presented, and tail parameters are adapted from Lipid14 alkane tails. 1,2-Dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers are simulated at a range of cholesterol contents. Experimental bilayer structural properties are compared with bilayer simulations and are found to be in good agreement. With this parameterization, another component of complex membranes is available for molecular dynamics with the Amber Lipid14 force field.

Diffusion and spectroscopy of water and lipids in fully hydrated dimyristoylphosphatidylcholine bilayer membranes

Microscopic structure and dynamics of water and lipids in a fully hydrated dimyristoylphosphatidylcholine phospholipid lipid bilayer membrane in the liquid-crystalline phase have been analyzed with all-atom molecular dynamics simulations based on the recently parameterized CHARMM36 force field. The diffusive dynamics of the membrane lipids and of its hydration water, their reorientational motions as well as their corresponding spectral densities, related to the absorption of radiation, have been considered for the first time using the present force field. In addition, structural properties such as density and pressure profiles, a deuterium-order parameter, surface tension, and the extent of water penetration in the membrane have been analyzed. Molecular self-diffusion, reorientational motions, and spectral densities of atomic species reveal a variety of time scales playing a role in membrane dynamics. The mechanisms of lipid motion strongly depend on the time scale considered, from fast ballistic translation at the scale of picoseconds (effective diffusion coefficients of the order of 10(-5) cm(2)/s) to diffusive flow of a few lipids forming nanodomains at the scale of hundreds of nanoseconds (diffusion coefficients of the order of 10(-8) cm(2)/s). In the intermediate regime of sub-diffusion, collisions with nearest neighbors prevent the lipids to achieve full diffusion. Lipid reorientations along selected directions agree well with reported nuclear magnetic resonance data and indicate two different time scales, one about 1 ns and a second one in the range of 2-8 ns. We associated the two time scales of reorientational motions with angular distributions of selected vectors. Calculated spectral densities corresponding to lipid and water reveal an overall good qualitative agreement with Fourier transform infrared spectroscopy experiments. Our simulations indicate a blue-shift of the low frequency spectral bands of hydration water as a result of its interaction with lipids. We have thoroughly analyzed the physical meaning of all spectral features from lipid atomic sites and correlated them with experimental data. Our findings include a "wagging of the tails" frequency around 30 cm(-1), which essentially corresponds to motions of the tail-group along the instantaneous plane formed by the two lipid tails, i.e., in-plane oscillations are clearly of bigger importance than those along the normal-to-the plane direction.

Detailed comparison of deuterium quadrupole profiles between sphingomyelin and phosphatidylcholine bilayers

Lipid rafts are microdomains rich in sphingomyelin (SM) and cholesterol (Chol). The essential question is why natural lipid rafts prefer SM rather than saturated diacyl glycerophosphocholine, although both form ordered membranes with Chol in model systems. Hence in this study, we synthesized site-specifically deuterated 1-palmitoyl-2-stearoyl-sn-glycero-3-phosphocholines that match the acyl chain length of stearoyl-SM (SSM), and compared their deuterium quadrupole coupling profiles in detail. The results suggest a deeper distribution of Chol in the SSM membranes, a lower entropic penalty upon accommodation of Chol in SSM membranes, and a higher thermal stability of acyl-chain orders in the SSM-Chol bilayers than in the 1-palmitoyl-2-stearoyl-sn-glycero-3-phosphocholine-Chol system at various Chol concentrations. The entropy effect and thermal stability should render SM a more preferred raft constituent than saturated diacyl glycerophosphocholine. Our data also demonstrate that the selective and comprehensive deuteration strategy is indispensable for accurate comparison of order profiles.

Sterols associated with small unilamellar vesicles (SUVs): intrinsic mobility role for 1H NMR detection

Small unilamellar vesicles (SUVs) of phospholipids are often used as a membrane model system for studying the interaction of molecules. When using NMR under the standard liquid-state conditions, SUV phospholipid proton spectra can be recorded, exhibiting sharp signals. This is not only because of the fast vesicular tumbling but also because of the combination of this tumbling with the individual motion of the lipids inside the bilayer. This appears evident because addition of cholesterol is responsible of broader resonances because of the slowing down of the lipid motion. On the other hand, no (1)H signal is detected for cholesterol in the bilayer. This lack of detection of the inserted molecules explains why generally SUVs are not considered as a good model for NMR studies under the standard liquid-state conditions. Here, we use two other sterols in order to demonstrate that an increase of the molecular mobility inside the bilayer could allow the detection of their proton resonances. For desmosterol and lanosterol, which show higher mobility inside the bilayer, with increasing lateral diffusion rates, (1)H sterol signals are detected in contrast to cholesterol. For the fast diffusing lanosterol, no significant improvement in detection is observed using deuterated lipids, demonstrating that homonuclear dipolar coupling is fully averaged out. Furthermore, in the case of low mobility such as for cholesterol, the use of a fast magic angle spinning probe is shown to be efficient to recover the full proton spectrum.

The effect of cholesterol on membrane dynamics on different timescales in lipid bilayers from fast field-cycling NMR relaxometry studies of unilamellar vesicles

The general applicability of fast field-cycling nuclear magnetic resonance relaxometry in the study of dynamics in lipid bilayers is demonstrated through analysis of binary unilamellar liposomes composed of 1,2-dioleoyl-sn-glycero-3-posphocholine (DOPC) and cholesterol. We extend an evidence-based method to simulating the NMR relaxation response, previously validated for single-component membranes, to evaluate the effect of the sterol molecule on local ordering and dynamics over multiple timescales. The relaxometric results are found to be most consistent with the partitioning of the lipid molecules into affected and unaffected portions, rather than a single averaged phase. Our analysis suggests that up to 25 mol%, each cholesterol molecule orders three DOPC molecules, providing experimental backup to the findings of many molecular dynamics studies. A methodology is established for studying dynamics on multiple timescales in unilamellar membranes of more complex compositions.

Solid-state NMR spectroscopy to study protein-lipid interactions

The appropriate lipid environment is crucial for the proper function of membrane proteins. There is a tremendous variety of lipid molecules in the membrane and so far it is often unclear which component of the lipid matrix is essential for the function of a respective protein. Lipid molecules and proteins mutually influence each other; parameters such as acyl chain order, membrane thickness, membrane elasticity, permeability, lipid-domain and annulus formation are strongly modulated by proteins. More recent data also indicates that the influence of proteins goes beyond a single annulus of next-neighbor boundary lipids. Therefore, a mesoscopic approach to membrane lipid-protein interactions in terms of elastic membrane deformations has been developed. Solid-state NMR has greatly contributed to the understanding of lipid-protein interactions and the modern view of biological membranes. Methods that detect the influence of proteins on the membrane as well as direct lipid-protein interactions have been developed and are reviewed here. Examples for solid-state NMR studies on the interaction of Ras proteins, the antimicrobial peptide protegrin-1, the G protein-coupled receptor rhodopsin, and the K(+) channel KcsA are discussed. This article is part of a Special Issue entitled Tools to study lipid functions.

Carotenoids, versatile components of oxygenic photosynthesis

Carotenoids (CARs) are a group of pigments that perform several important physiological functions in all kingdoms of living organisms. CARs serve as protective agents, which are essential structural components of photosynthetic complexes and membranes, and they play an important role in the light harvesting mechanism of photosynthesizing plants and cyanobacteria. The protection against reactive oxygen species, realized by quenching of singlet oxygen and the excited states of photosensitizing molecules, as well as by the scavenging of free radicals, is one of the main biological functions of CARs. X-ray crystallographic localization of CARs revealed that they are present at functionally and structurally important sites of both the PSI and PSII reaction centers. Characterization of a CAR-less cyanobacterial mutant revealed that while the absence of CARs prevents the formation of PSII complexes, it does not abolish the assembly and function of PSI. CAR molecules assist in the formation of protein subunits of the photosynthetic complexes by gluing together their protein components. In addition to their aforementioned indispensable functions, CARs have a substantial role in the formation and maintenance of proper cellular architecture, and potentially also in the protection of the translational machinery under stress conditions.

Phospholipid vesicle stability and temporal variations in acyl chain organization

We report on our investigation of the long-term dimensional stability and acyl chain organization of unilamellar vesicles containing 1,2-dimyristoyl-sn-phosphatidylcholine (DMPC) and DMPC with cholesterol. Vesicles of 100% DMPC, 80% DMPC/20% cholesterol and 70% DMPC/30% cholesterol, formed by extrusion in aqueous buffer solution (pH 8) were shown to remain dimensionally stable for periods in excess of 600 h by dynamic light scattering (DLS) measurements. The rotational diffusion dynamics of perylene confined in the vesicle acyl chain region revealed structural evolution that was dependent on vesicle composition. Re-extrusion of the vesicles caused no changes in the average diameter or size distribution, but did give rise to diminished organization in the lipid acyl chain region for DMPC vesicles. Cholesterol-containing vesicles exhibited somewhat less pronounced change in organization on re-extrusion, suggesting a structurally mediating role for cholesterol.

Another Piece of the Membrane Puzzle: Extending Slipids Further

To be able to model complex biological membranes in a more realistic manner, the force field Slipids (Stockholm lipids) has been extended to include parameters for sphingomyelin (SM), phosphatidylglycerol (PG), phosphatidylserine (PS) lipids, and cholesterol. Since the parametrization scheme was faithful to the scheme used in previous editions of Slipids, all parameters are consistent and fully compatible. The results of careful validation of a number of key structural properties for one and two component lipid bilayers are in excellent agreement with experiments. Potentials of mean force for transferring water across binary mixtures of lipids and cholesterol were also computed in order to compare water permeability rates to experiments. In agreement with experimental and simulation studies, it was found that the permeability and partitioning of water is affected by cholesterol in lipid bilayers made of saturated lipids to the largest extent. With the extensions of Slipids presented here, it is now possible to study complex systems containing many different lipids and proteins in a fully atomistic resolution in the isothermic-isobaric (NPT) ensemble, which is the proper ensemble for membrane simulations.

Molecular dynamics simulation of hydrated DPPC monolayers using charge equilibration force fields

We present results of molecular dynamics simulations of a model DPPC-water monolayer using charge equilibration (CHEQ) force fields, which explicitly account for electronic polarization in a classical treatment of intermolecular interactions. The surface pressure, determined as the difference between the monolayer and pure water surface tensions at 323 K, is predicted to be 22.92 ±1.29 dyne/cm, just slightly below the broad range of experimental values reported for this system. The surface tension for the DPPC-water monolayer is predicted to be 42.35 ±1.16 dyne/cm, in close agreement with the experimentally determined value of 40.9 dyne/cm. This surface tension is also consistent with the value obtained from DPPC monolayer simulations using state-of-the-art nonpolarizable force fields. The current results of simulations predict a monolayer-water potential difference relative to the pure water-air interface of 0.64 ±0.02 Volts, an improved prediction compared to the fixed-charge CHARMM27 force field, yet still overestimating the experimental range of 0.3 to 0.45 Volts. As the charge equilibration model is a purely charge-based model for polarization, the current results suggest that explicitly modeled polarization effects can offer improvements in describing interfacial electrostatics in such systems.

An NMR database for simulations of membrane dynamics

Computational methods are powerful in capturing the results of experimental studies in terms of force fields that both explain and predict biological structures. Validation of molecular simulations requires comparison with experimental data to test and confirm computational predictions. Here we report a comprehensive database of NMR results for membrane phospholipids with interpretations intended to be accessible by non-NMR specialists. Experimental ¹³C-¹H and ²H NMR segmental order parameters (S(CH) or S(CD)) and spin-lattice (Zeeman) relaxation times (T(1Z)) are summarized in convenient tabular form for various saturated, unsaturated, and biological membrane phospholipids. Segmental order parameters give direct information about bilayer structural properties, including the area per lipid and volumetric hydrocarbon thickness. In addition, relaxation rates provide complementary information about molecular dynamics. Particular attention is paid to the magnetic field dependence (frequency dispersion) of the NMR relaxation rates in terms of various simplified power laws. Model-free reduction of the T(1Z) studies in terms of a power-law formalism shows that the relaxation rates for saturated phosphatidylcholines follow a single frequency-dispersive trend within the MHz regime. We show how analytical models can guide the continued development of atomistic and coarse-grained force fields. Our interpretation suggests that lipid diffusion and collective order fluctuations are implicitly governed by the viscoelastic nature of the liquid-crystalline ensemble. Collective bilayer excitations are emergent over mesoscopic length scales that fall between the molecular and bilayer dimensions, and are important for lipid organization and lipid-protein interactions. Future conceptual advances and theoretical reductions will foster understanding of biomembrane structural dynamics through a synergy of NMR measurements and molecular simulations.

On the possibility of the amphotericin B-sterol complex formation in cholesterol- and ergosterol-containing lipid bilayers: a molecular dynamics study

Amphotericin B (AmB) is a well-known membrane-active antibiotic that has been used to treat systemic fungal infections for more than 45 years. Therapeutic application of AmB is based on the fact that it is more active against ergosterol-containing membranes of fungal cells than against mammalian membranes with cholesterol. In this paper, we examine the hypothesis according to which the selectivity of the AmB's membrane action originates from its different ability to form the binary complexes with the relevant sterols. To this end, molecular dynamics simulations were performed for systems containing the preformed models of AmB/sterol complexes embedded in lipid bilayers containing either cholesterol or ergosterol. The initial structures of the studied binary associates were selected on the basis of a systematic scan of all possible mutual positions and orientations of the two molecules. The results obtained demonstrate that in general the complexes with ergosterol are more stable on the 100 ns time scale. Furthermore, on the basis of motional correlation analysis, taking into account the effects of lipid environment, we propose that, within the sterol-enriched liquid-ordered membrane phases, AmB molecules exhibit a greater tendency to bind ergosterol than cholesterol. The analysis of the interactions suggests that this affinity difference is of enthalpic origin and may arise from the considerable difference in the energy of the van der Waals interactions between AmB and the two types of sterols. Thus, our current results: (i) support the hypothesis that binary AmB/sterol complexes form within a lipid membrane and (ii) suggest that the higher toxicity may at least partly be attributed to the higher affinity of AmB for ergosterol than for cholesterol within a lipid membrane environment.

Volumetric Properties, Compressibilities and Volume Fluctuations in Phospholipid-Cholesterol Bilayers

We conducted detailed measurements of the apparent specific volume of dipalmitoylphosphatidylcholine (DPPC)–cholesterol mixtures in excess water as a function of pressure up to 70MPa (700bar) at 20, 38 and 50°C. The volumetric properties and the isothermal compressibility κ T of the lipid vesicles were determined at cholesterol concentrations, χchol, ranging up to 50 mol. The thermodynamic data are compared with other physico-chemical properties of phospholipid–cholesterol mixtures. Furthermore, the thermodynamic properties of the system are discussed in the light of the various T, χchol–phase diagrams and computer simulation studies published in the literature.

Retinal dynamics during light activation of rhodopsin revealed by solid-state NMR spectroscopy

Rhodopsin is a canonical member of class A of the G protein-coupled receptors (GPCRs) that are implicated in many of the drug interventions in humans and are of great pharmaceutical interest. The molecular mechanism of rhodopsin activation remains unknown as atomistic structural information for the active metarhodopsin II state is currently lacking. Solid-state (2)H NMR constitutes a powerful approach to study atomic-level dynamics of membrane proteins. In the present application, we describe how information is obtained about interactions of the retinal cofactor with rhodopsin that change with light activation of the photoreceptor. The retinal methyl groups play an important role in rhodopsin function by directing conformational changes upon transition into the active state. Site-specific (2)H labels have been introduced into the methyl groups of retinal and solid-state (2)H NMR methods applied to obtain order parameters and correlation times that quantify the mobility of the cofactor in the inactive dark state, as well as the cryotrapped metarhodopsin I and metarhodopsin II states. Analysis of the angular-dependent (2)H NMR line shapes for selectively deuterated methyl groups of rhodopsin in aligned membranes enables determination of the average ligand conformation within the binding pocket. The relaxation data suggest that the beta-ionone ring is not expelled from its hydrophobic pocket in the transition from the pre-activated metarhodopsin I to the active metarhodopsin II state. Rather, the major structural changes of the retinal cofactor occur already at the metarhodopsin I state in the activation process. The metarhodopsin I to metarhodopsin II transition involves mainly conformational changes of the protein within the membrane lipid bilayer rather than the ligand. The dynamics of the retinylidene methyl groups upon isomerization are explained by an activation mechanism involving cooperative rearrangements of extracellular loop E2 together with transmembrane helices H5 and H6. These activating movements are triggered by steric clashes of the isomerized all-trans retinal with the beta4 strand of the E2 loop and the side chains of Glu(122) and Trp(265) within the binding pocket. The solid-state (2)H NMR data are discussed with regard to the pathway of the energy flow in the receptor activation mechanism.

Charge equilibration force fields for lipid environments: applications to fully hydrated DPPC bilayers and DMPC-embedded gramicidin A

Polarizable force fields for lipid and solvent environments are used for molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer and gramicidin A (gA) dimer embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. The lipid bilayer is modeled using the CHARMM charge equilibration (CHEQ) polarizable force field for lipids and the TIP4P-FQ force field to represent solvent. For the DPPC bilayer system, results are compared to the same system simulated using the nonpolarizable CHARMM27r (C27r) force field and TIP3P water. Calculated atomic and electron density profiles, head group orientations as measured by the phosphorus-nitrogen vector orientation, and deuterium order parameters are found to be consistent with previous simulations and with experiment. The CHEQ model exhibits greater water penetration into the bilayer interior, as demonstrated by the potential of mean force calculated from the water density profile. This is a result of the variation of the water molecular dipole from 2.55 D in the bulk to 1.88 D in the interior. We discuss this finding in the context of previous studies (both simulation and experiment) that have investigated the extent of penetration of water into DPPC bilayers. We also discuss the effects of including explicit polarization on the water dipole moment variation as a function of distance from the bilayer. We show distributions of atomic charges over the course of the simulation since the CHEQ model allows the charges to fluctuate. We have calculated the interfacial dipole potential, which the CHEQ model predicts to be 0.95 V compared to 0.86 V as predicted by the C27r model. We also discuss dielectric permittivity profiles and the differences arising between the two models. We obtain bulk values of 72.77 for the CHEQ model (TIP4P-FQ water) and 91.22 for C27r (TIP3P), and values approaching unity in the membrane interior. Finally, we present results of simulations of gA embedded in a DMPC bilayer using the CHEQ model and discuss structural properties.

A comparative study of the COX-1 and COX-2 isozymes bound to lipid membranes

The monotopic proteins COX-1 and -2 in dimeric form bound to lipid bilayer membranes are studied using molecular dynamics simulations within an aqueous environment. The 25-ns simulations are performed for both isozymes with arachidonic acid bound in the cyclooxygenase sites. The interactions between the enzymes and the lipids are analyzed, providing insight into the attachment mechanism of monotopic proteins to membranes. Our study reveals some key differences between the two isozymes that include the orientations at which they sit on the surface of the membranes and the depths to which they embed within the membranes. The differences in membrane association of the isozymes indicate that they may integrate distinctively with the same membrane, and/or with different membranes or their lipid components. Our results indicate that arachidonic acid can be bound in the cyclooxygenase active site in distinct catalytically competent conformations that lead to certain hydroperoxy acids; and the arachidonic acid and/or cyclooxygenase sites undergo a conformational change which makes only one subunit of each homodimer catalytically active.

Molecular dynamics simulations of a DMPC bilayer using nonadditive interaction models

We present a polarizable force field based on the charge-equilibration formalism for molecular dynamics simulations of phospholipid bilayers. We discuss refinement of headgroup dihedral potential parameters to reproduce ab initio conformational energies of dimethylphosphate calculated at the MP2/cc-pVTZ level of theory. We also address the refinement of electrostatic and Lennard-Jones (van der Waals) parameters to reproduce ab initio polarizabilities and water interaction energies of dimethylphosphate and tetramethylammonium. We present results of molecular dynamics simulations of a solvated dimyristoylphosphatidylcholine bilayer using this polarizable force field as well as the nonpolarizable, fixed-charge CHARMM27 and CHARMM27r force fields for comparison. Calculated atomic and electron-density profiles, deuterium order parameters, and headgroup orientations are found to be consistent with previous simulations and with experiment. Polarizable interaction models for solvent and lipid exhibit greater water penetration into the lipid interior; this is due to the variation of water molecular dipole moment from a bulk value of 2.6 Debye to a value of 1.9 Debye in the membrane interior. The reduction in the electrostatic component of the desolvation free-energy penalty allows for greater water density. The surface dipole potential predicted by the polarizable model is 0.95 V compared to the value of 0.8 V based on nonpolarizable force-field calculations. Effects of inclusion of explicit polarization are discussed in relation to water dipole moment and varying charge distributions. Dielectric permittivity profiles for polarizable and nonpolarizable interactions exhibit subtle differences arising from the nature of the individual component parameterizations; for the polarizable force field, we obtain a bulk dielectric permittivity of 79, whereas the nonpolarizable force field plateaus at 97 (the value for pure TIP3P water). In the membrane interior, both models predict unit permittivities, with the polarizable models contributing from one to two more units due to the optical dielectric (high-frequency dipole fluctuations). This contribution is a step toward the continuing development of a CHARMM (Chemistry at Harvard Molecular Mechanics) polarizable force field for simulations of biomacromolecular systems.

Phase equilibria in DOPC/DPPC-d62/cholesterol mixtures

There is broad interest in the question of fluid-fluid phase coexistence in membranes, in particular, whether evidence for liquid-disordered (l(d))-liquid-ordered (l(o)) two-phase regions or membrane "rafts" can be found in natural membranes. In model membrane systems, such phase behavior is observed, and we have used deuterium nuclear magnetic resonance spectroscopy to map the phase boundaries of ternary mixtures containing 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), chain-perdeuterated 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC-d(62)), and cholesterol. For both this ternary model system and the binary DPPC-d(62)/cholesterol system, we present clear evidence for l(d)-l(o) two-phase coexistence. We have selected sample compositions to focus on this region of fluid-fluid phase coexistence and to determine its temperature and composition ranges. The deuterium nuclear magnetic resonance spectra for compositions near the l(d)-l(o) phase boundary at high cholesterol concentrations show evidence of exchange broadening or critical fluctuations in composition, similar to that reported by Vist and Davis. There appears to be a line of critical compositions ranging from 48 degrees C for a DOPC/DPPC-d(62)/cholesterol composition of 0:75:25, to approximately -8 degrees C for the composition 57:14:29. At temperatures below this two-phase region, there is a region of three-phase coexistence (l(d)-l(o)-gel). These results are collected and presented in terms of a partial ternary phase diagram that is consistent with previously reported results of Vist and Davis.

Solubilization of hydrophobic guest molecules in the monoolein discontinuous QL cubic mesophase and its soft nanoparticles

Hydrophobic bioactive guest molecules were solubilized in the discontinuous cubic mesophase (QL) of monoolein. Their effects on the mesophase structure and thermal behavior, and on the formation of soft nanoparticles upon dispersion of the bulk mesophase were studied. Four additives were analyzed. They were classified into two types based on their presumed location within the lipid bilayer and their influence on the phase behavior and structure. Differential scanning calorimetry (DSC), small-angle X-ray scattering (SAXS), polarized light microscopy, cryogenic-transmission electron microscopy (cryo-TEM), and dynamic light scattering (DLS) were used for the analysis. We found that carbamazepine and cholesterol (type I molecules) likely localize in the hydrophobic domains, but close to the hydrophobic-hydrophilic region. They induce strong perturbation to the mesophase packing by influencing both the order of the lipid acyl chains and interactions between lipid headgroups. This results in significant reduction of the phase transition enthalpy, and phase separation into lamellar and cubic mesophases above the maximum loading capacity. The inclusion of type I molecules in the mesophase also prevents the formation of soft nanoparticles with long-range internal order upon dispersion. In their presence, only vesicles or sponge-like nanoparticles form. Phytosterols and coenzyme Q10 (type II molecules) present only moderate effects. These molecules reside in the hydrophobic domains, where they cannot alter the lipid curvature or transform the QL mesophase into another phase. Therefore, above maximum loading, excess solubilizate precipitates in crystal forms. Moreover, when type II-loaded QL is dispersed, nanoparticles with long-range order and cubic symmetry (i.e., cubosomes) do form. A model for the growth of the ordered nanoparticles was developed from a series of intermediate structures identified by cryo-TEM. It proposes the development of the internal structure by fusion events between bilayer segments.

Membrane binding of lipidated Ras peptides and proteins--the structural point of view

Biological membranes are interesting interfaces, at which important biological processes occur. In addition to integral membrane proteins, a number of proteins bind to the membrane surface and associate with it. Posttranslational lipid modification is one important mechanism, by which soluble molecules develop a propensity towards the membrane and reversibly bind to it. Membrane binding by insertion of hydrophobic lipid moieties is relevant for up to 10% of all cellular proteins. A particular interesting lipid-modified protein is the small GTPase Ras, which plays a key role in cellular signal transduction. Until recently, the structural basis for membrane binding of Ras was not well-defined. However, with the advent of new synthesis techniques and the advancement of several biophysical methods, a number of structural and dynamical features about membrane binding of Ras proteins have been revealed. This review will summarize the chemical biology of Ras and discuss in more detail the biophysical and structural features of the membrane bound C-terminus of the protein.

Fluid mechanical matching of H+-ATP synthase subunit c-ring with lipid membranes revealed by 2H solid-state NMR

The F(1)F(o)-ATP synthase utilizes the transmembrane H(+) gradient for the synthesis of ATP. F(o) subunit c-ring plays a key role in transporting H(+) through F(o) in the membrane. We investigated the interactions of Escherichia coli subunit c with dimyristoylphosphatidylcholine (DMPC-d(54)) at lipid/protein ratios of 50:1 and 20:1 by means of (2)H-solid-state NMR. In the liquid-crystalline state of DMPC, the (2)H-NMR moment values and the order parameter (S(CD)) profile were little affected by the presence of subunit c, suggesting that the bilayer thickness in the liquid-crystalline state is matched to the transmembrane hydrophobic surface of subunit c. On the other hand, hydrophobic mismatch of subunit c with the lipid bilayer was observed in the gel state of DMPC. Moreover, the viscoelasticity represented by a square-law function of the (2)H-NMR relaxation was also little influenced by subunit c in the fluid phase, in contrast with flexible nonionic detergents or rigid additives. Thus, the hydrophobic matching of the lipid bilayer to subunit c involves at least two factors, the hydrophobic length and the fluid mechanical property. These findings may be important for the torque generation in the rotary catalytic mechanism of the F(1)F(o)-ATPse molecular motor.

Solid-state 2H NMR spectroscopy of retinal proteins in aligned membranes

Solid-state 2H NMR spectroscopy gives a powerful avenue to investigating the structures of ligands and cofactors bound to integral membrane proteins. For bacteriorhodopsin (bR) and rhodopsin, retinal was site-specifically labeled by deuteration of the methyl groups followed by regeneration of the apoprotein. 2H NMR studies of aligned membrane samples were conducted under conditions where rotational and translational diffusion of the protein were absent on the NMR time scale. The theoretical lineshape treatment involved a static axial distribution of rotating C-C2H3 groups about the local membrane frame, together with the static axial distribution of the local normal relative to the average normal. Simulation of solid-state 2H NMR lineshapes gave both the methyl group orientations and the alignment disorder (mosaic spread) of the membrane stack. The methyl bond orientations provided the angular restraints for structural analysis. In the case of bR the retinal chromophore is nearly planar in the dark- and all-trans light-adapted states, as well upon isomerization to 13-cis in the M state. The C13-methyl group at the "business end" of the chromophore changes its orientation to the membrane upon photon absorption, moving towards W182 and thus driving the proton pump in energy conservation. Moreover, rhodopsin was studied as a prototype for G protein-coupled receptors (GPCRs) implicated in many biological responses in humans. In contrast to bR, the retinal chromophore of rhodopsin has an 11-cis conformation and is highly twisted in the dark state. Three sites of interaction affect the torsional deformation of retinal, viz. the protonated Schiff base with its carboxylate counterion; the C9-methyl group of the polyene; and the beta-ionone ring within its hydrophobic pocket. For rhodopsin, the strain energy and dynamics of retinal as established by 2H NMR are implicated in substituent control of activation. Retinal is locked in a conformation that is twisted in the direction of the photoisomerization, which explains the dark stability of rhodopsin and allows for ultra-fast isomerization upon absorption of a photon. Torsional strain is relaxed in the meta I state that precedes subsequent receptor activation. Comparison of the two retinal proteins using solid-state 2H NMR is thus illuminating in terms of their different biological functions.

Interaction of cholesterol-like molecules in polyunsaturated phosphatidylcholine lipid bilayers as revealed by a self-consistent field theory

Cholesterol is one of the most abundant components in biological membranes. In this paper we apply a detailed state-of-the-art self-consistent field (SCF) theory to predict the influence of cholesterol-look-alikes in the bilayer composed of 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphatidylcholine (18:022:6omega3cis PC) lipids with a polyunsaturated 22:6 and a fully saturated 18:0 tail. The cholesterol-like molecule is composed of a hydroxyl group, a rigid chain fragment with length n segments and a branched semiflexible moiety with methylene side groups. We vary both the length of the rigid fragment in the cholesterol-look-alikes and their mole fraction in the tensionless bilayers. We find that these additives significantly increase the order of the saturated tails, but influence the conformational properties of the unsaturated tail much less. With increasing loading the bilayer thickness and the area available per PC head group increase. The hydroxyl group anchors close to the membrane-water interface, but with increasing loading the distribution of this polar group widens. The orientational order of the rigid part is high and we conclude that the cholesterol has significant mobility in the normal direction in the hydrophobic region of the bilayer indicating that one singly hydroxyl group is giving only a weak anchoring to the water-interface. Cholesterol-look-alikes increase the fluctuation of the tail ends and decrease the interdigitation of the tails. Several of our predictions correspond to molecular dynamics (MD) simulation results, but there are also important differences. Most notably the cholesterol-look-alikes can visit the membrane symmetry-plane more easily in SCF than in MD. Possible reasons for this are discussed.

Influence of Perfluorinated Compounds on the Properties of Model Lipid Membranes

The influence of selected perfluorinated compounds (PFCs), perfluorooctanoic acid (PFOA) or perfluorooctanesulfonic acid (PFOS), on the structure and organization of lipid membranes was investigated using model membranes-lipid monolayers and bilayers. The simplest model--a lipid monolayer--was studied at the air-water interface using the Langmuir-Blodgett technique with surface pressure and surface potential measurements. Lipid bilayers were characterized by NMR techniques and molecular dynamics simulations. Two phospholipids, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), characterized by different surface properties have been chosen as components of the model membranes. For a DPPC monolayer, a phase transition from the liquid-expanded state to the liquid-condensed state can be observed upon compression at room temperature, while a DMPC monolayer under the same conditions remains in the liquid-expanded state. For each of the two lipids, the presence of both PFOA and PFOS leads to the formation of a more fluidic layer at the air-water interface. Pulsed field gradient NMR measurements of the lateral diffusion coefficient (DL) of DMPC and PFOA in oriented bilayers reveal that, upon addition of PFOA to DMPC bilayers, DL of DMPC decreases for small amounts of PFOA, while larger additions produce an increased DL. The DL values of PFOA were found to be slightly larger than those for DMPC, probably as a consequence of the water solubility of PFOA. Furthermore, 31P and 2H NMR showed that the gel-liquid crystalline phase transition temperature decreased by the addition of PFOA for concentrations of 5 mol % and above, indicating a destabilizing effect of PFOA on the membranes. Deuterium order parameters of deuterated DMPC were found to increase slightly upon increasing the PFOA concentration. The monolayer experiments reveal that PFOS also penetrates slowly into already preformed lipid layers, leading to a change of their properties with time. These experimental observations are in qualitative agreement with the computational results obtained from the molecular dynamics simulations showing a slow migration of PFCs from the surrounding water phase into DPPC and DMPC bilayers.