Kinetic Mechanism of the Dechlorinating Flavin-dependent Monooxygenase HadA*

The accumulation of chlorophenols (CPs) in the environment, due to their wide use as agrochemicals, has become a serious environmental problem. These organic halides can be degraded by aerobic microorganisms, where the initial steps of various biodegradation pathways include an oxidative dechlorinating process in which chloride is replaced by a hydroxyl substituent. Harnessing these dechlorinating processes could provide an opportunity for environmental remediation, but detailed catalytic mechanisms for these enzymes are not yet known. To close this gap, we now report transient kinetics and product analysis of the dechlorinating flavin-dependent monooxygenase, HadA, from the aerobic organism Ralstonia pickettii DTP0602, identifying several mechanistic properties that differ from other enzymes in the same class. We first overexpressed and purified HadA to homogeneity. Analyses of the products from single and multiple turnover reactions demonstrated that HadA prefers 4-CP and 2-CP over CPs with multiple substituents. Stopped-flow and rapid-quench flow experiments of HadA with 4-CP show the involvement of specific intermediates (C4a-hydroperoxy-FAD and C4a-hydroxy-FAD) in the reaction, define rate constants and the order of substrate binding, and demonstrate that the hydroxylation step occurs prior to chloride elimination. The data also identify the non-productive and productive paths of the HadA reactions and demonstrate that product formation is the rate-limiting step. This is the first elucidation of the kinetic mechanism of a two-component flavin-dependent monooxygenase that can catalyze oxidative dechlorination of various CPs, and as such it will serve as the basis for future investigation of enzyme variants that will be useful for applications in detoxifying chemicals hazardous to human health.

Chlorophenols (CPs) 2 such as 2,4,6-trichlorophenol (2,4,6-TCP), 2,4-dichlorophenol (2,, and pentachlorophenol (PCP) are commonly used as herbicides, insecticides, and fungicides in agricultural industries. CPs can also be derived from pesticide degradation. For example, 2,4,5-trichlorophenol (2,4,5-TCP) is derived from the degradation of 2,4,5-trichlorophenoxyacetate (1), and 2-chloro-4-bromo phenol is derived from the organophosphate insecticide, profenofos (2). These compounds are classified as persistent environmental pollutants that can accumulate in water and soil (3), and most of them are registered under the priority pollutants list (4) and/or the priority list of hazardous substances (5). High human exposure to these compounds can lead to severe health problems such as cancer, DNA damage, organ damage, or even acute death (3, 6 -8). As these compounds are found widespread in the environment due to their utilization in agriculture and in industrial processes, a biodegradation process that can provide a sustainable method to detoxify them is essential to the reduction of CP contamination in fresh water and soil.
Aerobic biodegradation of CPs generally involves replacement of a chloride group by a hydroxyl group in the initial steps of these biodegradation pathways. This oxidative dehalogenation is often catalyzed by flavin-dependent monooxygenases (20). These oxygenases include both types of monooxygenases, those that utilize a single protein as well as those that utilize two proteins. The most well studied single-component flavin-dependent dechlorinase is PcpB from S. chlorophenolicum. This enzyme binds to FAD as a cofactor and requires NADPH and O 2 as co-substrates for the conversion of PCP to 2,3,5,6-tetrachlorobenzoquinone in the initial step of PCP degradation. The quinone product can be chemically reduced to 2,3,5,6-tetrachloro-p-hydoquinone (18,19). However, turnover of PcpB using PCP as a substrate is very slow, implying that PCP may not be a native substrate of PcpB (18). Many oxidative dechlorination reactions are also catalyzed by two-component flavindependent enzymes, including TftD from B. cepacia AC1100 (9,10), TcpA from R. eutropha JMP134 (13), and HadA from R. pickettii DTP0602 (15). Although a few of these enzymes have been shown to catalyze oxidative dechlorination of CPs, their reaction mechanisms are not well understood. None of these enzymes have been explored with regard to their kinetic mechanisms, and the rate constants associated with the individual steps have yet to be measured. Only X-ray structures of TftD and TcpA and the overall bioconversion activities of these enzymes have been reported (9,10,13,15,21).
Based on previous reports, the substrate utilization of HadA (14,15) seems to be broader than that of TcpA (specific to 2,4,6-TCP; 86.7% identity to HadA) and TftD (specific for 2,4,5-TCP; 64.2% identity to HadA), making HadA attractive for development as a biocatalyst for CP detoxification. Previously, studies of the HadA reaction were only carried out with whole cells and lysates of cells overexpressing HadA (14,15) without measuring the products that were specifically generated from the flavin-bound enzyme. As study of the transient kinetics of a two-component flavin-dependent dechlorinase has not been reported, we have chosen HadA as a model for investigating the substrate utilization and rate constants associated with the individual steps under single turnover conditions to identify the flavin species and reaction steps and pathways involved in oxidative dechlorination.
Based on sequence similarity, HadA belongs to the family of two-component flavin-dependent monooxygenases that utilize two separate proteins in their catalysis (class D monooxygenase) (22). A reductase component of these enzymes catalyzes flavin reduction, and an oxygenase component catalyzes substrate oxygenation (22,23). The most well studied enzyme in this family is p-hydroxyphenylacetate hydroxylase (HPAH) from Acinetobacter baumannii. The mechanisms of reduced flavin generation by the reductase (C 1 ), as well as the mechanism of the oxygenase (C 2 ) reaction, have been determined for this enzyme. Furthermore, the individual rate constants associated with each step and the mechanism of reduced flavin transfer between the two components have been reported (24 -30). This knowledge has been useful for applying HPAH for the synthesis of bioactive phenolic acids and in the hydroxylation of aniline derivatives (31,32).
When the sequence of HadA was examined by sequence similarity network (SSN) analysis using the Enzyme Function Initiative (EFI) program to determine its relationship with any homologous enzymes known to date, the results (supplemental Fig. S1) showed that there are 2,839 enzymes in which the sequences have Ͼ34% identity with HadA, and only 11 of these 2,839 enzymes were reviewed by UniProt. The enzymes identified to be homologous to HadA were mostly flavin-dependent monooxygenases that incorporate a hydroxyl group into phenols such as CPs or nitrophenols (9,10,(33)(34)(35). As these enzymes showed potential applications in bioremediation of toxic compounds, the understanding of kinetic mechanisms associated with a flavin-dependent dechlorinase (HadA) reported herein should also be useful for giving an overall picture about the reaction mechanism of these enzymes in general.
In this work, the enzyme from R. pickettii DTP0602 (HadA) was overexpressed and purified to homogeneity, and its reaction with several CPs was explored. The reactions of HadA with 4-CP (the most optimal substrate) were investigated using stopped-flow and rapid-quench flow techniques to measure the rate constants associated with the individual steps. The results clearly determined the order of substrate binding and the rate constants associated with each of the individual steps and identified C4a-hydroperoxy-FAD and C4a-hydroxy-FAD as intermediates in the dechlorination process catalyzed by HadA.

Expression and Purification of HadA
Recombinant HadA was overexpressed and purified according to the purification methods described under "Experimental Procedures" to yield protein with Ͼ95% purity (supplemental Fig. S2). The protein preparation protocol resulted in ϳ0.6 g of purified HadA from 3.9 liters of culture. HadA had no cofactor bound, and the native molecular mass suggests that HadA is a tetrameric protein with a subunit molecular mass of ϳ59 kDa (supplemental Fig. S3).

Flavin Specificity of HadA
As HadA is the oxygenase component of a two-component flavin-dependent monooxygenase, we first needed to identify its flavin specificity and its ability to form C4a-hydroperoxyflavin. The enzymes in this class use specific reduced flavins as substrates and form C4a-hydroperoxyflavin as a reactive intermediate to perform monooxygenation plus other additional activities (22,23,36,37). Binding of FADH Ϫ or FMNH Ϫ to HadA was investigated to identify which flavin could bind to HadA and yield a binary complex able to react with oxygen. Solutions of FADH Ϫ and FMNH Ϫ were mixed with buffers containing oxygen (air saturation) in the presence or absence of HadA. The results indicate that only the oxidation of FADH Ϫ was affected by the presence of HadA, indicating that the enzyme can only bind to FADH Ϫ and not FMNH Ϫ (supplemental Fig. S4). The reaction of free FADH Ϫ with oxygen showed an increase in absorbance at 380 and 450 nm with an observed rate constant of ϳ2.3 s Ϫ1 (red lines, supplemental Fig.  S4A). In the presence of HadA, the oxidation of FADH Ϫ resulted in two observable phases. The first phase consisted of a large absorbance increase at 380 nm (blue solid line, supplemental Fig. S4A), and the second phase was characterized by a large and slow (observed rate constant of 0.005 s Ϫ1 ) increase in absorbance at 450 nm (blue dashed line, supplemental Fig.  S4A). These data indicate that in the presence of HadA, FADH Ϫ oxidation resulted in the formation of C4a-hydroperoxy-FAD, which can be detected by a large absorbance increase at 380 nm where C4a-hydroperoxy-FAD absorbs maximally (more analysis to follow). The mixing of FMNH Ϫ with air-saturated HadA (blue lines, supplemental Fig. S4B) yielded the same kinetic traces as those of the reaction without HadA (red lines, supplemental Fig. S4B). Both reactions showed the increase in absorbance at 380 and 446 nm with an observed rate constant of 3.2 s Ϫ1 . These data clearly suggest that HadA can only bind to FADH Ϫ and not FMNH Ϫ .

Identification of Product from HadA Single Turnover Reactions
To identify the substrate specificity of HadA, the products resulting from single turnover reactions of HadA (see "Experimental Procedures") were identified by comparing their HPLC and mass spectroscopic profiles with those of standard compounds. As most of the products were expected to be quinone derivatives, ascorbic acid (1 mM) was added into the reaction to convert them into their hydroquinone form, which is more stable. The results of the product identification studies (Table 1 and supplemental Fig. S5, A-C) show that HadA catalyzes the hydroxylation of various CPs, including 4-chlorophenol (4-CP), 2-chlorophenol (2-CP), 2,4-dichlorophenol (2,4-DCP), 2,4,5trichlorophenol (2,4,5-TCP), and 2,4,6-trichlorophenol (2,4,6-TCP) to produce hydroquinone (HQ), chlorohydroquinone (CHQ), 2,5-dichlorohydroquinone (2,5-DCHQ), and 2,6-dichlorohydroquinone (2,6-DCHQ), respectively. The product identification results are summarized in Table 1. Based on these reactions, HadA incorporates a hydroxyl group at the paraposition of the original hydroxyl group of all five substrates investigated. When the para-position contains a chloride substituent (in the case of 4-CP, 2,4-DCP, 2,4,5-TCP, and 2,4,6-TCP), HadA catalyzes hydroxylation plus Cl Ϫ elimination ( Fig.  1). Comparison of the percentage of product formation after single turnover reactions for the different substrates indicates that HadA prefers to use mono-substituted CPs over multiplesubstituted CPs. These results contradict previously reported data from cell lysate experiments (where HadA was not purified) in which 2,4,6-TCP was reported as the most preferred substrate for HadA (14,15).

Preference of Hydroxylation Site Identified by Multiple Turnover Reactions of HadA and 2,4,5-TCP
The results in Fig. 2C also show that the product from 2,4,5-TCP multiple turnover reactions, 2,5-DCHQ, decreased over time, indicating that the compound may be further hydroxylated. To characterize the nature of the second conversion step of 2,4,5-TCP, the reaction of HadA with 2,5-DCHQ was investigated. The results showed that 2,5-DCHQ was completely depleted within 4 h (supplemental Fig. S6A), with concomitant appearance of a hydroxylation product that has an m/z of ϩ16 atomic units compared with the value for 2,5-DCHQ. This compound also has mass spectroscopic and HPLC profiles consistent with those of 5-chlorohydroxyquinol (5-CHQL) (supplemental Fig. S6, A and B). All results suggest that HadA can convert 2,5-DCHQ to 5-CHQL. Therefore, the overall reaction of HadA and 2,4,5-TCP can be summarized as in Fig. 3 in which the first hydroxylation and Cl Ϫ elimination step occurs at position 4 and then the second hydroxylation and Cl Ϫ elimination step proceeds at position 2.

pH Optimization of HadA
To identify the optimum pH for the HadA reaction to be used in transient kinetics experiments, single turnover reactions of HadA and its best substrate, 4-CP, were carried out at various pH values, and the amount of product formation was analyzed according to the protocol described under "Experimental Procedures." An anaerobic solution of HadA/FADH Ϫ in 20 mM For chlorophenols with a para-substituent, HadA catalyzes hydroxylation plus Cl Ϫ elimination to yield quinone derivatives that can be reduced by ascorbic acid to form hydroquinone products (upper pathway). In the absence of a para-substituent, HadA catalyzes hydroxylation alone to yield hydroquinone (lower pathway). a Product was characterized in the presence of ascorbic acid. b Product formation from single turnover reactions was estimated based on the depletion of substrate. A dash represents reactions that were not performed. ND (not determined), t1 ⁄ 2 of 2,4,6-TCP could not be analyzed because the reaction did not finish after 600 min.
HEPES, pH 7.5, was mixed with an equal volume of an aerobic solution of 4-CP at various pH values from 6.5 to 10.5. The HQ production results (Fig. 4) indicated that the reaction of HadA has an optimum pH within the range of ϳ7.0 -8.2. Therefore, the transient kinetics experiments were all performed at pH 7.5. Transient kinetic experiments were also performed at various pH values to explore if pH affects the ability of HadA to form C4a-hydroperoxy-FAD (shown under Section F).

Effect of Ascorbic Acid on the HadA Reaction
As most of the reactions in this report contained ascorbic acid (1 mM) to stabilize the hydroquinone products, the effect of ascorbic acid on the reaction of HadA was investigated. The transient kinetics of the HadA reactions in the absence or presence of ascorbic acid was investigated using stopped-flow techniques. A solution of the HadA⅐FADH Ϫ binary complex was mixed with air-saturated 20 mM HEPES, pH 7.5, in the presence or absence of 1 mM ascorbic acid on the stopped-flow spectrophotometer. Kinetic traces of both experiments were almost identical (supplemental Fig. S7), suggesting that ascorbic acid does not interfere with HadA catalytic behavior and that it can be used for stabilizing the substrate and product in the HadA reaction.

Reaction Mechanism of HadA with 4-CP as a Substrate Investigated by Transient Kinetics
A: Binding of FADH Ϫ (Step 1 in Fig. 5; FADH Ϫ Is the First Substrate to Bind to HadA)-For two-component flavin-dependent monooxygenases that use reduced flavin as a substrate such as HadA, the enzyme must bind to FADH Ϫ to form an active binary complex HadA⅐FADH Ϫ before it can react with oxygen to form C4a-hydroperoxy-FAD.  Fig. 6). This is much slower than the oxidation rate of free FADH Ϫ (t1 ⁄ 2 ϳ0.3 s, dashed line in Fig. 6). When the binding of HadA and FADH Ϫ was not complete, the traces showed two observable phases that directly reflected the oxidation of two distinct species of free and HadAbound FADH Ϫ (Fig. 6). Therefore, at t ϭ 2 s, the difference between the absorbance at 450 nm of free FADH Ϫ and that of the HadA-bound FADH Ϫ oxidation represents the amount of HadA⅐FADH Ϫ binary complex formed. A plot of these values versus free HadA concentration represents a binding isotherm that is consistent with the K d value of 2.0 Ϯ 0.2 M (Fig. 6, inset). This result also established the experimental conditions required to achieve complete binding of FADH Ϫ by HadA (a ratio of HadA/FADH Ϫ as 3:1).

A.2. Kinetics of HadA and FADH Ϫ Binding (Step 1 in Fig. 5)-
The binding rate of FADH Ϫ to HadA was measured using double-mixing stopped-flow experiments by monitoring the progress of HadA⅐FADH Ϫ formation. An anaerobic solution of FADH Ϫ (25 M) was mixed with an anaerobic solution of HadA (75 M) in the first mixing step and was incubated for various periods of time (0.01-200 s) before oxygenated buffer (air saturation) was added in the second mixing step. Kinetic traces were obtained by monitoring the change of absorbance at 450 nm. A longer preincubation time prior to adding the oxygenated buffer resulted in slower kinetics of a 450-nm increase in absorbance, raising the t1 ⁄ 2 from ϳ0.3 s (free FADH Ϫ oxidation) to ϳ100 s (HadA-bound FADH Ϫ reaction). This is because formation of the HadA⅐FADH Ϫ complex is more complete (as explained above, the HadA⅐FADH Ϫ binary complex results in  formation of C4a-hydroperoxy-FAD, which eliminates H 2 O 2 slower than the oxidation of free FADH Ϫ (Section A.1)). Results (supplemental Fig. S8A) showed that the behavior of free FADH Ϫ oxidation only disappeared when the incubation time of the first mixing was longer than 20 s. A plot of the absorbance changes at 450 nm at 2 s (representing the amount of HadA⅐FADH Ϫ binary complex formed) versus the incubation age time was analyzed according to a single exponential equation to yield an observed rate constant of 0.33 Ϯ 0.04 s Ϫ1 , which is equivalent to a bimolecular rate constant for the binding of HadA and FADH Ϫ of 4.4 Ϯ 0.5 ϫ 10 3 M Ϫ1 ⅐s Ϫ1 (inset in supplemental Fig. S8A). Fig. 5)-To obtain a spectrum of C4a-hydroperoxy-FAD and to characterize the kinetics of the reaction, a solution of the HadA⅐FADH Ϫ binary complex was mixed with various concentrations of oxygen in 20 mM HEPES, pH 7.5. Analysis of the results (top inset in Fig. 7) indicates that kinetics of these reactions are composed of three phases. The first (0.01-0.2 s) and second (0.2-6 s) phases showed large increases in the absorbance at 380 nm without significant changes at 450 nm. Both phases were dependent on the oxygen concentration. Plots of the observed rate constants (k obs ) of both phases versus oxygen concentrations were linear, consistent with second-order rate constants of 6.4 Ϯ 0.4 ϫ 10 4 and 3.3 Ϯ 0.5 ϫ 10 3 M Ϫ1 ⅐s Ϫ1 (bottom inset in Fig. 7), respectively. The dependence of both phases on oxygen concentration suggests that the HadA⅐ FADH Ϫ complex probably exists as a mixture of fast and slow reacting species (Step 2 in Fig. 5). These two phases represent the reaction of the HadA⅐FADH Ϫ complex with oxygen to form C4a-hydroperoxy-FAD because the spectrum at the end of each phase (solid line with filled circles in Fig. 7) is similar to the spectrum of C4a-hydroperoxyflavin that is produced in the reactions of the oxygenase component of p-hydroxyphenylacetate hydroxylase (25) and luciferase (38). The last phase showed a decrease in absorbance at 380 nm and a large increase in absorbance at 450 nm. The observed rate constant was independent of oxygen concentration. This phase represents a slow intermediate decay by H 2 O 2 elimination with a rate constant of 5.4 Ϯ 0.2 ϫ 10 Ϫ3 s Ϫ1 to form the final species, oxidized FAD (Step 8 in Fig. 5). These data indicate that in the absence of an aromatic substrate, the overall reaction of the binary complex of HadA⅐FADH Ϫ with oxygen can be summarized as in Steps 3 and 8 in Fig. 5. Fig. 5)-We noted that a mixing of an aerobic HadA solution with an anaerobic FADH Ϫ solution resulted in C4a-hydroperoxy-FAD formation with a monophasic kinetic trace (supplemental Fig. S4A), which is different from the results in the top inset of Fig. 7 that shows biphasic formation of C4a-hydroperoxy-FAD. These results suggest that pre-incubation of HadA and FADH Ϫ may lead to a conformational change in the enzyme, in which one form of the HadA⅐FADH Ϫ complex reacts faster with oxygen than the other. To determine whether the pre-incubation of HadA with FADH Ϫ did indeed cause this effect, double-mixing stoppedflow experiments similar to those in Section A.2 were carried out in which the incubation time of the first mixing step to form the HadA⅐FADH Ϫ complex was varied. In the second mixing step, an oxygenated solution was added to form C4a-hydroperoxy-FAD, which could then be detected at 380 nm (supplemental Fig. S8B). The results in supplemental Fig. S8B indicate that at an incubation time of Ͻ20 s, the kinetics of C4a-hydroperoxy-FAD formation was monophasic with a bimolecular constant of 9.9 Ϯ 0.1 ϫ 10 4 M Ϫ1 ⅐s Ϫ1 . When the incubation time was increased to greater than 20 s, biphasic kinetics of C4a-hydroperoxy-FAD formation was observed. These results indicate that the biphasic kinetics observed when mixing a preformed binary complex of HadA⅐FADH Ϫ with oxygen was indeed due to a conformational change in the enzyme to result in two forms of the complex, HadA⅐FADH Ϫ (fast reacting species) and HadA*⅐FADH Ϫ (slow reacting species). The rate constant for this enzymatic conformational change was calculated from a plot of ⌬A 380 nm at t ϭ 2 s (the decrease in absorbance at 380 represents conversion of the enzyme to the slow reacting species) versus the incubation time as 0.14 Ϯ 0.01 s Ϫ1 (inset in supplemental Fig. S8B).

C: Binding of 4-CP and HadA (Step 10 in Fig. 5, 4-CP Can Bind to HadA but Does Not Lead to a Productive Path)
C.1. 4-CP Can Also Bind to ApoHadA and This Binding Prevents the Binding of FADH Ϫ -Although formation of the HadA⅐FADH Ϫ binary complex in a two-component flavin-dependent monooxygenase is required before the enzyme can react with oxygen, it is known that certain substrates can also bind to some enzymes to form a dead-end complex before the binding of reduced flavin (38 -40). To identify whether 4-CP can bind to apoHadA, solutions of air-saturated HadA (75 M) that were pre-incubated with various concentrations of 4-CP (0.1-10 mM) were mixed with an anaerobic solution of FADH Ϫ (25 M) to monitor the formation of C4a-hydroperoxy-FAD using a stopped-flow spectrophotometer. Final concentrations are given as after mixing. If a complex of HadA⅐4-CP can form, the reaction of enzyme-bound FADH Ϫ and oxygen would react differently from the reaction in the absence of 4-CP. Therefore, the kinetic traces of the reactions were measured at 380 nm (Fig. 8A) and 450 nm (Fig. 8B) to monitor whether or not the presence of 4-CP perturbs the formation of C4a-hydroperoxy-FAD. The results clearly show that upon increasing the 4-CP concentration, the amount of C4a-hydroperoxy-FAD formed decreased, and the kinetics of the reaction became more similar to that of the reaction of free FADH Ϫ with oxygen (not bound to HadA). At the highest 4-CP concentration used (orange line in Fig. 8, A and B), the reaction was practically the same with the reaction of free FADH Ϫ reacting with oxygen (dashed line in Fig. 8, A and B). These data clearly demonstrate that 4-CP can bind to HadA, and this binding prevents the binding of FADH Ϫ to HadA, i.e. 4-CP can bind to HadA to form a dead-end complex that cannot proceed through a produc-tive path (Step 10 in Fig. 5). In another case in which HadA was pre-incubated with FADH Ϫ to form the HadA⅐FADH Ϫ binary complex prior to mixing with 4-CP (0.1-10 mM), it was found that at concentrations of 4-CP Ͻ2 mM, C4a-hydroperoxy-FAD can still form to the same extent as in the reaction in the absence of 4-CP (supplemental Fig. S9A). This implies that once the HadA⅐FADH Ϫ complex is formed, 4-CP cannot compete off the FADH Ϫ binding on HadA at concentrations Ͻ2 mM.

C.2. Formation of 4-CP and HadA Leads to Formation of a Dead-end Complex That Does Not Result in Dechlorination-
To confirm that formation of the HadA⅐4-CP complex (in Section C.1) is not relevant to a productive path with dechlorination activity, the amounts of product (HQ) formed from two different mixing experiments, in which HadA⅐4-CP was preformed or not preformed (Fig. 8, A and B and supplemental Fig.  S9), were measured and compared. The results (Fig. 8A, inset) clearly show that higher amounts of 4-CP (Ͼ2 mM) reduced the amount of product formation in both mixing methods, indicating that 4-CP exerts a substrate inhibition effect. However, mixing of a preformed complex of HadA⅐4-CP plus oxygen with FADH Ϫ resulted in a much lower amount of product formation than the mixing of a binary complex of HadA⅐FADH Ϫ with oxygen plus 4-CP (Fig. 8A, inset). At 4-CP concentrations of Ͼ3.2 mM, the preformed complex of HadA⅐4-CP did not give any dechlorinated product. This result also agrees with the kinetic traces because at this concentration of 4-CP, FADH Ϫ could not bind to HadA and only free FADH Ϫ oxidation was observed (orange line in Fig. 8, A and B). All data suggest that 4-CP can bind to apoHadA to form a dead-end complex that does not lead to dechlorination (Step 10 in Fig. 5).
C.3. Binding Kinetics and Dissociation Constant of the HadA⅐4-CP Complex-The rate of 4-CP binding to HadA was investigated using double-mixing stopped-flow experiments. In the first mixing, an aerobic solution of HadA (75 M) was mixed with 4-CP (0.8, 1.6, 3.2, and 6.4 mM) and was incubated at various time periods before the solution was mixed against an anaerobic solution of FADH Ϫ (25 M) in the second mixing. Concentrations indicated are final concentrations. Longer incubation times and higher 4-CP concentrations allowed more complete formation of the HadA⅐4-CP complex, leading to the prevention of FADH Ϫ binding, and the resulting free FADH Ϫ is oxidized with t1 ⁄ 2 ϳ0.3 s. Therefore, the progress of 4-CP binding to HadA can be monitored by plotting ⌬A 450 nm at t ϭ 2 s versus incubation times of the first mixing. The plots were analyzed using an exponential equation to calculate the observed rate constant of HadA and 4-CP binding at each 4-CP concentration (supplemental Fig. S10A). The observed rate constant values were plotted against the 4-CP concentration to calculate k on , k off , and K d for the binding of 4-CP to HadA as 25 Ϯ 2 M Ϫ1 ⅐s Ϫ1 , 10 Ϯ 2 ϫ 10 Ϫ4 s Ϫ1 , and 40 Ϯ 10 M, respectively (supplemental Fig. S10B).
A dissociation constant (K d ) for the binding of HadA and 4-CP could also be calculated from the amount of free FADH Ϫ in the system because the amount of free FADH Ϫ is inversely proportional to the amount of 4-CP occupying HadA. Analysis of the K d value was carried out according to supplemental Equation S5. Upon increasing the 4-CP concentration, more free FADH Ϫ oxidation, which is oxidized faster than the HadAbound FADH Ϫ (Fig. 8B), was observed. Therefore, the increase of absorbance at 450 nm at 2 s (which is due to free FADH Ϫ oxidation) with increasing 4-CP was used for calculating the K d value for the binding of HadA and 4-CP as 29 Ϯ 4 M (inset in Fig. 8B). This value agrees well with the dissociation constant obtained from the kinetics experiments mentioned above. Fig. 5) The results in Sections A-C indicate that FADH Ϫ is the first substrate to bind to HadA. Although 4-CP can bind to HadA, the HadA⅐4-CP is a dead-end complex that does not lead to formation of C4a-hydroperoxy-FAD. To proceed with the productive path, 4-CP must bind the HadA⅐FADH Ϫ complex or the HadA⅐C4a-hydroperoxy-FAD intermediate. To identify whether 4-CP can bind to a binary complex of HadA⅐FADH Ϫ before it reacts with oxygen, two stopped-flow experiments were carried out to monitor rate constants of C4a-hydroperoxy-FAD formation and decay. In the first experiment, a solution of HadA⅐FADH Ϫ in the presence of 4-CP was mixed with aerobic buffer, and in the second experiment, a solution of HadA⅐FADH Ϫ binary complex alone was mixed with aerobic buffer plus 4-CP. Both experiments were monitored for absorbance changes at 380 and 450 nm to observe the formation of C4a-hydroperoxy-FAD. If the ternary complex can form, it likely reacts with oxygen differently from the reaction of the binary complex of HadA⅐FADH Ϫ with oxygen. Previous studies of C 2 oxygenase, which is also a two-component oxygenase, have shown that a ternary complex of C 2 ⅐FMNH Ϫ ⅐HPA reacts with oxygen ϳ50-fold slower than the reaction of a binary complex of C 2 ⅐FMNH Ϫ (25). The results (supplemental Fig. S11) show similar kinetic traces for both experiments, indicating that pre-incubation of 4-CP with the HadA⅐FADH Ϫ binary complex had no influence on the rate constant of C4a-hydroperoxy-FAD formation. These results suggest that the ternary complex of HadA⅐FADH Ϫ ⅐4-CP is not necessary and that 4-CP likely binds to the enzyme after C4a-hydroperoxy-FAD formation (Step 4 in Fig. 5).

E: Kinetics of HadA⅐FADH ؊ Binary Complex Reaction with O 2 in the Presence 4-CP
The results in Sections A-D have established the kinetic mechanism for the reaction of HadA from Steps 1-4 and 10 in Fig. 5. The knowledge obtained was used for setting up a mixing condition that allows the reaction of HadA to proceed in a productive manner by avoiding pre-mixing of HadA and 4-CP to prevent formation of the dead-end complex (Step 10 in  Fig. 9) showed three observable phases. The first two phases are the same as those observed for the reaction performed in the absence of 4-CP (Fig. 7, top inset), indicating that these two phases represent formation of C4a-hydroperoxy-FAD with a bimolecular rate constant of 7.0 Ϯ 0.2 ϫ 10 4 and 3.8 Ϯ 0.3 ϫ 10 3 M Ϫ1 ⅐s Ϫ1 , respectively, as explained previously in Section B.1. The last phase showed a decrease of absorbance at 380 nm and a large increase of absorbance at 450 nm that is independent of oxygen concentration. The kinetics of this phase are significantly different from those data of the reaction in the absence of substrate. The rate constant of the reaction in the presence of 4-CP was 6.9 Ϯ 0.2 ϫ 10 Ϫ2 s Ϫ1 , ϳ13-fold faster than that without substrate (5.4 Ϯ 0.2 ϫ 10 Ϫ3 s Ϫ1 ). Although the kinetics showed only a single observable phase, this third phase is the combination of the reactions of a coupling path in which dechlorination (hydroxylation and Cl Ϫ elimination) takes place (Steps 4 -7 in Fig. 5) and an uncoupling path in which oxidized FAD is produced by H 2 O 2 elimination (Steps 8 and 9 in Fig. 5). To obtain the individual rate constants associated with each reaction path, stopped-flow spectrofluorometric (Section E.2) and rapid-quench flow (Section E.3) experiments were carried out.

E.2. Detection of the Hydroxylation Rate Constant by Monitoring Formation of the C4a-Hydroxy-FAD Intermediate Using
Stopped-flow Spectrofluorometry-As the hydroxylation of 4-CP to form HQ requires the transfer of the hydroxyl group from C4a-hydroperoxy-FAD to the substrate to result in C4ahydroxy-FAD, monitoring the kinetics of the change between the two forms of the flavin intermediates would be a direct measurement of the rate constant associated with hydroxyl group transfer from C4a-hydroperoxy-FAD into 4-CP. Previous studies have shown that C4a-hydroperoxyflavin and C4ahydroxyflavin in some enzymes have different fluorescence  MARCH 24, 2017 • VOLUME 292 • NUMBER 12 characteristics (36,37,(41)(42)(43), permitting direct measurement of the hydroxylation step. Therefore, an experiment similar to Section E.1 was carried out, but the flavin fluorescence was monitored rather than absorbance. The fluorescence kinetics of the reactions in the presence of 4-CP was monitored by using excitation wavelengths at 380 nm (excitation 380) and 450 nm (excitation 450) to monitor emission of C4a-hydroxy-FAD and oxidized FAD, respectively. Both flavin species were detected by their fluorescence at emission wavelengths of Ͼ495 nm (emission Ͼ495). The results showed three observable phases (red lines in Fig. 9). The first phase (0.001-0.2 s) is a lag phase with a very small change of fluorescence signal at both excitation wavelengths. This phase occurs at the same time as the formation of HadA⅐C4a hydroperoxy-FAD monitored by absorbance. These data indicate that HadA⅐C4a-hydroperoxy-FAD is not highly fluorescent. The second phase (0.2-20 s) of the reaction showed a large increase in fluorescence at excitation 380, but only a small change in fluorescence at excitation 450. The observed rate constant of this phase was independent of oxygen concentration. This phase likely represents formation of C4a-hydroxy-FAD, which is known to be more fluorescent than C4a-hydroperoxy-flavin (36,37,(41)(42)(43). The last phase showed a large decrease in fluorescence at excitation 380 and a large increase in fluorescence at excitation 450, implying that this phase represents a dehydration reaction to eliminate H 2 O from C4a-hydroxy-FAD to form oxidized FAD (Step 7 in Fig. 5). A control reaction in the absence of 4-CP was also carried out, and the result did not show any formation or decay of excitation 380 species (supplemental Fig. S12); only the increase of fluorescence due to formation of oxidized FAD was observed. These results indicate that the fluorescence signal detected from excitation 380 and emission Ͼ495 in the presence of 4-CP is indeed derived from C4a-hydroxy-FAD. Therefore, the rate constant of the formation of fluorescence species with excitation 380 (second phase of the fluorescence trace, 0.36 Ϯ 0.02 s Ϫ1 , k obs, OH ) represents the observed rate constant of the incorporation of ϪOH from C4a-hydroperoxy-FAD to 4-CP (the hydroxylation step, Step 5 in Fig. 5).

Kinetic Mechanism of HadA Monooxygenase
As the product yield for 4-CP hydroxylation is not 100%, the observed hydroxylation rate (k obs, OH ) is the combination of the rate constants from the coupling (Step 5 in Fig. 5) and uncoupling paths after 4-CP binding (Steps 9 in Fig. 5) and can be described according to the kinetic model for a bifurcation reaction (Equation 1).
The true hydroxylation rate constant (k OH or k 5 ) that represents the hydroxylation step can be calculated based on the observed hydroxylation constant and product yield according to Equation 2.
The calculation yielded a hydroxylation rate constant for 4-CP (k OH or k 5 ) of 0.18 Ϯ 0.02 s Ϫ1 ).

E.3. Rate Constant of Product Formation Determined by Rapid-quenched Flow
Experiments-Although the previous stopped-flow spectrofluorometric experiments could measure a rate constant associated with the hydroxylation step, the measurements were extrapolated only from observing changes in flavins (C4a-hydroperoxy-FAD, C4a-hydroxy-FAD, and oxidized FAD) rather than from direct observation of the substrate 4-CP. Therefore, the rate constant associated with hydroxylation plus dechlorination of 4-CP was investigated using rapidquench flow (RQF) techniques. The RQF experiment was carried out in the same manner as the stopped-flow experiments, except that only oxygen at a final concentration of 0.13 mM was used. The reactions were carried out in 20 mM HEPES, pH 7.5, with 1 mM ascorbic acid at 25°C in the RQF machine that was placed inside an anaerobic glove box and were quenched at various time points from 0 to 100 s by adding an equal volume of 0.2 M HCl (final concentration). Samples were collected, prepared, and analyzed by HPLC/DAD to measure the amount of HQ formation at various times (see "Experimental Procedures"), and the results are shown by the green line with solid green circles (Fig. 9).
The observed rate constant (k obs, Pro ) of product formation was calculated using a single exponential equation as 0.24 Ϯ 0.01 s Ϫ1 . This value is in the same range as the observed rate constant of the hydroxylation step (Section E.2). As the reaction of HadA and 4-CP does not yield 100% product formation, the observed rate constant from RQF experiments (k obs, Pro ) also contains contributions from the uncoupling (H 2 O 2 elimination) path similar to the explanation in E.2. Therefore, a rate constant associated with the product formation steps (k Pro ) was calculated by Equations 3 and 4 as 0.12 Ϯ 0.01 s Ϫ1 . However, the rate constant of product formation is an observed value resulting from hydroxylation and chloride elimination reactions (Steps 5 and 6 in Fig. 5), which can be described according to the Equation 5.
As the hydroxylation rate constants (k OH ) can be calculated from measuring the kinetics of C4a-hydroxy-FAD formation (Section E.2), the rate constant associated with chloride elimination (k Eli or k 6 ) was calculated as 0.36 Ϯ 0.01 s Ϫ1 (Step 6 in Fig. 5).

F: Kinetics of the HadA⅐FADH ؊ Binary Complex Reaction with O 2 in the Presence 4-CP at Various pH Values
Reactions of HadA⅐FADH Ϫ and air-saturated 4-CP were also carried out at various pH values from 6.5 to 10.5 to explore whether the low product formation at pH Ͻ7.0 and Ͼ8.2 was due to the low yield of C4a-hydroperoxy-FAD in HadA (see "Experimental Procedures"). The kinetic traces obtained at absorbance of 380 nm showed that at different pH values the intermediates (C4a-hydroperoxy-FAD and C4a-hydroxy-FAD) form and decay at different rates (supplemental Fig. S13A). However, they all have the highest absorbance values (com-bined signals from the two intermediates) in a similar range (within 10% variation). This is in contrast to the yield of production formation at various pH values in which the yield variation is greater than 50% (Fig. 4). These data suggest that the variation of product formation at various pH values is not due to the impairment of HadA to form C4a-hydroperoxy-FAD but likely is due to the protonation status of substrate or active site residues that can affect the hydroxylation process.

Overall Reaction of HadA
The results from Sections A-E have established the kinetic mechanism of the HadA reaction (Fig. 5). HadA binds FADH Ϫ as the first substrate to form a binary complex of HadA⅐FADH Ϫ that can be equilibrated in two forms in which both forms can react with oxygen to form HadA⅐C4a-hydroperoxy-FAD. 4-CP is the third substrate to bind to the enzyme after the reactive intermediate C4a-hydroperoxy-FAD is formed. The fate of the ternary complex of HadA⅐C4a-hydroperoxy-FAD⅐4-CP bifurcates into coupling (hydroxylation plus Cl Ϫ elimination) and uncoupling paths of H 2 O 2 elimination.
In the absence of 4-CP, the HadA⅐FADH Ϫ binary complex reacts with oxygen to form HadA⅐C4a-hydroperoxy-FAD in which it eliminates H 2 O 2 slowly (0.005 s Ϫ1 at 25°C). We propose that in the absence of 4-CP in vivo, the enzyme predominantly accumulates in the HadA⅐C4a-hydroperoxy-FAD form. This way HadA can avoid formation of the HadA⅐4-CP deadend complex, and the enzyme is ready for dechlorination/hydroxylation activity. All kinetic and thermodynamic parameters of the individual steps are summarized in Table 2.

Discussion
Our results reported herein are the first elucidation of the kinetic mechanism of a two-component flavin monooxygenase that can catalyze oxidative dechlorination of various CPs. These enzymes are prevalent in nature especially in microbes that have potential uses for CP detoxification. Analysis of the products from single and multiple turnover reactions clearly demonstrate that HadA prefers to use 4-CP and 2-CP more than CPs with multiple substituents. Transient kinetics based on stopped-flow and rapid-quench flow experiments have shown the involvement of C4a-hydroperoxy-FAD and C4a-hydroxy-FAD in the oxidative dechlorination process of a two-compo-nent flavin-dependent dechlorinase. The data have also elucidated the order of substrate binding, the catalytic steps, and the rate constants involved and have identified non-productive and productive paths of the HadA reactions.
Data from the single turnover reactions (Table 1) have clearly elucidated the substrate specificity of HadA. The experiments directly measuring products resulting from the first catalytic cycle have shown that the hydroxylation efficiency of HadA with CPs are as follows: 4-CP ϳ 2-CP Ͼ 2,4-DCP Ͼ 2,4,5-TCP Ͼ 2,4,6-TCP, which is different from the previously published results of the HadA investigation using the lysate of Escherichia coli overexpressing the hadA gene, which showed that the multiply substituted CPs were metabolized by the lysate more efficiently than mono-substituted CPs (14,15). It can be envisaged that for the lysate experiments the overall bioconversion process can be influenced by other enzymes in the lysate that can further metabolize a first product. In single turnover experiments in which the purified HadA was used, the results can clearly state the substrate preference of HadA because other factors affecting the reaction efficiency, such as the ability of a reductase to generate FADH Ϫ and efficiency of the FADH Ϫ transfer between the two proteins versus the loss of FADH Ϫ during turnover, can be avoided. The results from single turnover reactions also agree well with the data from multiple turnover reactions in which purified HadA was reconstituted with C 1 reductase and G6PD. The data indicate that HadA prefers to use 4-CP Ͼ 2,4-DCP Ͼ 2,4,5-TCP Ͼ 2,4,6-TCP ( Fig. 2 and Table 1).

Summary of kinetic and thermodynamic parameters for the reaction of HadA⅐FADH ؊ with 4-CP
Step

Kinetic and thermodynamic parameters
HadA⅐FADH ؊ with 4-CP reactions, which was speculated to be a C4a-flavin adduct (9). With our current findings, it is very likely that the species previously detected in the multiple turnovers of TftD was either C4a-hydroperoxy-FAD or C4a-hydroxy-FAD. Studies of a single component pentachlorophenol hydroxylase (PcpB) from S. chlorophenolicum (18,19) also showed the existence of C4aflavin adducts.
The transient kinetic studies in this report elucidated the order of substrate binding and identified the rate-limiting step in the overall reaction of HadA. For the productive path, HadA binds to FADH Ϫ as the first substrate to form a binary complex of HadA⅐FADH Ϫ that reacts with oxygen to result in C4a-hydroperoxy-FAD. The detection of highly fluorescent flavin species only after 4-CP binds to the C4a-hydroperoxy-FAD demonstrates the involvement of C4a-hydroxy-FAD in the HadA dechlorination. C4a-hydroxy-FAD has been shown to be highly fluorescent in flavin-dependent enzymes (36,37,(41)(42)(43). As this step occurs before Cl Ϫ elimination to form hydroquinone, the data suggest that the transfer of a terminal ϪOH from C4ahydroperoxy-FAD to 4-CP results in C4a-hydroxy-FAD that occurs before Cl Ϫ elimination. Based on the values of all the rate constants involved in the reaction of HadA ( Fig. 5 and Table 2), the hydroxylation step is a major rate-limiting step of the HadA reaction. This result agrees with computational calculation models of the reaction of TftD in which the hydroxylation step was proposed to be the rate-determining step for dechlorination (50). The involvement of C4a-flavin adduct species in the HadA reaction makes its mechanism unique from the mechanisms of reductive dehalogenation performed by tyrosine deiodinase from human (51,52) and glutathione-dependent dehalogenase from S. chlorophenolicum (53).
HadA also has various mechanistic properties distinct from other enzymes that make up the class D monooxygenases. 4-CP can bind to HadA to form a dead-end complex of HadA⅐4-CP in which its binding obstructs the productive path of FADH Ϫ binding (Fig. 8). This result is different from those of the p-hydroxyphenylacetate hydroxylases (members of class D monooxygenase) from A. baumannii and Pseudomonas aeruginosa in which the substrate HPA cannot bind to the apoenzymes (25,26,54). However, formation of an enzyme⅐substrate dead-end complex was also observed in bacterial luciferase (a member of the class C monooxygenases) in which an aldehyde substrate can bind to the enzyme and prevent the binding of FMNH Ϫ that is required for proceeding to a productive path (39,40). However, as the rate constants of FADH Ϫ binding and formation of C4a-hydroperoxy-FAD ( Fig. 5 and Table 2) are much faster than the rate constant of H 2 O 2 elimination, it is possible that under physiological conditions and in the absence of CPs the enzyme mainly exists in the C4a-hydroperoxy-FAD form. This mode would be an efficient means for avoiding formation of a dead-end complex of HadA⅐4-CP such that the enzyme can be maintained in a state ready for the dechlorination process. When CPs are present, the ligand can bind to the HadA⅐C4a-hydroperoxy-FAD, and the reaction can promptly proceed to a productive dechlorination path. C4a-hydroperoxyflavin is also thought to be a major species that accumulates during the reactions of p-hydroxyphenylacetate hydroxylase from A. baumannii and Pseudomonas aeruginosa under physiological conditions (25,26,54).
HadA can bind to FADH Ϫ and exist in two populations consisting of two different forms of HadA⅐FADH Ϫ that can react with oxygen to form C4a-hydroperoxy-FAD with different bimolecular rate constants. A longer incubation time resulted in more of the species that reacts slower with oxygen. This property does not exist in other enzymes in the class D monooxygenases such as p-hydroxyphenylacetate hydroxylases (members of class D monooxygenase) from A. baumannii and P. aeruginosa (25,26,54). However, similar results have been observed in a few two-component monooxygenases in other classes. Similar properties were reported for a flavin-dependent halogenase or RebH (class F monooxygenase). Preincubation of RebH and FADH Ϫ resulted in a less active form of RebH⅐FADH Ϫ that has a lower product yield (42). For alkanesulfonate monooxygenase or SsuD (class C monooxygenase), when SsuD⅐FADH Ϫ was premixed prior to reacting with oxygenated substrate, the amount of product formation is 4-fold lower than that produced from the mixing of free FMNH Ϫ and a solution of SsuD containing oxygenated substrate (55). However, unlike in the HadA reaction, stopped-flow experiments of RebH and SsuD reactions in which complexes of enzyme and reduced flavins were preformed did not reveal formation of the C4a-hydroperoxyflavin intermediate (42,55). The results of the HadA reaction are similar to the stopped-flow results of luciferase from Photobacterium leiognathi TH1 (PlLuxAB) in which two phases of C4a-hydroperoxy-FMN formation were detected when mixing a solution of preformed PlLuxAB/FMNH Ϫ with O 2 , whereas only one phase was detected in the non-premixed reaction (56). It is not clear why these enzymes possess this conformational change property. The mixture of enzyme populations that react differently with oxygen may be involved with the control of oxygen diffusion into the active site or the ability of reduced flavin to form the C4a-hydroperoxyflavin intermediate. It is known that diffusion of oxygen into a flavin-dependent monooxygenase can influence the rate of C4a-hydroperoxyflavin formation (57). Recently, it has been shown that formation of C4a-hydroperoxyflavin in a flavin-dependent enzyme occurs via a mechanism of proton-coupled electron transfer (58,59). In the C 2 reaction, in which the formation of C4a-hydroperoxyflavin is almost barrier-less with a bimolecular rate constant at 4°C of 1.1 ϫ 10 6 M Ϫ1 s Ϫ1 , a His residue that is located near the flavin C4a-position is necessary to facilitate the rapid formation of the intermediate (25,28). When we examined the X-ray structure of apoTcpA (Protein Data Bank code 4G5E) (21) that is similar to HadA and aligned it with the C 2 structure (Protein Data Bank code 2JBT) (60), His-290 and Arg-101 could be identified as potential proton donors at the active site. It is possible that the conformational change of HadA to a slower reacting species may be due to the change in the positions of these residues with respect to the flavin C4aposition, resulting in a slower rate of oxygen reaction. However, the conformational change in this case does not affect the overall turnover of HadA because the rate-limiting step of HadA is not the oxygen reaction step (Fig. 5 and Table 2).
Overall, the results in this report provide a solid basic mechanistic understanding of the flavin-dependent oxidative dechlorinases that will be useful for future applications. The data can also be used for comparison with future site-directed mutagenesis studies. As HadA homologs are found widespread in nature, an understanding of the HadA reaction ( Fig. 5 and Table 2) also contributes to the general understanding of reactions of its homologs identified by SSN analysis (supplemental Fig. S1) (9,10,13,(33)(34)(35). These enzymes were identified in microbes that showed potential in bioremediation applications.
In conclusion, this study is the first report of the kinetic mechanism of a flavin-dependent monooxygenase (dechlorinase) HadA from R. pickettii DTP0062. Combined approaches of transient kinetics have elucidated the order of substrate binding and rate constants associated with the individual steps. The data have also revealed the involvement of C4a-hydroperoxy-FAD and C4a-hydroxy-FAD in the oxidative dechlorination process. As the usage of CPs is still widespread worldwide, understanding how nature uses flavin and oxygen to battle these toxic chemicals provides knowledge that may lead to future innovations for detoxifying chemicals hazardous to human health.

Experimental Procedures
Construction of the hadA Gene-A synthetic gene of hadA based on the sequence shown in GenBank TM accession no. BAA13105 (15) was synthesized by GenScript. The gene contains 1,554 nucleotides that encode for 517 amino acids. The hadA synthetic gene was digested with NdeI and BamHI restriction enzymes (Novagen) at the 5Ј and 3Ј ends, respectively, before it was ligated with the E. coli expression vector, pET-11a, to generate the hadA-pET-11a vector. The sequence of the hadA-pET-11a vector was analyzed (Macrogen) to verify its identity.
Protein Expression-The hadA-pET-11a plasmid carrying a recombinant hadA gene was transformed into E. coli BL21 (DE3) cells. The resulting transformant was grown on an LB agar plate containing 50 g/ml ampicillin at 37°C overnight. Auto-induction medium (61) was used to grow cells and induce protein expression. Briefly, a single colony of grown cells was inoculated into 100 ml of ZYM-5052 starter medium containing 50 g/ml ampicillin. The starter was incubated with shaking at 220 rpm, 37°C overnight. Starting culture cells (1% of v/v) were then inoculated into six flasks of 650 ml of ZYP-5052 broth containing 50 g/ml ampicillin. The culture was incubated with shaking at 220 rpm, 37°C, until the A 600 reached 1.0. The temperature was then switched to 25°C, and the cultures were continuously shaken for an additional 16 h. Cell paste was harvested by centrifugation and stored at Ϫ80°C until needed.
Protein Purification-The cell paste was suspended in 50 mM sodium phosphate (NaP i ), pH 7.0, containing 5 mM EDTA, 100 M phenylmethanesulfonyl fluoride (PMSF), and 1 mM dithiothreitol (DTT). Cells were disrupted by ultrasonication before being centrifuged at 27,000 ϫ g for 1 h to discard cell debris. Ammonium sulfate ((NH 4 ) 2 SO 4 ) precipitation (0 -20%) was used to remove the remaining cell debris and other contaminant proteins. The protein fraction containing the HadA enzyme was precipitated in 20 -40% (NH 4 ) 2 SO 4 . The resulting pellet was resuspended in ϳ50 ml of 30 mM NaP i , pH 6.5, containing 50 mM NaCl (buffer A) and was dialyzed against 4 liters of buffer A. Denatured proteins were removed by centrifugation. The remaining clear solution of proteins was loaded onto a DEAE-Sepharose TM fast-flow (GE healthcare) anion exchange column (volume ϳ200 ml), which was pre-equilibrated with buffer A. The column was washed with 10 column volumes of buffer A to get rid of unbound proteins. Elution of protein was performed using a linear gradient of 50 -400 mM NaCl in 30 mM NaP i , pH 6.5. The fractions were pooled and concentrated, and the protein was exchanged into 20 mM HEPES, pH 7.5, by passing it through a Sephadex TM G-25 column (GE Healthcare). Purified protein was stored at Ϫ80°C before use. All purification processes were performed at 4°C. The total protein amount was determined using a Bradford assay, and the purity of HadA was assessed by SDS-PAGE (15%) with Coomassie Brilliant Blue staining and compared with standard protein markers (Enzmart Biotech, Thailand). Concentrations of purified HadA were determined based on absorbance at 280 nm using the extinction coefficient value of 74.9 mM Ϫ1 cm Ϫ1 calculated from ProtParam tool of ExPASy.
Identification of Products from Single Turnover Reactions-A solution of HadA (75 M) and FAD (25 M) in 20 mM HEPES, pH 7.5, was placed inside an anaerobic glove box to remove oxygen. The anaerobic solution of HadA/FAD was reduced by adding ϳ5 mg/ml dithionite stock solution into an equivalent amount (or slightly higher due to trace amount of oxygen) of FAD to produce HadA⅐FADH Ϫ . Absorption spectra of the enzyme solution were monitored to ensure the complete reduction of the complex. The resulting mixture was mixed with an equal volume of air-saturated substrate (0.8 mM substrate and 0.13 mM oxygen) in tightly closed microcentrifuge tubes. Note that all reactions contained excess ascorbic acid (1 mM) to convert the resulting quinone product after Cl Ϫ elimination into hydroquinone compounds that are more stable and easier to quantitate by HPLC analysis. All concentrations are indicated as final concentrations after mixing. The effect of ascorbic acid on the catalytic cycle of HadA was tested, and it was found not to interfere with the HadA reaction. After mixing, the final concentrations of the reagents were half of the original concentrations. Under these conditions, FADH Ϫ was the limiting reagent, and the reaction only proceeds for a single turnover. Reactions were left 15 min at room temperature before being quenched by adding an equal volume of acetonitrile. The resulting white precipitate of denatured protein was removed by centrifugation at 13,800 ϫ g, 4°C, for 30 min. Supernatants were diluted 3-fold in 20 mM HEPES, pH 7.5, before being filtered through an Amicon Ultra-0.5 centrifugal filter device with a molecular mass cutoff of 10 kDa (Millipore). Filtrates were collected and analyzed by HPLC/DAD or HPLC/DAD/mass spectroscopic detector (MSD) (Agilent) (see below). Commercially available substrates and products were purchased from Sigma, Merck, and Tokyo Chemical Industry Co., Ltd. (TCI). Standard curves of each authentic compound were constructed to calculate the amount of product formed or substrate consumed. Coupling percentage was calculated by comparing the product yield with the limiting reagent, FADH Ϫ .
Bioconversion by Multiple Turnover Reactions-Multiple turnover reactions of HadA were carried out by coupling the reaction with those of G6PD (Sigma) and flavin reductase (C 1 from A. baumannii). Reactions (5 ml) contained substrate (100 M), Glc-6-P (1 mM), G6PD (0.5 unit/ml), C 1 (1 M), NAD ϩ (5 M), FAD (10 M), HadA (10 M), and ascorbic acid (1 mM) in 20 mM HEPES, pH 7.5. Reactions were stopped at various time points to separate products from proteins similar to the protocol described above. The amount of product formed was analyzed using HPLC/DAD or HPLC/DAD/MSD as described below.
Product Analysis-Products from single and multiple turnover reactions were analyzed using an HPLC (Agilent Technologies 1100 or 1260 Infinity series) equipped with a UV-visible DAD and/or electrospray ionization mass spectroscopic detector (ESI-MSD) (Agilent Technologies 6120 Quadrupole LC/MS). The stationary phase was a Nova-Pak (Waters) C18 reverse phase column with a 4-m particle size and a 3.9 ϫ 150-mm column size. Mobile phases were either a gradient of H 2 O/methanol containing 0.1% formic acid or a gradient of H 2 O/acetonitrile containing 0.1% formic acid. All HPLC grade solvents were purchased from Burdick & Jackson Scientific Co., Ltd. The HPLC was run at ambient temperature with a flow rate of 0.5 ml/min, and a sample injection volume of 20 l. The MSD was set in negative mode with a scan range of m/z from 10 to 500 or selected ion monitoring mode that was set at the specific m/z of the substrate and expected product. Products were analyzed by comparing their HPLC retention time and mass with those of standard compounds.
pH Effects on HadA Reactions-To investigate the effects of pH on HadA, single turnover reactions of HadA using 4-CP as a substrate were carried out at various pH values. A pH jump method in which an enzyme solution in a lower ionic strength buffer at pH 7.5 was mixed with buffers with higher ionic strength at various pH values was employed. An anaerobic mixture of HadA and FADH Ϫ (75 and 25 M, respectively) was prepared in 20 mM HEPES, pH 7.5, as described above in an anaerobic glove box. To the mixture of HadA⅐FADH Ϫ was added an equal volume of air-saturated 4-CP (0.8 and 0.13 mM oxygen) in a three buffer system consisting of N-(2-acetamido)-2-aminoethanesulfonic acid/Tris/ethanolamine (ionic strength ϭ 0.1 M) at various pH values with ascorbic acid (1 mM). This buffer system covered a pH range of 5.5-10.5 (62). Samples were processed and analyzed for product formation by HPLC/ DAD as described above.
Transient Kinetics of the HadA Reaction-Single mixing or double mixing stopped-flow spectrophotometer model SF-61SX and rapid-quench flow model RQF-63 from TgK Scientific (Bradford-on-Avon, UK) were used to monitor the kinetics of the HadA reaction. The flow path of the instrument was made anaerobic by flushing it with a solution of ϳ5 mg/ml sodium dithionite and leaving dithionite solution in the instrument overnight before use. The flow unit was rinsed three times with anaerobic buffer that was bubbled with nitrogen gas to remove excess sodium dithionite before performing the experiment. To generate a solution of the HadA⅐FADH Ϫ complex, an anaerobic solution of FAD (25 M) and HadA (75 M) in 20 mM HEPES, pH 7.5, was stoichiometrically reduced by adding an anaerobic solution of sodium dithionite (ϳ5 mg/ml) in an anaerobic glove box. The reduced binary complex was placed inside an anaerobic glass tonometer and was later loaded onto the stopped-flow spectrophotometer. A solution of the HadA⅐ FADH Ϫ complex was mixed with 20 mM HEPES, pH 7.5, buffer containing 4-CP and various concentrations of oxygen at 25°C in the stopped-flow instrument. Oxygenated buffers were prepared by bubbling the buffer with various mixtures of oxygen and nitrogen gases (final concentration of 0.13, 0.31, and 0.61 mM oxygen) (25). Kinetic traces were obtained by monitoring the absorbance at wavelengths of 380 and 450 nm. Observed rate constants of the reactions were analyzed using ProgramA software (developed by C. J. Chiu, R. Chang, J. Diverno, and D. P. Ballou, University of Michigan, Ann Arbor). The rate constants of the oxygen-dependent phases were plotted using KaleidaGraph software (Synergy Software) and analyzed as a bimolecular reaction with oxygen.
Analysis of Sequence Similarity Network (SSN)-A sequence similarity network (SSN) for HadA was carried out using the option A of Enzyme Function Initiative (EFI)-Enzyme Similarity Tool and sequences from UniProt database that had identity to HadA of 34% or higher (63,64). SSNs of HadA and its homologs were generated using an EFI alignment score of 70 (equivalent to grouping sequences of Ͼ27.74% in the same cluster). Networks of these sequences were visualized using Cytoscape version 3.2.1 (65).
Author Contributions-P. P. performed experiments, analyzed the results, and wrote the manuscript. K. T. performed dechlorinase sequence analysis. J. S. gave suggestions on data analysis. P. C. conceived the study, analyzed the results, and wrote the manuscript.