A new method is presented for docking molecular fragments to a rigid protein with evaluation of the binding energy. Polar fragments are docked with at least one hydrogen bond with the protein while apolar fragments are positioned in the hydrophobic pockets. The electrostatic contribution to the binding energy, which consists of screened intermolecular energy and protein and fragment desolvation terms, is evaluated efficiently by a numerical approach based on the continuum dielectric approximation. The latter is also used to predetermine the hydrophobic pockets of the protein by rolling a low dielectric sphere over the protein surface and calculating the electrostatic desolvation of the protein and van der Waals interaction energy. The method was implemented in the program SEED (solvation energy for exhaustive docking). The SEED continuum electrostatic approach has been successfully validated by a comparison with finite difference solutions of the Poisson equation for more than 2,500 complexes of small molecules with thrombin and the monomer of HIV-1 aspartic proteinase. The fragments docked by SEED in the active site of thrombin reproduce the structural features of the interaction patterns between known inhibitors and thrombin. Moreover, the combinatorial connection of these fragments yields a number of compounds that are very similar to potent inhibitors of thrombin.