Biocatalytic and structural properties of a highly engineered halohydrin dehalogenase

Origin of Chemoselectivity of Halohydrin Dehalogenase-Catalyzed Epoxide Ring-Opening Reactions

By performing molecular dynamics (MD), quantum mechanical/molecular mechanical (QM/MM) calculations, and QM cluster calculations, the origin of chemoselectivity of halohydrin dehalogenase (HHDH)-catalyzed ring-opening reactions of epoxide with the nucleophilic reagent NO2- has been explored. Four possible chemoselective pathways were considered, and the computed results indicate that the pathway associated with the nucleophilic attack on the Cα position of epoxide by NO2- is most energetically favorable and has an energy barrier of 12.9 kcal/mol, which is close to 14.1 kcal/mol derived from experimental kinetic data. A hydrogen bonding network formed by residues Ser140, Tyr153, and Arg157 can strengthen the electrophilicity of the active site of the epoxide substrate to affect chemoselectivity. To predict the energy barrier trends of the chemoselective transition states, multiple analyses including distortion analysis and electrophilic Parr function (Pk+) analysis were carried out with or without an enzyme environment. The obtained insights should be valuable for the rational design of enzyme-catalyzed and biomimetic organocatalytic epoxide ring-opening reactions with special chemoselectivity.

Biocatalytic Thionation of Epoxides for Enantioselective Synthesis of Thiiranes

Expanding the enzymatic toolbox for the green synthesis of valuable molecules is still of high interest in synthetic chemistry and the pharmaceutical industry. Chiral thiiranes are valuable sulfur-containing heterocyclic compounds, but relevant methods for their enantioselective synthesis are limited. Herein, we report a biocatalytic thionation strategy for the enantioselective synthesis of thiiranes, which was developed based on the halohydrin dehalogenase (HHDH)-catalyzed enantioselective ring-opening reaction of epoxides with thiocyanate and a subsequent nonenzymatic rearrangement process. A novel HHDH was identified and engineered for enantioselective biocatalytic thionation of various aryl- and alkyl-substituted epoxides on a preparative scale, affording the corresponding thiiranes in up to 43 % isolated yield and 98 % ee. Large-scale synthesis and useful transformations of chiral thiiranes were also performed to demonstrate the utility and scalability of the biocatalytic thionation strategy.

Biocatalytic approach to chiral fluoroaromatic scaffolds

Ten different fluorinated aromatic epoxides have been tested as potential substrates for halohydrin dehalogenase (HHDH) HheC. The majority of investigated epoxides are useful building blocks in synthetic chemistry applications, with a number of them being polysubstituted. Moderate to high enantioselectivities (E_R = 15 → 200) were observed in azidolysis, allowing the synthesis of enantioenriched (R)-azido alcohols containing fluorine in the molecule. In the case where a reaction runs over 50% conversion, enantiopure (S)-epoxides are also available. While o-F-styrene oxide was easily converted into a product, a sterically challenging o-CF₃-derivative was not accepted by HheC. In silico probing of the binding site indicates that, in order to accommodate an o-CF₃-derivative in the HheC active site, it is necessary to eliminate steric hindrance. Hence, we extended our research by probing several available HheC variants containing relevant modifications in the active site. The active mutant P84V/F86P/T134A/N176A (named HheC-M4) was identified, showing not only high activity towards o-CF₃-styrene oxide, but also inverted enantioselectivity (E_S = 27). Since (S)-enantioselective HHDHs are rare and therefore valuable for their synthetic application, this enzyme was screened on the initial panel of substrates. The observed (S)-enantioselectivity (E_S = 1-111) is ascribed to the formation of the additional space by introduced mutations in HheC-M4, which is also confirmed by classical MD simulations. Successive molecular docking demonstrated that this newly formed tunnel located close to the protein surface is a critical feature of HheC-M4, representing a novel binding site.

Insights into the molecular determinants of thermal stability in halohydrin dehalogenase HheD2

Halohydrin dehalogenases (HHDHs) are promising enzymes for application in biocatalysis due to their promiscuous epoxide ring-opening activity with various anionic nucleophiles. So far, seven different HHDH subtypes A to G have been reported with subtype D containing the by far largest number of enzymes. Moreover, several characterized members of subtype D have been reported to display outstanding characteristics such as high catalytic activity, broad substrate spectra or remarkable thermal stability. Yet, no structure of a D-type HHDH has been reported to date that could be used to investigate and understand those features on a molecular level. We therefore solved the crystal structure of HheD2 from gamma proteobacterium HTCC2207 at 1.6 Å resolution and used it as a starting point for targeted mutagenesis in combination with molecular dynamics (MD) simulation, in order to study the low thermal stability of HheD2 in comparison with other members of subtype D. This revealed a hydrogen bond between conserved residues Q160 and D198 to be connected with a high catalytic activity of this enzyme. Moreover, a flexible surface region containing two α-helices was identified to impact thermal stability of HheD2. Exchange of this surface region by residues of HheD3 yielded a variant with 10 °C higher melting temperature and reaction temperature optimum. Overall, our results provide important insights into the structure-function relationship of HheD2 and presumably for other D-type HHDHs. DATABASES: Structural data are available in PDB database under the accession number 7B73.

Conformational Landscapes of Halohydrin Dehalogenases and Their Accessible Active Site Tunnels

Halohydrin dehalogenases (HHDH) are industrially relevant biocatalysts exhibiting a promiscuous epoxide-ring opening reactivity in the presence of small nucleophiles, thus giving access to novel carbon–carbon, carbon–oxygen, carbon–nitrogen, and carbon–sulfur bonds. Recently, the repertoire of HHDH has been expanded, providing access to some novel HHDH subclasses exhibiting a broader epoxide substrate scope. In this work, we develop a computational approach based on the application of linear and non-linear dimensionality reduction techniques to long time-scale Molecular Dynamics (MD) simulations to study the HHDH conformational landscapes. We couple the analysis of the conformational landscapes to CAVER calculations to assess their impact on the active site tunnels and potential ability towards bulky epoxide ring opening reaction. Our study indicates that the analyzed HHDHs subclasses share a common breathing motion of the halide binding pocket, but present large deviations in the loops adjacent to the active site pocket and N-terminal regions. Such conformational differences affect the available tunnels for epoxide binding to the active site. The superior activity of the HHDH G subclass towards bulkier substrates is explained by the additional structural elements delimiting the active site region, its rich conformational heterogeneity, and the substantially wider and frequently observed active site tunnels. This study therefore provides key information for HHDH promiscuity and engineering.

Position 123 of halohydrin dehalogenase HheG plays an important role in stability, activity, and enantioselectivity

HheG from Ilumatobacter coccineus is a halohydrin dehalogenase with synthetically useful activity in the ring opening of cyclic epoxides with various small anionic nucleophiles. This enzyme provides access to chiral β-substituted alcohols that serve as building blocks in the pharmaceutical industry. Wild-type HheG suffers from low thermostability, which poses a significant drawback for potential applications. In an attempt to thermostabilize HheG by protein engineering, several single mutants at position 123 were identified which displayed up to 14 °C increased apparent melting temperatures and up to three-fold higher activity. Aromatic amino acids at position 123 resulted even in a slightly higher enantioselectivity. Crystal structures of variants T123W and T123G revealed a flexible loop opposite to amino acid 123. In variant T123G, this loop adopted two different positions resulting in an open or partially closed active site. Classical molecular dynamics simulations confirmed a high mobility of this loop. Moreover, in variant T123G this loop adopted a position much closer to residue 123 resulting in denser packing and increased buried surface area. Our results indicate an important role for position 123 in HheG and give first structural and mechanistic insight into the thermostabilizing effect of mutations T123W and T123G.

Biocatalyzed Synthesis of Statins: A Sustainable Strategy for the Preparation of Valuable Drugs

Statins, inhibitors of 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase, are the largest selling class of drugs prescribed for the pharmacological treatment of hypercholesterolemia and dyslipidaemia. Statins also possess other therapeutic effects, called pleiotropic, because the blockade of the conversion of HMG-CoA to (R)-mevalonate produces a concomitant inhibition of the biosynthesis of numerous isoprenoid metabolites (e.g., geranylgeranyl pyrophosphate (GGPP) or farnesyl pyrophosphate (FPP)). Thus, the prenylation of several cell signalling proteins (small GTPase family members: Ras, Rac, and Rho) is hampered, so that these molecular switches, controlling multiple pathways and cell functions (maintenance of cell shape, motility, factor secretion, differentiation, and proliferation) are regulated, leading to beneficial effects in cardiovascular health, regulation of the immune system, anti-inflammatory and immunosuppressive properties, prevention and treatment of sepsis, treatment of autoimmune diseases, osteoporosis, kidney and neurological disorders, or even in cancer therapy. Thus, there is a growing interest in developing more sustainable protocols for preparation of statins, and the introduction of biocatalyzed steps into the synthetic pathways is highly advantageous—synthetic routes are conducted under mild reaction conditions, at ambient temperature, and can use water as a reaction medium in many cases. Furthermore, their high selectivity avoids the need for functional group activation and protection/deprotection steps usually required in traditional organic synthesis. Therefore, biocatalysis provides shorter processes, produces less waste, and reduces manufacturing costs and environmental impact. In this review, we will comment on the pleiotropic effects of statins and will illustrate some biotransformations nowadays implemented for statin synthesis.

Molecular modification of a halohydrin dehalogenase for kinetic regulation to synthesize optically pure (S)-epichlorohydrin

Asymmetric synthesis of chiral epichlorohydrin (ECH) from 1,3-dichloro-2-propanol (1,3-DCP) using halohydrin dehalogenases (HHDHs) is of great value due to the 100% theoretical yield and high enantioselectivity. The vital problem in the asymmetric synthesis is to prepare optically pure ECH. In this study, key amino acid residues located at halide ion channels of HheC (P175S/W249P) (HheC_PS) were modified to regulate the kinetic parameters. HheC_PS I81W, F86N and V94R were constructed with the corresponding halide ion channels destroyed. The catalytically efficiencies (k_cat/K_m) of the three mutants exhibited 0.38-, 0.23- and 0.23-fold decrease toward (S)-ECH and the reverse reaction was significantly inhibited. As the results, (S)-ECH was synthesized with >99% enantiomeric excess (e.e.) and 63.42%, 67.08% and 57.01% yields, respectively, under 20 mM 1,3-DCP as substrate. To our knowledge, this is the first investigation of the molecule kinetic modification of HHDHs and also the first report for the biosynthesis of optically pure (S)-ECH from 1,3-DCP using HHDHs.

Identification and characterization of a highly S-enantioselective halohydrin dehalogenase from Tsukamurella sp. 1534 for kinetic resolution of halohydrins

Halohydrin dehalogenases are remarkable enzymes which possess promiscuous catalytic activity and serve as potential biocatalysts for the synthesis of chiral halohydrins, epoxides and β-substituted alcohols. The enzyme HheC exhibits a highly R enantioselectivity in the processes of dehalogenation of vicinal halohydrins and ring-opening of epoxides, which attracts more attentions in organic synthesis. Recently dozens of novel potential halohydrin dehalogenases have been identified by gene mining, however, most of the characterized enzymes showed low stereoselectivity. In this study, a novel halohydrin dehalogenase of HheA10 from Tsukamurella sp. 1534 has been heterologously expressed, purified and characterized. Substrate spectrum and kinetic resolution studies indicated the HheA10 was a highly S enantioselective enzyme toward several halohydrins, which produced the corresponding epoxides with the ee (enantiomeric excess) and E values up to >99% and >200 respectively. Our results revealed the HheA10 was a promising biocatalyst for the synthesis of enantiopure aromatic halohydrins and epoxides via enzymatic kinetic resolution of racemic halohydrins. What's more important, the HheA10 as the first individual halohydrin dehalogenase with the highly S enantioselectivity provides a complementary enantioselectivity to the HheC.

Halofunctionalization of alkenes by vanadium chloroperoxidase from Curvularia inaequalis

The vanadium-dependent chloroperoxidase from Curvularia inaequalis is a stable and efficient biocatalyst for the hydroxyhalogenation of a broad range of alkenes into halohydrins. Up to 1 200 000 TON with 69 s^-1 TOF were observed for the biocatalyst. A bienzymatic cascade to yield epoxides as reaction products is presented.

A robust cosolvent-compatible halohydrin dehalogenase by computational library design

To improve the applicability of halohydrin dehalogenase as a catalyst for reactions in the presence of organic cosolvents, we explored a computational library design strategy (Framework for Rapid Enzyme Stabilization by Computational libraries) that involves discovery and in silico evaluation of stabilizing mutations. Energy calculations, disulfide bond predictions and molecular dynamics simulations identified 218 point mutations and 35 disulfide bonds with predicted stabilizing effects. Experiments confirmed 29 stabilizing point mutations, most of which were located in two distinct regions, whereas introduction of disulfide bonds was not effective. Combining the best mutations resulted in a 12-fold mutant (HheC-H12) with a 28°C higher apparent melting temperature and a remarkable increase in resistance to cosolvents. This variant also showed a higher optimum temperature for catalysis while activity at low temperature was preserved. Mutant H12 was used as a template for the introduction of mutations that enhance enantioselectivity or activity. Crystal structures showed that the structural changes in the H12 mutant mostly agreed with the computational predictions and that the enhanced stability was mainly due to mutations that redistributed surface charges and improved interactions between subunits, the latter including better interactions of water molecules at the subunit interfaces.

Enzymatic Halogenation and Dehalogenation Reactions: Pervasive and Mechanistically Diverse

Naturally produced halogenated compounds are ubiquitous across all domains of life where they perform a multitude of biological functions and adopt a diversity of chemical structures. Accordingly, a diverse collection of enzyme catalysts to install and remove halogens from organic scaffolds has evolved in nature. Accounting for the different chemical properties of the four halogen atoms (fluorine, chlorine, bromine, and iodine) and the diversity and chemical reactivity of their organic substrates, enzymes performing biosynthetic and degradative halogenation chemistry utilize numerous mechanistic strategies involving oxidation, reduction, and substitution. Biosynthetic halogenation reactions range from simple aromatic substitutions to stereoselective C-H functionalizations on remote carbon centers and can initiate the formation of simple to complex ring structures. Dehalogenating enzymes, on the other hand, are best known for removing halogen atoms from man-made organohalogens, yet also function naturally, albeit rarely, in metabolic pathways. This review details the scope and mechanism of nature's halogenation and dehalogenation enzymatic strategies, highlights gaps in our understanding, and posits where new advances in the field might arise in the near future.

Exploring the thermostable properties of halohydrin dehalogenase from Agrobacterium radiobacter AD1 by a combinatorial directed evolution strategy

As a crucial factor for biocatalysts, protein thermostability often arises from a combination of factors that are often difficult to rationalize. In this work, the thermostable nature of halohydrin dehalogenase from Agrobacterium radiobacter AD1 (HheC) was systematically explored using a combinatorial directed evolution approach. For this, a mutagenesis library of HheC mutants was first constructed using error-prone PCR with low mutagenesis frequency. After screening approximately 2000 colonies, six mutants with eight mutation sites were obtained. Those mutation sites were subsequently combined by adopting several rounds of iterative saturation mutagenesis (ISM) approach. After four rounds of saturation mutagenesis, one best mutant ISM-4 with a 3400-fold improvement in half-life (t _1/2) inactivation at 65 °C, 18 °C increase in apparent T _m value, and 20 °C increase in optimum temperature was obtained, compared to wild-type HheC. To the best of our knowledge, the mutant represents the most thermostable HheC variant reported up to now. Moreover, the mutant was as active as wild-type enzyme for the substrate 1,3-dichloro-2-propanol, and they remained most enantioselectivity of wild-type enzyme in the kinetic resolution of rac-2-chloro-1-phenolethanol, exhibiting a great potential for industrial applications. Our structural investigation highlights that surface loop regions are hot spots for modulating the thermostability of HheC, the residues located at these regions contribute to the thermostability of HheC in a cooperative way, and protein rigidity and oligomeric interface connections contribute to the thermostability of HheC. All of these essential factors could be used for further design of an even more thermostable HheC, which, in turn, could greatly facilitate the application of the enzyme as a biocatalyst.

Discovery of a Xylooligosaccharide Oxidase from Myceliophthora thermophila C1

By inspection of the predicted proteome of the fungus Myceliophthora thermophila C1 for vanillyl-alcohol oxidase (VAO)-type flavoprotein oxidases, a putative oligosaccharide oxidase was identified. By homologous expression and subsequent purification, the respective protein could be obtained. The protein was found to contain a bicovalently bound FAD cofactor. By screening a large number of carbohydrates, several mono- and oligosaccharides could be identified as substrates. The enzyme exhibits a strong substrate preference toward xylooligosaccharides; hence it is named xylooligosaccharide oxidase (XylO). Chemical analyses of the product formed upon oxidation of xylobiose revealed that the oxidation occurs at C1, yielding xylobionate as product. By elucidation of several XylO crystal structures (in complex with a substrate mimic, xylose, and xylobiose), the residues that tune the unique substrate specificity and regioselectivity could be identified. The discovery of this novel oligosaccharide oxidase reveals that the VAO-type flavoprotein family harbors oxidases tuned for specific oligosaccharides. The unique substrate profile of XylO hints at a role in the degradation of xylan-derived oligosaccharides by the fungus M. thermophila C1.

Recent advances on halohydrin dehalogenases-from enzyme identification to novel biocatalytic applications

Halohydrin dehalogenases are industrially relevant enzymes that catalyze the reversible dehalogenation of vicinal haloalcohols with formation of the corresponding epoxides. In the reverse reaction, also other negatively charged nucleophiles such as azide, cyanide, or nitrite are accepted besides halides to open the epoxide ring. Thus, novel C-N, C-C, or C-O bonds can be formed by halohydrin dehalogenases, which makes them attractive biocatalysts for the production of various β-substituted alcohols. Despite the fact that only five individual halohydrin dehalogenase enzyme sequences have been known until recently enabling their heterologous production, a large number of different biocatalytic applications have been reported using these enzymes. The recent characterization of specific sequence motifs has facilitated the identification of novel halohydrin dehalogenase sequences available in public databases and has largely increased the number of recombinantly available enzymes. These will help to extend the biocatalytic repertoire of this enzyme family and to foster novel biotechnological applications and developments in the future. This review gives a general overview on the halohydrin dehalogenase enzyme family and their biochemical properties and further focuses on recent developments in halohydrin dehalogenase biocatalysis and protein engineering.

A high-throughput adrenaline test for the exploration of the catalytic potential of halohydrin dehalogenases in epoxide ring-opening reactions

The principle of the adrenaline test for enzymes is based on the quantification of periodate-sensitive reaction products with adrenaline to produce a chromogenic compound adrenochrome that can be easily detected. Here, a rapid whole-cell -based adrenaline assay for the activity measurement of halohydrin dehalogenases (HHDHs) in nucleophile-mediated epoxide ring-opening reactions is presented. The assay was validated using two types of model reactions (glycidol with nucleophiles and nitrite with epoxides). Moreover, the reliability of the assay was confirmed by gas chromatography analysis. Our results demonstrated that the developed assay is efficient in both library screening and the evaluation of catalytic diversity and specificity of HHDHs. Thus, the assay represents a valuable tool in the evolution of HHDHs for its industrial applications. Moreover, the adrenaline test exhibits a great potential for enzyme assay and could be easily adopted for other suitable enzymes.

MDC-Analyzer-facilitated combinatorial strategy for improving the activity and stability of halohydrin dehalogenase from Agrobacterium radiobacter AD1

Halohydrin dehalogenase from Agrobacterium radiobacter AD1 (HheC) displays a broad substrate range with high regio- and enantioselectivity of both ring-closure and ring-opening reactions, making the enzyme a useful catalyst for the production of optically pure epoxides and β-substituted alcohols. In this study, we report a novel method using an MDC-Analyzer-facilitated combinatorial strategy to improve the activity and stability of HheC by simultaneously randomizing multiple contiguous residues. Six contiguous active-site residues, which are the hotspots for improving the activity of HheC, were simultaneously selected and randomized using the MDC-Analyzer-facilitated combinatorial strategy, resulting in a high-quality mutagenesis library. After screening a total of 1152 clones, three positive mutants were obtained, which exhibited approximately 3.5-5.9-fold higher kcat values than the wild-type HheC toward 1,3-dichloro-2-propanol (1,3-DCP). However, the inactivation half-life of the best mutant (DG9) at 55 °C decreased 9-fold compared with that of the wild-type HheC. To improve the stability of mutant DG9, seven contiguous potential surface amino acids were revealed by using the B-FITTER tool. Two charged amino acids, Glu and Lys, which are more abundant in thermophilic proteins than in their mesophilic counterparts, were selected to substitute those seven amino acids and were combined together via an MDC-Analyzer-facilitated combinatorial strategy. Two mutants displaying 1.6- and 2.3-fold higher half-life τ1/2 (55 °C) values than their DG9 template were obtained after screening only 384 clones. The results indicated that an MDC-Analyzer-facilitated combinatorial strategy represents an efficient tool for the directed evolution of functional enzymes with multiple contiguous targeting sites.

X-ray crystallographic validation of structure predictions used in computational design for protein stabilization

Protein engineering aimed at enhancing enzyme stability is increasingly supported by computational methods for calculation of mutant folding energies and for the design of disulfide bonds. To examine the accuracy of mutant structure predictions underlying these computational methods, crystal structures of thermostable limonene epoxide hydrolase variants obtained by computational library design were determined. Four different predicted effects indeed contributed to the obtained stabilization: (i) enhanced interactions between a flexible loop close to the N-terminus and the rest of the protein; (ii) improved interactions at the dimer interface; (iii) removal of unsatisfied hydrogen bonding groups; and (iv) introduction of additional positively charged groups at the surface. The structures of an eightfold and an elevenfold mutant showed that most mutations introduced the intended stabilizing interactions, and side-chain conformations were correctly predicted for 72-88% of the point mutations. However, mutations that introduced a disulfide bond in a flexible region had a larger influence on the backbone conformation than predicted. The enzyme active sites were unaltered, in agreement with the observed preservation of catalytic activities. The structures also revealed how a c-Myc tag, which was introduced for facile detection and purification, can reduce access to the active site and thereby lower the catalytic activity. Finally, sequence analysis showed that comprehensive mutant energy calculations discovered stabilizing mutations that are not proposed by the consensus or B-FIT methods.

Exploring the enantioselective mechanism of halohydrin dehalogenase from Agrobacterium radiobacter AD1 by iterative saturation mutagenesis

Halohydrin dehalogenase from Agrobacterium radiobacter AD1 (HheC) shows great potential in producing valuable chiral epoxides and β-substituted alcohols. The wild-type (WT) enzyme displays a high R-enantiopreference toward most aromatic substrates, whereas no S-selective HheC has been reported to date. To obtain more enantioselective enzymes, seven noncatalytic active-site residues were subjected to iterative saturation mutagenesis (ISM). After two rounds of screening aspects of both activity and enantioselectivity (E), three outstanding mutants (Thr134Val/Leu142Met, Leu142Phe/Asn176His, and Pro84Val/Phe86Pro/Thr134Ala/Asn176Ala mutants) with divergent enantioselectivity were obtained. The two double mutants displayed approximately 2-fold improvement in R-enantioselectivity toward 2-chloro-1-phenylethanol (2-CPE) without a significant loss of enzyme activity compared with the WT enzyme. Strikingly, the Pro84Val/Phe86Pro/Thr134Ala/Asn176Ala mutant showed an inverted enantioselectivity (from an ER of 65 [WT] to an ES of 101) and approximately 100-fold-enhanced catalytic efficiency toward (S)-2-CPE. Molecular dynamic simulation and docking analysis revealed that the phenyl side chain of (S)-2-CPE bound at a different location than that of its R-counterpart; those mutations generated extra connections for the binding of the favored enantiomer, while the eliminated connections reduced binding of the nonfavored enantiomer, all of which could contribute to the observed inverted enantiopreference.

Enzymatic approaches to the preparation of chiral epichlorohydrin

Enantiomerically pure epichlorohydrin is a key chiral synthon in the preparation of 4-chloro-3-hydroxybutyrate, pheromones,l-carnitine, and β-adrenergic blockers.

Expanding the Halohydrin Dehalogenase Enzyme Family: Identification of Novel Enzymes by Database Mining

Halohydrin dehalogenases are very rare enzymes that are naturally involved in the mineralization of halogenated xenobiotics. Due to their catalytic potential and promiscuity, many biocatalytic reactions have been described that have led to several interesting and industrially important applications. Nevertheless, only a few of these enzymes have been made available through recombinant techniques; hence, it is of general interest to expand the repertoire of these enzymes so as to enable novel biocatalytic applications. After the identification of specific sequence motifs, 37 novel enzyme sequences were readily identified in public sequence databases. All enzymes that could be heterologously expressed also catalyzed typical halohydrin dehalogenase reactions. Phylogenetic inference for enzymes of the halohydrin dehalogenase enzyme family confirmed that all enzymes form a distinct monophyletic clade within the short-chain dehydrogenase/reductase superfamily. In addition, the majority of novel enzymes are substantially different from previously known phylogenetic subtypes. Consequently, four additional phylogenetic subtypes were defined, greatly expanding the halohydrin dehalogenase enzyme family. We show that the enormous wealth of environmental and genome sequences present in public databases can be tapped for in silico identification of very rare but biotechnologically important biocatalysts. Our findings help to readily identify halohydrin dehalogenases in ever-growing sequence databases and, as a consequence, make even more members of this interesting enzyme family available to the scientific and industrial community.