Evidence for a mechanism by which omega-3 polyunsaturated lipids may affect membrane protein function

We have calculated the lateral pressure profile from well-converged, experimentally validated, molecular dynamics simulations of hydrated lipid bilayer membranes containing highly polyunsaturated fatty acids. The three simulations, each 30 ns in length, contain omega-3 fatty acids, omega-6 fatty acids, and a mixture of omega-3 fatty acids and cholesterol and were continued from previously published simulations that demonstrated excellent agreement with a wide variety of experimental measurements. We find that the distribution of lateral stress within the hydrophobic core of the membrane is sensitively dependent on the degree of chain unsaturation and on the presence of cholesterol. Replacing omega-3 fatty acids with omega-6 chains, or incorporating cholesterol into the membrane, shifts the repulsive lateral chain pressure away from the lipid/water interface toward the bilayer interior. This may support a previously proposed mechanism by which lipid composition may affect conformational equilibrium for integral membrane proteins.     

Structure and Dynamics of Cholesterol-Containing Polyunsaturated Lipid Membranes Studied by Neutron Diffraction and NMR

A direct and quantitative analysis of the internal structure and dynamics of a polyunsaturated lipid bilayer composed of 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (18:0-22:6n3-PC) containing 29 mol% cholesterol was carried out by neutron diffraction, (2)H-NMR and (13)C-MAS NMR. Scattering length distribution functions of cholesterol segments as well as of the sn-1 and sn-2 hydrocarbon chains of 18:0-22:6n3-PC were obtained by conducting experiments with specifically deuterated cholesterol and lipids. Cholesterol orients parallel to the phospholipids, with the A-ring near the lipid glycerol and the terminal methyl groups 3 ? away from the bilayer center. Previously, we reported that the density of polyunsaturated docosahexaenoic acid (DHA, 22:6n3) chains was higher near the lipid-water interface. Addition of cholesterol partially redistributes DHA density from near the lipid-water interface to the center of the hydrocarbon region. Cholesterol raises chain-order parameters of both stearic acid and DHA chains. The fractional order increase for stearic acid methylene carbons C(8)-C(18) is larger, reflecting the redistribution of DHA chain density toward the bilayer center. The correlation times of DHA chain isomerization are short and mostly unperturbed by the presence of cholesterol. The uneven distribution of saturated and polyunsaturated chain densities and the cholesterol-induced balancing of chain distributions may have important implications for the function and integrity of membrane receptors, such as rhodopsin.     

Properties of unsaturated phospholipid bilayers: Effect of cholesterol

Properties of hydrated unsaturated phosphatidylcholine (PC) lipid bilayers containing 40 mol % cholesterol and of pure PC bilayers have been studied. Various methods were applied, including molecular dynamics simulations, self-consistent field calculations, and the pulsed field gradient nuclear magnetic resonance technique. Lipid bilayers were composed of 18:0/18:1(n-9)cis PC, 18:0/18:2(n-6)cis PC, 18:0/18:3(n-3)cis PC, 18:0/20:4(n-6)cis PC, and 18:0/22:6(n-3)cis PC molecules. Lateral self-diffusion coefficients of the lipids in all these bilayers, mass density distributions of atoms and atom groups with respect to the bilayer normal, the C-H and C-C bond order parameter profiles of each phospholipid hydrocarbon chain with respect to the bilayer normal were calculated. It was shown that the lateral self-diffusion coefficient of PC molecules of the lipid bilayer containing 40 mol % cholesterol is smaller than that for a corresponding pure PC bilayer; the diffusion coefficients increase with increasing the degree of unsaturation of one of the PC chains in bilayers of both types (i.e., in pure bilayers or in bilayers with cholesterol). The presence of cholesterol in a bilayer promoted the extension of saturated and polyunsaturated lipid chains. The condensing effect of cholesterol on the order parameters was more pronounced for the double C=C bonds of polyunsaturated chains than for single C-C bonds of saturated chains.     

Constant-pressure molecular dynamics investigation of cholesterol effects in a dipalmitoylphosphatidylcholine bilayer.

We report a 1.4-ns constant-pressure molecular dynamics simulation of cholesterol at 12.5 mol% in a dipalmitoylphosphatidylcholine (DPPC) bilayer at 50 degrees C and compare the results to our previous simulation of a pure DPPC bilayer. The interlamellar spacing was increased by 2.5 A in the cholesterol-containing bilayer, consistent with x-ray diffraction results, whereas the bilayer thickness was increased by only 1 A. The bilayer/water interface was more abrupt because the lipid headgroups lie flatter to fill spaces left by the cholesterol molecules. This leads to less compensation by the lipid headgroups of the oriented water contribution to the membrane dipole potential and could explain the experimentally observed increase in the magnitude of the dipole potential by cholesterol. Our calculations suggested that 12.5 mol% cholesterol does not significantly affect the conformations and packing of the hydrocarbon chains and produces only a slight reduction in the empty free volume. However, cholesterol has a significant influence on the subnanosecond time scale lipid dynamics: the diffusion constant for the center-of-mass "rattling" motion was reduced by a factor of 3, and the reorientational motion of the methylene groups was slowed along the entire length of the hydrocarbon chains.     

Structural properties of a highly polyunsaturated lipid bilayer from molecular dynamics simulations.

The structure of a fully hydrated mixed (saturated/polyunsaturated) chain lipid bilayer in the biologically relevant liquid crystalline phase has been examined by performing a molecular dynamics study. The model membrane, a 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (SDPC, 18:0/22:6 PC) lipid bilayer, was investigated at constant (room) temperature and (ambient) pressure, and the results obtained in the nanosecond time scale reproduced quite well the available experimental data. Polyunsaturated fatty acids are found in high concentrations in neuronal and retinal tissues and are essential for the development of human brain function. The docosahexaenoic fatty acid, in particular, is fundamental for the proper function of the visual receptor rhodopsin. The lipid bilayer order has been investigated through the orientational order parameters. The water-lipid interface has been explored thoroughly in terms of its dimensions and the organization of the different components. Several types of interactions occurring in the system have been analyzed, specifically, the water-hydrocarbon chain, lipid-lipid and lipid-water interactions. The distribution of dihedral angles along the chains and the molecular conformations of the polyunsaturated chain of the lipids have also been studied. Special attention has been focused on the microscopic (molecular) origin of the effects of polyunsaturations on the different physical properties of membranes.     

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

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

Biophysical and biochemical mechanisms by which dietary N-3 polyunsaturated fatty acids from fish oil disrupt membrane lipid rafts

N-3 polyunsaturated fatty acids (PUFAs) from fish oil exert their functional effects by targeting multiple mechanisms. One mechanism to emerge in the past decade is the ability of n-3 PUFA acyl chains to perturb the molecular organization of plasma membrane sphingolipid/cholesterol-enriched lipid raft domains. These domains are nanometer-scale assemblies that coalesce to compartmentalize select proteins for optimal function. Here we review recent evidence on how n-3 PUFAs modify lipid rafts from biophysical and biochemical experiments from several different model systems. A central theme emerges from these studies. N-3 PUFA acyl chains display tremendous conformational flexibility and a low affinity for cholesterol and saturated acyl chains. This unique flexibility of n-3 PUFA acyl chains impacts the organization of inner and outer leaflet lipid rafts by disrupting acyl chain packing and molecular order within rafts. Ultimately, the disruption in raft organization has consequences for protein clustering and thereby signaling. Overall, elucidating the complex mechanisms by which n-3 PUFA acyl chains reorganize membrane architecture will enhance the translation of these fatty acids into the clinic for treating several diseases.     

The structure of polyunsaturated lipid bilayers important for rhodopsin function: a neutron diffraction study

The structure of oriented 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine bilayers with perdeuterated stearoyl- or docosahexaenoyl hydrocarbon chains was investigated by neutron diffraction. Experiments were conducted at two different relative humidities, 66 and 86%. At both humidities we observed that the polyunsaturated docosahexaenoyl chain has a preference to reside near the lipid water interface. That leaves voids in the bilayer center that are occupied by saturated stearoyl chain segments. This uneven distribution of saturated- and polyunsaturated chain densities is likely to result in membrane elastic stress that modulates function of integral receptor proteins like rhodopsin.     

Interplay of unsaturated phospholipids and cholesterol in membranes: effect of the double-bond position

The structural and dynamical properties of lipid membranes rich in phospholipids and cholesterol are known to be strongly affected by the unsaturation of lipid acyl chains. We show that not only unsaturation but also the position of a double bond has a pronounced effect on membrane properties. We consider how cholesterol interacts with phosphatidylcholines comprising two 18-carbon long monounsaturated acyl chains, where the position of the double bond is varied systematically along the acyl chains. Atomistic molecular dynamics simulations indicate that when the double bond is not in contact with the cholesterol ring, and especially with the C18 group on its rough β-side, the membrane properties are closest to those of the saturated bilayer. However, any interaction between the double bond and the ring promotes membrane disorder and fluidity. Maximal disorder is found when the double bond is located in the middle of a lipid acyl chain, the case most commonly found in monounsaturated acyl chains of phospholipids. The results suggest a cholesterol-mediated lipid selection mechanism in eukaryotic cell membranes. With saturated lipids, cholesterol promotes the formation of highly ordered raft-like membrane domains, whereas domains rich in unsaturated lipids with a double bond in the middle remain highly fluid despite the presence of cholesterol.     

Comparative molecular dynamics study of lipid membranes containing cholesterol and ergosterol

Sterol molecules are essential for maintaining the proper structure and function of eukaryotic cell membranes. The influence of cholesterol (the principal sterol of higher animals) on the lipid bilayer properties was extensively studied by both experimental and simulation methods. In contrast, the effect of ergosterol (the principal fungal sterol) on the membrane structure and dynamics is much less recognized. This work presents the results of comparative molecular dynamics simulation of the hydrated dimyristoylphosphatidylcholine bilayer containing approximately 25 mol % of cholesterol or ergosterol. A detailed analysis of the molecular properties (e.g., bilayer thickness, lipid order, diffusion, intermolecular interactions, etc.) of both sterol-induced liquid-ordered membrane phases is presented. Presence of sterols in the membrane significantly changes its property, especially fluidity and molecular packing. Moreover, in accordance with the experiments, our calculations show that, compared to cholesterol, ergosterol has higher ordering effect on the phospholipid acyl chains. This different influence on the properties of the lipid bilayer stems from differences in conformational freedom of sterol side chains. Additionally, obtained models of lipid membranes containing human and fungal sterols, constituting the result of our work, can be also utilized in other chemotherapeutic studies on interaction of selected ligands (e.g., antifungal compounds) with membranes.     

Cholesterol hydroxyl group is found to reside in the center of a polyunsaturated lipid membrane

Cholesterol and saturated lipid species preferentially partition into liquid ordered microdomains, such as lipid rafts, away from unsaturated lipid species for which the sterol has less affinity in the surrounding liquid-disordered membrane. To observe how cholesterol interacts with unsaturated phospholipids, we have determined, from one-dimensional neutron scattering length density profiles, the depth of cholesterol in phosphatidylcholine (PC) bilayers with varying amounts of acyl chain unsaturation. Through the use of [2,2,3,4,4,6-(2)H(6)]-labeled cholesterol, we show that in 1-palmitoyl-2-oleoylphosphatidylcholine (16:0-18:1 PC), 1,2-dioleoylphosphatidylcholine (18:1-18:1 PC), and 1-stearoyl-2-arachidonylphosphatidylcholine (18:0-20:4 PC) bilayers the center of mass of the deuterated sites is approximately 16 A from the bilayer center. This location places the hydroxyl group of the sterol moiety at the hydrophobic/hydrophilic bilayer interface, which is the generally accepted position. In dramatic contrast, for 20:4-20:4 PC membranes the hydroxyl group is found, unequivocally, sequestered in the bilayer center. We attribute the change in location to the high disorder of polyunsaturated fatty acids (PUFA) that is incompatible with close proximity to the steroid moiety in its usual "upright" orientation.     

The Structure of Polyunsaturated Lipid Bilayers Important for Rhodopsin Function: A Neutron Diffraction Study

The structure of oriented 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine bilayers with perdeuterated stearoyl- or docosahexaenoyl hydrocarbon chains was investigated by neutron diffraction. Experiments were conducted at two different relative humidities, 66 and 86%. At both humidities we observed that the polyunsaturated docosahexaenoyl chain has a preference to reside near the lipid water interface. That leaves voids in the bilayer center that are occupied by saturated stearoyl chain segments. This uneven distribution of saturated- and polyunsaturated chain densities is likely to result in membrane elastic stress that modulates function of integral receptor proteins like rhodopsin.     

Movement of cholesterol between vesicles prepared with different phospholipids or sizes

The rates of exchange of [4-14C]cholesterol between lipid vesicles prepared with different phospholipids and with different sizes have been measured. The first-order rate constants were higher using vesicles prepared from phosphatidylcholines with highly branched or polyunsaturated fatty acyl chains than with saturated diacyl or di-O-alkyl chains. The rate measurements indicate that the affinity of cholesterol for phospholipid does not vary significantly on change of the type of linkage (ether or ester) in phosphatidylcholine (PC) or of the positions of the fatty acyl chains in 1,2-diacyl-PC bearing one saturated and one unsaturated chain; furthermore, egg phosphatidylglycerol and egg phosphatidylethanolamine appear to have comparable affinities for cholesterol. However, the molecular packing in the bilayer and nearest-neighbor interactions involving cholesterol appear tightened more by N-palmitoylsphingomyelin than by dipalmitoyl-PC; on incorporation of 44 mol % of these phospholipids (which have the same fatty acyl chain composition) into either small or large unilamellar vesicles prepared with egg phosphatidylglycerol, the exchange rates were strikingly slower when the donor species contained sphingomyelin compared with PC. The rate of cholesterol exchange was 100% faster with small unilamellar vesicles than with large unilamellar vesicles as donors, suggesting that the looser packing in the highly curved small vesicles facilitates cholesterol desorption. The cholesterol exchange rate did not vary with the size of the acceptor vesicles, which indicates that desorption is the rate-limiting step in the exchange process in the presence of excess acceptors.     

The effect of variations in growth temperature, fatty acid composition and cholesterol content on the lipid polar head-group composition of Acholeplasma laidlawii B membranes

We have systematically investigated the effect of variations in growth temperature, fatty acid composition and cholesterol content on the membrane lipid polar headgroup composition of Acholeplasma laidlawii B. Two important lipid compositional parameters have been determined from such an analysis. The first parameter studied was the ratio of the two major neutral glycolipids of this organism, monoglucosyldiacylglycerol (MGDG) and diglucosyldiacylglycerol (DGDG). As the former lipid prefers to exist in a reversed hexagonal phase at higher temperatures, with unsaturated fatty acyl chains or in the presence of cholesterol, the ratio of these two lipids reflects the phase state preference of the total A. laidlawii membrane lipids. Although we find that the MGDG/DGDG ratio is reduced in response to an increase in fatty acid unsaturation, increases in growth temperature or cholesterol content reduce this ratio only in cells enriched in a saturated but not an unsaturated fatty acid. The second parameter studied was the ratio of these neutral glycolipids to the only phosphatide in the A. laidlawii membrane, phosphatidylglycerol (PG); this parameter reflects the relative balance of uncharged and charged lipids in the membrane of this organism. We find that the MGDG + DGDG/PG ratio is lowest in cells enriched in the saturated fatty acid even though these cells already have the highest lipid bilayer surface charge density. Moreover, this ratio is not consistently related to growth temperature or changes in cholesterol levels, as expected. We therefore conclude that A. laidlawii strain B, apparently unlike strain A, does not possess coherent regulatory mechanisms for maintaining either the phase preference or the surface charge density of its membrane lipid constant in response to variations in growth temperature, fatty acid composition or cholesterol content.     

The Molecular Basis of Polyunsaturated Fatty Acid Interactions with the Shaker Voltage-Gated Potassium Channel

Voltage-gated potassium (K_V) channels are membrane proteins that respond to changes in membrane potential by enabling K⁺ ion flux across the membrane. Polyunsaturated fatty acids (PUFAs) induce channel opening by modulating the voltage-sensitivity, which can provide effective treatment against refractory epilepsy by means of a ketogenic diet. While PUFAs have been reported to influence the gating mechanism by electrostatic interactions to the voltage-sensor domain (VSD), the exact PUFA-protein interactions are still elusive. In this study, we report on the interactions between the Shaker K_V channel in open and closed states and a PUFA-enriched lipid bilayer using microsecond molecular dynamics simulations. We determined a putative PUFA binding site in the open state of the channel located at the protein-lipid interface in the vicinity of the extracellular halves of the S3 and S4 helices of the VSD. In particular, the lipophilic PUFA tail covered a wide range of non-specific hydrophobic interactions in the hydrophobic central core of the protein-lipid interface, while the carboxylic head group displayed more specific interactions to polar/charged residues at the extracellular regions of the S3 and S4 helices, encompassing the S3-S4 linker. Moreover, by studying the interactions between saturated fatty acids (SFA) and the Shaker K_V channel, our study confirmed an increased conformational flexibility in the polyunsaturated carbon tails compared to saturated carbon chains, which may explain the specificity of PUFA action on channel proteins.     

Fluid phase lipid areas and bilayer thicknesses of commonly used phosphatidylcholines as a function of temperature

The structural parameters of fluid phase bilayers composed of phosphatidylcholines with fully saturated, mixed, and branched fatty acid chains, at several temperatures, have been determined by simultaneously analyzing small-angle neutron and X-ray scattering data. Bilayer parameters, such as area per lipid and overall bilayer thickness have been obtained in conjunction with intrabilayer structural parameters (e.g. hydrocarbon region thickness). The results have allowed us to assess the effect of temperature and hydrocarbon chain composition on bilayer structure. For example, we found that for all lipids there is, not surprisingly, an increase in fatty acid chain trans-gauche isomerization with increasing temperature. Moreover, this increase in trans-gauche isomerization scales with fatty acid chain length in mixed chain lipids. However, in the case of lipids with saturated fatty acid chains, trans-gauche isomerization is increasingly tempered by attractive chain-chain van der Waals interactions with increasing chain length. Finally, our results confirm a strong dependence of lipid chain dynamics as a function of double bond position along fatty acid chains.     

Synthesis and polymorphic phase behaviour of polyunsaturated phosphatidylcholines and phosphatidylethanolamines

A series of phosphatidylcholines and phosphatidylethanolamines was synthesized containing two acyl chains of the following polyunsaturated fatty acids: linoleic acid (18:2), linolenic acid (18:3), arachidonic acid (20:4) and docosahexaenoic acid (22:6). In addition two phospholipids with mixed acid composition were synthesized: 16:0/18:1_c phosphatidylcholine and 16:0/18:1_c phosphatidylethanolamine. The structural properties of these lipids in aqueous dispersions in the absence and in the presence of equimolar cholesterol were studied using ³¹P-NMR, freeze fracturing and differential scanning calorimetry (DSC).The phosphatidylcholines adopt a bilayer configuration above 0°C. Incorporation of 50 mol% of cholesterol in polyunsaturated species induces a transition at elevated temperatures into structures with ³¹P-NMR characteristics typical of non-bilayer organizations. When the acyl chains contain three or more double bonds, this non-bilayer organization is most likely the hexagonal H_II phase, 16:0/15:1_c phosphatidylethanolamine shows a bilayer to hexagonal transition temperature of 75°C. The polyunsaturated phosphatidylethanolamines exhibit a bilayer to hexagonal transition temperature below 0°C which decreases with increasing unsaturation and which is lowered by approximately 10°C upon incorporation of 50 mol% of cholesterol. Finally, it was found that small amounts of polyunsaturated fatty acyl chains in a phosphatidylethanolamine disproportionally lower its bilayer to hexagonal transition temperature.     

Molecular dynamics simulation of dipalmitoylphosphatidylcholine membrane with cholesterol sulfate

Using the molecular dynamics simulation technique, we studied the changes occurring in a dipalmitoylphosphatidylcholine (DPPC):cholesterol (CH) membrane at 50 mol% sterol when cholesterol is replaced with cholesterol sulfate (CS). Our simulations were performed at constant pressure and temperature on a nanosecond time scale. We found that 1) the area per DPPC:CS heterodimer is greater than the area of the DPPC:CH heterodimer; 2) CS increases ordering of DPPC acyl chains, but to a lesser extent than CH; 3) the number of hydrogen bonds between DPPC and water is decreased in a CS-containing membrane, but CS forms more water hydrogen bonds than CH; and 4) the membrane dipole potential reverses its sign for a DPPC-CS membrane compared to a DPPC-CH bilayer. We also studied the changes occurring in lipid headgroup conformations and determined the location of CS molecules in the membrane. Our results are in good agreement with the data available from experiments.     

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

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