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.

Characterizing Polymer Hydration Shell Compressibilities with the Small-System Method.

The small-system method (SSM) exploits the unique feature of finite-sized open systems, whose thermodynamic quantities scale with the inverse system size. This scaling enables the calculation of properties in the thermodynamic limit of macroscopic systems based on computer simulations of finite-sized systems. We herein extend the SSM to characterize the hydration shell compressibility of a generic hydrophobic polymer in water. By systematically increasing the strength of polymer-water repulsion, we find that the excess inverse thermodynamic correction factor (Δ1/Γs∞) and compressibility (Δχs) of the first hydration shell change sign from negative to positive. This occurs with a concurrent decrease in water hydrogen bonding and local tetrahedral order of the hydration shell water. The crossover lengthscale corresponds to an effective polymer bead diameter of 0.7 nm and is consistent with previous works on hydration of small and large hydrophobic solutes. The crossover lengthscale in polymer hydration shell compressibility, herein identified with the SSM approach, relates to hydrophobic interactions and macromolecular conformational equilibria in aqueous solution. The SSM approach may further be applied to study thermodynamic properties of polymer solvation shells in mixed solvents.