Integral equation study of hydrophobic interaction: A comparison between the simple point charge model for water and a Lennard-Jones model for solvent

Pair Potential between Poloxamer 407 Micelles in 14 wt % Aqueous Solution Clarifying the Sol-Gel-Sol Transition by Temperature Rise

Poloxamer 407 (P407) is used as a safety-guaranteed, invaluable pharmaceutical nanocarrier. The aqueous solution of P407 exhibits sol-to-gel and gel-to-sol transitions, specifically during a temperature rise. Here, we develop a method to determine the pair potential between colloidal particles based primarily on experimental small-angle scattering data. Using this approach, the pair potential between the P407 micelles in the sol-gel-sol transition state is revealed without prelimiting any type of interaction forces. The results indicate that the increase in the attractive interaction contributes to enhancing the volume fraction, which is the decisive parameter for the gelation in terms of the Alder transition theory, i.e., the fluid-to-crystal phase transition of the hard-sphere model in the micelle system. This study demonstrates that one of the key mechanisms of the gel-to-sol transition upon heating is the enhancement of the structural fluctuation due to widening of the potential well in the intermicellar pair potential.

How do water-mediated interactions and osmotic second virial coefficients vary with particle size?

We examine quantitatively the solute-size dependences of the effective interactions between nonpolar solutes in water and in a simple liquid. The potential w(r) of mean force and the osmotic second virial coefficients B are calculated with high accuracy from molecular dynamics simulations. As the solute diameter increases from methane's to C60's with the solute-solute and solute-solvent attractive interaction parameters fixed to those for the methane-methane and methane-water interactions, the first minimum of w(r) lowers from -1.1 to -4.7 in units of the thermal energy kT. Correspondingly, the magnitude of B (<0) increases proportional to σα with some power close to 6 or 7, which reinforces the solute-size dependence of B found earlier for a smaller range of σ [H. Naito, R. Okamoto, T. Sumi and K. Koga, J. Chem. Phys., 2022, 156, 221104]. We also demonstrate that the strength of the attractive interactions between solute and solvent molecules can qualitatively change the characteristics of the effective pair interaction between solute particles, both in water and in a simple liquid. If the solute-solvent attractive force is set to be weaker (stronger) than a threshold, the effective interaction becomes increasingly attractive (repulsive) with increasing solute size.

Thermodynamic, structural, and dynamical properties of nano-confined water using SPC/E and TIP4P models by molecular dynamics simulations

Different morphologies of water molecules are confined between two parallel graphene surfaces.

Model-potential-free analysis of small angle scattering of proteins in solution: insights into solvent effects on protein–protein interaction

The model-potential-free method advances the analysis of small-angle scattering for protein solutions and provides new insights into protein–protein interactions.

Solvation free energies of alanine peptides: the effect of flexibility

The electrostatic (ΔGel), van der Waals cavity-formation (ΔGvdw), and total (ΔG) solvation free energies for 10 alanine peptides ranging in length (n) from 1 to 10 monomers were calculated. The free energies were computed both with fixed, extended conformations of the peptides and again for some of the peptides without constraints. The solvation free energies, ΔGel, and components ΔGvdw, and ΔG, were found to be linear in n, with the slopes of the best-fit lines being γel, γvdw, and γ, respectively. Both γel and γ were negative for fixed and flexible peptides, and γvdw was negative for fixed peptides. That γvdw was negative was surprising, as experimental data on alkanes, theoretical models, and MD computations on small molecules and model systems generally suggest that γvdw should be positive. A negative γvdw seemingly contradicts the notion that ΔGvdw drives the initial collapse of the protein when it folds by favoring conformations with small surface areas. When we computed ΔGvdw for the flexible peptides, thereby allowing the peptides to assume natural ensembles of more compact conformations, γvdw was positive. Because most proteins do not assume extended conformations, a ΔGvdw that increases with increasing surface area may be typical for globular proteins. An alternative hypothesis is that the collapse is driven by intramolecular interactions. We find few intramolecular H-bonds but show that the intramolecular van der Waals interaction energy is more favorable for the flexible than for the extended peptides, seemingly favoring this hypothesis. The large fluctuations in the vdw energy may make attributing the collapse of the peptide to this intramolecular energy difficult.

Soft collective fluctuations governing hydrophobic association

The interaction between two associating hydrophobic particles has traditionally been explained in terms of the release of entropically frustrated hydration shell water molecules. However, this picture cannot account for the kinetics of hydrophobic association and is therefore not capable of providing a microscopic description of the hydrophobic interaction (HI). Here, Monte Carlo simulations of a pair of molecular-scale apolar solutes in aqueous solution reveal the critical role of collective fluctuations in the hydrogen bond (HB) network for the microscopic picture of the HI. The main contribution to the HI is the relaxation of solute-water translational correlations. The existence of a heat capacity maximum at the desolvation barrier is shown to arise from softening of non-HB water fluctuations and the relaxation of many-body correlations in the labile HB network. The microscopic event governing the kinetics of hydrophobic association has turned out to be a relatively large critical collective fluctuation in hydration water displacing a substantial fraction of HB clusters from the inner to the outer region of the first hydration shell.