Quantifying the Hydrophobic Effect. 1. A Computer Simulation−Molecular-Thermodynamic Model for the Self-Assembly of Hydrophobic and Amphiphilic Solutes in Aqueous Solution

Determining interfacial tension and critical micelle concentrations of surfactants from atomistic molecular simulations

Hypothesis Atomistically-detailed models of surfactants provide quantitative information on the molecular interactions and spatial distributions at fluid interfaces. Hence, it should be possible to extract from this information, macroscopical thermophysical properties such as interfacial tension, critical micelle concentrations and the relationship between these properties and the bulk fluid surfactant concentrations. Simulations and Experiments Molecular-scale interfacial of systems containing n-dodecyl β-glucoside (APG12) are simulated using classical molecular dynamics. The bulk phases and the corresponding interfacial regions are all explicitly detailed using an all-atom force field (PCFF+). During the simulation, the behaviour of the interface is analyzed geometrically to obtain an approximated value of the critical micelle concentration (CMC) in terms of the surfactant area number density and the interfacial tension is assessed through the analysis of the forces amongst molecules. New experimental determinations are reported for the surface tension of APG12 at the water/air and at the water/n-decane interfaces. Findings We showcase the application of a thermodynamic framework that inter-relates interfacial tensions, surface densities, CMCs and bulk surfactant concentrations, which allows the in silico quantitative prediction of interfacial tension isotherms.

Perspectives on the Lindman Hypothesis and Cellulose Interactions

In the history of cellulose chemistry, hydrogen bonding has been the predominant explanation when discussing intermolecular interactions between cellulose polymers. This is the general consensus in scholarly textbooks and in many research articles, and it applies to several other biomacromolecules' interactions as well. This rather unbalanced description of cellulose has likely impacted the development of materials based on the processing of cellulose-for example, via dissolution in various solvent systems and regeneration into solid materials, such as films and fibers, and even traditional wood fiber handling and papermaking. In this review, we take as a starting point the questioning of the general description of the nature of cellulose and cellulose interactions initiated by Professor Björn Lindman, based on generic physicochemical reasoning about surfactants and polymers. This dispute, which became known as "the Lindman hypothesis", highlights the importance of hydrophobic interactions in cellulose systems and that cellulose is an amphiphilic polymer. This paper elaborates on Björn Lindman's contribution to the subject, which has caused the scientific community to revisit cellulose and reconsider certain phenomena from other perspectives.

Structural features of interfacial water predict the hydrophobicity of chemically heterogeneous surfaces

The hydrophobicity of an interface determines the magnitude of hydrophobic interactions that drive numerous biological and industrial processes. Chemically heterogeneous interfaces are abundant in these contexts; examples include the surfaces of proteins, functionalized nanomaterials, and polymeric materials. While the hydrophobicity of nonpolar solutes can be predicted and related to the structure of interfacial water molecules, predicting the hydrophobicity of chemically heterogeneous interfaces remains a challenge because of the complex, non-additive contributions to hydrophobicity that depend on the chemical identity and nanoscale spatial arrangements of polar and nonpolar groups. In this work, we utilize atomistic molecular dynamics simulations in conjunction with enhanced sampling and data-centric analysis techniques to quantitatively relate changes in interfacial water structure to the hydration free energy (a thermodynamically well-defined descriptor of hydrophobicity) of chemically heterogeneous interfaces. We analyze a large data set of 58 self-assembled monolayers (SAMs) composed of ligands with nonpolar and polar end groups of different chemical identity (amine, amide, and hydroxyl) in five mole fractions, two spatial patterns, and with scaled partial charges. We find that only five features of interfacial water structure are required to accurately predict hydration free energies. Examination of these features reveals mechanistic insights into the interfacial hydrogen bonding behaviors that distinguish different surface compositions and patterns. This analysis also identifies the probability of highly coordinated water structures as a unique signature of hydrophobicity. These insights provide a physical basis to understand the hydrophobicity of chemically heterogeneous interfaces and connect hydrophobicity to experimentally accessible perturbations of interfacial water structure.

Light scattering from mixtures of interacting, nonionic micelles with hydrophobic solutes

Model equations for the Rayleigh ratio and the electric field autocorrelation function are derived using thermodynamic fluctuation theory applied to crowded solute-containing micellar solutions and microemulsions with negligible molecular species and polydispersity. This theory invokes non-equilibrium thermodynamics and enforces local equilibrium between molecular solute, surfactant, and the various micellar species, in order to elucidate the influence of self-assembly on light scattering correlation functions. We find that self-assembly driven variations in the average micelle radius and aggregation number along gradients in concentration, which were previously shown to drive strong multicomponent diffusion effects expressed via the ternary diffusivity matrix [D], do not affect the scattering functions in the limit of zero local polydispersity. Hence, theoretical predictions for the Rayleigh ratio and the field autocorrelation function for ternary mixtures of solute-containing, locally monodisperse micellar solutions are identical to those developed for binary mixtures of monodisperse, colloidal hard spheres. However, self-assembly driven multicomponent diffusion phenomena are predicted to influence the thermodynamic driving forces for diffusion in these mixtures. In support of our theoretical results, measurements for the Rayleigh ratio and the field autocorrelation function for ternary aqueous solutions of decaethylene glycol monododecyl ether (C₁₂E₁₀) with either decane or limonene solute were performed for several molar ratios and volume fractions up to ϕ ≈ 0.25, and for binary mixtures of C₁₂E₁₀/water up to ϕ ≈ 0.5. Excellent agreement between our light scattering theory and experimental data is achieved for low to moderate volume fractions (ϕ < 0.3), and at higher concentrations when our theoretical results are corrected to account for micelle dehydration.

Prediction of Micellar Thermodynamics of Nonionic Surfactants Based on the Square Gradient Theory

A geometry-dependent contribution based on the square gradient theory of van der Waals is proposed as a predictive modification of the interfacial energy contribution for the micellar thermodynamic theory. The model has an analytic prediction for the spherical and cylindrical geometries. For ellipsoidal geometry, a simple yet physically meaningful approximation is proposed. The critical micelle concentration (CMC) and the surface tension isotherm under the new contribution are compared with the classical theory. The modified model describes qualitatively the available experimental data and the surface isotherm, showing an improvement in the predictions of the CMC.

Preparation of Bioderived and Biodegradable Surfactants Based on an Intrinsically Disordered Protein Sequence

Surfactants, block copolymers, and other types of micellar systems are used in a wide variety of biomedical and industrial processes. However, most commonly used surfactants are synthetically derived and pose environmental and toxicological concerns throughout their product life cycle. Because of this, bioderived and biodegradable surfactants are promising alternatives. For biosurfactants to be implemented industrially, they need to be produced on a large scale and also have tailorable properties that match those afforded by the polymerization of synthetic surfactants. In this paper, a scalable and versatile production method for biosurfactants based on a hydrophilic intrinsically disordered protein (IDP) sequence with a genetically engineered hydrophobic domain is used to study variables that impact their physicochemical and self-assembling properties. These amphiphilic sequences were found to self-assemble into micelles over a broad range of temperatures, pH values, and ionic strengths. To investigate the role of the IDP hydrophilic domain on self-assembly, variants with increased overall charges and systematically decreased IDP domain lengths were produced and examined for their sizes, morphologies, and critical micelle concentrations (CMCs). The results of these studies indicate that decreasing the length of the IDP domain and consequently the molecular weight and hydrophilic fraction leads to smaller micelles. In addition, significantly increasing the amount of charged residues in the hydrophilic IDP domain results in micelles of similar sizes but with higher CMC values. This represents an initial step in developing a quantitative model for the future engineering of biosurfactants based on this IDP sequence.

Identifying nonadditive contributions to the hydrophobicity of chemically heterogeneous surfaces via dual-loop active learning

Hydrophobic interactions drive numerous biological and synthetic processes. The materials used in these processes often possess chemically heterogeneous surfaces that are characterized by diverse chemical groups positioned in close proximity at the nanoscale; examples include functionalized nanomaterials and biomolecules, such as proteins and peptides. Nonadditive contributions to the hydrophobicity of such surfaces depend on the chemical identities and spatial patterns of polar and nonpolar groups in ways that remain poorly understood. Here, we develop a dual-loop active learning framework that combines a fast reduced-accuracy method (a convolutional neural network) with a slow higher-accuracy method (molecular dynamics simulations with enhanced sampling) to efficiently predict the hydration free energy, a thermodynamic descriptor of hydrophobicity, for nearly 200 000 chemically heterogeneous self-assembled monolayers (SAMs). Analysis of this dataset reveals that SAMs with distinct polar groups exhibit substantial variations in hydrophobicity as a function of their composition and patterning, but the clustering of nonpolar groups is a common signature of highly hydrophobic patterns. Further molecular dynamics analysis relates such clustering to the perturbation of interfacial water structure. These results provide new insight into the influence of chemical heterogeneity on hydrophobicity via quantitative analysis of a large set of surfaces, enabled by the active learning approach.

Hydrophobic interactions control the self-assembly of DNA and cellulose

Desoxyribosenucleic acid, DNA, and cellulose molecules self-assemble in aqueous systems. This aggregation is the basis of the important functions of these biological macromolecules. Both DNA and cellulose have significant polar and nonpolar parts and there is a delicate balance between hydrophilic and hydrophobic interactions. The hydrophilic interactions related to net charges have been thoroughly studied and are well understood. On the other hand, the detailed roles of hydrogen bonding and hydrophobic interactions have remained controversial. It is found that the contributions of hydrophobic interactions in driving important processes, like the double-helix formation of DNA and the aqueous dissolution of cellulose, are dominating whereas the net contribution from hydrogen bonding is small. In reviewing the roles of different interactions for DNA and cellulose it is useful to compare with the self-assembly features of surfactants, the simplest case of amphiphilic molecules. Pertinent information on the amphiphilic character of cellulose and DNA can be obtained from the association with surfactants, as well as on modifying the hydrophobic interactions by additives.

Stabilization of aqueous graphene dispersions utilizing a biocompatible dispersant: a molecular dynamics study

Investigation of the high efficiency of flavin mononucleotide sodium salt (FMNS) for the stabilization of aqueous graphene dispersions using all-atom molecular dynamics simulations.

Growth of wormlike micelles in nonionic surfactant solutions: Quantitative theory vs. experiment

Despite the considerable advances of molecular-thermodynamic theory of micelle growth, agreement between theory and experiment has been achieved only in isolated cases. A general theory that can provide self-consistent quantitative description of the growth of wormlike micelles in mixed surfactant solutions, including the experimentally observed high peaks in viscosity and aggregation number, is still missing. As a step toward the creation of such theory, here we consider the simplest system - nonionic wormlike surfactant micelles from polyoxyethylene alkyl ethers, C_iE_j. Our goal is to construct a molecular-thermodynamic model that is in agreement with the available experimental data. For this goal, we systematized data for the micelle mean mass aggregation number, from which the micelle growth parameter was determined at various temperatures. None of the available models can give a quantitative description of these data. We constructed a new model, which is based on theoretical expressions for the interfacial-tension, headgroup-steric and chain-conformation components of micelle free energy, along with appropriate expressions for the parameters of the model, including their temperature and curvature dependencies. Special attention was paid to the surfactant chain-conformation free energy, for which a new more general formula was derived. As a result, relatively simple theoretical expressions are obtained. All parameters that enter these expressions are known, which facilitates the theoretical modeling of micelle growth for various nonionic surfactants in excellent agreement with the experiment. The constructed model can serve as a basis that can be further upgraded to obtain quantitative description of micelle growth in more complicated systems, including binary and ternary mixtures of nonionic, ionic and zwitterionic surfactants, which determines the viscosity and stability of various formulations in personal-care and house-hold detergency.

Melatonin-directed micellization: a case for tryptophan metabolites and their classical bioisosteres as templates for the self-assembly of bipyridinium-based supramolecular amphiphiles in water

The bulk solution properties of amphiphilic formulations are derivative of their self-assembly into higher ordered supramolecular assemblies known as micelles and of their ordering at the air-water interface. Exerting control over the surface-active properties of amphiphiles and their propensity to aggregate in pure water is most often fine-tuned by covalent modification of their molecular structure. Nevertheless structural constraints which limit the performance of amphiphiles do emerge when trying to develop more sophisticated systems which undergo for example, shape-defined controlled assembly and/or respond to external stimuli. In this regard, the template-modulated assembly of the so-called "supramolecular amphiphiles" continues to make progress ordering molecules that otherwise have very little to no driving force to aggregate in a prescribed manner in aqueous solutions. Herein we describe the template-modulated micellization and ordering at the air-water interface of bipyridinium-based supramolecular amphiphiles triggered by host-guest interactions with high specificity for the neurotransmitter melatonin over its biosynthetic synthon l-tryptophan and the thermodynamic parameters governing the template-modulated micellization process. When bound to the bipyridinium units of micellized surfactant molecules, melatonin effectively serves as "molecular glue" capable of lowering the CMC by 52% as compared to untemplated solutions. Analysis of this system suggests that a hallmark of donor-acceptor template-modulated micellization in water is a strong positively correlated temperature dependence of the CMC and the absence of a U-shaped CMC-temperature curve. Our findings make a case for the incorporation of l-tryptophan-based metabolites and their classical synthetic pharmaceutical bioisosteres as potential targets/components of donor-acceptor CT-based supramolecular amphiphile systems/materials operating in water.

Ion dehydration controls adsorption at the micellar interface: hydrotropic ions

The properties of ionic micelles depend on the nature of the counterion, and these effects become more evident as the ion adsorption at the interface increases. Prediction of the relative extent of ion adsorption is required for rational design of ionic micellar aggregates. Unlike the well understood adsorption of monatomic ions, the adsorption of polyatomic ions is not easily predicted. We combined experimental and computational methods to evaluate the affinity of hydrotropic ions, i.e., ions with polar and apolar regions, to the surface of positively charged micelles. We analyzed cationic micelles of dodecyltrimethylammonium and six hydrotropic counterions: methanesulfonate, trifluoromethanesulfonate, benzenesulfonate, acetate, trifluoroacetate and benzoate. Our results demonstrated that the apolar region of hydrotropic ions had the largest influence on micellar properties. The dehydration of the apolar region of hydrotropic ions upon their adsorption at the micellar interface determined the ion adsorption extension, differently to what was expected based on Collins' law of matching affinities. These results may lead to more general models to describe the adsorption of ions, including polyatomic ions, at the micellar interface.

Driving Force for Spontaneous Perforation of Bilayers Formed by Ionic Amphiphiles in Aqueous Salt

Spontaneous perforation of amphiphilic membranes is important in both living matter and technology because of an impact on functions of biological membranes and shape transitions of self-assembling structures. Nevertheless, no definite molecular mechanism has been established so far even for simple ionic surfactant systems. We show that spontaneous perforation of a bilayer formed by an ionic amphiphile is driven by electrostatics. Creation of large pores with a concave-convex geometry of the rim is promoted by lower electrostatic free energy than that for a flat nonperforated bilayer. The opposite effect comes from the elasticity of the hydrocarbon tails of the amphiphile that prefer flat geometry of a nonperforated bilayer. The balance between electrostatics and tail deformation controls the appearance of pores; this balance is modulated by added salt that screens the electrostatic interactions. We illustrate the proposed mechanism with the aid of classical aggregation model that has been extended by including an analytical description of the electrostatic contribution for the toroidal rim of a pore. Numerical solution of the linearized Poisson-Boltzmann equation confirms the role of electrostatic forces in formation of pores. For the ionic surfactants of C_nTAB family, we predict shape transitions including bilayer perforations and formation of branched micellar networks induced by changing salinity or temperature and demonstrate the effect of surfactant's molecular parameters on these transitions.

Combined Molecular Dynamics Simulation-Molecular-Thermodynamic Theory Framework for Predicting Surface Tensions

A molecular modeling approach is presented with a focus on quantitative predictions of the surface tension of aqueous surfactant solutions. The approach combines classical Molecular Dynamics (MD) simulations with a molecular-thermodynamic theory (MTT) [ Y. J. Nikas, S. Puvvada, D. Blankschtein, Langmuir 1992 , 8 , 2680 ]. The MD component is used to calculate thermodynamic and molecular parameters that are needed in the MTT model to determine the surface tension isotherm. The MD/MTT approach provides the important link between the surfactant bulk concentration, the experimental control parameter, and the surfactant surface concentration, the MD control parameter. We demonstrate the capability of the MD/MTT modeling approach on nonionic alkyl polyethylene glycol surfactants at the air-water interface and observe reasonable agreement of the predicted surface tensions and the experimental surface tension data over a wide range of surfactant concentrations below the critical micelle concentration. Our modeling approach can be extended to ionic surfactants and their mixtures with both ionic and nonionic surfactants at liquid-liquid interfaces.

Complex Formations between Surfactants and Polyelectrolytes of the Same Charge on a Water Surface

The mechanism of complex formation between surfactants and polyelectrolytes with the same charge on the water surface was investigated using molecular dynamics simulations and phase-sensitive sum-frequency generation vibrational spectroscopy. Although complex formation between highly charged surfactants and polyelectrolytes of the same charge is generally expected to be prohibited by the electrostatic repulsive force, our study shows that it is possible to form thermodynamically stable complexes when excess ions are present in the solution. We found that anionic partially hydrolyzed polyacrylamide (HPAM) could interact with anionic sodium dodecyl sulfate (SDS) on a water surface in the presence of salts. With excess Na⁺ ions in the solution, the charge screening effect allows HPAM to weakly interact with SDS via hydrogen bonds. In the presence of divalent Ca²⁺ ions, the surfactant and the polymer are strongly coupled by forming Ca²⁺ ion bridges and hydrogen bonds. Our calculation shows that the presence of Ca²⁺ ions creates a steep binding energy of ∼30 kJ/mol near the water surface. These results were qualitatively verified using phase-sensitive sum-frequency generation vibrational spectroscopy.

The influence of poly(ethylene glycol) on the micelle formation of alkyl maltosides used in membrane protein crystallization

The influence of poly(ethylene glycol) on the micelle formation of alkyl maltosides under conditions of membrane protein crystallization is investigated.

Modeling of the effects of ion specificity on the onset and growth of ionic micelles in a solution of simple salts

A new version of the molecular thermodynamic model has been developed that takes into account the effect of ion specificity on the free energy of aggregation. The specificity of salt is reflected by differences in the bare ionic sizes and polarizabilities leading to the difference in the dispersion interaction of ions with the aggregate. The model also contains parameters that characterize the compactness of ionic pairs formed between a mobile ion and surfactant's headgroup. The values of these parameters show that more chaotropic heads form tighter pairs with chaotropic ions whereas more cosmotropic heads form more compact pairs with cosmotropic ions. The formation of compact pairs in the micelle corona diminishes the preferable curvature of the aggregates and promotes their growth. The model has been applied to aqueous solutions of cationic (alkyltrimethylammonium, alkyldimethylammonium, and alkylpyridinium) and anionic (alkylsulfate and alkylcarboxylate) surfactants in the presence of simple 1:1 salts. With a single set of parameter values, the model reproduces the critical micelle concentration-salinity curves and the sphere-to-rod transitions or the absence of thereof and describes the aggregate growth for different simple salts, in good agreement with experiment.

Study of the microcharacter of ultrastable aqueous foam stabilized by a kind of flexible connecting bipolar-headed surfactant with existence of magnesium ion

In this paper, ultrastable aqueous foam stabilized by a kind of flexible connecting bipolar-headed surfactant alkyl polyoxyethylene sulfate (AE3S) with coexisting Mg(2+) was reported. Detailed molecular behaviors of AE3S in foam film with coexisting divalent cationic Ca(2+) or Mg(2+) were investigated by molecular dynamic simulation, comparing with the traditional surfactant sodium dodecyl sulfate (SDS), to find out how the microcharacter and array behavior of molecules in the foam film determined by molecular interaction effect the foam stability. It was found that the ultrastable foam film obtained by the cooperation of magnesium ions and AE3S was driven from two aspects: one is the favorable arrangement of surfactant molecules, and the other is the increase of capacity of foam films for resolutely holding water molecules deduced by a dipolar pair formed by the flexible connecting head groups of AE3S and hydrated Mg(2+) via intermolecular coactions, both related to the presence of magnesium ions. Foam lamella stability measurement and foam decay method were both used to evaluate the stability of foam. Fourier transform infrared (FT-IR) was used to detect the composition variation of foam film in the drainage process; the vibration peak of OH for water molecule shifted from the 3390 cm(-1) (being assigned to the bulk water integrated by hydrogen bonds) to 3685 cm(-1) (being assigned to the vibration of isolated water molecules) for the ultrastable foam film after complete drainage, which agreed very well with the molecular simulation results.

Molecular dynamics study of the effect of calcium ions on the monolayer of SDC and SDSn surfactants at the vapor/liquid interface

The effect of Ca(2+) ions on the hydration shell of sodium dodecyl carboxylate (SDC) and sodium dodecyl sulfonate (SDSn) monolayer at vapor/liquid interfaces was studied using molecular dynamics simulations. For each surfactant, two different surface concentrations were used to perform the simulations, and the aggregation morphologies and structural details have been reported. The results showed that the aggregation structures relate to both the surface coverage and the calcium ions. The divalent ions can screen the interaction between the polar head and Na(+) ions. Thus, Ca(2+) ions locate near the vapor/liquid interface to bind to the headgroup, making the aggregations much more compact via the salt bridge. The potential of mean force (PMF) between Ca(2+) and the headgroups shows that the interaction is decided by a stabilizing solvent-separated minimum in the PMF. To bind to the headgroup, Ca(2+) should overcome the energy barrier. Among contributions to the PMF, the major repulsive interaction was due to the rearrangement of the hydration shell after the calcium ions entered into the hydration shell of the headgroup. The PMFs between the headgroup and Ca(2+) in the SDSn systems showed higher energy barriers than those in the SDC systems. This result indicated that SDSn binds the divalent ions with more difficulty compared with SDC, so the ions have a strong effect on the hydration shell of SDC. That is why sulfonate surfactants have better efficiency in salt solutions with Ca(2+) ions for enhanced oil recovery.