A coarse-grained model of dimethyl sulfoxide for molecular dynamics simulations with lipid membranes

Interface Water Assists in Dimethyl Sulfoxide Crossing and Poration in Model Bilayer

Understanding the mechanism of transport and pore formation by a commonly used cryoprotectant, dimethyl sulfoxide (DMSO), across cell membranes is fundamentally crucial for drug delivery and cryopreservation. To shed light on the mechanism and thermodynamics of pore formation and crossing behavior of DMSO, extensive all-atom molecular dynamics simulations of 1,2-dimyristoyl-rac-glycero-3-phosphocholine (DMPC) bilayers are performed at various concentrations of DMSO at a temperature above the physiological temperature. Our results unveil that DMSO partially depletes water from the interface and positions itself between lipid heads without full dehydration. This induces a larger area per headgroup, increased disorder, and enhanced fluidity without any disintegration even at the highest DMSO concentration studied. The enhanced disorder fosters local fluctuations at the interface that nucleate dynamic and transient pores. The potential of mean force (PMF) of DMSO crossing is derived from two types of biased simulations: a single DMSO pulling using the umbrella sampling technique and a cylindrical pore formation using the recently developed chain reaction coordinate method. In both cases, DMSO crossing encounters a barrier attributed to unfavorable polar nonpolar interactions between DMSO and lipid tails. As the DMSO concentration increases, the barrier height reduces along with the faster lateral and perpendicular diffusion of DMSO suggesting favorable permeation. Our findings suggest that the energy required for pore formation decreases when water assists in the formation of DMSO pores. Although DMSO displaces water from the interface toward the far interface region without complete dehydration, the presence of interface water diminishes pore formation free energy. The existence of interface water leads to the formation of a two-dimensional percolated water-DMSO structure at the interface, which is absent otherwise. Overall, these insights into the mechanism of DMSO crossing and pore formation in the bilayer will contribute to understanding cryoprotectant behavior under supercooled conditions in the future.

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.

Application of molecular dynamics simulation in self-assembled cancer nanomedicine

Self-assembled nanomedicine holds great potential in cancer theragnostic. The structures and dynamics of nanomedicine can be affected by a variety of non-covalent interactions, so it is essential to ensure the self-assembly process at atomic level. Molecular dynamics (MD) simulation is a key technology to link microcosm and macroscale. Along with the rapid development of computational power and simulation methods, scientists could simulate the specific process of intermolecular interactions. Thus, some experimental observations could be explained at microscopic level and the nanomedicine synthesis process would have traces to follow. This review not only outlines the concept, basic principle, and the parameter setting of MD simulation, but also highlights the recent progress in MD simulation for self-assembled cancer nanomedicine. In addition, the physicochemical parameters of self-assembly structure and interaction between various assembled molecules under MD simulation are also discussed. Therefore, this review will help advanced and novice researchers to quickly zoom in on fundamental information and gather some thought-provoking ideas to advance this subfield of self-assembled cancer nanomedicine.

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.

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.

How Do Urea and Trimethylamine N-Oxide Influence the Dehydration-Induced Phase Transition of a Lipid Membrane?

Living organisms are often exposed to extreme dehydration, which is detrimental to the structure and function of the cell membrane. The lipid membrane undergoes fluid-to-gel phase transition due to dehydration and thus loses fluidity and functionality. To protect the fluid phase of the bilayer these organisms adopt several strategies. Enhanced production of small polar organic solutes (also called osmolytes) is one such strategy. Urea and trimethylamine N-oxide (TMAO) are two osmolytes found in different organisms combating osmotic stress. Previous experiments have found that both these osmolytes have strong effects on lipid membrane under different hydration conditions. Urea prevents the dehydration-induced phase transition of the lipid membrane by directly interacting with the lipids, while TMAO does not inhibit the phase transition. To provide atomistic insights, we have carried out all-atom molecular dynamics (MD) simulation of a lipid membrane under varying hydration levels and studied the effect of these osmolytes on different structural and dynamic properties of the membrane. This study suggests that urea significantly inhibits the dehydration-induced fluid-to-gel phase transition by strongly interacting with the lipid membrane via hydrogen bonds, which balances the reduced lipid hydration due to the decreasing water content. In contrast, TMAO is excluded from the membrane surface due to unfavorable interaction with the lipids. This induces further dehydration of the lipids which reinforces the fluid-to-gel phase transition. We have also studied the counteractive role of TMAO on the effect of urea on lipid membrane when both the osmolytes are present. TMAO draws some urea molecules out of the membrane and thereby reduces the effect of urea on the lipid membrane at lower hydration levels. This is similar to the counteraction of urea's deleterious effects on protein by TMAO. All these observations are consistent with the experimental results and thus provide deep molecular insights into the role of these osmolytes in protecting the fluid phase of the membrane, the key survival strategy against osmotic-stress-induced dehydration.

A Comparative Study on DMSO-Induced Modulation of the Structural and Dynamical Properties of Model Bilayer Membranes

Modulating the structures and properties of biomembranes via permeation of small amphiphilic molecules is immensely important, having diverse applications in cell biology, biotechnology, and pharmaceuticals, because their physiochemical and biological interactions lead to new pathways for transdermal drug delivery and administration. In this work, we have elucidated the role of dimethyl sulfoxide (DMSO), broadly used as a penetration-enhancing agent and cryoprotective agent on model lipid membranes, using a combination of fluorescence microscopy and time-resolved fluorescence spectroscopy. Spatially resolved fluorescence lifetime imaging microscopy (FLIM) has been employed to unravel how the fluidity of the DMSO-induced bilayer regulates the structural alteration of the vesicles. Moreover, we have also shown that the dehydration effect of DMSO leads to weakening of the hydrogen bond between lipid headgroups and water molecules and results in faster solvation dynamics as demonstrated by femtosecond time-resolved fluorescence spectroscopy. It has been gleaned that the water dynamics becomes faster because bilayer rigidity decreases in the presence of DMSO, which is also supported by time-resolved rotational anisotropy measurements. The enhanced diffusivity and increased membrane fluidity in the presence of DMSO are further ratified at the single-molecule level through fluorescence correlation spectroscopy (FCS) measurements. Our results indicate that while the presence of DMSO significantly affects the 1,2-dimyristoyl-rac-glycero-3-phosphocholine (DMPC) and 1,2-dipalmitoyl-rac-glycero-3-phosphatidylcholine (DPPC) bilayers, it has a weak effect on 1,2-dimyristoyl-sn-glycero-3-phospho-rac-glycerol (DMPG) vesicles, which might explain the preferential interaction of DMSO with the positively charged choline group present in DMPC and DPPC vesicles. The experimental findings have also been further verified with molecular dynamics simulation studies. Moreover, it has been observed that DMSO is likely to have a differential effect on heterogeneous bilayer membranes depending on the structure and composition of their headgroups. Our results illuminate the importance of probing the lipid structure and composition of cellular membranes in determining the effects of cryoprotective agents.