In this Communication, we have further developed and applied our grand canonical nonequilibrium Monte Carlo method to sample the binding of small molecules to protein systems with the goal of providing a tool which can help in the discovery of new fragment leads. The method has been successfully applied to detect binding sites, elucidate binding geometries, and calculate binding affinities.