| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
,
Department of Biochemistry & Molecular Biology, University of Georgia, Athens, Georgia 30602
Received 7 January 2007/ Accepted 5 April 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
P. furiosus utilizes a range of both simple and complex carbohydrates and converts them to acetate, to CO2, to H2, and, if S0 is present, to H2S. Its glycolytic pathway has been extensively studied and served as one of the model systems for elucidating the modified Embden-Meyerhof pathway in archaea (39, 40, 51). The key feature of this pathway is that the classical enzymes glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and phosphoglycerate kinase are replaced by a single ferredoxin-linked enzyme, glyceraldehyde-3-phosphate ferredoxin oxidoreductase (GAPOR, EC 1.2.1.-). This converts glyceraldehyde-3-phosphate to glycerate-3-phosphate without the generation of ATP (51). Consequently, only ferredoxin serves as the electron acceptor in glycolysis and no NAD(P)H is formed. The reason for this became apparent upon the discovery in P. furiosus of a membrane-bound hydrogenase (MBH) which evolved H2 from reduced ferredoxin in an energy-conserving manner via a proton motive force (41). The other oxidation step in the conversion of glucose to acetate is also coupled to the reduction of ferredoxin, and this is catalyzed by pyruvate ferredoxin oxidoreductase (POR) (4). Hence, in the production of acetate from glucose, all of the reductant is generated as reduced ferredoxin. Its oxidation via the MBH is thought to result in the conservation of energy equivalent to 1.2 ATP per glucose (41). Two cytoplasmic hydrogenases from P. furiosus have also been characterized (29), but these use NADPH, rather than ferredoxin, as the electron carrier, and their functions in fermentative metabolism are unclear. It is possible that one or both serve to recycle the H2 produced by the membrane-bound enzyme to generate NADPH for biosynthesis (29, 47).
For P. furiosus and related heterotrophic Thermococcales, S0 reduction, like H2 production, was proposed to be a mechanism for disposing of excess reductant (9, 43). H2 production is now regarded as an energy-conserving process (41), but it is not known if this is also true of S0 reduction. The S0 reduction system of the mesophilic bacterium Wolinella succinogenes is generally accepted as a model system for anaerobic S0 respiration in which H2S production is coupled to energy conservation (17). W. succinogenes uses H2 or formate as the electron donor, and their oxidation is coupled through cytochrome b and quinones to a membrane-bound, molybdopterin-containing sulfur reductase (8). A similar S0-reducing respiratory system has been characterized for other autotrophs, including the hyperthermophilic bacterium Aquifex aeolicus (13) and the hyperthermoacidophilic archaeon A. ambivalens (27), and it appears to be present in the H2-oxidizing hyperthermophilic archaea Pyrodictium brockii (36) and Pyrodictium abyssi (24).
P. furiosus and the other heterotrophic hyperthermophilic archaea seem to have a mechanism for S0 reduction that is different from that found for the autotrophic species. The available genome sequences of three Pyrococcus and one Thermococcus species (7, 11, 23, 38) do not contain obvious homologs of the molybdenum-containing sulfur reductase of W. succinogenes or A. ambivalens (2). There are also no reports of the presence of quinones or cytochromes in these organisms. Three enzymes from P. furiosus, the two cytosolic hydrogenases (29) and a sulfide dehydrogenase (30), have been previously reported to possess S0 reductase activity in vitro. However, both the activity and the expression of the two hydrogenases dramatically decreased in cells grown in the presence of S0 (1, 46). Similarly, the sulfide dehydrogenase is now thought to function in vivo as a ferredoxin:NADPH oxidoreductase (28), and the expression of its genes is related to the carbon source rather than to S0 (44). Consequently, none of these three enzymes is likely to play a role in S0 reduction in vivo.
In a previous study with P. furiosus, the effect of S0 on the expression of a selected group of genes (271) was investigated using a targeted DNA microarray by comparing cells grown for many generations (as batch cultures) in the presence or absence of S0 (46). While a significant number of genes were affected, including those encoding the two cytoplasmic hydrogenases (which were down-regulated in S0-grown cells), this study was limited by (i) the small number of genes analyzed and (ii) the use of batch-grown cells, where regulated genes involved in the primary metabolism of S0 could not be distinguished from those causing secondary or other effects. In the present work, the primary response of P. furiosus to S0 has been investigated using a kinetic approach, where S0 is added to a log-phase culture and changes in gene expression are analyzed using a complete genome DNA microarray. In addition, the enzyme responsible for NAD(P)H-linked S0 reduction was characterized independently by biochemical approaches, and the up-regulation of the gene encoding it (PF1186) was shown to be a component of the primary response of P. furiosus to S0 addition.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Gas analyses. By use of the 100-ml cultures, samples (500 µl each) were taken from the liquid and headspace, and these were injected into 10-ml anaerobic vials containing an inner reaction vial (500-µl Eppendorf tube) surrounded by 0.1 M NaOH (1 ml) to capture H2S. Sulfuric acid (100 µl, 2.0 M) was added to the inner vial to release acid-labile sulfide from the liquid culture. H2S production was measured in the NaOH phase of the double-vial system with the methylene blue assay (6). H2 gas was detected in the headspace of the vials with a gas chromatograph (Shimadzu GC-8A). The Bradford method was used to estimate protein concentration for harvested cells with bovine serum albumin as a standard (5).
Enzyme assays.
Intact P. furiosus cells were harvested from fermentor mid-log-phase cultures by centrifugation (10,000 x g, 10 min) and were gently resuspended in fresh growth medium without S0 or maltose to a final protein concentration of
15 mg/ml. The production of H2 and H2S by intact cells using maltose (50 mM) as the source of reductant was measured with the 10-ml double-vial system described above containing 0.1 M NaOH (1 ml). The inner reaction vial (500 µl) contained the reaction mixture (100 µl), including 12.8 g/liter colloidal sulfur (Fluka, Milwaukee, WI) where indicated. NAD(P)H-dependent elemental sulfur oxidoreductase (NSR) activity was measured using the same double-vial system. The standard NSR assay mixture (50 µl) contained 50 mM phosphate buffer (pH 7.0), 10 mM NAD(P)H, 6.4 g/liter (wt/vol) colloidal sulfur, and 200 µM coenzyme A (CoA). The mixture was incubated for 5 min at 85°C, the reaction was stopped, and H2S was released by the addition to the reaction mixture of 100 µl 2 M H2SO4 and quantitated as described above. One unit of NSR catalyzed the production of 1 µmol of H2S per min under these conditions. Ferredoxin-linked H2S production was measured by the same method except that NADPH was omitted and the electron donor was ferredoxin reduced by the POR of P. furiosus. The 50-µl assay mixture contained 50 mM phosphate buffer (pH 7.0), 10 mM pyruvate, 1 mM CoA, 5 mM ADP, 10 µg ml1 POR, 10 µM ferredoxin, and 6.4 g/liter colloidal sulfur. The addition of ADP allows the acetyl-CoA that is generated by the POR reaction to be utilized by acetyl-CoA synthetase (present in the cell extract), thereby preventing accumulation of acetyl-CoA. The reaction mixture was incubated at 85°C for 10 min, and H2S formation was determined as described above. The kinetic analyses were carried out under the same conditions except that NADPH, NADH, CoA, and CoA disulfide (CoA-S-S-CoA) concentrations were varied as indicated. POR and ferredoxin were purified from P. furiosus as described previously (45). Polysulfide was prepared as described previously (30), and the final concentration in the assay mixture was 11 mM. NAD(P)H, CoA, CoA-S-S-CoA, coenzyme M, glutathione, and dephospho- and desulfo-CoA were purchased from Sigma (St. Louis, MO).
RNA extraction and DNA microarray analyses. Total RNA was extracted from cell extracts of P. furious by use of acid-phenol (46) and stored at 80°C until needed. The design and construction of DNA microarrays containing all of the predicted 2,192 open reading frames (ORFs) in the annotated genome of P. furiosus (37), preparation of cDNA from the RNA samples, and hybridization experiments were all performed as previously described (46). Fluorescently labeled cDNA was prepared using an ARES DNA labeling kit (Molecular Probes, Eugene, OR). The resulting amine-modified cDNA was purified using a QIAquick PCR purification kit (QIAGEN, Valencia, CA) according to the manufacturer's instructions except that the wash buffer was replaced with 75% (vol/vol) ethanol and the cDNA was eluted with 45 µl of distilled water and dried under vacuum. The amine-modified cDNA was labeled with Alexa dye 488, 546, 594, or 647 (Molecular Probes) according to the manufacturer's instructions. The labeled cDNA was purified using a QIAquick PCR purification kit (QIAGEN) and dried under vacuum. Differentially labeled cDNAs derived from P. furiosus cells grown in the absence of S0 or from cells harvested at various times after S0 addition (up to 60 min) were pooled and hybridized to the DNA microarray, and the fluorescence intensities for the Alexa dyes were measured as described previously (46). For the microarray experiments, each log2 value represents an average of two hybridization experiments performed in triplicate using cDNA derived from two different cultures of P. furiosus. The log2 ratios were subjected to an unpaired t test function with two-tailed distribution to get the raw P values. The raw P values were subjected to a family-wise error rate correction by use of the Holm's step-down (18) procedure to give final P values. All microarray data are deposited in the NCBI GEO database (reference number GPL4688).
QPCR. RNA was isolated as described above and further purified twice using the Absolutely RNA cleanup kit (Stratagene, La Jolla, CA) with an intermediate DNase (Ambion, Austin, TX) treatment (30 min, 37°C). cDNA was then prepared as described previously (53). The genes PF1186, PF2051, PF2052, PF1441-PF1453, PF1423-PF1436, and PF2025 were selected for study, and the constitutively expressed gene encoding the POR gamma subunit (PF0971) was selected as a control. Primers for the genes were designed using the program Array Designer v.1.16 (Premier Biosoft International, Palo Alto, CA). All QPCR experiments were carried out using an Mx3000P instrument (Stratagene) using the Brilliant SYBR green QPCR master mix (Stratagene). The comparative cycle threshold method was used to analyze the resulting data (2a), which are expressed as log2 changes (n-fold).
Purification of NSR from P. furiosus. Frozen P. furiosus cell paste (100 g) grown on peptides in the presence of S0 (50) was lysed anaerobically by osmotic shock in 200 ml of 50 mM Tris-HCl (pH 8.0) under argon followed by sonication (Branson sonifier, 10 min, power setting 4). The cell extract was centrifuged at 120,000 x g for 1 h to fractionate the soluble cytoplasmic fraction from the insoluble membrane fraction. NSR was purified from the cytoplasmic fraction by anaerobic multistep chromatography using an Akta Basic (GE Healthcare, Piscataway, NJ). Unless otherwise stated, 50 mM Tris-HCl (pH 8.0) buffer was used, and all column chromatography materials were obtained from GE Healthcare. The cytoplasmic fraction (340 ml, 5,243 mg, 5,774 units) was loaded onto a DEAE-Sepharose column (150 ml) at a flow rate of 15 ml min1. The column was washed with 2 column volumes (CV) of buffer and eluted with a NaCl gradient (0 to 1.0 M) over 20 CV. NSR eluted as 200 to 290 mM NaCl was applied to the column. Fractions with the highest specific activity were pooled (330 ml, 1,046 mg, 3,815 units) and, after being diluted with an equal volume of buffer, were loaded onto a Blue Sepharose column (40 ml) at a flow rate of 15 ml min1, washed with 2 CV, and eluted with a linear gradient of NaCl (0 to 2.0 M over 20 CV). NSR eluted as a broad peak between 0.5 and 1.9 M NaCl. Active fractions were pooled (275 ml, 198 mg, 2,163 units), diluted threefold with buffer, loaded onto a hydroxyapatite (Bio-Rad) column (40 ml) at a flow rate of 10 ml min1, and washed with 10 CV of buffer containing 5 mM phosphate (pH 7.4). Some of the activity (254 ml, 33 mg, 491 units) eluted during the wash and was purified separately, while the remainder (75 ml, 39 mg, 577 units) eluted when 129 to 162 mM phosphate was applied as part of a linear gradient of phosphate (5 to 500 mM in 18 CV). The active fractions from the wash and gradient elution steps were pooled separately, and each was diluted with an equal volume of 2.0 M (NH4)2SO4 in 50 mM Tris (pH 8.0) and loaded onto a Phenyl Sepharose HP column (20 ml) at a flow rate of 10 ml min1. The column was washed with 2 CV of 1.0 M (NH4)2SO4 in 50 mM Tris (pH 8.0) and eluted with decreasing concentrations of (NH4)2SO4 (1.0 to 0 M over 20 CV). NSR activity was eluted as 575 to 465 mM (NH4)2SO4 was applied. The two NSR samples (11 ml, 3 mg, 317 units, and 12 ml, 2.3 mg, 224 units) from the Phenyl Sepharose columns were concentrated separately using a Q-Sepharose FF column (5 ml), and each was loaded onto a Bioscale Q5 column (Bio-Rad, Hercules, CA). NSR activity eluted as 130 to 150 mM NaCl was applied. The two NSR samples were indistinguishable by sodium dodecyl sulfate (SDS) gel analysis and specific activity and were combined to yield 2.0 mg with a total activity of 221 units.
Purification of recombinant NSR. To clone the gene encoding NSR (PF1186), attB PCR primers were designed based on the Gateway cloning technology (Invitrogen) with a tobacco etch virus protease cleavage site two residues upstream of the start site on the N terminus. The forward and reverse primers were CTTACAAGTTTGTACAAAAAAGCAGGCTTAGAAAACCTGTATTTTCAGGGAGGAGAAAAGAAAAAGGTAGTCATAAT and CTTACCACTTTGTACAAGAAAGCTGGGTGTCACAAAACCCTGGCGAGGAC, respectively. Pfu polymerase (Stratagene) was used to amplify the gene of interest from P. furiosus genomic DNA, and the resulting PCR product was purified using a QIAquick PCR purification kit (QIAGEN). This was cloned into the destination vector pDEST C1 containing a six-His tag according to the manufacturer's protocols (Invitrogen). This destination vector was then transformed into the expression strain of Escherichia coli BL21(DE3)pRIL (Stratagene). These cells were grown aerobically on 2XYT medium in 2.8-liter Fernbach flasks (1 liter medium) at 37°C for 6 h with shaking at 200 rpm before induction with 1 mM isopropyl-ß-D-thiogalactopyranoside (IPTG). After 16 h at 16°C, cells were harvested (10,000 x g, 20 min), resuspended in 100 ml 50 mM Tris-HCl (pH 7.4), degassed under argon, and frozen at 20°C. All purification procedures were carried out under anaerobic conditions, and all chemicals were acquired from Sigma (St. Louis, MO). Frozen cells (28 g [wet weight]) were thawed in the presence of lysozyme (200 µg ml1), DNase I (5 µg ml1), and phenylmethylsulfonyl fluoride (1 mM) and incubated with shaking at 37°C for 1 h. The cell extract was sonicated (Branson sonifier, 10 min, power setting 4), incubated at 80°C for 30 min, and then centrifuged (40,000 x g for 45 min) to remove denatured proteins. The heat-treated cytoplasmic extract (108 ml, 206 mg) containing 1,650 units of NSR activity was loaded onto a 12-ml Ni-nitrilotriacetic acid drip column (HIS-Select nickel affinity gel; Sigma) equilibrated with 50 mM Tris-HCl (pH 8.0) containing 0.5 M NaCl. NSR was eluted with 300 mM imidazole in the same buffer. The eluted His-tagged NSR protein (20 ml, 31 mg, 1,352 units) was incubated at 23°C for 3 h with AcTEV protease according to the manufacturer's protocols (Invitrogen), diluted 30-fold in 50 mM Tris-HCl (pH 8.0), and loaded onto the second Ni-nitrilotriacetic acid column. Protease-cleaved NSR without the N-terminal His tag (97 ml, 14.6 mg, 1,443 units) did not bind to the column, while residual His-tagged protein (18.6 ml, 14.9 mg, 1,044 units) was eluted with imidazole as described above. Both proteins were concentrated separately using a Q-Sepharose HP column (5 ml), where purified recombinant nontagged NSR (10 mg, 1,672 units) and His-tagged NSR (12.5 mg, 1,400 units) were obtained.
Other methods.
SDS-polyacrylamide gel electrophoresis analysis of purified NSR was performed using 4 to 20% Long Life gels (Life Therapeutics, Australia) with a Tris-HEPES buffer system. Samples were heated at 100°C for 10 min prior to loading. Gel filtration chromatography was performed using a Superdex S-200 column (320 ml; GE Healthcare) equilibrated with 50 mM Tris-HCl buffer (pH 8.0) containing 300 mM NaCl. Matrix-assisted laser desorption ionization was performed on a Bruker Autoflex (time of flight) mass spectrometer. SDS-polyacrylamide gel electrophoresis gel bands of purified NSR were excised, destained, and dehydrated with 50% acetonitrile in 50 mM NH4HCO3 and then digested with 15 µl of 10 µg/ml trypsin for 16 h. Peptides were then extracted from the gel slice by three 15-min washes (once with 50 mM NH4HCO3 and twice with 75% acetonitrile, 0.5% trifluoroacetic acid). Peptides were purified using NuTip C18 tips (Glygen Corp., Columbia, MD) and spotted (1 µl, containing
-cyano-4-hydroxycinnamic acid) directly on a matrix-assisted laser desorption ionization plate. Data analysis was performed in Protein Prospector v 3.2.1 using MS-Fit (http://prospector.ucsf.edu/).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
In contrast to ferredoxin-linked activity, NADPH-dependent S0 reductase activity was readily measured in cell extracts and in the cytosolic fraction but not in the membranes. As shown in Fig. 3, the activity was greatly stimulated by the addition of CoA to give a specific activity in the cytoplasm of 0.7 units/mg. The enzyme responsible for the activity was purified from the cytoplasm by use of anaerobic column chromatography. Only one peak of this NADPH- and CoA-dependent activity eluted from the first chromatography column, together with a minor CoA-independent peak of NADPH-linked S0-reducing activity, which was due to ferredoxin:NADPH oxidoreductase (28). The primary S0-reducing enzyme was purified by three subsequent chromatography steps to a specific activity in S0 reduction of approximately 100 units/mg, which represented a greater-than-140-fold purification compared to the cell extract. This NADPH NSR has the highest activity reported for any sulfur reductase-type enzyme system (13, 30, 33). The purified NSR preparation gave rise to a single protein band when analyzed by SDS-polyacrylamide electrophoresis (see Fig. S3 in the supplemental material), and this corresponded to a molecular weight of approximately 50,000. The apparent mass of the holoenzyme from analytical gel filtration was approximately 100 kDa, suggesting that NSR is a homodimer (data not shown). Analysis of the SDS gel band by trypsin digestion/mass spectrometry revealed that NSR corresponds to PF1186 (43% coverage), which is predicted to encode a protein of 48,720 Da, in agreement with the biochemical analyses of NSR.
|
In the NSR assays, the rate of sulfide production from colloidal sulfur was linear (up to 10 min) suggesting that this is the true substrate for the enzyme. A lag phase in sulfide production would be expected if polysulfide, which is generated by the reaction of sulfide with elemental sulfur, was the substrate for NSR. Accordingly, a less-than-twofold increase in activity was observed, both at pH 7.0 and at pH 9.0, when polysulfide (11 mM) was used as the substrate compared to when elemental sulfur (6.4 g/liter) was used. Polysulfide is stable only at pH values of
8 and readily dissociates to colloidal sulfur and sulfide at neutral pH (14). A much greater stimulation of activity would be observed if polysulfide were the preferred substrate, particularly at the higher pH. Presumably, the colloidal sulfur generated from polysulfide is a better substrate for NSR than the elemental sulfur typically added to the assay mixture (leading to an
2-fold increase in activity). The pH optimum for sulfide production (pH 6.5 [data not shown]) is also consistent with elemental sulfur rather than polysulfide being the substrate of NSR. Given the CoA dependence of the reaction, it is possible that polysulfide derivatives of CoA are generated during catalysis, and this is currently under investigation.
PF1186 is a member of a large family of flavin adenine dinucleotide (FAD)-dependent pyridine nucleotide-disulfide oxidoreductase genes (InterPro IPR013027). Accordingly, both native and recombinant forms of NSR were yellow in color and exhibited a UV/visible spectrum characteristic of flavoprotein (A460/A280 = 0.13 for the native enzyme). PF1186 and its homolog (PH0572) from P. horikoshii were previously proposed to function as NAD(P)H-dependent CoA-S-S-CoA reductase (CoADR) genes, and the recombinant form of the P. horikoshii enzyme was characterized (15). The aerobically purified P. horikoshii CoADR apoprotein was reconstituted with FAD and had a specific activity for CoA-S-S-CoA reduction of approximately 8.3 µmol/min/mg (at 75°C) using NADPH as the electron donor. The reported Km value for CoA-S-S-CoA (30 µM) is comparable to what was found (10 µM) in the present study using P. furiosus NSR in the S0 reduction assay, although the Km value for NADPH and P. horikoshii CoADR (<9 µM [15]) is 3 orders of magnitude lower than that which was determined with P. furiosus NSR (8.5 mM). The reason for this discrepancy is discussed below. It should be noted that for the P. furiosus enzyme, the specific activity for CoA-S-S-CoA reduction (6.0 µmol CoA-S-S-CoA reduced/min/mg) is about 20-fold lower than the activity that this enzyme exhibits in the S0 reduction assay. A second discrepancy is that P. horikoshii CoADR was reported to be a homotetramer (198 kDa [15]). This is in contrast to results presented here for P. furiosus NSR (92% sequence identity), which indicate that it is a homodimer. It was also reported that CoADR activity could not be purified from P. furiosus (15), but this attempt utilized cells grown in the absence of S0. Such cells would be expected to have a much lower content of the product of PF1186, according to the data presented in Fig. 2, a conclusion supported by the molecular analyses described below.
These results therefore suggest that the previously reported CoADR activity of the PF1186 homolog (in P. horikoshii) represents only a partial reaction of its true physiological function, which is now proposed to be CoA-dependent S0 reduction. Purified P. furiosus NSR did not reduce S0 to H2S in detectable amounts in the absence of CoA (or CoA-S-S-CoA), and this cofactor could not be replaced by dephospho- or desulfo-CoA, nor could it be replaced by glutathione or the methanogenic cofactor coenzyme M. The Km value for CoA (10 µM) is much lower than the intracellular concentration in P. furiosus (860 µM [20]), indicating that NSR would normally be saturated. The affinity of P. furiosus NSR for NADPH (Km, 8.5 mM) is surprisingly low, however, given the intracellular nicotinamide nucleotide concentration in P. furiosus (0.5 mM [34]). This high apparent Km value for NADPH may be related to the proposed mechanism of NADPH oxidation by NSR. P. horikoshii CoADR (15) is thought to react with two NADPH molecules, one to reduce the active site Cys (from the sulfenic acid derivative) to give the reduced enzyme ("EH2," which could react with CoA-S-S-CoA, generating CoA) and one to produce the EH2NADPH active state (in which the FAD remains oxidized). Consequently, the prior kinetic analyses (15) might have measured only the first reaction, while the S0 reduction assay reported herein measures the second reaction. This would explain the difference in the kinetic constants for NADPH between this study and the earlier one (15), the nonlinear kinetics observed for S0 reduction (see Fig. S4 in the supplemental material), and the apparent high Km value for NADPH. A detailed study of the mechanism of S0 reduction by NSR is currently under way. Suffice it to say that since NSR contains only one cysteinyl residue, it is feasible that CoA provides the active site with a second thiol group to enable the two-electron reduction of S0, where the resulting disulfide is reduced by NADPH. Consequently, the ability of NSR to reduce CoA-S-S-CoA appears to be an artifactual side reaction of the catalytic mechanism and is not thought to have any physiological relevance.
Transcriptional analyses. Growth studies showed that P. furiosus has a rapid response to the addition of S0, with H2S detected within 10 min and a pause in growth for approximately 20 min (Fig. 1). Transcriptional analyses of the complete genome (2,192 ORFs [37]) were therefore conducted on RNA extracted from cells harvested 10, 20, 30, and 60 min after S0 addition. In spite of the apparent dramatic physiological response (Fig. 1), only 19 ORFs were significantly up-regulated more than threefold (P < 0.05) within 10 min of adding S0, and a total of 34 were down-regulated. These ORFs are shown in Tables 1 and 2, respectively, together with those that are part of corresponding and potentially regulated operons. The regulation of these 34 ORFs is proposed to represent the primary response to S0. QPCR analyses were carried out with selected key ORFs to assess the validity of the DNA microarray data, and the results (Fig. 4) are discussed below. In analyzing the data, advantage was taken of the availability of the genome sequences of three other members of the Thermococcales order, Pyrococcus abyssi, P. horikoshii, and Thermococcus kodakaraensis (7, 11, 23). All three utilize S0, and it is assumed that they use the same mechanism as that of P. furiosus. Therefore, if an S0-regulated ORF in P. furiosus is not conserved in the other three species, this calls into question a direct role in primary S0 metabolism.
|
|
|
MBH and MBX are part of the NADH dehydrogenase complex I/hydrogenase family, which includes NADH:quinone oxidoreductases (NUO or complex I), F420 quinone oxidoreductases, and energy-converting hydrogenases (Ech) (10, 16). These complexes consist of a core of six homologous subunits and through evolution and recruitment of additional subunits diverged into complexes with different physiological functions (10). For example, the core enzymes of MBH and MBX (MbxH, J-N and MbhH, J-N, respectively) are supplemented with six subunits of a ubiquitous family of cation/proton antiporters (49). Previous studies have shown that MBH is an energy-conserving complex in which the oxidation of ferredoxin and the reduction of protons are coupled to the generation of a proton motive force (41). MBH lacks homologs of the NADH- and flavin-binding subunits (NuoEFG) of the typical complex I (10, 16) and couples the oxidation of ferredoxin, rather than NADH, to proton pumping and in this case H2 production. The remarkable gene conservation and sequence similarity between MBH and MBX suggests that they have very similar functions.
As illustrated in Fig. 5, MBX is proposed, like MBH, to oxidize ferredoxin and to conserve energy by pumping protons. The question then becomes, what entity does MBX reduce, thereby ultimately leading to S0 reduction? The most logical explanation is that MBX reduces NADP and NADPH is reoxidized by NSR as S0 is reduced to H2S (Fig. 5). However, we were unable to measure ferredoxin-dependent reduction of NADP or of S0 using a membrane preparation obtained from P. furiosus cells harvested from 10 to 60 min after S0 addition or using cells grown for many generations with S0. Since the amount of S0 reduced by growing P. furiosus cells approximates the amount of H2 produced on a cellular protein basis (see Fig. S2 in the supplemental material), it would seem likely that MBX and MBH conserve comparable amounts of energy during growth with and without S0, respectively, and that most if not all of the reductant used to reduce S0 flows through MBX. We speculate that our inability to measure an activity associated with MBX is due to instability of the complex. In this regard, it should be noted that MBH readily loses its ability to use ferredoxin as an electron donor during membrane fractionation (42, 47). Similarly, bacterial complex I has yet to be isolated as an intact complex (54). Attempts to stabilize P. furiosus MBX are in progress.
|
In addition to those encoding NSR, MBX, and PDO, there are only four other ORFs that are significantly up-regulated within 10 min of S0 addition (Table 1), and these appear to be present as two operons. PF2051 and PF2052 are both annotated as transcriptional regulators, but only PF2051 contains regulatory domains (IPR001845), while PF2052 contains a nucleotide triphosphate pyrophosphohydrolase domain (IPR004518). Both ORFs most likely form an operon, as their sequences overlap by 22 nucleotides, and this synteny is conserved in the genome sequences of the other three Thermococcales. It is possible that these regulators play a key role in the primary response to S0, although the nature of the effector is unknown. The other two ORFs that are regulated, PF0261 and PF0262, overlap by 12 nucleotides. PF0262 shows homology to multidrug efflux systems (IPR001036), while PF0261 has a nucleic acid binding fold (IPR008994). Homologs of both ORFs are found in the genomes of the other two Pyrococcus species but not in the genome sequence of T. kodakaraensis, so their role in S0 metabolism is unclear.
Primary response to S0 (down-regulated ORFs). The most striking feature of the list of 34 ORFs whose expression is down-regulated within 10 min of S0 addition (Table 2) is that 22 of them contain the structural genes of the three hydrogenases and another involves hydrogenase maturation. It was already known that cells grown for many generations with S0 contain only very low hydrogenase activity (1), and it is now apparent that the biosynthesis of all three hydrogenases is rapidly curtailed within minutes of S0 addition, indicating a complete shutdown at the genetic level of H2 metabolism in the presence of S0 (Table 2 and Fig. 4). Cells continue to produce H2 for 2 h or so after S0 addition (Fig. 1) due to the existing hydrogenase protein in the cell, but by 1 h the rate is <10% of the rate of H2S production, and eventually no H2 is produced. The small amount of H2 subsequently consumed (Fig. 1) may reflect a differential stability between the MBH (less stable) and the cytoplasmic enzymes, which are proposed to consume H2 and reduce NADP (29). The hydrogenase genes are well conserved within the Thermococcales order, although hydrogenase II, whose exact function is unknown, is absent from P. horikoshii and T. kodakaraensis (29).
It is difficult to rationalize the roles of the remaining 11 ORFs that are part of the primary response to S0 (Table 2), particularly when not all of them are conserved in the other three Thermococcales. PF0450 encodes a putative glutamine synthetase but it was not regulated in cells grown for many generations with S0 (46), and the fact that it is affected so quickly after S0 addition is puzzling. While PF0450 has homologs in the other sequenced Thermococcales, this is not the case for the other S0-responsive ORFs, suggesting that they may not be playing essential roles in S0 metabolism. PF0528-PF0531 appears to form an operon and is annotated as a cobalt transporter. This operon is conserved in T. kodakaraensis but not in P. horikoshii or P. abyssi. Similarly, PF1621 contains a fibronectin-like fold (IPR008957) but has a homolog only in T. kodakaraensis. Since P. horikoshii and P. abyssi do not utilize sugars like P. furiosus and T. kodakaraensis, perhaps the latter two ORF systems are involved in sugar metabolism. PF0925-PF0926 is predicted to contain a radical SAM domain, but close homologs are absent from P. horikoshii, P. abyssi, and T. kodakaraensis. PF0736 and PF0736.1 are hypothetical ORFs on opposing strands and show little sequence similarity to any protein in the NCBI database (2). PF0913 is a homolog of subunit E of formylmethanofuran dehydrogenase, an enzyme found in methanogens, but this is not the catalytic subunit and its function is unknown (52).
Secondary response to S0. The ORFs that are up-regulated only 10 min after S0 addition are assumed to represent the primary response, and all remain up-regulated at 30 min. At this time, an additional 27 ORFs are up-regulated more than threefold (or are part of a potentially regulated operon), and these appear to represent a secondary response to S0 (Table 3). This is supported by the fact that most of them (15 of 27) are involved in metabolism of glutamate or branched-chain amino acids. PF0204-PF206 and PF1852 are potentially involved in glutamate biosynthesis (21, 44) and these might compensate for down-regulation (by 4.1-fold) of the ORF that encodes glutamate dehydrogenase (data not shown), although why is not clear. Three ORFs involved in iron metabolism are also up-regulated by S0 (Table 3). Presumably, the production of intracellular sulfide (by cytoplasmic NSR) might lead to insoluble iron sulfides, and ORFs involved in ferrous iron transport (PF0857) and iron-sulfur cluster biosynthesis (PF1285, PF1286) are up-regulated in response to the products of S0 reduction.
|
In conclusion, while the addition of S0 to a culture of P. furiosus cells causes growth to stall, indicating a large metabolic shift (Fig. 1), only two key enzymes that appear to be directly involved in S0 reduction, MBX and NSR, were identified. As shown in Fig. 5, MBX and NSR are proposed to be the key enzymes responsible for the reoxidation of ferredoxin and NAD(P)H, respectively. Another primary response to S0 availability is the concomitant shutdown of H2 metabolism, resulting in the preferential transfer of reducing equivalents to the reduction of S0 rather than protons (Fig. 5). This novel S0-reducing system involving NSR and MBX is so far unique to the heterotrophic Thermococcales and contrasts with the cytochrome- and quinone-based S0-reducing system in autotrophic archaea (and bacteria). Future research will focus on elucidating the precise role of MBX and the mechanism of S0 reduction by NSR.
| ACKNOWLEDGMENTS |
|---|
We thank Peter S. Horanyi for providing the pDEST C1 Gateway expression vector, Cindy Lim for preliminary QPCR studies, Scott D. Hamilton-Brehm for assistance in constructing the microarrays, Frank E. Jenney, Jr. and Angeli L. Menon for many helpful discussions, and Farris L. Poole II for bioinformatics analyses.
| FOOTNOTES |
|---|
Published ahead of print on 20 April 2007. ![]()
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles: