Enzymes from the Haloacid Dehalogenase (HAD) Superfamily Catalyse the Elusive Dephosphorylation Step of Riboflavin Biosynthesis

Discovery of a novel symbiotic lineage associated with a hematophagous leech from the genus Haementeria

Similarly to other strict blood feeders, leeches from the Haementeria genus (Hirudinida: Glossiphoniidae) have established a symbiotic association with bacteria harbored intracellularly in esophageal bacteriomes. Previous genome sequence analyses of these endosymbionts revealed co-divergence with their hosts, a strong genome reduction, and a simplified metabolism largely dedicated to the production of B vitamins, which are nutrients lacking from a blood diet. 'Candidatus Providencia siddallii' has been identified as the obligate nutritional endosymbiont of a monophyletic clade of Mexican and South American Haementeria spp. However, the Haementeria genus includes a sister clade of congeners from Central and South America, where the presence or absence of the aforementioned symbiont taxon remains unknown. In this work, we report on a novel bacterial endosymbiont found in a representative from this Haementeria clade. We found that this symbiont lineage has evolved from within the Pluralibacter genus, known mainly from clinical but also environmental strains. Similarly to Ca. Providencia siddallii, the Haementeria-associated Pluralibacter symbiont displays clear signs of genome reduction, accompanied by an A+T-biased sequence composition. Genomic analysis of its metabolic potential revealed a retention of pathways related to B vitamin biosynthesis, supporting its role as a nutritional endosymbiont. Finally, comparative genomics of both Haementeria symbiont lineages suggests that an ancient Providencia symbiont was likely replaced by the novel Pluralibacter one, thus constituting the first reported case of nutritional symbiont replacement in a leech without morphological changes in the bacteriome.

IMPORTANCE

Obligate symbiotic associations with a nutritional base have likely evolved more than once in strict blood-feeding leeches. Unlike those symbioses found in hematophagous arthropods, the nature, identity, and evolutionary history of these remains poorly studied. In this work, we further explored obligate nutritional associations between Haementeria leeches and their microbial symbionts, which led to the unexpected discovery of a novel symbiosis with a member of the Pluralibacter genus. When compared to Providencia siddallii, an obligate nutritional symbiont of other Haementeria leeches, this novel bacterial symbiont shows convergent retention of the metabolic pathways involved in B vitamin biosynthesis. Moreover, the genomic characteristics of this Pluralibacter symbiont suggest a more recent association than that of Pr. siddallii and Haementeria. We conclude that the once-thought stable associations between blood-feeding Glossiphoniidae and their symbionts (i.e., one bacteriome structure, one symbiont lineage) can break down, mirroring symbiont turnover observed in various arthropod lineages.

Systematic Analysis of the Effect of Genomic Knock-Out of Non-Essential Promiscuous HAD-Like Phosphatases YcsE, YitU and YwtE on Flavin and Adenylate Content in Bacillus Subtilis

Studying the metabolic role of non-essential promiscuous enzymes is a challenging task, as genetic manipulations usually do not reveal at which point(s) of the metabolic network the enzymatic activity of such protein is beneficial for the organism. Each of the HAD-like phosphatases YcsE, YitU and YwtE of Bacillus subtilis catalyzes the dephosphorylation of 5-amino-6-ribitylamino-uracil 5'-phosphate, which is essential in the biosynthesis of riboflavin. Using CRISPR technology, we have found that the deletion of these genes, individually or in all possible combinations failed to cause riboflavin auxotrophy and did not result in significant growth changes. Analysis of flavin and adenylate content in B. subtilis knockout mutants showed that (i) there must be one or several still unidentified phosphatases that can replace the deleted proteins; (ii) such replacements, however, cannot fully restore the intracellular content of any of three flavins studied (riboflavin, FMN, FAD); (iii) whereas bacterial fitness was not significantly compromised by mutations, the intracellular balance of flavins and adenylates did show some significant changes.

Arsenophonus symbiosis with louse flies: multiple origins, coevolutionary dynamics, and metabolic significance

IMPORTANCE

Insects that live exclusively on vertebrate blood utilize symbiotic bacteria as a source of essential compounds, e.g., B vitamins. In louse flies, the most frequent symbiont originated in genus Arsenophonus, known from a wide range of insects. Here, we analyze genomic traits, phylogenetic origins, and metabolic capacities of 11 Arsenophonus strains associated with louse flies. We show that in louse flies, Arsenophonus established symbiosis in at least four independent events, reaching different stages of symbiogenesis. This allowed for comparative genomic analysis, including convergence of metabolic capacities. The significance of the results is twofold. First, based on a comparison of independently originated Arsenophonus symbioses, it determines the importance of individual B vitamins for the insect host. This expands our theoretical insight into insect-bacteria symbiosis. The second outcome is of methodological significance. We show that the comparative approach reveals artifacts that would be difficult to identify based on a single-genome analysis.

Evolution of an Alternative Genetic Code in the Providencia Symbiont of the Hematophagous Leech Haementeria acuecueyetzin

Strict blood-feeding animals are confronted with a strong B-vitamin deficiency. Blood-feeding leeches from the Glossiphoniidae family, similarly to hematophagous insects, have evolved specialized organs called bacteriomes to harbor symbiotic bacteria. Leeches of the Haementeria genus have two pairs of globular bacteriomes attached to the esophagus which house intracellular "Candidatus Providencia siddallii" bacteria. Previous work analyzing a draft genome of the Providencia symbiont of the Mexican leech Haementeria officinalis showed that, in this species, the bacteria hold a reduced genome capable of synthesizing B vitamins. In this work, we aimed to expand our knowledge on the diversity and evolution of Providencia symbionts of Haementeria. For this purpose, we sequenced the symbiont genomes of three selected leech species. We found that all genomes are highly syntenic and have kept a stable genetic repertoire, mirroring ancient insect endosymbionts. Additionally, we found B-vitamin pathways to be conserved among these symbionts, pointing to a conserved symbiotic role. Lastly and most notably, we found that the symbiont of H. acuecueyetzin has evolved an alternative genetic code, affecting a portion of its proteome and showing evidence of a lineage-specific and likely intermediate stage of genetic code reassignment.

Enhancing the biosynthesis of riboflavin in the recombinant Escherichia coli BL21 strain by metabolic engineering

In this study, to construct the riboflavin-producing strain R1, five key genes, ribA, ribB, ribC, ribD, and ribE, were cloned and ligated to generate the plasmid pET-AE, which was overexpressed in Escherichia coli BL21. The R1 strain accumulated 182.65 ± 9.04 mg/l riboflavin. Subsequently, the R2 strain was constructed by the overexpression of zwf harboring the constructed plasmid pAC-Z in the R1 strain. Thus, the level of riboflavin in the R2 strain increased to 319.01 ± 20.65 mg/l (74.66% increase). To further enhance ribB transcript levels and riboflavin production, the FMN riboswitch was deleted from E. coli BL21 with CRISPR/Cas9 to generate the R3 strain. The R4 strain was constructed by cotransforming pET-AE and pAC-Z into the R3 strain. Compared to those of E. coli BL21, the ribB transcript levels of R2 and R4 improved 2.78 and 3.05-fold, respectively. The R4 strain accumulated 437.58 ± 14.36 mg/l riboflavin, increasing by 37.17% compared to the R2 strain. These results suggest that the deletion of the FMN riboswitch can improve the transcript level of ribB and facilitate riboflavin production. A riboflavin titer of 611.22 ± 11.25 mg/l was achieved under the optimal fermentation conditions. Ultimately, 1574.60 ± 109.32 mg/l riboflavin was produced through fed-batch fermentation with 40 g/l glucose. This study contributes to the industrial production of riboflavin by the recombinant E. coli BL21.

Phylogenomic analysis and metabolic role reconstruction of mutualistic Rhizobiales hindgut symbionts of Acromyrmex leaf-cutting ants

Rhizobiales are well-known plant-root nitrogen-fixing symbionts, but the functions of insect-associated Rhizobiales are poorly understood. We obtained genomes of three strains associated with Acromyrmex leaf-cutting ants and show that, in spite of being extracellular gut symbionts, they lost all pathways for essential amino acid biosynthesis, making them fully dependent on their hosts. Comparison with 54 Rhizobiales genomes showed that all insect-associated Rhizobiales lost the ability to fix nitrogen and that the Acromyrmex symbionts had exceptionally also lost the urease genes. However, the Acromyrmex strains share biosynthesis pathways for riboflavin vitamin, queuosine and a wide range of antioxidant enzymes likely to be beneficial for the ant fungus-farming symbiosis. We infer that the Rhizobiales symbionts catabolize excess of fungus-garden-derived arginine to urea, supplementing complementary Mollicutes symbionts that turn arginine into ammonia and infer that these combined symbiont activities stabilize the fungus-farming mutualism. Similar to the Mollicutes symbionts, the Rhizobiales species have fully functional CRISPR/Cas and R-M phage defenses, suggesting that these symbionts are important enough for the ant hosts to have precluded the evolution of metabolically cheaper defenseless strains.

High-yield production of coenzyme F420 in Escherichia coli by fluorescence-based screening of multi-dimensional gene expression space

Coenzyme F₄₂₀ is involved in bioprocesses such as biosynthesis of antibiotics by streptomycetes, prodrug activation in Mycobacterium tuberculosis, and methanogenesis in archaea. F₄₂₀-dependent enzymes also attract interest as biocatalysts in organic chemistry. However, as only low F₄₂₀ levels are produced in microorganisms, F₄₂₀ availability is a serious bottleneck for research and application. Recent advances in our understanding of the F₄₂₀ biosynthesis enabled heterologous overproduction of F₄₂₀ in Escherichia coli, but the yields remained moderate. To address this issue, we rationally designed a synthetic operon for F₄₂₀ biosynthesis in E. coli. However, it still led to the production of low amounts of F₄₂₀ and undesired side-products. In order to strongly improve yield and purity, a screening approach was chosen to interrogate the gene expression-space of a combinatorial library based on diversified promotors and ribosome binding sites. The whole pathway was encoded by a two-operon construct. The first module ("core") addressed parts of the riboflavin biosynthesis pathway and F_O synthase for the conversion of GTP to the stable F₄₂₀ intermediate F_O. The enzymes of the second module ("decoration") were chosen to turn F_O into F₄₂₀. The final construct included variations of T7 promoter strengths and ribosome binding site activity to vary the expression ratio for the eight genes involved in the pathway. Fluorescence-activated cell sorting was used to isolate clones of this library displaying strong F₄₂₀-derived fluorescence. This approach yielded the highest titer of coenzyme F₄₂₀ produced in the widely used organism E. coli so far. Production in standard LB medium offers a highly effective and simple production process that will facilitate basic research into unexplored F₄₂₀-dependent bioprocesses as well as applications of F₄₂₀-dependent enzymes in biocatalysis.

The Di-Symbiotic Systems in the Aphids Sipha maydis and Periphyllus lyropictus Provide a Contrasting Picture of Recent Co-Obligate Nutritional Endosymbiosis in Aphids

Dependence on multiple nutritional bacterial symbionts forming a metabolic unit has repeatedly evolved in many insect species that feed on nutritionally unbalanced diets such as plant sap. This is the case for aphids of the subfamilies Lachninae and Chaitophorinae, which have evolved di-symbiotic systems in which the ancient obligate nutritional symbiont Buchnera aphidicola is metabolically complemented by an additional nutritional symbiont acquired more recently. Deciphering how different symbionts integrate both metabolically and anatomically in such systems is crucial to understanding how complex nutritional symbiotic systems function and evolve. In this study, we sequenced and analyzed the genomes of the symbionts B. aphidicola and Serratia symbiotica associated with the Chaitophorinae aphids Sipha maydis and Periphyllus lyropictus. Our results show that, in these two species, B. aphidicola and S. symbiotica complement each other metabolically (and their hosts) for the biosynthesis of essential amino acids and vitamins, but with distinct metabolic reactions supported by each symbiont depending on the host species. Furthermore, the S. symbiotica symbiont associated with S. maydis appears to be strictly compartmentalized into the specialized host cells housing symbionts in aphids, the bacteriocytes, whereas the S. symbiotica symbiont associated with P. lyropictus exhibits a highly invasive phenotype, presumably because it is capable of expressing a larger set of virulence factors, including a complete flagellum for bacterial motility. Such contrasting levels of metabolic and anatomical integration for two S. symbiotica symbionts that were recently acquired as nutritional co-obligate partners reflect distinct coevolutionary processes specific to each association.

Microbial 5'-nucleotidases: their characteristics, roles in cellular metabolism, and possible practical applications

5′-Nucleotidases (EC 3.1.3.5) are enzymes that catalyze the hydrolytic dephosphorylation of 5′-ribonucleotides and 5′-deoxyribonucleotides to their respective nucleosides and phosphate. Most 5′-nucleotidases have broad substrate specificity and are multifunctional enzymes capable of cleaving phosphorus from not only mononucleotide phosphate molecules but also a variety of other phosphorylated metabolites. 5′-Nucleotidases are widely distributed throughout all kingdoms of life and found in different cellular locations. The well-studied vertebrate 5′-nucleotidases play an important role in cellular metabolism. These enzymes are involved in purine and pyrimidine salvage pathways, nucleic acid repair, cell-to-cell communication, signal transduction, control of the ribo- and deoxyribonucleotide pools, etc. Although the first evidence of microbial 5′-nucleotidases was obtained almost 60 years ago, active studies of genetic control and the functions of microbial 5′-nucleotidases started relatively recently. The present review summarizes the current knowledge about microbial 5′-nucleotidases with a focus on their diversity, cellular localizations, molecular structures, mechanisms of catalysis, physiological roles, and activity regulation and approaches to identify new 5′-nucleotidases. The possible applications of these enzymes in biotechnology are also discussed.

Key points

• Microbial 5′-nucleotidases differ in molecular structure, hydrolytic mechanism, and cellular localization.

• 5′-Nucleotidases play important and multifaceted roles in microbial cells.

• Microbial 5′-nucleotidases have wide range of practical applications.

The novel phosphatase RosC catalyzes the last unknown step of roseoflavin biosynthesis in Streptomyces davaonensis

The bacterium Streptomyces davaonensis produces the antibiotic roseoflavin, which is a riboflavin (vitamin B₂ ) analog. The key enzyme of roseoflavin biosynthesis is the 8-demethyl-8-amino-riboflavin-5'-phosphate (AFP) synthase RosB which synthesizes AFP from riboflavin-5'-phosphate. AFP is not a substrate for the last enzyme of roseoflavin biosynthesis the N, N-dimethyltransferase RosA, which generates roseoflavin from 8-demethyl-8-amino-riboflavin (AF). Consequently, the roseoflavin biosynthetic pathway depends on a phosphatase, which dephosphorylates AFP to AF. Here, we report on the identification and characterization of such an AFP phosphatase which we named RosC. The gene rosC is located immediately downstream of rosA and both genes are part of a cluster comprising 10 genes. Deletion of rosC from the chromosome of S. davaonensis led to reduced roseoflavin levels in the corresponding recombinant strain. In contrast to wild-type S. davaonensis, cell-free extracts of the rosC deletion strain did not catalyze dephosphorylation of AFP. RosC was purified from an overproducing Escherichia coli strain. RosC is the fastest enzyme of roseoflavin biosynthesis (k_cat 31.3 ± 1.4 min^-1 ). The apparent K_M for the substrate AFP was 34.5 µM. Roseoflavin biosynthesis is now completely understood--it takes three enzymes (RosB, RosC, and RosA) to convert the flavin cofactor riboflavin-5'-phosphate into a potent antibiotic.

Bacterial symbionts support larval sap feeding and adult folivory in (semi-)aquatic reed beetles

Symbiotic microbes can enable their host to access untapped nutritional resources but may also constrain niche space by promoting specialization. Here, we reconstruct functional changes in the evolutionary history of the symbiosis between a group of (semi-)aquatic herbivorous insects and mutualistic bacteria. Sequencing the symbiont genomes across 26 species of reed beetles (Chrysomelidae, Donaciinae) spanning four genera indicates that the genome-eroded mutualists provide life stage-specific benefits to larvae and adults, respectively. In the plant sap-feeding larvae, the symbionts are inferred to synthesize most of the essential amino acids as well as the B vitamin riboflavin. The adult reed beetles’ folivory is likely supported by symbiont-encoded pectinases that complement the host-encoded set of cellulases, as revealed by transcriptome sequencing. However, mapping the occurrence of the symbionts’ pectinase genes and the hosts’ food plant preferences onto the beetles’ phylogeny reveals multiple independent losses of pectinase genes in lineages that switched to feeding on pectin-poor plants, presumably constraining their hosts’ subsequent adaptive potential.

Symbiotic microbes in insects can enable their hosts to access untapped nutritional resources. Here, the authors show that symbiotic bacteria in reed beetles can provide essential amino acids to sap-feeding larvae and help leaf-feeding adults to degrade pectin, respectively.

Production of riboflavin and related cofactors by biotechnological processes

Riboflavin (RF) and its active forms, the cofactors flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD), have been extensively used in the food, feed and pharmaceutical industries. Modern commercial production of riboflavin is based on microbial fermentation, but the established genetically engineered production strains are facing new challenges due to safety concerns in the food and feed additives industry. High yields of flavin mononucleotide and flavin adenine dinucleotide have been obtained using whole-cell biocatalysis processes. However, the necessity of adding expensive precursors results in high production costs. Consequently, developing microbial cell factories that are capable of efficiently producing flavin nucleotides at low cost is an increasingly attractive approach. The biotechnological processes for the production of RF and its cognate cofactors are reviewed in this article.

Identification and characterization of the missing phosphatase on the riboflavin biosynthesis pathway in Arabidopsis thaliana

Despite the importance of riboflavin as the direct precursor of the cofactors flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN), the physiologically relevant catalyst dephosphorylating the riboflavin biosynthesis pathway intermediate 5-amino-6-ribitylamino-2,4(1H,3H) pyrimidinedione 5'-phosphate (ARPP) has not been characterized from any organism. By using as the query sequence a previously identified plastidial FMN hydrolase AtcpFHy1 (At1g79790), belonging to the haloacid dehalogenase (HAD) superfamily, seven candidates for the missing ARPP phosphatase were found, cloned, recombinantly expressed, and purified. Activity screening showed that the enzymes encoded by AtcpFHy1, At4g11570, and At4g25840 catalyze dephosphorylation of ARPP. AtcpFHy1 was renamed AtcpFHy/PyrP1, At4g11570 and At4g25840 were named AtPyrP2 and AtGpp1/PyrP3, respectively. Subcellular localization in planta indicated that AtPyrP2 was localized in plastids and AtGpp1/PyrP3 in mitochondria. Biochemical characterization of AtcpFHy/PyrP1 and AtPyrP2 showed that they have similar K_m values for the substrate ARPP, with AtcpFHy/PyrP1 having higher catalytic efficiency. Screening of 21 phosphorylated substrates showed that AtPyrP2 is specific for ARPP. Molecular weights of AtcpFHy/PyrP1 and AtPyrP2 were estimated at 46 and 72 kDa, suggesting dimers. pH and temperature optima for AtcpFHy/PyrP1 and AtPyrP2 were ~7.0-8.5 and 40-50°C. T-DNA knockout of AtcpFHy/PyrP1 did not affect the flavin profile of the transgenic plants, whereas silencing of AtPyrP2 decreased accumulation of riboflavin, FMN, and FAD. Our results strongly support AtPyrP2 as the missing phosphatase on the riboflavin biosynthesis pathway in Arabidopsis thaliana. The identification of this enzyme closes a long-standing gap in understanding of the riboflavin biosynthesis in plants.

Genomic approach to studying nutritional requirements of Clostridium tyrobutyricum and other Clostridia causing late blowing defects

Clostridium tyrobutyricum is the main microorganism responsible for the late blowing defect in hard and semi-hard cheeses, causing considerable economic losses to the cheese industry. Deeper knowledge of the metabolic requirements of this microorganism can lead to the development of more effective control approaches. In this work, the amino acids and B vitamins essential for sustaining the growth of C. tyrobutyricum were investigated using a genomic approach. As the first step, the genomes of four C. tyrobutyricum strains were analyzed for the presence of genes putatively involved in the biosynthesis of amino acids and B vitamins. Metabolic pathways could be reconstructed for all amino acids and B vitamins with the exception of biotin (vitamin B7) and folate (vitamin B9). The biotin pathway was missing the enzyme amino-7-oxononanoate synthase that catalyzes the condensation of pimeloyl-ACP and l-alanine to 8-amino-7-oxononanoate. In the folate pathway, the missing genes were those coding for para-aminobenzoate synthase and aminodeoxychorismate lyase enzymes. These enzymes are responsible for the conversion of chorismate into para-aminobenzoate (PABA). Two C. tyrobutyircum strains whose genome was analyzed in silico as well as other 10 strains isolated from cheese were tested in liquid media to confirm these observations. 11 strains showed growth in a defined liquid medium containing biotin and PABA after 6-8 days of incubation. No strain showed growth when only one or none of these compounds were added, confirming the observations obtained in silico. Furthermore, the genome analysis was extended to genomes of single strains of other Clostridium species potentially causing late blowing, namely Clostridium beijerinckii, Clostridium sporogenes and Clostridium butyricum. Only the biotin biosynthesis pathway was incomplete for C. butyricum and C. beijerincki. In contrast, C. sporogenes showed missing enzymes in biosynthesis pathways of several amino acids as well as biotin, folate, and cobalamin (vitamin B12). These observations agree with the results of growth experiments of these species in liquid media reported in the literature. The results of this study suggest that biotin and folate are potential targets for reducing late blowing in cheese and highlight the usefulness of genomic analysis for identifying essential nutrients in bacteria.

Reinventing the Wheel and Making It Round Again: Evolutionary Convergence in Buchnera-Serratia Symbiotic Consortia between the Distantly Related Lachninae Aphids Tuberolachnus salignus and Cinara cedri

Virtually all aphids (Aphididae) harbor Buchnera aphidicola as an obligate endosymbiont to compensate nutritional deficiencies arising from their phloem diet. Many species within the Lachninae subfamily seem to be consistently associated also with Serratia symbiotica We have previously shown that both Cinara (Cinara) cedri and Cinara (Cupressobium) tujafilina (Lachninae: Eulachnini tribe) have indeed established co-obligate associations with both Buchnera and S. symbiotica However, while Buchnera genomes of both Cinara species are similar, genome degradation differs greatly between the two S. symbiotica strains. To gain insight into the essentiality and degree of integration of S. symbiotica within the Lachninae, we sequenced the genome of both Buchnera and S. symbiotica endosymbionts from the distantly related aphid Tuberolachnus salignus (Lachninae: Tuberolachnini tribe). We found a striking level of similarity between the endosymbiotic system of this aphid and that of C. cedri In both aphid hosts, S. symbiotica possesses a highly reduced genome and is found exclusively intracellularly inside bacteriocytes. Interestingly, T. salignus' endosymbionts present the same tryptophan biosynthetic metabolic complementation as C. cedri's, which is not present in C. tujafilina's. Moreover, we corroborate the riboflavin-biosynthetic-role take-over/rescue by S. symbiotica in T. salignus, and therefore, provide further evidence for the previously proposed establishment of a secondary co-obligate endosymbiont in the common ancestor of the Lachninae aphids. Finally, we propose that the putative convergent split of the tryptophan biosynthetic role between Buchnera and S. symbiotica could be behind the establishment of S. symbiotica as an obligate intracellular symbiont and the triggering of further genome degradation.

Sequence-based Network Completion Reveals the Integrality of Missing Reactions in Metabolic Networks

Genome-scale metabolic models are central in connecting genotypes to metabolic phenotypes. However, even for well studied organisms, such as Escherichia coli, draft networks do not contain a complete biochemical network. Missing reactions are referred to as gaps. These gaps need to be filled to enable functional analysis, and gap-filling choices influence model predictions. To investigate whether functional networks existed where all gap-filling reactions were supported by sequence similarity to annotated enzymes, four draft networks were supplemented with all reactions from the Model SEED database for which minimal sequence similarity was found in their genomes. Quadratic programming revealed that the number of reactions that could partake in a gap-filling solution was vast: 3,270 in the case of E. coli, where 72% of the metabolites in the draft network could connect a gap-filling solution. Nonetheless, no network could be completed without the inclusion of orphaned enzymes, suggesting that parts of the biochemistry integral to biomass precursor formation are uncharacterized. However, many gap-filling reactions were well determined, and the resulting networks showed improved prediction of gene essentiality compared with networks generated through canonical gap filling. In addition, gene essentiality predictions that were sensitive to poorly determined gap-filling reactions were of poor quality, suggesting that damage to the network structure resulting from the inclusion of erroneous gap-filling reactions may be predictable.

Covalent docking predicts substrates for haloalkanoate dehalogenase superfamily phosphatases

Enzyme function prediction remains an important open problem. Though structure-based modeling, such as metabolite docking, can identify substrates of some enzymes, it is ill-suited to reactions that progress through a covalent intermediate. Here we investigated the ability of covalent docking to identify substrates that pass through such a covalent intermediate, focusing particularly on the haloalkanoate dehalogenase superfamily. In retrospective assessments, covalent docking recapitulated substrate binding modes of known cocrystal structures and identified experimental substrates from a set of putative phosphorylated metabolites. In comparison, noncovalent docking of high-energy intermediates yielded nonproductive poses. In prospective predictions against seven enzymes, a substrate was identified for five. For one of those cases, a covalent docking prediction, confirmed by empirical screening, and combined with genomic context analysis, suggested the identity of the enzyme that catalyzes the orphan phosphatase reaction in the riboflavin biosynthetic pathway of Bacteroides.