A set of algorithms were designed, implemented and evaluated in order to, first, identifyclusters of conserved waters in binding pockets, i.e. hydration sites. Then, their contributionto the free energy of binding in a ligand-protein association was quantified by calculatingtheir enthalpy and entropy. The information obtained by using these algorithms couldcontribute to the development of new drugs by generating new ligands that target specifichigh-energy, unfavorable waters. Evaluation tests show that our algorithms can indeedprovide relevant data about how hydration sites influence ligand-protein binding.