Hydrogen-bonding propensities of sphingomyelin in solution and in a bilayer assembly: a molecular dynamics study

Sphingomyelin is enriched within lipid microdomains of the cell membrane termed lipid rafts. These microdomains play a part in regulating a variety of cellular events. Computer simulations of the hydrogen-bonding properties of sphingolipids, believed to be central to the organization of these domains, can delineate the possible molecular interactions that underlie this lipid structure. We have therefore used molecular dynamics simulations to unravel the hydrogen-bonding behavior of palmitoylsphingomyelin (PSM). A series of eight simulations of 3 ns each of a single PSM molecule in water showed that the sphingosine OH and NH groups can form hydrogen bonds with the phosphate oxygens of their own polar head, in agreement with NMR data. Simulations of PSM in a bilayer assembly were carried out for 8 ns with three different force field parameterizations. The major physico-chemical parameters of the simulated bilayer agree with those established experimentally. The sphingosine OH group was mainly involved in intramolecular hydrogen bonds, in contrast to the almost exclusive intermolecular hydrogen bonds formed by the amide NH moiety. During the bilayer simulations the intermolecular hydrogen bonds among lipids formed a dynamic network characterized by the presence of hydrogen-bonded lipid clusters of up to nine PSM molecules.