An interplay of excluded-volume and polymer-(co)solvent attractive interactions regulates polymer collapse in mixed solvents

Interactions of adsorbing cosolutes with hydrophobic hydration shells

The analysis of water density fluctuations in the hydration shell of nonpolar solutes provides insights into water-mediated interactions, especially hydrophobic interactions. These fluctuations are sensitive to small perturbations due to changes in thermodynamic conditions, such as temperature and pressure, but also to the presence of cosolutes, such as salts or small organic molecules. Herein, we investigate the effect of two classes of adsorbing cosolutes, using urea and methanol as representatives, on the fluctuations in energy and solvent density within the solvation shell of a model extended hydrophobic solute. We focus on the interactions of the cosolutes with the hydrophobic hydration shell, rather than with the solute itself, which though important remain largely unexplored. We calculate and analyze the interfacial partial molar energy of the cosolute, using a methodology based on the small system method. This approach provides correlated solvent density and energy fluctuations and allows us to decompose them into contributions due to interactions between the different components present in the solvation shell of the solute. The results show that adsorbed urea molecules interact more favorably with water than nonadsorbed urea molecules, which leads to the attenuation of interfacial density fluctuations and thus to the stabilization of the solvation shell. By contrast, the adsorbed methanol molecules interact preferably with other methanol molecules in the solvation shell, leading to a nano-phase segregated structure, which enhances interfacial fluctuations.

Salt Effects on Caffeine across Concentration Regimes

Salts affect the solvation thermodynamics of molecules of all sizes; the Hofmeister series is a prime example in which different ions lead to salting-in or salting-out of aqueous proteins. Early work of Tanford led to the discovery that the solvation of molecular surface motifs is proportional to the solvent accessible surface area (SASA), and later studies have shown that the proportionality constant varies with the salt concentration and type. Using multiscale computer simulations combined with vapor-pressure osmometry on caffeine-salt solutions, we reveal that this SASA description captures a rich set of molecular driving forces in tertiary solutions at changing solute and osmolyte concentrations. Central to the theoretical work is a new potential energy function that depends on the instantaneous surface area, salt type, and concentration. Used in, e.g., Monte Carlo simulations, this allows for a highly efficient exploration of many-body interactions and the resulting thermodynamics at elevated solute and salt concentrations.

Phase Transitions of Associative Biomacromolecules

Multivalent proteins and nucleic acids, collectively referred to as multivalent associative biomacromolecules, provide the driving forces for the formation and compositional regulation of biomolecular condensates. Here, we review the key concepts of phase transitions of aqueous solutions of associative biomacromolecules, specifically proteins that include folded domains and intrinsically disordered regions. The phase transitions of these systems come under the rubric of coupled associative and segregative transitions. The concepts underlying these processes are presented, and their relevance to biomolecular condensates is discussed.

Molecular Crowders Can Induce Collapse in Hydrophilic Polymers via Soft Attractive Interactions

A comprehensive understanding of protein folding and biomolecular self-assembly in the intracellular environment requires obtaining a microscopic view of the crowding effects. The classical view of crowding explains biomolecular collapse in such an environment in terms of the entropic solvent excluded volume effects subjected to hard-core repulsions exerted by the inert crowders, neglecting their soft chemical interactions. In this study, the effects of nonspecific, soft interactions of molecular crowders in regulating the conformational equilibrium of hydrophilic (charged) polymers are examined. Using advanced molecular dynamics simulations, collapse free energies of an uncharged, a negatively charged, and a charge-neutral 32-mer generic polymer are computed. The strength of the polymer-crowder dispersion energy is modulated to examine its effect on polymer collapse. The results show that the crowders preferentially adsorb and drive the collapse of all three polymers. The uncharged polymer collapse is opposed by the change in solute-solvent interaction energy but is overcompensated by the favorable change in the solute-solvent entropy as observed in hydrophobic collapse. However, the negatively charged polymer collapses with a favorable change in solute-solvent interaction energy due to reduction in the dehydration energy penalty as the crowders partition to the polymer interface and shield the charged beads. The collapse of a charge-neutral polymer is opposed by the solute-solvent interaction energy but is overcompensated by the solute-solvent entropy change. However, for the strongly interacting crowders, the overall energetic penalty decreases since the crowders interact with polymer beads via cohesive bridging attractions to induce polymer collapse. These bridging attractions are found to be sensitive to the binding sites of the polymer, since they are absent in the negatively charged or uncharged polymers. These interesting differences in thermodynamic driving forces highlight the crucial role of the chemical nature of the macromolecule as well as of the crowder in determining the conformational equilibria in a crowded milieu. The results emphasize that the chemical interactions of the crowders should be explicitly considered to account for the crowding effects. The findings have implications in understanding the crowding effects on the protein free energy landscapes.

The conformational phase diagram of neutral polymers in the presence of attractive crowders

Extensive coarse-grained molecular dynamics simulations are performed to investigate the conformational phase diagram of a neutral polymer in the presence of attractive crowders. We show that, for low crowder densities, the polymer predominantly shows three phases as a function of both intra-polymer and polymer-crowder interactions: (1) weak intra-polymer and weak polymer-crowder attractive interactions induce extended or coil polymer conformations (phase E), (2) strong intra-polymer and relatively weak polymer-crowder attractive interactions induce collapsed or globular conformations (phase CI), and (3) strong polymer-crowder attractive interactions, regardless of intra-polymer interactions, induce a second collapsed or globular conformation that encloses bridging crowders (phase CB). The detailed phase diagram is obtained by determining the phase boundaries delineating the different phases based on an analysis of the radius of gyration as well as bridging crowders. The dependence of the phase diagram on strength of crowder-crowder attractive interactions and crowder density is clarified. We also show that when the crowder density is increased, a third collapsed phase of the polymer emerges for weak intra-polymer attractive interactions. This crowder density-induced compaction is shown to be enhanced by stronger crowder-crowder attraction and is different from the depletion-induced collapse mechanism, which is primarily driven by repulsive interactions. We also provide a unified explanation of the observed re-entrant swollen/extended conformations of the earlier simulations of weak and strongly self-interacting polymers in terms of crowder-crowder attractive interactions.

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.

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.