Impact of amphiphilic molecules on the structure and stability of homogeneous sphingomyelin bilayer: Insights from atomistic simulations

Mechanistic Insight on Ethanol Driven Swelling and Disruption of Cholesterol Containing Biomimetic Vesicles From Coarse-Grained Molecular Dynamics

We have performed coarse-grained (CG) molecular dynamics (MD) simulations to delineate the impact of ethanol (EtOH) on cholesterol (CHOL) containing biomimetic bilayer and vesicle composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipids. We have first deduced the missing interaction parameters for the POPC-CHOL-EtOH-water system within the SPICA/SDK CG force-field (CG-FF). By monitoring the electron density profiles, the orientational order parameter, and reproducing the all-atom MD-derived free energy for the insertion of ethanol from the bulk aqueous phase to the hydrophobic core of the POPC-CHOL lipid bilayer, we successfully determined all the missing non-bonding interaction parameters for the POPC-CHOL-EtOH-water system. The proposed force field was applied to investigate the effect of ethanol at various concentrations on unilamellar vesicles composed of POPC and cholesterol. It was found that 40 mol% or more concentration of ethanol is required to disintegrate or rupture the POPC-CHOL vesicle membranes. While cholesterol offers some resilience against the detrimental effects of ethanol, we still observe an increase in vesicle size (swelling) and a contraction in the bilayer thickness (thinning) as ethanol concentration rises from 0 to 30 mol%. At ethanol concentrations exceeding 30 mol%, the vesicles become increasingly susceptible to disintegration due to enhanced penetration of ethanol and water molecules into the hydrophobic core of the membranes.

Computational Insights on the Assembly of the Dengue Virus Membrane-Capsid-RNA Complex

Dengue virus, an arbovirus from the genus Flavivirus in the family Flaviviridae, forms a nucleocapsid structure through interactions between its genome and multiple copies of the capsid protein. Experimental studies have confirmed the interaction between the viral capsid protein and lipid droplets, indicating a protein-lipid interaction. Cryo-EM studies show that in immature viruses, the nucleocapsid is located close to the viral membrane. This study uses multiple MD simulations to explore the orientation of the capsid protein relative to the lipid membrane, focusing on how the protein's hydrophobic pocket interacts with the membrane. We also investigated the interaction between the capsid protein and RNA, considering the effects of sequence length and identity. Finally, we construct a model of the lipid-protein-RNA complex, demonstrating that the capsid protein's hydrophobic pocket interacts with the membrane, while the positively charged H4 helix interacts with the negatively charged RNA. This research may identify crucial interactions for immature virus particle formation and provide insights for future therapeutic interventions.

Tracking Cholesterol Flip-Flop in Mammalian Plasma Membrane through Coarse-Grained Molecular Dynamics Simulations

Plasma membrane (PM) simulations at longer length and time scales at nearly atomistic resolution can provide invaluable insights into cell signaling, apoptosis, lipid trafficking, and lipid raft formation. We propose a coarse-grained (CG) model of a mammalian PM considering major lipid head groups distributed asymmetrically across the membrane bilayer and validate the model against bilayer structural properties from atomistic simulation. Using the proposed CG model, we identify a recurring pattern in the passive collective cholesterol transbilayer motion and study the individual cholesterol flip-flop events and associated pathways along with lateral ordering in the bilayer during a flip-flop event. We identify two discrete cholesterol flip-flop pathways: (i) a systematic rototranslational pathway and (ii) intraleaflet inversion followed by interleaflet translation (or reverse). We observe a periodic cholesterol enrichment in the exoplasmic leaflet of the PM bilayer and examine the underlying cholesterol-lipid affinities. We observe closer association between cholesterol and palmitoylsphingomyelin (PSM) lipid, relative to other lipids, and conclude that the cholesterol enrichment in the exoplasmic leaflet can be attributed to higher PSM content in that leaflet, together leading to formation of short-lived PSM-cholesterol-rich domains.

Fatty Alcohol Membrane Model for Quantifying and Predicting Amphiphilicity

Amphiphilicity is an important property for drug development and self-assembly. This paper introduces a general approach based on a simple fatty alcohol (dodecanol) membrane model that can be used to quantify the amphiphilicity of small molecules that are in good agreement with experimental surface tension data. By applying the model to a systematic series of compounds, it was possible to elucidate the effect of different motifs on amphiphilicity. The results further indicate that amphiphilicity correlates strongly with water-octanol partition coefficients (logP) for the 29 organic molecules examined in the 0 < logP < 4 range. Importantly, the simulation of the model membrane is an order of magnitude faster than a phospholipid membrane (e.g., 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) simulation and offers a simple atomistic approach for quantifying and predicting amphiphilicity of small drug-like molecules that could be used in quantitative structure-activity relationship studies.

Ultrafast and Reversible Superwettability Switching of 3D Graphene Foams via Solvent-Exclusive Plasma Treatments

For highly active electron transfer and ion diffusion, controlling the surface wettability of electrically and thermally conductive 3D graphene foams (3D GFs) is required. Here, we present ultrasimple and rapid superwettability switching of 3D GFs in a reversible and reproducible manner, mediated by solvent-exclusive microwave arcs. As the 3D GFs are prepared with vapors of nonpolar acetone or polar water exclusively, short microwave radiation (≤10 s) leads to plasma hotspot-mediated production of methyl and hydroxyl radicals, respectively. Upon immediate radical chemisorption, the 3D surfaces become either superhydrophobic (water contact angle = ∼170°) or superhydrophilic (∼0°), and interestingly, the wettability transition can be repeated many times due to the facile exchange between previously chemisorbed and newly introduced radicals via the formation of methanol-like intermediates. When 3D GFs of different surficial polarities are incorporated into electric double-layer capacitors with nonpolar ionic liquids or polar aqueous electrolytes, the polarity matching between graphene surfaces and electrolytes results in ≥548.0 times higher capacitance compared to its mismatching at ≥0.5 A g-1, demonstrating the significance of wettability-controlled 3D GFs.

Modulation of bacterial membranes and cellular macromolecules by dimethyl sulfoxide: A dose-dependent study providing novel insights

Using Escherichia coli as a model, this manuscript delves into the intricate interactions between dimethyl sulfoxide (DMSO) and membranes, cellular macromolecules, and the effects on various aspects of bacterial physiology. Given DMSO's wide-ranging use as a solvent in microbiology, we investigate the impacts of both non-growth inhibitory (1.0 % and 2.5 % v/v) and slightly growth-inhibitory (5.0 % v/v) concentrations of DMSO. The results demonstrate that DMSO causes alterations in bacterial membrane potential, influences the electrochemical characteristics of the cell surface, and exerts substantial effects on the composition and structure of cellular biomolecules. Genome-wide gene expression data from DMSO-treated E. coli was used to further investigate and bolster the results. The findings of this study provide valuable insights into the complex relationship between DMSO and biological systems, with potential implications in drug delivery and cellular manipulation. However, it is essential to exercise caution when utilizing DMSO to enhance the solubility and delivery of bioactive compounds, as even at low concentrations, DMSO exerts non-inert effects on cellular macromolecules and processes.

A critical review of ethanol effects on neuronal firing: A metabolic perspective

Ethanol metabolism is relatively understudied in neurons, even though changes in neuronal metabolism are known to affect their activity. Recent work demonstrates that ethanol is preferentially metabolized over glucose as a source of carbon and energy, and it reprograms neurons to a state of reduced energy potential and diminished capacity to utilize glucose once ethanol is exhausted. Ethanol intake has been associated with changes in neuronal firing and specific brain activity (EEG) patterns have been linked with risk for alcohol use disorder (AUD). Furthermore, a haplotype of the inwardly rectifying potassium channel subunit, GIRK2, which plays a critical role in regulating excitability of neurons, has been linked with AUD and shown to be directly regulated by ethanol. At the same time, overexpression of GIRK2 prevents ethanol-induced metabolic changes. Based on the available evidence, we conclude that the mechanisms underlying the effects of ethanol on neuronal metabolism are a novel target for developing therapies for AUD.

Stability of Cytoplasmic Membrane of Escherichia coli Bacteria in Aqueous and Ethanolic Environment

The mechanism of action of any antibacterial agent or disinfectant depends largely on their interaction with the bacterial membrane. Herein, we use the SPICA (surface property fitting coarse graining) force-field and develop a coarse-grained (CG) model for the structure of the cytoplasmic membrane of Escherichia coli (E. coli) and its interaction with water and ethanol. We elucidate the impact of different concentrations of ethanol on the cytoplasmic membrane bilayers and vesicles of E. coli using the CG molecular dynamics (CG MD) simulations. Our modeling approach first focuses on the parametrization of the required force-field for POPG lipid and its interaction with water, ethanol, and POPE lipid. Subsequently, the structural stability of the E. coli bacterial membrane in the presence of high and low concentrations of ethanol is delineated. Both flat bilayers as well as vesicles of E. coli membrane were considered for the CG MD. Our results reveal that, at low ethanol concentrations (<30 mol %), the size of the E. coli vesicles increases with discernible deformations in their shapes. Because of ethanol-induced interdigitation, thinning of the E. coli vesicular membrane is also observed. However, at higher ethanol concentrations (>30 mol %), the integrity of the vesicles is lost because of deteriorating invasion of ethanol molecules into the vesicle bilayer and significant weakening of lipid-lipid interactions. At higher ethanol concentrations (40 and 70 mol %), both the multivesicle and single-vesicle bacterial membranes exhibit a similar rupturing pattern wherein the extraction of lipids from the membrane and formation of aggregates of the component lipids are observed. These aggregates consist of polar head groups of 3-5 POPE/POPG lipids with intertwined nonpolar tails.

Influence of docosahexaenoic acid on the interfacial behavior of cholesterol-containing lipid membranes: Interactions with small amphiphiles and hydration properties

Cholesterol is known to significantly modify both the structural and the dynamical properties of lipid membranes. On one side, the presence of free cholesterol molecules has been determined to stiffen the membrane bilayer by stretching the hydrophobic tails. Additionally, recent experimental and computational findings have made evident the fact that cholesterol also alters the dynamics and the hydration properties of the polar head groups of DPPC model lipid membranes. In turn, we have recently shown that the Omega-3 fatty acid docosahexaenoic acid, DHA, counteracts the effect of cholesterol on DPPC membrane's mechanical properties by fluidizing the bilayer. However, such behavior represents in fact a global outcome dominated by the larger lipid hydrophobic tails that neither discriminates between the different parts of the membrane nor elucidates the effect on membrane hydration and binding properties. Thus, we now perform molecular dynamics simulations to scrutinize the influence of DHA on the interfacial behavior of cholesterol-containing lipid membranes by characterizing their hydration properties and their binding to amphiphiles. We find that while cholesterol destabilizes interactions with amphiphiles and slightly weakens the lipid's hydration layer, the incorporation of DHA practically restores the interfacial behavior of pure DPPC.

Wrapping-Trapping versus Extraction Mechanism of Bactericidal Activity of MoS2 Nanosheets against Staphylococcus aureus Bacterial Membrane

The promising broad-spectrum antibacterial activity of two-dimensional molybdenum disulfide (2D MoS₂) has been widely recognized in the past decade. However, a comprehensive understanding of how the antibacterial pathways opted by the MoS₂ nanosheets varies with change in lipid compositions of different bacterial strains is imperative to harness their full antibacterial potential and remains unexplored thus far. Herein, we present an atomistic molecular dynamics (MD) study to investigate the distinct modes of antibacterial action of MoS₂ nanosheets against Staphylococcus aureus (S. aureus) under varying conditions. We observed that the freely dispersed nanosheets readily adhered to the bacterial membrane outer surface and opted for an unconventional surface directed "wrapping-trapping" mechanism at physiological temperature (i.e., 310 K). The adsorbed nanosheets mildly influenced the membrane structure by originating a compact packing of the lipid molecules present in its direct contact. Interestingly, these surface adsorbed nanosheets exhibited extensive phospholipid extraction to their surface, thereby inducing transmembrane water passage analogous to the cellular leakage, even at a slight increment of 20 K in the temperature. The strong van der Waals interactions between lipid fatty acyl tails and MoS₂ basal planes were primarily responsible for this destructive phospholipid extraction. In addition, the MoS₂ nanosheets bound to an imaginary substrate, controlling their vertical alignment, demonstrated a "nano-knives" action by spontaneously piercing inside the membrane core through their sharp corner, subsequently causing localized lipid ordering in their vicinity. The larger nanosheet produced a more profound deteriorating impact in all of the observed mechanisms. Keeping the existing knowledge about the bactericidal activity of 2D MoS₂ in view, our study concludes that their antibacterial activity is strongly governed by the lipid composition of the bacterial membrane and can be intensified either by controlling the nanosheet vertical alignment or by moderately warming up the systems.

Chronic alcohol metabolism results in DNA repair infidelity and cell cycle-induced senescence in neurons

Chronic binge-like drinking is a risk factor for age-related dementia, however, the lasting and irreversible effect of alcohol on the brain remains elusive. Transcriptomic changes in brain cortices revealed pro-ageing hallmarks upon chronic ethanol exposure and these changes predominantly occur in neurons. The changes are attributed to a prioritized ethyl alcohol oxidation in these cells via the NADPH-dependent cytochrome pathway. This hijacks the folate metabolism of the 1-carbon network which supports the pathway choice of DNA repair via the non-cell cycle-dependent mismatch repair networks. The lost-in-function of such results in the de-inactivation of the less preferred cell cycle-dependent homologous recombination (HR) repair, forcing these post-mitotic cells to re-engage in a cell cycle-like process. However, mature neurons are post-mitotic. Therefore, instead of successfully completing a full round of cell cycle which is necessary for the completion of HR-mediated repair; these cells are arrested at checkpoints. The resulting persistence of repair intermediates induces and promotes the nuclear accumulation of p21 and cyclin B-a trigger for permanent cell cycle exits and irreversible senescence response. Supplementation of bioactive 5-methyl tetrahydrofolate simultaneously at times with ethyl alcohol exposure supports the fidelity of the 1-carbon network and hence the activity of the mismatch repair. This prevents aberrant and irreversible cell cycle re-entry and senescence events of neurons. Together, our findings offer a direct connection between binge-drinking behaviour and its irreversible impact on the brain, which makes it a potential risk factor for dementia.

Estrogenic mixtures induce alterations in lipidomic profiles in the gonads of female oysters

This study aimed to reveal possible alterations to lipidomic profiles in Sydney rock oysters, Saccostrea glomerata, exposed to estrogenic mixtures (i.e., estrone, E1; 17β-estradiol, E2; estriol, E3; 17α-ethinylestradiol, EE2; bisphenol A, BPA; 4-t-octylphenol, 4-t-OP; and 4-nonylphenol, 4-NP) at "low" and "high" concentrations, typical of those detected in Australian and global receiving waters. A seven-day acute exposure window exhibited significantly lower abundances of many non-polar metabolites in digestive gland, gills, and gonads. Overall, there was a strong effect of the carrier solvent ethanol (despite a low exposure of 0.0002%), with all solvent containing treatments exhibiting lower abundances of lipidic metabolites, especially in the gill and digestive gland. No significant changes of the lipidome were exhibited in the male gonad by estrogenic exposure. However, in the female gonad, significant reductions of phospholipids and phosphatidylcholine were associated with exposure to high estrogenic mixtures. We hypothesise that the decreases in these phospholipids in the female gonad may be attributable to 1) lower algal consumption and thus lower uptake of lipidic building blocks; 2) a reduction of available substrates for phospholipid and phosphatidylcholine synthesis; and/or 3) induction of reactive oxygen species via estrogen metabolism, which may cause lipid peroxidation and lower abundance of phospholipids.

Deciphering Ethanol-Driven Swelling, Rupturing, Aggregation, and Fusion of Lipid Vesicles Using Coarse-Grained Molecular Dynamics Simulations

Traditionally, liquid ethanol is known to enhance the permeability of lipid membranes and causes vesicle aggregation and fusion. However, how the amphiphilic ethanol molecules perturb the lipid vesicles to facilitate their aggregation or fusion has not been addressed at any level of molecular simulations. Herein, not only have we developed a coarse-grained (CG) model for liquid ethanol, its aqueous mixture, and hydrated lipid membranes for molecular dynamics (MD) simulations, but also utilized it to delineate the aggregation and fusion of lipid vesicles using CG-MD simulations with multimillion particles. We have systematically parametrized the force-field for pure ethanol and its interactions with hydrated POPC and POPE model lipid membranes. In this process, we have successfully reproduced the bulk ethanol structure and concentration-dependent density of aqueous ethanol. To quantify the interaction of ethanol with lipid membranes, we have reproduced the transfer free energy of the ethanol molecule across the hydrated bilayers, and the concentration-dependent distribution of ethanol molecules across the lipid bilayers. After having acceptable force-field parameters for ethanol-membrane interactions, we have checked the effect of ethanol toward the vesicles comprising POPC lipids. We observe a rapid increase in the size of the POPC lipid vesicles with increasing amounts of ethanol up to 30 mol %. We unambiguously observe swelling and decrease in the thickness of the POPC vesicles with increasing amounts of ethanol up to 30 mol %, beyond which the vesicles begin to lose their integrity and rupture at higher mol % of ethanol. The fusion study of two vesicles demonstrates that fused vesicles can be obtained from 20 to 30 mol % of ethanol provided that they are brought closer than a critical distance at a particular mol %. The multivesicle simulations show that along with the increase in the sizes of vesicles the propensity of vesicle aggregation increases as the mol % of ethanol increases.

Mechanistic Insight on BioIL-Induced Structural Alterations in DMPC Lipid Bilayer

The emerging application risks of traditional ionic liquids (ILs) toward the ecosystem have changed the perception regarding their greenness. This resulted in the exploration of their more biocompatible alternatives known as biocompatible ILs (BioILs). Here, we have investigated the impact of two such biocompatible cholinium amino acid-based ILs on the structural behavior of model homogeneous DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) lipid bilayer using all-atom molecular dynamics simulation technique. Two classic cholinium-amino acid-based ILs, cholinium glycinate ([Ch][Gly]) and cholinium phenylalaninate ([Ch][Phe]), which differ only by the side chain lengths and hydrophobicity of the anions, have been utilized in the present work. Simultaneous analysis of the bilayer structural properties reveals that the existence of [Ch][Gly] BioIL above a particular concentration induces phase transition from fluid phase to gel phase in the DMPC lipid bilayer. Such a freezing of lipid bilayer upon the exposure to concentrated aqueous solution of [Ch][Gly] BioIL indicates the harmfulness of this BioIL toward the cell membranes majorly containing DMPC lipids, as the cell freezing can negatively affect its stability and functionality. Despite having a more hydrophobic amino acid side chain of [Phe]^- anion in [Ch][Phe], in the case of bilayer-[Ch][Phe] systems we observe the minimal impact of [Ch][Phe] BioIL on the DMPC bilayer properties up to 10 mol % concentration. In the presence of these BioIL, we observe the thickening of the bilayer and accumulation of the cations and anions of the BioILs at the interface of DMPC lipid heads and tails. The transfer free-energy profile of a [Phe]^- anion from aqueous phase to membrane center also indicates the anion partitioning at lipid head-tail interface and its inability to penetrate in the lipid membrane tail region. In contrast, the free-energy profile for a [Gly]^- anion offers a very high energy barrier to the insertion of [Gly]^- into the membrane interior, leading to accumulation of [Gly]^- anions at the lipid head-water region.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

The role of the envelope protein in the stability of a coronavirus model membrane against an ethanolic disinfectant

Ethanol is highly effective against various enveloped viruses and can disable the virus by disintegrating the protective envelope surrounding it. The interactions between the coronavirus envelope (E) protein and its membrane environment play key roles in the stability and function of the viral envelope. By using molecular dynamics simulation, we explore the underlying mechanism of ethanol-induced disruption of a model coronavirus membrane and, in detail, interactions of the E-protein and lipids. We model the membrane bilayer as N-palmitoyl-sphingomyelin and 1-palmitoyl-2-oleoylphosphatidylcholine lipids and the coronavirus E-protein. The study reveals that ethanol causes an increase in the lateral area of the bilayer along with thinning of the bilayer membrane and orientational disordering of lipid tails. Ethanol resides at the head-tail region of the membrane and enhances bilayer permeability. We found an envelope-protein-mediated increase in the ordering of lipid tails. Our simulations also provide important insights into the orientation of the envelope protein in a model membrane environment. At ∼25 mol. % of ethanol in the surrounding ethanol-water phase, we observe disintegration of the lipid bilayer and dislocation of the E-protein from the membrane environment.

A Glimpse of the Structural Biology of the Metabolism of Sphingosine-1-Phosphate

As a key sphingolipid metabolite, sphingosine-1-phosphate (S1P) plays crucial roles in vascular and immune systems. It regulates angiogenesis, vascular integrity and homeostasis, allergic responses, and lymphocyte trafficking. S1P is interconverted with sphingosine, which is also derived from the deacylation of ceramide. S1P levels and the ratio to ceramide in cells are tightly regulated by its metabolic pathways. Abnormal S1P production causes the occurrence and progression of numerous severe diseases, such as metabolic syndrome, cancers, autoimmune disorders such as multiple sclerosis, and kidney and cardiovascular diseases. In recent years, huge advances on the structure of S1P metabolic pathways have been accomplished. In this review, we have got a glimpse of S1P metabolism through structural and biochemical studies of: sphingosine kinases, S1P transporters and S1P receptors, and the development of therapeutics targeting S1P signaling. The progress we summarize here could provide fresh perspectives to further the exploration of S1P functions and facilitate the development of therapeutic molecules targeting S1P signaling with improved specificity and therapeutic effects.

Molecular Dynamics Evaluation of the Effect of Cholinium Phenylalaninate Biocompatible Ionic Liquid on Biomimetic Membranes

Toward the search of sustainable green solvents, choline amino acid ([Ch][AA]) ionic liquids (ILs), mainly derived from renewable feedstocks, have emerged as a promising atoxic alternative to the conventional solvents. Recent studies have shown the remarkably benign nature of cholinium-based ILs against biomimetic phospholipid membranes. However, few of the contemporaneous experimental studies have contradicted the aforesaid ecofriendly nature of these ILs with anions comprising longer alkyl or aromatic tails. Aiming to understand the influence of amino acid side-chain variation in a particular bio-IL on biomembranes, herein, we have evaluated the effect of cholinium phenylalaninate ([Ch][Phe]) IL on the structural stability of homogeneous biomimetic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) lipid bilayers using atomistic molecular dynamics simulations. Although we observe spontaneous intercalation of aromatic rings of [Phe]^- anions in the hydrophobic region of the bilayers, the polar backbone of the anion remains coordinated with the lipid polar part through strong electrostatic and H-bonding interactions. Besides, the [Ch]⁺ cations get accumulated at the lipid-water interface to counter the excess negative charge density. The intercalation of the anionic rings causes significant perturbations in the lipid structural arrangement while still maintaining the bilayer integrity. The quantitative evaluation to probe the deteriorating effect of this bio-IL application establishes anions as the principal component causing the observed structural perturbations. The analysis of the structural properties along with the free energy assessment reveals the higher efficacy of [Ch][Phe] bio-IL to perturb the POPE bilayer structure than the POPC bilayer.

Molecular Structure of Sphingomyelin in Fluid Phase Bilayers Determined by the Joint Analysis of Small-Angle Neutron and X-ray Scattering Data

We have determined the fluid bilayer structure of palmitoyl sphingomyelin (PSM) and stearoyl sphingomyelin (SSM) by simultaneously analyzing small-angle neutron and X-ray scattering data. Using a newly developed scattering density profile (SDP) model for sphingomyelin lipids, we report structural parameters including the area per lipid, total bilayer thickness, and hydrocarbon thickness, in addition to lipid volumes determined by densitometry. Unconstrained all-atom simulations of PSM bilayers at 55 °C using the C36 CHARMM force field produced a lipid area of 56 Ų, a value that is 10% lower than the one determined experimentally by SDP analysis (61.9 Ų). Furthermore, scattering form factors calculated from the unconstrained simulations were in poor agreement with experimental form factors, even though segmental order parameter (S_CD) profiles calculated from the simulations were in relatively good agreement with S_CD profiles obtained from NMR experiments. Conversely, constrained area simulations at 61.9 Ų resulted in good agreement between the simulation and experimental scattering form factors, but not with S_CD profiles from NMR. We discuss possible reasons for the discrepancies between these two types of data that are frequently used as validation metrics for molecular dynamics force fields.

Molecular Dynamics Simulations of Membrane Permeability

This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.