Dissecting the low catalytic capability of flavin-dependent halogenases
Although flavin-dependent halogenases (FDHs) are attractive biocatalysts, their practical applications are limited because of their low catalytic efficiency. Here, we investigated the reaction mechanisms and structures of tryptophan 6-halogenase (Thal) from Streptomyces albogriseolus using stopped-flow, rapid-quench flow, quantum/mechanics molecular mechanics calculations, crystallography, and detection of intermediate (hypohalous acid [HOX]) liberation. We found that the key flavin intermediate, C4a-hydroperoxyflavin (C4aOOH-FAD), formed by Thal and other FDHs (tryptophan 7-halogenase [PrnA] and tryptophan 5-halogenase [PyrH]), can react with I−, Br−, and Cl− but not F− to form C4a-hydroxyflavin and HOX. Our experiments revealed that I− reacts with C4aOOH-FAD the fastest with the lowest energy barrier and have shown for the first time that a significant amount of the HOX formed leaks out as free HOX. This leakage is probably a major cause of low product coupling ratios in all FDHs. Site-saturation mutagenesis of Lys79 showed that changing Lys79 to any other amino acid resulted in an inactive enzyme. However, the levels of liberated HOX of these variants are all similar, implying that Lys79 probably does not form a chloramine or bromamine intermediate as previously proposed. Computational calculations revealed that Lys79 has an abnormally lower pKa compared with other Lys residues, implying that the catalytic Lys may act as a proton donor in catalysis. Analysis of new X-ray structures of Thal also explains why premixing of FDHs with reduced flavin adenine dinucleotide generally results in abolishment of C4aOOH-FAD formation. These findings reveal the hidden factors restricting FDHs capability which should be useful for future development of FDHs applications.
<!>Halogenation activity of Thal<!><!>Determination of the hydrogen peroxide elimination ratio from the uncoupling (nonproductive) path<!>Kinetics of hypohalous acid (HOCl, HOBr, and hypoiodous acid) formation catalyzed by Thal, PrnA, and PyrH<!><!>Kinetics of hypohalous acid (HOCl, HOBr, and hypoiodous acid) formation catalyzed by Thal, PrnA, and PyrH<!><!>Investigation of the functional role of Lys79 using site-saturation mutagenesis<!><!>Computational calculations to investigate the pKa of Lys79<!><!>Kinetics of Thal:FADH-∗ inactive complex formation<!><!>Crystal structures of Thal: Implications for limitation of noniodination activity<!><!>Discussion<!>Single turnover reaction of Thal<!>Transient kinetics experiments investigated by stopped-flow spectrometry<!>QM/MM calculations<!>Investigation of HOX leakage from Thal<!>Characterization of brominated D-luciferin products using HPLC-QTOF<!>Site-saturation mutagenesis of Thal at Lys79<!>Crystallization of Thal<!>Data availability<!>Conflict of interest<!>Supplementary data
<p>Edited by F. Peter Guengerich</p><p>Halogenation is extremely important for the agricultural, pharmaceutical, and chemical industries (1). It is especially of importance for drug synthesis because one-third of currently produced drugs, particularly new drugs currently in clinical trials are halogenated compounds (2, 3, 4, 5, 6, 7, 8, 9). However, by traditional chemical synthesis, halogenations are generally not environmentally friendly reactions because most of these reactions need to be carried out under harsh conditions and often involve the use of toxic reagents, catalysts, and solvents. Chemical halogenation reactions also lack regioselectivity, resulting in many by-products (10, 11, 12, 13). The use of enzymes, i.e., halogenases, to catalyze halogenation is an attractive green chemistry method because the reactions can be done under mild conditions. Halogenases discovered to date can be classified into four groups: haloperoxidases, fluorinases, α-ketoglutarate-dependent halogenases, and flavin-dependent halogenases (FDHs) (14).</p><p>FDHs belong to the class of two-component flavin-dependent monooxygenases which require a reductase component to generate reduced FAD (FADH−) for a halogenase (monooxygenase) component to catalyze halogenation (14, 15, 16, 17, 18, 19). Several X-ray structures of FDHs have been solved (18) which showed that all FDHs have two separated binding sites. While the FAD binding site is well conserved among different FDHs, a significant difference lies in the substrate binding site which allows FDHs to utilize a broad range of substrates and catalyze halogenation at various positions (14). Therefore, most of the previous studies have been focused on enzyme engineering to alter substrate scope and regio-specificity (14, 20, 21, 22, 23, 24, 25, 26), whereas less attention has been paid to mechanistic investigation.</p><p>To date, mechanistic investigation of FDHs has only been carried out with the reactions of tryptophan 7-halogenases (RebH (15, 27, 28) and PrnA (29)) with only one investigation of transient kinetics of RebH reported (27, 28). The results indicate that RebH can form a C4a-hydroperoxyflavin intermediate (C4aOOH-FAD) to react with Br− and Cl− (the reactions with I− and F− were not explored) to generate a reactive halogenation species (28), hypohalous acid (HOX) (hypobromous acid [HOBr] and hypochlorous acid [HOCl]). HOX is thought to be transferred to a second active site to react with tryptophan (27, 28, 29). Because site-directed mutagenesis of the conserved Lys to Ala (15, 30, 31) and to Arg (31) resulted in inactive tryptophan halogenases, two models to explain the mechanism of halogenation have been proposed. In the first model, HOX is thought to react with Lys to form the chloramine intermediate which acts as a halogenating species (27), whereas in the second model, HOX is thought to be a direct halogenation reagent with the conserved Lys facilitating and controlling the binding of HOX in the proper geometry necessary for halogenation (15, 32). Interestingly, none of the studies since then have investigated the function of the conserved Lys. None of the previous studies has also investigated the mechanism of HOX formation and determined the factors controlling HOX production. Most FDHs are inefficient because their turnovers are quite slow (kcat = 0.5–3 min−1), and they can be easily inactivated (26). Despite gram scale synthesis of halogenated tryptophan product was reported for the reaction of RebH, full conversion was obtained only after 8 days (33), which is too slow to be accepted for a real industrial process. Therefore, kinetic and structural investigation to understand mechanistic features controlling their halogenation are necessary for future improvement of FDHs for industrial applications.</p><p>Tryptophan 6-halogenase or Thal from Streptomyces albogriseolus (34, 35) represents an interesting FDH system for investigating the mechanisms of HOX formation and factors controlling halogenation. As Thal catalyzes tryptophan halogenation at the 6-position, not 7-position as in RebH or PrnA, the results would also contribute to the understanding of a mechanism related to halogenation of tryptophan at the 6-position. Recently, crystal structures of Thal were solved and the results showed that FAD and tryptophan substrates could not be co-crystallized (34). This feature is different from other tryptophan halogenases such as RebH (36), tryptophan 5-halogenase (PyrH) (37) and PrnA (15) which could be co-crystallized with both FAD and tryptophan. This suggests that protein dynamics between the flavin and substrate binding sites of Thal are strongly linked. Therefore, we set out to investigate the mechanism of HOX formation and its role in regulating the Thal reaction and used Thal as a model for exploring the factors contributing to the catalytic inefficiency of FDHs in general.</p><p>In this report, the reaction of Thal was investigated with regard to four main mechanistic issues to unravel why the reactions of FDHs are not very efficient. (A) The mechanism of HOX formation was investigated using transient kinetics (stopped-flow absorbance and fluorescence and rapid-quench flow techniques) in combination with quantum/mechanics molecular mechanics (QM/MM) computational calculations. (B) The cause of the low product coupling ratio was investigated by measuring the leakage of the HOX intermediate using an organic substrate (D-luciferin) as a halogenating target to detect free HOX. (C) The functional role of the conserved Lys was investigated using site-saturation mutagenesis studies along with detection of free HOX formation and theoretical pKa calculations. Unlike the mechanism previously proposed, our results do not support formation of chloramine or bromamine intermediates but rather suggest that the conserved Lys likely protonates the HOX in the halogenation reaction. (D) Results from transient kinetics, newly crystallized structures of Thal:FADH− and Thal:FAD:AMP and molecular dynamics (MD) simulations could be used to explain the cause for formation of the dead-end inactive Thal:FADH-∗ complex.</p><!><p>To explore whether Thal represents a halogenase which can catalyze a reasonable halogenation reaction, multiple turnover reactions of Thal were carried out and compared with the reactions of PrnA and PyrH. The products from the Thal reaction were characterized by NMR spectroscopy (Fig. S2C). The results indicated that among the three halogenases investigated, Thal was the fastest halogenase (Fig. S2A). These data also clearly imply that most of the halogenases are not very efficient enzymes which agrees with previous data reporting that FDHs suffer from having low kcat values (26). Therefore, understanding the cause of catalytic inefficiency in Thal as a representative of FDHs would contribute to identification of the bottleneck in the reaction of FDHs in general.</p><!><p>Reaction of Thal with tryptophan. Thal catalyzes halogenation of tryptophan to yield halogenated tryptophan as the product. FAD, oxidized flavin adenine dinucleotide; FADH−, reduced flavin adenine dinucleotide.</p><p>Kinetic traces of the Thal reaction in the presence of halide ions. An air-saturated solution of Thal premixed with NaCl (blue lines) or NaBr (green lines) was mixed with an anaerobic solution of FADH−. Black lines are kinetic traces without addition of halide ions. Flavin intermediates were monitored by absorption changes at 380 nm (solid lines) and 450 nm (dashed lines). C4aOOH-FAD, C4a-hydroperoxyflavin; FADH−, reduced flavin adenine dinucleotide.</p><!><p>The kinetics of C4aOOH-FAD formation in the presence of 10 mM Br− or Cl− (the condition in which HOBr or HOCl was formed maximally, as explained in detail in the next section) was the same as that in the absence of Br− or Cl− (Fig. 2), but the kinetics of the reaction after C4aOOH-FAD formation was different. The reaction in the presence of halide ions showed three phases. The first phase (0.02–0.15 s) showed an increase in absorbance at 380 nm. The second (0.15–2.5 s) and third phases (2.5–100 s) showed a decrease in absorbance at 380 nm and an increase in absorbance at 450 nm. The kinetics of the second phase of both reactions (kobs of 1.6 s−1) were similar to that in the absence of halide, indicating that this phase was the uncoupling path to elimination of H2O2 without forming HOBr or HOCl. As the amplitude of this phase (absorbance 450 nm) in the presence of Br− and Cl− was about 20% and 55%, respectively, the results indicate that the uncoupling path is more prominent in the presence of Cl− compared with Br− and that C4aOOH-FAD reacts more efficiently with Br− than Cl−.</p><p>The kinetics of the third phase was clearly different from the reaction in the absence of the halide ions. This phase was later identified as the step of C4a-hydroxyflavin (C4aOH-FAD) decay after HOX formation (see the following results in the next section). These data indicate that even under saturating amounts of halide, the reaction of Thal still proceeds through an uncoupling path (no formation of HOX) with a frequency of about 20% and 55% in the case of Br− and Cl−, respectively. The data also indicate that 80% of the C4aOOH-FAD intermediate can react with Br− to form HOBr. However, our results in the first section (halogenation activity of Thal) indicate that only 30% of the tryptophans were brominated, implying that not all of the HOBr formed reacts with tryptophan. We therefore investigated the mechanisms of HOX formation by FDHs and examined their control of HOX reactivity.</p><!><p>Because the reaction of halide ions with C4aOOH-FAD is crucial for the catalysis of FDHs and this reaction is not well understood, we investigated the reactivity of halide ions with Thal and other tryptophan halogenases (PrnA and PyrH) by monitoring the kinetics of its formation.</p><!><p>Kinetic mechanisms of the Thal reaction with tryptophan. FAD refers to oxidized FAD. FADH− refers to reduced FAD. Thal:FADH-∗ refers to inactive complex which cannot form C4aOOH-FAD. TRP refers to tryptophan. C4aOOH-FAD refers to C4a-hydroperoxyflavin. C4aOH-FAD refers to C4a-hydroxyflavin. X− represents the halide ions. HOX represents hypohalous acid. In this model, TRP can bind to Thal together with a flavin.</p><p>Kinetic traces of the Thal reaction with various halide ions. Formation of C4aOH-FAD was monitored at Em > 495 nm using excitation wavelengths of 380 nm (solid lines) and 450 nm (solid lines with circles) in the presence of various halide ions (Cl− (blue), Br− (orange), I− (green), and F− (red)). Gray lines are reactions without addition of halide ion. Note that the difference in fluorescence intensities was because of the difference of photomultiplier tube voltage used for the various halide reactions. Inset shows a plot of the observed rate constants of C4aOH-FAD formation resulting from the reaction of C4aOOH-FAD with halide ions (Cl− (circle), Br− (square), and I− (triangular)) at various concentrations. C4aOOH-FAD, C4a-hydroperoxyflavin; C4aOH-FAD, C4a-hydroxyflavin; Em, emission wavelength.</p><!><p>To explore whether other halide ions (F−, I−, and Br−) can also react with C4aOOH-FAD in the Thal reaction, reactions of Thal with various halide ions were carried out. Results in Figure 4 indicate that C4aOOH-FAD in Thal can react with I−, Br−, and Cl− to form a highly fluorescent C4aOH-FAD intermediate. The exception was found for F− in which a highly fluorescent flavin could not be detected. These results clearly indicate that Thal can catalyze the reaction of C4aOOH-FAD with I−, Br−, and Cl− to form C4aOH-FAD and hypoiodous acid (HOI), HOBr, and HOCl, respectively, and it is the first evidence to show formation of HOI in the reaction of an FDH. Previously, the reaction of RebH was only investigated in regard to its ability to form HOBr and HOCl.</p><p>To explore whether the reactivity of Thal to form HOI is only specific to this enzyme, we also examined the reactivity of other halogenases such as PrnA and PyrH with I−. Similar experiments to those in Figure 4 were carried out. The results (Fig. S7) clearly showed that PrnA and PyrH could also react with I− to form HOI and C4aOH-FAD. However, product analysis indicates that Thal could not catalyze iodination of tryptophan (Fig. S8), indicating that the factor limiting iodination of tryptophan is not on the flavin reaction site (described in more detail using crystal structures of Thal in the last section).</p><!><p>Reactivity of halide ions with C4aOOH-FAD to form C4aOH-FAD</p><p>C4aOOH-FAD, C4a-hydroperoxyflavin; C4aOH-FAD, C4a-hydroxyflavin.</p><p>Potential energies of various states for HOX (X = I, Br, or Cl) formation catalyzed by Thal</p><p>HOBr, hypobromous acid; HOCl, hypochlorous acid; HOI, hypoiodous acid; HOX, hypohalous acid.</p><p>Detection of HOBr leakage from Thal (wild-type). HPLC-MS/MS chromatograms of the reaction without Thal (black) and the reaction in the presence of Thal (red). Multiple turnover reactions were performed to generate HOBr. A reaction solution after removal of enzymes was incubated with D-luciferin. Reactions were analyzed using HPLC-MS/MS. Mass spectra of bromo-D-luciferin and dibromo-D-luciferin were analyzed using HPLC-QTOF (inset). HOBr, hypobromous acid; QTOF, Quadrupole Time-of-Flight.</p><!><p>As the conserved catalytic Lys (Lys79Thal) is present in all FDHs reported and only two variants of the catalytic Lys (substitution of Lys to Ala (15, 27) or Arg (31)) were reported in the past, we subsequently mutated Lys79Thal to each of the other 19 amino acids using a site-saturation mutagenesis approach (Experimental Procedures) to investigate whether halogenation can also be achieved by other types of amino acids in this position. As the results in the previous section indicate that HOX can leak out from the wild-type enzyme, we also wanted to explore whether changing Lys79 to other residues could change the level of free HOX leakage. If Lys79 can form the chloramine or bromamine as previously proposed, exchanging this residue with other amino acids would generate significantly higher amounts of free HOX.</p><p>Plasmid libraries for site-saturation mutagenesis were constructed and evaluated to confirm that they generated all 20 amino acids in the Lys79 position, and a total 168 clones of Lys79 variants were tested for their ability to halogenate tryptophan (Experimental Procedures). The activity of the variants was analyzed by detecting the bromotryptophan product using HPLC-MS/MS. The protein production level of each variant was also measured using a GFP split system (47, 48). The results showed that all of the variants were completely inactive for tryptophan bromination (data not shown). These data confirm that Lys79 is crucial for the Thal reaction and cannot be replaced by another amino acid. One of the inactive variants, Lys79Thr, which expressed well was selected for further mechanistic investigation to understand the catalytic role of Lys79 (Fig. S18).</p><!><p>Kinetics of Lys79Thr reacting with oxygen. An air-saturated solution of Lys79Thr (blue lines) or wild-type (black lines) (30 μM) was mixed with a solution of FADH− (15 μM) using a single mixing stopped-flow spectrometer. Kinetic traces were monitored at 380 nm (solid lines) and 450 nm (dashed lines). FADH−, reduced flavin adenine dinucleotide.</p><p>ComparingD-luciferin bromination acitivity of Lys79Thr with wild-type Thal. Peak areas of brominated D-luciferin product of Lys79Thr (orange), wild-type (green), and the reaction without Thal (blue).</p><p>Proposed reaction mechanisms of tryptophan halogenation. The reaction was proposed to occur via electrophilic aromatic substitution. A, in tryptophan 7-halogenase (RebH), the Walsh group proposed that tryptophan can be halogenated by a covalent intermediate (Lys-NH-X) (27). B, alternatively, HOX can directly interact with a fully protonated Lys via hydrogen bonding (29, 32) to halogenate tryptophan. X represents I, Br, and Cl. Our findings support mechanism B and suggest the role of Lys in protonation of HOX. HOX, hypohalous acid.</p><!><p>As our results in the previous section imply that Lys79 does not form a chloramine or bromamine intermediate, this suggests that their role may be to form hydrogen bonding with HOX as originally proposed by van Pée et al. (29). However, the participation of Lys79 in the reaction of Thal is crucial because results from site-saturation mutagenesis experiments showed that Lys79 cannot be replaced by any other amino acids. To investigate why all FDHs contain a strictly conserved Lys as a catalytic residue, propka 3.0 (49) was used for calculating the pKa of Lys residues of Thal (PDB code 7CU1) at pH 7.2. The calculation results showed that the pKa values of Lys79 in both subunits are abnormally lower, 7.48 (Chain A) and 7.59 (Chain B), than other Lys residues which have typical pKa values around 9.0 to 11.0 (Table S4). This pKa value allows Lys79 to be in the fully protonated state and can also act as a proton donor during turnovers (Fig. 8B). Although Arg has a similar chain length compared with Lys, its higher pKa compared with Lys would not promote it to act as a proton donor. Therefore, our study implies that the conserved Lys in FDHs may participate in the reaction by providing a proton for HOX to increase the electrophilicity of X (Fig. 8B). Our data do not support formation of a covalent intermediate (chloramine) which was previously proposed for the RebH reaction (27). Results reported here agree with the model proposed in the theoretical study of PrnA (32) and provides more insights into the functional roles of the conserved Lys.</p><!><p>Kinetics of Thal:FADH-∗ inactive complex formation. An anaerobic solution of FADH− was premixed with an anaerobic solution of Thal at various age times (incubation times) (0.01–30 s) in the first mixing. Aerobic buffer was then added into the second mixing step. Final concentrations were 128 μM oxygen, 15 μM FADH−, and 30 μM Thal. Kinetic traces were monitored at 380 nm (solid lines) and 450 nm (solid lines with circles). The inset was a plot of absorbance decrease at 380 nm versus incubation time. FADH−, reduced flavin adenine dinucleotide.</p><!><p>The results in Figure 9 indicate that the amount of absorbance increase at 380 nm during the first phase diminished upon prolonged incubation. These data indicate that a binary complex of active Thal:FADH− changes over time to form an inactive complex, Thal:FADH-∗, which does not allow formation of C4aOOH-FAD (Fig. 3). A plot of the decrease in absorbance at 380 nm versus incubation time represents the kinetics of the formation of the inactive complex (Fig. 9). The data can be fitted with a single exponential decay (inset of Fig. 9), which corresponds to a rate of isomerization to form an inactive complex as 0.41 s−1. Results in Figure 9 showed that the shortest preincubation (0.01 s) gave the highest yield of C4aOOH-FAD formation.</p><!><p>Data collection and refinement statistics for Thal structures</p><p>FAD, oxidized flavin adenine dinucleotide; FADH−, reduced flavin adenine dinucleotide.</p><p>Alignment of various conformations of the flavin binding loop of Thal structures. Superposition of conformations of Chain A (green) and Chain B (white) observed in the Thal:FADH− complex (PDB code 7CU2) defined as "loose" and "tight" conformations, respectively, and Chain A (blue) of Thal:FAD:AMP (PDB code 7CU1) defined as "open" conformation. FAD, oxidized flavin adenine dinucleotide; FADH−, reduced flavin adenine dinucleotide.</p><p>Water networks observed in the open conformation of Chain A of Thal:FAD:AMP (PDB code 7CU1).A, residues within 4 Å of water networks (green sphere). B, spherical display of water molecules shows water networks connecting bulk solvent to the FAD binding site. FAD, oxidized flavin adenine dinucleotide.</p><p>Rearrangement of water networks resulting from protein dynamic changes.A, water network arrangements around the flavin site in Chain A of the Thal:FAD:AMP complex (PDB code 7CU1) were analyzed by MD simulations. Arrangement of water networks observed in the structure of the Thal:FAD:AMP complex (PDB code 7CU1). B, rearrangement of water molecules and the green loop after 4 ns of the production phase performed by MD simulation. C, close-up views of interactions of water molecules around the isoalloxazine ring after 4 ns of the production phase by MD simulations. MD, molecular dynamics; FAD, oxidized flavin adenine dinucleotide.</p><p>Mechanism of inactive complex formation in FDHs.A, water networks located around the flavin were identified in the crystal structure of Thal. Preincubation of Thal with FADH− under anaerobic conditions promotes conformational changes resulting in water rearrangement around the flavin. These water molecules may promote proton transfer to N1 of FADH− to produce FADH2 which is less reactive with O2 to form C4aOOH-FAD (inactive complex). B, freshly mixing Thal with FADH− under aerobic conditions can avoid inactive complex formation. C4aOOH-FAD, C4a-hydroperoxyflavin; FAD, oxidized flavin adenine dinucleotide; FDHs, flavin-dependent halogenases; FADH−, reduced flavin adenine dinucleotide; HOX, hypohalous acid.</p><!><p>Because the data in the HOX formation section indicated that Thal could form HOI, we further investigated whether this enzyme could catalyze iodination of tryptophan. Only chlorinated and brominated but not iodinated tryptophan could be detected (Fig. S8). To explore factors controlling halogenation by Thal, we solved the structure of the Thal:TRP complex (PBD code 7CU0, 2.05–1.95 Å) which has a better resolution compared with the previously reported structure (PDB code 6H44, 2.62–2.55 Å) (35).</p><p>Structures of Thal and various FDHs including PrnA and PyrH (15, 37) were overlaid and analyzed. Superposition of PyrH, Thal, and PrnA structures showed that the three residues (Lys79, Glu358, and Ser359 in Thal) are found at the same position in all of these enzymes (Fig. S20). Moreover, the halogenation sites of the substrate are all placed next to the side-chain of Lys79. We noticed that the space between the tryptophan substrate and Lys79 was confined by the surrounding residues in the pocket and were not large enough to accommodate the larger sphere of the iodide species. We modeled HOI into the active site and found that the van der Waals sphere of the tryptophan does not allow this putative HOI to fit (data not shown). The tight space between the halogenation site of tryptophan and the residues lining the binding pockets would not allow HOI to be accommodated. These data clearly explain the results in the HOX formation section which show that although HOI can be formed at the flavin site, iodination of tryptophan cannot occur because of the limited space in the tryptophan binding site.</p><!><p>Comparison of structures of Thal (PDB code 7CU0) and VirX1 (PDB code 6QGM).A, alignment of VirX1 (purple) and Thal (gray) structures identified the distinct α-helical (yellow) structure encompassing the Trp substrate (green stick) in Thal. B, surface representation of VirX1 aligned with Thal.</p><!><p>Our results reported herein highlight four new key mechanistic insights on FDHs: (A) C4aOOH-FAD formed by tryptophan halogenases can react with I−, Br−, and Cl− to form HOI, HOBr, and HOCl, respectively. However, iodination of tryptophan cannot occur because there is not enough space to accommodate the HOI in the tryptophan binding site. (B) Using D-luciferin for detection of free HOX, we observed leakage of HOX from Thal, and this feature probably contributes to the low product coupling ratio of the enzyme. (C) Site-saturation mutagenesis showed that the active site Lys is indeed crucial because none of the other 19 amino acids can replace the function of Lys. Computational results indicate that the active site Lys has an abnormally lower pKa compared with other Lys residues, possibly allowing it to act as a proton donor for HOX during halogenation. (D) Stopped-flow experiments have shown that prolonged incubation of Thal with FADH− leads to formation of an inactive enzyme species in which C4aOOH-FAD cannot be detected. Various conformations of newly solved Thal:FADH− (PDB code 7CU2) and Thal:FAD:AMP (PDB code 7CU1) structures combined with MD simulations revealed that protein conformational changes and water molecule arrangement may lead to inactive complex formation.</p><p>Our stopped-flow results clearly showed that all tryptophan halognases (Thal, PrnA, and PyrH) can react with all halide ions (I−, Br−, and Cl−) except F− to yield the C4aOH-FAD intermediate and a halogenating species, i.e., HOX. Rate constants of the formation of the halogenating species showed that I− gave the highest rate constant. These data are well correlated with the QM/MM calculations, indicating that the potential energy barrier to form HOX is the lowest for the reaction of C4aOOH-FAD with I−. Although HOI can be formed at the flavin site, an iodinated tryptophan could not be detected as a product. This is likely because of the constrained space in the tryptophan binding site that does not accommodate the larger iodide species.</p><p>The site-saturation mutagenesis and computational studies reported here can explain the mechanistic role of the catalytic Lys. Although all FDHs discovered to date have a strictly conserved Lys at their active sites, it is still unclear whether the residue can form a chloramine intermediate (27) or hydrogen bonds with HOX (29, 32). Only site-directed mutation to Ala (15, 27) or Arg (31) were investigated in the past. Our site-saturation mutagenesis results show that all variants are inactive, indicating that the catalytic Lys is not dispensible. We found that the catalytic Lys is unique for its abnormal pKa value and has the lowest pKa compared with other Lys residues. We propose that this catalytic Lys may serve as a proton donor to HOX (Fig. 8B). Our findings agree with the mechanism previously proposed that Lys does not form a chloramine intermediate (Fig. 8B) (29, 32).</p><p>Kinetics and crystal structures of Thal have also shown that the Thal structure is highly dynamic. Upon prolonged incubation, the Thal:FADH− complex can isomerize to an inactive form which C4aOOH-FAD cannot be detected (Figs. S5 and S6). This is similar to the property of RebH which could also form the inactive complex after prolonged incubation (28). Similar phenomena were also found in several two-component flavin-dependent monooxygenases. A dechlorinating/denitrating flavin-dependent monooxygenase (HadA) from R. pickettii DTP0062 also shifts to form a species which is less reactive with oxygen upon longer preincubation of FADH− and HadA (38). FMNH−-dependent alkanesulfonate monooxygenase (SsuD) from E. coli has also been reported to form an inactive complex when SsuD is preincubated with FMNH− in the absence of oxygen (55). However, some flavin-dependent monooxygenases, i.e., the monooxygenase component (C2) of p-hydroxyphenylacetate hydroxylase, does not have any problem when the enzyme is preincubated with FMNH− (39). As shown in Figure 10, Figure 11, Figure 12, Figure 13 and mentioned in the Results section, Thal protein conformational changes can cause water rearrangements which may lead to an improper protonation status of FADH− or disrupt an arrangement which facilitates the PCET process for C4aOOH-FAD formation. It is interesting to note that enzymes which are robust in formation of C4aOOH-flavin such as C2 use an active site residue (His396 in C2) not water networks as a proton donor for the PCET process. We speculate that this arrangement can be more precise in controlling the protonation status of reduced flavin and the PCET process than Thal which relies on water networks to transfer protons.</p><p>In conclusion, our data have expanded the understanding of the mechanism of halide ion reactivity of tryptophan halogenases. Although tryptophan halogenases can form reactive halogenating species, HOI, HOBr, and HOCl, the enzyme can only catalyze chlorination and bromination but not iodination because the space at the tryptophan binding site cannot accommodate the halogenating iodine species. We believe that future efforts to widen the substrate binding site may allow some tryptophan halogenases such as Thal to catalyze the iodination reaction. We also found HOX leakage which might be the cause for the low catalytic inefficiency of FDHs in general. Future efforts in enzyme engineering to decrease HOX leakage may improve the halogenation activity of FDHs. The catalytic Lys has an abnormally low pKa which allows it to serve as a proton donor during product formation. The fundamental catalytic properties reported should be useful for improvement of halogenases to serve as green biocatalysts in the future.</p><!><p>Single turnover of tryptophan bromination activity of Thal was carried out in a Hi-Tech Scientific Model RQF-63 rapid-quench flow apparatus at 25 °C. An air-saturated solution containing Thal, NaBr, and D,L-tryptophan in 20 mM potassium dihydrogen phosphate buffer pH 7.2 was mixed with an anaerobic solution of FADH−, and the reactions were stopped at different time points by quenching the reaction with 0.4 M HCl. Final concentrations after quenching were 15 μM Thal, 7.5 μM FADH−, 15 μM D,L-tryptophan, and 5 mM NaBr. Reactions were analyzed by HPLC-MS/MS using an Eclipse Plus C18 column (3.5 μm, 4.6 × 100 mm) with a gradient mobile phase of 95:5 to 20:80 of H2O:acetonitrile containing 0.1% (v/v) formic acid as an eluent at a flow rate of 0.5 ml/min and monitoring for mass/charge (m/z) of tryptophan and bromotryptophan. Multiple reaction monitoring mode of HPLC-MS/MS with positive ion mode was set for the dectection.</p><!><p>Rapid kinetics experiments were performed using a TgK Scientific model SF-61DX stopped-flow in single- or double-mixing mode or with a model SF-61SX stopped-flow spectrophotometer in single-mixing mode. An oxygen-scrubbing solution, as previously described (39), was used for scrubbing trace amounts of oxygen from the apparatus to make the system anaerobic. Reactions were carried out in 20 mM potassium dihydrogen phosphate buffer, pH 7.2, 25 °C. Trace amounts of oxygen were removed from the glass tonometer by evacuating air and purging the system with nitrogen gas. FAD was reduced in a glove box (Belle Technology, Weymouth, UK) to obtain FADH− by titration with sodium dithionite. Flavin intermediates were monitored by UV-Visible or fluorescence spectrometry. Apparent rate constants (kobs) were analyzed using Program A software (written by R. Chang, J. Chiu, J. Dinverno, and D. P. Ballou, University of Michigan). More details for each experiment are described in Supplemental Information.</p><!><p>The reactivity of C4aOOH-FAD with different halide ions to form C4aOH-FAD and HOX was investigated using QM/MM. The crystal structure of the Thal:FAD:AMP complex (PDB code 7CU1) was used for these calculations. The model setup and the QM/MM calculations are described in the Supporting Information. We applied protocols similar to those used previously in QM/MM investigations of other enzyme-catalyzed reaction mechanisms (43, 44). The positions of hydrogen atoms in the enzyme were located using the CHARMM procedure HBUILD (56). The protonation states of titratable amino acid residues were chosen based on predictions from PROPKA (57). The atom types in the topology files were assigned based on the CHARMM27 parameter set (58). The system was divided into QM and MM parts for QM/MM simulations. The QM part consisted of Cl− and C4aOOH-FAD with a net charge of 0 e. The model was divided into the reaction zone which is an area of 21 Å radius, centered on the carbon atom C4a of C4aOOH-FAD molecule, and the buffer zone is the rest of the protein system. First, the energy system was minimized using Adopted Basis Newton-Raphson minimization with the AM1/CHARMM27 QM/MM method. The system was then equilibrated with QM/MM molecular dynamics at the AM1/CHARMM27 level. MP2/6-31G(d),I=lanl2dz/CHARMM27 calculations on the reaction (59, 60, 61) were performed with the QoMMMa program (62). The optimized structure from this ab initio QM/MM reaction modeling calculations was used for higher level ab initio QM/MM single point calculations using the Molpro program (63). Energy profiles were calculated at the SCS-MP2/aug-cc-pvdz,I=auc-cc-pwcvdz-pp/CHARMM27 level; the SCS-MP2 method has been shown to give accurate results, close to coupled cluster calculations for other enzyme-catalyzed reactions (43, 44, 45, 46). More details of the calculations are given in the Supporting Information.</p><!><p>Experiments were done by performing multiple turnover reactions in the absence of substrate tryptophan to produce HOX. In this study, HOBr was generated after 1.5 h of incubation with 20 mM potassium dihydrogen phosphate buffer pH 7.2, 20 μM FAD, 20 μM NAD+, 100 mM NaBr, 0.4 μM flavin reductase (C1), 25 mM HCOONa, 0.2 μM formate dehydrogenase from Pseudomonas sp. (psFDH), and 20 μM Thal at room temperature (22–25 °C). Enzymes in the reaction were removed by centrifugation through Microcon (10 kDa cut-off). D-luciferin (2 μM final concentration) was added into the reaction solution. After incubation for 2 h, reactions were analyzed using HPLC-diode array detector/MS and HPLC-MS/MS.</p><!><p>A large scale reaction (6.3 ml) of D-luciferin bromination was set up to characterize the products. The reaction was comprised of 20 μM FAD, 20 μM NAD+, 100 mM NaBr, 0.3 μM C1, 25 mM HCOONa, 0.2 μM psFDH, 5 μM D-luciferin, and 20 μM Thal (wild-type) in 20 mM potassium dihydrogen phosphate buffer pH 7.2. Reactions were incubated for 48 h at room temperature (22–25 °C) in a dark room. Fresh enzymes including Thal, C1, and psFDH were added into the reaction mixtures every 4 h. Enzymes in the reaction were removed by centrifugation through a Microcon (10 kDa cut-off). Reactions were analyzed using HPLC-QTOF.</p><!><p>For the site-saturation mutagenesis at the Lys79 position, NNK-degenerate oligonucleotides were designed as showed in Table S2. Using PCR amplification, Thal (wild-type) containing GFP11 tag in the pET28b vector was used as a template (47). PCR products were purified from an agarose gel, and the products were treated with DpnI and cutsmart buffer at 37 °C for 16 h to remove the Thal (wild-type)-pET28b plasmid. Plasmids were transformed into BL21(DE3) competent cells by the heat shock method. LB broth (900 μl) was added to transformed BL21(DE3) cells. Cell cultures were grown at 37 °C and 220 rpm for 2 h before centrifugation and discarding 900 μl of supernatant. The remaining 100 μl were resuspended and spreaded onto 15 ml of LB agar containing 25 μg/ml kanamycin. Cultures were further grown at 37 °C for 16 h. Cells were harvested by resuspending colonies using 3 ml of LB, and extracted DNA was used for sequence analysis to evaluate site-saturation libraries (64). Libraries which had Qpool values of at least 0.7 were further cultured on LB agar plates containing 50 μg/ml kanamycin for 16 h. A total of 168 colonies of Lys79 variants, 12 colonies of wild-type, six colonies of BL21(DE3) as negative controls, respectively, were cultured in 96-well plates containing TB medium with and without kanamycin. Cells were cultured at 37 °C with shaking at 50 rpm for 16 h. Cells were divided into two parts. The first part was kept as a library stock by mixing each cell cuture with glycerol to make up a final concentration of 25% (v/v) and stored at −80 °C. The second part was used for assay of protein level and tryptophan bromination.</p><p>Gene expression was done by transferring 7 μl of a starter culture into 200 μl of TB medium (3% v/v inoculation). Cells were cultured at 37 °C with shaking at 50 rpm for 7 h. Then, isopropyl β-D-1-thiogalactopyranoside with a final concentration of 1 mM was added into the culture, and the temperature was switched from 37 to 25 °C. Cells were further cultured at 25 °C with shaking at 50 rpm for 16 h. Cells were harvested and lysed. Lysates were centrifuged for 30 min at 4000 rpm 4 °C. The supernatants were taken for activity assay and protein level analysis.</p><p>Gene expression yields were determined using high throughput GFP split system (47, 48). Thal was tagged with a small portion of GFP sequence called a GFP11 at the C-terminal (Thal-GFP11). The purified GFP1-10 (150 μl) was incubated with 20 μl of supernatant at 4 °C for 1 day. The protein levels were determined by measuring fluorescence intensity at an emission wavelength of 530 nm using an excitation wavelength of 488 nm which resulted from the interaction of GFP1-10 and Thal-GFP11.</p><p>The supernatant (40 μl) was transferred to new 96-well plates. A solution of 12.5 mM NaBr, 10 μM NAD+, 10 μM FAD, 22 mM HCOONa, 0.1 mM tryptophan, and 0.1 μM C1 in 20 mM potassium dihydrogen phosphate buffer pH 7.2 was added (50 μl) into each well. The reaction was initiated by the addition of 0.165 μM psFDH to each well, and the reactions were incubated at 60 rpm, 25 °C. After the reactions were incubated for 3 h, they were quenched with 0.2 M HCl and analyzed using HPLC-MS/MS with the same conditions as described in the single turnover reaction section.</p><!><p>Crystals of Thal apoenzyme and tryptophan complexes were obtained from a microbatch crystallization (65) in a solution containing 0.28 to 0.38 M NaF, 23% to 31% (w/v) PEG3350. Thal:FAD:AMP and Thal:TRP crystals were obtained from soaking of apo Thal crystals into a well containing crystallizing solution and FAD and Trp at a ratio of 2.2:5 mM for 7 days and 0.8:1.6 mM for 2 days, respectively at 25 °C. Thal:FADH− co-complex crystals were produced by placing Thal:Trp complex crystals in a well containing crystallizing solution containing 5 mM tryptophan, 4 mM FAD, and 20% (v/v) glycerol at 25 °C followed by the addition of Na2S2O4 powder. The crystals were harvested by flash-freezing in liquid nitrogen. Diffraction data were collected from the flash frozen crystals on a single crystal X-ray diffraction D8 venture at the Characterization and Testing Center, National Science and Technology Development Agency, Thailand. Data processing was performed with Proteum3 software. Structural phases were solved by molecular replacement using PHASER (66) in the CCP4 suite (67) with SttH (PDB code 5HY5) as a search model. Structure refinement and building were performed with REFMAC5 (68) and COOT (69, 70). Parameter files of FADH− for refinement was generated with LIBCHECK in the CCP4 suite. The structures were validated in PROCHECK (71). LigPlot (72) and PyMOL (https://pymol.org/) was used for molecular graphics.</p><p>Figures were created with ChemBioDraw Ultra (Version 15.1), BioRender.com, and the PyMOL Molecular graphic System (Version 2.3.4).</p><!><p>All data are contained within the manuscript and Supporting Information. Data of X-ray structures are available at Protein Data Bank under PDB codes indicated.</p><!><p>The authors declare that they have no conflict of interest with the contents of this article.</p><!><p>Supplementary Figure and Tables</p>