A kinetic model for quantitative evaluation of the effect of hydrogen and osmolarity on hydrogen production by Caldicellulosiruptor saccharolyticus
Biotechnology for Biofuels volume 4, Article number: 31 (2011)
Caldicellulosiruptor saccharolyticus has attracted increased interest as an industrial hydrogen (H2) producer. The aim of the present study was to develop a kinetic growth model for this extreme thermophile. The model is based on Monod kinetics supplemented with the inhibitory effects of H2 and osmotic pressure, as well as the liquid-to-gas mass transfer of H2.
Mathematical expressions were developed to enable the simulation of microbial growth, substrate consumption and product formation. The model parameters were determined by fitting them to experimental data. The derived model corresponded well with experimental data from batch fermentations in which the stripping rates and substrate concentrations were varied. The model was used to simulate the inhibition of growth by H2 and solute concentrations, giving a critical dissolved H2 concentration of 2.2 mmol/L and an osmolarity of 0.27 to 29 mol/L. The inhibition by H2, being a function of the dissolved H2 concentration, was demonstrated to be mainly dependent on H2 productivity and mass transfer rate. The latter can be improved by increasing the stripping rate, thereby allowing higher H2 productivity. The experimentally determined degree of oversaturation of dissolved H2 was 12 to 34 times the equilibrium concentration and was comparable to the values given by the model.
The derived model is the first mechanistically based model for fermentative H2 production and provides useful information to improve the understanding of the growth behavior of C. saccharolyticus. The model can be used to determine optimal operating conditions for H2 production regarding the substrate concentration and the stripping rate.
The question of climate change, together with the increasing scarcity and cost of fossil fuels, has triggered research on the sustainable production of energy carriers, such as biofuels. Although a wide spectrum of alternative fuels and processes are available, it is still not known which of them will succeed in replacing fossil fuels in the long term. Processes for producing some biofuels, such as ethanol and methane, are already highly developed for commercial use, while others, including biohydrogen, require more basic research before they can be produced in an economically feasible way .
The current major drawbacks of fermentative hydrogen (H2) production are its low yield  and the requirement of gas stripping to remove H2 from the liquid . Dark fermentation can be performed at either moderate or elevated temperatures: the productivity is generally higher in the former case , whereas the latter provides higher yields . However, for an economically sustainable biohydrogen process, high productivity should be accompanied by high yields . One solution is to focus on the process design, such as the optimization of H2 removal and substrate concentration . A higher substrate concentration is required to avoid excess amounts of water, since a surplus of process water has a negative effect on both the environmental and economic aspects of the process . It is important to find alternatives to nitrogen (N2) sparging, as this dilutes the H2, leading to a lower-grade product or the need for expensive gas upgrading to remove the N2.
The extreme thermophile Caldicellulosiruptor saccharolyticus may be of industrial interest because of its ability to produce high yields of H2 from a wide variety of sugars [8, 9], ranging from pentose and hexose monomers  to complex (hemi)cellulosic materials [8, 11, 12]. However, the growth of this microorganism, as well as the production of H2 and acetic acid, is subject to substrate and product inhibition . In addition, it is sensitive to increased osmotic pressure , making it necessary to dilute the substrate. Moreover, C. saccharolyticus is inhibited by H2 , which is a general trait of H2-producing organisms. A high dissolved H2 concentration (H2aq) inhibits H2-generating hydrogenases [1, 2], leading to increased nicotinamide adenine dinucleotide, reduced/nicotinamide adenine dinucleotide (NADH/NAD) ratios, resulting in a metabolic shift toward reduced products such as lactate and ethanol [1, 3, 4] and thus decreased H2 productivity.
The critical partial H2 pressure () is usually the parameter coupled to growth inhibition and lactate formation [5, 13, 15–17], although it has been shown that there are mass transfer limitations that cause the fermentation medium to be supersaturated with H2 [18, 19]. Therefore, H2aq should be considered rather than when discussing inhibition by H2. Stripping with N2, which decreases H2aq , can significantly improve the H2 yield . Indeed, optimizing the process should consist of a compromise between minimizing the stripping gas flow rate, henceforth called the "stripping rate," and maintaining a suitable H2 yield and productivity . However, optimization must take H2 productivity into account, since it influences H2aq greatly; otherwise, the optimization will be relevant only to the specific system studied and cannot be applied to other biohydrogen systems. Furthermore, better insight into H2 productivity and liquid-to-gas mass transfer will help considerably in finding technical solutions allowing an economical and feasible fermentation process to be designed.
This paper describes the development of a kinetic Monod-based model for cell growth and product formation by C. saccharolyticus, which takes the liquid-to-gas mass transfer into consideration. The derived model successfully simulates batch fermentations of C. saccharolyticus and predicts the influence of H2aq and osmotic pressure on growth as well as the metabolic shift to lactate production. To the best of our knowledge, this is the first time the inhibitory effects and liquid-to-gas mass transfer of H2 and carbon dioxide (CO2) have been included in a model for fermentative H2 production. The results obtained with the model are relevant for biohydrogen systems other than that described herein using C. saccharolyticus.
Cultivation and fermentation
C. saccharolyticus was cultivated in a N2 atmosphere in modified DSM640 medium as described previously . Cultures for inocula were grown overnight in 250-mL serum flasks containing 50 mL of modified DSM640 medium, together with a 0.4% carbon source and 0.02% cysteine.
Fermentations were performed in a jacketed 3-L reactor (Applikon, Schiedam, the Netherlands) at a working volume of 1 L. The pH was monitored using an ADI 1025 Bio Console Bio Controller (Applikon) and maintained at pH 6.6 (corresponding to neutral pH at 70°C) by the addition of sodium hydroxide (NaOH). The temperature was thermostatically maintained at 70°C ± 1°C, and the stirring rate was set to 350 rpm. Prior to inoculation, the medium was reduced by the addition of 0.1% cysteine. The fermentor was continuously stripped with N2.
A total of eight experiments were performed, with glucose concentrations of 5 or 10 g/L and stripping rates from 0.78 to 6 L of N2/hour, to estimate the parameters of the kinetic model. The experiments were performed on four different occasions. An overview of the experiments is given in Table 1. The H2 productivity and the cumulative H2 production were determined as described previously .
Measurement of dissolved hydrogen concentration
The dissolved H2 concentration, H2aq, was determined in samples of cell suspensions from the bioreactor in serum flasks. The pressure in the flasks was reduced using a vacuum pump to withdraw the fermentation medium (30 to 40 mL) from the fermentor. The sealed flasks were placed in a water bath (20°C) with magnetic stirring to ensure good mixing. After one hour in the water bath, the pressure in the flasks was rapidly raised to atmospheric pressure by inserting a small needle through the cap, allowing air to enter the flasks. After 12 hours, the H2 concentration in the headspace of the serum flasks was measured as described below (Analytical methods). The volume of the collected liquid and the total volume of the individual serum flasks were measured to estimate the original H2aq in the sample.
Headspace samples were analyzed for CO2 and H2 using a Varian CP-4900 Micro GC gas chromatograph (Varian, Inc., Middelburg, the Netherlands) equipped with a thermal conductivity detector (100 mA). The results were analyzed using Galaxie Chromatography Workstation version 188.8.131.52 software (Varian, Inc.). The optical density of aerobic samples was measured at 620 nm using a spectrophotometer (Hitachi U-1100; Hitachi High-Technologies Corp., Tokyo, Japan). Concentrations of acetic acid, lactic acid and ethanol were analyzed using high-pressure liquid chromatography (Waters Corp., Milford, MA, USA). An Aminex HPX-87H ion exchange column (Bio-Rad Laboratories, Hercules, CA, USA) was used at 45°C with a mobile phase of 5 mM H2SO4 at a flow rate of 0.6 mL/minute. The chromatograph was equipped with a refractive index detector (RID-6A; Shimadzu, Kyoto, Japan).
The model employed in this study considers the kinetics of growth and product formation, including inhibition, liquid-to-gas mass transfer and the chemical equilibrium of CO2 and carbonates. The constants used in the model are listed in Table 2.
Growth and product kinetics in C. saccharolyticus
The kinetic model for C. saccharolyticus is based on simultaneous solution of the mass balance equations. The following reactions occurring in C. saccharolyticus were taken into consideration in the model:
The elemental composition of C. saccharolyticus was previously determined by de Vrije et al. . The value of the stoichiometric parameter k in equation 3 is not known and was therefore estimated (described in detail below, Experimental design and estimation of kinetic parameters). Furthermore, reaction 3 is not balanced, since elements are available in the fermentation medium which were not considered in the model (for example, yeast extract). C. saccharolyticus is subject to several inhibiting agents affecting cell growth as well as the production of acetic acid and H2. These include (1) elevated concentrations of dissolved H2 and H2aq, and (2) high osmolarity. Osmolarity (named "OSM" in the model) is calculated by summing the molar can be described by the following expression:
where G, Ac and Lac are the concentrations of glucose, acetic acid and lactic acid, respectively, and CO2sol is the concentration of bicarbonate and carbonate, resulting from the CO2 produced by the bacterium. The stoichiometric factor 2 arises from the assumption that for each mole of acid produced, one mole of NaOH is added to maintain the pH. The background osmolarity, 0.10 mol/L, resulting from yeast extract and other nutrients was calculated based on osmolarity measurements presented previously .
The mass balance for the growth of the cells (in batch fermentation) can be described by the following expression:
where μ, X and rcd are the specific growth rate, the cell mass concentration and the cell death rate, respectively.
The specific growth rate, μ, is given by Monod kinetics:
Here μmax is the maximum specific growth rate, G is the glucose concentration and KG is the glucose affinity parameter. OSMcrit, H2aqcrit, and nμ are the critical osmolarity, critical dissolved H2 concentration and the exponential parameters describing the level of inhibition, respectively. The critical parameters represent the values at which inhibition is 100%. The kinetic inhibition expressions, as presented by Han and Levenspiel , are valid only for concentrations of the inhibiting agents (OSM or H2aq) lower than the respective critical concentrations. Elevated osmolarity inhibits the growth of C. saccharolyticus , which may lead to a metabolic shift toward the production of lactic acid .
The product formation rates were expressed using an adapted form of the Luedeking-Piret equation :
Here, α is a growth-associated constant, and β is a non-growth-associated constant. All products were assumed to be only growth-associated, hence β is zero. This simplification is based on the assumption that other nutrients are in excess. It is further assumed that H2 inhibits its own production and thus also that of acetate in both a direct and an indirect way: indirectly through the inhibition of growth, and directly through specifically inhibiting the hydrogenases, which results in an increase in the NADH/NAD ratio in the cell . This in turn induces the activity of the enzyme lactate dehydrogenase, initiating a metabolic shift toward lactate . Besides being activated through H2, lactic acid production is also activated via osmolarity, which at high concentrations inhibits growth . The growth-associated constants (equation 7) for the metabolic products are as follows:
Here RAcF, and RLacFOSM are the factors relating the production rates of acetic acid and lactic acid to the growth rate; INHIB and ACT are the inhibition and activation factors, respectively (see below, Product formation) YGAc and are the stoichiometric yield coefficients of the reactions described above (equation 1); and is the estimated H2 yield coefficient (mol H2/mol glucose). The following inhibition mechanisms were studied.
Inhibition and activation by dissolved H2:
Activation of lactic acid production by osmolarity:
The mass balances of the products in the liquid phase consist of an accumulation, time-dependent, and a product formation rate, concentration-dependent, term, and, when applicable, a liquid-to-gas mass transfer term.
Here is the overall volumetric mass transfer coefficient for H2 and is the H2 saturation concentration.
where is the overall volumetric mass transfer coefficient for CO2; CO2aq and are the concentration of dissolved CO2 and the saturation concentration of CO2, respectively; K1 and K2 are the dissociation constants; and CO2sol is the total concentration of carbonates. The last two expressions in the CO2 mass balance (equation 18) are related to the reaction of CO2 with water to form carbonates (see Additional file: Liquid-to-gas mass transfer and chemistry of CO2__Theory_experiments_and_results). The expressions and in equations 17 and 18 describe the liquid-to-gas mass transfer of H2 and CO2, respectively (see Additional file: Liquid-to-gas mass transfer and chemistry of CO2__Theory_experiments_and_results).
The mass balance of glucose can be written as follows:
where YGX is the cell mass yield coefficient parameter (which is the model notation of the constant k in equation 3) and YGAc and YGLac are the stoichiometric yield coefficients (mol/mol substrate) describing the theoretical yields of cell mass, acetic acid and lactic acid, respectively.
Gas phase balances
The gaseous products are produced in the liquid phase, after which they are transported to the gas phase. The mass balances of the various gaseous compounds comprise an accumulation term and a gas-to-liquid mass transfer term as well as an inflow and an outflow.
Here V l and Vg are the liquid and gas volume; k l a is the volumetric mass transfer coefficient for the various compounds; H2aq and CO2aq are the dissolved concentrations of H2 and CO2; and are the dissolved concentrations of the compounds at equilibrium; is the flow rate of N2 into the fermentor; and , and are the effluent flow rates of N2, CO2 and H2, respectively.
Other equations used
The equilibrium between CO2, bicarbonate and carbonate was taken into consideration as follows:
Since the pH was controlled during fermentation, it was kept constant at 6.6 in the model. For more details, see Additional file 1: Liquid-to-gas mass transfer and chemistry of CO2__Theory_experiments_and_results.
Furthermore, the effect of the stripping rate on the volumetric mass transfer coefficients was determined by using the following equation:
where k l a i0 is the volumetric mass transfer coefficient at a stripping rate of , k l a i is the volumetric mass transfer coefficient at a stripping rate of and γ is an experimentally determined exponential coefficient.
Experimental design and estimation of kinetic parameters
The mass transfer parameters ( and ) and the exponential coefficient (γ) were determined before the parameters for the microbial kinetics were estimated. The mass transfer parameters were estimated using a methodology described by Hill . The volumetric mass transfer coefficient for CO2, , was determined and the volumetric mass transfer coefficient for H2 () was calculated using the following relation between k l a values and the diffusion coefficients for the compounds in water :
Here the diffusion coefficients for CO2 () and H2 () are 1.98 × 10-5 and 4.65 × 10-5 cm2/second, respectively, according to Pauss et al. . Only the most important results in the estimation of the mass transfer parameters are given here; the detailed results and the methodology can be found in Additional file: Liquid-to-gas mass transfer and chemistry of CO2__Theory_experiments_and_results.
To determine the kinetic parameters, fermentations were performed at two substrate concentrations and various stripping rates (Table 1). The kinetic parameters studied and estimated were YGX, , μmax, RAcF, , H2aqcrit, , RLacFOSM, OSMcrit and nμ. Since the experiments were characterized by different growth rates as well as H2 and biomass yields, all parameters except the three related to osmolarity (RLacFOSM, OSMcrit and nμ) were estimated simultaneously for each of the 5 g/L experiments (that is, experiments 1 and 3 through 7), giving five parameter sets (Table 3). Experiments 6 and 7 were carried out under almost the same conditions with the same inocula and were thus used together to estimate one parameter set. The effect of osmolarity was assumed to be negligible in these experiments (experiments 1 and 3 through 7). The validity of this assumption was checked and is discussed below (Kinetic parameters: inhibition by osmolarity). An average of all parameters except μmax was then calculated and used to simulate all six experiments. The parameters related to osmolarity were estimated using the 10 g/L experiments (that is, experiments 2 and 8).
Finally, the cell death rate (rcd; equation 5), was calculated directly from the experimental data by determining, for each experiment, the slope of the biomass curve as it declined and calculating an average value of the slope based on all experiments. The overall affinity constant (KG; equation 6) was taken from the literature on bacteria with the same type of glucose transporter system, that is, the ATP-binding cassette transport system (values are given in Table 2) [10, 27, 28].
The parameters in the model were estimated using a nonlinear least squares method which minimizes the objective function:
where LSQ is the least squares quadratic, that is, the sum of the squared residuals, Yexp,i,tand Ymodel,i,t, are the experimental and model values, respectively, at time t for each compound i, that is, glucose, acetic acid, lactic acid, cell mass and H2. The parameter estimations were carried out and the 95% confidence intervals were calculated using the MATLAB functions nlinfit and nlparci, respectively (MathWorks, Natick, MA, USA).
Results and Discussion
The order in which the parameters are estimated is important. Since the mass transfer influences fermentation, and thus the corresponding kinetic parameters, the mass transfer parameters were estimated prior to determining the fermentation-related kinetic parameters.
Estimation of the volumetric mass transfer coefficient
The estimated overall volumetric mass transfer coefficients for CO2 () in the fermentation medium were in the range of 3.4 to 7.2/hour for stripping rates ranging from 2.0 to 9.8 L/hour (33 to 164 mL/minute). The values obtained with their respective stripping rates were used to calculate the exponential parameter γ (equation 24), giving a value of 0.46 ± 0.02, which agrees well with previous findings [20, 29]. On the basis of the mass transfer results of CO2 and equation 25 was determined to be 9.0 ± 0.1/hour at a stripping rate of 6 L/hour (100 mL/minute). This value is consistent with previous results for similar systems with stripping rates between 5 and 2,000 mL/minute, where the k l a values were in the range of 1 to 120/hour [18, 20]. Hence, the expression calculating (equation 24) is:
where is the stripping rate. For more detailed results on mass transfer, see Additional file: Liquid-to-gas mass transfer and chemistry of CO2__Theory_experiments_and_results.
In fermentation experiments 1 through 7, the glucose was completely consumed and acetic acid, H2 and lactic acid were the main products (Figures 1, 2, 3, 4, 5, 6 and 7). In experiment 8, complete glucose consumption was not achieved within 40 hours (Figure 8).
The six 5 g/L experiments (that is, experiments 1 and 3 through 7) differed with respect to growth rate and H2 yield (Figures 1, 2, 3, 4, 5 and 6). The reason for the variations is not clear. However, the two sets of experiments, carried out in parallel using inocula from the same preculture (experiments 3 and 4 as well as experiments 5 through 8), showed a strong resemblance. Batch experiments are generally characterized by poor reproducibility, as has also been observed with Escherichia coli . Owing to the variations between the experiments, each of the 5 g/L experiments was used separately to estimate all parameters except OSMcrit, nμ and RLacFOSM, which were estimated using the 10 g/L experiments (experiments 2 and 8 and Figures 7 and 8). In addition, acetic acid and H2 were not produced at a stoichiometric ratio (1:2) in any of the experiments, which is in agreement with prior observations in continuous culture [15, 21] where the H2/acetic acid ratio depended on the dilution rate. Instead, less acetate was produced in relation to H2, and the H2 yield was thus considered a parameter in the model. A possible explanation of the nonstoichiometric yield is that acetic acid is used as a building block by the cells.
The estimated maximum growth rate of C. saccharolyticus varied between the different experiments; however, the values were similar for experiments carried out in parallel (that is, experiments 3 and 4 and 5 through 8) (Table 3). The μmax determined for glucose ranged from 0.21 to 0.36/hour and corresponded well with earlier findings, which varied from 0.2 to 0.5/hour [14, 21]. Vrije et al.  found different values of μmax obtained in washout experiments with C. saccharolyticus, depending on the growth history of the culture. The reason for the variation in growth rates has not yet been elucidated, but the quality of the inocula used seems to depend on its growth history. The variation in growth rate was considered in the present model by not using an average value of μmax, while average values were used for all the other parameters (described in more detail below, Kinteic parameters: inhibition by dissolved hydrogen).
Together with changes in the product yields, the variations in growth rates strongly suggest that an as yet unknown factor influences the metabolism. Also, it is likely that the highest value of μmax determined in this study (0.36/hour) is not the true maximum, that is to say, the real μmax that can be achieved for glucose by this organism. This could mean that the growth of C. saccharolyticus is normally suppressed, as has been observed for E. coli . The cause of the suppression of the growth of C. saccharolyticus is not understood and must be further studied to gain a better understanding of the physiology of the microorganism and to improve the model .
Kinetic parameters: inhibition by dissolved hydrogen
As described above (Growth and product kinetics in C. saccharolyticus), the effects of two environmental conditions, that is to say, H2aq and osmolarity, on the growth, the production of acetic acid and H2, and the activation of lactic acid production of C. saccharolyticus were studied. All parameters except and RLacFOSM, OSMcrit and nμ were estimated for each of the experiments with 5 g/L glucose (Table 3), from which their average values were calculated (Table 4). The experimental results are presented in Figures 1, 2, 3, 4, 5, and 6, together with the predictions of the model using the average values of the parameters (except for μmax). The average value of μmax was not used because of the large variation in this parameter; instead, the average within each set of parallel experiments was used. The model then agreed well with the experimental results and also reflected the effect of H2aq on the microorganism. Interestingly, the parameter has quite a high value (4.5), which indicates a steplike response to the inhibiting agent, H2; that is, at lower concentrations, the inhibition increases only a little with the concentration, whereas at higher concentrations the inhibition increases dramatically with increased concentration. The value of the H2 yield coefficient parameter, , ranged from 3.7 to 6.2 mol H2/mol glucose. The upper value exceeds the Thauer limit, which is 4 mol H2/mol glucose . However, the observed H2 yield obtained in the simulations and in the experiments never exceeded the Thauer limit; that is, the molar ratio of H2 produced to glucose consumed is always below 4. The high values of are an effect of the nonstoichiometric ratio between H2 and acetic acid seen during the experiments, described as in the model (equation 10). The ratios are often above 2 mol H2/mol acetic acid and, in the model, this comes out as a higher H2 yield coefficient parameter. The reason for the H2 and acetic acid ratio not being equal to 2, as the stoichiometry of reaction 1 shows, might be that some of the acetic acid produced is taken up and utilized by the cell.
Dissolved hydrogen: degree of supersaturation and its effect
Supersaturation of dissolved H2 occurred during fermentation, as measured experimentally and predicted by the model (Table 5), as reported previously [18, 19]. Oversaturation ranged from 4.9 to 52 times the saturation concentration, which corresponds well with the results of other studies [18, 19]. The large deviation in the measured values of H2aq (Table 5) is an effect of the large uncertainties in the measurements. The oversaturation of H2aq can also be clearly seen in Figure 9, where values of H2aq calculated with the model for a few experiments are shown together with the saturation concentration of H2 for experiment 6.
The experiment with the lowest stripping rate (experiment 4) showed the highest production of lactic acid (Figure 3). The other two experiments with low stripping rates (experiments 6 and 7) (Figures 5 and 6) showed only a slight increase in lactic acid production compared to the experiments with high stripping rates (experiments 1, 3 and 5) (Figures 1, 2 and 4) and very little lactic acid production compared to experiment 4. This can be explained by the lower growth rates resulting in lower H2 productivity, and hence lower values of H2aq, in experiments 5 through 7 (Figures 4, 5 and 6). Experiments 3 and 4 both showed high initial growth rates characterized by a high value of μmax and hence high H2 productivity, leading to higher values of H2aq. In experiment 4, with a low stripping rate (0.78 L/hour), H2aq reached a plateau faster than it did in the other experiments (Figure 9). At that stage, the culture reached maximum H2 productivity (data not shown), which resulted in linear rates of glucose consumption and H2 and acetic acid production (Figure 4). This behavior was reflected by the model (Figure 4) and can be explained by growth inhibition by H2. When H2 productivity reaches a critical value in the fermentation system, the cells respond with a reduction in the specific growth rate, and hence a plateau is observed. This plateau was reached much faster in experiment 4 than in the other two low-stripping-rate experiments (experiments 6 and 7), which was due to a combination of the higher initial growth rate and H2 productivity and the lower stripping rate (Table 1).
The model with estimated parameters for inhibition by H2aq can be used to indicate how H2 productivity, , H2aq and the overall mass transfer coefficient are related (Figure 10). The dissolved H2 concentration in experiment 4 reached almost 80% of the value of H2aqcrit (indicated by the red line in Figure 10) after 10 hours, and, Dissolved hydrogen: degree of supersaturation and its effect), H2 productivity and hence H2aq were constant until glucose was depleted (Figures 3 and 9). The dissolved H2 concentration in experiment 5, on the other hand, was always well above the red line in Figure 10 and also showed little lactic acid production. In experiments 6 and 7, H2aq also ended up near the red line (Figure 10); however, both reached it much later (at 17 hours; that is, only 4 hours before the glucose was depleted). This was due to lower growth rates and slightly higher stripping rates. Finally, H2aq in experiment 1 crossed the red line after 16 hours with little production of lactic acid, which may be due to an exponential increase in H2 productivity, and the line was crossed when very little glucose remained; hence almost no substrate was left for lactic acid production (Figure 1).
The organism thus seems to respond to an increase in H2aq by adjusting its growth rate and directing its metabolism toward the production of lactate, thus fine-tuning H2 productivity to prevent H2aq values detrimental to the cell. On the basis of both the experimental data and the model, it is clear that dissolved H2 has a negative effect on C. saccharolyticus and that H2 productivity and the stripping rate have a considerable impact on H2aq. This can explain why fermentation with high H2 productivity commonly results in low H2 yields [15, 21, 32, 33] and why an increase in the stripping rate tends to increase H2 yield and productivity .
Predicting the effect of dissolved hydrogen concentration
The dissolved H2 concentration can be lowered by stripping with N2 gas (Figure 11), which also results in higher H2 yields [19, 34–36]. Kraemer and Bagley  published details of a study in which they optimized the stripping rate with respect to H2 yield and productivity. They concluded that a stripping rate of 12 mL/minute, resulting in a volumetric mass transfer rate of 5/hour, was optimal. However, this was only valid for their particular system and not for systems with higher H2 productivity and different liquid-to-gas mass transfer properties. Figure 11 shows H2aq calculated by our model at various stripping rates and H2 productivities for the reactor system used herein. It is evident that H2aq is a function of both variables.
Furthermore, the current model cannot simulate a system without sparging, since k l a is based on the stripping rate (equation 24). However, Figure 10 provides an excellent basis for further discussions, since the relation between , and H2 productivity still holds, even without stripping. For instance, with C. saccharolyticus at a total pressure of 1 atm, no stripping, a productivity of 10 mmol H2/L/hour, and a gas composition of 60% H2 and 40% CO2, a value of of at least 7.5/hour is required to avoid serious H2 inhibition (Figure 10).
Kinetic parameters: inhibition by osmolarity
Experiments 2 and 8 were carried out with almost identical initial conditions, but they progressed quite differently, probably due mainly to different inocula and a slight difference in the initial glucose concentration. It was evident that H2aq as the sole inhibitor could not adequately describe the system at elevated substrate concentrations (Figures 7 and 8). Indeed, the model fitted the data from experiments 2 and 8 very well after introducing growth inhibition by osmotic pressure.
The average value of μmax in experiments 5 through 7, which were carried out using the same inocula as experiment 8, was used for experiment 8. However, since no experiments were carried out using the same inocula as experiment 2, together with a low glucose concentration, the value of μmax for experiment 2 is not known. Instead, the average value of μmax for the three sets of 5 g/L experiments (experiment 1, experiments 3 and 4, and experiments 5 through 7) was assumed (0.27/hour) when estimating the inhibition parameters for osmolarity, that is, RLacFOSM, OSMcrit and nμ (Table 6). Applying this extended model to the experiments with lower glucose concentrations confirmed that the osmolarity had a minor influence on the outcomes of the fermentations under those conditions. To estimate the three parameters describing the inhibition by osmolarity, μmax for experiment 2 had to be assumed. This assumption mainly affects the value of the exponential factor, nμ, and not RLacFOSM or OSMcrit (Table 5). This strongly indicates that OSMcrit for C. saccharolyticus is in the range of 0.27 to 0.29 mol/L, a range that is in agreement with previous results of a study by Willquist et al. .
Although this rather simple model describes the effect of osmolarity well, some aspects of it cannot be reproduced by this type of unstructured, mechanistic model. It is evident from Figure 8 that the growth rate, and hence the glucose consumption and acetic acid production rates, are not affected at the initiation of growth, since the model excluding inhibition by osmolarity fits the first three points perfectly. However, at an osmolarity of about 0.2 mol/L, there is a shift toward a reduction in the growth rate, glucose consumption, and acetic acid production rate, as well as an increase in lactic acid production, rapidly at first, followed by a slow decline of all of them. It is not possible to model this effect, that is, a rapid shift followed by a slow gradual decrease, using a simple inhibition expression. One could assume a very high exponential factor, but that would predict only the beginning of the experiment; however, as inhibition came into effect, the decrease in the growth rate to zero would not be slow and gradual, but very considerable and rapid.
Batch fermentations of C. saccharolyticus on glucose were successfully simulated using Monod kinetics extended to include liquid-to-gas mass transfer and inhibition by H2aq and osmolarity. In agreement with previous measurements [18, 19], the model predicted high oversaturation of H2 in the liquid, which was also confirmed experimentally.
To the best of our knowledge, we have demonstrated for the first time the possibility of predicting H2aq and its effect on the fermentation in a laboratory-scale bioreactor as a function of stripping rate and productivity using the derived model. With this tool, it was possible to conclude that the inhibition of growth and H2 production depend on both H2aq and osmolarity. The model described the inhibition by both H2aq and osmolarity satisfactorily and can hence be used to predict the optimal stripping rate and substrate concentration. It also shows how C. saccharolyticus responds to increased osmolarity and H2aq. In addition, the modeled relationship between mass transfer, H2 productivity and H2aq is not specific to C. saccharolyticus, but can also be applied to other systems. To widen the applicability of the model even further, future studies should aim at incorporating the effect of the stirrer speed and the absence of stripping on the liquid-to-gas mass transfer.
Ac: acetic acid concentration (mol/L); CO2aq: CO2 concentration in the liquid phase (mol/L); : saturation concentration (mol/L); CO2sol: solubilized CO2 (mol/L); : flow rate of CO2 out of the fermentor (L/hour); : flow rate of H2 out of the fermentor (L/hour); : flow rate of N2 into the fermentor (L/hour); : flow rate of N2 out of the fermentor (L/hour); G: substrate concentration (mol/L); H2aq: H2 concentration in the liquid phase (mol/L); H2aqcrit: critical H2 concentration in the liquid phase (mol/L); : saturation concentration (mol/L); : Henry's constant for H2 (L atm/mol); : Henry's constant for CO2 (L atm/mol); K1: dissociation constant for CO2 into bicarbonate; K2: dissociation constant for bicarbonate into carbonate; : mass transfer coefficient for H2 (per hour); : mass transfer coefficient for CO2 (per hour); KG: saturation constant (mol/L); Lac: lactic acid concentration (mol/L); nμ: inhibition coefficient for inhibition by osmolarity; : inhibition coefficient for inhibition by H2; OSM: osmolarity parameter (mol/L); OSMcrit: critical osmolarity parameter (mol/L); : partial pressure of H2 (Pa); : critical partial pressure of H2 (Pa); Ptot: total pressure (Pa); PPi: pyrophosphate; LSQ: least squares quadratic; R: gas constant (L atm/K/mol); RAcF: maximum Ac production rate (mol/g cell mass/hour); : maximum Lac production rate caused by H2 (mol/g cell mass/hour); RLacFOSM: maximum Lac production rate caused by osmotic pressure (mol/g cell mass/hour); rcd: death rate of cells (per hour); T: temperature (Kelvin); Vg: gas volume in liquid (L); Vl: liquid volume in reactor (L); X: cell mass concentration (mol/L); YGAc: stoichiometric acetic acid yield (mol acetic acid/mol substrate); : stoichiometric CO2 yield (mol CO2/mol substrate); : stoichiometric H2 yield (mol H2/mol substrate); YGLac: stoichiometric lactic acid yield (mol lactic acid/mol substrate); YGX: biomass yield (mol cell mass/mol substrate); α: cell-growth-associated constant; αAc: cell-growth-associated constant for Ac production rate (mol Ac/mol cells); : cell growth-associated constant for CO2 production rate (mol CO2/mol cells); : cell growth-associated constant for H2 production rate (mol H2/mol cells); αLac: cell-growth-associated constant for Lac production rate (mol Lac/mol cells); β: non-cell-growth-associated constant; γ: exponential coefficient for mass transfer correlation; μ: specific growth rate (per hour); μmax: maximum specific growth rate (per hour).
Hallenbeck PC, Benemann JR: Biological hydrogen production: fundamentals and limiting processes. Int J Hydrogen Energy. 2002, 27: 1185-1193. 10.1016/S0360-3199(02)00131-3.
Hallenbeck PC, Ghosh D: Advances in fermentative biohydrogen production: the way forward?. Trends Biotechnol. 2009, 27: 287-297. 10.1016/j.tibtech.2009.02.004.
Kraemer JT, Bagley DM: Improving the yield from fermentative hydrogen production. Biotechnol Lett. 2007, 29: 685-695. 10.1007/s10529-006-9299-9.
Chou CJ, Jenney FE, Adams MWW, Kelly RM: Hydrogenesis in hyperthermophilic microorganisms: implications for biofuels. Metab Eng. 2008, 10: 394-404. 10.1016/j.ymben.2008.06.007.
Kengen SWM, Goorissen HP, Verhaart M, Stams AJM, van Niel EWJ, Claassen PAM: Biological hydrogen production by anaerobic microorganisms. Biofuels. Edited by: Soetaert W, Vandamme EJ. 2009, Chichester, UK: John Wiley & Sons, 197-221.
van Groenestijn JW, Hazewinkel JHO, Nienoord M, Bussmann PJT: Energy aspects of biological hydrogen production in high rate bioreactors operated in the thermophilic temperature range. Int J Hydrogen Energy. 2002, 27: 1141-1147. 10.1016/S0360-3199(02)00096-4.
Ljunggren M, Zacchi G: Techno-economic evaluation of a two-step biological process for hydrogen production. Biotechnol Prog. 2010, 26: 496-504.
de Vrije T, Bakker R, Budde MAW, Lai MH, Mars AE, Claassen PAM: Efficient hydrogen production from lignocellulosic energy crop Miscanthus by the extreme thermophilic bacteria Caldicellulosiruptor saccharolyticus and Thermotoga neapolitana. Biotechnol Biofuels. 2009, 2: 12-10.1186/1754-6834-2-12.
van Niel EWJ, Budde MAW, de Haas GG, van der Wal FJ, Claassen PAM, Stams AJM: Distinctive properties of high hydrogen producing extreme thermophiles, Caldicellulosiruptor saccharolyticus and Thermotoga elfii. Int J Hydrogen Energy. 2002, 27: 1391-1398. 10.1016/S0360-3199(02)00115-5.
van der Werken HJG, Verhaart MRA, VanFossen AL, Willquist K, Lewis DL, Nichols JD, Goorissen HP, Mongodin EF, Nelson KE, van Niel EWJ, Stams AJM, Ward DE, de Vos WM, van der Oost J, Kelly RM, Kengen SWM: Hydrogenomics of the extremely thermophilic bacterium Caldicellulosiruptor saccharolyticus. Appl Environ Microbiol. 2008, 74: 6720-6729. 10.1128/AEM.00968-08.
Ivanova G, Rakhely G, Kovacs KL: Hydrogen production from biopolymers by Caldicellulosiruptor saccharolyticus and stabilization of the system by immobilization. Int J Hydrogen Energy. 2008, 33: 6953-6961. 10.1016/j.ijhydene.2008.08.058.
Kádár Z, de Vrije T, Budde MAW, Szengyel Z, Réczey K, Claassen PAM: Hydrogen production from paper sludge hydrolysate. Appl Biochem Biotechnol. 2003, 105-108: 557-566.
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. 10.1002/bit.10463.
Willquist K, Claassen PAM, van Niel EWJ: Evaluation of the influence of CO2 on hydrogen production by Caldicellulosiruptor saccharolyticus. Int J Hydrogen Energy. 2009, 34: 4718-4726. 10.1016/j.ijhydene.2009.03.056.
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.
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. 10.1099/mic.0.27159-0.
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.
Pauss A, Andre G, Perrier M, Guiot SR: Liquid-to-gas mass transfer in anaerobic processes: inevitable transfer limitations of methane and hydrogen in the biomethanation process. Appl Environ Microbiol. 1990, 56: 1636-1644.
Kraemer JT, Bagley DM: Supersaturation of dissolved H2 and CO2 during fermentative hydrogen production with N2 sparging. Biotechnol Lett. 2006, 28: 1485-1491. 10.1007/s10529-006-9114-7.
Kraemer JT, Bagley DM: Optimisation and design of nitrogen-sparged fermentative hydrogen production. Int J Hydrogen Energy. 2008, 33: 6558-6565. 10.1016/j.ijhydene.2008.08.033.
de Vrije T, Mars AE, Budde MAW, Lai MH, Dijkema C, de Waard P, Claassen PAM: Glycolytic pathway and hydrogen yield studies of the extreme thermophile Caldicellulosiruptor saccharolyticus. Appl Microbiol Biotechnol. 2007, 74: 1358-1367. 10.1007/s00253-006-0783-x.
Han K, Levenspiel O: Extended Monod kinetics for substrate, product, and cell inhibition. Biotechnol Bioeng. 1988, 32: 430-437. 10.1002/bit.260320404.
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. 10.1016/j.ymben.2010.01.001.
Luedeking R, Piret EL: A kinetic study of the lactic acid fermentation. J Biochem Microbiol Technol Eng. 1959, 1: 393-412. 10.1002/jbmte.390010406.
Willquist K: Physiology of Caldicellulosiruptor saccharolyticus: a hydrogen cell factory. PhD thesis. 2010, Lund University, Department of Applied Microbiology
Hill GA: Measurement of overall volumetric mass transfer coefficients for carbon dioxide in a well-mixed reactor using a pH probe. Ind Eng Chem Res. 2006, 45: 5796-5800. 10.1021/ie060242t.
Linton KJ, Higgins CF: Structure and function of ABC transporters: the ATP switch provides flexible control. Pflügers Arch Eur J Physiol. 2007, 453: 555-567. 10.1007/s00424-006-0126-x.
Chevance FF, Erhardt M, Lengsfeld C, Lee SJ, Boos W: Mlc of Thermus thermophilus: a glucose-specific regulator for a glucose/mannose ABC transporter in the absence of the phosphotransferase system. J Bacteriol. 2006, 188: 6561-6571. 10.1128/JB.00715-06.
Nielsen JH, Villadsen J, Lidén G: Bioreaction Engineering Principles. 2003, New York: Kluwer Academic/Plenum Publishers, 2
Koch AL: Diffussion: the crucial process in many aspects of the biology of bacteria. Adv Microb Ecol. 1990, 2: 37-70.
Thauer RK, Jungermann K, Decker K: Energy conservation in chemotrophic anaerobic bacteria. Bacteriol Rev. 1977, 41: 100-180.
Zeidan AA, van Niel EWJ: Developing a thermophilic hydrogen-producing co-culture for efficient utilization of mixed sugars. Int J Hydrogen Energy. 2009, 34: 4524-4528. 10.1016/j.ijhydene.2008.07.092.
Koskinen PEP, Lay C, Puhakka JA, Lin P, Wu S, Örlygsson J, Lin CY: High-efficiency hydrogen production by an anaerobic thermophilic culture from an Icelandic hot spring. Biotechnol Bioeng. 2008, 101: 665-678. 10.1002/bit.21948.
Hussy I, Hawkes FR, Dinsdale R, Hawkes DL: Continuous fermentative hydrogen production from a wheat starch co-product by mixed microflora. Biotechnol Bioeng. 2003, 84: 619-626. 10.1002/bit.10785.
Mizuno O, Dinsdale R, Hawkes FR, Hawkes DL, Noike T: Enhancement of hydrogen production from glucose by nitrogen gas sparging. Bioresour Technol. 2000, 73: 59-65. 10.1016/S0960-8524(99)00130-3.
Kyazze G, Martinez-Perez N, Dinsdale R, Premier GC, Hawkes FR, Guwy AJ, Hawkes DL: Influence of substrate concentration on the stability and yield of continuous biohydrogen production. Biotechnol Bioeng. 2006, 93: 971-979. 10.1002/bit.20802.
The authors gratefully acknowledge the Commission of the European Communities, Sixth Framework Programme, Priority 6, Sustainable Energy Systems (019825 HYVOLUTION), as well as the Swedish Energy Agency for their financial support.
The authors declare that they have no competing interests.
ML developed the model, carried out the parameter estimations and planned as well as performed the mass transfer experiments. ML also wrote the main part of the manuscript and took part in the planning and execution of the fermentation experiments. KW took part in the development of the model, planned and carried out the main part of the fermentation experiments, analyzed the results and assisted in the mass transfer experiments. KW also wrote parts of the manuscript. GZ and EVN participated in the coordination of the study and reviewed the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Liquid-to-gas mass transfer and chemistry of CO2__Theory_experiments_and_results. This file contains information on the estimation of the mass transfer parameters as well as the set-up of the mass transfer experiments. The file also contains more detailed information concerning the theory behind the chemistry of carbon dioxide and liquid-to-gas mass transfer. (DOCX 40 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Ljunggren, M., Willquist, K., Zacchi, G. et al. A kinetic model for quantitative evaluation of the effect of hydrogen and osmolarity on hydrogen production by Caldicellulosiruptor saccharolyticus. Biotechnol Biofuels 4, 31 (2011). https://doi.org/10.1186/1754-6834-4-31