A new method for the representation of interactions between solutes and their hydration waters is presented. The algorithm is based on a generalization of the molecular interaction potential. The method, after a careful parametrization of the van der Waals interactions between the solute and the waters, provides a very accurate quasi-quantum mechanical picture of the hydration characteristics of the solute. Some applications of the method to the study of hydration patterns of relevant molecules are reported.