A Novel F420-dependent Thioredoxin Reductase Gated by Low Potential FAD

A recent report suggested that the thioredoxin-dependent metabolic regulation, which is widespread in all domains of life, existed in methanogenic archaea about 3.5 billion years ago. We now show that the respective electron delivery enzyme (thioredoxin reductase, TrxR), although structurally similar to flavin-containing NADPH-dependent TrxRs (NTR), lacked an NADPH-binding site and was dependent on reduced coenzyme F420 (F420H2), a stronger reductant with a mid-point redox potential (E′0) of −360 mV; E′0 of NAD(P)H is −320 mV. Because F420 is a deazaflavin, this enzyme was named deazaflavin-dependent flavin-containing thioredoxin reductase (DFTR). It transferred electrons from F420H2 to thioredoxin via protein-bound flavin; Km values for thioredoxin and F420H2 were 6.3 and 28.6 μm, respectively. The E′0 of DFTR-bound flavin was approximately −389 mV, making electron transfer from NAD(P)H or F420H2 to flavin endergonic. However, under high partial pressures of hydrogen prevailing on early Earth and present day deep-sea volcanoes, the potential for the F420/F420H2 pair could be as low as −425 mV, making DFTR efficient. The presence of DFTR exclusively in ancient methanogens and mostly in the early Earth environment of deep-sea volcanoes and DFTR's characteristics suggest that the enzyme developed on early Earth and gave rise to NTR. A phylogenetic analysis revealed six more novel-type TrxR groups and suggested that the broader flavin-containing disulfide oxidoreductase family is more diverse than previously considered. The unprecedented structural similarities between an F420-dependent enzyme (DFTR) and an NADPH-dependent enzyme (NTR) brought new thoughts to investigations on F420 systems involved in microbial pathogenesis and antibiotic production.

In this report, we show that the electron delivery enzyme of the above-mentioned redox regulation system of M. jannaschii is tuned for operation under highly reduced conditions. This archaeon carries two Trxs, Mj-Trx1 which is a typical Trx, and Mj-Trx2 with unknown activity, and a TrxR (Mj-TrxR) that is an NTR homolog yet does not use NAD(P)H (20). We demonstrate that Mj-TrxR uses coenzyme F 420 , a deazaflavin derivative present in all methanogens, as electron carrier. We reveal novel redox properties of this deazaflavin-dependent flavincontaining thioredoxin reductase (DFTR) and present the respective ecological implications. Also, the structural and phylogenetic analyses of DFTR homologs bring up new possibilities with the long studied and widespread disulfide oxidoreductase family of flavoenzymes and the expanding field of F 420 -dependent enzymes. Similarly, the properties and distribution of DFTR and the conditions prevailing in M. jannaschii's habitat lead to the proposal that this enzyme developed in the low redox environment of early Earth and gave rise to NTR. This implication fits the hypothesis that the Trx-based redox control developed on Earth long before oxygen appeared (20).

Discovery of a Deazaflavin-dependent Thioredoxin Reductase
The thioredoxin reductase of M. jannaschii (Mj-TrxR) was expressed recombinantly in E. coli with an NH 2 -terminal His 6 tag, and it was soluble. The recombinant protein was purified to homogeneity via Ni 2ϩ -affinity chromatography, and the His 6 tag was removed via TEV-protease digestion. The UV-visible spectra of the purified protein with absorbance maxima at 380 and 460 nm indicated that it contained flavin (Fig. 1A) (25). However, it had a partial incorporation of this cofactor, and incubation with FAD led to full incorporation. The denatured molecular mass of the protein as determined via SDS-PAGE was found to be 33 kDa (Fig. 1B), and the value of the apparent native molecular mass as estimated via size exclusion chromatography was 63 kDa. Thus, the native Mj-TrxR was likely a homodimer of 33-kDa subunits, which is a typical property of low molecular weight flavin-containing thioredoxin reductases (1)(2)(3)(4)25).
The initial test for the Trx reductase activity in Mj-TrxR was via the commonly used insulin disulfide reduction assay (26), and it was performed under anaerobic conditions (20,27). It utilized the following electron flow: electron donor 3 Mj-TrxR 3 Mj-Trx 3 insulin. In this assay, Mj-TrxR did not reduce either Mj-Trx1 or Mj-Trx2 with NADH, NADPH, coenzyme M, or reduced ferredoxin as electron donor. However, with dithionite, an effective artificial electron donor for NTR (28), Mj-TrxR reduced Mj-Trx1 but not Mj-Trx2. Thus, Mj-TrxR was a bona fide Trx reductase and recognized Mj-Trx1 as a substrate. UV-visible spectra of Mj-DFTR mixtures incubated with various electron donors are shown in Fig. 2. Then, we found that with reduced coenzyme F 420 (F 420 H 2 ), Mj-TrxR reduced Mj-Trx1. The details for this assay appear below. Coenzyme F 420 is a 7,8-didemethyl-8-hydroxy-5-deazaflavin derivative, and accordingly, we called the Mj-TrxR a DFTR (Fig. 3).

Kinetic Properties of Mj-DFTR
The enzyme activity was assayed by following the oxidation of F 420 H 2 spectrophotometrically at 400 nm; F 420 , but not F 420 H 2 , has a detectable absorbance at the 360 -450 nm range with a maximum at 420 nm, and 401 nm is the isosbestic point with respect to pH (29). Concentration of Trx was held constant by coupling the assay to the reduction of oxidized glutathione or GSSG (30). The flow of electrons in the assay was as follows: F 420 H 2 3 DFTR 3 Trx 3 GSSG. The DFTR activity was not observed in the absence of Trx or TrxR in the assay. The activity of the enzyme increased with a rise in temperature up to 90°C; the assays at higher temperatures could not be performed due to technical limitations. The activation energy of the DFTR-catalyzed reaction as determined from the slope of the linear section of the Arrhenius plot (50 -90°C) was 36.14 kJ mol Ϫ1 (Fig. 4A). The optimal pH for activity was 8.0 ( Fig. 4B).
At a fixed F 420 H 2 concentration of 60 M and varying concentrations of Mj-Trx1 (0.5-40 M), the apparent K m for Mj-Trx1 was found to be 6.3 Ϯ 2 M (Fig. 5A). Under similar conditions and with Mj-Trx2 concentrations of 25-500 M, the value of K m for Mj-Trx2 was 371 Ϯ 211 M (Fig. 5B). Assays with 15 M Mj-Trx1 and 4 -60 M F 420 H 2 , yielded an apparent K m value of 28.6 Ϯ 2.5 M for F 420 H 2 (Fig. 5C). The specific activity of the enzyme with 60 M F 420 H 2 and 30 M Mj-Trx1 was 29.7 Ϯ 1.3 mol/min/mg. The reactivity of DFTR toward Mj-Trx2 was poor, and the respective V max value was 3-fold lower than that observed with Mj-Trx1. Also, unlike most Trxs, Mj-Trx2 does not reduce insulin, and in contrast to Mj-Trx1, it is not recognized by Escherichia coli NTR (20). Thus, it is unclear whether this protein is a true Trx and how it is reduced or oxidized.  (32). F 420 showed absorbance maxima at 420 nm (32). FIGURE 3. Electron flow in the deazaflavin-dependent thioredoxin reductase system of M. jannaschii. F 420 H 2 -generating systems: F 420 -reducing hydrogenase (Frh) (93); Hmd-Mtd cycle, H 2 -dependent methylenetetrahydromethanopterin dehydrogenase (Hmd) ϩ F 420 H 2 -dependent methylenetetrahydromethanopterin dehydrogenase (Mtd) (94). Inset, structures of coenzyme F 420 and F 420 H 2 . *, the Trx targets and Trx controlled systems were identified in proteomics study (20).
value of DFTR for DTNB of 182.2 Ϯ 70 M was obtained (Fig.  5D). The specific activity of DTNB reduction at F 420 H 2 and DTNB concentrations of 60 and 1000 M, respectively, was 26.8 Ϯ 0.4 mol/min/mg. Unlike the high molecular mass NTRs, DFTR did not reduce lipoate and selenite, and low molecular mass NTRs in general behave similarly (2,3).

Redox Properties of Mj-DFTR, E 0 Values for the Redox Centers
Similar to prokaryotic NTRs that are composed of two identical approximate 35-kDa subunits (32), the DFTR is expected to have two redox centers as follows: a protein-bound FAD and a catalytically active Cys pair residing within a CXXC motif (Fig.  6). We determined the mid-point redox potential (EЈ 0 ) values for these two centers via redox titration as follows.
(i) For the EЈ 0 for the Cys-disulfide/Cys-SH pair of the catalytically active CTMC motif of Mj-DFTR, this value was determined via the measurement of the fluorescence emission intensities of mBBr-labeled forms of the enzyme preparations that were subjected to incubation in a series of redox buffers offering a potential range of Ϫ350 to Ϫ200 mV (Fig. 7A). The data fit the Nernst equation for a 2-electron reduction reaction and yielded an EЈ 0 value of Ϫ279 Ϯ 7 mV at pH 7.0 ( Fig. 7A) which com- pared well to the reported values of 230 -277 mV for the EЈ 0 of Trxs from a variety of sources (33).
(ii) For the EЈ 0 of the Mj-DFTR-bound flavin cofactor, our first attempt to obtain this value involved redox titration of Mj-DFTR with F 420 H 2 , which revealed a complex picture for the interaction between this substrate and the enzyme that has been elaborated under the "Discussion." The second attempt involved an artificial reductant and provided the EЈ 0 value. We describe both below.
Titration with F 420 H 2 -Because Mj-DFTR carried only one catalytically active Cys pair, it was expected to follow the mechanism reported for E. coli NTR (34) and allow the determination of the EЈ 0 value from the equilibrium constant (K eq ) for the Mj-DFTR reaction (E ox ϩ F 420 H 2 3 EH 2 ϩ F 420 , where E is Mj-DFTR) and the EЈ 0 of the F 420 /F 420 H 2 pair, which is Ϫ360 mV (36) (E ox , oxidized enzyme; EH 2 , 2-electron reduced form of the enzyme). We attempted to obtain the K eq value from the changes in the UV-visible spectrum of Mj-DFTR upon redox titration with F 420 H 2 (Fig. 7B). To avoid the wavelengths where all or most of the participating cofactors exhibit significant absorbance values (29,37), the A 480 and A 540 values (where A is absorbance) were chosen for the determination of the concen-trations of FAD and the FAD-thiolate charge transfer species from which similar information on other species could be calculated. Fig. 7C presents the resulting data in terms of respective molar absorbance values, ⑀ 480, EH2/Eox and ⌬⑀ 540, EH2/Eox . It was observed that both values attained saturation only after the addition of 2 eq of F 420 H 2 , and as judged from the 480 nm data, even at this stage only 1 eq of FAD was reduced. As elaborated under the "Discussion," although these observations brought up interesting aspects of the enzyme, it also highlighted the need for a method that will involve a non-physiological reductant and an artificial reporting system and thereby would allow direct determination of the EЈ 0 value of Mj-DFTRbound FAD. Accordingly, an alternative method, as presented below, was used.
Titration with Dithionite in the Presence of Methyl Viologen (MV) as Redox Reference-The primary basis for this measurement was that at any point of the titration, the EЈ of Mj-DFTR-FAD/Mj-DFTR-FADH 2 (EЈ Mj-DFTR-FAD/Mj-DFTR-FADH2 ) and the MV ϩ /MV 0 (EЈ MV ϩ /MV 0) would be the same. The reduction of FAD and methyl viologen (MV ϩ ) to FADH 2 and MV 0 (1-electron reduced methyl viologen) with various amounts of dithionite was followed spectrophotometrically; the former caused a decrease in A 460 and the latter raised A 604 . The concentrations of reduced MV 0 and enzyme-bound FAD were determined from A 604 and A 460 values, respectively, employing the following extinction coefficient values: MV 0 , 12 mM Ϫ1 cm Ϫ1 at 604 nm and pH 7 (a value determined experimentally using the published value of 13.5 mM Ϫ1 cm Ϫ1 at pH 8 (38)); enzyme-bound FAD, 12.7 mM Ϫ1 cm Ϫ1 at 460 nm and pH 7, determined with oxidized enzyme carrying one FAD/subunit. With the concentrations of MV ϩ and MV 0 , the EЈ values of the titrated mixtures were calculated from the Nernst equation for a 1-electron transfer reaction (Equation 1), where EЈ 0 of the MV ϩ /MV 0 pair was Ϫ446 mV (38). Then the equilibrium concentrations of FADH 2 versus EЈ data were fitted to Equation 2, which considers a 2-electron transfer reaction (n ϭ 2), yielding an EЈ 0 value of Ϫ389 Ϯ 24 mV for FAD-bound Mj-DFTR (Fig.  7D).

Primary Structure Features of DFTR
Mj-DFTR was composed of 301 amino acids with theoretical molecular mass and pI values of 33 kDa and 6.86, respectively, which are typical of low molecular mass NTRs (2,3). At the amino acid sequence level, it was more similar to flavin containing TrxRs of Thermoplasma acidophilum (19), a facultative anaerobic archaeon (39), and Thermotoga maritima (16), an anaerobic bacterium (40), than to E. coli NTR (41) and ferredoxin-dependent and flavin containing TrxR of C. pasteurianum, an anaerobic fermentative bacterium (8) (Fig. 6); the respective percent identity and similarity values appear in the legend of Fig. 6. The nature of the reductant used by TrxR of T. acidophilum is unknown (19), and the enzyme of T. maritima is an NTR of dual specificity for NADH and NADPH (16); both organisms are thermophiles (39,40). Fig. 6 presents an amino acid sequence alignment of Mj-DFTR with the abovementioned flavin containing TrxRs (8,16,19,41) that leveraged the x-ray crystallographic structures of the E. coli and T. acidophilum enzymes (PDB codes 1CL0 and 3CTY, respectively (19,41)).
The defining features of NTR are the FAD-and NADPHbinding elements and an active site CXXC motif with two redox-active cysteine residues (X, non-conserved amino acid residue), and Mj-DFTR mimicked NTR in the FAD-binding elements but not in the rest (Fig. 6). The NTR interacts with FAD via a dinucleotide-binding motif (DBM FAD ) encompassing three conserved elements (GXGXXG, GG, and ATG) and two additional COOH-terminal binding sequences (GD and a G-helix) (42). Mj-DFTR contained all these features (Fig. 6). The NADPH-binding elements of NTR (DBM NADPH ) are GXGXX(G/A), GGGXXA, and HRR motifs, and the last ele-ment is for binding the 2Ј-phosphate group of the coenzyme. In Mj-DFTR, the sequence elements equivalent to GXGXX(G/A) and GGGXXA of NTR were found to be GRGISY and GRDTPA, respectively, representing a replacement of the glycine/alanine residue by a larger aromatic and hydrophobic amino acid for the first and a substitution of two glycine residues with two charged residues for the second (Fig. 6). In addition, the 2Ј-phosphate-binding HRR residues of NTRs were absent in Mj-DFTR. Another special feature was found at the active site of Mj-DFTR, which carried a CTMC motif, whereas E. coli NTR carries CATC at this location (41); the respective element in T. acidophilum TrxR is CSTC (19).
As true for the majority of F 420 -dependent enzymes (43), the coenzyme F 420 -binding amino acid residues of Mj-DFTR could not be identified from primary and secondary sequence alignments. It is known that the features of the F 420 -binding sites of different groups of F 420 -dependent enzymes are different (43)(44)(45).

Distribution and Phylogenetic Clades of DFTR Homologs
Because of high sequence similarities between flavin-containing thioredoxin reductases and other pyridine dinucleotide oxidoreductase family proteins (46), especially the alkyl hydroperoxide reductases, we employed a two-step method involving the following for the identification of DFTR homologs: identification of candidates via Psi-BLAST searches (47) using Mj-DFTR (locus number, MJ1536) as a query and with stringent parameters (low e-value, 1e Ϫ25 ; 10 iterations); phylogenetic analysis-based screening of the candidates. Such an analysis of the extracted methanogen proteins suggested that within this group the DFTR homologs were limited to the phylogenetically deeply rooted methanogens belonging to the class of Methanococci under the euryarchaeal phylum, and these organisms lacked NTR and ferredoxin-thioredoxin reductase (FTR) homologs ( Fig. 8 and supplemental Table S1). Other methanogens carried NTR, and some of these harbored one or more FTR homologs as well ( Fig. 8 and supplemental Table S1). A close homolog of DFTR was not found in the genomes of bacteria and non-methanogenic archaea. However, as discussed below, the study suggested that some of the yet to be studied TrxRs have the potential of using electron donors other than NADPH. Fig. 9 presents the results from a phylogenetic analysis of the NTR homologs (including DFTR) of all methanogens, selected non-methanogenic archaea, known non-canonical (non-NA-D(P)H utilizing) flavin containing TrxRs, and selected TrxRs of bacteria and archaea that have been studied via direct assays; some information about these proteins appear in supplemental Table S2. In this analysis the homologs of DFTR and ferredox in-dependent clostridial TrxRs formed two tightly clustered groups (marked with F 420 H 2 and Ferredoxin, respectively, in Fig. 9) that seemed to have diverged recently from a common ancestor. The homologs of the non-canonical TrxR of T. acidophilum with unknown electron donor specificity (19) also formed a cohesive clade (Non-NAD(P)H I, Fig. 9). However, it was closely related to an NTR group, which included NADPHdependent TrxR of Lactobacillus casei (protein accession number, WP_012491289.1), which is a facultative anaerobic bacte-rium ( Fig. 9). The study revealed four more clades of potentially non-canonical enzymes. The non-NAD(P)H II group contained TrxR of D. vulgaris (D.vulgaris_3 in Fig. 9; protein accession number AAS95935.1) that does not utilize NAD(P)H and uncharacterized TrxRs of halophilic archaea, whereas the non-NAD(P)H III clade included experimentally validated Sulfolobus solfataricus TrxR-1 (S. solfa_1, protein accession number WP_010923960.1) and certain unstudied enzymes of methanogens from the Methanomicrobia class (Fig. 9). Two additional groups of methanogen TrxRs of non-canonical types (Unknown I and Unknown II, Fig. 9) were also identified, and none of these has been examined experimentally. However, their relationships with the phylogenetically nearest neighbors suggested that the TrxRs of the Unknown I group are likely to be NTRs. The T. maritima TrxR that can use both NADH and NADPH (16) belonged to an NTR group containing bona fide enzymes from D. vulgaris (D.vulgaris_2, accession number AAS94860) (6) and Methanosarcina acetivorans (accession number, AAM04784) (9).
All known and putative non-canonical TrxRs lacked the phosphate-binding HRR motif, and most of them contained non-E. coli NTR-type CXXC motifs. However, the candidate non-canonical TrxRs could not be analyzed for their potentials of utilizing F 420 as electron carrier, because the F 420 -binding residues of DFTR are not known, and the F 420 -binding residues do not show broad conservation across F 420 -dependent enzymes families (43).
The case of Methanopyrus kandleri, a strictly hydrogenotrophic methanogen that lives in deep-sea hydrothermal vents and grows at temperatures as high as 110°C (48), raised an intriguing question. It does not carry a Trx or DFTR but an NTR-type TrxR (Fig. 8), offering an opportunity of discovering a new TrxR reductase substrate. For evolutionary deductions, this enzyme presents a problematic case, because in terms of 16S rRNA sequence M. kandleri appears closely related to the Methanococci (Fig. 8), and the whole genome-based analysis contradicts this position (49). FIGURE 8. Thioredoxins and thioredoxin reductases in representative methanogens. The information is presented on a maximum likelihood inferencebased 16S-rRNA phylogenetic tree. Dash, not detectable in homology based searches; black bullet at the branches, a confidence value Ͼ70%; scale bar, number of base substitution per site; Trx, thioredoxin; NTR, FTR, and DFTR, NADPH-, ferredoxin-, and deazaflavin-dependent thioredoxin reductase, respectively. The NCBI accession numbers for Trxs and Trx reductases are listed in supplemental Table S1.

Discussion
The results described above revealed a novel type of thioredoxin reductase and provided the mechanistic basis for the dependence of the enzyme on coenzyme F 420 and the respective ecological relevance. A comparative study with this enzyme provided new information about the catalytic properties and evolutionary biology of a long-studied and broadly important family of enzymes and the evolutionary development of the widely distributed NTR. It also presents leads for the expanding field of F 420 -dependent enzymes in a broad range of organisms. We elaborate these outcomes below.
Our study began with a bioinformatic analysis that showed that every methanogen carried at least one NTR homolog, and FTR was present only in late evolving and nutritionally diverse species that live in close association with anaerobic bacteria (Fig. 8); only three species contained more than one NTR homolog (Fig. 8). FTR is thought to have originated in bacteria and acquired by the late evolving methanogens through horizontal gene transfer (50). Thus, NTR homologs could be considered a fundamental and ancient component of methanogens. An analysis of the structural features identified the NTR homologs of late evolving methanogens as potentially true NTRs, whereas those of Methanococci, such as M. jannaschii, were non-canonical types (Figs. 6 and 9); as mentioned under the "Results," Methanococci lacked FTR (Fig. 8). M. jannaschii and all other Methanocaldococcus species carried one non-canonical TrxR FIGURE 9. Maximum likelihood phylogenetic analysis of DFTR and selected flavin-containing TrxRs. The phylogenetic tree was constructed as described previously (91). Name in boldface, bona fide TrxR validated by direct activity assay; details are in supplemental Table S1. Black and white bullets near branches: bootstrap values Ͼ70 (calculated from 100 replicates). Each label appearing on the outer side of the tree; electron donor is utilized by the respective clade. Unknown, TrxR with unknown reductant. Non-NAD(P)H, inability to use NAD(P)H as electron donor as determined via direct activity assay of representative members (supplemental Table S1). Symbols for the departures in the amino acid sequences at select conserved elements of TrxRs as described in Fig. 6: filled black and gray circles, unrecognizable HRR and GGG motifs, respectively; open circle, absence of the GG motif. Ec-DLD, dihydrolipoamide dehydrogenase of E. coli used as an outgroup. Details for the abbreviations for the host organism names and the accession numbers of respective TrxRs are listed in supplemental Table S2. and 2-4 Trxs, presenting a minimal system among the methanogens ( Fig. 8; supplemental Table S1), and these organisms live in the deep-sea hydrothermal vents presenting early Earthlike environments that are characterized by low redox potential values and high temperatures (51,52). Therefore, we hypothesized that Methanocaldococcus species carry the history of the development of NTR, and to explore this possibility further, we have investigated the catalytic properties of M. jannaschii TrxR (Mj-TrxR).
Matching the predictions from the structural features (Fig.  6), Mj-TrxR was found to contain bound flavin and did not utilize NAD(P)H. It utilized coenzyme F 420 H 2 , a deazaflavin derivative, an electron donor, and this is the first report of a DFTR. F 420 is present in high abundance in all methanogenic archaea tested (29,53), and the gene for F 0 -synthase, the key F 420 -biosynthesis enzyme (54,55), is present in all methanogen genomes (54,55). The coenzyme is present in trace amounts in Halobacteria and Thermoplasma (56), which are non-methanogenic archaea, and among the bacteria, it is found mostly in the Actinobacteria (57). In methanogens, F 420 is the electron carrier for 12 previously described enzymes (36), and six of these utilize F 420 H 2 (43,44,58,59). The observed K m value for F 420 H 2 (28.6 M) was similar to that recorded for several F 420 H 2 -oxidizing enzymes (5.4 -150 M) (58,59), and a similar agreement was seen between the K m values of DFTR for Mj-Trx1 DFTR (6.3 M) and that of other thioredoxin reductases for their cognate Trxs (0.6 -86 M) (9,15). These values were also physiologically meaningful; the F 420 contents of methanogen cells fall in the range 200 -4000 pmol/mg dry cell weight (36), which roughly translates to an intercellular concentrations of 60 -1200 M, assuming that 70% of the mass of a cell is due to water.
The chemistry of the other analogous F 420 -dependent enzymatic reactions suggested that DFTR would use F 420 for hydride transfer, a role also served by NAD(P) ϩ (60). Because the Mj-DFTR-bound FAD could be reduced with F 420 H 2 (Fig.  3) and reduced enzyme could be oxidized with Mj-Trx1, the enzyme was expected to exhibit a Ping-Pong mechanism involving the above-mentioned two independent half-reactions. A preliminary bi-substrate kinetic analysis (Fig. 10) supports this expectation as the reciprocal plots for the initial velocity data collected with varied concentrations of F 420 H 2 (5-60 M) and Mj-Trx1 at two fixed concentrations (2 and 15 M) produced two parallel lines. Because there were indications of complications arising from inhibition of the enzyme by F 420 H 2 and the generation of single-electron reduced deazaflavin species (37) in the presence of EDTA and visible light used in spectrophotometry, albeit in very minor amounts, the detailed mechanistic studies will be pursued and reported in the future. It should be noted that to avoid such a problem, all redox titrations were performed in the absence of EDTA.
The observed conservation of several key residues of NTR in DFTR (Fig. 6) brought up the following question. Why can't DFTR utilize NAD(P)H and need a low potential electron carrier such as F 420 for activity? To answer this question, we determined the EЈ 0 values for the Cys-disulfide/Cys-SH pair at the catalytic 130 CXXC 133 unit and the protein-bound FAD/FADH 2 system. The former value (EЈ 0 , Ϫ279 mV) did not yield a sur-prise, as it was similar to the EЈ 0 of the Cys-disulfide/Cys pair of E. coli NTR (Ϫ270 mV) (41).
In contrast, the results from the efforts to measure the EЈ 0 value of Mj-DFTR-bound FAD provided the illuminating information about the redox properties of DFTR and the ecophysiological implications for the organisms carrying this enzyme. The first attempt to measure the EЈ 0 value of Mj-DFTR-bound FAD involved redox titration of the oxidized form of the enzyme with F 420 H 2 . It presented several challenges. One was due to overlapping UV-visible spectra of the flavin and deazaflavin species that made measurement at the most responsive wavelengths (400 -470 mm) rather difficult (Fig. 7B). When we opted for analysis at less sensitive wavelengths of 480 and 540 nm, which represented oxidized form of FAD and the charge transfer between FAD and Cys-thiolate at the active site, the results presented a complex picture. The observations at both 480 and 540 nm showed that the full reduction of the enzymebound flavin required 2 eq of F 420 H 2 (Fig. 7C). This was surprising, as even the eukaryotic NTR from Drosophila melanogaster that carries two catalytically active Cys pairs attains the EH 2 state with only 1 eq of NADH, and the consumption of 2 eq amounts of NADH converts this enzyme to an EH 4 form (35). The D. melanogaster NTR cycles between EH 2 and EH 4 states during catalysis, and the UV-visible spectra gathered after redox titration of the protein with NADH exhibit four isosbestic points (35). In contrast, the titration of Mj-DFTR with F 420 H 2 produced two isosbestic points (at 445 and 500 nm; Fig.  7B), which is indicative of an EH 2 species (35). The consumption of 2 eq of F 420 H 2 in part was certainly due to the reduction of the enzyme-bound FAD and the active site Cys-disulfide and was most likely a part of the added F 420 H 2 bound to the enzyme. This is because the enzyme will be loaded with reducing equivalents at all sites, the F 420 H 2 -binding site, FAD, and the active site Cys pair, and as the equilibrium principle dictates, none of these sites will be saturated with the addition of 2 eq of F 420 H 2 . Even if it were possible to measure the concentration of these species, the analysis of the resultant data would be further challenged by the possibility that in the enzyme-bound form the deazaflavin would present an EЈ 0 value that is different from that measured with the free forms of F 420 and F 420 H 2 (Ϫ360 mV) (36). We also found that the changes in the flavin spectrum at 540 nm during the titration with F 420 H 2 (Fig. 7B) are not characteristics of a charge transfer between FAD and Cys thiolate (35). For these reasons, in our second attempt, we used dithionite as the reductant in the titration on Mj-DFTR, and the MV (the MV ϩ /MV 0 redox pair) was the redox reference. Both dithionite and methyl viologen are not expected to bind the enzyme. This analysis provided an EЈ 0 value of Ϫ389 mV for the Mj-DFTR-bound FAD. This property makes the transfer of electrons from NAD(P)H to Mj-DFTR-bound FAD highly endergonic, as EЈ 0 of NAD(P) ϩ /NAD(P)H pair is Ϫ320 mV (60). Even with F 420 H 2 as the electron source, a similar but less severe limitation is expected as the EЈ 0 of F 420 /F 420 H 2 is Ϫ360 mV (36). As discussed below, it would restrict the operation of the enzyme in an environment represented by the habitat of M. jannaschii. However, the FADH 2 so generated would readily reduce the Cys-disulfide formed at the 130 CMTC 133 element that has an EЈ 0 value of Ϫ279 mV. This design is consistent with the life style of M. jannaschii as described above and discussed below (51,52,61).
Although much remains to be learned about the redox properties of Mj-DFTR, including the structural elements that determine the redox properties of protein-bound flavins, the above-described redox characteristics could be rationalized in light of the conditions prevailing in the habitat of M. jannaschii. Within the deep-sea hydrothermal vents, the partial pressure of hydrogen varies from 4 Pa to 233 kPa (52), which translated to EЈ values of Ϫ284 to Ϫ425 mV for the H ϩ /H 2 pair (calculated from the EЈ 0 value of Ϫ414 mV of H ϩ /H 2 (62)). Because the reduction of F 420 in a methanogen cell is well coordinated with the environmental partial pressure of hydrogen (63,64), EЈ of F 420 /F 420 H 2 will also range between Ϫ284 and Ϫ425 mV and allow facile reduction of Mj-DFTR-bound FAD and efficient catalysis at high partial pressures of hydrogen. These observations call for a detailed mechanistic study on this novel thioredoxin reductase.
To examine the distribution of DFTR, we performed phylogenetic analysis of DFTR homologs and identified six groups of potentially non-canonical TrxRs (Fig. 9). These proteins carried non-NTR type CXXC motifs, did not contain a fully conserved nicotinamide cofactor-binding site, and fully lacked the 2Ј-phosphate-specific HRR motif (Fig. 6). In almost all these TrxRs, the CXXC center was preceded not by a GXGXX(G/A) motif as in an NTR but by GXGXXY (Fig. 9). The presence of Tyr, a large aromatic and hydrophobic residue, instead of a smaller Gly/Ala residue, at a position intervening an NADPH binding-type site and the CXXC motif could influence either the redox potential of the catalytic Cys residues or the electron donor specificity or both. In few cases, a Thr residue appeared in place of the Tyr (data not shown), suggesting a hydroxyl group and hydrophobicity as the common denominator for a residue at this location in non-canonical TrxRs. However, these departures from the NTR structure did not suggest the nature of the electron donors used by the non-canonical enzymes. In fact, the structure-guided comparative analysis presented in Fig. 6 indicated that this aspect would likely be determined by a limited number of changes within rather well conserved active site architecture. The difficulty of identifying these signatures is underscored by the DFTR as well as the following additional examples. The non-canonical TrxR of C. pasteurianum (9), an anaerobic bacterium, that uses ferredoxin with EЈ 0 of approximately Ϫ400 mV (48), is a close relative of DFTR (Fig. 11), and it is likely that similar to DFTR the EЈ 0 of the FAD bound to this clostridial enzyme is low enough to exclude the participation of NAD(P)H as an electron source and to justify the use of ferredoxin as electron carrier. However, none of these features could have been predicted from the primary structure of this enzyme. This case is further complicated by the fact that the sequence of the CXXC motif of this enzyme is same as that of E. coli NTR (Fig. 6). An opposite situation is presented by the T. maritima TrxR, which, despite having significant structural differences with the E. coli NTR (Fig. 6), is still a nicotinamide-dependent enzyme, using both NADH and NADPH (16). Another surprise is offered by the flavin-containing TrxR of T. acidophilum, an anaerobic archaeon, where the EЈ 0 of the protein-bound FAD is Ϫ305 mV, and yet the enzyme cannot use NAD(P)H as a reductant (65). A recent report describes a TrxR from an aerobic cyanobacterium that is a close structural homolog of NTR and cannot use nicotinamides (66). The overall situation is reminiscent of the enolase superfamily, where few changes within a generally conserved active site have allowed transformation of a diversity of substrates (67). Therefore, it would be rather difficult if not impossible to identify new non-NTR flavin-dependent TrxRs and predict the natures of their electron donors solely via computational analysis. This caution aside, the report opens up the long studied flavin-dependent thioredoxin reductase family for further exploration. These enzymes belong to a broader disulfide oxidoreductase family of flavoenzymes (lipoamide dehydrogenase, glutathione reductase, trypanothione reductase, mercuric reductase, and NADH peroxidase), which are known for the diversity of their substrates (46). However, the members that have been studied thus far are nicotinamidedependent, and the possibility that some of the uncharacterized members use other types of electron carriers has not been explored.
Considering the possibility that the electron donor choices for the non-NAD(P)H-dependent enzymes could be diverse, we have proposed a naming strategy for this group where the identity of the electron donor (for example, deazaflavin or D) and the presence of flavin (F) are indicated. The four letter abbreviations, such as DFTR for deazaflavin-dependent flavincontaining thioredoxin reductase, will distinguish these enzymes from FTR (ferredoxin thioredoxin reductase), which contains iron-sulfur cluster and utilize ferredoxin (F) as electron carrier (5).
The DFTR brought a major surprise in the area of the structure-function relationships in deazaflavin-dependent enzymes. Although both F 420 and NAD ϩ (P) catalyze hydride transfer and facilitate chemically similar reactions, none of the previously described F 420 -dependent enzymes show meaningful similarities to functionally analogous or other type NAD ϩ (P)-dependent enzymes (27,(43)(44)(45). In contrast, at the primary structural level DFTR exhibits 33% identities and 50% similarities to E. coli NTR, and they share several conserved sequence features (Fig.  6). This is highly significant, as F 420 is key to biological production of methane by methanogens (36), antibiotic production by Actinobacteria (36,57), and diseases such as tuberculosis (68 -70). Thus, for the research on F 420 metabolism in these areas and others, one now needs to consider the possibility of encountering deazaflavin-dependent enzymes that are significantly homologous to cognate nicotinamide-dependent enzymes. The discovery of DFTR provided a new perspective on the evolution of NTR. The distribution of DFTR homologs led to the hypothesis that the development of NTR occurred in concert with the environmental changes on Earth. It is striking that although coenzyme F 420 is present in all characterized methanogens, DFTR homologs were limited to the most phylogenetically deeply rooted group euryarchaeal class of Methanococci (Figs. 8 and 9), which utilize primarily hydrogen as energy source (51,71), carry limited numbers (2-4) of Trx homologs, and lack NTR and FTR. The more ancient members of this group, the Methanocaldococcus species live in the early Earth environment of deep-sea hydrothermal vents (51,52). In contrast, late evolving methanogens, such as Methanobacteria and Methanomicrobia, were devoid of DFTR homologs, carried 1-9 Trxs, 1-3 NTRs, and 0 -3 FTRs (9,20). The Methanococci are more sensitive to oxygen than the late evolving methanogens (72). Thus, DFTR seems to be part of the simpler Trx systems that exist in more fastidious and phylogenetically and deeply rooted methanogens inhabiting niches presenting early Earth type reduced environment with an abundant supply of hydrogen, and a transition to the NTR and FTR systems brought about more tolerance to oxygen exposure or it is a consequence of an adaptation to niches that experience occasional oxygen exposure and/or less reduced conditions. Hence, DFTR likely represents one of the early forms of flavin-containing thioredoxin reductases from which NTR developed.
The Trx-based metabolic regulation has been studied extensively in aerobic organisms and was considered less relevant to the strictly anaerobic microorganisms (20). This is likely because the focus has been on the catalytic Cys residues (EЈ 0 , Ϫ279 mV), which are unlikely to be oxidized and consequently of less use in metabolic regulation in the anaerobic microorganisms living not only in the absence of oxygen but also in a highly reduced environment. The Trx reductase with a flavin-gating system described in this report changes this situation. By virtue of carrying a low potential flavin, the enzyme would function only under highly reduced conditions, leaving target proteins with cysteine disulfides even under moderate anaerobic conditions. This control on catalysis will allow the thioredoxin system to be responsive to redox changes within reduced environments. Because flavins provide a wider range of possibilities in terms of redox potentials (EЈ 0 , ϩ150 to Ϫ500 mV) (73), the gating system would also offer a diversity of redox potential ranges for the thioredoxin system to operate. Thus, our findings open up a new area of research on methanogens and other anaerobes, which are highly relevant in bioenergy production, bioremediation, climate change, and gut metabolism in relation to beef and dairy production and the development of type 2 diabetes and obesity in human (20).

Experimental Procedures
Materials-E. coli NiCo21(DE3) competent cells (74) were obtained from New England Biolabs (Ipswich, MA). F 420 was purified from Methanothermobacter thermoautotrophicus and reduced to generate F 420 H 2 employing previously described methods (59,75). C. pasteurianum ferredoxin and bidirectional hydrogenase were generous gifts from Dr. J. H. Chen of Virginia Tech (76). All other chemicals were purchased from standard suppliers.
Generation of Recombinant and Homogeneous Preparation of M. jannaschii Thioredoxins and DFTR-Heterologous expression and purification of Mj-Trx1 and Mj-Trx2 were performed as described previously (20). An overexpression vector for Mj-DFTR, pUL206, was constructed by amplifying the respective open reading frame (locus tag number MJ1536) from M. jannaschii chromosomal DNA and cloning it into the NdeI and BamHI sites of pTev5, a T7-based expression vector (77). It was designed to generate the recombinant protein with an NH 2 -terminal His 6 tag (MSYYHHHHHHDYDIPTSENLYFQ-GASH). The plasmid was then transformed into E. coli NiCo21(DE3). The resulting strain was grown at 37°C in 20 liters of Luria Bertani media containing 100 and 34 g/ml ampicillin and chloramphenicol, respectively. After 5 h of growth, cells were harvested at 10,000 ϫ g using a Beckman Coulter Avanti J-E centrifuge (Beckman Coulter, Brea, CA) for 20 min at 4°C and stored at Ϫ20°C until used.
DFTR purification was performed under air at 4°C as described previously (20). Briefly, 20 g of frozen cell pellets were thawed on ice and resuspended in 50 ml of lysis solution containing 50 mM sodium phosphate buffer, pH 7.0, 300 mM NaCl, 10 mM imidazole, and a cOmplete TM EDTA-free protease inhibitor mixture tablet (Roche Applied Science). The cells in the suspension were lysed by three passages through a French pressure cell operating at a pressure of 1.28 ϫ 10 8 Pa and the cell lysate was centrifuged at 20,000 ϫ g at 4°C for 30 min. From the resulting supernatant, recombinant Mj-DFTR was purified via Ni 2ϩ -nitrilotriacetic acid chromatography (20) where the protein was eluted at an imidazole concentration of 150 mM. The eluted fractions that were found to contain homogeneous DFTR via SDS-PAGE were pooled and concentrated using an Amicon ultracentrifugal filter, 3-kDa molecular mass cutoff (EMD Millipore, Bedford, MA). The NH 2 -terminal His 6 tag of the purified protein was cleaved by incubation at 4°C with recombinant His 6 -TEV (recombinant TEV-protease carrying a His 6 tag) (78), at an Mj-DFTR to protease ratio of 20:1 under overnight dialysis against a solution containing Tris-Cl buffer, pH 8.0, NaCl, and DTT at final concentrations of 50, 200, and 5 mM, respectively. DTT was then removed by three more rounds of dialysis against the same solution but without DTT and containing a higher final concentration of NaCl (400 mM). From this preparation the recombinant His 6 -TEV, Mj-DFTR carrying uncleaved His 6 tag, and the cleaved His 6 tag were removed via a passage through an Ni 2ϩ -nitrilotriacetic acid column. The protein concentrations in Mj-TrxR and Trx solutions were determined via the Bradford assay (79) using a protein assay kit (Bio-Rad) following the manufacturer's protocol.
The yield of Mj-DFTR was 0.8 mg of purified protein/liter of culture. This low protein yield was due to two factors. First, the protein was expressed in E. coli deliberately at a basal level in the absence of isopropyl 1-thio-␤-D-galactopyranoside-driven induction of the T7lac-promoter. Second, the culture was harvested at a relatively low optical density (0.8 at 600 nm, as measured with a DU-800 UV-visible spectrophotometer, Beckman Coulter, Inc., Brea, CA). These cultivation and expression conditions were chosen to increase the purity of Mj-TrxR protein preparation. The results from early trials showed that Mj-TrxR purified from E. coli cells overexpressing the protein under isopropyl 1-thio-␤-D-galactopyranoside induction contained several tightly bound proteins.
Size Exclusion Chromatography-The size exclusion chromatography was performed employing a 7.8-mm ϫ 30-cm TSK-GEL G3000SWXL column (TosoHaas, Montgomeryville, PA), a 6-mm ϫ 4-cm SWxl guard column (TosoHaas), and a Shimadzu Prominence HPLC system consisting of LC-20AD dual pumps, SIL-20A autosampler, SPD-M20A diode array detector, and CBM 20A controller system (Shimadzu Scientific Instruments, Columbia, MD) as described previously (80). An aqueous mobile phase composed of 100 mM potassium phosphate buffer, pH 7, and 150 mM NaCl was used. The flow rate was 1 ml/min. Before applying a sample, the column was equilibrated with the mobile phase for 1 h. The elution was monitored at 280 nm. A mixture of protein standards (Bio-Rad) containing the following components (protein or molecule, estimated molecular mass) was used for the calibration of the column: vitamin B 12 , 1,350 Da; myoglobin, 17,000 Da; ovalbumin, 44,000 Da; ␥-globulin, 158,000; thyroglobulin, 670,000 Da. Then, 160 g of Mj-DFTR dissolved in a 75-l mobile phase was analyzed.
Reconstitution of Purified DFTR with Flavin Adenine Dinucleotide-The purified Mj-DFTR was yellowish in color, and the UV-visible spectra of the protein exhibited absorbance maxima at 380 and 460 nm, which are typical characteristics of flavin-containing proteins (25). The ratio of the absorbance of freshly purified Mj-TrxR at 280 and 460 nm (A 280 /A 460 ) was typically 11, suggesting a partial incorporation of flavin into the protein; full incorporation of FAD in E. coli NTR leads to an A 280 /A 460 value of 4 (32). To achieve full incorporation of FAD, the purified preparation of DFTR was incubated in a solution containing 100 mM potassium phosphate buffer, pH 7, 100 mM KCl, and 1 mM FAD at 4°C for an hour. The excess FAD was removed by extensively washing the protein with the same solution but without FAD on the membrane of an Amicon ultracentrifugal filter, 3-kDa molecular mass cutoff (EMD Millipore, Bedford, MA). This final preparation exhibited an A 280 /A 460 value of 4.1.
Insulin Disulfide Reduction Assay-The insulin disulfide reduction assay was performed as described previously (20,26), except it occurred under anaerobic conditions (20). The assay temperature was 25°C. The following electron flow scheme was used: electron donor 3 Mj-TrxR 3 Mj-Trx 3 insulin. NADH, NADPH, dithionite, coenzyme M, and ferredoxin at final concentrations of 0.1, 0.1, 1, and 6 mM and 5 M, respectively, were examined as potential electron donors. The assays were performed in round glass cuvettes sealed with cutoff butyl rubber stoppers with N 2 (1.3 ϫ 10 5 Pa) in the headspace (27), except with ferredoxin as electron carrier H 2 (1.3 ϫ 10 5 Pa) was used. Ferredoxin was from C. pasteurianum, and it was reduced in situ with bidirectional hydrogenase from the same organism, and the respective final concentrations in the assay were 5 M and 50 nM, respectively (76).
Oxidation of Mj-DFTR with 2-Aldrithiol-This task was accomplished with the addition of 2-aldrithiol to a final concentration of 1 mM to a solution (total volume of 4 ml) containing 100 mM potassium phosphate, pH 7, and 50 M enzyme and incubation of the mixture at 4°C for 90 min. The excess aldrithiol was removed by three rounds of filtration in an Amicon Ultra-0.5 centrifugal filter device 3-kDa molecular mass cutoff (Millipore, Billerica, MA). Each round involved centrifugation at 10,000 ϫ g for 10 min followed by resuspension of the retentate in 0.5 ml of 100 mM MOPS buffer, pH 7.
Reduction of DFTR-bound Flavin with Various Reductants-A solution containing 30 M oxidized Mj-DFTR, 100 mM potassium phosphate buffer, pH 6.8, and one of the following reductants (final concentration) was incubated at 70°C for 30 min: NADPH (0.1 mM), dithionite (1 mM), and F 420 H 2 (0.1 mM). Then the UV-visible spectrum of the mixture was recorded using a Beckman Coulter DU800 spectrophotometer (Brea, CA). In each case, a solution containing all components but without Mj-DFTR was used as controls.
Determination of the Midpoint Redox Potential of the Catalytic Cys-Disulfide/Cys-SH Pair of the Mj-DFTR-A series of solutions containing 30 M oxidized DFTR, varying concentrations of DTT (threo-1,4-dimercapto-2,3-butanediol), and oxidized DTT (trans-1,2-dithiane-4,5-diol) giving a total concentration of 5 mM for these two compounds and 100 mM MOPS buffer, pH 7.0, were incubated inside an anaerobic chamber at 25°C for 2 h. Then to each of these mixtures mBBr, a thiolspecific labeling reagent, was added to a final concentration of 4 mM, and the incubation was continued for another 45 min but in the dark. The labeling reaction was stopped by the addition of trichloroacetic acid to a final concentration of 10% (v/v), and the mixtures were incubated on ice for 1 h. The resulting pro-tein pellet was recovered by centrifugation at 10,000 ϫ g for 10 min, washed with ice-cold acetone, and air-dried and dissolved in a 200-l solution containing 100 mM Tris-HCl buffer, pH 7, and 1% SDS. The solution was then analyzed via fluorometry by the use of an Infinite M200 plate reader (Tecan Group Ltd., Männedorf, Switzerland) with the excitation and emission wavelengths set to 380 and 450 nm, respectively. Three mixtures of the following compositions served as controls: lacking both the DTT and oxidized DTT (providing background fluorescence for subtraction from all other values); only 5 mM DTT (providing a fully reduced Mj-DFTR preparation); only 5 mM oxidized DTT or ox-DTT (providing a fully oxidized Mj-DFTR preparation). From the collected fluorescence emission intensity data, the relative amounts (%) of the DFTR molecules with fully reduced thiol groups were calculated by setting the values for the fully oxidized and fully reduced preparations to 0 and 100%, respectively, and these data were plotted against the redox potential values for the incubation mixtures that were obtained from the following equation, where R and F are universal gas and Faraday constants, respectively; T and n are temperature and number of electrons involved in the reaction, and EЈ 0 for the DTT/DTT ox is Ϫ332 mV (81).Then the data in the plot were fitted to the Nernst equation for a 2-electron transfer reaction (n ϭ 2), see Equation 3.
Determination of the EЈ 0 of the Mj-DFTR-bound Flavin Cofactor-This task involved redox titration of Mj-DFTRbound FAD with F 420 H 2 or dithionite as reductant. A series of aqueous solutions containing oxidized Mj-DFTR (25 M), F 420 H 2 (5-115 M), or dithionite (0 -500 M) and 100 mM potassium phosphate buffer, pH 7.0 were prepared from anaerobic stocks inside an anaerobic chamber (Coy Lab Products, Grass Lake, MI) and incubated there for 2 h at 25°C. The internal gas atmosphere of the chamber was composed of N 2 and H 2 (ratio 98:2, v/v); hydrogen did not have any influence on the spectral characteristic of the mixtures (data not shown). A solution lacking the reductant (F 420 H 2 or dithionite)was used as control. After incubation, each mixture was transferred into a UV-transparent cuvette with a 1-cm light path (BrandTech Scientific, Inc., Essex, CT), and the cuvette was sealed with a butyl rubber stopper (size 00 with 1 ⁄5 of the bottom cut off) and brought outside from the anaerobic chamber. Then the respective absorption spectrum (200 -800 nm) was collected on a Beckman Coulter DU800 spectrophotometer (Brea, CA) and analyzed as described under "Results." With F 420 H 2 as the reductant, the F 420 /F 420 H 2 system served as the reference for calculating the redox potential (EЈ) values for the mixtures. For the titration with dithionite, MV, a redox dye, was added as the reference, and the resultant MV ϩ /MV 0 system allowed calculations for system EЈ values.
Kinetic Assay of F 420 H 2 -dependent Trx Reductase (DFTR)-The assay was performed anaerobically with nitrogen (1.3 ϫ 10 5 Pa) in the headspace (20,82). The assay mixture (total volume, 1000 l) contained 100 mM potassium phosphate buffer, pH 7.0, 2 mM EDTA, 1 mM oxidized glutathione, F 420 H 2 , and Mj-Trx or DTNB at desired levels. The reaction was initiated by the addition of the enzyme and monitored using a Beckman Coulter DU800 spectrophotometer by following an increase in absorbance at 400 nm. For assay with DTNB, the increase in absorbance at 400 nm was due to the formation of both F 420 and TNB. Accordingly, the initial rate of TNB production was calculated using composite extinction coefficient values of 24.4 mM Ϫ1 cm Ϫ1 (⑀ 400 nm, TNB ϩ 0.5 ϫ ⑀ 400 nm, F420 , ⑀ 400 nm, TNB ϭ 11.9 mM Ϫ1 cm Ϫ1 ). A reaction without Mj-DFTR served as a control. The initial velocity was calculated using an extinction coefficient value of 25 mM Ϫ1 cm Ϫ1 for F 420 at 400 nm (83). For pH studies, the potassium phosphate buffer was replaced with constant ionic strength buffers composed of 60 mM MES, 120 mM Tris, and 60 mM glacial acetic acid and adjusted to the desired pH values (4 -9.5) with HCl or NaOH (84).
For kinetic analysis, the assay mixture was supplemented with oxidized glutathione (1 mM) to maintain a constant level of Mj-Trx (30). Each assay was performed in triplicate. The apparent kinetic constants were calculated by fitting the initial velocity data to the Henri-Michaelis-Menten's equation using Solver function in the Microsoft excel (85).
A multiple sequence alignment was performed by the use of PROMALS3D (92). Three-dimensional structures of E. coli NTR and T. acidophilum TrxR (PDB codes 1CL0 and 3CTY, respectively) served as the references.