INTRODUCTION
Dyes that contain an azo bond (−N=N−) contribute significantly to world dye production (
Stolz, 2001). Almost half of the dye used for dying fabrics does not bind to the fabric, and as a result, is lost in the effluent (
Shore, 1995). As these dyes can be carcinogenic (
Alves de Lima et al., 2007) their removal from the effluent is essential. This has driven investigations into both electrochemical reduction of azo dyes (
Wang et al., 2010b), as well as the use of both aerobic and anaerobic bacteria for bioremediation (
dos Santos et al., 2007;
You et al., 2007).
Azo bonds are also found in some drugs including those used for the treatment of inflammatory bowel disease (IBD). Azo compounds have also been shown to be effective antibacterial (
Farghaly and Abdalla, 2009) and antitumor (
El-Shafei et al., 2009) agents. Azo drugs for the treatment of IBD consist of a carrier compound linked via an azo bond to 5-amino salycilic acid (5-ASA). 5-ASA if administered orally is rapidly removed from the digestive system and so it is modified with an inert carrier to prolong its action (
Lichtenstein and Kamm, 2008). The azo bond is cleaved by intestinal bacteria (
Peppercorn and Goldman, 1972) releasing 5-ASA.
Enzymes able to reduce an azo bond have been found in a number of species including
Escherichia coli (
Nakanishi et al., 2001),
Pseudomonas aeruginosa (
Wang et al., 2007),
Enterococcus faecalis (
Chen et al., 2004),
Sinorhizobium meliloti (
Ye et al., 2007),
Saccharomyces cerevisiae (
Sollner et al., 2007) and humans (
Cui et al., 1995). These enzymes are flavin dependent oxidoreductases. Despite some of them sharing only minimal sequence identity, they have a common three dimensional fold (Fig. 1B) (
Li et al., 1995;
Ito et al., 2006;
Wang et al., 2007;
Ye et al., 2007;
Binter et al., 2009).
The azoreductase reaction is proposed to be catalyzed via a Bi-Bi ping-pong mechanism (
Nakanishi et al., 2001;
Liu et al., 2008;
Wang et al., 2010a), and we have indicated that the reaction is likely to occur via reduction of the hydrazone tautomer of the substrate, as identified from structural studies (Fig. 1A) (
Ryan et al., 2010). The first reduction reaction is proposed to be followed by cleavage of the hydrazine intermediate via an acid or base catalyzed reaction (Fig. 1C). This cleavage releases a quinoneimine that is believed to be the substrate for the second reduction reaction (Fig. 1D).
A number of bacterial enzymes, classified as azoreductases, can also reduce quinones (
Nakanishi et al., 2001;
Chen et al., 2004;
Liu et al., 2008;
Binter et al., 2009). The mechanism for azo reduction that has been proposed (
Ryan et al., 2010) explains the ability of azoreductases to reduce both azo and quinone substrates. It has been suggested that the major biological role of azoreductases is in the detoxification of quinones (
Liu et al., 2008,
2009). These data suggest that azoreductases should be reclassified as NAD(P)H quinone oxidoreductases (NQOs). NQOs have also been shown to be important for bacterial metabolism of some substrates, e.g.,
p-nitrophenol (
Zhang et al., 2009). Even though mammalian NQOs have very low sequence identity with the bacterial enzymes (less than 15% between paAzoR1 and either human NQO1 or NQO2), they also reduce quinone and azo compounds (
Cui et al., 1995). Recent evidence implicates the eukaryotic NQOs in the regulation of proteasomal degradation of tumor suppressor proteins (
Sollner et al., 2009b;
Alard et al., 2010). Therefore, further understanding of the relationship between azoreductases and NAD(P)H quinone oxidoreductases has wide implications in biology.
Three azoreductases have been identified in the genome of
P. aeruginosa: these are termed paAzoR1 (PA0785), paAzoR2 (PA1962) and paAzoR3 (PA3223 (
Wang et al., 2007)). The three enzymes have significantly different substrate specificities (
Ryan et al., 2010), and we describe here a detailed enzymic characterization of these proteins, including the NADH quinone oxidoreductase activity of paAzoR3.
RESULTS
Cofactor and hydride donor preference
Absorbance spectra for each purified enzyme (Fig. 2) show that all three proteins have a typical flavoprotein signature: absorbance maxima shifted 7–10 nm bathochromically compared to that of flavin mononucleotide (FMN). The spectra each have a distinct shoulder at λ
max of 486 nm as an indication of binding of the flavin cofactor to the protein (Fig. 2) (
Duurkens et al., 2007). The flavins are released as indicated by disappearance of the shoulder in the spectra upon denaturation of the protein and were identified as FMN via TLC (Supplemental Fig. 1).
paAzoR1 is unusual among azoreductases. Unlike many of these enzymes (
Nakanishi et al., 2001;
Deller et al., 2006), it can utilize both NADH and NADPH as electron donors with a preference for NADPH (
Wang et al., 2007). Interestingly, this preference was reversed on mutation of Tyr131 to phenylalanine in the enzyme’s active site (
Wang et al., 2010a). In contrast to paAzoR1, which shows a preference for NADPH as the hydride donor, both paAzoR2 and paAzoR3 show a preference for NADH (Fig. 3A). paAzoR1 shows the same hydride donor preference when reducing the quinoneimine 2,6-dichloroindophenol (DCIP) as it does when reducing the azo substrate methyl red (Fig. 3B); however, the difference in quinone reductase activity with each of the two nicotinamide cofactors is smaller than what is observed when reducing azo substrates. In contrast paAzoR3 utilizes both NADH and NADPH with a preference for NADH of less than twofold when reducing methyl red, while the rate of reduction of DCIP with NADH is almost 20 times greater than with NADPH (Fig. 3C).
The reasons behind the nicotinamide cofactor preference remain unclear. Two crystal structures have been solved showing nicotinamide cofactors bound to flavodoxin-like proteins. The first structure showed rat NQO binding NADP
+ and has been withdrawn from the PDB (
Li et al., 1995). The second structure shows NADH binding to EmoB from
Mesorhizobium sp. BNC1 (PDB: 2VZJ) (
Nissen et al., 2008). NADH binds EmoB in a solvent exposed cleft that extends from the solvent channel occupied by methyl red in the structure bound to paAzoR1. When the structure of EmoB is aligned with the paAzoR1, however, there are significant clashes between residue side chains and NADH; as a result, it is unlikely that NADH will bind in this cleft. The most possible binding position would be similar to what was modeled by
Sollner et al. (2009a), where NADH binds in the same pocket of the active site as balsalazide (Fig. 1A).
Enzyme stability
Although
P. aeruginosa lives in relatively mild conditions in soil and on human skin and in the intestinal flora, paAzoR1 has previously been shown to be thermostable with a melting temperature (T
m) of approximately 55°C (
Wang et al., 2007). Both paAzoR2 and paAzoR3 show significantly lower melting temperatures when compared to wild type paAzoR1 (40°C and 50°C, respectively, Fig. 4). Thermostability is not uncommon among azoreductases, e.g.,
B. subtilis (T
m 86.5°C) and
S. cerevisiae (T
m 60.2°C) determined via circular dichroism (
Deller et al., 2006). When Tyr131 was mutated to Phe within the active site of paAzoR1 (equivalent residues are prolines in paAzoR2 and paAzoR3), it caused a significant shift in T
m to approximately 42°C (
Wang et al., 2010a). Oligomerisation state has previously been suggested as a possible reason for high thermostability (
Deller et al., 2006), although later work on
B. subtilis azoreductase disagreed (
Binter et al., 2009). Preliminary data from analytical ultracentrifugation experiments on paAzoR2 and paAzoR3 (
Wang, 2008) indicate they form dimeric and trimeric assemblies, respectively, which supports the proposal that oligomerisation state may play a role in determining thermostability.
paAzoR1, paAzoR2 and paAzoR3 are tolerant of high salt concentrations showing greater than 40% maximal activity even at NaCl concentrations up to 2.5 M (Fig. 5). In contrast to paAzoR2 and paAzoR3 whose peak activity is at 0.1 M NaCl (Fig. 5B and 5C), peak activity for paAzoR1 is around 1 M NaCl (Fig. 5A). This high tolerance to salt may be due to stronger interactions at the predominantly hydrophobic oligmer interfaces. This high salt tolerance is advantageous for degradation of azo dyes, since salt concentrations can be as high as 2–3 M in textile effluents (
Environmental Protection Agency, 1997).
Azo substrate kinetics
The true values for the kinetic constants for all three azoreductases are shown in Table 1. One notable result is the variation in
Km values for the preferred hydride donor. In the cases of paAzoR2 and paAzoR3 the
Km for NADH is 525 μM and 298 μM, respectively; in contrast, the
Km of paAzoR1 for NADPH is over 10 mM. This suggests that the activity of the enzymes may be tied to the redox state of the cell. It has also been proposed that the
in vivo role of human NQO1 is in regulating the intracellular redox state (
Long et al., 2002).
During assays for paAzoR1 activity against methyl red, substrate inhibition was observed at methyl red concentrations of > 80 μM. The reason for this is believed to be the low affinity of paAzoR1 for NADPH, as the Km for methyl red is ~30 folds lower than for NADPH. No substrate inhibition was observed with balsalazide; thus, although the binding of methyl red is significantly tighter than for balsalazide (greater than 9:1 difference), balsalazide shows a markedly higher Vmax (more than 4:1).
Azobenzene was tested as a substrate for paAzoR1 and paAzoR3. In contrast to the other azo substrates, azobenzene has neither hydroxyl nor amine groups as substituents on the aromatic ring (Supplemental Fig. 2A); thus, it cannot form a hydrazone tautomer. Even using up to 20 μg of either paAzoR1 or paAzoR3 in the reaction, there was no reduction of azobenzene with either NADH or NADPH as the hydride donor. Azobenzene is significantly smaller than balsalazide (Fig. 1A); therefore, it can easily be accommodated by the active site. These data support the proposal that formation of a hydrazone tautomer is essential for azo compounds to be reduced by this class of enzymes.
Expression of azoreductases during growth on azo compounds
Colonies of P. aeruginosa were grown in LB media containing one of the following drugs: balsalazide, sulfasalazine, olsalazine and 5-ASA, or dyes methyl red and ponceau BS. mRNA levels were monitored for each of the azoreductases via Q-PCR (Fig. 6).
In media containing methyl red, the levels of paAzoR2 mRNA increased by a factor of approximately 50. Growth on methyl red also caused a significant (3.5 fold) increase in paAzoR3 mRNA levels. Interestingly growth on methyl red had the opposite effect on paAzoR1 whose mRNA levels were reduced by more than 50%.
The other substrate where a significant increase was observed was growth on balsalazide, where again paAzoR2 mRNA levels showed the most significant increase (2.5 fold). mRNA levels for paAzoR1 also increased but to a less extent (1.5 fold).
More significant increases in azoreductase mRNA levels have been described in
E. coli and
B. subtilis in the presence of quinones (
Töwe et al., 2007;
Liu et al., 2009), which is compatible with the proposal that quinones are the primary physiologic substrate for azoreductases.
NAD(P)H quinone oxidoreductase activity
We previously demonstrated that paAzoR1 is able to reduce the quinoneimine substrate DCIP (Supplemental Fig. 2B) (
Ryan et al., 2010). This reduction has now been studied in more detail in the cases of both paAzoR1 and paAzoR3. The ability of these enzymes to reduce the quinone substrate menadione (2-methyl-1,4-napthoquinone, Supplemental Fig. 2C) was also evaluated.
The
Km value for DCIP with paAzoR1 is similar to those obtained for azo compounds: 85.3 ± 3.6 μM for DCIP compared to 98.6 μM for balsalazide (
Wang et al., 2010a). The
Vmax value is similar to those observed for azo substrates:
Vmax = 0.53 ± 0.01 μM/s for DCIP compared to 0.81 μM/s for balsalazide (
Wang et al., 2010a). The specific activity for DCIP was however 3.5 fold higher than for balsalazide: 106 μM⋅s
-1⋅mg
-1 for DCIP compared to 29.7 μM⋅s
-1⋅mg
-1 for balsalazide (
Ryan et al., 2010).
In the case of paAzoR3, DCIP is reduced at least 30 times faster than methyl red. The maximum activity, which could be measured was at 18.75 μM DCIP, and is equivalent to a specific activity ~3220 μM⋅s
-1⋅mg
-1 (> 30 times that observed for methyl red the best azo substrate for paAzoR3) (
Ryan et al., 2010). Due to substrate inhibition of the enzyme, it was not possible to determine the
Km of DCIP (Fig. 7); however, the
Km is likely to be less than 20 μM.
Menadione was also chosen as a test substrate as it is a substrate of azoreductases from
E. coli (
Nakanishi et al., 2001),
E. faecalis (
Chen et al., 2004) and
R. sphaeroides (
Liu et al., 2008). The structure of menadione bound to human NQO2 (PDB 2QR2) (
Foster et al., 1999) is shown in Fig. 1B overlaid onto the structure of paAzoR1. The proposed mechanism for menadione reduction, based on data from azo compounds, is shown in Fig. 1E. paAzoR1 did not show activity with menadione even when up to 20 μg of enzyme was used in the assay using either NADH or NADPH as hydride donor. These data reinforce the concept that paAzoR1 has a strong preference for reduction of single ring substrates, e.g., methyl red (Supplemental Fig. 2D), rather than naphthol containing substrates, e.g., ponceau BS (Supplemental Fig. 2E) (
Ryan et al., 2010).
On the other hand, paAzoR3, which reduced naphthol containing azo substrates, showed a high affinity for menadione (Km<7.5 μM) with NADH as a proton donor. The specific activity for the reduction of menadione was calculated as 910 μM⋅s-1⋅mg-1 (ten times the specific activity observed with methyl red). When NADPH was used as the proton donor no activity was observed. This shows that paAzoR3 is an NADH quinone oxidoreductase as well as an azoreductase.
DISCUSSION
The physiologic role of bacterial azoreductases has remained an unanswered question although recent work has suggested they may be involved in the detoxification of quinones (
Liu et al., 2009). The accepted mechanism for azoreduction was unable to explain how both azo and quinone substrates are reduced by the same enzyme. However, we have resolved the apparent contradiction by proposing a novel mechanism for azoreduction (Fig. 1C and 1D) (
Ryan et al., 2010). The inability of paAzoR1 and paAzoR3 to reduce azobenzene, a substrate unable to form the hydrazone tautomer required for the proposed mechanism, lends strong support to this mechanism. The ability of both paAzoR1 and paAzoR3 to reduce DCIP, an analog of the proposed quinoneimine substrates, also supports the proposed mechanism. The high affinity of paAzoR3 for the quinone, menadione and its ability to reduce it at a higher rate than azo compounds further support its role as an NADH quinone oxidoreductase.
A number of therapeutic compounds, including antitumor drugs such as mitomycin C (
Tomasz and Palom, 1997), are quinone based agents that are activated in humans upon reduction by azoreductase-like proteins. A better understanding of the enzymes that activate these drugs in the body may allow the design of more selective compounds. There is an extensive history of using information from bacterial enzymes to understand the interactions of human enzymes with xenobiotics, e.g., cytochrome P450s (
Poulos et al., 1985) and more recently the arylamine N-acetyl transferase, enzymes (
Sinclair et al., 2000). The use of a series of bacterial enzymes has proved particularly instructive (
Lewis, 2001;
Sandy et al., 2005).
This manuscript describes the first detailed characterization of a series of azoreductases from the same organism. Due to the sequence diversity which is observed among proteins which adopt the flavodoxin-like fold there may be further azoreductase-like enzymes in
P. aeruginosa that have yet to be characterized, e.g., pa0949 from
P. aeruginosa shares only 11% sequence identity with paAzoR1, but shares approximately 40% sequence identity with tryptophan repressor binding protein A (WrbA) from
E. coli which has an azoreductase activity (
Patridge and Ferry, 2006). The number of azoreductase-like enzymes and the diversity of their substrates may suggest that they play more than one physiologic role and provides insight into the interactions of the human structural and functional homologs with xenobiotics. This is important for their interactions with not only drugs but also environmental toxins.
The sequence diversity among flavodoxin-like proteins is possible in part due to the non-sequence specific manner in which the flavin binding cradle operates (
Ryan et al., 2010). Diverse enzymes belonging to this class of enzymes have been identified within organisms but also among different kingdoms. These are connected by common activities that link proteins from archaea, eukaryotes and prokaryotes. A phylogenetic tree has been generated using the neighbor-joining method for the amino acid sequence of flavodoxin like proteins (Fig. 8) (
Armougom et al., 2006). This tree clearly shows the division of the enzymes into three groups: bacterial azoreductases, enzymes that are related to the mammalian NQO enzymes and heavy metal reducing enzymes. What is particularly interesting is that these groups cross kingdoms, especially in the case of the heavy metal enzymes, including plants, yeast, and bacterial and archaeal members. This indicates that evolution does not follow the standard species lines and may suggest convergent evolution.
MATERIALS AND METHODS
Chemicals were purchased from Sigma-Aldrich, unless otherwise stated. DNA primers were also synthesized by Sigma-Aldrich (Supplemental Table 1). The P. aeruginosa PAO1 Washington strain was kindly provided by Dr Gail Preston, Department of Plant Sciences, University of Oxford, UK.
paAzoR1, paAzoR2 and paAzoR3 were expressed and purified as described previously (
Wang et al., 2007). Enzyme characterization was carried out as detailed previously (
Wang et al., 2007). 2,6-dichlorindophenol (DCIP) reduction was measured via its absorbance at 595 nm and an extinction coefficient of 14,690 M
-1⋅cm
-1 was used. Reduction of menadione was measured, as described previously (
Nakanishi et al., 2001), via monitoring the oxidation of NAD(P)H (absorbance at 340 nm). To measure accurately rates for menadione and azobenzene reduction 0.5 μg of enzyme was used per well, 150 μM NAD(P)H and 30–7.5 μM substrate in a total volume of 100 µL of 20 mM Tris-HCl, pH 8.0, containing 100 mM NaCl per well of a 96-well plate. Thermostability assays were carried out as described previously (
Wang et al., 2010a), via measurement of enzyme activity versus methyl red after ten minutes incubation at defined temperatures.
Quantitative RT-PCR
Cultures of
P. aeruginosa PAO1 were grown in LB medium containing 1 mM azo compound or 5-ASA to an absorbance at 600 nm of 1.0. A sample (5 mL) of each culture was removed RNA was extracted using an RNeasy Mini kit (Qiagen) followed by treatment with additional RNase-free DNase (Qiagen) to remove contaminating DNA. The concentration of RNA from each sample was adjusted to the same value for cDNA synthesis using a QuantiTect Reverse Transcription kit (Qiagen). Primers used in this study (RQ-785, RQ-962, RQ-223, rpoD, forward and reverse) were specifically designed to give PCR products around 150 base pair (Supplemental Table 1). Expression levels of
paazor1,
paazor2 and
paazor3 were detected by quantitative RT-PCR analysis with a LightCycler (Roche) in comparison with expression of a
P. aeruginosa σ factor rpoD (
Savli et al., 2003). Each amplification mixture (10 μL QuantiTect 2×SYBR Green PCR master mix (Qiagen), 0.5 μg cDNA, 0.5 μM suitable primers and sterile water in a total 20 μL) was subjected to the following thermo-cycling program: one cycle of 95°C for 15 min to activate the HotStart Taq DNA polymerase; 40 cycles of denaturation (94°C for 15 s), annealing (60°C for 20 s), extension (72°C for 15 s) and data acquisition (65°C for 5 s). All reactions were performed in duplicate.
In all PCR studies, RNA concentration was determined by measuring the absorbance at 260 nm and the integrity was confirmed by 1% agarose gel electrophoresis. Absence of DNA contamination was confirmed by running duplicate samples without reverse transcriptase as a negative control.
Thin layer chromatography
Thin layer chromatography (TLC) analysis (
Sherma, 2003) was performed on a 0.1 mm thick silica gel-coated glass plate (Merck) using the mobile solvent system: tertiary amyl alcohol/formic acid/water (3:1:1,
v/v) (
Fetzner et al., 1992). Flavin co-factor was extracted from the protein component by denaturation (95°C, 10 min) and centrifugation (10,000 g, 5 min). Flavins were detected by direct visual observation of their yellow color. Authentic FMN, FAD and riboflavin were prepared in 20 mM Tris-HCl pH 8.0 buffer as standards.
Higher Education Press and Springer-Verlag Berlin Heidelberg 2010