A new method is presented to quantitatively estimate and graphically display the propensity of nonpolar groups to bind at the surface of proteins. It is based on the calculation of the binding energy, i.e., van der Waals interaction plus protein electrostatic desolvation, of a nonpolar probe sphere rolled over the protein surface, and on the color coding of this quantity on a smooth molecular surface (hydrophobicity map). The method is validated on ten protein-ligand complexes and is shown to distinguish precisely where polar and nonpolar groups preferentially bind. Comparisons with existing approaches, like the display of the electrostatic potential or the curvature, illustrate the advantages and the better predictive power of the present method. Hydrophobicity maps will play an important role in the characterization of binding sites for the large number of proteins emerging from the genome projects and structure modeling approaches.