Skip to main content

Assessment of metabolic flux distribution in the thermophilic hydrogen producer Caloramator celer as affected by external pH and hydrogen partial pressure



Caloramator celer is a strict anaerobic, alkalitolerant, thermophilic bacterium capable of converting glucose to hydrogen (H2), carbon dioxide, acetate, ethanol and formate by a mixed acid fermentation. Depending on the growth conditions C. celer can produce H2 at high yields. For a biotechnological exploitation of this bacterium for H2 production it is crucial to understand the factors that regulate carbon and electron fluxes and therefore the final distribution of metabolites to channel the metabolic flux towards the desired product.


Combining experimental results from batch fermentations with genome analysis, reconstruction of central carbon metabolism and metabolic flux analysis (MFA), this study shed light on glucose catabolism of the thermophilic alkalitolerant bacterium C. celer. Two innate factors pertaining to culture conditions have been identified to significantly affect the metabolic flux distribution: culture pH and partial pressures of H2 (P H2). Overall, at alkaline to neutral pH the rate of biomass synthesis was maximized, whereas at acidic pH the lower growth rate and the less efficient biomass formation are accompanied with more efficient energy recovery from the substrate indicating high cell maintenance possibly to sustain intracellular pH homeostasis. Higher H2 yields were associated with fermentation at acidic pH as a consequence of the lower synthesis of other reduced by-products such as formate and ethanol. In contrast, P H2 did not affect the growth of C. celer on glucose. At high P H2 the cellular redox state was balanced by rerouting the flow of carbon and electrons to ethanol and formate production allowing unaltered glycolytic flux and growth rate, but resulting in a decreased H2 synthesis.


C. celer possesses a flexible fermentative metabolism that allows redistribution of fluxes at key metabolic nodes to simultaneously control redox state and efficiently harvest energy from substrate even under unfavorable conditions (i.e. low pH and high P H2). With the H2 production in mind, acidic pH and low P H2 should be preferred for a high yield-oriented process, while a high productivity-oriented process can be achieved at alkaline pH and high P H2.


The need to circumvent environmental and social issues concerning the depletion of fossil fuels and the greenhouse gas emissions has led to exploration of alternative sources of energy. Hydrogen (H2) is considered as a promising energy carrier for the future because of its high energy content and non-polluting properties [1]. Moreover, it can be employed as a non-fuel commodity in a variety of industrial chemo-physical processes for which the demand is increasing [2].

Dark fermentation is a potential carbon neutral process for production of H2 from organic substrates by mesophilic or thermophilic anaerobic microorganisms. In order to establish an economically viable biological process for H2 production the yield needs to be maximized [3]. Mesophilic microorganisms are not capable of producing H2 at high yield being it reported in the range of 1-2 mol H2/mol hexose [4]. On the other hand, thermophiles have shown great potential for H2 generation, mainly because of the more favorable thermodynamics of the reaction at elevated temperatures which allows to generate a limited variety of by-products and to reach the theoretical yield of 4 mol H2/mol hexose [5].

Caloramator celer, formerly known as Thermobrachium celere [6], is a strict anaerobic, alkalitolerant, thermophilic bacterium capable of converting C6 sugars to H2, CO2, acetate, ethanol and formate as major metabolites. Previous studies have shown that C. celer produces H2 at high yields both in a naturally occurring microbial community and in pure culture [79]. However, the distribution of the end metabolites can vary depending on the growth conditions allowing the accumulation of significant amounts of ethanol and formate with consequent reduction of H2 yield.

Two innate factors pertaining to culture conditions have been identified to have a prominent role in the distribution of metabolic fluxes during anaerobic fermentation in C. celer: partial pressures of H2 (P H2) and culture pH [8, 9]. The inhibitory effect on microbial growth and H2 production caused by H2 accumulation both in liquid and gaseous phases during the fermentation is a well known phenomenon and it has been demonstrated for several H2-producing organisms [913]. The fermentation medium can be easily supersaturated with H2 due to liquid-to-gas mass transfer limitations and this seems to be inevitable also at elevated temperatures [14, 15]. In these conditions H2 synthesis becomes thermodynamically unfavorable and consequently the disposal of accumulated reducing equivalents in the cell is mediated by a metabolic shift towards production of more reduced metabolites, such as lactate, ethanol, acetone, butanol, or alanine [5]. As a consequence of the changes in the fermentation profile both H2 and ATP yields decrease. In previous studies C. celer achieved high H2 yields ( ≥ 3.3 mol H2/mol glucose) at low culture-to-headspace volume ratio and when H2 was periodically removed from the headspace [8, 9]. These observations suggest that in order to maximize the H2 yield in C. celer P H2 in the fermentation vessel should be maintained at low level.

Culture pH is another factor affecting the distribution of metabolic fluxes during glucose fermentation. However, there exists a certain disagreement on the optimal pH to employ for H2 production [16]. Certainly, the optimal value needs to be studied case by case as it depends on the metabolic and physiological properties of the microorganism under investigation. Moreover, the ideal pH for cell growth may not be the same as the one for maximal H2 production [17].

Recently, a draft whole-genome sequence of C. celer has been obtained allowing to gain insight into the metabolic potential of this organism [18]. In particular, genomic analysis revealed pathways involved in pyruvate catabolism and end-product synthesis. To get a better understanding of the fermentative metabolism of C. celer under controlled conditions, a combination of experimental data from batch cultures in continuous stirred-tank reactor (CSTR), recently available genomic data and methods of metabolic flux analysis (MFA) were employed to analyze flux distribution towards end-products, taking in special consideration the effect of two innate factors, pH and P H2, on hydrogen production.


Effect of culture pH on fermentative metabolism of C. celer

The effect of culture pH on growth and product formation in C. celer was investigated for the first time under pH-controlled conditions. Four different pH levels (8, 7, 6, 5.5) were tested during fermentation in a CSTR. At pH 8 and 7, C. celer showed about 3-fold higher maximum growth rates (μ max) than those obtained at pH 6 and 5.5 (Table 1). Similarly, higher glucose consumption rates were observed at pH 8 and 7 compared to pH 6 and 5.5. The highest biomass formation and biomass yield were similar at pH 8 and 7 and slightly lower at pH 6, while at pH 5.5 a significant reduction was observed (Table 1). Again, cultures at pH 8 and 7 displayed higher yield of biomass per ATP (YX/ATP) than the one calculated at acidic pHs, whereas the ATP yield (YATP/S) increased at acidic conditions (Table 1).

Table 1 pH-controlled batch fermentations of C. celer on 5 g/l of glucose at different culture pH

Besides growth, culture pH affected also the metabolite profile of C. celer during glucose fermentation. The highest H2 accumulation (106 mmol H2/l) and H2 yield (2.97 mol H2/mol glucose) were observed at pH 6, while at higher pH their values were about 33% lower (Table 1). In contrast, higher volumetric H2 productivity (Q H2) was obtained at alkaline conditions (25.2 mmol H2/l/h) and gradually decreased as the culture pH was reduced. The lower H2 production observed at pH 8 and 7 was accompanied by an increase of formate and ethanol yields, both being about 2-fold higher compared to those observed at pH 6 (Table 1). The lower RVEP (total molar reduction values of reduced end-products) [19] and the acetate-to-ethanol ratio calculated for fermentation at pH 8 and 7 reflected the increased production of formate and ethanol (Table 1).

Combined effect of culture pH and P H2 on fermentative metabolism of C. celer

Combined effect of pH and P H2 on growth and metabolite production profiles

To investigate the combined effect of pH and P H2 on growth, product formation and carbon flux distribution in C. celer, four different experimental conditions were tested using pH and application of N2 sparging as variables (Case I-IV). Sparging the reactor with N2 did not significantly affect the growth of C. celer at pH 7 and 6 (Figure 1A). In fact, despite the remarkable difference in the P H2max in the reactor headspace with and without N2 sparging, similar growth rates and biomass yields were achieved within the same pH level (Table 2). On the other hand, at same sparging conditions the pH had a higher impact on growth with about 2.5-fold reduction of μ max (Table 2) and a slower overall glucose consumption rate observed at pH 6 (Figure 1B). Also, YX/ATP was lower at pH 6, especially under sparging conditions (Table 2). Nevertheless, despite a slower growth at pH 6, the highest biomass formation was comparable to the one reached at pH 7 (about 1 gCDW/l) (Figure 1A).

Figure 1
figure 1

Growth and product formation by C. celer in pH-controlled batch fermentations performed at pH 7 and 6, with and without N 2 sparging. Growth (A), glucose consumption (B) and end product accumulation (C-F) profiles of C. celer cultured with 10 g/l glucose at pH 7 (straight lines) with N2 sparging (Case I, filled squares) and without N2 sparging (Case II, filled circles), and at pH 6 (dashed lines) with N2 sparging (Case III, open squares) and without N2 sparging (Case IV, open circles). Larger symbols in panel A represent sampling points for ATP, NADH and NAD+ measurements. Data are from one representative fermentation per condition.

Table 2 pH-controlled batch fermentations of C. celer on 10 g/l of glucose at pH 7 and 6, with and without N 2 sparging

Both culture pH and P H2 influenced the carbon flux distribution and the end-product profile in C. celer during glucose fermentation. The beneficial effect of H2 removal from the reactor on H2 production was more evident at pH 6 (Case III and IV). In fact, at this pH, when N2 sparging was applied, H2 accumulation increased by about 26% and H2 yield by 19% (Case III) (Figure 1C; Table 2). On the other hand, at pH 7 (Case I and II) the effect of P H2 on H2 production was marginal with a 9% and 13% increase at low P H2 for H2 accumulation and H2 yield, respectively (Case I). In addition, applying N2 sparging did not affect Q H2 within the same pH level. At the same sparging condition acidic pH promoted higher H2 accumulation and H2 yield, while neutral pH favored higher Q H2 (Figure 1C; Table 2). In fact, when N2 sparging was applied (Case I and III) H2 accumulation increased by 32% and H2 yield by 33% as pH set-point was lowered from 7 to 6. Even so, the Q H2 at pH 7 was about 36% higher than at pH 6. A similar trend could be observed without N2 sparging (Case II and IV), although the differences between H2 accumulation and H2 yield at pH 6 and pH 7 were significant but less pronounced.

Carbon and electron flow at the pyruvate and acetyl-CoA nodes was affected by both pH and P H2. While acetate accumulation profile and yield were minimally influenced by the P H2, formate and ethanol synthesis significantly increased at both pH levels when H2 concentration was allowed to build up in the reactor (Case II and IV) (Figure 1D, 1E, 1F; Table 2). Specifically, formate accumulation and yield were about 2-fold higher without N2 sparging regardless of the culture pH (Case II and IV), whereas the effect of P H2 on ethanol production was more drastic at pH 6 (Case IV) with an increase of accumulation and yield of 2.5- and 4-fold, respectively. The ability of C. celer to produce reduced end-products (RVEP), i.e. H2 and ethanol, remained unchanged under different P H2 (Table 2). On the other hand, acetate-to-ethanol ratio showed that at the acetyl-CoA node the redirection of carbon and electron flow toward ethanol synthesis is less favorable at low P H2, especially at pH 6 (Case III) where the highest acetate-to-ethanol ratio of 14.2 was observed (Table 2).

Acidic pH (Case III and IV) promoted a slight increase of acetate accumulation and yield, whereas formate accumulation and yield were reduced by more than 2-fold at pH 6 both with and without gas stripping (Case III and IV) (Figure 1D, 1E; Table 2). A drastic reduction of ethanol synthesis occurred only at acidic pH with N2 sparging (Case III), while at high P H2 both ethanol accumulation and yield remained at high level regardless of the culture pH (Case II and IV) (Figure 1F; Table 2). The increase of RVEP from 2.8-2.9 at pH 7 (Case I and II) to 3.3-3.4 at pH 6 (Case III and IV) indicates that pH influenced the conversion of glucose to reduced end-products (H2 and ethanol). The flux at the acetyl-CoA node was affected by pH only when sparging was applied resulting in a 3-fold increase of the acetate-to-ethanol ratio at pH 6 (Case III), while no variation in the ratio was observed when H2 concentration was allowed to build up in the reactor.

Combined effect of pH and P H2 on ATP and redox levels

To evaluate the effect of pH and P H2 on energy and redox metabolism of C. celer, intracellular ATP, NADH and NAD+ were measured in the exponential and stationary phase of growth (Figure 1A). At the same culture pH ATP levels were not significantly different both in exponential and stationary phase regardless of the application of N2 sparging, suggesting that intracellular ATP was not affected by the P H2 (Table 3). In contrast, culture pH showed a more significant influence on the ATP levels. In the exponential phase ATP was at least 2-fold higher at pH 7 (Case I and II) compared to pH 6 (Case III and IV). In the stationary phase intracellular ATP at pH 7 doubled (Case I and II), while at pH 6 was reduced by one-third (Case III and IV). As a consequence the difference between the ATP levels observed at pH 7 and 6 increased from 2-fold in exponential phase to about 7- to 8-fold in the stationary phase (Table 3).

Table 3 Intracellular ATP, NADH and NAD + concentrations, and NADH/NAD + ratio

In the exponential phase the NADH/NAD+ ratio measured with and without N2 sparging resulted to be similar at both pH levels (Table 3). However, at pH 7 the NADH/NAD+ ratio increased by 2-fold (from 0.06 to 0.13) during the stationary phase when P H2 reached about 56 kPa (Case II), while it remained unaffected when P H2 was kept low by N2 sparging (Case I). Surprisingly, a dramatic increase of the NADH/NAD+ ratio was observed during the stationary phase at pH 6 (Case III and IV). This was caused by an expected depletion of NAD+ pool rather than an increase in intracellular NADH concentration (Table 3).

Metabolic flux analysis under different pH and P H2

To investigate more in detail the effect of pH and N2 sparging on the fermentative metabolism of C. celer, metabolic flux analysis (MFA) was carried out for Case I to IV at different growth phases (Figure 2). A node analysis was performed to estimate the distribution of the fluxes at key metabolic branch points (Figure 3). In the proposed metabolic network of C. celer (Figure 2A; Additional file 1: Table S1), phosphoenolpyruvate (PEP) is the first important intermediate being potentially used by three reactions (v1, v6, v11) (Figure 3A). Generally, most of the PEP was converted to pyruvate by reaction v1 and v6 especially as the fermentation proceeded. In the exponential phase only a fraction of the PEP was converted to oxaloacetate (v11) and almost none in the stationary phase. The distribution at this metabolic node during exponential phase was minimally affected by the difference in P H2 at both pH. On the other hand, the conversion of PEP to oxaloacetate (v11) was affected by pH. In fact, at pH 7 (Case I and II) the PEP directed to the “malate shunt” (v11-v13) in the exponential phase was almost 2-fold higher than at pH 6 (Case III and IV). Given the higher growth rate observed at pH 7 (Figure 1A, Table 2), a higher flux through the “malate shunt” was not surprising since it served as the only source of NADPH for biomass synthesis.

Figure 2
figure 2

Metabolic flux distribution of C. celer during glucose fermentation at different stages of growth. Reconstructed metabolic network of central carbon metabolism (A) and heatmap of metabolic fluxes normalized with respect to glucose uptake of 100 mmol (B). Stoichiometry for each reaction is defined in Table S1 (Additional file 1). Abbreviations: 3PG, 3-Phosphoglycerate; AcCoA, Acetyl-CoA; ACE, Acetate; ADP, Adenosine diphosphate; ATP, Adenosine triphosphate; BIO, Biomass; E4P, Erythrose-4-phosphate; ETH, Ethanol; F6P, Fructose-6-phosphate; Fdred, Reduced ferredoxin; FOR, Formate; G3P, Glyceraldehyde-3-phosphate; G6P, Glucose-6-phosphate; GLC, Glucose; GLN, Glutamine; GLU, Glutamate; H2, Hydrogen; MAL, Malate; NADH, Reduced nicotinamide adenine dinucleotide; OAA, Oxaloacetate; Pi, Orthophosphate; PPi, Pyrophosphate; PEP, Phosphoenolpyruvate; PYR, Pyruvate; Ri5P, Ribose-5-phosphate; S7P, Sedoheptulose-7-phosphate; X5P, Xylulose-5-phosphate.

Figure 3
figure 3

Node analysis of metabolic fluxes in C. celer at different stages of growth. Fluxes at the phosphoenolpyruvate node (A), pyruvate node (B), acetyl-CoA node (C) and H2 flux (D) at pH 7 with N2 sparging (Case I, black) and without N2 sparging (Case II, dark grey), and at pH 6 with N2 sparging (Case III, light grey) and without N2 sparging (Case IV, white). The values represent the percentage of the flux considering 100 as the incoming flux at the node.

The distribution of the flux at the pyruvate node determines the size of the reduced ferredoxin (Fdred) available for H2 production since only conversion of pyruvate to acetyl-CoA by pyruvate:ferredoxin oxidoreductase (PFOR) generates reducing equivalents. Both pH and P H2, as well as the growth phase, influenced the fluxes at this node (Figure 3B). When H2 was stripped from the reactor, the v15/v14 ratio in the exponential phase almost doubled both at pH 7 (Case I) and pH 6 (Case III) suggesting that lower P H2 directs the flux through PFOR. A more dramatic effect on flux distribution at this node was exerted by the culture pH. In the exponential phase neutral pH favored acetyl-CoA formation through pyruvate formate lyase (PFL) (v15/v14 < 1 in Case I and II), whereas at acidic pH pyruvate was mainly directed through PFOR (v15/v14 > 1 in Case III and IV). Although decreased in the later stage of growth, a significant formate flux was measured at pH 7 when sparging was not applied (Case II), whereas it almost ceased in all other conditions (Case I, III and IV). In general, a direct correlation was observed between formate flux (v14) and the growth rate (Figure 4A). Routing the carbon and electron flow through PFL has negative implication for H2 production since this reaction does not supply Fdred. This is confirmed by the negative correlation existing between the formate flux (v14) and the Fd-dependent H2 flux (v19) estimated at different growth phases (Figure 5A). In addition, the importance of the flux distribution at the pyruvate node for H2 production in C. celer is also highlighted by the correlation between the fraction of the flux through PFOR in the exponential phase and the overall H2 yield (Figure 4B).

Figure 4
figure 4

Correlation between distribution of fluxes at the pyruvate node and growth rate, H 2 yield and pyruvate flux. Correlation between formate flux (v14) and growth rate (A) in exponential phase (circles), transition phase (diamonds) and stationary phase (squares) at pH 7 with N2 sparging (Case I, black) and without N2 sparging (Case II, dark grey), and at pH 6 with N2 sparging (Case III, light grey) and without N2 sparging (Case IV, white). Correlation between fraction of flux through PFOR during exponential phase and pyruvate flux during exponential phase (open symbols), and overall H2 yield (filled symbols) (B) at pH 7 with N2 sparging (Case I, diamonds) and without N2 sparging (Case II, triangles), and at pH 6 with N2 sparging (Case III, circles) and without N2 sparging (Case IV, squares).

Figure 5
figure 5

Correlation between formate and ethanol fluxes, and H 2 fluxes. Correlation between formate flux (v14) and Fd-dependent H2 flux (v19) (A), between ethanol flux (v17) and NADH-dependent H2 flux (v18) (B), and between the sum of formate and ethanol fluxes (v14+v17) and overall H2 flux (v18+v19) (C) in exponential phase (circles), transition phase (diamonds) and stationary phase (squares) at pH 7 with N2 sparging (Case I, black) and without N2 sparging (Case II, dark grey), and at pH 6 with N2 sparging (Case III, light grey) and without N2 sparging (Case IV, white).

Acetyl-CoA is another critical intermediate in C. celer metabolism. Conversion of acetyl-CoA to acetate (v16) generates one extra ATP, while ethanol production is used as an electron sink for the oxidation of NADH (v17). Although production of acetate was generally favored, ethanol flux varied depending of the experimental conditions (Figure 3C). Applying N2 sparging to prevent H2 building up in the reactor promoted a lower ethanol flux in the exponential phase (Case I and III). In particular, the flux directed to ethanol was 3- and 2-fold lower for Case I and III respectively, and similarly the v16/v17 ratio was 3- and 2.5-fold higher suggesting that low P H2 favored acetate over ethanol production in the exponential phase. Ethanol flux in the exponential phase was lower at pH 7 compared to pH 6 at both sparging conditions. However, while at pH 6 (Case III and IV) the ethanol production completely ceased in the stationary phase, at pH 7 (Case I and II) it started increasing at the end of the exponential phase resulting in about 25% of the flux to be directed to ethanol regardless of the P H2 in the reactor. Since conversion of acetyl-CoA to ethanol requires NADH as reduced cofactor, it competes with NADH-dependent H2 production for the use of reducing equivalents. Figure 5B supports this scenario showing that in C. celer NADH-dependent H2 flux (v18) estimated at different growth phases was negatively correlated with ethanol flux (v17).

Data from MFA was also used to analyze the H2 flux and how its two components, the NADH-dependent (v18) and the Fd-dependent (v19) H2 production, contribute to it (Figure 3D). In the exponential phase N2 sparging had a great effect on the NADH-dependent H2 flux (v19) that increased from 23.4 (Case II) to 66.5% (Case I) of the total H2 flux at pH 7 and from 20.9 (Case IV) to 40.2% (Case III) at pH 6. However, no evident trend can be identified in the distribution of H2 flux with respect of the culture pH in this growth phase. In general, with the exception of Case I, during the exponential phase Fd-dependent proton reduction (v19) was the favorite reaction for H2 synthesis. As the fermentation entered in the stationary phase, at pH 7 (Case I and II) Fd-dependent H2 flux (v19) was almost twice the NADH-dependent H2 flux (v18), while at pH 6 (Case III and IV) the H2 flux was equally distributed between the Fd-dependent (v19) and the NADH-dependent reaction (v18). When comparing the overall H2 flux (v18 + v19) normalized by glucose uptake (v1) it is clear that low P H2 achieved by N2 sparging (Case I and III) promoted higher H2 production in C. celer, whereas when H2 was allowed to accumulate in the reactor (Case II and IV) H2 evolution suffered (Figure 3D). Similarly, running the fermentation at acidic pH (Case III and IV) instead of at neutral pH (Case I and II) increased the disposal of reducing equivalent through H2 production, especially in the stationary phase when fluxes of the competing reactions (v14 and v17) were minimal (Figure 5C).


Generation and consumption of reducing equivalents in C. celer

During glucose catabolism reducing equivalents are generated in the form of NADH or Fdred. These cofactors need to be reoxidized allowing the glycolytic flux to proceed. During fermentation this can be achieved through different reactions involving the synthesis of reduced molecules such as H2, ethanol, lactate, butyrate and alanine. Specifically, C. celer can produce NADH during the oxidation of glyceraldehyde-3-phosphate by an NAD+-dependent glyceraldehyde-3-phosphate dehydrogenase (GAPDH, TCEL_00702) and Fdred in the oxidation of pyruvate to acetyl-CoA by pyruvate:ferredoxin oxidoreductases (PFOR, TCEL_02202-02206/TCEL_01566) (Figure 2A; Additional file 1: Table S1). While in C. celer glycolysis (Embden-Meyerhof pathway) yields 2 NADH per one glucose metabolized, generation of Fdred is dependent on flux distribution at the pyruvate node. At this node pyruvate is non-oxidatively dissimilated to acetyl-CoA and formate by pyruvate formate lyase (PFL, TCEL_00503) or oxidatively decarboxylated to acetyl-CoA by PFOR. Only the latter generates Fdred and when 100% of the flux goes through PFOR, generation of reducing equivalents is maximized resulting in the formation of 2 NADH and 4 Fdred. If all reduced electron carriers were recycled by hydrogenases, the complete oxidation of one molecule of glucose would yield 4 H2 molecules. The genome of C. celer contains two clusters coding for putative [FeFe]-hydrogenases (TCEL_00581-00584/TCEL_01273-01277) and one cluster coding for a putative [NiFe]-hydrogenase (TCEL_00187-00205) (Figure 2A; Additional file 1: Table S1). Based on similarity of the gene organization within the clusters and the high homology with the cytosolic enzyme complex of Caldanaerobacter subterraneus subsp. tengcongensis [12], the two [FeFe]-hydrogenases most likely utilize NADH as an electron donor for proton reduction. The [NiFe]-hydrogenase is predicted to be a multimeric Fd-dependent membrane-bound enzyme similar to the MBH complex characterized in Pyrococcus furiosus [20]. In this species MBH couples the oxidation of Fdred with energy conservation by forming an ion gradient that can be used to generate ATP via a membrane-bound ATP synthase.

However, recycling reducing equivalent through H2 production is not always thermodynamically favorable. At a given temperature the oxidation of reduced cofactors (particularly of NADH) by hydrogenases is a function of the H2 concentration [5]. Supersaturation of the aqueous phase with H2 can easily occur due to limitation of liquid-to-gas mass transfer rate [14, 15]. As a consequence, at high H2 concentrations the disposal of accumulated reducing equivalents in the cell is mediated by a metabolic shift towards production of more reduced metabolites such as ethanol, lactate, and alanine. The genome analysis revealed that in C. celer the only alternative for NADH oxidation is ethanol production, since no lactate dehydrogenase nor NAD:ferredoxin oxidoreductase (NFOR) were identified (Figure 2A; Additional file 1: Table S1). In addition, although C. celer encodes a complete pathway for butyrate synthesis that would consume NADH, no butyrate was detected as end-product in this study. Ethanol production mediated by putative bifunctional acetaldehyde/alcohol dehydrogenase (TCEL_01373) and alcohol dehydrogenase (TCEL_00064) has been shown to increase at high P H2 validating the hypothesis that this pathway serves as an alternative to hydrogenase for NADH reoxidation [8, 9].

As in C. celer the only possibility for Fdred oxidation relies on the activity of the ferredoxin-dependent MBH complex, the distribution of carbon and electron flow at the pyruvate nodes dictates the size of the Fdred pool available for H2 production. Despite the similar standard Gibbs energy (ΔGo´) [19], the reaction catalyzed by PFOR is expected to be less favorable at high H2 concentration due to the increased Gibbs energy (ΔG′) of the ferredoxin oxidation by H2 production, whereas the thermodynamics of the reaction catalyzed by PFL does not change. Therefore, in C. celer the branched pyruvate node can serve as a safety valve relieving the cell from the burden of ferredoxin reoxidation in unfavorable conditions and thus avoiding a decrease in the metabolic flux.

Effect of culture pH on growth and fermentative metabolism of C. celer

In previous studies performed under non-controlled conditions the optimal initial pH for H2 accumulation in C. celer was observed to be 8.2 at 67°C (approximately 9.0 at room temperature) [8, 9] in accordance with the optimal pH for growth [21]. However, in non-pH-controlled conditions the dynamic pH profile caused by the production of organic acids masked the real effect of pH on the metabolism of C. celer. In the current study, pH-controlled fermentations using both 5 and 10 g/l of glucose showed that alkaline to neutral pH (8 and 7 at room temperature) clearly favor high growth rates in C. celer confirming the alkalitolerant nature of this organism. On the other hand, at moderately acidic pH (6 and 5.5 at room temperature) the μ max was reduced by 2 to 3-fold (Table 1, 2) regardless of the P H2 in the system. However, with 10 g/l of glucose in the medium, growth ceased before complete glucose depletion and the biomass yield was identical for all the conditions (Table 2). This could be caused by nutrient limitation or by inhibiting conditions such as increased osmolality and accumulation of by-products in the culture [13, 22]. When C. celer grew at high rates in alkaline/neutral conditions, biomass formation was more energy-efficient despite the lower overall ATP recovery from glucose (Table 1, 3). In contrast, at acidic pH the growth rates were reduced and inefficient biomass formation was observed despite the more effective energy recovery. Given the alkalitolerant nature of C. celer, it is likely that at suboptimal extracellular pH this organism needs to maintain the intracellular pH within the physiological range. This goal can be achieved by H+-ATPase which exports protons at the expense of ATP hydrolysis. The activity of this enzyme is known to increase at low intracellular pH, thus consuming a substantial portion of the intracellular ATP produced via substrate level phosphorylation [23, 24]. Based on the intracellular ATP levels measured at different culture pH conditions (Table 3) and the flux of ATP hydrolysis (v21) (Figure 2B), it can be hypothesized that in C. celer at lower pH non-growth-associated ATP consumption increases most probably to achieve cytoplasmic pH homeostasis.

An analysis of the end-product distribution clearly indicated that C. celer exhibited different metabolic patterns depending on the culture pH. Surprisingly, at moderately acidic pH conversion of glucose to H2 was more favorable ( > 2 mol H2/mol glucose) than at alkaline to neutral pHs ( < 2 mol H2/mol glucose) as a consequence of the lower synthesis of other reduced by-products, i.e. formate and ethanol (Tables 1, 2). On the other hand, volumetric H2 productivity was directly correlated with μ max, both increasing as pH was shifted to alkaline pH. Notably, a transition in the metabolic behavior of C. celer was observed between pH 7 and 6 (Tables 1, 2). In particular, both the formate yield and the formate flux during exponential phase were more than 2-fold higher at pH 7 regardless of the P H2 in the reactor (Table 2; Figures 2B, 3B). A widespread tendency among PFL-encoding organisms to accumulate formate during glucose fermentation as culture pH becomes more neutral/alkaline was observed, albeit with different magnitude, both in pure [2532] and mixed cultures [33, 34]. In absence of a formate hydrogen lyase (FHL) in C. celer and other strictly anaerobic bacteria the synthesis of formate by PFL competes with generation of H2 [19].

Although at this stage the exact mechanism of this metabolic shift triggered by culture pH in C. celer is unknown, several hypotheses can be proposed: i) a simultaneous anabolic and catabolic role of PFL, as reported in closely related clostridia [35], would justify the higher flux through PFL observed at elevated growth rates in this (Figure 4A) and previous studies [9] possibly to meet energy and anabolic demands at optimal pH for growth; ii) since pyruvate can trigger its own utilization by PFL participating as allosteric effector in the activation of this enzyme [36, 37] and acetyl-CoA can inhibit PFOR [3840], at elevated growth rates observed at neutral/alkaline pH higher pyruvate and acetyl-CoA fluxes favored the reaction catalyzed by PFL, whereas pyruvate was mainly directed to PFOR under reduced pyruvate and acetyl-CoA fluxes observed at lower growth rates in acidic conditions and in the transition from the exponential to the stationary phase (Figures 2B, 3B, 4B); iii) the activity of PFL is optimal at pH slightly above 7, but reduced under acidic conditions [37, 41]; iv) while the ΔGo´ of the reaction catalyzed by PFL is not affected by pH, the reoxidation of Fdred via proton reduction necessary to drive the oxidative pyruvate dissimilation by PFOR, becomes less favorable as the pH increases [27]. Overall, given the direct correlation between flux distribution at the pyruvate node during exponential phase and H2 yield (Figure 4B), carbon and electron flow should be channeled through PFOR instead of PFL to maximize H2 synthesis.

Unlike formate, ethanol production does not seem to be dependent on pH culture. In fact, at 5 g/l of glucose with N2 sparging the highest ethanol yield (0.48 mol/mol glucose) was achieved at pH 7 and the lowest (0.20 mol/mol glucose) at pH 6 (Table 1). However, ethanol yields at pH 8 and 5.5 were similar. Moreover, at high P H2 no difference was observed between pH 7 and 6 (Table 2). This suggests that regulation of carbon flux through the acetyl-CoA branch point could potentially serve as a means of controlling the disposal of reducing equivalents necessary to maintain an internal redox balance.

Effect of P H2 on growth and fermentative metabolism of C. celer

The results obtained in this study show that under controlled conditions P H2 does not affect the growth of C. celer. Both biomass formation and bioenergetics parameters (YX/ATP and YATP/S) were the same regardless the P H2 (Figure 1A; Table 2). Also growth rate and consequently volumetric H2 productivity were not dependent on P H2 (Table 2). Under given conditions, the growth rate is dependent on the glycolytic flux which, in turn, is a function of the activity of GAPDH. The activity of GAPDH is negatively affected by an increase of the NADH/NAD+ ratio [42, 43]. During batch cultivations the NADH/NAD+ in C. celer was low at exponential phase at both pH (Table 3), suggesting a complete NADH reoxidation regardless of the P H2. As a consequence, the flux through GAPDH (v4) was not inhibited by high P H2 (Figure 2B). Unlike other anaerobic bacteria which employ alternative pathways to H2 production for the reoxidation of NADH (e.g. ethanol and lactate production) only in response to high NADH/NAD+ [43, 44], in C. celer the ethanol production started already in the exponential phase (Figure 1F) preventing the NADH/NAD+ to increase and inhibit GAPDH. As a result, in the conditions at which NADH-dependent H2 formation is supposed to be inhibited by high P H2 (Case II and IV) ethanol yield is higher (Table 2). NADH was maintained at low levels also in the stationary phase (Table 3). At this stage reoxidation of NADH through ethanol synthesis followed a pH-dependent pattern as at pH 7 ethanol flux (v17) and accumulation were still significant, while at pH 6 almost ceased (Figures 1F, 3C).

The pyruvate node is another important metabolic point for the management of reducing equivalents. Indeed, at high P H2 the overall formate yield and its flux (v14) in the exponential phase were higher suggesting that the non-oxidative reaction catalyzed by PFL serves as an alternative route for pyruvate dissimilation at high H2 concentration. In fact, while acetyl-CoA production by PFOR produces one pair of reducing equivalents (Fdred), in the PFL reaction the reducing equivalents remain with the product formate. Thus, regulation of carbon flux at the pyruvate node could potentially serve as an effective means of controlling the disposal of reducing equivalents necessary to maintain an internal redox balance. Additionally, this electron rerouting was pH-dependent since in Case II (pH 7 without N2 sparging) formate was accumulated also in the stationary phase and 25.9% of the flux at the pyruvate node was directed through PFL, whereas in all the other conditions formate accumulation and flux almost ceased at this stage of growth. Increased formate yield in response to high P H2 was reported also for Clostridium thermocellum [45, 46]. Formate production has never been previously proposed as a pathway to control redox balance in the cell, probably due to the scarce occurrence of a similarly branched pyruvate node in other strict anaerobes.

The majority of the flux being directed to acetate production in all conditions (Figures 2B, 3C) seems to suggest that C. celer aims for maximization of ATP production to sustain high growth rate at alkaline/neutral pH or mechanisms for intracellular pH homeostasis at acidic pH. According to the metabolic reconstruction (Figure 2A; Additional file 1: Table S1), four combinations of pathways for glucose catabolism are possible in C. celer (Eq. (1-4)):

Glucose + 4 H 2 O 2 Acetate + 2 HCO 3 + 4 H + + 4 H 2 4 ATP
Glucose + 3 H 2 O 1 Acetate + 1 Ethanol + 2 HCO 3 + 3 H + + 2 H 2 3 ATP
Glucose + 2 H 2 O 2 Acetate + 2 Formate + 4 H + + 2 H 2 4 ATP
Glucose + 1 H 2 O 1 Acetate + 1 Ethanol + 2 Formate + 3 H + 3 ATP

Eqs. (1-3) direct 100% of the carbon flux through PFOR (v15), whereas Eqs. (3-4) direct the flux entirely through PFL (v14). Conversion of glucose to acetate by Eq. (1) allows for maximal ATP generation (4 ATP) by completely relying on hydrogenase activity to regenerate reduced cofactors (NADH and Fdred). However, at high H2 concentration hydrogenases are inhibited and less-energy efficient pathways for disposal of reducing equivalents need to be activated to sustain the catabolic flux. Indeed, at 70°C and 1 M of dissolved H2 this reaction is the least thermodynamically favorable (ΔG′70°C = -153.1 kJ), but becomes more favorable as H2 concentration decreases (Additional file 2: Figure S1). The metabolic setup of C. celer allows for a reaction yielding 3 ATP (Eq. (4)) whose ΔG′70°C is completely independent of the H2 concentration (Additional file 2: Figure S1). Conversely, reactions thermodynamically independent of the H2 concentration used by other fermentative bacteria for balancing high NADH/NAD+ (e.g. ethanol and lactate synthesis) yield only 2 ATP [12, 43, 44]. Thus Eq. (4) ensures C. celer a higher energy-recovery from glucose breakdown and unaltered glycolytic flux even under unfavorable conditions for H2 production. In addition, C. celer can potentially still obtain 4 ATP through Eq. (3) which is slightly more favorable at high H2 concentrations than Eq. (1) (Additional file 2: Figure S1).

As shown by the inverse correlation between ethanol and formate fluxes with the two components of H2 flux (Figure 5A, B, C), ethanol and formate production serves as an alternative to H2 production for maintaining the redox balance in the cell when hydrogenases are inhibited: ADH assists NADH-dependent hydrogenase to keep low NADH levels, while PFL simultaneously supplies acetyl-CoA and stores reducing equivalents in formate thus regulating the pool of Fdred available for Fd-dependent hydrogenase. Given the stoichiometry of the fermentations and the distribution of the fluxes, C. celer utilizes a combination of the putative pathways to balance anabolic and catabolic requirements as well as intracellular redox state depending on the conditions. One exception can, however, be observed in Case III (pH 6 with N2 sparging) when C. celer entered the stationary phase. Under these conditions an equal contribution by the two components of the H2 flux and a normalized H2 flux close to the theoretical stoichiometric coefficient of 4 were observed (Figure 3D). Since no ethanol and formate production as well as biomass synthesis were detected, the fermentation proceeded according to Eq. (1). As previously observed [9], in C. celer the complete reoxidation of both Fdred and NADH through H2 production occurred at low pH and P H2 when glucose was used for non-growth associated cell maintenance.


Combining experimental results from batch fermentations, with genome analysis, reconstruction of central carbon metabolism and metabolic flux analysis (MFA), this study shed light on glucose catabolism of the thermophilic alkalitolerant bacterium C. celer. This organism possesses a flexible fermentative metabolism that allows efficient energy harvesting from substrate even under unfavorable conditions (i.e. low pH and high P H2). Two innate factors pertaining to culture conditions have been identified to significantly affect the metabolic flux distribution: culture pH and P H2.

Overall, at alkaline to neutral pH the rate of biomass synthesis of C. celer is maximized and the flux at the pyruvate node mainly directed to PFL suggesting a higher activity of this enzyme as well as a possible role in anabolic metabolism. On the other hand, at acidic pH the lower growth rate and the less efficient biomass formation are accompanied with a more efficient energy recovery from the substrate indicating a high cell maintenance requirement possibly to sustain intracellular pH homeostasis. Higher H2 yields were associated with fermentation at acidic pH as a consequence of lower synthesis of other reduced by-products such as formate and ethanol. In contrast, P H2 does not affect growth of C. celer on glucose. At high P H2 the cellular redox state was balanced by rerouting the flow of carbon and electrons to ethanol and formate production allowing unaltered glycolytic flux and growth rate, but resulting in a decreased H2 synthesis. Overall, with the optimization of H2 production in mind, C. celer offers the flexibility of shifting from a high yield-oriented process (at acidic pH with N2 sparging) to a high productivity-oriented process (at alkaline pH without N2 sparging). In particular, the tolerance exhibited by C. celer towards H2 build-up in the reactor is attractive because H2 production should be achieved preferably without the need for sparging gas to reduce capital costs for the gas-upgrading process [47]. The metabolic control and regulation of the pyruvate node as well as the genetic, metabolic and physiological traits that allow C. celer to withstand high P H2 merit further studies.

Materials and methods

Microorganism and medium

Caloramator celer strain JW/YL-NZ35, former Thermobrachium celere (equivalent to DSMZ 8682 and ATCC 700318), was obtained from the Deutsche Sammlung von Mikroorganismen und Zellkulturen (Braunschweig, Germany). C. celer was cultivated in a modified ATCC 2072 medium.

( containing (per liter): KH2PO4 0.75 g; Na2HPO4 · 2H2O 1.53 g; KCl 1 g; (NH4)2SO4 0.5 g; NH4Cl 0.5 g; MgCl2 · 6H2O 0.1 g; CaCl2 · 6H2O 0.11 g; FeSO4 · 7H2O 0.2 g; cystein-HCl 0.2 g; resazurin 0.001 g; trace element solution 10 ml; vitamin solution 10 ml; yeast extract 2 g; tryptone 2 g; glucose 5 or 10 g. Stock solutions under nitrogen atmosphere containing glucose, cystein-HCl, MgCl2 · 6H2O, CaCl2 · 6H2O, FeSO4 · 7H20 and vitamins were sterilized separately and added anaerobically to autoclaved medium at the required concentrations. Routine subcultures and inoculum development were conducted in 250 ml serum bottles containing 50 ml of medium at 67°C.

Growth conditions

Cultures were grown in a jacketed, 3-liter bioreactor equipped with an ADI 1025 Bio-Console and an ADI 1010 Bio-Controller (Applikon, Schiedam, The Netherlands) at a working volume of 1 l. The initial pH was maintained at 67°C by automatic titration with 4 M NaOH. The temperature was thermostatically kept at 67 ± 1°C and the stirring rate was set to 250 rpm. A condenser with 5°C cooling water was fitted to the bioreactor's headplate. Prior to inoculation, the medium was sparged with N2 and supplemented with an anoxic solution of cysteine-HCl at a final concentration of 0.2 g/l to make the medium completely anaerobic along with the other stock solutions. Glucose was used as a primary substrate in all experiments at an initial concentration of 5 or 10 g/l. The medium was inoculated with 50 ml (5% v/v) of a culture in the exponential phase.

The effect of culture pH on the fermentative metabolism of C. celer was studied at four different pH values (8, 7, 6 and 5.5) measured at room temperature. In this experiment the bioreactor was constantly sparged with N2 at 100 ml/min. Glucose was added as substrate in concentration of 5 g/l. The effect of P H2 at two pH levels (7 and 6) was investigated by applying two experimental conditions: continuous N2 sparging at 100 ml/min for continuous removal of produced H2 (Case I at pH 7 and Case III at pH 6) and in absence of N2 sparging with the bioreactor's gas outlet open to allow accumulation of H2 in the headspace at atmospheric pressure (Case II at pH 7 and Case IV at pH 6). In the latter experimental condition the overall gas production was obtained by measuring the volume of displaced NaHCO3-saturated solution in a graduated cylinder. In this experiment glucose was added as substrate in concentration of 10 g/l. Gas samples from the headspace for H2 and CO2 determination and culture samples for monitoring growth, substrate consumption and metabolite formation were regularly withdrawn during growth. Samples for the measurement of intracellular ATP, NADH and NAD+ concentrations were collected as described earlier [48] during the exponential and stationary phase of the cultivations supplemented with 10 g/l glucose. All the experimental conditions were performed in biological duplicates and included negative controls omitting glucose.

Analytical methods

Headspace samples were analyzed for H2 and CO2 concentration by gas chromatography, using a dual channel CP-4900 Micro-GC (Varian, Middelburg, The Netherlands) as described earlier [49]. The results were analyzed with a Galaxie Chromatography Workstation (v. Cell concentrations were determined by measuring the absorbance at 620 nm using a spectrophotometer (Ultrospec 2100 pro, Amesham Biosciences, UK). Cell dry weight (CDW) was determined as described earlier [44]. Glucose, acetate, formate, ethanol and butyrate were analyzed by HPLC (Waters, Milford, USA) on an Aminex HPX-87H ion exchange column (Bio-Rad, Hercules, USA) at 45°C, with 5 mM H2SO4 (0.6 ml min-1) as the mobile phase. The column was equipped with a refractive index detector (RID-6A, Shimadzu, Japan).

NAD(H) assay

Cell culture samples (1 ml) were collected in the exponential and stationary phase of growth for NADH and NAD+ determination. Samples were immediately quenched by transferring them to a microcentrifuge tube containing 1 ml of ice and centrifuged for 1 min at 12,100 g. The pellets were immediately frozen and stored at -80°C until further analysis. Intracellular concentrations of NADH and NAD+ were determined by a cyclic assay as described earlier [48]. Intracellular levels of NADH and NAD+ were expressed per 1 g of CDW.

Measurement of ATP

Cell culture samples (1 ml) were collected in the exponential and stationary phase of growth for ATP determinations. Samples were collected in screw-cap microcentrifuge tubes containing 1 ml of ice-cold chloroform and immediately frozen into liquid nitrogen. Samples were stored at -80°C until further analysis. Intracellular concentrations of ATP were measured with an ATP Bioluminescence assay kit HSII (Roche Molecular Biochemicals, Germany) as described earlier [48]. Intracellular levels of ATP were expressed per 1 g of CDW.


H2 accumulation (mmol H2/l) was calculated in two different ways depending on the experimental design. When N2 sparging was applied the calculations were based on the flow rate of the influent N2 gas and the percentages of H2 and CO2 in the effluent gas, whereas when N2 sparging was not applied the flow rate of the effluent gas was measured by the liquid displacement method using a NaCO3-saturated solution to avoid any further CO2 solubilization. Molar H2 and CO2 were calculated using the ideal gas law based on their concentration in the effluent gas. Volumetric H2 productivity (Q H2; mmol H2/l/h) was determined from the slope of the straight line obtained by plotting the cumulative H2 (mmol H2/l) against the time (h) during exponential growth.

Carbon balances were calculated from the total amount of carbon-containing products formed (in C-mol) and the amount of sugar consumed (in C-mol). Electron balances were calculated after multiplying the amount of each metabolite and sugar by the corresponding degree of reduction (in mol electrons per C-mol). The chemical formula of biomass was assumed to be CH1.8O0.5 N0.2. Metabolite yields as well as carbon and electron balances were calculated by subtracting the background metabolite production in control cultivations (i.e. cultivation without the substrate) and the precultures’ carryover metabolites from the results.

The yield of biomass per ATP (YX/ATP; gCDW/mol ATP) and the ATP yield (YATP/S; mol ATP/mol glucose) were calculated based on the Eq. (5) and Eq. (6), respectively.

Y X / ATP = biomass 2 * acetate + 1 * ethanol
Y ATP / S = 2 * acetate + 1 * ethanol glucose consumed

The calculation assumes that C. celer uses a PEP-dependent phosphotransferase system (PTS) for glucose uptake. Specific rates of substrate consumption uptake and metabolite formation (r) were calculated according to

r = dC dt * X

where r is the specific production rate (mmol/gCDW/h), C is the substrate or metabolite concentration (mM) and X is the biomass (gCDW).

In silico model construction and metabolic flux analysis

The genome of C. celer strain JW/YL-NZ35 has been recently sequenced and a high quality draft sequence has been annotated as described earlier [18]. The available genomic data was used to reconstruct the central carbon metabolism of C. celer and to build a stoichiometric model including glycolytic pathway, pentose phosphate pathway, end product synthesis pathways, transport mechanisms and biomass synthesis. In order to define each reaction and metabolite in the network, genes encoding for putative enzymes involved in the aforementioned pathways were identified by manual curation of the data from the annotated genome by a combination of sequence alignments [50], gene context analysis, database and literature searches. In addition, physiological evidences were used either when association between genes and reactions could not be established or when a reaction was removed from the network despite a corresponding gene being identified in the genome. The final model consisted of 24 reactions and 20 intracellular metabolites (Additional file 1: Table S1).

Metabolic flux analysis (MFA) was employed to calculate the unknown intracellular fluxes at different stages of C. celer´s growth (exponential, transition and stationary phase) under the tested conditions. A metabolic matrix was constructed based on the model according to the law of mass conservation and on the pseudo-steady state hypothesis for the intracellular metabolites. Given the overdetermined nature of the system, the mass balance equations for all the metabolites were expressed in matrix form as:

V c = - G c * × G m × V m

where Vc is the calculated flux vector, Gc * is pseudo-inverse matrix of the calculated reactions, Gm is the matrix of the measured reactions and Vm the measured flux vector. The matrix of MFA model was solved with MATLAB R2013a (The MathWorks, Inc., Natick, USA).


  1. Kapdan IK, Kargi F: Bio-hydrogen production from waste materials. Enzyme Microb Technol. 2006, 38: 569-582. 10.1016/j.enzmictec.2005.09.015.

    Article  CAS  Google Scholar 

  2. Lattin WC, Utgikar VP: Transition to hydrogen economy in the United States: A 2006 status report. Int J Hydrogen Energy. 2007, 32: 3230-3237. 10.1016/j.ijhydene.2007.02.004.

    Article  CAS  Google Scholar 

  3. Hallenbeck PC, Ghosh D: Advances in fermentative biohydrogen production: the way forward?. Trends Biotechnol. 2009, 27: 287-297.

    Article  CAS  Google Scholar 

  4. Rittmann S, Herwig C: A comprehensive and quantitative review of dark fermentative biohydrogen production. Microb Cell Fact. 2012, 11: 115-

    Article  CAS  Google Scholar 

  5. Verhaart MRA, Bielen AAM, van der Oost J, Stams AJM, Kengen SWM: Hydrogen production by hyperthermophilic and extremely thermophilic bacteria and archaea: mechanisms for reductant disposal. Environ Technol. 2010, 31: 993-1003.

    Article  CAS  Google Scholar 

  6. Baena S, Patel B: Genus V: Caloramator Collins, Lawson, Willems, Cordoba, Fernández-Garayzábal, Garcia, Cai, Hippe and Farrow 1994, 812VP emend. Chrisostomos, Patel, Dwivedi and Denman 1996, 497. Bergey's manual of systematic bacteriology. Volume 3: The Firmicutes. Edited by: De Vos P, Garrity G, Jones D, Krieg N, Ludwig W, Rainey F, Schleifer K, Whitman W. 2009, 834-838. New York, NY: Springer-Verlag, 2,

    Google Scholar 

  7. Koskinen PEP, Lay C, Puhakka JA, Lin P, Wu S, Örlygsson J, Lin C: High-efficiency hydrogen production by an anaerobic, thermophilic enrichment culture from an icelandic hot spring. Biotechnol Bioeng. 2008, 101: 665-678.

    Article  CAS  Google Scholar 

  8. Ciranna A, Santala V, Karp M: Biohydrogen production in alkalithermophilic conditions: Thermobrachium celere as a case study. Bioresour Technol. 2011, 102: 8714-8722.

    Article  CAS  Google Scholar 

  9. Ciranna A, Santala V, Karp M: Enhancing biohydrogen production of the alkalithermophile Thermobrachium celere. Int J Hydrogen Energy. 2012, 37: 5550-5558. 10.1016/j.ijhydene.2011.12.105.

    Article  CAS  Google Scholar 

  10. Kengen SM, Stams AJM, de Vos WM: Sugar metabolism of hyperthermophiles. FEMS Microbiol Rev. 1996, 18: 119-137. 10.1111/j.1574-6976.1996.tb00231.x.

    Article  CAS  Google Scholar 

  11. Schröder C, Selig M, Schönheit P: Glucose fermentation to acetate, CO2 and H2 in the anaerobic hyperthermophilic eubacterium Thermotoga maritima: involvement of the Embden-Meyerhof pathway. Arch Microbiol. 1994, 161: 460-470.

    Google Scholar 

  12. Soboh B, Linder D, Hedderich R: A multisubunit membrane-bound [NiFe] hydrogenase and an NADH-dependent Fe-only hydrogenase in the fermenting bacterium Thermoanaerobacter tengcongensis. Microbiology. 2004, 150: 2451-2463.

    Article  CAS  Google Scholar 

  13. van Niel EWJ, Claassen PAM, Stams AJM: Substrate and product inhibition of hydrogen production by the extreme thermophile, Caldicellulosiruptor saccharolyticus. Biotechnol Bioeng. 2003, 81: 255-262.

    Article  CAS  Google Scholar 

  14. Ljunggren M, Willquist K, Zacchi G, Van Niel EWJ: A kinetic model for quantitative evaluation of the effect of hydrogen and osmolarity on hydrogen production by Caldicellulosiruptor saccharolyticus. Biotechnol Biofuels. 2011, 4: 31-

    Article  CAS  Google Scholar 

  15. Zhang F, Zhang Y, Chen M, Zeng RJ: Hydrogen supersaturation in thermophilic mixed culture fermentation. Int J Hydrogen Energy. 2012, 37: 17809-17816. 10.1016/j.ijhydene.2012.09.019.

    Article  CAS  Google Scholar 

  16. Wang J, Wan W: Factors influencing fermentative hydrogen production: A review. Int J Hydrogen Energy. 2009, 34: 799-811. 10.1016/j.ijhydene.2008.11.015.

    Article  CAS  Google Scholar 

  17. Hawkes FR, Dinsdale R, Hawkes DL, Hussy I: Sustainable fermentative hydrogen production: Challenges for process optimisation. Int J Hydrogen Energy. 2002, 27: 1339-1347. 10.1016/S0360-3199(02)00090-3.

    Article  CAS  Google Scholar 

  18. Ciranna A, Larjo A, Kivistö A, Santala V, Roos C, Karp M: Draft Genome Sequence of the Hydrogen- and Ethanol-Producing Anaerobic Alkalithermophilic Bacterium Caloramator celer. Genome Announc. 2013, 1: e00471-13.

    Article  Google Scholar 

  19. Carere CR, Rydzak T, Verbeke TJ, Cicek N, Levin DB, Sparling R: Linking genome content to biofuel production yields: A meta-analysis of major catabolic pathways among select H2 and ethanol-producing bacteria. BMC Microbiology. 2012, 12: 295-

    Article  CAS  Google Scholar 

  20. Sapra R, Bagramyan K, Adams MWW: A simple energy-conserving system: Proton reduction coupled to proton translocation. Proc Natl Acad Sci USA. 2003, 100: 7545-7550.

    Article  CAS  Google Scholar 

  21. Engle M, Li Y, Rainey F, DeBlois S, Mai V, Reichert A, Mayer F, Messner P, Wiegel J: Thermobrachium celere gen. nov., sp. nov., a rapidly growing thermophilic, alkalitolerant, and proteolytic obligate anaerobe. Int J Syst Bacteriol. 1996, 46: 1025-1033.

    Article  CAS  Google Scholar 

  22. Ciranna A, Ferrari R, Santala V, Karp M: Inhibitory effects of substrate and soluble end products on biohydrogen production of the alkalithermophile Caloramator celer: kinetic, metabolic and transcription analyses. Int J Hydrogen Energy. 2014, 39: 6391-6401. 10.1016/j.ijhydene.2014.02.047.

    Article  CAS  Google Scholar 

  23. Kobayashi H: Regulation of cytoplasmic pH in streptococci. Sugar transport and metabolism in Gram-positive bacteria. Edited by: Reizer J, Peterkofsky A. 1987, London, United Kingdom: Ellis Harwood

    Google Scholar 

  24. O'Sullivan E, Condon S: Relationship between acid tolerance, cytoplasmic pH, and ATP and H+-ATPase levels in chemostat cultures of Lactococcus lactis. Appl Environ Microbiol. 1999, 65: 2287-2293.

    Google Scholar 

  25. Kessler D, Knappe J: Anaerobic dissimilation of pyruvate. Escherichia coli and Salmonella: cellular and molecular biology. Edited by: Neidhardt FC, Curtiss RIII, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schaechter M, Umbarger HE. 1996, 199-204. Washington, DC: ASM Press, 2

    Google Scholar 

  26. Sridhar J, Eiteman MA: Metabolic flux analysis of Clostridium thermosuccinogenes: Effects of pH and culture redox potential. Appl Biochem Biotechnol. 2001, 94: 51-69.

    Article  CAS  Google Scholar 

  27. Liu IC, Whang LM, Ren WJ, Lin PY: The effect of pH on the production of biohydrogen by clostridia: Thermodynamic and metabolic considerations. Int J Hydrogen Energy. 2011, 36: 439-449. 10.1016/j.ijhydene.2010.10.045.

    Article  Google Scholar 

  28. Combet-Blanc Y, Kalamba KK, Kergoat PY: Effect of pH on Bacillus thermoamylovorans growth and glucose fermentation. Appl Environ Microbiol. 1995, 61: 656-659.

    CAS  Google Scholar 

  29. Su Y, Rhee MS, Ingram LO, Shanmugam KT: Physiological and fermentation properties of Bacillus coagulans and a mutant lacking fermentative lactate dehydrogenase activity. J Ind Microbiol Biotechnol. 2011, 38: 441-450.

    Article  CAS  Google Scholar 

  30. Koussémon M, Combet-Blanc Y, Ollivier B: Glucose fermentation by Propionibacterium microaerophilum: Effect of pH on metabolism and bioenergetic. Curr Microbiol. 2003, 46: 141-145.

    Article  Google Scholar 

  31. Asanuma N, Iwamoto M, Hino T: Structure and transcriptional regulation of the gene encoding pyruvate formate-lyase of a ruminal bacterium, Streptococcus bovis. Microbiology. 1999, 145: 151-157.

    Article  CAS  Google Scholar 

  32. Even S, Lindley ND, Cocaign-Bousquet M: Transcriptional, translational and metabolic regulation of glycolysis in Lactococcus lactis subsp. cremoris MG 1363 grown in continuous acidic cultures. Microbiology. 2003, 149: 1935-1944.

    Article  CAS  Google Scholar 

  33. Temudo MF, Kleerebezem R, Van Loosdrecht M: Influence of the pH on (Open) mixed culture fermentation of glucose: A chemostat study. Biotechnol Bioeng. 2007, 98: 69-79.

    Article  CAS  Google Scholar 

  34. Lee H, Rittmann BE: Evaluation of metabolism using stoichiometry in fermentative biohydrogen. Biotechnol Bioeng. 2009, 102: 749-758.

    Article  CAS  Google Scholar 

  35. Thauer RK, Kirchniawy FH, Jungermann KA: Properties and function of the pyruvate-formate-lyase reaction in clostridiae. Eur J Biochem. 1972, 27: 282-290.

    Article  CAS  Google Scholar 

  36. Takahashi-Abbe S, Abe K, Takahashi N: Biochemical and functional properties of a pyruvate formate-lyase (PFL)-activating system in Streptococcus mutans. Oral Microbiol Immunol. 2003, 18: 293-297.

    Article  CAS  Google Scholar 

  37. Knappe J, Blaschkowski HP, Gröbner P, Schmitt T: Pyruvate Formate-Lyase of Escherichia coli: the Acetyl-Enzyme Intermediate. Eur J Biochem. 1974, 50: 253-263.

    Article  CAS  Google Scholar 

  38. Kerscher L, Oesterhelt D: The catalytic mechanism of 2-oxoacid: Ferredoxin oxidoreductases from Halobacterium halobium. One electron transfer at two distinct steps of the catalytic cycle. Eur J Biochem. 1981, 116: 595-600.

    Article  CAS  Google Scholar 

  39. Williams K, Lowe PN, Leadlay PF: Purification and characterization of pyruvate: ferredoxin oxidoreductase from the anaerobic protozoon Trichomonas vaginalis. Biochem J. 1987, 246: 529-536.

    Article  CAS  Google Scholar 

  40. Pineda E, Encalada R, Rodríguez-Zavala JS, Olivos-García A, Moreno-Sánchez R, Saavedra E: Pyruvate:ferredoxin oxidoreductase and bifunctional aldehyde-alcohol dehydrogenase are essential for energy metabolism under oxidative stress in Entamoeba histolytica. FEBS Journal. 2010, 277: 3382-3395.

    Article  CAS  Google Scholar 

  41. Asanuma N, Hino T: Effects of pH and energy supply on activity and amount of pyruvate formate-lyase in Streptococcus bovis. Appl Environ Microbiol. 2000, 66: 3773-3777.

    Article  CAS  Google Scholar 

  42. Payot S, Guedon E, Cailliez C, Gelhaye E, Petitdemange H: Metabolism of cellobiose by Clostridium cellulolyticum growing in continuous culture: Evidence for decreased NADH reoxidation as a factor limiting growth. Microbiology. 1998, 144: 375-384. 10.1099/00221287-144-2-375.

    Article  CAS  Google Scholar 

  43. Payot S, Guedon E, Gelhaye E, Petitdemange H: Induction of lactate production associated with a decrease in NADH cell content enables growth resumption of Clostridium cellulolyticum in batch cultures on cellobiose. Res Microbiol. 1999, 150: 465-473.

    Article  CAS  Google Scholar 

  44. Willquist K, van Niel EWJ: Lactate formation in Caldicellulosiruptor saccharolyticus is regulated by the energy carriers pyrophosphate and ATP. Metab Eng. 2010, 12: 282-290.

    Article  CAS  Google Scholar 

  45. Rydzak T, Levin DB, Cicek N, Sparling R: End-product induced metabolic shifts in Clostridium thermocellum ATCC 27405. Appl Microbiol Biotechnol. 2011, 92: 199-209.

    Article  CAS  Google Scholar 

  46. Carere C, Rydzak T, Cicek N, Levin D, Sparling R: Role of transcription and enzyme activities in redistribution of carbon and electron flux in response to N2 and H2 sparging of open-batch cultures of Clostridium thermocellum ATCC 27405. Appl Microbiol Biotechnol. 2014, 98: 2829-2840.

    Article  CAS  Google Scholar 

  47. Ljunggren M, Zacchi G: Techno-economic analysis of a two-step biological process producing hydrogen and methane. Bioresour Technol. 2010, 101: 7780-7788.

    Article  CAS  Google Scholar 

  48. Willquist K, Pawar S, Van Niel E: Reassessment of hydrogen tolerance in Caldicellulosiruptor saccharolyticus. Microb Cell Fact. 2011, 10: 111-

    Article  CAS  Google Scholar 

  49. Zeidan AA, van Niel EWJ: A quantitative analysis of hydrogen production efficiency of the extreme thermophile Caldicellulosiruptor owensensis OLT. Int J Hydrogen Energy. 2010, 35: 1128-1137. 10.1016/j.ijhydene.2009.11.082.

    Article  CAS  Google Scholar 

  50. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402.

    Article  CAS  Google Scholar 

Download references


Alessandra Radici is kindly acknowledged for language editing of the manuscript. AC acknowledges support from the Academy of Finland (grants no.126974, 139830, 272602), the Finnish Doctoral Programme in Environmental Science and Technology (EnSTe) and Tampere University of Technology Graduate School. SSP acknowledges support from the Swedish Research Council (VR).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Alessandro Ciranna.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AC, SSP and EWJvN designed the study. AC and SSP planned and preformed the batch experiments. AC planned and preformed NAD(H) and ATP quantifications and metabolite analysis. AC and SSP performed metabolic flux analysis. AC wrote the manuscript. SSP, VS, MK and EWJvN participated in manuscript drafting. VS, MK and EWJvN supervised and coordinated the study. All authors have read and approved the manuscript.

Electronic supplementary material


Additional file 1: Table S1: Reconstructed genomic model of the central carbon metabolism in Caloramator celer. (DOCX 32 KB)


Additional file 2: Figure S1: Effect of the H2 concentration on the Gibbs energy change (ΔG´) at 70°C of four predicted reactions involved in glucose fermentation in C. celer. (DOCX 912 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Ciranna, A., Pawar, S.S., Santala, V. et al. Assessment of metabolic flux distribution in the thermophilic hydrogen producer Caloramator celer as affected by external pH and hydrogen partial pressure. Microb Cell Fact 13, 48 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: