Structural properties of docosahexaenoyl phospholipid bilayers investigated by solid-state 2H NMR spectroscopy.

Polyunsaturated lipids in cellular membranes are known to play key roles in such diverse biological processes as vision, neuronal signaling, and apoptosis. One hypothesis is that polyunsaturated lipids are involved in second messenger functions in biological signaling. Another current hypothesis affirms that the functional role of polyunsaturated lipids relies on their ability to modulate physical properties of the lipid bilayer. The present research has employed solid-state 2H NMR spectroscopy to acquire knowledge of the molecular organization and material properties of polyunsaturated lipid bilayers. We report measurements for a homologous series of mixed-chain phosphatidylcholines containing a perdeuterated, saturated acyl chain (n:0) at the sn-1 position, adjacent to docosahexaenoic acid (DHA, 22:6omega3) at the sn-2 position. Measurements have been performed on fluid (L(alpha))-state multilamellar dispersions as a function of temperature for saturated acyl chain lengths of n = 12, 14, 16, and 18 carbons. The saturated sn-1 chains are therefore used as an intrinsic probe with site-specific resolution of the polyunsaturated bilayer structure. The 2H NMR order parameters as a function of acyl position (order profiles) have been analyzed using a mean-torque potential model for the chain segments, and the results are discussed in comparison with the homologous series of disaturated lipid bilayers. At a given absolute temperature, as the sn-1 acyl length adjacent to the sn-2 DHA chain is greater, the order of the initial chain segments increases, whereas that of the end segments decreases, in marked contrast with the corresponding disaturated series. For the latter, the order of the end segments is practically constant with acyl length, thus revealing a universal chain packing profile. We find that the DHA-containing series, while more complex, is still characterized by a universal chain packing profile, which is shifted relative to the homologous saturated series. Moreover, we show how introduction of DHA chains translates the order profile along the saturated chains, making more disordered states accessible within the bilayer central region. As a result, the area per lipid headgroup is increased as compared to disaturated bilayers. The systematic analysis of the 2H NMR data provides a basis for studies of lipid interactions with integral membrane proteins, for instance in relation to characteristic biological functions of highly unsaturated lipid membranes.

Area per lipid and acyl length distributions in fluid phosphatidylcholines determined by (2)H NMR spectroscopy.

Deuterium ((2)H) NMR spectroscopy provides detailed information regarding the structural fluctuations of lipid bilayers, including both the equilibrium properties and dynamics. Experimental (2)H NMR measurements for the homologous series of 1, 2-diacyl-sn-glycero-3-phosphocholines with perdeuterated saturated chains (from C12:0 to C18:0) have been performed on randomly oriented, fully hydrated multilamellar samples. For each lipid, the C-D bond order parameters have been calculated from de-Paked (2)H NMR spectra as a function of temperature. The experimental order parameters were analyzed using a mean-torque potential model for the acyl chain segment distributions, and comparison was made with the conventional diamond lattice approach. Statistical mechanical principles were used to relate the measured order parameters to the lipid bilayer structural parameters: the hydrocarbon thickness and the mean interfacial area per lipid. At fixed temperature, the area decreases with increasing acyl length, indicating increased van der Waals attraction for longer lipid chains. However, the main effect of increasing the acyl chain length is on the hydrocarbon thickness rather than on the area per lipid. Expansion coefficients of the structural parameters are reported and interpreted using an empirical free energy function that describes the force balance in fluid bilayers. At the same absolute temperature, the phosphatidylcholine (PC) series exhibits a universal chain packing profile that differs from that of phosphatidylethanolamines (PE). Hence, the lateral packing of phospholipids is more sensitive to the headgroup methylation than to the acyl chain length. A fit to the area per lipid for the PC series using the empirical free energy function shows that the PE area represents a limiting value for the packing of fluid acyl chains.

Modulation of glycophorin A transmembrane helix interactions by lipid bilayers: molecular dynamics calculations.

Starting from the glycophorin A dimer structure determined by NMR, we performed simulations of both dimer and monomer forms in explicit lipid bilayers with constant normal pressure, lateral area, and temperature using the CHARMM potential. Analysis of the trajectories in four different lipids reveals how lipid chain length and saturation modulate the structural and energetic properties of transmembrane helices. Helix tilt, helix-helix crossing angle, and helix accessible volume depend on lipid type in a manner consistent with hydrophobic matching concepts: the most relevant lipid property appears to be the bilayer thickness. Although the net helix-helix interaction enthalpy is strongly attractive, analysis of residue-residue interactions reveals significant unfavorable electrostatic repulsion between interfacial glycine residues previously shown to be critical for dimerization. Peptide volume is nearly conserved upon dimerization in all lipid types, indicating that the monomeric helices pack equally well with lipid as dimer helices do with one another. Enthalpy calculations indicate that the helix-environment interaction energy is lower in the dimer than in the monomer form, when solvated by unsaturated lipids. In all lipid environments there is a marked preference for lipids to interact with peptide predominantly through one rather than both acyl chains. Although our trajectories are not long enough to allow a full thermodynamic treatment, these results demonstrate that molecular dynamics simulations are a powerful method for investigating the protein-protein, protein-lipid, and lipid-lipid interactions that determine the structure, stability and dynamics of transmembrane alpha-helices in membranes.

Analysis of simulated NMR order parameters for lipid bilayer structure determination.

The conventional formula for relating CD2 average order parameters to average methylenic travel is flawed when compared to molecular dynamics simulations of dipalmitoylphosphatidylcholine. Inspired by the simulated probability distribution functions, a new formula is derived that satisfactorily relates these quantities. This formula is used to obtain the average chain length , and the result agrees with the direct simulation result for . The simulation also yields a hydrocarbon thickness 2. The result = is consistent with a model of chain packing with both early chain termination and partial interdigitation of chains from opposing monolayers. The actual simulated area per lipid is easily obtained from the order parameters. However, when this method is applied to NMR order parameter data from dimyristoylphosphatidylcholine, the resulting is 10% larger than the currently accepted value.

Fluid phase structure of EPC and DMPC bilayers.

X-ray diffraction data taken at high instrumental resolution were obtained for EPC and DMPC under various osmotic pressures, primarily at T = 30 degrees C. The headgroup thickness DHH was obtained from relative electron density profiles. By using volumetric results and by comparing to gel phase DPPC we obtain areas AEPCF = 69.4 +/- 1.1 A2 and ADMPCF = 59.7 +/- 0.2 A2. The analysis also gives estimates for the areal compressibility KA. The AF results lead to other structural results regarding membrane thickness and associated waters. Using the recently determined absolute electrons density profile of DPPC, the AF results also lead to absolute electron density profiles and absolute continuous transforms magnitude of F(q) for EPC and DMPC. Limited measurements of temperature dependence show directly that fluctuations increase with increasing temperature and that a small decrease in bending modulus KC accounts for the increased water spacing reported by Simon et al. (1995) Biophys. J. 69, 1473-1483.

Structure and interactions of fully hydrated dioleoylphosphatidylcholine bilayers.

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

DMSO produces a new subgel phase in DPPC: DSC and X-ray diffraction study.

Equilibrium phases and the kinetics of subgel phase transformation of dipalmitoylphosphatidylcholine (DPPC) hydrated with mixtures of dimethylsulfoxide (DMSO)/water have been studied using differential scanning calorimetry (DSC) and X-ray diffraction (XRD). The rate of gel-to-subgel transformation is decreased with a small increase in X, the DMSO/water mole fraction, but then speeds up and becomes faster than in pure water by X = 0.16. The DSC scans show multiple subgel peaks, some of which can be attributed to impacted domain growth. For X greater than 0.10, XRD shows that there is a new, stable subgel phase, S, which also accounts for some of the multiplicity of DSC peaks. Our electron density profiles show that the thickness of the bilayer in the S phase is greater than in the usual C subgel phase. We suggest that the S subgel phase is characterized by different headgroup ordering and smaller chain tilt angle than in the C subgel phase. Electron density profiles show that increasing X decreases the water space between bilayers in all phases, subgel, gel and fluid (L alpha). For X = 0.20, a different gel phase is also observed that may be due to subtle changes in the orientation of chain tilt first observed in partially dehydrated DMPC. The dehydrating effect of DMSO explains the results of a previous study, confirmed in this study, that increasing the concentration of DMSO raises the main transition temperature and eliminates the ripple phase.

Effect of substrate roughness on D spacing supports theoretical resolution of vapor pressure paradox.

The lamellar D spacing has been measured for oriented stacks of lecithin bilayers prepared on a variety of solid substrates and hydrated from the vapor. We find that, when the bilayers are in the L(alpha) phase near 100% relative humidity, the D spacing is consistently larger when the substrate is rougher than when it is smooth. The differences become smaller as the relative humidity is decreased to 80% and negligible differences are seen in the L(beta') phase. Our interpretation is that rough substrates frustrate the bilayer stack energetically, thereby increasing the fluctuations, the fluctuational repulsive forces, and the water spacing compared with stacks on smooth surfaces. This interpretation is consistent with and provides experimental support for a recently proposed theoretical resolution of the vapor pressure paradox.

Determination of component volumes of lipid bilayers from simulations.

An efficient method for extracting volumetric data from simulations is developed. The method is illustrated using a recent atomic-level molecular dynamics simulation of L alpha phase 1,2-dipalmitoyl-sn-glycero-3-phosphocholine bilayer. Results from this simulation are obtained for the volumes of water (VW), lipid (V1), chain methylenes (V2), chain terminal methyls (V3), and lipid headgroups (VH), including separate volumes for carboxyl (Vcoo), glyceryl (Vgl), phosphoryl (VPO4), and choline (Vchol) groups. The method assumes only that each group has the same average volume regardless of its location in the bilayer, and this assumption is then tested with the current simulation. The volumes obtained agree well with the values VW and VL that have been obtained directly from experiment, as well as with the volumes VH, V2, and V3 that require certain assumptions in addition to the experimental data. This method should help to support and refine some assumptions that are necessary when interpreting experimental data.

X-ray structure determination of fully hydrated L alpha phase dipalmitoylphosphatidylcholine bilayers.

Bilayer form factors obtained from x-ray scattering data taken with high instrumental resolution are reported for multilamellar vesicles of L alpha phase lipid bilayers of dipalmitoylphosphatidylcholine at 50 degrees C under varying osmotic pressure. Artifacts in the magnitudes of the form factors due to liquid crystalline fluctuations have been eliminated by using modified Caillé theory. The Caillé fluctuation parameter eta 1 increases systematically with increasing lamellar D spacing and this explains why some higher order peaks are unobservable for the larger D spacings. The corrected form factors fall on one smooth continuous transform F(q); this shows that the bilayer does not change shape as D decreases from 67.2 A (fully hydrated) to 60.9 A. The distance between headgroup peaks is obtained from Fourier reconstruction of samples with four orders of diffraction and from electron density models that use 38 independent form factors. By combining these results with previous gel phase results, area AF per lipid molecule and other structural quantities are obtained for the fluid L alpha phase. Comparison with results that we derived from previous neutron diffraction data is excellent, and we conclude from diffraction studies that AF = 62.9 +/- 1.3 A2, which is in excellent agreement with a previous estimate from NMR data.