From homogeneous dispersion to micelles-a molecular dynamics simulation on the compromise of the hydrophilic and hydrophobic effects of sodium dodecyl sulfate in aqueous solution

Revealing the Formation Dynamics of Janus Polymer Particles: Insights from Experiments and Molecular Dynamics

Seeded emulsion polymerization is one of the best-known methods for preparing polymer particles with a controlled size, composition, and shape. It first requires the preparation of seed particles, which are then swollen with additional monomer (the same as the one used for the seed or a different one), to either increase the seed's size or change its morphology. The use of surfactants plays a central role in guaranteeing the required colloidal stability and contributing to the final shape and structure of the particles by lowering the interfacial energy between the polymer of the seed and the added monomer. We here study the polymerization of methyl methacrylate in the presence of polystyrene seed particles at various surfactant concentrations in the presence and absence of a surfactant (sodium dodecyl sulfate). We first show experimentally that the morphology of the colloidal particles can be tuned from Janus to core-shell, depending on the presence or absence of surfactant on the seeds particles' surface. Furthermore, using classical molecular dynamics simulations, we investigate the mechanism and behavior of the surfactants during the first stages of the polymerization process. We use a newly developed approach based on contact statistical analysis to confirm the critical role played by the organization of surfactant molecules on the surface of the seed particles in dictating the final particle morphology.

Molecular Dynamics Study on the Aggregation Behavior of Triton X Micelles with Different PEO Chain Lengths in Aqueous Solution

The aggregation structure of Triton X (TX) amphiphilic molecules in aqueous solution plays an important role in determining the various properties and applications of surfactant solutions. In this paper, the properties of micelles formed by TX-5, TX-114, and TX-100 molecules with different poly(ethylene oxide) (PEO) chain lengths in TX series of nonionic surfactants were studied via molecular dynamics (MD) simulation. The structural characteristics of three micelles were analyzed at the molecular level, including the shape and size of micelles, the solvent accessible surface area, the radial distribution function, the micelle configuration, and the hydration numbers. With the increase of PEO chain length, the micelle size and solvent accessible surface area also increase. The distribution probability of the polar head oxygen atoms on the surface of the TX-100 micelle is higher than that in the TX-5 or TX-114 micelle. In particular, the tail quaternary carbon atoms in the hydrophobic region are mainly located at the micelle exterior. For TX-5, TX-114, and TX-100 micelles, the interactions between micelles and water molecules are also quite different. These structures and comparisons at the molecular level contribute to the further understanding of the aggregation and applications of TX series surfactants.

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.

A novel method for constructing continuous intrinsic surfaces of nanoparticles

In recent years, the field of nanotechnology has become increasingly prevalent in the disciplines of science and engineering due to it’s abundance of application areas. Therefore, the ability to study and characterize these materials is more relevant than ever. Despite the wealth of simulation and modeling studies of nanoparticles reported in the literature, a rigorous description of the interface of such materials is rarely included in analyses which are pivotal to understanding interfacial behavior. We propose a novel method for constructing the continuous intrinsic surface of nanoparticles, which has been applied to a model system consisting of a sodium dodecyl sulfate micelle in the presence of testosterone propionate. We demonstrate the advantages of using our continuous intrinsic surface definition as a means to elucidate the true interfacial structure of the micelle, the interfacial properties of the hydrating water molecules, and the position of the drug (testosterone propionate) within the micelle. Additionally, we discuss the implications of this algorithm for future work in the simulation of nanoparticles. Graphical Abstract

Complexation of a Globular Protein, β-Lactoglobulin, with an Anionic Surfactant in Aqueous Solution

The complexation of a globular protein, β-lactoglobulin (BLG), and an anionic surfactant sodium dodecyl sulfate (SDS) in aqueous media was investigated using capillary zone electrophoresis, electrophoretic, static, and dynamic light scattering, and small-angle X-ray scattering in a considerably high protein concentration range (0.27 mM < C_P < 3 mM). On increasing the molar concentration C_R of the surfactant, cooperative binding of SDS to BLG starts at C_R/C_P ≈ 1; the BLG-SDS complex consists mainly of the BLG dimer and approximately 20 SDS molecules, where BLG takes a compact conformation similar to that of the native BLG up to C_R/C_P ≈ 20. At C_R/C_P higher than 30, the BLG dimer in the BLG-SDS complex dissociates into a unimer, but the dissociated BLG unimer still takes a compact conformation at least at 30 < C_R/C_P < 65.

Molecular dynamics study on the aggregation behaviour of different positional isomers of sodium dodecyl benzenesulphonate

Schematic diagram of anionic structures (1Φ12 and 6Φ12).

Tetraalkylammonium interactions with dodecyl sulfate micelles: a molecular dynamics study

2D number-density maps for S atoms of the polar head in the x–y plane.

CHARMM-GUI Martini Maker for Coarse-Grained Simulations with the Martini Force Field

Coarse-grained simulations are widely used to study large biological systems. Nonetheless, building such simulation systems becomes nontrivial, especially when membranes with various lipid types are involved. Taking advantage of the frameworks in all-atom CHARMM-GUI modules, we have developed CHARMM-GUI Martini Maker for building solution, micelle, bilayer, and vesicle systems as well as systems with randomly distributed lipids using the Martini force field. Martini Maker supports 82 lipid types and different flavors of the Martini force field, including polar and nonpolar Martini, Dry Martini, and ElNeDyn (an elastic network model for proteins). The qualities of the systems generated by Martini Maker are validated by simulations of various examples involving proteins and lipids. We expect Martini Maker to be a useful tool for modeling large, complicated biomolecular systems in a user-friendly way.

Coarse-grained molecular dynamics simulation of self-assembly and surface adsorption of ionic surfactants using an implicit water model

We perform coarse-grained molecular dynamics simulations for sodium dodecyl sulfate (SDS) surfactant using a modification of the Dry Martini force field (Arnarez et al. 2014) with implicit water. After inclusion of particle mesh Ewald (PME) electrostatics, an artificially high dielectric constant for water (ε(r) = 150), and reparameterization, we obtain structural and thermodynamic properties of SDS micelles that are close to those obtained from the standard Martini force field with explicit water, which in turn match those of atomistic simulations. The gains in computational efficiency obtained by removing explicit water allow direct simulations of the self-assembly of SDS in solution. We observe surfactant exchange among micelles and micelle fission and fusion and obtain realistic, equilibrated micelle size distributions at modest computational cost, as well as a transition to cylindrical micelles at high surfactant concentration or with added salt. We further apply this parametrized force field to study the adsorption of SDS onto hydrophobic surfaces and calculate the adsorption kinetics and equilibrium adsorption isotherm. The greatly increased speed of computation of surfactant self-assembly made possible by this Dry Martini method should allow future simulation of competitive adsorption of multiple surfactant species to surfaces, as well as simulation of micellar shape transitions.

Molecular dynamics simulation of SDS and CTAB micellization and prediction of partition equilibria with COSMOmic

Molecular dynamics (MD) simulations of the self-assembly of different ionic surfactants have been performed in order to obtain representative micellar structures. Subsequently, these structures were used to predict the partition behavior of various solutes in these micelles with COSMOmic, an extension of COSMO-RS. This paper includes multiple self-assembled micelles of SDS (sodium dodecyl sulfate, anionic surfactant) and CTAB (cetyltrimethylammoniumbromide, cationic surfactant) at different concentrations. Micellar size, density profiles, and shape (eccentricity) have been investigated. However, the size strongly depends on the functional definition of a micelle. For this reason, we present a method based on the free monomer concentration in aqueous solution as an optimization criterion for the micelle definition. The combination of MD with COSMOmic has the benefit of combining detailed atomistic information from MD with fast calculations of COSMOmic. For the first time the influence of micelle structure on pratition equilibria, predicted with COSMOmic, were investigated. In case of SDS more than 4600 and for CTAB more than 800 single micelles have been studied. The predictions of the partition coefficients with COSMOmic are in good agreement with experimental data. Additionally, the most favorable locations of selected molecules in the micelles as well as probable energy barriers are determined even for complex solutes such as toluene, propanolol, ephedrine, acetone, phenol, lidocaine, syringic acid, coumarin, isovanillin, ferulic acid, and vanillic acid. This method can therefore be applied as a potential screening tool for solutes (e.g., drugs) to find the optimal solute-surfactant combination.

Insight into α-synuclein plasticity and misfolding from differential micelle binding

Misfolded species of the 140-residue protein α-synuclein (αS) are implicated in the demise of dopaminergic neurons, resulting in fatal neurodegeneration. The intrinsically unstructured protein binds curved synaptic vesicle membranes in helical conformations but misfolds into amyloid fibrils via β-sheet interactions. Breaks in helical αS conformation may offer a pathway to transition from helical to sheet conformation. Here, we explore the evolution of broken αS helix conformations formed in complex with SDS and SLAS micelles by molecular dynamics simulations. The population distribution of experimentally observed αS conformations is related to the spatial concentration of intrinsic micelle shape perturbations. For the success of micelle-induced αS folding, we posit the length of the first helical segment formed, which controls micelle ellipticity, to be a key determinant. The degree of micelle curvature relates to the arrangement and segmental motions of helical secondary structure elements. A criterion for assessing the reproduction of such intermediate time scale protein dynamics is introduced by comparing the sampling of experimental and simulated spin label distributions. Finally, at the sites of breaks in the elongated, marginally stable αS helix, vulnerability to forming a transient, intramolecular β-sheet is identified. Upon subsequent intermolecular β-sheet pairing, pathological αS amyloid formation from initial helical conformation is thus achievable.

CHARMM-GUI micelle builder for pure/mixed micelle and protein/micelle complex systems

Micelle Builder in CHARMM-GUI, http://www.charmm-gui.org/input/micelle , is a web-based graphical user interface to build pure/mixed micelle and protein/micelle complex systems for molecular dynamics (MD) simulation. The robustness of Micelle Builder is tested by simulating four detergent-only homogeneous micelles of DHPC (dihexanoylphosphatidylcholine), DPC (dodecylphosphocholine), TPC (tetradecylphosphocholine), and SDS (sodium dodecyl sulfate) and comparing the calculated micelle properties with experiments and previous simulations. As a representative protein/micelle model, Pf1 coat protein is modeled and simulated in DHPC micelles with three different numbers of DHPC molecules. While the number of DHPC molecules in direct contact with Pf1 protein converges during the simulation, distinct behavior and geometry of micelles lead to different protein conformations in comparison to that in bilayers. It is our hope that CHARMM-GUI Micelle Builder can be used for simulation studies of various protein/micelle systems to better understand the protein structure and dynamics in micelles as well as distribution of detergents and their dynamics around proteins.

Investigation of Chiral Molecular Micelles by NMR Spectroscopy and Molecular Dynamics Simulation

NMR spectroscopy and molecular dynamics (MD) simulation analyses of the chiral molecular micelles poly-(sodium undecyl-(L,L)-leucine-valine) (poly-SULV) and poly-(sodium undecyl-(L,L)- valine-leucine) (poly-(SUVL)) are reported. Both molecular micelles are used as chiral selectors in electrokinetic chromatography and each consists of covalently linked surfactant chains with chiral dipeptide headgroups. To provide experimental support for the structures from MD simulations, NOESY spectra were used to identify protons in close spatial proximity. Results from the NOESY analyses were then compared to radial distribution functions from MD simulations. In addition, the hydrodynamic radii of both molecular micelles were calculated from NMR-derived diffusion coefficients. Corresponding radii from the MD simulations were found to be in agreement with these experimental results. NMR diffusion experiments were also used to measure association constants for polar and non-polar binaphthyl analytes binding to both molecular micelles. Poly(SUVL) was found to bind the non-polar analyte enantiomers more strongly, while the more polar analyte enantiomers interacted more strongly with poly(SULV). MD simulations in tum showed that poly(SUL V) had a more open structure that gave greater access for water molecules to the dipeptide headgroup region.

Flower micelle of amphiphilic random copolymers in aqueous media

The structure of the flower micelle formed by an amphiphilic random copolymer, sodium (2-acrylamido)-2-methylpropanesulfonate and N-dodecylmethacrylamide p(AMPS/C12), in 0.05 M aqueous NaCl was investigated by fully atomistic molecular dynamics simulation as well as by light scattering, and the results were compared with the flower micelle model of the minimum loop size, recently proposed by Kawata et al. [Macromolecules 2007, 40, 1174-1180]. After a sufficiently long simulation time, simulated p(AMPS/C12) chain with the degree of polymerization of 200 and C12 content of 50 mol % formed a unicore micelle, of which radius of gyration was much smaller than the AMPS homopolymer with the same degree of polymerization. The simulated micellar structure was analyzed in terms of density distribution functions for dodecyl groups, the main chain, and sulfonate groups as functions of the radial distance r from the center of mass of dodecyl groups. Only dodecyl groups exist at r less, similar 1.5 nm, and the main chain and sulfonate groups distribute in the range of r between 1.5 and 3.5 nm, but there were dodecyl groups coexisting with the main chain and sulfonate groups beyond r = 1.5 nm. All these structural features, as well as hydrodynamic radius data for p(AMPS/C12) with C12 contents higher than ca. 20 mol % obtained by light scattering, agreed with the predictions of the flower micelle model of the minimum loop size.

Phenyl-ring-bearing cationic surfactants: effect of ring location on the micellar structure

A series of isomeric cationic surfactants (S1-S5) bearing a long alkyl chain that carries a 1,4-phenylene unit and a trimethyl ammonium headgroup was synthesized; the location of the phenyl ring within the alkyl tail was varied in an effort to understand its influence on the amphiphilic properties of the surfactants. The cmc's of the surfactants were estimated using ionic conductivity measurements and isothermal calorimetric titrations (ITC); the values obtained by the two methods were found to be in excellent agreement. The ITC measurements provided additional insight into the various thermodynamic parameters associated with the micellization process. Although all five surfactants have exactly the same molecular formula, their micellar properties were seen to vary dramatically depending on the location of the phenyl ring; the cmc was seen to decrease by almost an order of magnitude when the phenyl ring was moved from the tail end (cmc of S1 is 23 mM) to the headgroup region (cmc of S5 is 3 mM). In all cases, the enthalpy of micellization was negative but the entropy of micellization was positive, suggesting that in all of these systems the formation of micelles is both enthalpically and entropically favored. As expected, the decrease in cmc values upon moving the phenyl ring from the tail end to the headgroup region is accompanied by an increase in the thermodynamic driving force (ΔG) for micellization. To understand further the differences in the micellar structure of these surfactants, small-angle neutron scattering (SANS) measurements were carried out; these measurements reveal that the aggregation number of the micelles increases as the cmc decreases. This increase in the aggregation number is also accompanied by an increase in the asphericity of the micellar aggregate and a decrease in the fractional charge. Geometric packing arguments are presented to account for these changes in aggregation behavior as a function of phenyl ring location.

Explanation of experimental results of mixed micelles of homologous surfactants through a MM2 bidimensional modeling

A computational modeling (in gas phase) to study the disposition of the homologous surfactants in a bidimensional simple model of mixed and homogeneous micelles was performed for the case of R-trimethylammonium bromide surfactants with different linear R lengths from R = C(5) to C(17). First, the bidimensional homogeneous (one component) micelle was modeled, and as a second step, heterogeneous (two components) bidimensional micelles were modeled. The difference in the number of carbon atoms between hydrocarbon chains of the surfactants in the heterogeneous micelles, Δn(C), ranged from 2 to 8. Results were contrasted with experimental data obtained at our own laboratory. The exothermic values of the steric energy changes showed strong attraction between components of homologous surfactants mixture, especially when one of the surfactants has a long chain. It may be argued that the inclusion of a shorter surfactant in the mixture and the twisting of the longer surfactant makes the bidimensional arrangement formation more exothermic. All predictions were in agreement with previous experimental results.

Coarse-grained molecular dynamics simulation of the aggregation properties of multiheaded cationic surfactants in water

The aggregation property of multiheaded surfactants has been investigated by constant pressure molecular dynamics (MD) simulation in aqueous medium. The model multiheaded surfactants contain more than one headgroup (x = 2, 3, and 4) for a single tail group. This increases the hydrophilic charge progressively over the hydrophobic tail which has dramatic consequences in the aggregation behavior. In particular, we have looked at the change in the aggregation property such as critical micellar concentration (cmc), aggregation number, and size of the micelles for the multiheaded surfactants in water. We find with increasing number of headgroups of the multiheaded surfactants that the cmc values increase and the aggregation numbers as well as the size of the micelles decrease. These trends are in agreement with the experimental findings as reported earlier with x = 1, 2, and 3. We also predict the aggregation properties of multiheaded surfactant with four headgroups (x = 4) for which no experimental studies exist yet.

Synthesis and micellar properties of surface-active ionic liquids: 1-Alkyl-3-methylimidazolium chlorides

A series of surface-active ionic liquids, RMeImCl, has been synthesized by the reaction of purified 1-methylimidazole and 1-chloroalkanes, RCl, R=C(10),C(12),C(14), and C(16), respectively. Adsorption and aggregation of these surfactants in water have been studied by surface tension measurement. Additionally, solution conductivity, electromotive force, fluorescence quenching of micelle-solubilized pyrene, and static light scattering have been employed to investigate micelle formation. The following changes resulted from an increase in the length of R: an increase of micelle aggregation number; a decrease of: minimum area/surfactant molecule at solution/air interface; critical micelle concentration, and degree of counter-ion dissociation. Theoretically-calculated aggregation numbers and those based on quenching of pyrene are in good agreement. Gibbs free energies of adsorption at solution/air interface, DeltaG(ads)(0), and micelle formation in water, DeltaG(mic)(0), were calculated, and compared to those of three surfactant series, alkylpyridinium chlorides, RPyCl, alkylbenzyldimethylammonium chlorides, RBzMe(2)Cl, and benzyl(3-acylaminoethyl)dimethylammonium chlorides, R(')AEtBzMe(2)Cl, respectively. Contributions to the above-mentioned Gibbs free energies from surfactant methylene groups (in the hydrophobic tail) and the head-group were calculated. For RMeImCl, the former energy is similar to that of other cationic surfactants. The corresponding free energy contribution of the head-group to DeltaG(mic)(0) showed the following order: RPyCl approximately RBzMe(2)Cl>RMeImCl>R(')AEtBzMe(2)Cl. The head-groups of the first two surfactant series are more hydrophobic than the imidazolium ring of RMeImCl, this should favor their aggregation. Micellization of RMeImCl, however, is driven by a relatively strong hydrogen-bonding between the chloride ion and the hydrogens in the imidazolium ring, in particular the relatively acidic H2. This interaction more than compensates for the relative hydrophilic character of the diazolium ring. As indicated by the corresponding DeltaG(mic)(0), micellization of R(')AEtBzMe(2)Cl is more favorable than that of RMeImCl because the CONH group of the former surfactant series forms hydrogen bonds to both the counter-ion and the neighboring molecules in the micelle.