Molecular Modeling of Triton X Micelles: Force Field Parameters, Self-Assembly, and Partition Equilibria

Molecular Dynamic Simulations of Aqueous Micellar Organometallic Catalysis: Methane Functionalization as a Case Study

Molecular Dynamics (MD) simulations constitute a powerful tool that provides a 3D perspective of the dynamical behavior of chemical systems. Herein the first MD study of the dynamics of a catalytic organometallic system, in micellar media, is presented. The challenging methane catalytic functionalization into ethyl propionate through a silver-catalyzed process has been targeted as the case study. The intimate nature of the micelles formed with the surfactants sodium dodecylsulfate (SDS) and potassium perfluorooctane sulfonate (PFOS) has been ascertained, as well as the relative distribution of the main actors in this transformation, namely methane, the diazo reagent and the silver catalyst, the latter in two different forms: the initial compound and a silver-carbene intermediate. Catalyst deactivation occurs with halide containing surfactants dodecyltrimethylammonium chloride (DTAC) and Triton X-100. Computed simulations allow explaining the experimental results, indicating that micelles behave differently regarding the degree of accumulation and the local distribution of the reactants and their effect in the molecular collisions leading to net reaction.

Collection of Partition Coefficients in Hexadecyltrimethylammonium Bromide, Sodium Cholate, and Lithium Perfluorooctanesulfonate Micellar Solutions: Experimental Determination and Computational Predictions

This study focuses on determining the partition coefficients (logP) of a diverse set of 63 molecules in three distinct micellar systems: hexadecyltrimethylammonium bromide (HTAB), sodium cholate (SC), and lithium perfluorooctanesulfonate (LPFOS). The experimental log p values were obtained through micellar electrokinetic chromatography (MEKC) experiments, conducted under controlled pH conditions. Then, Quantum Mechanics (QM) and machine learning approaches are proposed for the prediction of the partition coefficients in these three micellar systems. In the applied QM approach, the experimentally obtained partition coefficients were correlated with the calculated values for the case of the 15 solvent mixtures. Using Density Function Theory (DFT) with the B3LYP functional, we calculated the solvation free energies of 63 molecules in these 16 solvents. The combined data from the experimental partition coefficients in the three micellar formulations showed that the 1-propanol/water combination demonstrated the best agreement with the experimental partition coefficients for the SC and HTAB micelles. Moreover, we employed the SVM approach and k-means clustering based on the generation of the chemical descriptor space. The analysis revealed distinct partitioning patterns associated with specific characteristic features within each identified class. These results indicate the utility of the combined techniques when we want an efficient and quicker model for predicting partition coefficients in diverse micelles.

Highly-Selective Harvesting of (6,4) SWCNTs Using the Aqueous Two-Phase Extraction Method and Nonionic Surfactants

Monochiral single-walled carbon nanotubes (SWCNTs) are indispensable for advancing the technology readiness level of nanocarbon-based concepts. In recent times, many separation techniques have been developed to obtain specific SWCNTs from raw unsorted materials to catalyze the development in this area. This work presents how the aqueous two-phase extraction (ATPE) method can be enhanced for the straightforward isolation of (6,4) SWCNTs in one step. Introducing nonionic surfactant into the typically employed mixture of anionic surfactants, which drive the partitioning, is essential to increasing the ATPE system's resolution. A thorough analysis of the parameter space by experiments and modeling reveals the underlying interactions between SWCNTs, surfactants, and phase-forming agents, which drive the partitioning. Based on new insight gained on this front, a separation mechanism is proposed. Notably, the developed method is highly robust, which is proven by isolating (6,4) SWCNTs from several raw SWCNT materials, including SWCNT waste generated over the years in the laboratory.

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.

Partition of the Reactive Species of the Suzuki-Miyaura Reaction between Aqueous and Micellar Environments

The Suzuki-Miyaura reaction between the aryl halide (1) and the phenyl boronic acid (2), in the presence of the palladium(0) complex (3) as catalyst, gives the cross-coupling product (4) in quantitative yield when performed in basic aqueous solution of the nonionic surfactant Kolliphor-EL (K-EL). The partition between the aqueous and micellar environments of the species of this reaction has been investigated by means of Molecular Dynamics (MD) simulations. Starting from the K-EL molecules dispersed in water, a micelle model has been generated by MD simulations, adopting the 2016H66 force field. Reagent and product species have been described with the same force field, once the reliability of this force field has been tested comparing the n-octanol/water partition free energies calculated from the MD and Free Energy Perturbation (FEP) method with those obtained from the quantum-mechanical SMD method. The potential of mean force for the transfer process between water and the micellar phase of the different species has been calculated by the MD simulations and the Umbrella Sampling (US) method. The overall picture that emerges from these results confirms that the molecular species involved in this reaction prefers the micellar environment and concentrates in different but close zones of the micelle. This supports the experimental evidence that the use of suitable surfactant agents promotes reactivity, allowing micelles to behave as nanoreactors in which reactive species are solubilized and enhance their local concentration.

Impact of drug aggregation on the structural and dynamic properties of Triton X-100 micelles

Surfactants are used in a wide range of chemical and biological applications, and for pharmaceutical purposes are frequently employed to enhance the solubility of poorly water soluble drugs. In this study, all-atom molecular dynamics (MD) simulations and small-angle neutron scattering (SANS) experiments have been used to investigate the drug solubilisation capabilities of the micelles that result from 10 wt% aqueous solutions of the non-ionic surfactant, Triton X-100 (TX-100). Specifically, we have investigated the solubilisation of saturation amounts of the sodium salts of two nonsteroidal anti-inflammatory drugs: ibuprofen and indomethacin. We find that the ibuprofen-loaded micelles are more non-spherical than the indomethacin-loaded micelles which are in turn even more non-spherical than the TX-100 micelles that form in the absence of any drug. Our simulations show that the TX-100 micelles are able to solubilise twice as many indomethacin molecules as ibuprofen molecules, and the indomethacin molecules form larger aggregates in the core of the micelle than ibuprofen. These large indomethacin aggregates result in the destabilisation of the TX-100 micelle, which leads to an increase in the amount of water inside of the core of the micelle. These combined effects cause the eventual division of the indomethacin-loaded micelle into two daughter micelles. These results provide a mechanistic description of how drug interactions can affect the stability of the resulting nanoparticles.

Structural Insights into Self-Assembled Aerosol-OT Aggregates in Aqueous Media Using Atomistic Molecular Dynamics

In water, the surfactant dioctyl sulfosuccinate (Aerosol-OT or AOT) exhibits diverse aggregate structures, ranging from micelles to lamella. An atomic-level understanding, however, of the formation and structure of these aggregates is lacking. Herein, using atomistic molecular dynamics (MD) with microsecond-long simulations, self-assembly of AOT in water is studied for concentrations of 1, 7.2, and 20 wt % at 293 K and for 7.2 wt % at 353 K. Assembly proceeds through stepwise association and dissociation of single AOT molecules, and the fusion and fission of AOT clusters. At 293 K, AOT self-assembles into either (i) spherical micelles (1 wt %), (ii) biphasic systems consisting of rod-like and prolate spheroidal micelles (7.2 wt %), or (iii) bilayers (20 wt %). We hypothesize that the observed rod-like structure is a precursor to lamellar microdomains found experimentally in biphasic dispersions. Increasing temperature to 353 K at 7.2 wt % results in a system consisting of prolate micelles but no rod-like micelles. Simulated phase behavior agrees with previously published experimental observations. Individual aggregates formed during self-assembly are identified using graph theory. Structural metrics of these aggregates like the radius of gyration, shape anisotropy, and prolateness are presented. Trends in structural metrics quantitatively reflect how shapes and sizes of AOT aggregates vary with surfactant concentration and temperature. These simulations provide deeper insight into open questions in the scientific community and demonstrate a method to generate physics-based micelle structures that can be used to rationalize experimental observations.

COSMO-RS Analysis of CO2 Solubility in N-Methyldiethanolamine, Sulfolane, and 1-Butyl-3-methyl-imidazolium Acetate Activated by 2-Methylpiperazine for Postcombustion Carbon Capture

Novel aqueous (aq) blends of N-methyldiethanolamine (MDEA), sulfolane (TMSO₂), and 1-butyl-3-methyl-imidazolium acetate ([bmim][Ac]) with amine activator 2-methylpiperazine (2-MPZ) are analyzed through conductor-like screening model for real solvents (COSMO-RS) for possible application in the chemisorption of CO₂. The molecules associated are analyzed for their ground-state energy, σ potential, and σ surface. Thermodynamic and physicochemical properties have been assessed and paralleled with the experimental data. Vapor pressure of the blended systems and pure component density and viscosity have been compared successfully with the experimental data. Important binary interaction parameters for the aqueous blends over a wide temperature, pressure, and concentration range have been estimated for NRTL, WILSON, and UNIQUAC 4 models. The COSMO-RS theory is further applied in calculating the expected CO₂ solubility over a pressure range of 1.0-3.0 bar and temperature range of 303.15-323.15 K. Henry's constant and free energy of solvation to realize the physical absorption through intermolecular interaction offered by the proposed solvents. Perceptive molecular learning from the behavior of chemical constituents involved indicated that the best suitable solvent is aq (MDEA + 2-MPZ).

Accessing lipophilicity of drugs with biomimetic models: A comparative study using liposomes and micelles

Lipophilicity is a physicochemical property of crucial importance in drug discovery and drug design. Biomimetic models, such as liposomes and micelles, constitute a valuable tool for the assessment of lipophilicity through the determination of partition coefficients (log K_p). However, the lack of standardization hampers the judgment about which model or method has the best and broadest passive drug permeation predictive capacity. This work provides a comparative analysis between the methodologies based on biomimetic models to determine the partition coefficient (log K_p). For that purpose, a set of reference substances preconized by the Organization for Economic Cooperation and Development (OECD) guidelines was used. The biomimetic models employed were liposomes and micelles composed by 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) and hexadecylphosphocholine (HePC), respectively. Both lipids were used as representative phospholipids of natural membranes. The partition coefficients between biomimetic models and aqueous phases were determined by derivative spectroscopy at physiological conditions (37 °C and pH 7.4). The partition coefficients obtained using biomimetic models are quite different and more reliable than the ones obtained using an octanol/water system. Comparing the performance of the two biomimetic models, micelles revealed to be suitable only for substances with high molar absorption coefficient and log K_p > 3, but in general liposomes are the best model for accessing lipophilicity of drugs. Furthermore, a comparison between experimental data and the partition coefficients determined by the computational method COSMOmic is also provided and discussed. As a final summarizing result, a decision tree is provided in order to guide the selection of a tool for assessing the lipophilicity of drugs.

Micellization and Partition Equilibria in Mixed Nonionic/Ionic Micellar Systems: Predictions with Molecular Models

In practical applications, surfactant solutions are mostly used in mixtures of nonionic and ionic surfactants because they have improved characteristics compared to those of single surfactant solutions. By adjusting the composition of the micelles and the pH value, the solubilization of solutes can be enhanced. Nevertheless, the partitioning of solutes between nonionic/ionic mixed micelles and the aqueous phase is studied to a much lesser extent than for single surfactant solutions. Theoretical methods to predict partition equilibria in mixed micelles are of interest for screening studies. For those, the composition of the mixed micelle has to be known. Here we investigate mixtures of TX-114 (Triton X-114), Brij35 (C12E23), SDS (sodium dodecyl sulfate), and CTAB (cetyltrimethylammonium bromide). First, to investigate the surfactant compositions in the micelles, molecular dynamics (MD) self-assembly simulations were applied. Thereafter, the predictive COSMO-RS model, which applies the pseudophase approach, and its extension to anisotropic systems termed COSMOmic were compared for the prediction of partition equilibria in mixed micelles, where various molar ratios of the surfactants were considered. It could be demonstrated that both methods are applicable and lead to reasonable predictions for neutral molecules. However, taking into account the three-dimensional structure of the micelle is beneficial because the calculations with COSMOmic are in better agreement with experimental results. Because the partitioning behavior of ionizable molecules in mixed micelles is of particular interest, the partitioning of ionized isovanillin in mixed Brij35/CTAB micelles at different micelle compositions was calculated with COSMOmic. Using a thermodynamic cycle, the position-dependent pK_a of isovanillin within the micelle is calculated on the basis of COSMOmic free energy profiles. As a result, the protolytic equilibrium of isovanillin within the micelles can be taken into account, which is crucial for the reliable prediction of partition coefficients.

MARTINI Coarse-Grained Model of Triton TX-100 in Pure DPPC Monolayer and Bilayer Interfaces

The coarse-grained MARTINI model of Triton TX-100 has been validated by direct comparison of the experimental and calculated area increase in pure DPPC lipid bilayers and monolayers at water/air interfaces in the presence of surfactant and by comparison of electron density profiles calculated with more detailed atomistic models based on the CHARMM force field. Bilayer simulations have been performed and compared with monolayers and with atomistic models. The validated CG model has been employed to study the phase separation of TX-100 molecules in lipid bilayers and the effect of the lipid bilayer curvature.