Examination of Mixtures Containing Sphingomyelin and Cholesterol

Apr 11, 2017 - Examination of Mixtures Containing Sphingomyelin and Cholesterol by Molecular Dynamics Simulations. Eric Wang† and Jeffery B. Klaudaâ...
68 downloads 16 Views 3MB Size
Article pubs.acs.org/JPCB

Examination of Mixtures Containing Sphingomyelin and Cholesterol by Molecular Dynamics Simulations Eric Wang† and Jeffery B. Klauda*,†,‡ †

Department of Chemical and Biomolecular Engineering and ‡Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, United States S Supporting Information *

ABSTRACT: The all-atom CHARMM36 (C36) force field is used to simulate bilayers of pure palmitoylsphingomyelin (PSM) as well as binary mixtures of PSM and stearoylsphingomyelin (SSM) at various cholesterol concentrations (XC) and temperatures. C36 simulation data is in good agreement with experimental deuterium order parameters and previous computational results, providing evidence of the utility of the force field for potentially studying more complex membranes. The area compressibility modulus is shown to achieve a large value of 2.82 ± 0.08 N/m in cholesterol-rich membranes (XC = 0.50). Surface area per lipid (SA/lip), tilt angle, membrane thicknesses, and acyl chain ordering are shown to have strong dependencies on cholesterol concentration. Relaxation times also indicate cholesterol dependence and show a strong preference for rotational axial motion over wobbling motion. Radial distribution functions and lipid clustering indicate strong relationships between lateral ordering and hydrogen bonding, which is long lived in SM membranes. These interactions lead to strong self-association of cholesterol at high concentrations, causing shielding from further SM−cholesterol interactions. The importance of a ternary component on SM−SM hydrogen bonds is revealed in light of previous results and is consequential in the modeling of lipid rafts.

1. INTRODUCTION Sphingolipids are a class of membrane lipids that are pervasive in eukaryotes and carry out a wide variety of essential cell functions such as lipid homeostasis, metabolism, antigen presentation, and cell trafficking.1 Many lipids contain saturated tails that aid in the membrane’s ability to form tightly packed structures that help membranes endure mechanical stress.2 A particular subset of sphingolipids of great interest in recent years is sphingomyelin (SM). The discovered functions for SM are numerous and continually being revealed. SM is abundant in brain and nervous tissue, and the sphingomyelinases which break down SMs to ceramides are believed to be connected to various neurological disorders.3 The hydrolysis of SM leads to the formation of critical second messengers such as ceramides and its derivatives, affecting properties of membranes such as lipid phases, cholesterol concentration and esterification, and protein recognition.2 Most strikingly, SM and its interactions with cholesterol are known to have important roles in the formation of lipid rafts, which are ordered microdomains within the cell membrane that emerged as a challenge to the prevailing fluid mosaic model and carry profound biological significance.4 Since the initial proposal of the lipid raft hypothesis,5 extensive research has led to many possible explanations for raft function and formation. Rafts have been suggested to play roles in signal transduction, protein sorting, microbial infection,6 and nerve terminal function.7 The details of lipid raft formation are © 2017 American Chemical Society

currently unknown, but it has been found that the properties of rafts are similar to those of liquid-ordered phases and that rafts prefer SM over other lipids such as phosphatidylcholines (PC) because of the deeper intrusion of cholesterol into SM membranes.6 Analysis of in vivo rafts has proven to be quite difficult, so studies often use systems of SM and cholesterol as a model. The rapid advance of computational power has allowed for molecular dynamics (MD) simulations to emerge as a powerful method of studying membranes. An early study in 1993 spent 6 months performing a 100 ps simulation of a simple bilayer,8 whereas present-day simulations can run in the millisecond range using advanced supercomputing architecture.9 MD simulations can be limited by time scales (submicroseconds for most systems on standard computational resources) and potential inaccuracies in the force fields but are cost efficient and provide access to a level of detail that is difficult to impossible with experimental methods. Systems of SM and cholesterol have previously been studied using MD, the most common SM being palmitoylsphingomyelin (PSM) and stearoylsphingomyelin (SSM). The structure of these two lipids are very similar, containing saturated tails that differ only by two carbons. Previous simulation studies of PSM/ Received: February 24, 2017 Revised: April 5, 2017 Published: April 11, 2017 4833

DOI: 10.1021/acs.jpcb.7b01832 J. Phys. Chem. B 2017, 121, 4833−4844

Article

The Journal of Physical Chemistry B

Figure 1. (a) Chemical structure of PSM, SSM, and cholesterol. Terminal tail carbons on SM lipids and the C3 and C17 carbons on cholesterol are labeled. (b) Snapshot of SSM at 50 °C and XC = 0.20 at the end of the simulation. SSM lipids are shown in cyan with colored headgroups, cholesterol molecules are in red, and water is shown as purple dots. (c) Snapshot of SSM at 50 °C and XC = 0.50 at the end of simulation. Color scheme is the same as in (b).

Table 1. System Parameters and General Bilayer Propertiesa lipid

T (°C)

XC

no. of lipids

no. of waters

time (ns)

PSM

45

0.00 0.20 0.33 0.00 0.20 0.33 0.20 0.33 0.50

80 100 100 80 100 100 100 100 100

2400 3500 3500 2400 3500 3500 3500 3500 3500

300 400 400 200 600 400 600 500 400

60

SSM

a

50

SA/lip (Å2) 54.8 43.9 41.7 57.5 46.8 43.1 44.5 41.7 40.3

(0.2) (0.2) (0.2) (0.2) (0.2) (0.1) (0.1) (0.1) (0.1)

KA (N/m) 0.45 1.08 2.01 0.45 0.61 1.14 0.84 2.02 2.82

(0.04) (0.06) (0.10) (0.04) (0.04) (0.06) (0.05) (0.09) (0.08)

Standard errors are represented in parentheses. SA/lip is the average for all lipids (SM and cholesterol).

2. METHODS 2.1. System Setup and Simulation Protocol. Experimental and computational studies have produced results on SM/cholesterol mixtures at cholesterol concentrations of 0.00, 0.20, 0.33, and 0.50,6,13,16,17 so mixtures were built to mimic these systems for simple comparison. PSM systems were built using cholesterol concentrations of 0.00, 0.20, and 0.33 at both 45 and 60 °C, while SSM systems were built using cholesterol concentrations of 0.20, 0.33, and 0.50 only at 50 °C (Figure 1). Pure SM bilayers contained 40 lipids per leaflet with 30 water molecules per lipid, and SM/cholesterol mixtures contained 50 lipids per leaflet with 35 water molecules per lipid. Triplicate replicas of each membrane system were generated in rectangular boxes and tetragonal cells (X = Y ≠ Z) using CHARMM-GUI Membrane Builder and were equilibrated for 385 ps with a 1 fs time step using the Membrane Builder18−21 six-step equilibration process to slowly remove restraints. After this, production runs were ran for 200−600 ns, until every replica had at least 150 ns of stable equilibrated data, which was determined by examination of surface area per lipid (SA/lip) plots. The total production time for each system can be found in Table 1. Most systems assumed a symmetric distribution across leaflets. However, for systems with XC = 0.33, the lipids were asymmetrically distributed between leaflets such that one leaflet would have 16 cholesterol and 34 SM while the other would have 17 cholesterol and 33 SM. Concerns about mismatch changing the lateral pressure of each leaflet are mitigated by

cholesterol and SSM/cholesterol mixtures were run often run 0.09), only becoming significant for the initial introduction of cholesterol (p = 0.004) or for large disparities in concentration such as between XC = 0.20 and XC = 0.50 (p = 0.03). Additionally, greater cholesterol concentration causes a tendency for cluster composition to reflect membrane composition which can be seen as the difference (RC − RN) approaches 0 for both cholesterol and SM. The tendency for cholesterol to participate in lipid clusters in higher proportions and form segregated domains (Figure 1b and 1c, Figure S5) at high concentration possibly explains the previously mentioned increase in SM SA/lip since this causes SM to be excluded from interactions with cholesterol. Decreases in SM−cholesterol H bonding are in agreement with this trend (Table 5). In contrast to cholesterol concentration, temperature has a minimal effect and produces statistically identical values for all cluster properties (p > 0.20). 3.9. Relaxation Times. Details of SM intramolecular motion can be explored by the calculation of three fit exponentials to the correlation functions of carbon−hydrogen and cross-chain carbon−carbon vectors as previously described. The cross-chain vector is expected to limit local torsional motion and only represent axial motion, giving insight into the motion of C2−H, which is an unknown combination of both axial motion and lipid wobble.50 The equilibration times, fit coefficients, and relaxation times are reported in Table S7. However, some systems require considerable time to relax, and the corresponding correlation functions did not plateau within the equilibrated time range, often >300 ns (Figure S13). As a result, the resulting fit coefficients and relaxation times are extracted from incomplete data and sometimes exhibit wide variability between replicas. Additionally, one particular case, PSM at 45 °C and XC = 0.20, was equilibrated long enough to plateau, but one replica plateaued at a significantly different value, producing different fit coefficients. Despite these obstacles, the calculated relaxation times remain useful for providing a general estimate and illustrating trends over temperature and cholesterol concentration. The fast and intermediate relaxation times (τ1 and τ2) do not show a consistent trend with either cholesterol concentration or temperature. However, the slow relaxation time (τ3) which is associated with lipid wobble51 is of most interest. In general, τ3 shows the expected trend of increasing time with cholesterol concentration and decreasing time with temperature (Table 7). The strongest shift occurs when cholesterol is first added in which τ3 increases 2−3-fold, and this occurs for both crosschain and C2−H relaxation times. The C2−H slow relaxation time of pure PSM was previously calculated to be 94 ns,13 so this work’s τ3 values of 150 ns in the 45 °C and XC = 0.20 case is consistent. Discrepancies appear when comparing against this work’s pure PSM τ3 value of 56 to 94 ns from the previous simulation,13 which may be the result of shorter time scales in our current work when compared to the past simulations of 400 ns. It is interesting to note that the τ3 relaxation times for the cross-chain and C2−H vectors are quite similar, implying the preferred mode is likely axial rotation and not lipid wobble. This bias for rotation is consistent with the tendency of cholesterol and SM to remain upright in the membrane which restricts space. Cholesterol’s orientation is clear from the tilt angle (Table 3), and SM can be inferred from the relatively high ordering (Tables S2−4, Figure S14) as well as the high proportion of amide−carbonyl inter H bond relative to amide−

Table 7. Slow Relaxation Time (τ3) Calculated from ThirdOrder Exponential Fits to Reorientational Correlation Functions for Both the Cross-Chain (C4S−C2F) and the C2−H Vector in Each Systema decay

lipid

T (°C)

XC

τ3 (ns)

cross chain

PSM

45

0.00 0.20b 0.33b 0.00 0.20 0.33 0.20b 0.33 0.50b 0.00 0.20b 0.33b 0.00 0.20 0.33 0.20b 0.33 0.50b

56 (7) 117 (9) 306 (66) 24 (1) 63 (6) 62 (3) 132 (21) 101 (17) 194 (18) 56 (4) 150 (25) 255 (10) 30 (3) 58 (6) 65 (2) 131 (25) 88 (5) 183 (15)

60

C2−H

SSM

50

PSM

45

60

SSM

50

a

Standard errors are represented in parentheses. bThe reorientational correlation function did not plateau within the time range, and fits are based on extrapolations of the curve.

phosphate bonds (Table S5). In addition, SM has a tendency to form long-lasting H bonds,11 which itself is known to slow even rotational motion of SM compared to other lipids,52 so the absence of lipid wobble is reasonable.

4. DISCUSSION AND CONCLUSION Many previous simulation studies with SM and cholesterol used relatively short time scales due to computational cost.10−12 Zidar et al.10 performed 30 ns simulations of SSM with 0−40% cholesterol using a force field for SSM that was transferred from analogous functional groups in the CHARMM force field. The unoptimized force field used by Zidar et al.10 resulted in a bilayer that was too dense, and thus, the SCD was higher when compared to that published here with the updated C36 lipid force field for sphingomyelin. Although our work agrees with Zidar et al.10 in that there is a general decrease in SM− cholesterol H bonding with cholesterol concentration, the values in our work are around 30% lower. Ultimately, it is hard to fully compare our work with Zidar et al. because 30 ns is too short to reach lateral equilibration in the bilayer, especially in multicomponent membranes. Other studies using various force fields have also been performed on SM/cholesterol mixtures,11,12 and the work by Zhang et al.12 using the Niemelä et al.53 force field for SM was run long enough to make meaningful comparisons (100 ns). These past results for the EDP of SSM with 34% cholesterol lack the secondary peaks caused by cholesterol at ±10 Å that we observe in our simulation (Figure 3b). Moreover, the DHH is higher with the past work (48.6 Å) compared to this work (47.9 Å). The SCD of the sphingosine chain by Zhang et al.12 lacks the decrease in order for the double-bonded carbons (C4 and C5) in comparison with neighboring carbons, and the order is higher than our work. This is consistent with the SA/lip being lower for the past work (39.4 ± 0.4 Å2) in comparison to our work (41.7 ± 0.1 Å2). Overall, the C36 force field predicts a slightly 4841

DOI: 10.1021/acs.jpcb.7b01832 J. Phys. Chem. B 2017, 121, 4833−4844

Article

The Journal of Physical Chemistry B

the amide plane through MD simulations and ultimately discover what level of lipid diversity is needed for strong lipid raft modeling. The mixtures presented in this work are either pure bilayers of SM or binary mixtures of SM/cholesterol at varying concentrations and temperatures. Comparison of deuterium order parameters with experimental data showed good agreement, and consistency was demonstrated with past C36 simulations based on order parameters, H bonding, relaxation times, SA/lip, and KA. An inaccuracy in the force field is demonstrated based on order parameter data, which shows an early phase transition to a liquid-ordered state at XC = 0.20 and minimal change from XC = 0.33 to XC = 0.50. However, order parameters for the liquid-ordered state are in agreement with experiment. In general, this study is further evidence that the C36 force field is an overall accurate model of molecular interactions and can be used to study more complex structures. Unexpected breaks from trends at high cholesterol concentrations are explained by the self-association of cholesterol which lead to a shielding effect. H-bonding results (Table 5) provide computational support for the decrease in SM−SM H bonding with cholesterol concentration, which has been experimentally noted and has implications in the study of liquid-ordered domains and consequently lipid rafts.

less dense membrane when compared to past simulations with other force fields and in agreement with simulations of pure SM bilayers.13 Since MD simulations with the C36 force field on pure and SM mixtures with cholesterol agree with experiment better than published work with other force fields, the C36 force field is well tuned for looking at mixtures with SM lipids. Experimental studies of SM/cholesterol binary mixtures by Raman spectroscopy54 and computational studies of ternary mixtures involving SM and cholesterol44 using the C36 force field have shed light on possible factors leading to the formation of liquid-ordered phases. In particular, they highlight the importance of hydrogen-bonding networks which are particularly high in SM/cholesterol mixtures. The formation of these liquid-ordered phases, occurring near a cholesterol concentration of 30%, could be similar to the formation of lipid rafts at the nanoscale.10 Notably, not all SM are relevant as oleoylsphingomyelin has been eliminated as a factor in raft formation because of its inability to form segregated domains with cholesterol, a result of unsaturation in the tail.55 Raman spectroscopy on SSM/cholesterol mixtures has suggested a decrease in SM−SM H bonding with cholesterol concentration which is apparently contradicted by the mechanism of amide plane tilting caused by SM−cholesterol interactions proposed by Sodt et al.44 in a ternary system containing SM, cholesterol, and PC. From these studies it appears that the presence of a ternary PC component has a profound impact on SM−SM H bonding. The results of the work published here are in agreement with a past binary SM/cholesterol simulation in the OPLS system56 and Raman spectroscopy as a clear decrease in SM−SM inter H bonding is demonstrated, which is especially prominent in the amide−carbonyl pair. This drop in H bonding can be explained by the tendency of cholesterol to disrupt groups of associating SM which then spread out while maintaining order, also resulting in an increase in SA/lip. Thus, H bonding decreases because the distance between Hbonding groups decreases. This is more powerful than the modest increases which occur in ternary SM/PC/cholesterol mixtures,14 suggesting that amide plane tilting is not always significant. Nevertheless, while binary mixtures remain suitable models of lipid rafts based on the high proportion of SM and cholesterol, future studies should take lipid diversity into account because the full significance of ternary components on lateral organization has yet to be explored. Despite the difficulties associated with applying experimental methods of studying binary systems to ternary systems and beyond, efforts are being made to take greater lipid diversity into account. The use of unsaturated PCs is relatively common because cholesterol promotes phase separation between the low main transition temperature (TM) PC and the high TM SM.57 Results from experimental studies of ternary mixtures found the presence of SM/cholesterol complexes with partitioning into the PC phase58 which is similar to the existence of such complexes in binary systems.59 Additionally, it was found that the use of different PC lipids such as 1-palmitoyl-2-oleoyl-snglycero-3-phosphocholine versus 1-stearoyl-2-oleoyl-sn-glycero3-phosphocholine yielded similar results, and both could be used to model rafts.60 Even greater diversity has been probed with the inclusion of ganglioside in SM/PC/cholesterol membranes to observe ganglioside-rich microdomains.61 Based on the importance of ternary components on SM−SM H bonding revealed in past work and differences with our observations for binary lipids, it would be interesting to probe how an additional component such as ganglioside would affect



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcb.7b01832. Figures of SA/lip vs time, deuterium order parameters, component EDPs, RDFs, XY-plane snapshots, reorientational correlation functions, and lipid simulation snapshots; tables of interpolated SA/lip values, order parameters, hydrogen bonding, clustering fractions, and fitted coefficients of reorientational correlation functions (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Jeffery B. Klauda: 0000-0001-8725-1870 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This research was supported by the NSF grant MCB-1149187. The high-performance computational resources used for this research are Deepthought and Deepthought2, which are maintained by the Division of Information Technology at the University of Maryland. Production runs used the Extreme Science and Engineering Discovery Environment (XSEDE) allocation on Stampede (MCB-100139). We especially thank Xiaohong Zhuang for sharing data analysis scripts and Hojin Kang for running simulations of SM/cholesterol mixtures.



REFERENCES

(1) Hla, T.; Dannenberg, A. J. Sphingolipid Signaling in Metabolic Disorders. Cell Metab. 2012, 16 (4), 420−434. (2) van Meer, G.; Voelker, D. R.; Feigenson, G. W. Membrane Lipids: Where They Are and How They Behave. Nat. Rev. Mol. Cell Biol. 2008, 9 (2), 112−124.

4842

DOI: 10.1021/acs.jpcb.7b01832 J. Phys. Chem. B 2017, 121, 4833−4844

Article

The Journal of Physical Chemistry B

Biological Membrane Simulations. J. Comput. Chem. 2014, 35 (27), 1997−2004. (21) Jo, S.; Lim, J. B.; Klauda, J. B.; Im, W. CHARMM-GUI Membrane Builder for Mixed Bilayers and Its Application to Yeast Membranes. Biophys. J. 2009, 97 (1), 50−58. (22) Park, S.; Beaven, A. H.; Klauda, J. B.; Im, W. How Tolerant are Membrane Simulations with Mismatch in Area per Lipid between Leaflets? J. Chem. Theory Comput. 2015, 11 (7), 3466−3477. (23) Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kale, L.; Schulten, K. Scalable Molecular Dynamics with NAMD. J. Comput. Chem. 2005, 26 (16), 1781−1802. (24) Feller, S. E.; Zhang, Y. H.; Pastor, R. W.; Brooks, B. R. Constant-Pressure Molecular-Dynamics Simulation - the Langevin Piston Method. J. Chem. Phys. 1995, 103 (11), 4613−4621. (25) Martyna, G. J.; Tobias, D. J.; Klein, M. L. Constant Pressure Molecular Dynamics Algorithms. J. Chem. Phys. 1994, 101 (5), 4177− 4189. (26) Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald - an NLog(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98 (12), 10089−10092. (27) Steinbach, P. J.; Brooks, B. R. New Spherical-cutoff Methods for Long-range Forces in Macromolecular Simulation. J. Comput. Chem. 1994, 15 (7), 667−683. (28) Lim, J. B.; Rogaski, B.; Klauda, J. B. Update of the Cholesterol Force Field Parameters in CHARMM. J. Phys. Chem. B 2012, 116 (1), 203−210. (29) Durell, S. R.; Brooks, B. R.; Bennaim, A. Solvent-Induced Forces between Two Hydrophilic Groups. J. Phys. Chem. 1994, 98 (8), 2198− 2202. (30) Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79 (2), 926−935. (31) Andersen, H. C. Rattle - a Velocity Version of the Shake Algorithm for Molecular-Dynamics Calculations. J. Comput. Phys. 1983, 52 (1), 24−34. (32) Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graphics 1996, 14 (1), 33−38. (33) Racine, J. Gnuplot 4.0: A Portable Interactive Plotting Utility. Journal of Applied Econometrics 2006, 21 (1), 133−141. (34) Brooks, B. R.; Brooks, C. L.; Mackerell, A. D.; Nilsson, L.; Petrella, R. J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, S.; et al. CHARMM: The Biomolecular Simulation Program. J. Comput. Chem. 2009, 30 (10), 1545−1614. (35) Brooks, B. R.; Bruccoleri, R. E.; Olafson, B. D. States, D. J.; Swaminathan, S.; Karplus, M., CHARMM - a Program for Macromolecular Energy, Minimization, and Dynamics Calculations. J. Comput. Chem. 1983, 4 (2), 187−217. (36) Dufourc, E. J.; Parish, E. J.; Chitrakorn, S.; Smith, I. C. P. Structural and Dynamical Detials of Cholesterol Lipid Interaction as Revealed by Deuterium NMR. Biochemistry 1984, 23 (25), 6062− 6071. (37) Jacobsen, J. P.; Schaumburg, K. Determination of Deuterium Quadrupole Coupling-constants in Methyl-groups by High-resolution NMR-Spectroscopy. J. Magn. Reson. 1977, 28 (2), 191−201. (38) Barber, C. B.; Dobkin, D. P.; Huhdanpaa, H. The Quickhull Algorithm for Convex Hulls. Acm Transactions on Mathematical Software 1996, 22 (4), 469−483. (39) Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825−2830. (40) Ester, M.; Kriegel, H.-P.; Sander, J.; Xiaowei, Xu A Densitybased Algorithm for Discovering Clusters in Large Spatial Databases with Noise. Proceedings of 2nd International Conference on Knowledge Discovery and Data Mining; AAAI, 1996; Vol. 96 (34), pp 226−231. (41) Hogberg, C. J.; Lyubartsev, A. P. A Molecular Dynamics Investigation of the Influence of Hydration and Temperature on

(3) Horres, C. R.; Hannun, Y. A. The Roles of Neutral Sphingomyelinases in Neurological Pathologies. Neurochem. Res. 2012, 37 (6), 1137−1149. (4) Simons, K.; Toomre, D. Lipid Rafts and Signal Transduction. Nat. Rev. Mol. Cell Biol. 2000, 1 (1), 31−39. (5) Simons, K.; Ikonen, E. Functional Rafts in Cell Membranes. Nature 1997, 387 (6633), 569−572. (6) Yasuda, T.; Kinoshita, M.; Murata, M.; Matsumori, N. Detailed Comparison of Deuterium Quadrupole Profiles between Sphingomyelin and Phosphatidylcholine Bilayers. Biophys. J. 2014, 106 (3), 631−638. (7) Saghy, E.; Szoke, E.; Payrits, M.; Helyes, Z.; Borzsei, R.; Erostyak, J.; Janosi, T. Z.; Setalo, G., Jr.; Szolcsanyi, J. Evidence for the Role of Lipid Rafts and Sphingomyelin in Ca2+-Gating of Transient Receptor Potential Channels in Trigeminal Sensory Neurons and Peripheral Nerve Terminals. Pharmacol. Res. 2015, 100, 101−116. (8) Venable, R. M.; Zhang, Y. H.; Hardy, B. J.; Pastor, R. W. Molecular-Dynamics Simulations of a Lipid Bilayer and of Hexadecane - an Investigation of Membrane Fluidity. Science 1993, 262 (5131), 223−226. (9) Shaw, D. E.; Dror, R. O.; Salmon, J. K.; Grossman, J. P.; Mackenzie, K. M.; Bank, J. A.; Young, C.; Deneroff, M. M.; Batson, B.; Bowers, K. J. Millisecond-Scale Molecular Dynamics Simulations on Anton. Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis, Portland, OR, Nov 14−20, 2009; ACM: New York, 2009,; p 11. (10) Zidar, J.; Merzel, F.; Hodoscek, M.; Rebolj, K.; Sepcic, K.; Macek, P.; Janezic, D. Liquid-Ordered Phase Formation in Cholesterol/Sphingomyelin Bilayers: All-Atom Molecular Dynamics Simulations. J. Phys. Chem. B 2009, 113 (48), 15795−15802. (11) Mombelli, E.; Morris, R.; Taylor, W.; Fraternali, F. HydrogenBonding Propensities of Sphingomyelin in Solution and in a Bilayer Assembly: A Molecular Dynamics Study. Biophys. J. 2003, 84 (3), 1507−1517. (12) Zhang, Z.; Bhide, S. Y.; Berkowitz, M. L. Molecular Dynamics Simulations of Bilayers Containing Mixtures of Sphingomyelin with Cholesterol and Phosphatidylcholine with Cholesterol. J. Phys. Chem. B 2007, 111 (44), 12888−12897. (13) Venable, R. M.; Sodt, A. J.; Rogaski, B.; Rui, H.; Hatcher, E.; MacKerell, A. D.; Pastor, R. W.; Klauda, J. B. CHARMM All-Atom Additive Force Field for Sphingomyelin: Elucidation of Hydrogen Bonding and of Positive Curvature. Biophys. J. 2014, 107 (1), 134− 145. (14) Bera, I.; Klauda, J. Molecular Simulations of Mixed Lipid Bilayers with Sphingomyelin, Glycerophospholipids and Cholesterol. Submitted for publication, 2017. (15) Klauda, J. B.; Venable, R. M.; Freites, J. A.; O’Connor, J. W.; Tobias, D. J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A. D.; Pastor, R. W. Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types. J. Phys. Chem. B 2010, 114 (23), 7830−7843. (16) Bartels, T.; Lankalapalli, R. S.; Bittman, R.; Beyer, K.; Brown, M. F. Raftlike Mixtures of Sphingomyelin and Cholesterol Investigated by Solid-State H-2 NMR Spectroscopy. J. Am. Chem. Soc. 2008, 130 (44), 14521−14532. (17) Matsumori, N.; Yasuda, T.; Okazaki, H.; Suzuki, T.; Yamaguchi, T.; Tsuchikawa, H.; Doi, M.; Oishi, T.; Murata, M. Comprehensive Molecular Motion Capture for Sphingomyelin by Site-Specific Deuterium Labeling. Biochemistry 2012, 51 (42), 8363−8370. (18) Jo, S.; Kim, T.; Iyer, V. G.; Im, W. CHARMM-GUI: a Webbased Graphical User Interface for CHARMM. J. Comput. Chem. 2008, 29 (11), 1859−1865. (19) Jo, S.; Kim, T.; Im, W. Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations. PLoS One 2007, 2 (9), e880. (20) Wu, E. L.; Cheng, X.; Jo, S.; Rui, H.; Song, K. C.; DavilaContreras, E. M.; Qi, Y. F.; Lee, J. M.; Monje-Galvan, V.; Venable, R. M.; et al. CHARMM-GUI Membrane Builder Toward Realistic 4843

DOI: 10.1021/acs.jpcb.7b01832 J. Phys. Chem. B 2017, 121, 4833−4844

Article

The Journal of Physical Chemistry B Structural and Dynamical Properties of a Dimyristoylphosphatidylcholine Bilayer. J. Phys. Chem. B 2006, 110 (29), 14326−14336. (42) Nagle, J. F.; Tristram-Nagle, S. Structure of Lipid Bilayers. Biochim. Biophys. Acta, Rev. Biomembr. 2000, 1469 (3), 159−195. (43) Zhuang, X. H.; Makover, J. R.; Im, W.; Klauda, J. B. A Systematic Molecular Dynamics Simulation Study of Temperature Dependent Bilayer Structural Properties. Biochim. Biophys. Acta, Biomembr. 2014, 1838 (10), 2520−2529. (44) Sodt, A. J.; Pastor, R. W.; Lyman, E. Hexagonal Substructure and Hydrogen Bonding in Liquid-Ordered Phases Containing Palmitoyl Sphingomyelin. Biophys. J. 2015, 109 (5), 948−955. (45) Boughter, C. T.; Monje-Galvan, V.; Im, W.; Klauda, J. B. Influence of Cholesterol on Phospholipid Bilayer Structure and Dynamics. J. Phys. Chem. B 2016, 120 (45), 11761−11772. (46) MacDermaid, C. M.; Kashyap, H. K.; DeVane, R. H.; Shinoda, W.; Klauda, J. B.; Klein, M. L.; Fiorin, G. Molecular Dynamics Simulations of Cholesterol-rich Membranes using a Coarse-grained Force Field for Cyclic Alkanes. J. Chem. Phys. 2015, 143 (24), 243144. (47) O’Connor, J. W.; Klauda, J. B. Lipid Membranes with a Majority of Cholesterol: Applications to the Ocular Lens and Aquaporin 0. J. Phys. Chem. B 2011, 115 (20), 6455−6464. (48) Yasuda, T.; Al Sazzad, M. A.; Jantti, N. Z.; Pentikainen, O. T.; Slotte, J. P. The Influence of Hydrogen Bonding on Sphingomyelin/ Colipid Interactions in Bilayer Membranes. Biophys. J. 2016, 110 (2), 431−440. (49) Zhuang, X.; Ou, A.; Klauda, J. Simulation of Simple Linoleic Acid-containing Lipid Membranes to Models of Soybean Plasma Membranes. Biophys. J. 2017, 112, 82a. (50) Klauda, J. B.; Roberts, M. F.; Redfield, A. G.; Brooks, B. R.; Pastor, R. W. Rotation of Lipids in Membranes: Molecular Dynamics Simulation, P-31 Spin-lattice Relaxation, and Rigid-body Dynamics. Biophys. J. 2008, 94 (8), 3074−3083. (51) Klauda, J. B.; Eldho, N. V.; Gawrisch, K.; Brooks, B. R.; Pastor, R. W. Collective and Noncollective Models of NMR Relaxation in Lipid Vesicles and Multilayers. J. Phys. Chem. B 2008, 112 (19), 5924− 5929. (52) Niemela, P.; Hyvonen, M. T.; Vattulainen, I. Structure and Dynamics of Sphingomyelin Bilayer: Insight Gained through Systematic Comparison to Phosphatidylcholine. Biophys. J. 2004, 87 (5), 2976−2989. (53) Niemela, P. S.; Hyvonen, M. T.; Vattulainen, I. Influence of Chain Length and Unsaturation on Sphingomyelin Bilayers. Biophys. J. 2006, 90 (3), 851−863. (54) Shirota, K.; Yagi, K.; Inaba, T.; Li, P. C.; Murata, M.; Sugita, Y.; Kobayashi, T. Detection of Sphingomyelin Clusters by Raman Spectroscopy. Biophys. J. 2016, 111 (5), 999−1007. (55) Epand, R. M.; Epand, R. F. Non-Raft Forming SphingomyelinCholesterol Mixtures. Chem. Phys. Lipids 2004, 132 (1), 37−46. (56) Rog, T.; Pasenkiewicz-Gierula, M. Cholesterol-sphingomyelin Interactions: A Molecular Dynamics Simulation Study. Biophys. J. 2006, 91 (10), 3756−3767. (57) de Almeida, R. F. M.; Fedorov, A.; Prieto, M. Sphingomyelin/ Phosphatidylcholine/Cholesterol Phase Diagram: Boundaries and Composition of Lipid Rafts. Biophys. J. 2003, 85 (4), 2406−2416. (58) Quinn, P. J.; Wolf, C. Egg-Sphingomyelin and Cholesterol Form a Stoichiometric Molecular Complex in Bilayers of Egg-Phosphatidylcholine. J. Phys. Chem. B 2010, 114 (47), 15536−15545. (59) Quinn, P. J. Structure of Sphingomyelin Bilayers and Complexes with Cholesterol Forming Membrane Rafts. Langmuir 2013, 29 (30), 9447−9456. (60) Wydro, P. Sphingomyelin/Phosphatidylcholine/Cholesterol Monolayers - Analysis of the Interactions in Model Membranes and Brewster Angle Microscopy Experiments. Colloids Surf., B 2012, 93, 174−179. (61) Yuan, C. B.; Furlong, J.; Burgos, P.; Johnston, L. J. The Size of Lipid Rafts: An Atomic Force Microscopy Study of Ganglioside GM1 Domains in Sphingomyelin/DOPC/Cholesterol Membranes. Biophys. J. 2002, 82 (5), 2526−2535.

4844

DOI: 10.1021/acs.jpcb.7b01832 J. Phys. Chem. B 2017, 121, 4833−4844