Raman spectroscopic study of polycrystalline mono- and polyunsaturated 1-eicosanoyl-d(39)-2-eicosenoyl-sn-glycero-3-phosphocholines: bilayer lipid clustering and acyl chain order and disorder characteristics

Polycrystalline lipid samples of a series of mono- and polyunsaturated, double bond positional isomers of 1-eicosanoyl-d(39)-2-eicosenoyl-sn-glycero-3-phosphocholines [C(20-d(39)):C(20:1 Delta(j))PC, with j = 5, 8, 11, or 13; C(20-d(39)):C(20:2 Delta(11,14))PC; and C(20-d(39)):C(20:3 Delta(11, 14,17))PC] were investigated using vibrational Raman spectroscopy to assess the acyl chain packing order-disorder characteristics and putative bilayer cluster formation of the isotopically differentiated acyl chains. Perdeuteration of specifically the saturated sn-1 acyl chains for these bilayer systems enables each chain's intra- and intermolecular conformational and organizational properties to be evaluated separately. Various saturated chain methylene CD(2) and carbon-carbon (C&bond;C) stretching mode peak height intensity ratios and line width parameters for the polycrystalline samples demonstrate a high degree of sn-1 chain order that is unaffected by either the double bond placement or number of unsaturated bonds within the sn-2 chain. In contrast, the unsaturated sn-2 chain spectral signatures reflect increasing acyl chain conformational disorder as either the cis double bond is generally repositioned toward the chain terminus or the number of double bonds increases from one to three. The lipid bilayer chain packing differences observed between the sn-1 and sn-2 chains of this series of monounsaturated and polyunsaturated 20 carbon chain lipids suggest the existence of laterally distributed microdomains predicated on the formation of highly ordered, saturated sn-1 chain clusters.     

Molecular dynamics simulations of unsaturated lipid bilayers: effects of varying the numbers of double bonds

The importance of unsaturated, and especially polyunsaturated phosphatidylcholine molecules for the functional properties of biological membranes is widely accepted. Here, the effects of unsaturation on the nanosecond-scale structural and dynamic properties of the phosphatidylcholine bilayer were elucidated by performance of multinanosecond molecular dynamics simulations of all-atom bilayer models. Bilayers of dipalmitoylphosphatidylcholine and its mono-, di-, and tetraunsaturated counterparts were simulated, containing, respectively, oleoyl, linoleoyl, or arachidonoyl chains in the sn-2 position. Analysis of the simulations focused on comparison of the structural properties, especially the ordering of the chains in the membranes. Although the results suggest some problems in the CHARMM force field of the lipids when applied in a constant pressure ensemble, the features appearing in the ordering of the unsaturated chains are consistent with the behaviour known from ²H NMR experiments. The rigidity of the double bonds is compensated by the flexibility of skew state single bonds juxtaposed with double bonds. The presence of double bonds in the sn-2 chains considerably reduces the order parameters of the CH bonds. Moreover, the double bond region of tetraunsaturated chains is shown to span all the way from the bilayer centre to the head group region.     

Implicit solvent model estimates of the stability of model structures of the alamethicin channel

Alamethicin is a hydrophobic helical peptide of 20 residues, which oligomerizes to form ion-conducting channels in membranes. The behavior of an intact alamethicin channel in POPC bilayers was recently studied, using 2 ns molecular dynamics (MD) simulations of a model hexameric channel. These simulations produced numerous conformations of the channel. In the present study, we used 11 of these channel conformations and carried out continuum-solvent model calculations, similar to those used for the monomers in our previous studies, to investigate the energetics of the channel inside the lipid bilayer. Our results suggest that, out of the 11 channel conformations produced by the MD simulations, only four are stable inside the lipid bilayer, with water-to-membrane free energies of transfer ranging from ~–6 to ~–10 kcal/mol. Analysis of the results suggests two causes for the apparent instability of the remainder of the structures inside the lipid bilayer, both resulting from the desolvation of channel polar groups (i.e. their transfer from the aqueous phase into the bilayer). The first is specific, uncompensated backbone hydrogen bonds, which exist in the region of the channel exposed to the hydrocarbon of the lipid bilayer. The second is exposure of intra-pore water molecules to the surrounding lipid. Thus, the association of these structures with the membrane involves a large electrostatic desolvation free-energy penalty. The apparent conflict between continuum-solvent and MD calculations, and its significance for the interpretation of membrane proteins simulations, are discussed.     

Molecular motions and dynamics of a diunsaturated acyl chain in a lipid bilayer: implications for the role of polyunsaturation in biological membranes.

The nature and dynamics of the motions of a diunsaturated fatty acyl chain in a lipid bilayer were examined using a comprehensive simulation program for 2H NMR line shapes developed by Wittebort et al. [Wittebort, R. J., Olejniczak, E. T., & Griffin, R. G. (1987) J. Chem. Phys. 36, 5411-5420]. A motional model in which the isolinoleoyl chain (18:2 delta 6,9) adopts two conformations consistent with the low energy structures proposed for 1,4-pentadiene [Applegate, K. R., & Glomset, J. A. (1986) J. Lipid Res. 27, 658-680], but undergoes a rapid jump between these states, is sufficient to account for the experimentally observed quadrupolar couplings, the 2H-2H and 1H-2H dipolar couplings, the longitudinal relaxation times, and the changes in the average conformation of the chain that occur with a variation in temperature. The jump motion originates via rotations about the C7-C8 and the C8-C9 carbon bonds and leads to the low order parameters assigned to the C8 methylene segment (0.18) and the C9-C10 double bond (0.11). In contrast, the C6-C7 double bond, which is not involved in the two-site jump, characterized by a relatively large order parameter (0.56). Fatty acyl chains containing three or more double bonds likely cannot undergo the same jump motion and consequently will be highly ordered structures. Correlation times for diffusion of the molecular long axis of the diunsaturated acyl chain about the bilayer normal (approximately 10(-10) s) and for the local jump motion (approximately 10(-10) s) were calculated.(ABSTRACT TRUNCATED AT 250 WORDS)     

Effects of various numbers and positions of cis double bonds in the sn-2 acyl chain of phosphatidylethanolamine on the chain-melting temperature

In an attempt to investigate systematically the effects of various single and multiple cis carbon-carbon double bonds in the sn-2 acyl chains of natural phospholipids on membrane properties, we have de novo synthesized unsaturated C20 fatty acids comprised of single or multiple methylene-interrupted cis double bonds. Subsequently, 15 molecular species of phosphatidylethanolamine (PE) with sn-1 C20-saturated and sn-2 C20-unsaturated acyl chains were semi-synthesized by acylation of C20-lysophosphatidylcholine with unsaturated C20 fatty acids followed by phospholipase D-catalyzed base-exchange reaction in the presence of excess ethanolamine. The gel-to-liquid crystalline phase transitions of these 15 mixed-chain PE, in excess H2O, were investigated by high resolution differential scanning calorimetry. In addition, the energy-minimized structures of these sn-1 C20-saturated/sn-2 C20-unsaturated PE were simulated by molecular mechanics calculations. It is shown that the successive introduction of cis double bonds into the sn-2 acyl chain of C(20):C(20)PE can affect the gel-to-liquid crystalline phase transition temperature, Tm, of the lipid bilayer in some characteristic ways; moreover, the effect depends critically on the position of cis double bonds in the sn-2 acyl chain. Specifically, we have constructed a novel Tm diagram for the 15 species of unsaturated PE, from which the effects of the number and the position of cis double bonds on Tm can be examined simultaneously in a simple, direct, and unifying manner. Interestingly, the characteristic Tm profiles exhibited by different series of mixed-chain PE with increasing degree of unsaturation can be interpreted in terms of structural changes associated with acyl chain unsaturation.     

Thermodynamic analysis of chain-melting transition temperatures for monounsaturated phospholipid membranes: dependence on cis-monoenoic double bond position.

Unsaturated phospholipid is the membrane component that is essential to the dynamic environment needed for biomembrane function. The dependence of the chain-melting transition temperature, T(t), of phospholipid bilayer membranes on the position, n(u), of the cis double bond in the glycerophospholipid sn-2 chain can be described by an expression of the form T(t) = T(t)(c)(1 + h'(c)|n(u) - n(c)|)/(1 + s'(c)|n(u) - n(c)|), where n(c) is the chain position of the double bond corresponding to the minimum transition temperature, T(t)(c), for constant diacyl lipid chain lengths. This implies that the incremental transition enthalpy (and entropy) contributed by the sn-2 chain is greater for whichever of the chain segments, above or below the double-bond position, is the longer. The critical position, n(c), of the double bond is offset from the center of the sn-2 chain by an approximately constant amount, deltan(c) approximately 1. 5 C-atom units. The dependence of the parameters T(t)(c), h'(c), and s'(c) on sn-1 and sn-2 chain lengths can be interpreted consistently when allowance is made for the chain packing mismatch between the sn-1 and sn-2 chains. The length of the sn-2 chain is reduced by approximately 0.8 C-atom units by the cis double bond, in addition to a shortening by approximately 1.3 C-atom units by the bent configuration at the C-2 position. Based on this analysis, a general thermodynamic expression is proposed for the dependence of the chain-melting transition temperature on the position of the cis double bond and on the sn-1 and sn-2 chain lengths. The above treatment is restricted mostly to double-bond positions close to the center of the sn-2 chain. For double bonds positioned closer to the carboxyl or terminal methyl ends of the sn-2 chain, the effects on transition enthalpy can be considerably larger. They may be interpreted by the same formalism, but with different characteristic parameters, h'(c) and s'(c), such that the shorter of the chain segments makes a considerably smaller contribution to the calorimetric properties of the chain-melting transition.     

Average structural and motional properties of a diunsaturated acyl chain in a lipid bilayer: effects of two cis-unsaturated double bonds

Isolinoleic acid (18:2 delta 6,9) deuterated at 10 different positions was esterified to form 1-palmitoyl-2-isolinoleoyl-sn-glycero-3-phosphocholine (PiLPC), and the average structural and motional properties of the diunsaturated chain, in aqueous dispersions of PiLPC, were examined by 2H NMR spectroscopy. For each sample, 2H spectra were acquired over a temperature range of 1-40 degrees C and the quadrupolar splittings interpreted in terms of carbon-deuterium bond order parameters, SCD. Furthermore, definition of the average orientation of the C8 methylene unit with respect to the bilayer normal [Baenziger, J. E., Smith, I. C. P., Hill, R. J., & Jarrell, H. C. (1988) J. Am. Chem. Soc. 110, 8229-8231] provided sufficient information to calculate both the average orientations and the molecular order parameters, Smol (which reflects the amplitudes of motion), for the C6-C7 and the C9-C10 double bonds. The results indicate that both the motional freedom (reflected in the order profile) and the average structure (reflected in the orientation of carbon segments with respect to the bilayer normal) are strongly affected by the presence of two cis-unsaturated double bonds. The data were interpreted in terms of two possible models whereby, in each case, the chain adopts a conformation consistent with the low-energy conformation of 1,4-pentadiene [Applegate, K. R., & Glomset, J. A. (1986) J. Lipid Res. 27, 658-680] but undergoes a two-site jump between the conformations. The jump motion arises mainly from rotations about the C7-C8 and the C8-C9 single bonds that disorder the C8 and the C9-C10 segments (Smol = 0.15 and 0.08, respectively) but leave the C6-C7 double bond relatively immobile (Smol = 0.55; all at 40 degrees C). It is suggested that acyl chains containing three or more double bonds could not undergo a similar jump motion and therefore would be highly ordered and not "fluid" as is generally thought.     

Calorimetric and molecular mechanics studies of the thermotropic phase behavior of membrane phospholipids.

In this review, we summarize the results of recent studies on the main phase transition behavior of phospholipid bilayers using the combined approaches of molecular mechanics simulations and high-resolution differential scanning calorimetry. Following a brief overview of the phase transition phenomenon exhibited by the lipid bilayer, we begin with the review by showing how several structural parameters underlying various phospholipids including phosphatidylcholine, phosphatidylethanolamine, and phosphatidylglycerol are defined and determined. Specifically, these structural parameters are obtained with saturated lipids packed in the gel-state bilayer using computer-based molecular mechanics calculations. Then we proceed to present the calorimetric data obtained with the lipid bilayer composed of saturated phospholipids as it undergoes the gel-to-liquid-crystalline phase transition in excess water. The general equations that can correlate the gel-to-liquid-crystalline phase transition temperature (T(m)) of the lipid bilayer with the structural parameters of the lipid molecule constituting the lipid bilayer are subsequently presented. From these equations, two tables of predicated T(m) values for well over 400 molecular species of saturated phosphatidylcholine and saturated phosphatidylethanolamine are generated. We further review the structure and chain-melting behavior of a large number of sn-1 saturated/sn-2 unsaturated phospholipids. Two T(m)-diagrams are shown, from which the effects of the number and the position of one to five cis carbon-carbon double bonds on T(m) can be viewed simultaneously. Finally, in the last part of this review, simple molecular models that have been invoked to interpret the characteristic T(m) trends exhibited by lipid bilayers composed of unsaturated lipids with different numbers and positions of cis carbon-carbon double bonds as seen in the T(m)-diagram are presented.     

Changes in a phospholipid bilayer induced by the hydrolysis of a phospholipase A2 enzyme: a molecular dynamics simulation study

Phospholipase A2 (PLA2) enzymes are important in numerous physiological processes. Their function at lipid-water interfaces is also used as a biophysical model for protein-membrane interactions. These enzymes catalyze the hydrolysis of the sn-2 bonds of various phospholipids and the hydrolysis products are known to increase the activity of the enzymes. Here, we have applied molecular dynamics (MD) simulations to study the membrane properties in three compositionally different systems that relate to PLA2 enzyme action. One-nanosecond simulations were performed for a 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) bilayer and for two of its PLA2-hydrolyzed versions, i.e., bilayers consisting of lysophospholipids and of either free charged linoleate or free uncharged linoleic acid molecules. The results revealed loosening of the structure in the hydrolyzed bilayer due to increased mobility of the molecules in the direction normal to the bilayer. This loss of integrity due to the hydrolysis products is in accord with observations that not only the presence of hydrolysis products, but also a variety of other perturbations of the membrane may activate PLA2. Additionally, changes were observed in other structural parameters and in the electrostatic potential across the membrane-water interface. These changes are discussed in relation to the simulation methodology and the experimental observations of PLA2-hydrolyzed membranes.     

Distribution of Pentachlorophenol in Phospholipid Bilayers: A Molecular Dynamics Study

Molecular dynamics computer simulations of pentachlorophenol (PCP) in palmitoyl-oleoyl-phosphatidylethanolamine and palmitoyl-oleoyl-phosphatidylcholine lipid bilayers were carried out to investigate the distribution of PCP and the effects of PCP on the phospholipid bilayer structure. Starting from two extreme starting structures, including PCP molecules outside the lipid bilayer, the PCP distribution converges in simulations of up to 50 ns. PCP preferentially occupies the region between the carbonyl groups and the double bonds in the acyl chains of the lipid molecules in the bilayer. In the presence of PCP, the lipid chain order increases somewhat in both chains, and the average tilt angle of the lipid chains decreases. The increase in the lipid chain order in the presence of PCP was more pronounced in the palmitoyl-oleoyl-phosphatidylcholine bilayer compared to the palmitoyl-oleoyl-phosphatidylethanolamine bilayer. The number of trans conformations of lipid chain dihedrals does not change significantly. PCP aligns parallel to the alkyl chains of the lipid to optimize the packing in the dense ordered chain region of the bilayer. The hydroxyl group of PCP forms hydrogen bonds with both water and lipid oxygen atoms in the water/lipid interface region.     

On the rate of genetic adaptation under natural selection.

Unsaturated fatty acids are constituents of nearly all biological membranes. They are always present in membranes which possess transmembrane potentials. Two completely different biosynthetic routes have evolved (aerobic and anaerobic) for placing cis double bonds in the 9 position on the fatty acids of membrane lipids. Bacterial membranes contain primarily monounsaturated fatty acids, whereas eukaryote membranes contain a significant fraction of polyunsaturated fatty acids. The polyunsaturated fatty acids are concentrated in organelles, such as chloroplasts and mitochondria that are known to manipulate transmembrane potentials. I propose that the function of the unsaturated fatty acids is to facilitate the transmission of a local compaction of the membrane (in response to a transmembrane potential) laterally through the membrane. The role of the cis double bond at position 9 is twofold: first to create a kink in the chains of a large fraction of membrane fatty acids enhancing the separation of two regions in the membrane and second to enhance the rigidity of the membrane in the region between the head group and the 9 double bond. This ordered region contains those carbons proximal to the 9 carbon and which are in a regular array of trans conformations. The presence of a reasonable proportion of cis double bonds at position 9 will tend to maintain these trans conformations utilizing pi-pi (van der Waals) interactions between adjacent hydrocarbon chains at position 9. The disordered region contains the carbons distal to the 9 carbon. These have greater degrees of freedom and considerable gauche conformations. The role of the double bonds in the polyunsaturated fatty acids distal to carbon 9 is to facilitate trans bilayer pi-pi (van der Waals) interactions enhancing compaction of the bilayer during the electrostriction. I further propose that it is the function of the ionic headgroups to form an interlocking polyionic network which constitutes an elastic sheet. These ionic interactions would serve as the restoring force converting the compaction into a wave. The facilitation of the compaction of the bilayer together with the polyionic restoring force permits the membrane to transmit conformational changes from one transmembrane protein to another. Since transmembrane potentials are created and responded to by proteins each in a single location, it is thus proposed that a potential compaction wave emanates from the first protein in all directions in the plane of the membrane. The proposed wave would have both physical and electrical components. The electrohydrodynamic wave would require that the compaction oscillations be coupled to an oscillating electrical field. These proposals are applied to mitochondrial oxidative phosphorylation, and to transport across biological membranes.     

Orientation of specifically 13C=O labeled phosphatidylcholine multilayers from polarized attenuated total reflection FT-IR spectroscopy.

Oriented multilayers of 1-myristoyl-2(1-13C)-myristoyl-sn-glycero-3-phosphatidylcholine (2[1-13C]DMPC) and 1-palmitoyl-2(1-13C)-palmitoyl-sn-glycero-3-phosphatidylcholine (2[1-13C]DPPC) were investigated by use of attenuated total reflection infrared spectroscopy with polarized light. Experiments were performed with the aim to determine the orientation of the two ester groups in these phospholipids in the solid state and in the hydrated state at temperatures below and above the respective gel to liquid-crystalline phase transitions. Substitution of the naturally occurring 12C carbonyl carbon atom by 13C in the ester group of the sn-2 chain of DMPC and DPPC shifts the infrared absorption of the carbonyl double bond stretching vibration to lower frequency. This results in two well-resolved ester C=O bands which can be assigned unequivocally to the sn-1 and sn-2 chains as they are separated by more than 40 cm-1. The two ester CO-O single bond stretching vibrations of the molecular fragments-CH2CO-OC-are also affected and the corresponding infrared absorption band shifts by 20 cm-1 on 13C-labeling of the carbonyl carbon atom. From the dichroic ratios of the individual ester bands in 2(1-13C)DMPC and 2(1-13C)DPPC we were able to demonstrate that the sn-1 and sn-2 ester C=O groups are similarly oriented with respect to the bilayer plane, with an angle greater than or equal to 60 degrees relative to the bilayer normal. The two CO-O single bonds on the other hand have very different orientations. The CH2CO-OC fragment of the sn-1 chain is oriented along the direction of the all-trans methylene chain, whereas the same molecular segment of the sn-2 carbon chain is directed toward the bilayer plane. This orientation of the ester groups is retained in the liquid-crystalline phase. The tilt angle of the hydrocarbon all-trans chains, relative to the membrane normal, is 25 degrees in the solid state of DMPC and DPPC multibilayers. In the hydrated gel state this angle varies between 26 degrees and 30 degrees, depending on temperature. Neither the orientation of the phosphate group, nor that of the choline group varies significantly in the different physical states of these phospholipids.     

Molecular Dynamics Simulations of Phospholipid Bilayers

Abstract

Molecular dynamics (MD) simulations at 37°C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a_t/a_h), correctly predict the aggregation preferences of the three lipids observed experimentally at 37°C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.     

2H-NMR study and molecular dynamics simulation of the location, alignment, and mobility of pyrene in POPC bilayers

The alignment of pyrene in a 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer was investigated using two different approaches, namely solid-state (2)H-NMR spectroscopy and molecular dynamics (MD) simulations. Quadrupolar splittings from (2)H-NMR spectra of deuterated pyrene-d(10) in an oriented lipid bilayer give information about the orientation of C-D bonds with respect to the membrane normal. From MD simulations, geometric information is accessible via trajectories. By defining molecular and bond order parameters, the data from MD trajectories and NMR spectra can be compared straightforwardly. To ensure that the results from both methods are comparable, parameters of the experimental and the simulation setup were chosen to be as similar as possible. From simulations, we saw that pyrene prefers a position inside the lipid membrane near the headgroups and has no tendency to diffuse from one monolayer of the membrane to the other. The results from simulation and NMR show that the normal of the molecular plane is aligned nearly perpendicular to the bilayer normal. The long axis of pyrene lies preferentially parallel to the bilayer normal within a range of +/-30 degrees . The results from the two different methods are remarkably consistent. The good agreement can be explained by the fact that the different kind of motions of a pyrene molecule are already averaged within a few nanoseconds, which is the timescale covered by the MD simulation.     

Effect of acyl chain unsaturation on the packing of model diacylglycerols in simulated monolayers.

In a companion study of the effects of acyl chain unsaturation on a series of model sn-1,2-diacylglycerols (DGs) we showed that individual DGs could adopt one of three energy-minimized conformations depending on the number and location of cis double bonds in the sn-2 chain. Here we show that each of these conformations promotes a distinct type of packing arrangement in a simulated DG monolayer. One conformation, shown by sn-1-18:0 DGs containing an sn-2 22:6(n-3)-, 20:4(n-6)-, or 20:3(n-9)- group, determines a regular packing that resembles a known hybrid subcell, HS2, of crystalline hydrocarbon chains. The second conformation, shown by DGs containing an sn-2 18:0-, 18:2(n-6)-, or 18:3(n-3)- group, determines a regular packing that resembles a second known, distinct hydrocarbon subcell, HS1. The third conformation, that of 18:0/18:1(n-9) DG, determines a much looser, less energetically favorable packing. Stable heterogeneous packings are possible for DGs that have similar conformations, but mixed packings of DGs that have dissimilar conformations are less stable. These results raise the possibility that differences in sn-2 acyl chain unsaturation among membrane sn-1,2-diacylglycerophospholipids may promote the formation of different domains.     

Influence of cis double bonds in the sn-2 acyl chain of phosphatidylethanolamine on the gel-to-liquid crystalline phase transition.

We have semisynthesized 19 species of mixed-chain phosphatidylethanolamines (PEs) in which the sn-1 acyl chain is derived from saturated fatty acids with varying chain lengths and the sn-2 acyl chain has different chain lengths but contains 0, 1, and 2 cis double bond(s). The gel-to-liquid crystalline phase transition temperatures (Tm) of lipid bilayers prepared from these 19 mixed-chain PEs were determined calorimetrically. When the Tm values are compared with those of saturated and monounsaturated counterparts, a common Tm profile is observed in the plot of Tm versus the number of cis double bonds. Specifically, a marked stepwise decrease in Tm is detected as the number of cis double bonds in the sn-2 acyl chain of the mixed-chain PE is successively increased from 0 to 1 and then to 2. The large Tm-lowering effect of the acyl chain unsaturation can be attributed to the increase in Gibbs free energy of the gel-state bilayer as a result of weaker lateral chain-chain interactions. In addition, we have applied molecular mechanics calculations to simulate the molecular structure of dienoic mixed-chain C(X):C(Y:2 delta n,n+3)PE in the gel-state bilayer, thus enabling the three independent structural parameters (N, delta C, and LS) to be calculated in terms of X, Y, and n, which are intrinsic quantities of C(X):C(Y:2 delta n,n+3)PE. When the Tm values and the corresponding N and delta C values of all dienoic mixed-chain PEs under study are first codified and then analyzed statistically by multiple regressions, the dependence of Tm on the structural parameters can be described quantitatively by a simple and general equation. The physical meaning and the usefulness of this simple and general equation are explained.     

Phosphatidylcholine fluidity and structure affect lecithin:cholesterol acyltransferase activity

The purpose of this study was to test the hypothesis that lipid fluidity regulates lecithin:cholesterol acyltransferase (LCAT) activity. Phosphatidylcholine (PC) species were synthesized that varied in fluidity by changing the number, type (cis vs. trans), or position of the double bonds in 18 or 20 carbon sn-2 fatty acyl chains and recombined with [(3)H]cholesterol and apolipoprotein A-I to form recombinant high density lipoprotein (rHDL) substrate particles. The activity of purified human plasma LCAT decreased with PC sn-2 fatty acyl chains containing trans versus cis double bonds and as double bonds were moved towards the methyl terminus of the sn-2 fatty acyl chain. The decrease in LCAT activity was significantly correlated with a decrease in rHDL fluidity (measured by diphenylhexatriene fluorescence polarization) for PC species containing 18 carbon (r(2) = 0.61, n = 18) and 20 carbon (r(2) = 0.93, n = 5) sn-2 fatty acyl chains. rHDL were also made containing 10% of the 18 carbon sn-2 fatty acyl chain PC species and 90% of an inert PC ether matrix (sn-1 18:1, sn-2 16:0 PC ether) to normalize rHDL fluidity. Even though fluidity was similar among the PC ether-containing rHDL, the order of PC reactivity with LCAT was significantly correlated (r(2) = 0.71) with that of 100% PC rHDL containing the same 18 carbon sn-2 fatty acyl chain species, suggesting that PC structure in the active site of LCAT determines reactivity in the absence of measurable differences in bilayer fluidity. We conclude that PC fluidity and structure are major regulators of LCAT activity when fatty acyl chain length is constant.     

Acyl chain dynamics of phosphatidylethanolamines containing oleic acid and dihydrosterculic acid: 2H NMR relaxation studies

The dynamical behavior of the acyl chains of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and 1-palmitoyl-2-dihydrosterculoyl-sn-glycero-3-phosphoethanolamine has been investigated by using 2H T1 and T2 relaxation times. Lipids were labeled at the 5-,9-,10-, and 16-positions of the sn-2 acyl chain. The profile of deuterium spin-lattice relaxation rate (T1(-1) vs. chain position is characterized in all systems by a marked discontinuity at the positions of the carbon-carbon double bond and the cyclopropane ring; the deuterons at these positions have relaxation rates which are greater than at any other labeled position of the sn-2 chain. For both types of sn-2 acyl chain, assuming a single-exponential correlation time and that the motion is within the rapid regime, the phosphatidylcholine lipid systems are less mobile than their phosphatidylethanolamine analogues. Systems containing an oleoyl chain are more dynamic than their analogues containing a dihydrosterculoyl chain. The rates of motion of the sn-2 acyl chains of phosphatidylethanolamine in a bilayer structure are slower than those of the lipid in an inverted hexagonal structure. In the hexagonal phase, the motional rates of a dihydrosterculoyl chain are slower than those of the corresponding positions of an oleoyl chain.     

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.     

Effects of phospholipid unsaturation on the membrane/water interface: a molecular simulation study.

Molecular dynamics (MD) simulations of fully hydrated bilayers in the liquid-crystalline state made of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) or 1-palmitoyl-2-elaidoyl-phosphatidylcholine (PEPC) were carried out to investigate the effect of the incorporation of a double bond in the phosphatidylcholine (PC) beta-chain (cis or trans) on the membrane/water interface. The bilayers reached thermal equilibrium after 3 and 1 ns of MD simulations, respectively, and productive runs were carried out for 3 ns for each bilayer. As reference systems, the 1,2-dimyristoyl-phosphatidylcholine (DMPC) bilayer (M. Pasenkiewicz-Gierula, Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi, 1999, Biophys. J. 76:1228-1240) and DMPC-cholesterol (Chol) bilayer containing 22 mol % Chol (M. Pasenkiewicz-Gierula, T. Róg, K. Kitamura, A. and Kusumi, 2000, Biophys. J. 78:1376-1389) were used. The study shows that at the interface of POPC, PEPC, and DMPC-Chol bilayers, average numbers of PC-water and PC-PC interactions are similar and, respectively, greater and smaller than in the DMPC bilayer. The average area/PC in mono-unsaturated bilayers is approximately 4 A(2) larger than in the DMPC bilayer; nevertheless, a strong correlation was found between a single molecular area (SMA) of a PC and the number of interactions this PC makes; i.e., PCs (either saturated or unsaturated) with the same SMA form similar numbers of intermolecular links. The numbers and corresponding SMAs are distributed about averages pertinent to each bilayer. No significant difference between cis and trans bonds was found.