Molecular conformation and bilayer pores in a nonionic surfactant lamellar phase studied with 1H-13C solid-state NMR and molecular dynamics simulations

Unravelling Guest Dynamics in Crystalline Molecular Organics Using 2H Solid-State NMR and Molecular Dynamics Simulation

2H solid-state NMR and atomistic molecular dynamics (MD) simulations are used to understand the disorder of guest solvent molecules in two cocrystal solvates of the pharmaceutical furosemide. Traditional approaches to interpreting the NMR data fail to provide a coherent model of molecular behavior and indeed give misleading kinetic data. In contrast, the direct prediction of the NMR properties from MD simulation trajectories allows the NMR data to be correctly interpreted in terms of combined jump-type and libration-type motions. Time-independent component analysis of the MD trajectories provides additional insights, particularly for motions that are invisible to NMR. This allows a coherent picture of the dynamics of molecules restricted in molecular-sized cavities to be determined.

Phase Behavior of Alkyl Ethoxylate Surfactants in a Dissipative Particle Dynamics Model

We present a dissipative particle dynamics (DPD) model capable of capturing the liquid state phase behavior of nonionic surfactants from the alkyl ethoxylate (C_nE_m) family. The model is based upon our recent work [Anderson et al. J. Chem. Phys. 2017, 147, 094503] but adopts tighter control of the molecular structure by setting the bond angles with guidance from molecular dynamics simulations. Changes to the geometry of the surfactants were shown to have little effect on the predicted micelle properties of sampled surfactants, or the water-octanol partition coefficients of small molecules, when compared to the original work. With these modifications the model is capable of reproducing the binary water-surfactant phase behavior of nine surfactants (C₈E₄, C₈E₅, C₈E₆, C₁₀E₄, C₁₀E₆, C₁₀E₈, C₁₂E₆, C₁₂E₈, and C₁₂E₁₂) with a good degree of accuracy.

Unveiling the phase behavior of CiEj non-ionic surfactants in water through coarse-grained molecular dynamics simulations

Poly(oxyethylene) alkyl ethers, usually denoted by CiEj surfactants, exhibit a rich phase behavior in water, self-assembling to form a variety of 3-D structures with a controllable morphology that find multiple applications across different industrial segments. Hence, being able to describe and understand the effect of molecular structure on the phase behavior of these systems is highly relevant for the efficient design of new materials and their applications. Considering the promising results obtained over the last decade using the MARTINI model to describe ethylene-oxide containing compounds, an extensive assessment of the ability of such a model to describe the phase behavior of CiEj in water was carried out and results are presented here. Given the overall poor temperature transferability of the MARTINI model, mostly due to the lack of an accurate representation of hydrogen bonding, simulations were carried out at a single temperature of 333 K, where most phases are expected to occur according to experiments. Different chain lengths of both the hydrophobic and hydrophilic moieties, spanning a wide range of hydrophilic-lipophilic balance values, were investigated and the phase diagrams of various CiEj surfactants explored over a wide concentration range. The model was able to satisfactorily describe the effect of surfactant structure and concentration on mesophase formation. The stability and dimensions of the obtained phases, and the prediction of some unique features such as the characterization of a singular lamellar phase are presented. The results obtained in this work highlight both the predictive ability and the transferability of the MARTINI forcefield in the description of such systems. Moreover, the model was shown to provide adequate descriptions of the micellar phase in terms of micelle dimensions, critical micelle concentration, and average aggregation number.

Probing n-Octyl α-d-Glycosides Using Deuterated Water in the Lyotropic Phase by Deuterium NMR

The lyotropic phase behavior of four common and easily accessible glycosides, n-octyl α-d-glycosides, namely, α-Glc-OC₈, α-Man-OC₈, α-Gal-OC₈, and α-Xyl-OC₈, was investigated. The presence of normal hexagonal (H_I), bicontinuous cubic (V_I), and lamellar (L_α) phases in α-Glc-OC₈ and α-Man-OC₈ including their phase diagrams in water reported previously was verified by deuterium nuclear magnetic resonance (²H NMR), via monitoring the D₂O spectra. Additionally, the partial binary phase diagrams and the liquid crystal structures formed by α-Gal-OC₈ and α-Xyl-OC₈ in D₂O were constructed and confirmed using small- and wide-angle X-ray scattering and ²H NMR. The average number of bound water molecules (n_b) per headgroup in the L_α phase was determined by the systematic measurement of the quadrupolar splitting of D₂O over a wide range of molar ratio values (glycoside/D₂O), especially at high glucoside composition. The number of bound water molecules bound to the headgroup was found to be around 1.5-2.0 for glucoside, mannoside, and galactoside, all of which possesses four OH groups. In the case of xyloside, which has only three OH groups, the bound water content is ∼2.0. Our findings confirmed that the bound water content of all n-octyl α-d-glycosides studied is lower compared to the number of possible hydrogen bonding sites possibly due to the fact that most of the OH groups are involved in intralayer interaction that holds the lipid assembly together.

Aqueous phase behavior of the PEO-containing non-ionic surfactant C12E6: A molecular dynamics simulation study

HYPOTHESIS

Non-ionic surfactants containing polyethylene oxide (PEO) chains are widely used in drug formulations, cosmetics, paints, textiles and detergents. High quality molecular dynamics models for PEO surfactants can give us detailed, atomic-scale information about the behavior of surfactant/water mixtures.

SIMULATIONS

We used two molecular dynamics force fields (FFs), 2016H66 and 53A6_DBW, to model the simple non-ionic PEO surfactant, hexaoxyethylene dodecyl ether (C₁₂E₆). We investigated surfactant/water mixtures that span the phase diagram of starting from randomly distributed arrangements. In some cases, we also started with prebuilt, approximate models. The simulations results were compared with the experimentally observed phase behavior.

FINDINGS

Overall, this study shows that the spontaneous self-assembly of PEO non-ionic surfactants into different colloidal structures can be accurately modeled with MD simulations using the 2016H66 FF although transitions to well-formed hexagonal phase are slow. Of the two FFs investigated, the 2016H66 FF better reproduces the experimental phase behavior across all regions of the C₁₂E_6/water phase diagram.

Spherical Micelles with Nonspherical Cores: Effect of Chain Packing on the Micellar Shape

Self-assembly of amphiphilic polymers into micelles is an archetypical example of a "self-confined" system due to the formation of micellar cores with dimensions of a few nanometers. In this work, we investigate the chain packing and resulting shape of C n -PEOx micelles with semicrystalline cores using small/wide-angle X-ray scattering (SAXS/WAXS), contrast-variation small-angle neutron scattering (SANS), and nuclear magnetic resonance spectroscopy (NMR). Interestingly, the n-alkyl chains adopt a rotator-like conformation and pack into prolate ellipses (axial ratio ϵ ≈ 0.5) in the "crystalline" region and abruptly arrange into a more spheroidal shape (ϵ ≈ 0.7) above the melting point. We attribute the distorted spherical shape above the melting point to thermal fluctuations and intrinsic rigidity of the n-alkyl blocks. We also find evidence for a thin dehydrated PEO layer (≤1 nm) close to the micellar core. The results provide substantial insight into the interplay between crystallinity and molecular packing in confinement and the resulting overall micellar shape.

Structure and Dynamics of Spherical and Rodlike Alkyl Ethoxylate Surfactant Micelles Investigated Using NMR Relaxation and Atomistic Molecular Dynamics Simulations

Predicting and controlling the properties of amphiphile aggregate mixtures require understanding the arrangements and dynamics of the constituent molecules. To explore these topics, we study molecular arrangements and dynamics in alkyl ethoxylate nonionic surfactant micelles by combining NMR relaxation measurements with large-scale atomistic molecular dynamics simulations. We calculate parameters that determine relaxation rates directly from simulated trajectories, without introducing specific functional forms to describe the dynamics. NMR relaxation rates, which depend on relative motions of interacting atom pairs, are influenced by wide distributions of dynamic time scales. We find that relative motions of neighboring atom pairs are rapid and liquidlike but are subject to structural constraints imposed by micelle morphology. Relative motions of distant atom pairs are slower than nearby atom pairs because changes in distances and angles are smaller when the moving atoms are further apart. Large numbers of atom pairs undergoing these slow relative motions contribute to predominantly negative cross-relaxation rates. For spherical micelles, but not for cylindrical micelles, cross-relaxation rates are positive only for surfactant tail atoms connected to the hydrophilic headgroup. This effect is related to the lower packing density of these atoms at the hydrophilic-hydrophobic boundary in spherical vs cylindrical arrangements, with correspondingly rapid and less constrained motion of atoms at the boundary.

Dynamics of Dislocations in Smectic A Liquid Crystals Doped with Nanoparticles

Edge dislocations are linear defects that locally break the positional order of the layers in smectic A liquid crystals. As in usual solids, these defects play a central role for explaining the plastic properties of the smectic A phase. This work focuses on the dynamical properties of dislocations in bulk samples prepared between two glass plates and in free-standing films. The emphasis will be put on the measurement of the mobility of edge dislocations in liquid crystals either pure or doped with nanoparticles. The experimental results will be compared to the existing models.

Transferable MARTINI Model of Poly(ethylene Oxide)

Motivated by the deficiencies of the previous MARTINI models of poly(ethylene oxide) (PEO), we present a new model featuring a high degree of transferability. The model is parametrized on (a) a set of 8 free energies of transfer of dimethoxyethane (PEO dimer) from water to solvents of varying polarity; (b) the radius of gyration in water at high dilution; and (c) matching angle and dihedral distributions from atomistic simulations. We demonstrate that our model behaves well in five different areas of application: (1) it produces accurate densities and phase behavior or small PEO oligomers and water mixtures; (2) it yields chain dimensions in good agreement with the experiment in three different solvents (water, diglyme, and benzene) over a broad range of molecular weights (∼1.2 kg/mol to 21 kg/mol); (3) it reproduces qualitatively the structural features of lipid bilayers containing PEGylated lipids in the brush and mushroom regime; (4) it is able to reproduce the phase behavior of several PEO-based nonionic surfactants in water; and (5) it can be combined with the existing MARTINI PS to model PS-PEO block copolymers. Overall, the new PEO model outperforms previous models and features a high degree of transferability.

Simulating Bilayers of Nonionic Surfactants with the GROMOS-Compatible 2016H66 Force Field

Polyoxyethylene glycol alkyl ether amphiphiles (C_iE_j) are important nonionic surfactants, often used for biophysical and membrane protein studies. In this work, we extensively test the GROMOS-compatible 2016H66 force field in molecular dynamics simulations involving the lamellar phase of a series of C_iE_j surfactants, namely C₁₂E₂, C₁₂E₃, C₁₂E₄, C₁₂E₅, and C₁₄E₄. The simulations reproduce qualitatively well the monitored structural properties and their experimental trends along the surfactant series, although some discrepancies remain, in particular in terms of the area per surfactant, the equilibrium phase of C₁₂E₅, and the order parameters of C₁₂E₃, C₁₂E₄, and C₁₂E₅. The polar head of the C_iE_j surfactants is highly hydrated, almost like a single polyethyleneoxide (PEO) molecule at full hydration, resulting in very compact conformations. Within the bilayer, all C_iE_j surfactants flip-flop spontaneously within tens of nanoseconds. Water-permeation is facilitated, and the bending rigidity is 4 to 5 times lower than that of typical phospholipid bilayers. In line with another recent theoretical study, the simulations show that the lamellar phase of C_iE_j contains large hydrophilic pores. These pores should be abundant in order to reproduce the comparatively low NMR order parameters. We show that their contour length is directly correlated to the order parameters, and we estimate that they should occupy approximately 7-10% of the total membrane area. Due to their highly dynamic nature (rapid flip-flops, high water permeability, observed pore formation), C_iE_j surfactant bilayers are found to represent surprisingly challenging systems in terms of modeling. Given this difficulty, the results presented here show that the 2016H66 parameters, optimized independently considering pure-liquid as well as polar and nonpolar solvation properties of small organic molecules, represent a good starting point for simulating these systems.

Structure and Dynamics of Nonionic Surfactant Aggregates in Layered Materials

The aggregation of surfactants on solid surfaces as they are adsorbed from solution is the basis of numerous technological applications such as colloidal stabilization, ore flotation, and floor cleaning. The understanding of both the structure and the dynamics of surfactant aggregates applies to the development of alternative ways of preparing hybrid layered materials. For this purpose, we study the adsorption of the triethylene glycol mono n-decyl ether (C₁₀E₃) nonionic surfactant onto a synthetic montmorillonite (Mt), an aluminosilicate clay mineral for organoclay preparation with important applications in materials sciences, catalysis, wastewater treatment, or as drug delivery. The aggregation mechanisms follow those observed in an analogous natural Mt, with the condensation of C₁₀E₃ in a bilayer arrangement once the surfactant self-assembles in a lamellar phase beyond the critical micelle concentration, underlining the importance of the surfactant state in solution. Solid-state ¹H nuclear magnetic resonance (NMR) at fast magic-angle spinning (MAS) and high magnetic field combined with¹H-¹³C correlation experiments and different types of ¹³C NMR experiments selectively probes mobile or rigid moieties of C₁₀E₃ in three different aggregate organizations: (i) a lateral monolayer, (ii) a lateral bilayer, and (iii) a normal bilayer. High-resolution ¹H{²⁷Al} CP-¹H-¹H spin diffusion experiments shed light on the proximities and dynamics of the different fragments and fractions of the intercalated surfactant molecules with respect to the Mt surface. ²³Na and ¹H NMR measurements combined with complementary NMR data, at both molecular and nanometer scales, precisely pointed out the location of the C₁₀E₃ ethylene oxide hydrophilic group in close contact with the Mt surface interacting through ion-dipole or van der Waals interactions.

Acyl Chain Disorder and Azelaoyl Orientation in Lipid Membranes Containing Oxidized Lipids

Oxidized phospholipids occur naturally in conditions of oxidative stress and have been suggested to play an important role in a number of pathological conditions due to their effects on a lipid membrane acyl chain orientation, ordering, and permeability. Here we investigate the effect of the oxidized phospholipid 1-palmitoyl-2-azelaoyl-sn-glycero-3-phosphocholine (PazePC) on a model membrane of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) using a combination of (13)C-(1)H dipolar-recoupling nuclear magnetic resonance (NMR) experiments and united-atom molecular dynamics (MD) simulations. The obtained experimental order parameter SCH profiles show that the presence of 30 mol % PazePC in the bilayer significantly increases the gauche content of the POPC acyl chains, therefore decreasing the thickness of the bilayer, although with no stable bilayer pore formation. The MD simulations reproduce the disordering effect and indicate that the orientation of the azelaoyl chain is highly dependent on its protonation state with acyl chain reversal for fully deprotonated states and a parallel orientation along the interfacial plane for fully protonated states, deprotonated and protonated azelaoyl chains having negative and positive SCH profiles, respectively. Only fully or nearly fully protonated azelaoyl chain are observed in the (13)C-(1)H dipolar-recoupling NMR experiments. The experiments show positive SCH values for the azelaoyl segments confirming for the first time that oxidized chains with polar termini adopt a parallel orientation to the bilayer plane as predicted in MD simulations.

Atomistic resolution structure and dynamics of lipid bilayers in simulations and experiments

Accurate details on the sampled atomistic resolution structures of lipid bilayers can be experimentally obtained by measuring C-H bond order parameters, spin relaxation rates and scattering form factors. These parameters can be also directly calculated from the classical atomistic resolution molecular dynamics simulations (MD) and compared to the experimentally achieved results. This comparison measures the simulation model quality with respect to 'reality'. If agreement is sufficient, the simulation model gives an atomistic structural interpretation of the acquired experimental data. Significant advance of MD models is made by jointly interpreting different experiments using the same structural model. Here we focus on phosphatidylcholine lipid bilayers, which out of all model membranes have been studied mostly by experiments and simulations, leading to the largest available dataset. From the applied comparisons we conclude that the acyl chain region structure and rotational dynamics are generally well described in simulation models. Also changes with temperature, dehydration and cholesterol concentration are qualitatively correctly reproduced. However, the quality of the underlying atomistic resolution structural changes is uncertain. Even worse, when focusing on the lipid bilayer properties at the interfacial region, e.g. glycerol backbone and choline structures, and cation binding, many simulation models produce an inaccurate description of experimental data. Thus extreme care must be applied when simulations are applied to understand phenomena where the interfacial region plays a significant role. This work is done by the NMRlipids Open Collaboration project running at https://nmrlipids.blogspot.fi and https://github.com/NMRLipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions

Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C-H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P-N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files ( https://zenodo.org/collection/user-nmrlipids ) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.

Stability of C(12)E(j) Bilayers Probed with Adhesive Droplets

The stability of model surfactant bilayers from the poly(ethylene glycol) mono-n-dodecyl ether (C12Ej) family was probed. The surfactant bilayers were formed by the adhesion of emulsion droplets. We generated C12Ej bilayers by forming water-in-oil (w/o) emulsions with saline water droplets, covered by the surfactant, in a silicone and octane oil mixture. Using microfluidics, we studied the stability of those bilayers. C12E1 allowed only short-lived bilayers whereas C12E2 bilayers were stable over a wide range of oil mixtures. At high C12E2 concentration, a two-phase region was displayed in the phase diagram: bilayers formed by the adhesion of two water droplets and Janus-like particles consisting of adhering aqueous and amphiphilic droplets. C12E8 and C12E25 did not mediate bilayer formation and caused phase inversion leading to o/w emulsion. With intermediate C12E4 and C12E5 surfactants, both w/o and o/w emulsions were unstable. We provided the titration of the C12E2 bilayer with C12E4 and C12E5 to study and predict their stability behavior.

NMR diffusion-encoding with axial symmetry and variable anisotropy: Distinguishing between prolate and oblate microscopic diffusion tensors with unknown orientation distribution

We introduce a nuclear magnetic resonance method for quantifying the shape of axially symmetric microscopic diffusion tensors in terms of a new diffusion anisotropy metric, DΔ, which has unique values for oblate, spherical, and prolate tensor shapes. The pulse sequence includes a series of equal-amplitude magnetic field gradient pulse pairs, the directions of which are tailored to give an axially symmetric diffusion-encoding tensor b with variable anisotropy bΔ. Averaging of data acquired for a range of orientations of the symmetry axis of the tensor b renders the method insensitive to the orientation distribution function of the microscopic diffusion tensors. Proof-of-principle experiments are performed on water in polydomain lyotropic liquid crystals with geometries that give rise to microscopic diffusion tensors with oblate, spherical, and prolate shapes. The method could be useful for characterizing the geometry of fluid-filled compartments in porous solids, soft matter, and biological tissues.