Here we apply the computational solvent mapping (CS-Map) algorithm toward the in silico identification of hot spots, that is, regions of protein binding sites that are major contributors to the binding energy and, hence, are prime targets in drug design. The CS-Map algorithm, developed for binding site characterization, moves small organic functional groups around the protein surface and determines their most energetically favorable binding positions. The utility of CS-Map algorithm toward the prediction of hot spot regions in druggable binding pockets is illustrated by three test systems: (1) renin aspartic protease, (2) a set of previously characterized druggable proteins, and (3) E. coli ketopantoate reductase. In each of the three studies, existing literature was used to verify our results. Based on our analyses, we conclude that the information provided by CS-Map can contribute substantially to the identification of hot spots, a necessary predecessor of fragment-based drug discovery efforts.