Molecular docking uses the three-dimensional structure of a receptor to screen databases of small molecules for potential ligands, often based on energetic complementarity. For many docking scoring functions, which calculate nonbonded interactions, metalloenzymes are challenging because of the partial covalent nature of metal-ligand interactions. To investigate how well molecular docking can identify potential ligands of metalloenzymes using a "standard" scoring function, we have docked the MDL Drug Data Report (MDDR), a functionally annotated database of 95 000 small molecules, against the X-ray crystal structures of five metalloenzymes. These enzymes included three zinc proteases, the nickel analogue of an iron enzyme, and a molybdenum metalloenzyme. The ability of the docking program to retrospectively enrich the annotated ligands as high-scoring hits for each enzyme and to calculate proper geometries was evaluated. In all five systems, the annotated ligands within the MDDR were enriched at least 20 times over random. To test the approach prospectively, a sixth target, the zinc -lactamase from Bacteroides fragilis, was screened against the fragment-like subset of the ZINC database. We purchased and tested 15 compounds from among the top 50 top-ranked ligands from docking, and found 5 inhibitors with apparent K i values less than 120 µM, the best of which was 2 µM. A more ambitious test still was predicting actual substrates for a seventh target, a Zn-dependent phosphotriesterase from Pseudomonas diminuta. Screening the Available Chemicals Directory (ACD) identified 25 thiophosphate esters as potential substrates within the top 100 ranked compounds. Eight of these, all previously uncharacterized for this enzyme, were acquired and tested, and all were confirmed experimentally as substrates. These results suggest that a simple, noncovalent scoring function may be used to identify inhibitors of at least some metalloenzymes.Metalloenzymes are important therapeutic targets, having attracted interest for diseases such as cancer (1, 2), arthritis (3), glaucoma (4), and infectious disease (5), among others. There has been considerable interest in the structural bases for their actions, and the structures of many metalloenzyme targets have been determined. These structures have provided templates for atomic resolution modeling of reaction mechanism and for inhibitor discovery.A problem faced in modeling metalloenzymes is the key role played by the metal. Molecular mechanics potential functions, which underlie most modeling approaches, have been parametrized to treat molecules as interacting either covalently or noncovalently. Metal-ligand interactions fit uncomfortably into either of these categories. Whereas vibrational spectroscopy indicates the presence of a metal ligand bond consistent with a force field treatment, the dative covalent "bond" is nearly an order of magnitude weaker than a covalent bond between first and second row elements. Further complicating the calculations is the unusually high charge ...