Biophysical Journal Volume 99 November 2010 2957–2966
2957
Importance of the Sphingosine Base Double-Bond Geometry for the Structural and Thermodynamic Properties of Sphingomyelin Bilayers Lorant Janosi* and Alemayehu Gorfe* Department of Integrative Biology and Pharmacology, University of Texas-Health Science Center, Houston, Texas
ABSTRACT The precise role of the sphingosine base trans double bond for the unique properties of sphingomyelins (SMs), one of the main lipid components in raftlike structures of biological membranes, has not been fully explored. Several reports comparing the hydration, lipid packing, and hydrogen-bonding behaviors of SM and glycerophospholipid bilayers found remarkable differences overall. However, the atomic interactions linking the double-bond geometry with these thermodynamic and structural changes remained elusive. A recent report on ceramides, which differ from SMs only by their hydroxyl headgroup, has shown that replacing the trans double bond of the sphingosine base by cis weakens the hydrogen-bonding potential of these lipids and thereby alters their biological activity. Based on data from extensive (a total 0.75 ms) atomistic molecular dynamics simulations of bilayers composed of all-trans, all-cis, and a trans/cis (4:1 ratio) racemic mixture of sphingomyelin lipids, here we show that the trans configuration allows for the formation of significantly more hydrogen bonds than the cis. The extra hydrogen bonds enabled tighter packing of lipids in the all-trans and trans/cis bilayers, thus reducing the average area per lipid while increasing the chain order and the bilayer thickness. Moreover, fewer water molecules access the lipid-water interface of the all-trans bilayer than of the all-cis bilayer. These results provide the atomic basis for the importance of the natural sphingomyelin trans double-bond conformation for the formation of ordered membrane domains.
INTRODUCTION Sphingomyelins (SMs) are one of the most abundant phospholipids (1,2) important for the formation of ordered lipid domains and rafts in the cell membranes of vertebrates (3–9). Ordered membrane domains play critical roles in many signaling processes (10–13). Moreover, SM metabolic products, such as ceramides, sphingosines, and sphingosine 1-phosphates, are involved in the regulation of cell proliferation, differentiation, and apoptosis (14–17), as well as in endocytosis and inflammation (16–20). SM has a phosphocholine polar headgroup, an N-acyl chain, and a base. The most common base in mammalian SM is sphingosine (1,3-dihydroxy-2-amino-4-octadecene), which contains a trans double bond between C4 and C5 (Fig. 1). Recent studies found that changing the natural trans double bond to cis or to a saturated bond has strong influence on the structural properties of sphingomyelin and ceramide (Cer) bilayers (21–23). (Cer differs from SM only by its hydroxyl headgroup.) Although similar to diacyl-glycero-phosphatidylcholines (diacyl-PCs from here on) in many regards, SMs exhibit a more complex behavior (24). For example, four distinct gel phases have been identified for 18:0-SM (25) whereas SMs with large mismatch in length between their N-acyl and sphingosine chains have multiple transitions (26). Most natural SMs exhibit a nonlinear dependence on the acyl-chain length (27), with their main phase transition temperature being close to the physiological temperature. For instance, for pure palmitoyl Submitted June 1, 2010, and accepted for publication September 13, 2010. *Correspondence:
[email protected] or alemayehu.g.abebe@uth. tmc.edu Editor: Thomas J. McIntosh. Ó 2010 by the Biophysical Society 0006-3495/10/11/2957/10 $2.00
SM (16:0-SM), the transition is ~314 K (26), which is similar to that of DPPC (di-16:0-PC). Compared with diacyl-PCs, insertion of cis double bonds in the acyl chain of SMs has little impact on their high melting temperatures. This is likely due to the intermolecular hydrogen-bonding networks that stabilize the sphingomyelin bilayers (27). Furthermore, despite the same headgroup, the hydration behavior of SM bilayers is significantly different from that of diacyl-PCs (28,29). The source of this difference is believed to involve hydrogen-bonding interactions, which in turn are affected by the stereochemistry and geometry of the double bond at the sphingosine base (30,31). For instance, studies on Cers found that saturation of the sphingosine trans double bond or its replacement with the cis isomer leads to the loss of the H-bond network between Cer’s hydroxyl groups (both the headgroup and C3OH) and the bound water molecules (21,23,32). To gain atomistic insight into the influence of the sphingosine double-bond geometry on the structural, dynamic, and thermodynamic properties of SM bilayers, we performed extensive molecular dynamics simulations of palmitoyl sphingomyelin (PSM) bilayers of pure trans (all-trans), pure cis (all-cis), and mixed trans and cis isomers of the sphingosine double bond. We found that the cis isomer significantly reduces the formation of intermolecular hydrogen bonds. This disrupts the high ordering of the sphingosine chain, resulting in reduced packing of the bilayer. The C4–C5 trans and cis unsaturations had opposing effects on the bilayer structures, similar to effects expected for trans/cis isomerization in the midsection of a tail: the trans double bond maintains large values of the order doi: 10.1016/j.bpj.2010.09.020
2958
Janosi and Gorfe
z
hydrophilic region 1 2 4
3
trans
interfacial region
cis
5
6
6 7
8
2
3 4 5 7
9 10 12
hydrophobic region
11 13
PA
14 15
SPH
16 18
PA
SPH
17
FIGURE 1 Atomistic structures of trans (left) and cis (right) isomers of the PSM lipid. PSM is composed of a phosphocholine headgroup, the sphingosine base (SPH), and a palmitoyl chain (PA). VMD (57) was used to render the images. The representation used for the PSM is the GROMOS96 united-atom model (36); i.e., only the hydrogen atoms of donors (like the O–H and N–H) are explicit.
parameter (i.e., close to the value of a saturated tail), whereas the cis isomer lowers it dramatically (33). We also found that the number of waters solvating the all-trans lipid-water interface is reduced by ~1.13 waters per lipid compared with the all-cis bilayer. These results are discussed in terms of their impact onto the elastic properties of SM bilayers. MODEL AND METHODS The structures of PSM trans and cis isomers are shown in Fig. 1. Each molecule consists of three main parts: the phosphocholine headgroup, the sphingosine (SPH) base, and the palmitoyl (PA) chain. The enantiomeric configuration of both molecules is D-erythro (2S, 3R). The first part of the SPH base (C1–C5) bridges the headgroup to the hydrophobic tails of the lipid, and contains a double bond between C4 and C5 as well as two groups capable of hydrogen bonding: the NH acting as donor and the OH both as donor and acceptor. Bilayers made up of PSM lipids with
the C4–C5 double bond in the trans and the cis conformations, as well as a racemic mixture of 80% trans and 20% cis, were investigated. To speed up system construction, we began from a previously fully equilibrated DOPC bilayer at 323 K (in the liquid-disordered phase). The cis and trans configurations were generated by imposing an improper dihedral potential around the double bond with minimum at 0 and 180 , respectively. The dihedral was used only during a standard steepest descent energy minimization procedure of 2500 steps. At the end of the minimization, we rigorously tested that the three systems have 1), all-trans, 2), all-cis, and 3), a mixture of 80% trans and 20% cis geometries of the sphingosine double bond. All three systems contained 128 PSM lipids fully solvated with 3735 water molecules, summing up to a total of 17,605 atoms. Each system was then subjected to 250 ns of free equilibration simulations in the NpT (constant temperature and pressure) ensemble at 323 K, and having a proper dihedral potential around the sphingosine double bond as described by Chiu et al. (34). The equilibrated all-trans and all-cis systems are shown in Fig. 2. All simulations were performed with GROMACS 4.0 (35) and the Gromos96 ffG43A1 force field (36). Parameters proposed by Chiu et al. (34) were used to describe interactions of the PSM lipids with themselves and the surrounding environment. The van der Waals interactions were evaluated up to 1.2 nm, with twin-range cutoff neighbor search approach, with a short-range neighbor-list radius of 1.0 nm. Coulomb interactions were evaluated using the particle-mesh Ewald method with short-range interaction cutoff of 1.0 nm. The bilayer and the water were separately coupled to a Nose´-Hoover heat bath. All systems were semiisotropically (in the membrane plane and perpendicular to it) coupled to a reference pressure of 1 bar with compressibility constant 4.6 105 bar1.
RESULTS AND DISCUSSION Double-bond isomerization alters lipid packing of PSM bilayer The structure of a lipid bilayer can be characterized by a number of equilibrium quantities, including the area per lipid (AL) and the bilayer thickness (dPP). Because both AL and dPP are directly related to lipid packing, we used them to examine the effect of the SM double-bond isomerization on the lateral and transverse organization of the alltrans, all-cis, and trans/cis bilayers. The time evolution of the areas per lipid AL(t), calculated as the ratio between the (fluctuating) area of the simulation box and the number of lipids in one leaflet, are shown in Fig. 3 a. Fig. 3 b displays the corresponding dPP (t), calculated as the average
FIGURE 2 Equilibrated all-trans (a) and all-cis (b) systems at 323 K. PSM lipids are represented in licorice green, phosphorus atoms in van der Waals tan spheres, and water in cyan continuous surface. Snapshots were rendered using VMD (57).
Biophysical Journal 99(9) 2957–2966
Role of Sphingosine Double Bond
2959 TABLE 1 List of simulated PSM bilayers with the reported equilibrated area per lipid molecule AL, simulation temperature T, and time t
a 0.58
80% trans + 20% cis pure trans pure cis
2
area / lipid [nm ]
0.56 0.54 0.52 0.50
50
100
150
200
4.2 4.1 P-P distance [nm]
T [K]
t [ns]
(40) (40) (40) (41) (42) (44)
0.565 0.573 0.599 0.515 0.520 0.525
325 325 325 310 323 323
8.0 8.0 8.0 5.4 50.0 50.0
250
b
4.0 3.9 3.8 3.7
80% trans + 20% cis pure trans pure cis
3.6 3.5 0
hALi [nm2]
Although not always explicitly stated, we presume that all of these simulations were done with the sphingosine base double bond in the trans conformation.
0.48 0
Ref.
50
100
150
200
250
time [ns]
FIGURE 3 Time evolution throughout the three simulations of the (a) area per lipid molecule and of the (b) lipid bilayer thickness measured as the average distance between the phosphorus atoms.
distance between the centers of mass of the phosphorous atoms at the two leaflets. Both AL(t) and dPP (t) equilibrated after ~40 ns, 110 ns, and 190 ns in the all-cis, alltrans, and mixed bilayer systems, respectively. The well-equilibrated portions of the trajectories (last 50 ns) yield average areas per lipid of hALitrans ¼ 0.485 nm2, hALicis ¼ 0.531 nm2, and hALitrans/cis ¼ 0.506 nm2. The experimental area per lipid for PSM bilayers varies between 0.643 nm2 at 323 K (37) and 0.47 nm2 at 328 K (38). This variation is in part due to differences in the experimental methods as well as in the models used to interpret the data (39). For example, when the area had to be derived from molecular volume and bilayer thickness, the specific volume of DPPC was used to approximate the unavailable volume of PSM (37,38). Several computer simulation studies of sphingomyelin bilayers have been reported (34,40–44). Different force-field parameters for SM and water led to variations in the calculated hALi (40–42,44) (Table 1), though all are within the broad range of the experimental results (37,38). In fact, the variations in the simulated values (including the current ones) are smaller than the variations in the experimental area per lipid.
A more direct comparison between the simulated and experimental data can be made for the bilayer thickness, because it is directly measured in experiments (37,38). The equilibrated portions of our trajectories yield thicktrans=cis cis ¼ nesses of dPtrans -P ¼ 4.16 nm, dP-P ¼ 3.84 nm, and dP-P 4.02 nm. The calculated bilayer thickness for the all-trans system (dPP ¼ 4.16 nm) matches reasonably well with the most recent experimental value reported by Maulik and Shipley (4.44 nm) (38). To evaluate the significance of the differences in the area and thickness of the trans and cis systems to the bilayers’ elastic properties, we computed the isothermal area compressibility modulus, KA, from the average area per lipid hALi and its variance s2AL (45,46), KA ¼
1 hAL ikB T ; NL s2AL
(1)
where NL is the number of lipids per leaflet. We obtained KAtrans ~2700 mN/m, which is more than twice larger than KAcis ~1200 mN/m. Although no experimental data is found for sphingomyelin, this result is in qualitative agreement with previously reported value for 18:0-SM simulated bilayer of KAtrans ¼ 4400 mN/m (34). However, insufficient fluctuations in the simulation box (possibly due to limitations of the force field) led to area compressibility moduli that are likely overestimated (47) for both the cis and trans isomers. Nonetheless, their relative values, which are most relevant for the purposes of this work, should hold. The data described above clearly show that the conformation of the sphingosine base double bond plays a key role in determining the equilibrium structure of the resulting bilayers. In the following sections, therefore, we will discuss the details of the atomic interactions that underlie this role. Organization of PSM lipids in bilayers The differences in dPP and AL indicate that the trans PSM forms a more compact and thicker bilayer than the cis. Much more detailed structural information can be obtained from the distribution of atoms along the bilayer transverse Biophysical Journal 99(9) 2957–2966
2960
Janosi and Gorfe
dimension. Fig. 4 shows the number density distributions for the most important lipid atoms (or groups of atoms) at the lipid-water interface, with solid lines for the all-trans and dashed lines for the all-cis. The reference is set to be the average position of the C2 carbon of the sphingosine base. For ease of comparison, we divided the lipid-water interface into four regions:
TABLE 2 Bilayer thicknesses measured as distances between the P atoms and C2 atoms and between the double bonds (average position of the C4 and C5 atoms)
I. Bulk water, defined by the bulk density of the water molecules. II. Headgroup region, between region I and the average position of the phosphate group. III. Interfacial region, between the phosphate group and the end of the wet region. IV. Dry lipid tails region, beyond the sphingosine double bond.
Symbol D is the thickness difference between the trans and cis systems.
Region III makes the transition from the hydrophilic regions I and II to the hydrophobic region IV. The C2 carbon of the sphingosine base lies approximately in the middle of region III and can be used to calculate the hydrophobic thickness of the bilayer. The z positions of the peaks for the distributions in Fig. 4 suggest that the thickness difference at the interfacial region between all-trans and all-cis is small, with ~80% of the total difference coming from the hydrophobic thickness (Table 2). The trans-PSM lipids are more stretched than the cis, a conclusion also supported by analysis of the chain order parameters (discussed later). The cause of the small difference in the thickness of the interfacial region, however, is not as straightforward. The compactness of the trans system allows for less water penetration into the interfacial region (see Ow in Fig. 4) and repositions the phosphate groups such that the C2P vectors are oriented ~6 closer to the OW C2 PO4 O1 PO4 O2 CO oxygen OH NH C4 & C5
water bulk density
-3
Number density [nm ]
30
20
I.
II.
III.
IV.
10
0 -1.2
-1.0
-0.8
-0.6 -0.4 -0.2 0.0 0.2 0.4 z relative to C2 avg. position [nm]
0.6
0.8
FIGURE 4 Comparison of number densities in the all-trans (solid) and all-cis (dashed) systems for selected atoms along the bilayer normal. The reference is the average position of the C2 carbon of the sphingosine base (see Fig. 1). Regions I and II are hydrophilic, corresponding to bulk water and lipid headgroup region; III and IV are the interfacial and hydrophobic regions. Biophysical Journal 99(9) 2957–2966
Thickness [nm]
trans cis D
dP–P
dC2–C2
ddb
4.16 3.84 0.32
3.52 3.25 0.27
3.07 2.81 0.26
bilayer normal. Hence, the phosphate atoms move further from the interface, slightly increasing the interfacial thickness of the lipid bilayer (see Fig. 1). Fig. 4 also indicates that the way the water density decays from the bulk value of >32 nm3 to zero in the hydrophobic region differs from a sigmoidal shape typically observed in diacyl-PCs (e.g., DPPC) (48); it contains a relatively flat region between z ¼ 0.5 nm and z ¼ 0.0 nm. As a result, the number of water molecules in this region (0.5 nm < z < 0.0 nm) is constant, with the value for the all-trans system being ~1 water molecule (per nm3) less than for the all-cis system. This remarkable result indicates that, due to the better packing of the trans lipids, the overall number of waters that access the all-trans lipid-water interface is reduced by ~1.13 waters per lipid compared with the all-cis bilayer. Hence, the total number of perturbed (interfacial) waters interacting with the trans lipids is also reduced. Consequently, the number and distribution of lipid-water and lipid-lipid hydrogen bonds in region III of the all-trans bilayer differs from that of the all-cis bilayer, as discussed below. Hydrogen bonding Evolution with time of the total number of hydrogen bonds (HBs) is reported in Fig. 5 for the all-trans and all-cis systems. This number, clearly larger for the trans, is a summation of inter- and intramolecular hydrogen bonds listed in Table 3. The most important intramolecular HBs in both bilayers involve the hydroxyl group as a donor and the carbonyl and phosphoryl oxygens as acceptors. The intermolecular HBs, on the other hand, are dominated by the amide group acting as donor to the carbonyl and hydroxyl groups of the neighboring lipids, in agreement with deuterium-NMR experimental findings (49). It is worth noting that the amide group participates in intermolecular HBs only, strengthening the interlipid interactions and thereby playing a crucial role in the tight packing of the PSM bilayers. The most significant difference between the trans and cis systems in terms of intramolecular HBs is in the hydroxylcarbonyl interaction, where there are six more HBs in the trans bilayer. In terms of intermolecular HBs, hydroxyl-amide and hydroxyl-hydroxyl interactions
Role of Sphingosine Double Bond
2961
with water molecules for 10 inter-PSM HBs and 7.5 intraPSM (of which six involve carbonyl groups). Finally, as previously noted, the amide groups intensify the intermolecular interactions in the trans bilayer by making 5.5 more HBs, but they also form eight more HBs with water molecules. These data clearly demonstrate that one of the major consequences of the presence of the trans double bond is to enable tighter lipid packing through the formation of larger number of HBs within the bilayer.
140
# PSM H-bonds
130 120 110 100 pure trans
90
Water order parameter
pure cis 80 0
50
100
150
200
250
time [ns]
FIGURE 5 Time evolution of the total number of hydrogen bonds in the trans and cis systems.
together are ~50% more numerous in the trans bilayer than in the cis. In total, ~10 more inter-PSM and ~7 more intraPSM (accounting for a total of additional 15%) HBs are formed in the all-trans bilayer than in the all-cis (see Fig. 6). The significance of these differences can be underscored by considering their potential contributions to the elastic strength of bilayers. With a conservative estimate of 2 kcal/mol energetic contribution from a single hydrogen bond (50), the extra ~10 interlipid hydrogen bonds would make the all-trans bilayer stiffer by ~20 kcal/mol than the all-cis, assuming that all the HB energy is translated into mechanical energy. The stronger interlipid interactions in the all-trans system (than in the all-cis) is accompanied by the loss of lipid-water HBs (dark-colored curve in the inset of Fig. 6). This is partially caused by the elimination of ~140 water molecules from the interfacial region of the trans bilayer (as calculated from Fig. 4), but also due to the exchange of lipid-water HBs for lipid HBs (light gray curve in the inset of Fig. 6). As a consequence, phosphoryl oxygens in the trans system become less hydrated and lose at least 12 HBs with water; carbonyl oxygens lose 16 HBs with water but gain six more intra-PSM HBs; hydroxyl groups exchange 12 HBs TABLE 3 Number of the most important intra- and interlipid hydrogen bonds, averaged over the last 50 ns of simulation for the pure trans and pure cis systems Nintra
Ninter
Donor
Acceptor
trans
cis
trans
cis
O–H N–H N–H O–H O–H N–H Total
O¼C O¼C O–H O–H O1/2-P O1/2-P
18.9 — — — 9.2 0.4 28.5
12.9 — — — 7.7 1.0 21.6
28.8 34.6 12.4 15.9 1.5 1.0 94.2
30.9 33.4 8.8 10.4 0.9 0.2 84.6
Hydrogen bonding was defined by an acceptor-hydrogen distance cutoff of ˚ and acceptor-hydrogen-donor angle cutoff of 90 (40,42). 2.5 A
Structural properties of the interfacial water may also be characterized by means of orientational order parameters (OPs). Orientational OPs can be defined in a general way using the so-called Wigner rotational D-matrices (see the Supporting Material). Rank 1 orientational OP (OP1) is equal to hcosbi, where b is the angle between the Ow / H vector and the bilayer normal (z axis). Plot of OP1 as a function of the Ow’s distance from the average z position of the C2 atom (open symbols in Fig. 7) shows practically no difference in the orientation of the interfacial waters in the all-trans and all-cis. As water approaches the bilayer, it tries to align its dipole toward the phosphoryl oxygens of the lipid headgroup (51). Past the phosphate group, this water orientation is maintained by donating hydrogen to the C¼O and OH groups of the lipids (see Fig. 6). Further into the bilayer (z >0.1 nm), water molecules start to accept HBs from OH and NH groups. The latter is located closer to the hydrophobic core and therefore contributes strongly to the water molecule’s dipole reorientation. As a consequence, the average dipole points toward the bilayer core for z < 0.1 nm, and away from it for z > 0.1 nm (see distributions in Fig. 4). Rank 2 OPs are defined similarly to rank one OPs and yield nonzero values only for 0
S ¼
1 2 3cos b 1 : 2
This OP (OP2) is particularly important because it is directly related to the quadrupole splitting of a deuterium NMR spectrum (52). Our calculations show that, although the PN vector population distribution is slightly different in the two systems (data not shown), its average orientation is about the same: ~81 for trans and ~79 for cis. Because OP2 in the two systems is very similar (see Fig. 7), our findings are in agreement with Aman et al. (48) that a notable change in the water OP2 value would imply a change in the population of the PN vector orientation along the membrane normal. Overall, the analysis of the water orientational OPs indicate that, despite differences in the number and interaction of the interfacial water molecules, the average orientation, Biophysical Journal 99(9) 2957–2966
2962
Janosi and Gorfe )
(1)
)
(3 (4 (2) (6) (5) -H O1-P O2-P -H w -H w O ----O w -O w -H w -H w -H w -H w -H-- -H-N Ow Ow Ow Ow O
=C
---O
Δintra: 6.9 Δinter: 9.6 Δtotal: 16.5 Δwater: 33.0
Δ # H-bonds
Number of H-bonds with PSM
16
80
60
8
100
0
80
-8 1 2
3
4
5
6
7
8
trans: intra trans: inter 60 cis: intra cis: inter trans: water 40 cis: water
40
20
0
(
(
(
(5
(6
(7
(8
and hence polarization, remains unaffected by trans/cis isomerization. PSM chain order parameter The orientational OP of the hydrocarbon chains is best described by the rank 2 OP, Smol ¼
FIGURE 6 Number of hydrogen bonds in the all-trans and all-cis systems. (Inset) Number of lipid-water hydrogen bonds that were lost when switching from cis to trans, and the number of inter- and intra-PSM H-bonds, as well as their total, gained at the same time. Hydrogen bonding was defined by an ˚ and acceptoracceptor-hydrogen distance cutoff of 2.5 A hydrogen-donor angle cutoff of 90 (40,42).
20
=C O=C O-H O-H 1-P 2-P 1-P 2-P ---O H--- -H--- -H--- H---O H---O ---O ---O H H H N- 3) N 4) O )O )N )N )O )O 2)
(1
120
Number of H-bonds with water
100
1 3cos2 f 1 ; 2
(2)
calculated at each nonterminal carbon position Cn in terms of f, the angle between the long molecular axis r, defined ! along the vector Cn1 Cnþ1, and the normal to the bilayer.
According to Eq. 2, Smol will be large (high-ordering) for large hcos2 fi, i.e., small f, and small but positive (lowerordering) for small hcos2 fi, i.e., f less than but close to the magic-angle of 54.74 . Smol for the PA and SPH chains in the pure bilayer simulations are plotted in Fig. 8, parts a and b, respectively. SPA mol is relatively constant throughout the chain in both cases, with lower values toward the tail ends due to their increased mobility at the interface between the two layers. Compared to the cis system, the higher packing of the trans lipids increase the overall ordering of the PA chain. Because all the PA tail bonds are saturated, SPA mol can be directly related
PA
0.8
0.2 trans: OP1 cis: OP1 trans: OP2 cis: OP2
Smol
0.1
0.6
a
0.4
total trans cis
0.8 0.6 Smol
-0.1
0.4 0.2
-0.2
0 0
-0.8 -0.4 0 0.4 z relative to avg. C2 glycerol position [nm]
0.8
FIGURE 7 First rank (open symbols) and second rank (solid symbols) water order parameters in the all-trans and all-cis systems. Biophysical Journal 99(9) 2957–2966
b
0.2
0
-0.3 -1.2
SPH pure bilayers
c
d
mixture 5
10 15 carbon atom
20
5
10 15 carbon atom
20
FIGURE 8 Molecular order parameters for the palmitoyl chain of the pure bilayers (a) and mixture (c) and for the sphingosine chain of the pure bilayers (b) and mixture (d). For the mixed bilayer, the average Smol calculated over the entire bilayer is also plotted.
Role of Sphingosine Double Bond
2963
PA to the deuterium order parameter SPA CD x 0:5 Smol (the average orientation of CH bonds is normal to the molec315 ¼ 0.249, calcuular axis) (52). Thus, the average hSPA CD i lated over carbon positions 3–15, matches very well the exp ¼ 0.221 5 0.047 reported experimental value of hSPA CD i at 321 K (53). In the case of SPH, the first two carbon atoms have low ordering in both systems, mainly because they have a large f-angle with respect to the bilayer normal inasmuch as they provide a direct link between the hydrophobic tails and the lipid headgroup (see Fig. 1). The double bond (carbons 4 and 5) and the following C6 are highly ordered in the alltrans system. In the all-cis system, however, C4 is ordered whereas C5 and C6 are highly disordered. Finally, the rest of the chain (C7–C15) exhibits similar ordering to the PA tails. Overall, the trans SPH tails are generally more ordered than the cis SPH tails due to the tighter packing of lipids. Because hcos2fi is directly related to the f-angle probability distribution p(f) (see Fig. 9) through the relation
2 cos f ¼
Z
fmax 0
cos2 fpðfÞdf;
the behavior of Smol around the double bond can be explained as follows. For trans double bond, the molecular axes r on positions 4 and 5 are expected to be almost parallel to each other, and therefore yield very similar p(f) and hence very similar values of S4mol and S5mol . The orientation of the molecular axis on these bonds influences the nearest neighbor, yielding S6mol values that are in-between the high ordering of the double bond and the average ordering of the rest of the chain. The cis double bond allows the same molecular orientation on position 4 due to the similar interfacial interactions. However, its different conformation enforces very different p(f) on positions 5 and 6. Because hf5i ~46 is close to the magic-angle and p(f5) is almost
hSmol itrans x0:52 > hSmol imixed x0:47 > hSmol icis x0:42; are larger for the better-packed bilayers: hAL itrans < hAL imixed < hAL icis : Such a relationship between chain order and area per lipid was found generally valid for previously well-studied lipids (47,54). Lipid dynamics We studied the dynamics of the lipids in terms of their lateral diffusion and rotational motions. The lateral diffusion coefficients were calculated from the slopes of the b
a C3-C5: 〈φ4〉 = 23.3°
0.03
C4=C5: 〈ψ〉 = 41.8°
C4=C5: 〈ψ〉 = 27.6°
C4-C6: 〈φ5〉 = 23.3°
C4-C6: 〈φ5〉 = 46.4°
C5-C7: 〈φ6〉 = 28.6°
0.02
Ci-1-Ci+1: 〈φi〉 = 33.3° C16-C18: 〈φ17〉 = 44.9°
C5-C7: 〈φ6〉 = 46.9°
0.02
Ci-1-Ci+1: 〈φi〉 = 33.2° C16-C18: 〈φ17〉 = 41.8° 0.01
0.01
0 0
C3-C5: 〈φ4〉 = 26.1°
0.03
p(φ), p(ψ)
p(φ), p(ψ)
symmetric, S5mol has a very low value. The change in the orientation of the chain due to the cis double-bond kink is almost always conserved by the next dihedral angle, maintaining a small value for S6mol (33,47,54). Fig. 8, c and d, shows Smol calculated for PA and SPH chains in the mixed trans/cis simulation, respectively. Clearly, the Smol values around the SPH double bond extracted separately for each of the two conformations are very similar to Smol from the pure bilayer simulations. Nonetheless, the total Smol computed for the entire bilayer, with trans/cis ratio of 80:20, is dictated by the dominating conformation of the system (see dark curve around the double bond in Fig. 8 d). Because the saturated parts of the chains are less rigid, their profiles are approximately the same for both conformations. In both pure bilayers and the mixture, the coupling between the tails maintains approximately the same average value for Ssat mol ; calculated over positions 7–15 for SPH and 4–12 for PA. Note that the actual positioning of the PA chain (with respect to the bilayer surface) relative to the SPH chain is shifted by 2–3 positions in the trans conformation and by 3–4 positions in the cis (see Fig. 1). These averages,
30
60 90 120 angle with z-axis [deg]
150
180
0 0
30
60 90 120 angle with z-axis [deg]
150
180
FIGURE 9 Angular distributions of f (see text), and of the angle between the bilayer normal and the double bond, j, for the sphingosine chain in the (a) all-trans and (b) all-cis systems. For clarity, f1 and f2 were omitted and an average angle distribution of the tail where Smol is relatively constant (see Fig. 8) was calculated for carbons 7–11. Biophysical Journal 99(9) 2957–2966
2964
mean-square displacement curves (see Fig. S1 a) at long times: 1 ½~ rðtÞ2 : (3) DL ¼ lim t/N 4t The obtained values are Dtrans ¼ 0.5 107 cm2/s and Dcis ¼ 1.6 107 cm2/s for all-trans and all-cis bilayers, respectively. The former is in good agreement with the experimental value of ~0.6 107 cm2/s obtained by the pulsed field gradient NMR technique (55). Clearly, the alltrans bilayer is not only more compact and ordered, but it is also less fluid than the all-cis. The autocorrelation functions calculated for an interfacial PA vector between the two chains (uniting CSPH 3 and C1 ) and for a vector along the headgroup (between P and N atoms), shown in Fig. S1 b, indicate that the trans configuration is slower in changing its orientation than the cis. Because of the intermolecular hydrogen bonds in the interfacial region, PA the decay of the interfacial reorientation ðCSPH 3 /C1 Þ autocorrelation function is extremely slow. This slow rotational motion is partially transmitted to the headgroup (P / N) due to the intramolecular contacts between the oxygens of the phosphate group and the hydroxyl/amide groups. CONCLUSIONS The conformation of the sphingosine double bond has been reported to play an important role in cholesterol binding to sphingomyelins (22), and to have a strong effect on the structure and hydrogen bonding (23), as well as elastic properties and packing behavior of ceramides (21). In this work we have performed an extensive computational study on the influence of the double-bond geometry on the structural properties of the palmitoyl sphingomyelin lipid bilayer. We found that the sphingosine base double-bond conformation has little impact on the orientation and dynamics of the solvent molecules because the ordering of the interfacial water is mainly regulated by the phosphate group. On the other hand, we observed tighter packing of the trans lipids, with the average area per lipid being ~10% larger for the all-cis system. This behavior is induced by the larger number of both intra- and intermolecular hydrogen bonds that the trans lipids can form. Additionally, the positioning of the trans phosphate group slightly farther from the hydrophobic core allows the lipids to move closer to each other. Hence, less water enters the interfacial region and hydrogen bonds with the excluded solvent molecules are exchanged for intermolecular hydrogen bonds. The tight hydrogen-bond network in the all-trans bilayer keeps the lipids closely packed, thus reducing the average area per lipid and increasing the ordering and length of the tails, as well as the overall thickness of the bilayer. As a consequence of the strong hydrogen-bond network the all-trans bilayer is less fluid than the all-cis. Finally, the high order around the Biophysical Journal 99(9) 2957–2966
Janosi and Gorfe
trans double bond and the disorder around the cis double bond found in the pure bilayers are maintained in the mixed bilayers as well: the ordering of the saturated part of the sphingosine base and of the palmitoyl chain are highest for all-trans, lowest for all-cis, and intermediate for the 4:1 trans/cis mixture. This emphasizes the importance of the double-bond conformation on the overall bilayer behavior with a direct impact on the average thermodynamic properties of bilayers. An important consequence of these effects of the trans double bond is the resultant changes in bilayer mechanical properties. The isothermal area compressibility modulus is much larger in the all-trans than in the all-cis system, suggesting that stretching bilayers of the naturally occurring trans PSM requires roughly twice the energy needed to stretch the corresponding cis bilayer. Furthermore, using the relationship between the bending rigidity, compressibility modulus, and hydrophobic thickness of the bilayer, KC ¼ KA dC2 2 C2 =24 (56), we can estimate from our data that KCtrans x2:7 KCcis : This shows that the trans double bond confers almost threetimes’ more resistance to bending of the bilayer than the cis double bond. These mechanical properties may explain the PSM-enrichment of functional lipid rafts in cellular membranes. It is worth noting that while most membrane lipids that are associated with phase separation and lipid rafts have intrinsic curvature, PSM does not. Instead, it appears that its involvement in phase separation arises from a fine balance of intra- and intermolecular hydrogen bonds that are tightly regulated by the trans double bond. SUPPORTING MATERIAL Description of orientational order parameters and a figure of mean square displacement and rotational correlation figure are available at http://www. biophysj.org/biophysj/supplemental/S0006-3495(10)01167-7. We gratefully acknowledge the Medical School at the University of Texas Health Science Center in Houston for financial support and the Texas Advanced Computing Center for computational resources.
REFERENCES 1. Karlsson, K. A. 1970. Sphingolipid long chain bases. Lipids. 5: 878–891. 2. Pruett, S. T., A. Bushnev, ., A. H. Merrill, Jr. 2008. Biodiversity of sphingoid bases (‘‘sphingosines’’) and related amino alcohols. J. Lipid Res. 49:1621–1639. 3. Ramstedt, B., and J. P. Slotte. 2006. Sphingolipids and the formation of sterol-enriched ordered membrane domains. Biochim. Biophys. Acta Biomembr. 1758:1945–1956.
Role of Sphingosine Double Bond 4. Samsonov, A. V., I. Mihalyov, and F. S. Cohen. 2001. Characterization of cholesterol-sphingomyelin domains and their dynamics in bilayer membranes. Biophys. J. 81:1486–1500. 5. Sonnino, S., A. Prinetti, ., G. Tettamanti. 2006. Dynamic and structural properties of sphingolipids as driving forces for the formation of membrane domains. Chem. Rev. 106:2111–2125. 6. Wang, T. Y., and J. R. Silvius. 2003. Sphingolipid partitioning into ordered domains in cholesterol-free and cholesterol-containing lipid bilayers. Biophys. J. 84:367–378. 7. Quinn, P. J., and C. Wolf. 2009. The liquid-ordered phase in membranes. Biochim. Biophys. Acta Biomembr. 1788:33–46. 8. Lingwood, D., and K. Simons. 2010. Lipid rafts as a membrane-organizing principle. Science. 327:46–50. 9. Mayor, S., and M. Rao. 2004. Rafts: scale-dependent, active lipid organization at the cell surface. Traffic. 5:231–240. 10. Brown, D. A. 2006. Lipid rafts, detergent-resistant membranes, and raft targeting signals. Physiology (Bethesda). 21:430–439. 11. Patra, S. K. 2008. Dissecting lipid raft facilitated cell signaling pathways in cancer. Biochim. Biophys. Acta-Rev. Cancer. 1785:182–206. 12. Simons, K., and D. Toomre. 2000. Lipid rafts and signal transduction. Nat. Rev. Mol. Cell Biol. 1:31–39. 13. Brown, D. A., and E. London. 1998. Functions of lipid rafts in biological membranes. Annu. Rev. Cell Dev. Biol. 14:111–136. 14. Obeid, L. M., C. M. Linardic, ., Y. A. Hannun. 1993. Programmed cell death induced by ceramide. Science. 259:1769–1771. 15. Ogretmen, B., and Y. A. Hannun. 2004. Biologically active sphingolipids in cancer pathogenesis and treatment. Nat. Rev. Cancer. 4: 604–616. 16. Hannun, Y. A., and L. M. Obeid. 2008. Principles of bioactive lipid signaling: lessons from sphingolipids. Nat. Rev. Mol. Cell Biol. 9:139–150. 17. Ohanian, J., and V. Ohanian. 2001. Sphingolipids in mammalian cell signaling. Cell. Mol. Life Sci. 58:2053–2068. 18. Hla, T. 2004. Physiological and pathological actions of sphingosine 1-phosphate. Semin. Cell Dev. Biol. 15:513–520. 19. Chalfant, C. E., and S. Spiegel. 2005. Sphingosine 1-phosphate and ceramide 1-phosphate: expanding roles in cell signaling. J. Cell Sci. 118:4605–4612. 20. Spiegel, S., and S. Milstien. 2003. Sphingosine-1-phosphate: an enigmatic signaling lipid. Nat. Rev. Mol. Cell Biol. 4:397–407. 21. Brockman, H. L., M. M. Momsen, ., R. Bittman. 2004. The 4,5double bond of ceramide regulates its dipole potential, elastic properties, and packing behavior. Biophys. J. 87:1722–1731. 22. Subbaiah, P. V., D. Sircar, ., R. Bittman. 2009. Effect of double bond geometry in sphingosine base on the antioxidant function of sphingomyelin. Arch. Biochem. Biophys. 481:72–79. 23. Phillips, S. C., G. Triola, ., M. C. Yappert. 2009. cis- versus trans-ceramides: effects of the double bond on conformation and H-bonding interactions. J. Phys. Chem. B. 113:15249–15255. 24. Ramstedt, B., and J. P. Slotte. 2002. Membrane properties of sphingomyelins. FEBS Lett. 531:33–37. 25. Bruzik, K. S., B. Sobon, and G. M. Salamonczyk. 1990. Nuclear magnetic resonance study of sphingomyelin bilayers. Biochemistry. 29:4017–4021. 26. Sripada, P. K., P. R. Maulik, ., G. G. Shipley. 1987. Partial synthesis and properties of a series of n-acyl sphingomyelins. J. Lipid Res. 28:710–718. 27. Koynova, R., and M. Caffrey. 1995. Phases and phase-transitions of the sphingolipids. Biochim. Biophys. Acta Lipids Lipid Metab. 1255: 213–236. 28. McIntosh, T. J., S. A. Simon, ., C. H. Huang. 1992. Interbilayer interactions between sphingomyelin and sphingomyelin/cholesterol bilayers. Biochemistry. 31:2020–2024.
2965 29. Jendrasiak, G. L., and R. L. Smith. 2001. The effect of the choline head group on phospholipid hydration. Chem. Phys. Lipids. 113:55–66. 30. Talbott, C. M., I. Vorobyov, ., M. C. Yappert. 2000. Conformational studies of sphingolipids by NMR spectroscopy. II. Sphingomyelin. Biochim. Biophys. Acta Biomembr. 1467:326–337. 31. Ferguson-Yankey, S. R., D. Borchman, ., M. C. Yappert. 2000. Conformational studies of sphingolipids by NMR spectroscopy. I. Dihydrosphingon-myelin. Biochim. Biophys. Acta Biomembr. 1467: 307–325. 32. Li, L., X. P. Tang, ., M. C. Yappert. 2002. Conformational characterization of ceramides by nuclear magnetic resonance spectroscopy. Biophys. J. 82:2067–2080. 33. Seelig, J., and N. Waespe-Sarcevic. 1978. Molecular order in cis and trans unsaturated phospholipid bilayers. Biochemistry. 17:3310–3315. 34. Chiu, S. W., S. Vasudevan, ., H. L. Scott. 2003. Structure of sphingomyelin bilayers: a simulation study. Biophys. J. 85:3624–3635. 35. Hess, B., C. Kutzner, ., E. Lindahl. 2008. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 4:435–447. 36. van Gunsteren, W., S. Billeter, ., I. Tironi. 1996. Biomolecular Simulation: The GROMOS 96 Manual and User Guide. Vdf, Zu¨rich, Switzerland. 37. Maulik, P. R., D. Atkinson, and G. G. Shipley. 1986. X-ray scattering of vesicles of n-acyl sphingomyelins. Determination of bilayer thickness. Biophys. J. 50:1071–1077. 38. Maulik, P. R., and G. G. Shipley. 1996. N-palmitoyl sphingomyelin bilayers: structure and interactions with cholesterol and dipalmitoylphosphatidylcholine. Biochemistry. 35:8025–8034. 39. Nagle, J. F., and S. Tristram-Nagle. 2000. Structure of lipid bilayers. Biochim. Biophys. Acta. 1469:159–195. 40. Mombelli, E., R. Morris, ., F. Fraternali. 2003. Hydrogen-bonding propensities of sphingomyelin in solution and in a bilayer assembly: a molecular dynamics study. Biophys. J. 84:1507–1517. 41. Hyvonen, M. T., and P. T. Kovanen. 2003. Molecular dynamics simulation of sphingomyelin bilayer. J. Phys. Chem. B. 107:9102–9108. 42. Niemela¨, P., M. T. Hyvo¨nen, and I. Vattulainen. 2004. Structure and dynamics of sphingomyelin bilayer: insight gained through systematic comparison to phosphatidylcholine. Biophys. J. 87:2976–2989. 43. Pandit, S. A., E. Jakobsson, and H. L. Scott. 2004. Simulation of the early stages of nano-domain formation in mixed bilayers of sphingomyelin, cholesterol, and dioleylphosphatidylcholine. Biophys. J. 87:3312–3322. 44. Niemela¨, P. S., M. T. Hyvo¨nen, and I. Vattulainen. 2006. Influence of chain length and unsaturation on sphingomyelin bilayers. Biophys. J. 90:851–863. 45. Lindahl, E., and O. Edholm. 2000. Mesoscopic undulations and thickness fluctuations in lipid bilayers from molecular dynamics simulations. Biophys. J. 79:426–433. 46. Feller, S. E., and R. W. Pastor. 1999. Constant surface tension simulations of lipid bilayers: the sensitivity of surface areas and compressibilities. J. Chem. Phys. 111:1281–1287. 47. Poger, D., and A. E. Mark. 2010. On the validation of molecular dynamics simulations of saturated and cis-monounsaturated phosphatidylcholine lipid bilayers: a comparison with experiment. J. Chem. Theory Comput. 6:325–336. 48. Aman, K., E. Lindahl, ., P. O. Westlund. 2003. Structure and dynamics of interfacial water in an La phase lipid bilayer from molecular dynamics simulations. Biophys. J. 84:102–115. 49. Steinbauer, B., T. Mehnert, and K. Beyer. 2003. Hydration and lateral organization in phospholipid bilayers containing sphingomyelin: a 2HNMR study. Biophys. J. 85:1013–1024. 50. Vorobyov, I., M. C. Yappert, and D. B. DuPre. 2002. Energetic and topological analyses of cooperative sigma H- and p-H-bonding interactions. J. Phys. Chem. A. 106:10691–10699. Biophysical Journal 99(9) 2957–2966
2966 51. Hyvo¨nen, M. T., T. T. Rantala, and M. Ala-Korpela. 1997. Structure and dynamic properties of diunsaturated 1-palmitoyl-2-linoleoyl-snglycero-3-phosphatidylcholine lipid bilayer from molecular dynamics simulation. Biophys. J. 73:2907–2923. 52. Seelig, J. 1977. Deuterium magnetic resonance: theory and application to lipid membranes. Q. Rev. Biophys. 10:353–418. 53. Mehnert, T., K. Jacob, ., K. Beyer. 2006. Structure and lipid interaction of n-palmitoylsphingomyelin in bilayer membranes as revealed by 2 H-NMR spectroscopy. Biophys. J. 90:939–946.
Biophysical Journal 99(9) 2957–2966
Janosi and Gorfe 54. Ro´g, T., K. Murzyn, ., M. Pasenkiewicz-Gierula. 2004. Effects of phospholipid unsaturation on the bilayer nonpolar region: a molecular simulation study. J. Lipid Res. 45:326–336. 55. Filippov, A., G. Ora¨dd, and G. Lindblom. 2003. The effect of cholesterol on the lateral diffusion of phospholipids in oriented bilayers. Biophys. J. 84:3079–3086. 56. Rawicz, W., K. C. Olbrich, ., E. Evans. 2000. Effect of chain length and unsaturation on elasticity of lipid bilayers. Biophys. J. 79:328–339. 57. Humphrey, W., A. Dalke, and K. Schulten. 1996. VMD: visual molecular dynamics. J. Mol. Graph. 14:33–38, 27–28.
