2H NMR studies of isomeric omega 3 and omega 6 polyunsaturated phospholipid 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.

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.

Unperturbed hydrocarbon chains and liquid phase bilayer lipid chains: a computer simulation study

In this work, the properties of saturated and unsaturated fatty acid acyl chains 16:0, 18:0, 18:1(n-9)cis, 18:2(n-6)cis, 18:3(n-3)cis, 18:4(n-3)cis, 18:5(n-3)cis, 20:4(n-6)cis, 20:5(n-3)cis and 22:6(n-3)cis in a bilayer liquid crystalline state and similar hydrocarbon chains (with CH[Formula: see text] terminal groups instead of C=O groups) in the unperturbed state characterised by a lack of long-range interaction were investigated. The unperturbed hydrocarbon chains were modelled by Monte Carlo simulations at temperature [Formula: see text] K; sixteen fully hydrated homogeneous liquid crystalline phosphatidylcholine bilayers containing these chains were studied by molecular dynamics simulations at the same temperature. To eliminate effects of the simulation parameters, the molecular dynamics and Monte Carlo simulations were carried out using the same structural data and force field coefficients. From these computer simulations, the average distances between terminal carbon atoms of the chains (end-to-end distances) were calculated and compared. The trends in the end-to-end distances obtained for the unperturbed chains were found to be qualitatively similar to those obtained for the same lipid chains in the bilayers. So, for understanding of a number of processes in biological membranes (e.g., changes in fatty acid composition caused by environmental changes such as temperature and pressure), it is possible to use, at least as a first approximation, the relationships between the structure and properties for unperturbed or isolated hydrocarbon chains.

Effects of lipid interactions on model vesicle engulfment by alveolar macrophages

The engulfment function of macrophages relies on complex molecular interactions involving both lipids and proteins. In particular, the clearance of apoptotic bodies (efferocytosis) is enabled by externalization on the cell target of phosphatidylserine lipids, which activate receptors on macrophages, suggesting that (local) specific lipid-protein interactions are required at least for the initiation of efferocytosis. However, in addition to apoptotic cells, macrophages can engulf foreign bodies that vary substantially in size from a few nanometers to microns, suggesting that nonspecific interactions over a wide range of length scales could be relevant. Here, we use model lipid membranes (made of phosphatidylcholine, phosphatidylserine, and ceramide) and rat alveolar macrophages to show how lipid bilayer properties probed by small-angle x-ray scattering and solid-state (2)H NMR correlate with engulfment rates measured by flow cytometry. We find that engulfment of protein-free model lipid vesicles is promoted by the presence of phosphatidylserine lipids but inhibited by ceramide, in accord with a previous study of apoptotic cells. We conclude that the roles of phosphatidylserine and ceramide in phagocytosis is based, at least in part, on lipid-mediated modification of membrane physical properties, including interactions at large length scales as well as local lipid ordering and possible domain formation.

Oleic and docosahexaenoic acid differentially phase separate from lipid raft molecules: a comparative NMR, DSC, AFM, and detergent extraction study

We have previously suggested that the omega-3 polyunsaturated fatty acid, docosahexaenoic acid (DHA) may in part function by enhancing membrane lipid phase separation into lipid rafts. Here we further tested for differences in the molecular interactions of an oleic (OA) versus DHA-containing phospholipid with sphingomyelin (SM) and cholesterol (CHOL) utilizing (2)H NMR spectroscopy, differential scanning calorimetry, atomic force microscopy, and detergent extractions in model bilayer membranes. (2)H NMR and DSC (differential scanning calorimetry) established the phase behavior of the OA-containing 1-[(2)H(31)]palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (16:0-18:1PE-d(31))/SM (1:1) and the DHA-containing 1-[(2)H(31)]palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphoethanolamine (16:0-22:6PE-d(31))/SM (1:1) in the absence and presence of equimolar CHOL. CHOL was observed to affect the OA-containing phosphatidylethanolamine (PE) more than the DHA-containing PE, as exemplified by >2 x greater increase in order measured for the perdeuterated palmitic chain in 16:0-18:1PE-d(31)/SM (1:1) compared to 16:0-22:6PE-d(31)/SM (1:1) bilayers in the liquid crystalline phase. Atomic force microscopy (AFM) experiments showed less lateral phase separation between 16:0-18:1PE-rich and SM/CHOL-rich raft domains in 16:0-18:1PE/SM/CHOL (1:1:1) bilayers than was observed when 16:0-22:6PE replaced 16:0-18:1PE. Differences in the molecular interaction of 16:0-18:1PE and 16:0-22:6PE with SM/CHOL were also found using biochemical detergent extractions. In the presence of equimolar SM/CHOL, 16:0-18:1PE showed decreased solubilization in comparison to 16:0-22:6PE, indicating greater phase separation with the DHA-PE. Detergent experiments were also conducted with cardiomyocytes fed radiolabeled OA or DHA. Although both OA and DHA were found to be largely detergent solubilized, the amount of OA that was found to be associated with raft-rich detergent-resistant membranes exceeded DHA by almost a factor of 2. We conclude that the OA-PE phase separates from rafts far less than DHA-PE, which may have implications for cellular signaling.

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.

Docosahexaenoic acid: membrane properties of a unique fatty acid

Docosahexaenoic acid (DHA) with 22-carbons and 6 double bonds is the extreme example of an omega-3 polyunsaturated fatty acid (PUFA). DHA has strong medical implications since its dietary presence has been positively linked to the prevention of numerous human afflictions including cancer and heart disease. The PUFA, moreover, is essential to neurological function. It is remarkable that one simple molecule has been reported to affect so many seemingly unrelated biological processes. Although details of a molecular mode of action remain elusive, DHA must be acting at a fundamental level common to many tissues that is related to the high degree of conformational flexibility that the multiple double bonds have been identified to confer. One likely target for DHA action is at the cell membrane where the fatty acid is known to readily incorporate into membrane phospholipids. Once esterified into phospholipids DHA has been demonstrated to significantly alter many basic properties of membranes including acyl chain order and "fluidity", phase behavior, elastic compressibility, permeability, fusion, flip-flop and protein activity. It is concluded that DHA's interaction with other membrane lipids, particularly cholesterol, may play a prominent role in modulating the local structure and function of cell membranes.

Combined Monte Carlo and molecular dynamics simulation of hydrated lipid-cholesterol lipid bilayers at low cholesterol concentration

We have applied a hybrid equilibration and sampling procedure for the atomic level simulation of a hydrated lipid bilayer to systems consisting of dipalmitoyl phosphatidylcholine (DPPC) and cholesterol, and palmitoyl-oleyl phosphatidylcholine (POPC) at low (approximately 6%) cholesterol concentration. The procedure is applied to bilayers of 94 molecules of DPPC, 6 molecules of cholesterol, and 3205 water molecules, and to bilayers of 120 molecules of POPC, 8 molecules of cholesterol, and 4268 water molecules, at a temperature of 325 K. After equilibration, three separate 400-ps continuous molecular dynamics runs, separated by 10,000 configurational bias Monte Carlo steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of DPPC and cholesterol, and of pure DPPC, and pure POPC. Certain properties of the bilayers are indistinguishable from cholesterol-free bilayers, including lateral diffusion and electron density. Other properties, most notably the order parameter profile, show the effect of cholesterol even at low concentrations.

Rapid deconvolution of NMR powder spectra by weighted fast Fourier transformation

The resolution enhancement conferred by numerical deconvolution of powder pattern spectra to spectra characteristic of a single alignment greatly simplifies solid state NMR spectral analysis. This is especially beneficial when the spectrum is a superposition of signals from multiple environments or sites of labelling. We have developed an innovative method to deconvolute (depake) spectra governed by axially symmetric second rank tensor interactions which possess a P2(cos theta) dependence upon orientation, where theta is the angle between the symmetry axis and the external magnetic field. Our approach differs substantially from previously published procedures which are iterative or require matrix inversion and, hence, are slow. The new method, instead, utilizes weighting functions in time and frequency domains to facilitate a rapidly executed solution based upon fast Fourier transformation (FFT). Its efficacy is demonstrated with 2H and 31P NMR data for model membranes.

Structure and dynamic properties of diunsaturated 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine lipid bilayer from molecular dynamics simulation

Unsaturated fatty acid chains are known to be an essential structural part of biomembranes, but only monounsaturated chains have been included in the molecular dynamics (MD) simulations of membrane systems. Here we present a 1-ns MD simulation for a diunsaturated 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC; 16:0/18:2[delta9,12]) bilayer. The structural behavior of the phosphatidylcholine headgroup, the glycerol backbone, and the hydrating water were assessed and found to be consistent with the existing information about similar systems from both experimental and computational studies. Further analysis was focused on the structure of the double bond region and the effects of the diunsaturation on the bilayer interior. The behavior of the diunsaturated sn-2 chains is affected by the tilted beginning of the chain and the four main conformations of the double bond region. The double bonds of the sn-2 chains also influenced the characteristics of the saturated chains in the sn-1 position. Furthermore, extreme conformations of the sn-2 chains existed that are likely to be related to the functional role of the double bonds. The results here point out the importance of polyunsaturation for the biological interpretations deduced from the membrane MD simulations.

Influence of the positions of cis double bonds in the sn-2-acyl chain of phosphatidylethanolamine on the bilayer's melting behavior

In an attempt to examine the effects of different numbers and positions of cis double bonds in the sn-2-acyl chain of phosphatidylethanolamine (PE) on the bilayer's melting behavior, 21 molecular species of PE were first semisynthesized, and their Tm and DeltaH values were subsequently determined by high resolution differential scanning calorimetry. In the plot of Tm versus the number of the cis double bond, some characteristic profiles were observed for the various series of PEs. For instance, if the cis double bond was first introduced into the sn-2-acyl chain of C(20):C(20)PE at the Delta5-position, the Tm was observed to reduce drastically. Subsequent stepwise additions of up to five cis double bonds at the methylene-interrupted positions toward the methyl end resulted in a progressive yet smaller decrease in Tm. If, on the other hand, the cis double bonds were introduced sequentially at the Delta11-, Delta11,14-, and Delta11,14,17-positions along the sn-2-acyl chain of C(20):C(20)PE, the Tm profile in the Tm versus the number of the cis double bond showed a down-and-up trend. Most interestingly, for positional isomers of C(20):C(20:3Delta5,8,11)PE, C(20):C(20:3Delta8,11,14)PE, and C(20):C(20:3Delta11,14,17)PE, an inverted bell-shaped Tm profile was detected in the plot of Tm against the position of the omega-carbon for these isomers. Similar Tm profiles were also observed for C(18):C(20)PE, C(20):C(18)PE, and their unsaturated derivatives. This work thus demonstrated that both the positions and the numbers of cis double bonds in the sn-2 acyl chain could exert noticeable influence on the gel-to-liquid crystalline phase transition behavior of the lipid bilayer. Finally, a molecular model was presented, with which the behavior of the gel-to-liquid crystalline phase transition observed for lipid bilayers composed of various sn-1-saturated/sn-2-unsaturated lipids can be rationalized.

Effect of unsaturation on the chain order of phosphatidylcholines in a dioleoylphosphatidylethanolamine matrix

The properties of phosphatidylcholines (PCs) having a perdeuterated stearic acid, 18:0d35, in the sn-1 position and the fatty acid 18:0, 18:1 omega 9, 18:2 omega 6, 18:3 omega 3, 20:4 omega 6, 20:5 omega 3, or 22:6 omega 3 at the sn-2 position were investigated in a matrix of dioleoylphosphatidylethanolamine (DOPE) by 2H and 31P NMR spectroscopy. At a mole ratio of DOPE/PC = 5:1, the lipids form liquid crystalline lamellar phases below 40 degrees C and coexisting lamellar, inverse hexagonal (Hll), and cubic phases at higher temperatures. The sn-1 chain of the PCs in a DOPE matrix is appreciably more ordered than in pure PCs, corresponding to an increase in the hydrophobic bilayer thickness of approximately 1 A. Distearoylphosphatidylcholine in the DOPE matrix has a higher sn-1 chain order than the unsaturated PCs. We observed distinct differences in the lipid order of upper and lower sections of the hydrocarbon chains caused by changes of temperature, unsaturation, headgroups, and ethanol. Unsaturation lowers chain order, mostly in the lower third of the hydrocarbon chains. By contrast, the increase in chain order caused by the DOPE matrix and the decrease in order with increasing temperature have a constant magnitude for the upper two-thirds of the chain and are smaller for the lower third. Addition of 2 M ethanol reduced order parameters, in effect reversing the increase in chain order caused by the DOPE matrix.

Incorporation of surface tension into molecular dynamics simulation of an interface: a fluid phase lipid bilayer membrane

In this paper we report on the molecular dynamics simulation of a fluid phase hydrated dimyristoylphosphatidylcholine bilayer. The initial configuration of the lipid was the x-ray crystal structure. A distinctive feature of this simulation is that, upon heating the system, the fluid phase emerged from parameters, initial conditions, and boundary conditions determined independently of the collective properties of the fluid phase. The initial conditions did not include chain disorder characteristic of the fluid phase. The partial charges on the lipids were determined by ab initio self-consistent field calculations and required no adjustment to produce a fluid phase. The boundary conditions were constant pressure and temperature. Thus the membrane was not explicitly required to assume an area/phospholipid molecule thought to be characteristic of the fluid phase, as is the case in constant volume simulations. Normal to the membrane plane, the pressure was 1 atmosphere, corresponding to the normal laboratory situation. Parallel to the membrane plane a negative pressure of -100 atmospheres was applied, derived from the measured surface tension of a monolayer at an air-water interface. The measured features of the computed membrane are generally in close agreement with experiment. Our results confirm the concept that, for appropriately matched temperature and surface pressure, a monolayer is a close approximation to one-half of a bilayer. Our results suggest that the surface area per phospholipid molecule for fluid phosphatidylcholine bilayer membranes is smaller than has generally been assumed in computational studies at constant volume. Our results confirm that the basis of the measured dipole potential is primarily water orientations and also suggest the presence of potential barriers for the movement of positive charges across the water-headgroup interfacial region of the phospholipid.

2H nuclear magnetic resonance order parameter profiles suggest a change of molecular shape for phosphatidylcholines containing a polyunsaturated acyl chain

Solid-state 2H nuclear magnetic resonance spectroscopy was used to determine the orientational order parameter profiles for a series of phosphatidylcholines with perdeuterated stearic acid, 18:0d35, in position sn-1 and 18:1 omega 9, 18:2 omega 6, 18:3 omega 3, 20:4 omega 6, 20:5 omega 3, or 22:6 omega 3 in position sn-2. The main phase transition temperatures were derived from a first moment analysis, and order parameter profiles of sn-1 chains were calculated from dePaked nuclear magnetic resonance powder patterns. Comparison of the profiles at 37 degrees C showed that unsaturation causes an inhomogenous disordering along the sn-1 chain. Increasing sn-2 chain unsaturation from one to six double bonds resulted in a 1.6-kHz decrease in quadrupolar splittings of the sn-1 chain in the upper half of the chain (or plateau region) and maximum splitting difference of 4.4 kHz at methylene carbon 14. The change in chain order corresponds to a decrease in the 18:0 chain length of 0.4 +/- 0.2 A with 18:2 omega 6 versus 18:1 omega 9 in position sn-2. Fatty acids containing three or more double bonds in sn-2 showed a decrease in sn-1 chain length of 0.7 +/- 0.2 A compared with 18:1 omega 9. The chain length of all lipids decreased with increasing temperature. Highly unsaturated phosphatidylcholines (three or more double bonds in sn-2) had shorter sn-1 chains, but the chain length was somewhat less sensitive to temperature. The profiles reveal that the sn-1 chain exhibits a selective increase in motional freedom in a region located toward the bottom half of the chain as sn-2 unsaturation is increased. This corresponds to an area increase around carbon atom number 14 that is three to four times greater than the increase for the top part of the chain. A similar asymmetric decrease in order, largest toward the methyl end of the chain, was observed when 1 -palmitoyl-2-oleoylphosphatidylethanolamine goes from a lamellar to an inverse hexagonal (H,,) phase. This is consistent with a change to a more wedge-shaped space available for the acyl chain.

Cholesterol condensation of alpha-linolenic and gamma-linolenic acid-containing phosphatidylcholine monolayers and bilayers

Cholesterol is demonstrated to condense phosphatidylcholine (PC) monolayers and bilayers containing stearic acid in the sn-1 position and alpha-linolenic acid in the sn-2 position (18:0, alpha-18:3 PC) but has no effect when gamma-linolenic acid occupies the sn-2 position (18:0,gamma-18:3 PC). Cholesterol-induced condensation is measured by area/molecule determinations made on monolayers using a Langmuir trough, while condensation in bilayers is followed by the fluorescent dyes merocyanine (MC540) and dansyllysine. Permeability to erythritol is also demonstrated to be diminished by cholesterol for the condensable 18:0,alpha-18:3 PC bilayer membranes but not the 18:0,gamma-18:3 PC membranes. alpha- and gamma-linolenic acid are isomers containing 18 carbons and three unsaturations. Both fatty acids have unsaturations at positions 9 and 12 and differ only in the location of the third unsaturation, at either position 6 for gamma-linolenic acid (an omega-6 fatty acid) and at position 15 for alpha-linolenic acid (an omega-3 fatty acid). Here lipid-cholesterol interaction is used to distinguish the effect of position of unsaturation on membrane structure.