[Paper Review] Structure and Thermodynamics of Molecular Hydration via Grid Inhomogeneous Solvation Theory

Crystal N. Nguyen, Michael K. Gilson+1 more|arXiv (Cornell University)|Aug 24, 2011

Protein Structure and Dynamics2 references21 citations

TL;DR

This paper introduces Grid Inhomogeneous Solvation Theory (GIST), a computational method that maps solvation thermodynamics on a 3D grid to quantify local solvent contributions to free energy, entropy, and enthalpy. Applied to cucurbit[7]uril, GIST reveals a disfavored toroidal water region in the nonpolar cavity, whose entropic penalty contributes to high binding affinity, offering a rigorous framework for analyzing hydration in biomolecular recognition.

ABSTRACT

Changes in hydration are central to the phenomenon of biomolecular recognition, but it has been difficult to properly frame and answer questions about their precise thermodynamic role. We address this problem by introducing Grid Inhomogeneous Solvation Theory (GIST), which discretizes the equations of Inhomogeneous Solvation Theory on a 3D grid in a volume of interest. Here, the solvent volume is divided into small grid boxes and localized thermodynamic entropies, energies and free energies are defined for each grid box. Thermodynamic solvation quantities are defined in such a manner that summing the quantities over all the grid boxes yields the desired total quantity for the system. This approach smoothly accounts for the thermodynamics of not only highly occupied water sites but also partly occupied and water depleted regions of the solvent, without the need for ad hoc terms drawn from other theories. The GIST method has the further advantage of allowing a rigorous end-states analysis that, for example in the problem of molecular recognition, can account for not only the thermodynamics of displacing water from the surface but also for the thermodynamics of solvent reorganization around the bound complex. As a preliminary application, we present GIST calculations at the 1-body level for the host cucurbit[7]uril, a low molecular weight receptor molecule which represents a tractable model for biomolecular recognition. One of the most striking results is the observation of a toroidal region of water density, at the center of the host's nonpolar cavity, which is significantly disfavored entropically, and hence may contribute to the ability of this small receptor to bind guest molecules with unusually high affinities.

Motivation & Objective

• To develop a rigorous, spatially resolved method for calculating solvation thermodynamics in complex molecular environments.
• To address the challenge of quantifying hydration contributions—especially in partially occupied or depleted solvent regions—without ad hoc corrections.
• To enable end-states analysis of solvation, including solvent reorganization during binding, for accurate free energy prediction.
• To provide a framework that captures both enthalpic and entropic contributions to hydration at atomic resolution.
• To apply the method to a model system, cucurbit[7]uril, to reveal insights into high-affinity molecular recognition.

Proposed method

• The method discretizes Inhomogeneous Solvation Theory (IST) on a 3D grid, dividing the solvent volume into small cubic boxes.
• For each grid box, local thermodynamic quantities—free energy, entropy, and enthalpy—are computed using statistical mechanical averaging over solvent configurations.
• Total system-wide solvation properties are obtained by summing contributions from all grid boxes, preserving thermodynamic consistency.
• The approach accounts for both fully occupied and low-occupancy solvent regions, avoiding empirical corrections.
• The method enables end-states analysis by comparing hydration states before and after binding, capturing solvent reorganization effects.
• The framework is implemented at the 1-body level for initial validation on cucurbit[7]uril.

Experimental results

Research questions

• RQ1 How can solvation thermodynamics be accurately quantified in regions of low or variable water occupancy?
• RQ2 What is the contribution of local solvent entropy and enthalpy to the overall binding free energy of a molecular receptor?
• RQ3 Can a spatially resolved method detect unfavorable solvent structures that may drive high-affinity binding?
• RQ4 How does solvent reorganization during binding affect the thermodynamic profile of hydration?
• RQ5 To what extent can grid-based solvation theory capture the full thermodynamic landscape of hydration in nonpolar cavities?

Key findings

• A toroidal region of water density forms at the center of cucurbit[7]uril’s nonpolar cavity, indicating a structurally constrained and thermodynamically unfavorable hydration state.
• This central water region exhibits a strong entropic penalty, which contributes to the high binding affinity of the receptor by destabilizing the hydrated state.
• The GIST method successfully captures both highly occupied and low-occupancy solvent regions without requiring ad hoc corrections.
• The approach enables a complete thermodynamic analysis of hydration, including solvent reorganization effects during binding.
• The 1-body GIST calculation for cucurbit[7]uril demonstrates the method’s capability to reveal non-intuitive hydration features linked to binding thermodynamics.
• The results suggest that entropic destabilization of internal water can be a key driving force in molecular recognition.