Dynamics of lipid bilayers from comparative analysis of 2H and 13C nuclear magnetic resonance relaxation data as a function of frequency and temperature

Efficient calculation of orientation-dependent lipid dynamics from membrane simulations

Molecular dynamics simulations of lipid membranes have become increasingly impactful in biophysics because they offer atomistic resolution of structural fluctuations in relation to their functional outputs. Yet quantitative characterization of multiscale processes is a formidable challenge due to the distribution of motions that evade analysis of discrete simulation data. Here we investigate the efficient calculation of CH bond relaxation rates from membrane simulations. Widely used computational approaches offer numerical simplicity but fall short of capturing crucial aspects of the orientation dependence of the dynamics. To circumvent this problem, we introduced a robust framework based on liquid crystal theory which considers explicitly the CH bond motions with respect to the director axis (bilayer normal). Analysis of the orientation dependence of the dynamics shows excellent agreement with experiment, illustrating how the ordering potential affects the calculated relaxation rates. Furthermore, a fit-based resampling of the autocorrelation function of the bond fluctuations validates the new approach for low-temporal resolution data. The recovered relaxation rates indicate that at short timescales, both with and without cholesterol, the local motions of CH bonds describe the bilayer microviscosity and resemble liquid hydrocarbons. Our results establish the critical role of the orientational anisotropy in analysis of membrane simulations, explain fundamental aspects of lipid dynamics, and provide guidelines for extracting information that can be compared to experimental data.

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.

Dynamic Sorting of Mobile and Rigid Molecules in Biomembranes by Magic-Angle Spinning 13C NMR

Understanding the membrane dynamics of complex systems is essential to follow their function. As molecules in membranes can be in a rigid or mobile state depending on external (temperature, pressure) or internal (pH, domains, etc.) conditions, we propose an in-depth examination of NMR methods to filter highly mobile molecular parts from others that are in more restricted environments. We have thus developed a quantitative magic-angle spinning (MAS) ¹³C NMR approach coupled with cross-polarization (CP) and/or Insensitive Nuclei Enhanced by Polarization Transfer (INEPT) on rigid and fluid unlabeled model membranes. We demonstrate that INEPT can detect only very mobile lipid headgroups in gel (solid-ordered) phases; the remaining rigid parts are only detected with CP. A direct correlation is established between the normalized line intensity as obtained by CP and the C-H (C-D) order parameters measured by wide-line ²H NMR or extracted from molecular dynamics: I_CP/I_DP^eq ≈ 5|S_CH|, indicating that when the order is greater than 0.2-0.3 (maximum value of 0.5 for chain CH₂), only rigid parts can be filtered and detected using CP techniques. In very fluid (liquid-disordered) membranes, where there are many more active motions, both INEPT and CP detect resonances, with, however, a clear propensity of each technique to detect mobile and restricted molecular parts, respectively. Interestingly, the ¹³C NMR chemical shift of lipid hydrocarbon chains can be used to monitor order-disorder phase transitions and calculate the fraction of chain defects (rotamers) and the part of the transition enthalpy due to bond rotations (6-7 kJ·mol^-1 for dimyristolphosphatidylcholine, DMPC). Cholesterol-containing membranes (liquid-ordered phases) can be dynamically contrasted as the rigid-body sterol is mainly detected by the CP technique, with a contact time of 1 ms, and the phospholipid by INEPT. Our work opens up a straightforward, robust, and cost-effective route for the determination of membrane dynamics by taking advantage of well-resolved conventional ¹³C NMR experiments without the need of isotopic labeling.

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.

A solution NMR view of Lipidic Cubic Phases: Structure, dynamics, and beyond

Nuclear magnetic resonance (NMR) spectroscopy is well-established nowadays for the elucidation of the 3D structures of proteins and protein complexes, the evaluation of biomolecular dynamics with atomistic resolution across a range of time scales, the screening of drug candidates with site specificity, and for the quantitation of molecular translational diffusion. Lyotropic lipidic cubic phases (LCPs) are lipid bilayer-based materials with a complex geometry, formed via the spontaneous self-assembly of certain lipids in an aqueous environment at specific temperature ranges. LCPs have been successfully applied to the in meso crystallization of membrane proteins for structural studies by X-ray crystallography, and have also shown promising potential for serving as matrices for drug and nutrient delivery/release in vivo. The characterization of the structural and dynamics properties of LCPs is of significant interest for the application of these materials. Here we present a systematic review detailing the characterization of LCPs by solution NMR. Using LCPs formed by monoolein (MO) as an example, various aspects of LCPs readily accessible by solution NMR are covered, including spectral perturbation in the presence of additives, quantification of hydration levels, ¹³C relaxation-based measurements for studying atom-specific dynamics along the MO hydrocarbon chain, PGSE NMR measurement of translational diffusion and its correlation with release profiles, and the encapsulation of soluble proteins in LCPs. A brief discussion of future perspectives for the characterization of LCPs by solution NMR is also presented.

Phospholipid headgroups govern area per lipid and emergent elastic properties of bilayers

Phospholipid bilayers are liquid-crystalline materials whose intermolecular interactions at mesoscopic length scales have key roles in the emergence of membrane physical properties. Here we investigated the combined effects of phospholipid polar headgroups and acyl chains on biophysical functions of membranes with solid-state 2H NMR spectroscopy. We compared the structural and dynamic properties of phosphatidylethanolamine and phosphatidylcholine with perdeuterated acyl chains in the solid-ordered (so) and liquid-disordered (ld) phases. Our analysis of spectral lineshapes of 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE-d62) and 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphocholine (DPPC-d62) in the so (gel) phase indicated an all-trans rotating chain structure for both lipids. Greater segmental order parameters (SCD) were observed in the ld (liquid-crystalline) phase for DPPE-d62 than for DPPC-d62 membranes, while their mixtures had intermediate values irrespective of the deuterated lipid type. Our results suggest the SCD profiles of the acyl chains are governed by methylation of the headgroups and are averaged over the entire system. Variations in the acyl chain molecular dynamics were further investigated by spin-lattice (R1Z) and quadrupolar-order relaxation (R1Q) measurements. The two acyl-perdeuterated lipids showed distinct differences in relaxation behavior as a function of the order parameter. The R1Z rates had a square-law dependence on SCD, implying collective mesoscopic dynamics, with a higher bending rigidity for DPPE-d62 than for DPPC-d62 lipids. Remodeling of lipid average and dynamic properties by methylation of the headgroups thus provides a mechanism to control the actions of peptides and proteins in biomembranes.

The dynamical Matryoshka model: 3. Diffusive nature of the atomic motions contained in a new dynamical model for deciphering local lipid dynamics

In accompanying papers [Bicout et al., BioRxiv https://doi.org/10.1101/2021.09.21.461198 (2021); Cissé et al., BioRxiv https://doi.org/10.1101/2022.03.30.486370 (2022)], a new model called Matryoshka model has been proposed to describe the geometry of atomic motions in phospholipid molecules in bilayers and multilamellar vesicles based on their quasielastic neutron scattering (QENS) spectra. Here, in order to characterize the relaxational aspects of this model, the energy widths of the QENS spectra of the samples were analyzed first in a model-free way. The spectra were decomposed into three Lorentzian functions, which are classified as slow, intermediate, and fast motions depending on their widths. The analysis provides the diffusion coefficients, residence times, and geometrical parameters for the three classes of motions. The results corroborate the parameter values such as the amplitudes and the mobile fractions of atomic motions obtained by the application of the Matryoshka model to the same samples. Since the current analysis was carried out independently of the development of the Matryoshka model, the present results enhance the validity of the model. The model will serve as a powerful tool to decipher the dynamics of lipid molecules not only in model systems, but also in more complex systems such as mixtures of different kinds of lipids or natural cell membranes.

A method to construct the dynamic landscape of a bio-membrane with experiment and simulation

Biomolecular function is based on a complex hierarchy of molecular motions. While biophysical methods can reveal details of specific motions, a concept for the comprehensive description of molecular dynamics over a wide range of correlation times has been unattainable. Here, we report an approach to construct the dynamic landscape of biomolecules, which describes the aggregate influence of multiple motions acting on various timescales and on multiple positions in the molecule. To this end, we use 13C NMR relaxation and molecular dynamics simulation data for the characterization of fully hydrated palmitoyl-oleoyl-phosphatidylcholine bilayers. We combine dynamics detector methodology with a new frame analysis of motion that yields site-specific amplitudes of motion, separated both by type and timescale of motion. In this study, we show that this separation allows the detailed description of the dynamic landscape, which yields vast differences in motional amplitudes and correlation times depending on molecular position.

Measuring the bending modulus of lipid bilayers with cholesterol

The effect of cholesterol on the bending modulus of DOPC lipid bilayers has become a controversial topic that has implications for methods of measuring the bending modulus. Recent results using neutron spin echo and nuclear magnetic resonance relaxation methods that involve linear transport properties have conflicted with earlier results from purely equilibrium experiments that do not involve linear transport properties. A general discussion indicates how one can be misled by data obtained by methods that involve linear transport properties. It is then shown specifically how the recent neutron spin echo results can be interpreted to agree with the earlier purely equilibrium experimental results, thereby resolving that conflict. Regarding the nuclear magnetic resonance relaxation method, it is noted that current interpretation of the data is unclear regarding the identity of the modulus that is involved, and an alternative interpretation is explored that does not disagree with the results of the equilibrium experiments.

How cholesterol stiffens unsaturated lipid membranes

Cholesterol is an integral component of eukaryotic cell membranes and a key molecule in controlling membrane fluidity, organization, and other physicochemical parameters. It also plays a regulatory function in antibiotic drug resistance and the immune response of cells against viruses, by stabilizing the membrane against structural damage. While it is well understood that, structurally, cholesterol exhibits a densification effect on fluid lipid membranes, its effects on membrane bending rigidity are assumed to be nonuniversal; i.e., cholesterol stiffens saturated lipid membranes, but has no stiffening effect on membranes populated by unsaturated lipids, such as 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). This observation presents a clear challenge to structure-property relationships and to our understanding of cholesterol-mediated biological functions. Here, using a comprehensive approach-combining neutron spin-echo (NSE) spectroscopy, solid-state deuterium NMR (²H NMR) spectroscopy, and molecular dynamics (MD) simulations-we report that cholesterol locally increases the bending rigidity of DOPC membranes, similar to saturated membranes, by increasing the bilayer's packing density. All three techniques, inherently sensitive to mesoscale bending fluctuations, show up to a threefold increase in effective bending rigidity with increasing cholesterol content approaching a mole fraction of 50%. Our observations are in good agreement with the known effects of cholesterol on the area-compressibility modulus and membrane structure, reaffirming membrane structure-property relationships. The current findings point to a scale-dependent manifestation of membrane properties, highlighting the need to reassess cholesterol's role in controlling membrane bending rigidity over mesoscopic length and time scales of important biological functions, such as viral budding and lipid-protein interactions.

Quasi-Elastic Neutron Scattering Studies on Hydration Water in Phospholipid Membranes

The dynamic behavior of hydration water in phospholipid membranes has been investigated to understand the relationship between water and biological molecules using various experimental techniques. Quasi-elastic neutron scattering (QENS) is an effective method for this purpose because the dynamic behaviors of both water and lipid molecules could be identified by using selective deuteration. In addition, the measurable ranges from the 10⁻¹² to 10⁻⁹ s time scale and the 10⁻¹¹ to 10⁻⁸ m length scale are suitable to investigate the slowing down of water molecules due to their interaction with lipid membranes. In this mini-review, QENS experiments on the dynamic behavior of hydration water molecules in neighboring phospholipid membranes are summarized.

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.

Phospholipid Bilayer Softening Due to Hydrophobic Gold Nanoparticle Inclusions

Liposome-nanoparticle assemblies (LNAs) are vital in the context of novel targeted drug-delivery systems, in addition to investigating nanoparticle-lipid bilayer interactions. Quantifying membrane structural properties and dynamics in presence of nanoparticle inclusions provides a simple model to elucidate nanoparticle effects on membrane biophysical properties. We present experimental evidences of bilayer softening due to small hydrophobic gold nanoparticle inclusions. LNA structure has been investigated by a combination of cryo-transmission electron microscopy, dynamic light scattering, and small-angle neutron scattering. Neutron spin echo spectroscopy demonstrated a remarkable ∼15% bending modulus decrease for LNAs relative to pure liposomes. Clear dependence of bending modulus on gold nanoparticle diameter and concentration was observed from our observations. Our findings point toward local bilayer fluidization by nanoparticle inclusions leading to an overall bilayer softening. These findings add valuable information to liposomal drug-delivery vehicle design and membrane biophysics research.

Lipid- and Cholesterol-Mediated Time-Scale-Specific Modulation of the Outer Membrane Protein X Dynamics in Lipid Bilayers

Membrane protein function fundamentally depends on lipid-bilayer fluidity and the composition of the biological membrane. Although dynamic interdependencies between membrane proteins and the surrounding lipids are suspected, a detailed description is still missing. To uncover lipid-modulated membrane protein backbone dynamics, time-scale-specific NMR relaxation experiments with residue-resolution were recorded. The data revealed that lipid order, modified either biochemically or biophysically, changes the dynamics of the immersed membrane protein in a specific and time-scale-dependent manner. A temperature-dependent dynamics analysis furthermore suggests a direct coupling between lipid and protein dynamics in the picosecond-nanosecond, microsecond, and millisecond time scales, caused by the lipid's trans-gauche isomerization, the segmental and rotational motion of lipids, and the fluidity of the lipid phase, respectively. These observations provide evidence of a direct modulatory capability of the membrane to regulate protein function through lipid dynamics ranging from picoseconds to milliseconds.

Relaxation dynamics of saturated and unsaturated oriented lipid bilayers

We present experimental measurements and analysis of the dynamics and the phase behaviour of saturated DMPC and unsaturated DOPC oriented multi-lamellar bilayers. Elastic and inelastic neutron scattering were used to directly probe the dynamical processes of these membrane systems on time and length scales relevant to the internal and localized motion of lipid monomers. Mobility in this regime can be informative in elucidating the local interactions responsible for material properties of these fluid lipid systems. DMPC and DOPC are structurally similar in terms of their membrane hydrophobic thickness; however, they exhibit different mechanical properties in terms of both elastic compressibility and bending moduli. The analyses suggest that the constraint imposed by the double bonds in DOPC acyl chains restricts atomic motion in both liquid and gel phases compared to DMPC. We discuss applications of molecular dynamics to further elucidate the atomic details of the dynamical processes. Such an understanding may suggest how membrane properties can be tuned using a variety of different lipid species.

Thermodynamics of lipid multi-lamellar vesicles in presence of sterols at high hydrostatic pressure

We compared the effect of cholesterol at different concentration on the phase behaviour of DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) multilamellar vesicles. We used pressure perturbation differential scanning calorimetry (PPC) that studies a system on the whole by giving access to relevant thermodynamic quantities, and elastic incoherent neutron scattering (EINS) that probes local motions of a system at the atomic level by allowing extraction of dynamical parameters. PPC revealed that the volume expansion coefficient of DMPC and DMPC/Cholesterol samples with 13 and 25 mol% cholesterol is a linear function of the heat capacity measured by differential scanning calorimetry. Neutron backscattering spectroscopy showed that the mean square displacements of H atoms do exhibit an increase with temperature and a decrease under increasing pressure. Cholesterol added at concentrations of 25 and 50 mol% led to suppression of the main phase transition. Taking advantage of these results, the present study aims (i) to show that calorimetry and EINS using the Bicout and Zaccai model equally permit to get access to thermodynamic quantities characterizing pure DMPC and DMPC/cholesterol mixtures, thus directly confirming the theoretical method, and (ii) to validate our approach as function of temperature and of pressure, as both are equally important and complementary thermodynamic variables.

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.

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.

Area per lipid and cholesterol interactions in membranes from separated local-field (13)C NMR spectroscopy

Investigations of lipid membranes using NMR spectroscopy generally require isotopic labeling, often precluding structural studies of complex lipid systems. Solid-state (13)C magic-angle spinning NMR spectroscopy at natural isotopic abundance gives site-specific structural information that can aid in the characterization of complex biomembranes. Using the separated local-field experiment DROSS, we resolved (13)C-(1)H residual dipolar couplings that were interpreted with a statistical mean-torque model. Liquid-disordered and liquid-ordered phases were characterized according to membrane thickness and average cross-sectional area per lipid. Knowledge of such structural parameters is vital for molecular dynamics simulations, and provides information about the balance of forces in membrane lipid bilayers. Experiments were conducted with both phosphatidylcholine (dimyristoylphosphatidylcholine (DMPC) and palmitoyloleoylphosphatidylcholine (POPC)) and egg-yolk sphingomyelin (EYSM) lipids, and allowed us to extract segmental order parameters from the (13)C-(1)H residual dipolar couplings. Order parameters were used to calculate membrane structural quantities, including the area per lipid and bilayer thickness. Relative to POPC, EYSM is more ordered in the ld phase and experiences less structural perturbation upon adding 50% cholesterol to form the lo phase. The loss of configurational entropy is smaller for EYSM than for POPC, thus favoring its interaction with cholesterol in raftlike lipid systems. Our studies show that solid-state (13)C NMR spectroscopy is applicable to investigations of complex lipids and makes it possible to obtain structural parameters for biomembrane systems where isotope labeling may be prohibitive.

Multi-component modeling of quasielastic neutron scattering from phospholipid membranes

We investigated molecular motions in the 0.3-350 ps time range of D2O-hydrated bilayers of 1-palmitoyl-oleoyl-sn-glycero-phosphocholine and 1,2-dimyristoyl-sn-glycero-phosphocholine in the liquid phase by quasielastic neutron scattering. Model analysis of sets of spectra covering scale lengths from 4.8 to 30 Å revealed the presence of three types of motion taking place on well-separated time scales: (i) slow diffusion of the whole phospholipid molecules in a confined cylindrical region; (ii) conformational motion of the phospholipid chains; and (iii) fast uniaxial rotation of the hydrogen atoms around their carbon atoms. Based on theoretical models for the hydrogen dynamics in phospholipids, the spatial extent of these motions was analysed in detail and the results were compared with existing literature data. The complex dynamics of protons was described in terms of elemental dynamical processes involving different parts of the phospholipid chain on whose motions the hydrogen atoms ride.

Solid state NMR of lipid model membranes

In this chapter we describe the use of solid state nuclear magnetic spectroscopy to study the structure of lyotropic phases and lipid model membranes and show its ability to probe, site specifically, at a sub-Ångstrom resolution. Here, we demonstrate the immense versatility of the technique and its ability to provide information on the different liquid crystalline phases present. A multinuclear for example (31)P, (1)H, and (13)C approach is able to elucidate both the structure and dynamics over a wide variety of timescales. This coupled with a non-perturbing label (2)H is able to provide information such as the order parameters for a wide variety of different liquid phases.

Investigation of rhodopsin dynamics in its signaling state by solid-state deuterium NMR spectroscopy

Site-directed deuterium NMR spectroscopy is a valuable tool to study the structural dynamics of biomolecules in cases where solution NMR is inapplicable. Solid-state (2)H NMR spectral studies of aligned membrane samples of rhodopsin with selectively labeled retinal provide information on structural changes of the chromophore in different protein states. Moreover (2)H NMR relaxation time measurements allow one to study the dynamics of the ligand during the transition from the inactive to the active state. Here we describe the methodological aspects of solid-state (2)H NMR spectroscopy for functional studies of rhodopsin, with an emphasis on the dynamics of the retinal cofactor. We provide complete protocols for the preparation of NMR samples of rhodopsin with 11-cis-retinal selectively deuterated at the methyl groups in aligned membranes. In addition we review optimized conditions for trapping the rhodopsin photointermediates; and we address the challenging problem of trapping the signaling state of rhodopsin in aligned membrane films.

High-resolution NMR reveals secondary structure and folding of amino acid transporter from outer chloroplast membrane

Solving high-resolution structures for membrane proteins continues to be a daunting challenge in the structural biology community. In this study we report our high-resolution NMR results for a transmembrane protein, outer envelope protein of molar mass 16 kDa (OEP16), an amino acid transporter from the outer membrane of chloroplasts. Three-dimensional, high-resolution NMR experiments on the ¹³C, ¹⁵N, ²H-triply-labeled protein were used to assign protein backbone resonances and to obtain secondary structure information. The results yield over 95% assignment of N, H_N, CO, C_α, and C_β chemical shifts, which is essential for obtaining a high resolution structure from NMR data. Chemical shift analysis from the assignment data reveals experimental evidence for the first time on the location of the secondary structure elements on a per residue basis. In addition T_1Z and T₂ relaxation experiments were performed in order to better understand the protein dynamics. Arginine titration experiments yield an insight into the amino acid residues responsible for protein transporter function. The results provide the necessary basis for high-resolution structural determination of this important plant membrane protein.

Curvature forces in membrane lipid-protein interactions

Membrane biochemists are becoming increasingly aware of the role of lipid-protein interactions in diverse cellular functions. This review describes how conformational changes in membrane proteins, involving folding, stability, and membrane shape transitions, potentially involve elastic remodeling of the lipid bilayer. Evidence suggests that membrane lipids affect proteins through interactions of a relatively long-range nature, extending beyond a single annulus of next-neighbor boundary lipids. It is assumed the distance scale of the forces is large compared to the molecular range of action. Application of the theory of elasticity to flexible soft surfaces derives from classical physics and explains the polymorphism of both detergents and membrane phospholipids. A flexible surface model (FSM) describes the balance of curvature and hydrophobic forces in lipid-protein interactions. Chemically nonspecific properties of the lipid bilayer modulate the conformational energetics of membrane proteins. The new biomembrane model challenges the standard model (the fluid mosaic model) found in biochemistry texts. The idea of a curvature force field based on data first introduced for rhodopsin gives a bridge between theory and experiment. Influences of bilayer thickness, nonlamellar-forming lipids, detergents, and osmotic stress are all explained by the FSM. An increased awareness of curvature forces suggests that research will accelerate as structural biology becomes more closely entwined with the physical chemistry of lipids in explaining membrane structure and function.

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.

Applications of neutron and X-ray scattering to the study of biologically relevant model membranes

Scattering techniques, in particular electron, neutron and X-ray scattering have played a major role in elucidating the static and dynamic structure of biologically relevant membranes. Importantly, neutron and X-ray scattering have evolved to address new sample preparations that better mimic biological membranes. In this review, we will report on some of the latest model membrane results, and the neutron and X-ray techniques that were used to obtain them.

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.

Collective and noncollective models of NMR relaxation in lipid vesicles and multilayers

NMR (13)C spin lattice relaxation (1/T(1)) rates of dipalmitoylphosphatidylcholine (DPPC) bilayers obtained from molecular dynamics simulations of 72 and 288 lipids are compared with each other, with experimental values from large liposomes obtained by magic angle spinning, and with previously published experimental data from small vesicles. The experimental results for multilayers and vesicles at the same frequencies differ only slightly. The simulation results indicate that T(1) relaxation in the 15.1 to 201.2 MHz carbon frequency range and up to 100 A length scale is dominated by fast isomerizations and slower lipid wobble (D perpendicular approximately 2.5 x 10(8) s(-1)). Rotational diffusion about the lipid long axis (described by D(parallel)) does not make a substantial contribution to the T(1). Modifications to the acyl chain torsional potential energy function used for the simulations substantially improve agreement with experiment.

X-ray scattering and solid-state deuterium nuclear magnetic resonance probes of structural fluctuations in lipid membranes

Molecular fluctuations are a dominant feature of biomembranes. Cellular functions might rely on these properties in ways yet to be determined. This expectation is suggested by the fact that membrane deformation and rigidity, which govern molecular fluctuations, have been implicated in a number of cellular functions. However, fluctuations are more challenging to measure than average structures, which partially explain the small number of dedicated studies. Here, it is shown that two accessible laboratory methods, small-angle X-ray scattering and solid-state deuterium nuclear magnetic resonance (NMR), can be used as complementary probes of structural fluctuations in lipid membranes. In the case of X-ray scattering, membrane undulations give rise to logarithmically varying positional correlations that generate scattering peaks with long (power-law) tails. In the case of 2H NMR spectroscopy, fluctuations in the magnetic-coupling energies resulting from molecular motions cause relaxation among the various spin energy levels, and yield a powerful probe of orientational fluctuations of the lipid molecules. A unified interpretation of the combined scattering and 2H NMR data is provided by a liquid-crystalline membrane deformation model. The importance of this approach is that it is possible to utilize a microscopic model for positional and orientational observables to calculate bulk material properties of liquid-crystalline systems.

Flexibility of ras lipid modifications studied by 2H solid-state NMR and molecular dynamics simulations

Human posttranslationally modified N-ras oncogenes are known to be implicated in numerous human cancers. Here, we applied a combination of experimental and computational techniques to determine structural and dynamical details of the lipid chain modifications of an N-ras heptapeptide in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes. Experimentally, 2H NMR spectroscopy was used to study oriented membranes that incorporated ras heptapeptides with two covalently attached perdeuterated hexadecyl chains. Atomistic molecular dynamics simulations of the same system were carried out over 100 ns including 60 DMPC and 4 ras molecules. Several structural and dynamical experimental parameters could be directly compared to the simulation. Experimental and simulated 2H NMR order parameters for the methylene groups of the ras lipid chains exhibited a systematic difference attributable to the absence of collective motions in the simulation and to geometrical effects. In contrast, experimental 2H NMR spin-lattice relaxation rates for Zeeman order were well reproduced in the simulation. The lack of slower collective motions in the simulation did not appreciably influence the relaxation rates at a Larmor frequency of 115.1 MHz. The experimental angular dependence of the 2H NMR relaxation rates with respect to the external magnetic field was also relatively well simulated. These relaxation rates showed a weak angular dependence, suggesting that the lipid modifications of ras are very flexible and highly mobile in agreement with the low order parameters. To quantify these results, the angular dependence of the 2H relaxation rates was calculated by an analytical model considering both molecular and collective motions. Peptide dynamics in the membrane could be modeled by an anisotropic diffusion tensor with principal values of Dparallel=2.1x10(9) s(-1) and Dperpendicular=4.5x10(5) s(-1). A viscoelastic fitting parameter describing the membrane elasticity, viscosity, and temperature was found to be relatively similar for the ras peptide and the DMPC host matrix. Large motional amplitudes and relatively short correlation times facilitate mixing and dispersal with the lipid bilayer matrix, with implications for the role of the full-length ras protein in signal transduction and oncogenesis.

Nanosecond molecular relaxations in lipid bilayers studied by high energy-resolution neutron scattering and in situ diffraction

We report a high energy-resolution neutron backscattering study to investigate slow motions on nanosecond time scales in highly oriented solid-supported phospholipid bilayers of the model system deuterated 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine, hydrated with heavy water. Wave-vector-resolved quasielastic neutron scattering is used to determine relaxation times tau , which can be associated with different molecular components, i.e., the lipid acyl chains and the interstitial water molecules in the different phases of the model membrane system. The inelastic data are complemented by both energy-resolved and energy-integrated in situ diffraction. From a combined analysis of the inelastic data in the energy and time domains, the corresponding character of the relaxation, i.e., the exponent of the exponential decay, is also determined. From this analysis we quantify two relaxation processes. We associate the fast relaxation with translational diffusion of lipid and water molecules while the slow process likely stems from collective dynamics.

Dynamics in atomistic simulations of phospholipid membranes: Nuclear magnetic resonance relaxation rates and lateral diffusion

It is shown that a long, near microsecond, atomistic simulation can shed some light upon the dynamical processes occurring in a lipid bilayer. The analysis focuses on reorientational dynamics of the chains and lateral diffusion of lipids. It is shown that the reorientational correlation functions exhibits an algebraic decay (rather than exponential) for several orders of magnitude in time. The calculated nuclear magnetic resonance relaxation rates agree with experiments for carbons at the C7 position while there are some differences for C3. Lateral diffusion can be divided into two stages. In a first stage occurring at short times, t<5 ns, the center of mass of the lipid moves due to conformational changes of the chains while the headgroup position remains relatively fixed. In this stage, the center of mass can move up to approximately 0.8 nm. The fitted short-time diffusion coefficient is D(1)=13 x 10(-7) cm(2) s(-1) On a longer time scale, the diffusion coefficient becomes D(2)=0.79 x 10(-7) cm(2) s(-1).

Dynamics of Alkyl Ammonium Intercalants within Organically Modified Montmorillonite: Dielectric Relaxation and Ionic Conductivity

The low-frequency (0.01 Hz-10 MHz) dynamic characteristics of alkyl quaternary ammonium exchanged montmorillonite (SC20A) were investigated to determine the correlation between temperature-dependent changes in the interlayer structure and collective mobility of the surfactant. From 25 to 165 degrees C, SC20A exhibits two interlayer transitions, one ascribed to the melting of the intercalated alkyl chains of the surfactant (20-40 degrees C) and another associated with an abrupt decrease in the interlayer's coefficient of thermal expansion (100 degrees C). For this temperature range, the excess surfactant and residual electrolytes present in commercially manufactured SC20A enhance the direct current conductivity and increase low-frequency space-charge polarization, which is believed to occur across percolation paths established by the surfaces of the SC20A crystallites. In contrast, a higher-frequency relaxation, which was less sensitive to process history and impurity content, is ascribed to relaxation within the interlayer at the surfactant-aluminosilicate interface electrostatic couple. The temperature dependence of these dielectric relaxations indicated a drastic increase in mobility as the interlayer organic phase transitions from static and glasslike into molten and mobile. Overall, SC20A displayed features of alternating current universality, including time-temperature superposition, common in other types of disordered ion-conducting media. The presence of long-range transport and its sensitivity to low amounts of impurities imply that from a dynamic perspective the local environment of the surfactants are substantially diverse and a minority fraction, such as at the edge of the crystallite (gallery and aluminosilicate layer), may dominate the lower-frequency dielectric response.

Picosecond dynamics of a membrane protein revealed by 2D IR

Fast protein dynamics can be missed with techniques that have relatively slow observation times. Using 2D IR spectroscopy and isotope labeling, we have probed the rapid, picosecond dynamics of a membrane protein in its native environment. By measuring the homogeneous and inhomogeneous IR linewidths of 11 amide I modes (backbone carbonyl stretch), we have captured the structural distributions and dynamics of the CD3zeta protein along its transmembrane segment that are lost with slower time-scale techniques. We find that the homogeneous lifetimes and population relaxation times are the same for almost all of the residues. In contrast, the inhomogeneous linewidths vary significantly with the largest inhomogeneous distribution occurring for residues near the N terminus and the narrowest near the center. This behavior is highly consistent with a recently reported experimental model of the protein and water accessibility as observed by molecular dynamics simulations. The data support the proposed CD3zeta peptide structure, and the simulations point to the structural disorder of water and lipid head-groups as the main source of inhomogeneous broadening. Taken together, this rigorous analysis of the vibrational dynamics of a membrane peptide provides experimental insight into a time regime of motions that has so far been largely unexplored.

Molecular motions in lipid bilayers studied by the neutron backscattering technique

We report a high energy-resolution neutron backscattering study to investigate slow motions on nanosecond time scales in highly oriented solid supported phospholipid bilayers of the model system DMPC-d54 (deuterated 1,2-dimyristoyl-sn-glycero-3-phoshatidylcholine), hydrated with heavy water. This technique allows to discriminate the onset of mobility at different length scales for the different molecular components, as, e.g., the lipid acyl-chains and the hydration water in between the membrane stacks, respectively, and provides a benchmark test regarding the feasibility of neutron backscattering investigations on these sample systems. We discuss freezing of the lipid acyl-chains, as observed by this technique, and observe a second freezing transition which we attribute to the hydration water.

Packing and Viscoelasticity of Polyunsaturated ω-3 and ω-6 Lipid Bilayers as Seen by2H NMR and X-ray Diffraction

Polyunsaturated phospholipids of the omega-3 and omega-6 classes play key roles in cellular functions, yet their mechanisms of biological action are still a matter of debate. Using deuterium ((2)H) NMR spectroscopy and small-angle X-ray diffraction, we show how membrane properties are modified by docosahexaenoic (DHA; 22:6) and arachidonic (AA; 20:4) acyl chains of the omega-3 and the omega-6 families, respectively. Structural and dynamical differences due to polyunsaturation are evident in both the ordered and disordered phases of mixed-chain (16:0)(22:6)PC and (16:0)(20:4)PC bilayers. Due to the lower chain melting temperature, the omega-6 AA bilayer is more disordered in the fluid (L(alpha)) state than the omega-3 DHA bilayer; it is thinner with a larger area per lipid. The thermal hysteresis observed for the DHA bilayer may represent the influences of angle-iron conformers in the gel state and back-bended, hairpinlike conformers in the fluid state, consistent with molecular dynamics studies. Interpretation of the (2)H NMR order profiles of (16:0-d(31))(22:6)PC and (16:0-d(31))(20:4)PC together with X-ray electron density profiles reveals an uneven distribution of mass; i.e., the sn-1 saturated chain is displaced toward the membrane center, whereas the sn-2 polyunsaturated chain is shifted toward the bilayer aqueous interface. Moreover, the (2)H NMR relaxation rates are increased by the presence of omega-6 AA chains compared to omega-3 DHA chains. When evaluated at the same amplitude of motion, relaxation parameters give a naturally calibrated scale for comparison of fluid lipid bilayers. Within this framework, polyunsaturated bilayers are relatively soft to bending and area fluctuations on the mesoscale approaching molecular dimensions. Significant differences are evident in the viscoelastic properties of the omega-3 and omega-6 bilayers, a possibly biologically relevant feature that distinguishes between the two phospholipid classes.

Collective dynamics of lipid membranes studied by inelastic neutron scattering

We have studied the collective short wavelength dynamics in deuterated 1,2-dimyristoyl-sn-glycero-3-phoshatidylcholine (DMPC) bilayers by inelastic neutron scattering. The corresponding dispersion relation variant Planck's over 2pi omega(Q) is presented for the gel and the fluid phase of this model system. The temperature dependence of the inelastic excitations indicates a phase coexistence between the two phases over a broad range and leads to a different assignment of excitations from that reported in a preceding inelastic x-ray scattering study [Phys. Rev. Lett. 86, 740 (2001)]]. As a consequence, we find that the minimum in the dispersion relation is actually deeper in the gel than in the fluid phase. Finally, we can clearly identify an additional nondispersive (optical) mode predicted by molecular dynamics simulations [Phys. Rev. Lett. 87, 238101 (2001)]].

Structure and dynamics of sphingomyelin bilayer: insight gained through systematic comparison to phosphatidylcholine

Sphingomyelin, one of the main lipid components of biological membranes, is actively involved in various cellular processes such as protein trafficking and signal transduction. In particular, specific lateral domains enriched in sphingomyelin and cholesterol have been proposed to play an important functional role in biomembranes, although their precise characteristics have remained unclear. A thorough understanding of the functional role of membranes requires detailed knowledge of their individual lipid components. Here, we employ molecular dynamics simulations to conduct a systematic comparison of a palmitoylsphingomyelin (PSM, 16:0-SM) bilayer with a membrane that comprises dipalmitoylphosphatidylcholine (DPPC) above the main phase transition temperature. We clarify atomic-scale properties that are specific to sphingomyelin due to its sphingosine moiety, and further discuss their implications for SM-rich membranes. We find that PSM bilayers, and in particular the dynamics of PSM systems, are distinctly different from those of a DPPC bilayer. When compared with DPPC, the strong hydrogen bonding properties characteristic to PSM are observed to lead to considerable structural changes in the polar headgroup and interface regions. The strong ordering of PSM acyl chains and specific ordering effects in the vicinity of a PSM-water interface reflect this issue clearly. The sphingosine moiety and related hydrogen bonding further play a crucial role in the dynamics of PSM bilayers, as most dynamic properties, such as lateral and rotational diffusion, are strongly suppressed. This is most evident in the rotational motion characterized by spin-lattice relaxation times and the decay of hydrogen bond autocorrelation functions that are expected to be important in complexation of SM with other lipids in many-component bilayers. A thorough understanding of SM bilayers would greatly benefit from nuclear magnetic resonance experiments for acyl chain ordering and dynamics, allowing full comparison of these simulations to experiments.

The simulation of 31P NMR line shapes of lipid bilayers using an analytically soluble model

Van Faassen's method for obtaining an explicit solution to a first order stochastic differential equation is applied to the simulation of 31P NMR line shapes of unoriented phospholipid bilayers in the Lalpha phase and of oriented bilayers in both the Lalpha and Lbeta' phases. The effects of the two slowest motions on the density matrix are described by the stochastic Liouville equation (SLE) which is solved analytically using the method of van Faassen. These two slowest motions are assumed to be a rotational re-orientation about the long molecular axis and a uniform wobble of this axis within a conical volume with re-orientation rates characterised by correlation times tau(parallel) and tau(perpendicular) respectively. In the present work the Hamiltonian contains the intramolecular dipole-dipole interaction between the phosphorous nucleus and the four closest methylene protons of the choline headgroup, as well as the anisotropic chemical shielding interaction. Hence the contribution to relaxation from cross correlation between the dipole-dipole and anisotropic chemical shielding interactions is included. The reorientation of the headgroup is assumed to be a rotation sufficiently fast to lead to complete axially symmetric averaging of the Hamiltonian about the rotational axis (the P-O11 bond axis). Evaluation of the line shape in the present work involves only numerical integration and is therefore less computationally demanding than the large matrix inversions involved in the approaches of Campbell, Freed et al. The present theory also uses fewer parameters than that of Dufourc et al. but nevertheless results in good agreement with these authors' measurements on DMPC bilayers, using a fixed value of 10 for the ratio tau(perpendicular)/tau(parallel) in the case of the Lalpha phase. However, in contrast to Dufourc et al., we find that these correlation times are equal for the Lbeta' phase. Finally, we have simulated the decoupled powder line shapes obtained from the Lbeta' phase of DPPC by Campbell and Meirovitch. Again, we get good agreement providing tau(perpendicular)=tau(parallel).

Elastic deformation of membrane bilayers probed by deuterium NMR relaxation

In deuterium ((2)H) NMR spectroscopy of fluid lipid bilayers, the average structure is manifested in the segmental order parameters (S(CD)) of the flexible molecules. The corresponding spin-lattice relaxation rates (R(1Z) depend on both the amplitudes and the rates of the segmental fluctuations, and indicate the types of lipid motions. By combining (2)H NMR order parameter measurements with relaxation studies, we have obtained a more comprehensive picture of lipids in the liquid-crystalline (L(alpha)) state than formerly possible. Our data suggest that a lipid bilayer constitutes an ordered fluid, in which the phospholipids are grafted to the aqueous interface via their polar headgroups, whereas the fatty acyl chains are in effect liquid hydrocarbon. Studies of (2)H-labeled saturated lipids indicate their R(1Z) rates and S(CD) order parameters are correlated by a model-free, square-law functional dependence, signifying the presence of relatively slow bilayer fluctuations. A new composite membrane deformation model explains simultaneously the frequency (magnetic field) dependence and the angular anisotropy of the relaxation. The results imply the R(1Z) rates are due to a broad spectrum of 3-D collective bilayer excitations, together with effective axial rotations of the lipids. For the first time, NMR relaxation studies show that the viscoelastic properties of membrane lipids at megahertz frequencies are modulated by the lipid acyl length (bilayer thickness), polar headgroups (bilayer interfacial area), inclusion of a nonionic detergent (C(12)E(8)), and the presence of cholesterol, leading to a range of bilayer softness. Our findings imply the concept of elastic deformation is relevant on lengths approaching the bilayer thickness and less (the mesoscopic scale), and suggest that application of combined R(1Z) and S(CD) studies of phospholipids can be used as a simple membrane elastometer. Heuristic estimates of the bilayer bending rigidity kappa and the area elastic modulus K(a) enable comparison to other biophysical studies, involving macroscopic deformation of thin membrane lipid films. Finally, the bilayer softness may be correlated with the lipid diversity of biomembranes, for example, with regard to membrane curvature, repulsive interactions between bilayers, and lipid-protein interactions.

Low-frequency collective chain dynamics in a model biomembrane

Proton NMR was employed as a probe for the collective hydrocarbon chain dynamics in decylammonium chloride (C10H21NH3Cl), a model biomembrane undergoing an irreversible structural phase transition sequence. Our rotating frame spin-lattice relaxation measurements revealed a low-frequency critical collective chain dynamics in the kHz regime, which is associated with the interdigitated to noninterdigitated chain configurational phase transition.