Length-Scale Effects in Hydrophobic Polymer Collapse Transitions

Concentration Fluctuation/Microheterogeneity Duality Illustrated with Aqueous 1,4-Dioxane Mixtures

The structural properties of aqueous 1-4 dioxane mixtures are studied by computer simulations of different water and dioxane force field models, from the perspective of illustrating the link between structural properties at the molecular level and measurable properties such as radiation scattering intensities and Kirkwood-Buff integrals (KBIs). A strategy to consistently correct the KBI obtained from simulations is proposed, which allows us to obtain the genuine KBI corresponding to a given pair of molecular species, in the entire concentration range, and without necessitating excessively large system sizes. The application of this method to the aqueous dioxane mixtures, with an all-atom CHARMM dioxane model and 2 water models, namely, SPC/E and TIP3P, allows one to understand the differences in the structure of the corresponding mixtures at the molecular level, particularly concerning the role of the water aggregates and its model dependence. This study allows us to characterize the dual role played by the concentration fluctuations and the domain segregation, particularly in what concerns the calculated X-ray spectra.

Distinctly different solvation behaviors of poly(N,N-diethylacrylamide) gels in water/acetone and water/DMSO mixtures

The solvation behaviors and intermolecular interactions of a poly(N,N-diethylacrylamide) (PDEA) gel network in water/DMSO and water/acetone mixtures have been investigated by variable-temperature high-resolution ¹H MAS NMR. Unlike decreasing volume phase transition temperature (VPTT) of the typical thermosensitive poly(N-isopropylacrylamide) (PNIPAM) gel induced by both acetone and DMSO in a water-rich region, distinct phase transition behaviors are revealed for the PDEA gel. That is, acetone is found to increase the VPTT of PDEA directly in the water-rich region while DMSO is also found to increase the VPTT of PDEA at a very low concentration but then decrease the VPTT as the concentration further increases. The above distinctly different VPTT shifts of PDEA are attributed to the different polymer-cosolvent interactions in water/acetone and water/DMSO systems. DMSO molecules with a strong water affinity are preferentially excluded by the PDEA gel network, and can promote the volume phase transition by favoring the collapse of the PDEA gel network, while acetone molecules preferentially adsorbed on the polymer interact with PDEA via non-specific van der Waals interaction, which favors the swollen state of the PDEA gel.

Connecting Non-Gaussian Water Density Fluctuations to the Lengthscale Dependent Crossover in Hydrophobic Hydration

We connect density fluctuations in liquid water to lengthscale dependent crossover in hydrophobic hydration. Specifically, we employ indirect umbrella sampling (INDUS) simulations to characterize density fluctuations in observation volumes of various sizes and shapes in water and as a function of temperature and salt concentration. Consistent with previous observations, density fluctuations are Gaussian in small molecular scale volumes, but they display non-Gaussian "low-density fat tails" in larger volumes. These non-Gaussian tails are indicative of the proximity of water to its liquid to vapor phase transition and have implications on biomolecular interactions and function. We show that the onset of non-Gaussian fluctuations in large volumes is accompanied by the formation of a cavity in the observation volume. We develop a model that uses the physics of cavity-water interface formation as a key ingredient and show that it captures the nature of non-Gaussian density fluctuations over a broad region in water and in salt solutions. We discuss the limitations of this model in the very low density region of the distribution. Our calculations provide new insights into the origins of non-Gaussian density fluctuations in water and their connections to lengthscale dependent crossover in hydrophobic hydration.

Equilibrium Swelling of Thermo-Responsive Gels in Mixtures of Solvents

Thermo-responsive (TR) gels of the LCST (lower critical solution temperature) type swell in water at temperatures below their volume phase transition temperature Tc and collapse above the critical temperature. When water is partially replaced with an organic liquid, these materials demonstrate three different types of equilibrium solvent uptake diagrams at temperatures below, above, in the close vicinity of Tc. A model is developed for equilibrium swelling of TR gels in binary mixtures of solvents. It takes into account three types of phase transitions in TR gels driven by (i) aggregation of hydrophobic side groups into clusters from which solvent molecules are expelled, (ii) replacement of water with cosolvent molecules in cage-like structures surrounding these groups, and (iii) replacement of water with cosolvent as the main element of hydration shells around backbone chains. The model involves a relatively small number of material constants that are found by matching observations on covalently cross-linked poly(N-isopropylacrylamide) macroscopic gels and microgels. Good agreement is demonstrated between the experimental data and results of numerical analysis. Classification is provided of the phase transition points on equilibrium swelling diagrams.

Solvation shell thermodynamics of extended hydrophobic solutes in mixed solvents

The ability of various cosolutes and cosolvents to enhance or quench solvent density fluctuations at solute–water interfaces has crucial implications on the conformational equilibrium of macromolecules such as polymers and proteins. Herein, we use an extended hydrophobic solute as a model system to study the effect of urea and methanol on the density fluctuations in the solute’s solvation shell and the resulting thermodynamics. On strengthening the solute–water/cosolute repulsive interaction, we observe distinct trends in the mutual affinities between various species in, and the thermodynamic properties of, the solvation shell. These trends strongly follow the respective trends in the preferential adsorption of urea and methanol: solute–water/cosolute repulsion strengthens, urea accumulation decreases, and methanol accumulation increases. Preferential accumulation of urea is found to quench the density fluctuations around the extended solute, leading to a decrease in the compressibility of the solvation shell. In contrast, methanol accumulation enhances the density fluctuations, leading to an increase in the compressibility. The mode of action of urea and methanol seems to be strongly coupled to their hydration behavior. The observations from this simple model is discussed in relation to urea driven swelling and methanol induced collapse of some well-known thermo-responsive polymers.

Cononsolvency of thermoresponsive polymers: where we are now and where we are going

Cononsolvency is an intriguing phenomenon where a polymer collapses in a mixture of good solvents. This cosolvent-induced modulation of the polymer solubility has been observed in solutions of several polymers and biomacromolecules, and finds application in areas such as hydrogel actuators, drug delivery, compound detection and catalysis. In the past decade, there has been a renewed interest in understanding the molecular mechanisms which drive cononsolvency with a predominant emphasis on its connection to the preferential adsorption of the cosolvent. Significant efforts have also been made to understand cononsolvency in complex systems such as micelles, block copolymers and thin films. In this review, we will discuss some of the recent developments from the experimental, simulation and theoretical fronts, and provide an outlook on the problems and challenges which are yet to be addressed.

Small-to-large length scale transition of TMAO interaction with hydrophobic solutes

We report the effect of trimethylamine N-oxide (TMAO) on the solvation of nonpolar solutes in water studied with molecular dynamics (MD) simulations and free-energy calculations. The simulation data indicate the occurrence of a length scale crossover in the TMAO interaction with repulsive Weeks-Chandler-Andersen (WCA) solutes: while TMAO is depleted from the hydration shell of a small WCA solute (methane) and increases the free-energy cost of solute-cavity formation, it preferentially binds to a large WCA solute (α-helical polyalanine), reducing the free-energy cost of solute-cavity formation via a surfactant-like mechanism. Significantly, we show that this surfactant-like behaviour of TMAO reinforces the solvent-mediated attraction between large WCA solutes by means of an entropic force linked to the interfacial accumulation of TMAO. Specifically, this entropic force arises from the natural tendency of adsorbed TMAO molecules to mix back into the bulk. It therefore favours solute-solute contact states that minimise the surface area exposed to the solvent and have a small overall number of TMAO molecules adsorbed. In contrast to the well-known depletion force, its effect is compensated by enthalpic solute-solvent interactions. Correspondingly, the hydrophobic association free energy of the large α-helical solutes passes through a minimum at low TMAO concentration when cohesive solute-solvent van der Waals interactions are considered. The observations reported herein are reminiscent to cosolvent effects on hydrophobic polymer coil-globule collapse free energies (Bharadwaj et al., Commun. Chem. 2020, 3, 165) and may be of general significance in systems whose properties are determined by hydrophobic self-assembly.

Direct Calculation of Entropic Components in Cohesive Interaction Free Energies

Cohesive interaction free energies entail an entropic component related to fluctuations of the energy associated with the attractive portion of the solute-solvent potential. The corresponding "fluctuation entropy" is fundamental in the solvation thermodynamics of macromolecular solutes and is linked to interfacial solvent density fluctuations and hydrophobic effects. Since the direct calculation of fluctuation entropy in molecular simulations is hampered by the poor sampling of high-energy tails in the solute-solvent energy distribution, indirect, and often approximate, routes for the calculation of fluctuation entropy are usually required, involving the modeling of geometrically frozen repulsive solute cavities in thermodynamic integration approaches. Herein, we propose a method to directly compute the fluctuation entropy by employing indirect umbrella sampling (INDUS). To validate the method, we consider model systems consisting of subnanometer oil droplets in water for which the fluctuation entropy can be computed exactly using indirect methods. The fluctuation entropy calculated with the newly proposed direct method agrees with the indirect reference calculations. We also observe that the solvation free energy and the contribution of the fluctuation entropy to it are of comparable magnitudes, particularly for larger oil droplets (∼1 nm). The proposed method can readily be employed for flexible macromolecular solutes and systems with extended hydrophobic surfaces or in the vicinity of a dewetting transition.