Evaluation of Host-Guest Binding Thermodynamics of Model Cavities with Grid Cell Theory

A previously developed cell theory model of liquid water was used to evaluate the excess thermodynamic properties of confined clusters of water molecules. The results are in good agreement with reference thermodynamic integration calculations, suggesting that the model is adequate to probe the thermodynamic properties of water at interfaces or in cavities. Next, the grid cell theory (GCT) method was applied to elucidate the thermodynamic signature of nonpolar association for a range of idealized host-guest systems. Polarity and geometry of the host cavities were systematically varied, and enthalpic and entropic solvent components were spatially resolved for detailed graphical analyses. Perturbations in the thermodynamic properties of water molecules upon guest binding are restricted to the immediate vicinity of the guest in solvent-exposed cavities, whereas longer-ranged perturbations are observed in buried cavities. Depending on the polarity and geometry of the host, water displacement by a nonpolar guest makes a small or large enthalpic or entropic contribution to the free energy of binding. Thus, no assumptions about the thermodynamic signature of the hydrophobic effect can be made in general. Overall the results warrant further applications of GCT to more complex systems such as protein-ligand complexes.