Multiple solvent crystal structures (MSCS) of porcine pancreatic elastase were used to map the binding surface the enzyme. Crystal structures of elastase in neat acetonitrile, 95% acetone, 55% dimethylformamide, 80% 5-hexene-1,2-diol, 80% isopropanol, 80% ethanol and 40% trifluoroethanol showed that the organic solvent molecules clustered in the active site, were found mostly unclustered in crystal contacts and in general did not bind elsewhere on the surface of elastase. Mixtures of 40% benzene or 40% cyclohexane in 50% isopropanol and 10% water showed no bound benzene or cyclohexane molecules, but did reveal bound isopropanol. The clusters of organic solvent probe molecules coincide with pockets occupied by known inhibitors. MSCS also reveal the areas of plasticity within the elastase binding site and allow for the visualization of a nearly complete first hydration shell. The pattern of organic solvent clusters determined by MSCS for elastase is consistent with patterns for hot spots in protein-ligand interactions determined from database analysis in general. The MSCS method allows probing of hot spots, plasticity and hydration simultaneously, providing a powerful complementary strategy to guide computational methods currently in development for binding site determination, ligand docking and design.