Multiscale molecular dynamics simulations of sodium dodecyl sulfate micelles: from coarse-grained to all-atom resolution

DPD simulations of anionic surfactant micelles: a critical role for polarisable water models

We investigate the effects of polarisable water models in dissipative particle dynamics (DPD) simulations, focussing on the influence these models have on the aggregation behaviour of sodium dodecyl sulfate solutions. Studies in the literature commonly report that DPD approaches underpredict the micellar aggregation number of ionic surfactants compared to experimental values. One of the proposed reasons for this discrepancy is that existing water models are insufficient to accurately model micellar solutions, as they fail to account for structural changes in water close to micellar surfaces. We show that polarisable DPD water models lead to more realistic counterion behaviour in micellar solutions, including the degree of counterion disassociation. These water models can also accurately reproduce changes in the dielectric constant of surfactant solutions as a function of concentration. We find evidence that polarisable water leads to the formation of more stable micelles at higher aggregation numbers. However, we also show that the choice of water model is not responsible for the underestimated aggregation numbers observed in DPD simulations. This finding addresses a key question in the literature surrounding the importance of water models in DPD simulations of ionic micellar solutions.

Revealing the Influence of Salts on the Hydration Structure of Ionic SDS Micelles by Contrast-Variation Small-Angle Neutron Scattering

The influence of lithium chloride (LiCl) on the hydration structure of anionic micelles of sodium dodecyl sulfate (SDS) in water was studied using the contrast-variation small-angle neutron scattering (SANS) technique. In the past, extensive computational studies have shown that the distribution of invasive water plays a critical role in the self-organization of SDS molecules and the stability of the assemblies. However, in past scattering studies the degree of the hydration level was not examined explicitly. Here, a series of contrast-variation SANS data was analyzed to extract the intramicellar radial distributions of invasive water and SDS molecules from the evolving spectral lineshapes caused by the varying isotopic ratios of water. By addressing the intramicellar inhomogeneous distributions of water and SDS molecules, a detailed description of how the counterion association influences the micellization behavior of SDS molecules is provided. The extension of our method can be used to provide an in-depth insight into the micellization phenomenon, which is commonly found in many soft matter systems.

Revisiting OPLS-AA Force Field for the Simulation of Anionic Surfactants in Concentrated Electrolyte Solutions

Hereby, we developed a set of nonbonded parameters within all-atom optimized potentials for liquid simulations (OPLS-AA) force field for the simulation of concentrated electrolyte solutions of anionic surfactants. More specifically, the aim of this paper is to assess the performance of five sets of atomic charges calculated using different population analyses (DDEC6, CHelpG, CHelpG-SMD, RESP, and CM5), as well as the original set of charges used in the literature for sodium dodecyl sulfate (SDS) simulation. Recently, Farafonov et al. have revised the SDS OPLS-AA force field; however, we were unable to obtain the experimental rodlike micelles using this parameter set on long time scale. In fact, the initial SDS bilayer micelle adopted a rodlike shape transiently and then broke down into spherical micelles. Updating OPLS-AA force field with DDEC6, CHelpG, and CHelpG-SMD charges resulted in stable rod micelles for a long simulation time (1 μs). The atomic charges of Farafonov (taken from Shelley et al.), RESP, and CM5 could not correctly describe SDS in concentrated electrolyte solutions. Analysis of the interaction of SDS with the counterions and solvent highlights the role of a balance of the intermolecular forces that must be met to describe adequately the anionic surfactant electrolyte solutions. Further, the optimization of the SDS Lennard-Jones parameters enabled the Farafonov set to properly reproduce the experimental rod micelle. In addition, we have examined the performance of different parameters of sodium ions: the first developed based on the Kirkwood-Buff integrals (KBI) and the second developed by Joung et al. The excessive ion pairing caused by KBI parameters screens significantly SDS-water interactions, which stabilize the rod micelle. Further, a tight interaction of the Na⁺-SDS head group resulted in stabilization of the bilayer micelle as observed in the case of Na⁺ parameters developed by Joung et al.

Hydration of diblock copolymer micelles: Effects of hydrophobicity and co-solvent

Diblock polymer micelles dispersed in an aqueous environment are being actively investigated for various applications, but there is only a qualitative understanding of the effect of the chemical structure on the micelle hydration and water dynamics as these properties are difficult to assess experimentally. Using all-atom molecular dynamics simulations, we investigate aqueous solutions of three comparable in size diblock copolymer micelles with core-forming blocks of different hydrophobicity: polybutadiene (PB), polycaprolactone (PCL), and polytetrahydrofuran (pTHF) with the same hydrophilic block, polyethylene oxide (PEO). We found that core-block hydrophobicity and ability to form hydrogen bonds with water strongly affect the water dynamics near the core: water molecules spend considerably less time in contact with the PB block than with PCL and pTHF blocks. We obtained polymer and solvent volume fraction profiles and determined that the interfacial width systematically increases with a decrease of core block hydrophobicity with water penetration into the core being negligible for PB-PEO and PCL-PEO micelles, while for pTHF-PEO micelles the interface is more diffuse and there is a noticeable penetration of water (17% by volume). For PCL-PEO micelles, which are commonly used in biomedical applications, we also investigated tetrahydrofuran (THF) penetration into the micelles from mixed THF:water solution at early stages of micelle dissolution. We found an inhomogeneous solvent distribution with a maximum of THF volume fraction in the interfacial core-corona region and partial exclusion from the PEO corona, which slows down micelle dissolution. These results can have important implications for micelle stability and use in biomedical applications.

Protobiotic Systems Chemistry Analyzed by Molecular Dynamics

Systems chemistry has been a key component of origin of life research, invoking models of life's inception based on evolving molecular networks. One such model is the graded autocatalysis replication domain (GARD) formalism embodied in a lipid world scenario, which offers rigorous computer simulation based on defined chemical kinetics equations. GARD suggests that the first pre-RNA life-like entities could have been homeostatically-growing assemblies of amphiphiles, undergoing compositional replication and mutations, as well as rudimentary selection and evolution. Recent progress in molecular dynamics has provided an experimental tool to study complex biological phenomena such as protein folding, ligand-receptor interactions, and micellar formation, growth, and fission. The detailed molecular definition of GARD and its inter-molecular catalytic interactions make it highly compatible with molecular dynamics analyses. We present a roadmap for simulating GARD's kinetic and thermodynamic behavior using various molecular dynamics methodologies. We review different approaches for testing the validity of the GARD model by following micellar accretion and fission events and examining compositional changes over time. Near-future computational advances could provide empirical delineation for further system complexification, from simple compositional non-covalent assemblies towards more life-like protocellular entities with covalent chemistry that underlies metabolism and genetic encoding.

The Fluidity of Phosphocholine and Maltoside Micelles and the Effect of CHAPS

The dynamics of phosphocholine and maltoside micelles, detergents frequently used for membrane protein structure determination, were investigated using electron paramagnetic resonance of spin probes doped into the micelles. Specifically, phosphocholines are frequently used detergents in NMR studies, and maltosides are frequently used in x-ray crystallography structure determination. Beyond the structural and electrostatic differences, this study aimed to determine whether there are differences in the local chain dynamics (i.e., fluidity). The nitroxide probe rotational dynamics in longer chain detergents is more restricted than in shorter chain detergents, and maltoside micelles are more restricted than phosphocholine micelles. Furthermore, the micelle microviscosity can be modulated with mixtures, as demonstrated with mixtures of 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate with n-dodecylphosphocholine, n-tetradecylphosphocholine, n-decyl-β-D-maltoside, or n-dodecyl-β-D-maltoside. These results indicate that observed differences in membrane protein stability in these detergents could be due to fluidity in addition to the already determined structural differences.

Peptide-surfactant interactions: A combined spectroscopic and molecular dynamics simulation approach

In the present contribution, we report a combined spectroscopic and computational approach aiming to unravel at atomic resolution the effect of the anionic SDS detergent on the structure of two model peptides, the α-helix TrpCage and the β-stranded TrpZip. A detailed characterization of the specific amino acids involved is performed. Monomeric (single molecules) and micellar SDS species differently interact with the α-helix and β-stranded peptides, emphasizing the different mechanisms occurring below and above the critical aggregation concentration (CAC). Below the CAC, the α-helix peptide is fully unfolded, losing its hydrophobic core and its Asp-Arg salt bridge, while the β-stranded peptide keeps its native structure with its four Trp well oriented. Above the CAC, the SDS micelles have the same effect on both peptides, that is, destabilizing the tertiary structure while keeping their secondary structure. Our studies will be helpful to deepen our understanding of the action of the denaturant SDS on peptides and proteins.