Computing Spatially Resolved Rotational Hydration Entropies from Atomistic Simulations

For a first-principles understanding of macromolecular processes, a quantitative understanding of the underlying free energy landscape and in particular its entropy contribution is crucial. The stability of biomolecules, such as proteins, is governed by the hydrophobic effect, which arises from competing enthalpic and entropic contributions to the free energy of the solvent shell. While the statistical mechanics of liquids, as well as molecular dynamics simulations, have provided much insight, solvation shell entropies remain notoriously difficult to calculate, especially when spatial resolution is required. Here, we present a method that allows for the computation of spatially resolved rotational solvent entropies via a nonparametric k-nearest-neighbor density estimator. We validated our method using analytic test distributions and applied it to atomistic simulations of a water box. With an accuracy of better than 9.6%, the obtained spatial resolution should shed new light on the hydrophobic effect and the thermodynamics of solvation in general.