Computational prediction of molecular hydration entropy with hybrid scaled particle theory and free-energy perturbation method.

Despite the importance of the knowledge of molecular hydration entropy (ΔShyd) in chemical and biological processes, the exact calculation of ΔShyd is very difficult, because of the complexity in solute-water interactions. Although free-energy perturbation (FEP) methods have been employed quite widely in the literature, the poor convergent behavior of the van der Waals interaction term in the potential function limited the accuracy and robustness. In this study, we propose a new method for estimating ΔShyd by means of combining the FEP approach and the scaled particle theory (or information theory) to separately calculate the electrostatic solute-water interaction term (ΔSelec) and the hydrophobic contribution approximated by the cavity formation entropy (ΔScav), respectively. Decomposition of ΔShyd into ΔScav and ΔSelec terms is found to be very effective with a substantial accuracy enhancement in ΔShyd estimation, when compared to the conventional full FEP calculations. ΔScav appears to dominate over ΔSelec in magnitude, even in the case of polar solutes, implying that the major contribution to the entropic cost for hydration comes from the formation of a solvent-excluded volume. Our hybrid scaled particle theory and FEP method is thus found to enhance the accuracy of ΔShyd prediction by effectively complementing the conventional full FEP method.