Glyceraldehyde-3-phosphate ferredoxin oxidoreductase (GAPOR) catalyzes a dual substrate electron-transfer reaction in which oxidation of glyceraldehyde-3-phosphate (G3P) yields 3-phosphoglycerate (3PG) concomitant with reduction of oxidized ferredoxin. GAPOR and its homologs are found mainly in archaea including both euryarchaeota and crenarcheaota and especially in hyperthermophiles (
26,
38,
42). Homologs can also be found in draft sequences of deltaproteobacterial species (M. A. Aklujkar, J. E. Butler, G. Reguera, E. S. Shelobolina, A. L. Lapidus, D. Bruce, and P. M. Richardson, presented at the 106th General Meeting American Society for Microbiology, Orlando, FL, 21 to 25 May 2006). It has been proposed that the role of GAPOR is to replace the conventional two-step reduction of G3P to 3PG catalyzed by glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and 3-phosphoglyceric acid phosphokinase (PGK) in only the glycolytic direction of the Embden-Meyerhof-Parnas pathway (
26). GAPOR uses ferredoxin as an electron acceptor while GAPDH uses NAD(P)
+. A glycolytic pathway dependent on GAPOR is supposed to be a primordial pathway in early evolution (
38) since at least some ferredoxins are more stable and more thermodynamically preferable for electron transport than pyridine nucleotides at high temperatures (
36).
M. maripaludis is the sole example so far of a mesophile that contains a gene annotated as GAPOR (GAPOR
Mm) in its genome sequence (
43).
M. maripaludis is a facultative autotroph (
18) which is capable of autotrophic growth in mineral medium in the presence of hydrogen and carbon dioxide. On the other hand,
M. maripaludis readily incorporates acetate and a variety of amino acids when they are present (
44,
45). The amino acids are incorporated into protein but are not subject to further catabolism (
45).
M. maripaludis has also been reported to be capable of synthesis of intracellular glycogen reserves (
46). This was suggested from the result that activities of glycolytic enzymes were enhanced while the intracellular glycogen content increased (
46). However,
M. maripaludis has never been reported to be able to assimilate any extracellular sugars. It is therefore surprising that
M. maripaludis possesses a homolog of
P. furiosus GAPOR (GAPOR
Pf) in its genome despite limited opportunities for glycolysis. Furthermore,
M. maripaludis also possesses two homologs of putative NAD(P)-dependent glyceraldehyde-3-phosphate dehydrogenases (phosphorylating [GAPDH] and nonphosphorylating [GAPN]), which all in theory are capable of catalyzing the G3P-dependent step of the glycolytic pathway. Yu et al. (
46) earlier performed GAPDH assays using crude extracts of
M. maripaludis; however, no genome sequence was available, and GAPOR had not yet been discovered at that time.
The present study was initiated with an effort toward understanding the physiological role of a gene product which at first instance would not be expected to play an important role in an organism best known for its autotrophic growth capability. Toward this aim, we verified the catalytic function of GAPORMm through synthesis and purification of the recombinant protein in Escherichia coli and investigated its biochemical and biophysical characteristics with a particular focus on metal specificity. Subsequently, we investigated the relationship between the three enzymes (GAPDH, GAPN, and GAPOR) that are all theoretically capable of catalyzing G3P metabolism in M. maripaludis.
MATERIALS AND METHODS
Media and growth conditions of M. maripaludis.
Techniques used for growing methanogens were those of Miller et al. (
25); a serum bottle capped with a butyl rubber inner cap and aluminum sealer was used for cultivation.
M. maripaludis ATCC 43000, a strain which is not identical with the S2 strain used for genomic sequencing (
17), was grown at 37°C in MPl (DSM medium 141 [
http://www.dsmz.de/microorganisms/html/media/medium000141.html ] with omission of acetate, yeast extract, and Trypticase) which was sparged (H
2-CO
2 [ratio of 80/20], liquid for 5 min; headspace for 30 s) prior to autoclaving. While the medium was hot, 17% cysteine · H
2O · HCl (wt/vol) was added through a 0.2-μm-pore-size filter until the medium lost the color of the redox indicator. The inocula were 10% (vol/vol) of cells in the mid-exponential phase. After the cells were inoculated, the headspace was overpressurized to ∼2 atmospheres with the same gas filtered through a 0.2-μm-pore-size filter (Millex; Millipore). For GAPOR assays, cells were harvested on a 0.2-μm-pore-size filter from 5 ml of culture in an N
2-CO
2-H
2 (90/5/5 ratio) atmosphere and lysed by the addition of BugBuster Reagent (Novagen) containing nonionic surfactants under the same atmosphere, with or without the addition of 4% (vol/vol) Triton X-100. The protein concentration in the lysate was determined using a detergent-compatible protein assay (DC Protein Assay; Bio-Rad). The cell lysate was stored at −80°C under the same atmosphere until use. For transcript assays, cells were harvested by centrifugation from 1 ml of culture under the same atmosphere. Three individual cultures were treated the same way.
Transcript analysis.
Total RNA was prepared using an EZNA Fungal RNA Miniprep Kit (Omega Bio-Tek, Inc.) according to the manufacturer's instructions. For real-time PCR assays, 30 ng of total RNA was used as a template. Primers GdF (5′-GTAGCTAAACAGGACGATATGAAAG) and GdR (5′-TCAACTACAATATCTGCATCTTCAA) for GAPDH (locus MMP0325), GnF (5′-GAATCAATATCTAAAAATGCGAAAA) and GnR (5′-ATGTTTCAACAACTTCTTCTTCAAC) for GAPN (locus MMP1487), and GpF (5′-GTAGCTAAACAGGACGATATGAAAG) and GpR (5′-TCAACTACAATATCTGCATCTTCAA) for GAPOR (locus MMP0945) were designed to yield fragments of approximately equal size (200 bp). Primers 16SF (5′-ATGCGAGTCTATGGTTTCGGCCATG) and 16SR (5′-GACGCTTTAGGCCCAATAAAAGTGG) were designed to amplify 16S rRNA as an internal control. Real-time PCR was performed with an AccessQuick reverse transcription-PCR System (Promega). Following the reverse transcription reaction (at 45°C for 45 min), PCR was performed according to the following conditions: 95°C for 2 min and 40 cycles of 95°C for 30 s, 58°C for 30 s, and 68°C for 1 min. The PCR products were analyzed by agarose gel electrophoresis. Reverse transcriptase was excluded from negative controls. Relative expression levels of target genes were quantitated by the comparative cycle threshold method (
4) using 16S rRNA as an internal control.
Cloning of GAPOR.
Genomic DNA was extracted from actively growing
M. maripaludis using a DNeasy Tissue Kit (QIAGEN) after repeated cycles of freezing and thawing (
3). The entire open reading frame (accession number CAF30501; locus MMP0945) of the gene encoding GAPOR
Mm was amplified by PCR using the following oligonucleotide primers: 5′-GACGACGACAAGATGAACATTTTGATTGATGG and 5′-GAGGAGAAGCCCGGTTTATTCTTTTAATTTCCAG in the presence of 8% (vol/vol) dimethyl sulfoxide and cloned into pET46 Ek-LIC vector (Novagen) using T4 DNA polymerase (Novagen). As the reverse primer included a stop codon, the pET46-Ek/LIC-GAPOR vector was expected to generate a fusion protein encoding only an N-terminal His tag. The DNA was used to transform
E. coli strains BL21(DE3) or Rosetta-gami 2(DE3) (Novagen).
Preparation and purification of protein samples.
E. coli was cultured (37°C; aerobic) in M9 minimal medium (0.8 g/liter of NH4Cl, 0.5 g/liter of NaCl, 7.5 g/liter of Na2HPO4 · 2H2O, 3.0 g/liter of KH2PO4), to which the following separately sterilized components were added (per liter of medium): 2 ml of 1 M MgSO4 · 7H2O, 1 ml of 0.1 M CaCl2, 0.3 ml of 1 mM filter-sterilized thiamine-HCl, and 10 ml of a trace element solution containing (per liter) 1 g of FeCl3 · 6H2O, 0.18 g of ZnSO4 · 7H2O, 0.12 g of CuCl2 · 2H2O, 0.12 g of MnSO4 · H2O, and 0.18 g of CoCl2 · 6H2O. Sterilized glucose was added to a final concentration of 2 g per liter. Any of the following was added further (per liter of medium): 10 ml of 10 mM Na2WO4 · 2H2O, 10 ml of 10 mM Na2MoO4 · 2H2O, 5 ml of 10 mM Na2WO4 · 2H2O and 5 ml of 10 mM Na2MoO4 · 2H2O, or 10 ml of deionized ultrafiltered water. For M9 medium precultures, the last condition was used. A total of 40 mg/liter of l-leucine was added to all media for cultivation of Rosetta-gami 2(DE3) cells. Large-scale expression cultures were conducted using either a 10-liter fermentor (MBF-1000; Eyela, Japan) with constant aeration of 0.1 MPa (at 120 rpm and 37°C) or 1- to 5-liter shake flasks (at 200 rpm and 37°C). Expression was induced at an optical density at 600 nm of 0.5 to 0.6 or of 0.1 for the 1- or 8-liter scale culture, respectively, with isopropyl-β-d-thiogalactopyranoside (IPTG), and bacteria was grown for 12 h. Alternatively, bacteria was grown in the same culture sparged with N2 (4 ml/min) after the addition of IPTG for anaerobic expression using a Minifor system (Lambda). The bacteria were harvested by centrifugation (8,000 × g for 20 min at 4°C), washed, and resuspended in cold phosphate buffer (10 mM; pH 8.0), and then lysed with BugBuster Reagent (Novagen) containing benzonase nuclease (Novagen) and 1 mM phenylmethylsulfonyl fluoride. Anaerobically induced samples were harvested and lysed under an N2 or N2-CO2-H2 atmosphere. The lysate was centrifuged (20,000 × g for 20 min at 4°C), and the pellet was washed in buffer A (10 mM phosphate buffer [pH 8.0], 1 mM EDTA, 4% [wt/vol] Triton X-100), centrifuged (20,000 × g for 20 min at 4°C), and resuspended in buffer A containing 8 M urea. The resuspended pellet was dialyzed (for 3 h at 4°C) against 20 mM buffer B (buffer A with the addition of 10 mM urea and 100 mM NaCl) and subsequently, under oxygen-limited conditions, against buffer B with the addition of 1 mM reduced glutathione and 0.2 mM oxidized glutathione (for 3 h at 4°C). Finally, the protein preparation was centrifuged (25,000 × g for 1 h at 4°C), and the supernatant was recovered. Protein preparation using GAPORMm expressed in Rosetta-gami 2(DE3) differed from that of BL21(DE3), as the supernatant obtained after the first centrifugation following lysis contained the majority of the activity. In this case, the addition of Triton X-100 and urea was not required for solubilization of GAPORMm. GAPORMm was purified to >75% homogeneity through a HisTrap HP (GE healthcare) column and gel filtration using Superdex 200 10/300 (GE healthcare) calibrated with low-molecular-weight and high-molecular-weight gel filtration calibration kits (GE healthcare). Protein preparations were routinely monitored with an Agilent 2100 Bioanalyzer.
Cofactor analysis.
The purified protein was subjected to metal content analysis using inductively coupled plasma-mass spectrometry ([ICP-MS] SPQ-9000 with SII nanotechnology) after wet-sample digestion by 5 ml of HNO
3 and 1 ml of HClO
4 (
11). Controls were conducted using [2Fe-2S]-ferredoxin from spinach (Sigma) and bovine serum albumin (Sigma).
Pterin content was measured using 0.5 mg of GAPOR by the previously described methods (
24,
32) with fluorescence detection (lambda excitation at 363 nm; lambda emission at 442 nm). The pterin content was calculated by comparison with pterin-6-carboxylic acid. Xanthine oxidase from buttermilk (Sigma) was used as a positive control.
Acid-labile sulfide was measured by a method described elsewhere (
5,
29). The positive control was Na
2S · 9H
2O, and 1 nmol of S
2− gave an
A 670 of 0.0014.
Enzyme assay.
All assays were conducted in screw-cap cuvettes (Sigma) with butyl rubber caps (Voigt Global Distribution, LLC, Kansas City, KS) at room temperature. All assay mixtures (excluding protein) and headspace were sparged with N
2 (99.9999%) for 2 min and 15 s, respectively; protein was added under a N
2 atmosphere, and thereafter cuvettes were capped. GAPOR activity was routinely measured as previously reported by Mukund et al. (
26) by following the reduction of benzyl viologen (BV) at
A 600 (ε
BV of 7,400 M
−1 cm
−1) (
1). The reaction mixture contained 30 μM G3P, 3 mM BV, 56 μM Na
2MoO
4, 50 mM EPPS (
N-[2-hydroxyethyl]piperazine-
N′-[3-propanesulfonic acid]) buffer, pH 8.4. G3P-dependent reduction of ferredoxin was monitored by following the reduction of metronidazole (MNZ) at
A 320 (ε
MNZ of 9.3 mM
−1 cm
−1) at room temperature, as modified from the procedure of Soboh et al. (
40). The reaction mixture contained 30 μM G3P, 2 mM
Clostridium pasteurianum ferredoxin, 0.1 mM MNZ, 56 μM Na
2MoO
4, and 50 mM EPPS buffer, pH 8.4.
C. pasteurianum ferredoxin was used instead of
M. maripaludis ferredoxin, which is not available from commercial resources. The assays were started with the addition of 63 mM cysteine-HCl, which is added to ensure completely anoxic conditions, using syringes. Extracts prepared aerobically from aerobically grown cells yielded a reversibly inactivated protein preparation as GAPOR activities measured in anoxic buffer under N
2 atmosphere could be recovered after the addition of
l-cysteine. Without the reducing agent, GAPOR reacted very slowly (38-fold decrease with 100 μM G3P). The presence or absence of oxygen during expression and protein purification did not influence total yield of GAPOR activity or the requirement for
l-cysteine-dependent activation.
Cofactor inhibition was tested using 0.1 mM G3P by the addition of either AMP, NADH, NADPH (all at 10 μM), or ATP (1 and 10 μM) in the reaction mixture using GAPORMm synthesized with Rosetta-gami 2(DE3). To maintain pH above 8.0 in reaction mixtures to which ATP was added, as confirmed using a pH indicator paper (Whatman), the concentration of EPPS buffer in each reaction mixture was increased to a 100 mM final concentration, including control reactions to which no cofactor was added.
Reverse GAPOR activity, i.e., catalysis of G3P formation, was measured by following the oxidation of prereduced BV at A 600. The reaction mixture contained 100 μM 3PG, 100 μM BV, 56 μM Na2MoO4, and 50 mM EPPS buffer, pH 8.4. Assays were started by the addition of 100 μM sodium dithionite using syringes.
GAPDH activity was determined in both glycolytic and gluconeogenic directions as described previously (
47) with minor modifications. The assay mixtures contained 2 mM NADP
+, 0.3 mM G3P, 50 mM Tris-Cl (pH 8.0), with or without 5 mM K
2HPO
4 (glycolytic direction) or 0.2 mM NADPH, 30 nM 3-phospho-
d-glyceroyl phosphate, 30 nM ATP, and 50 mM Tris-Cl (pH 8.0) (gluconeogenic direction). Each assay was started by the addition of protein. The reduction of NADP
+ and oxidation of NADPH were monitored by fluorescence spectrometry. GAPN activity was determined as for GAPDH activity above, except with the omission of ATP.
Negative control reaction assays were conducted using reaction mixtures in which either protein, G3P, or electron acceptor was omitted. No activity was observed in any control reaction. Omission of IPTG during expression of recombinant GAPORMm resulted in complete loss of measurable GAPOR activity in crude lysates. One unit of activity is defined as 1 μmol of G3P oxidized per min.
Flux balance analysis.
Flux balance analysis was performed using the condensed stoichiometric model of
M. maripaludis and
Desulfovibrio vulgaris syntrophic coculture developed by Stolyar et al. (
41). The input files of the published model were imported without modification and analyzed using the software CellNetAnalyzer (
21), Matlab (version 7.0), and the Matlab Optimization module (Mathworks). The following constraints were fixed prior to flux optimization: lactate uptake (D86), 10; ethanol synthesis (D78), 0; glycerol synthesis (D83), 0; SO
4 uptake (D01), 0; alanine uptake (M79), 0. If GAPDH (M32) is constrained to zero prior to flux optimization, all values reported in Fig.
2 of Stolyar et al. (
41) are replicated with the exception of reaction M16. In a modified model, changes were incorporated in the reactions catalyzed by GAPDH (M32, 1 ZM_GAP + 1 ZM_ADP = 1 ZM_3PG + 1 ZM_ATP + 1 ZM_NADH; change to 1 ZM_GAP + 1 ZM_ADP = 1 ZM_3PG + 1 ZM_ATP + 1 ZM_NADPH) and GAPOR (M33, 1 ZM_GAP + 1 ZM_ADP = 1 ZM_3PG + 1 ZM_ATP + 1 ZM_Fdx(red); change to 1 ZM_GAP = 1 ZM_3PG + 1 ZM_Fdx(red), an irreversible reaction). To simulate the effect of variation in H
2 availability, reaction M15 (uptake of H
2 by
M. maripaludis) was constrained in a range between 0 to 16.2 μmol/h prior to flux optimization.
DISCUSSION
GAPOR
Pf contains W (
26), and it appears specific for W because even in medium containing excess Mo or vanadium (V) and trace amounts of W, no GAPOR
Pf activity is detected (
27). The factors which affect the choice of metals used by enzymes containing molybdenum cofactor (Mo-Co), a cofactor consisting of (i) organic molybdopterin and (ii) inorganic molybdenum or tungsten, remains a question of great interest (
20). Kletzin and Adams suggested that W is obligatory in enzyme reactions which require very low redox potentials (≤−420 mV), such as the ferredoxin-dependent oxidoreductases AOR, GAPOR, and formaldehyde ferredoxin oxidoreductase (
22). Analysis of recombinant GAPOR
Mm expressed in minimal medium clearly indicated that GAPOR containing Mo and pterin is catalytically active and that the reaction catalyzed by GAPOR is not dependent on W. The accumulation of active GAPOR
Mm in
M. maripaludis cultured in chemically defined medium containing Mo but not W further supports the idea that GAPOR
Mm containing Mo and not W is catalytically active. This verifies that GAPOR
Mm encodes an enzyme with said activity and that the enzyme is actively synthesized under autotrophic growth conditions. The synthesis of active GAPOR containing only Mo suggests that the lack of GAPOR activity in extracts of
P. furiosus cultured in the absence of W and presence of Mo or V (
27) does not result from a chemical incompatibility related to the redox potential of the substrates and end products (
22).
The addition of tungstate to the medium did not result in the formation of active Mo-Co containing recombinant GAPOR
Mm, even in the absence of added Mo. This is unlikely to be due to a lack of capability by
E. coli to synthesize W-containing Mo-Co. It has previously been demonstrated that
E. coli can both incorporate tungstate using the same mechanism for molybdate uptake (
33) and synthesize Mo-Co containing both W and Mo (
7,
33). In addition, Mo-Co biosynthesis is an ancient, ubiquitous, and highly conserved pathway (
31), and
E. coli is no exception (
35). Therefore, that active GAPOR
Mm containing a tungsten cofactor was not synthesized is most likely explained, at least in this case, by the metal cofactor specificity of GAPOR
Mm. An incompatibility between
E. coli and GAPOR
Mm with respect to maturation factor requirements or the presence or absence of a nucleotide, or a type of nucleotide, would be expected to influence equally the formation of both molybdenum- and tungsten-containing Mo-Co.
Mo-GAPOR is the first Mo-containing AOR member for which the corresponding gene sequence has been verified. Since the [4Fe-4S] cluster and molybdopterin ligands identified in sequences encoding other W-containing AORs (
9) are conserved, the root of metal specificity must be determined by other as yet unidentified residues, if not dictated by the organism and/or environments as a whole.
The simultaneous addition of an equal concentration of Mo and W resulted in completely inactive MoW-GAPOR, suggesting that the presence of equimolar amounts of W inhibits the formation of Mo-containing Mo-Co, at least in the special case of GAPOR
Mm. However, this should not pose limitations on the synthesis of Mo-dependent enzymes in natural environments since Mo is far more abundant than W in most geographical locations (
22). GAPOR
Mm therefore most likely evolved to utilize Mo in accordance with the metal abundance of the habitat of
M. maripaludis, thereby differentiating from “ancestral” hyperthermophilic GAPOR at least with respect to metal specificity.
It has been suggested that GAPOR operates solely in the glycolytic direction (
26) while GAPDH and PGK support a gluconeogenic pathway (
26,
37) in sugar-utilizing
P. furiosus. Mukund and Adams (
26) also suggested that
P. furiosus produces 2 mol less ATP per mol of glucose and disposes of excess electrons as molecular hydrogen as a result of exploiting GAPOR instead of GAPDH and PGK (
26). This apparently inefficient mechanism may represent an ancestral form of glycolysis at the early stages of the evolutionary process of thermophilic organisms and can be explained by the fact that certain ferredoxins may be more stable than NAD(P)H at high temperatures (
38). However, it is puzzling why GAPOR, which catalyzes the reaction in only the glycolytic direction, is found in a facultative autotroph that is not known to assimilate sugars. Furthermore, as the genome of
M. maripaludis also contains two more candidates potentially capable of catalyzing G3P metabolism (GAPDH and GAPN), it is difficult to understand why ATP generation would be sacrificed if the organism is genetically capable of replacing the GAPOR reaction with ATP-generating GAPDH. We therefore performed a time course study to verify the presence of each of the three genes and encoded enzymes and to verify the directionality of the reactions that they catalyze (Table
2). Transcript and enzyme assays showed that only GAPOR
Mm operates in the glycolytic direction, while only GAPDH operates in the gluconeogenic direction under autotrophic conditions, as the reactions catalyzed by GAPOR
Mm and GAPDH were unidirectional with opposing directionality. This contrasts to earlier results (
46) in which phosphorylating NADP
+-dependent G3P oxidation was reported using crude extracts of
M. maripaludis cultured under conditions similar to those used in the present study. Notably, the glycolytic GAPDH activity reported earlier was >90 and >2,000 times lower than the gluconeogenic GAPDH and glycolytic GAPOR activities, respectively, that were observed in the present study. While catalytically irreversible GAPDH isoenzymes that oxidize G3P (
14) or catalytically reversible GAPDH isoenzymes with proposed physiological unidirectional roles (
12) have previously been characterized, to our knowledge there are no previous reports of any phosphorylating GAPDH which catalyzes only formation, but not oxidation, of G3P. Further study is required to elucidate the cause of such apparent unidirectionality.
No sign of GAPN-encoding transcript or GAPN activity was observed at any time point of the time course study. It was recently demonstrated that a gene annotated as GAPN in
Methanocaldococcus jannaschii that shares 56% amino acid identity with
M. maripaludis GAPN (locus tag MMP1487) in fact encodes an aldehyde dehydrogenase with no specificity for G3P (
15). The three residues proposed to be essential for binding of the phosphate moiety of G3P are also missing in the
M. maripaludis GAPN sequence. Even if the
M. maripaludis gene annotated as GAPN were expressed, it is likely that the encoded protein in any case would not catalyze nonphosphorylating NAD(P)
+-dependent G3P oxidation.
Flux balance analysis, using a model corrected with respect to the stoichiometry and directionality of the reactions catalyzed by GAPOR and GAPDH, predicted that glycolytic flux would be strictly limited, at least under steady-state conditions of autotrophic growth, in order to achieve theoretically optimum biomass formation. As only GAPOR
Mm is capable of catalyzing G3P catabolism and assuming that no other pathways of glycogen catabolism exist in
M. maripaludis, this suggests that glycogen catabolism may also be limited under non-steady-state autotrophic growth conditions. Furthermore, although the GAPDH-GAPOR pair, together with other enzymes catalyzing reversible or a combination of irreversible steps of glycolysis/gluconeogenesis, would allow a cycle of glycogen storage and turnover as previously suggested (
46), such a cycle of glycogen turnover, or a cycle between pyruvate and G3P alone, that is dependent upon both GAPDH and GAPOR in concert would constitute, in effect, an ATP-dependent NADPH:ferredoxin-oxidoreductase reaction. This raises the questions of how such a reaction cycle and/or glycogen turnover would be regulated and what would be the role of the metabolic reactions catalyzed by GAPDH-GAPOR in
M. maripaludis.
Both the transcript and enzyme analysis suggested that GAPOR was constitutive throughout the growth period, while GAPDH activity was observed mainly in the early growth stages of the cultures. This contrasts with previous results with
P. furiosus in which GAPOR
Pf was suggested to be transcriptionally controlled and induced by the presence of extracellular sugar (
41). As methanogens in general do not assimilate extracellular sugar, a similar regulation mechanism is not possible with
M. maripaludis, which may explain why GAPOR activity was constitutive. Regulation may instead be posttranscriptional. Both GAPOR
Pf (
23) and GAPOR
Mm (Fig.
2) are strongly inhibited by G3P; however, no other apparent allosteric interactions were previously noted with GAPOR
Pf (
23). It is difficult to imagine how the apparently irreversible substrate inhibition exhibited by GAPOR
Mm alone would be able to exert any control over glycolytic flux. We therefore tested whether any cofactors would be able to affect the GAPOR-catalyzed reaction. Interestingly, ATP completely inhibited the reaction even at concentrations as low as 1 μM. As this concentration is lower than intracellular levels typically reported in other microorganisms (
30), it is likely that an additional, as yet unknown, positive effector also influences GAPOR activity in vivo. Nevertheless, it is clear that, in contrast to GAPOR
Pf, GAPOR
Mm is posttranscriptionally regulated.
The results of flux balance analysis and assays determining the molecular regulation of GAPOR
Mm and enzyme directionality of GAPOR and GAPDH together allow us to formulate a tentative hypothesis regarding the metabolic role of GAPOR
Mm, which is summarized in Fig.
3B. Given that only GAPOR
Mm is able to catalyze G3P catabolism and that the enzyme is subject to posttranscriptional regulation, it would appear that GAPOR
Mm plays a central role in regulation of glycogen catabolism. The suggested roles of glycogen metabolism in microorganisms vary depending on the species, including temporary conditions of stress imposed by limitations in growth substrate (i.e., starvation) or chemical warfare or undefined continuous cycling (
2). As flux balance analysis indicated no role for GAPOR
Mm under steady-state autotrophic growth conditions, arguing against a role in continuous cycling, GAPOR
Mm and glycogen catabolism are therefore more likely to play a metabolic role under non-steady-state growth conditions (for example, stationary phase or starvation) or under growth conditions dependent on alternative substrates. The fact that ATP is a potent negative regulator of GAPOR
Mm further supports the notion that glycogen degradation in
M. maripaludis may be limited to conditions of starvation, as the concentration of ATP would be expected to be lower under such conditions. Alternatively, the role of ATP as a negative regulator of GAPOR
Mm may also be important in order to limit continuous ATP-driven and NADPH-dependent reduction of ferredoxin, a reaction that should be necessary only in the special case where both glycogen and H
2, the major source of electrons for reduction of ferredoxin (
41), are limited. In order to test this, we simulated the effect of reduced H
2 availability on predicted reaction rates using flux balance analysis (Fig.
3A). At least stoichiometrically, a cyclic GAPOR-GAPDH reaction was predicted to maximize biomass formation in the event that the uptake of only H
2 and no other external metabolites is limited. Whether such a putative reaction cycle takes place in vivo or not will require further investigation. Nevertheless, it can be concluded that both the proposed role of GAPOR and GAPDH in glycogen catabolism, as well as a putative ATP- and NADPH-dependent reduction of ferredoxin, most likely are important only for maintaining metabolic homeostasis under nonoptimal growth conditions.