We map ligand binding sites on protein surfaces in molecular dynamics simulations using chlorobenzene as a probe molecule. The method was validated on four proteins. Two types of affinity maps that identified halogen and hydrophobic binding sites on proteins were obtained. Our method could prove useful for the discovery and development of halogenated inhibitors.