A non-linear model of hydrogen production by Caldicellulosiruptor saccharolyticus for diauxic-like consumption of lignocellulosic sugar mixtures

Background Caldicellulosiruptor saccharolyticus is an attractive hydrogen producer suitable for growth on various lignocellulosic substrates. The aim of this study was to quantify uptake of pentose and hexose monosaccharides in an industrial substrate and to present a kinetic growth model of C. saccharolyticus that includes sugar uptake on defined and industrial media. The model is based on Monod and Hill kinetics extended with gas-to-liquid mass transfer and a cybernetic approach to describe diauxic-like growth. Results Mathematical expressions were developed to describe hydrogen production by C. saccharolyticus consuming glucose, xylose, and arabinose. The model parameters were calibrated against batch fermentation data. The experimental data included four different cases: glucose, xylose, sugar mixture, and wheat straw hydrolysate (WSH) fermentations. The fermentations were performed without yeast extract. The substrate uptake rate of C. saccharolyticus on single sugar-defined media was higher on glucose compared to xylose. In contrast, in the defined sugar mixture and WSH, the pentoses were consumed faster than glucose. Subsequently, the cultures entered a lag phase when all pentoses were consumed after which glucose uptake rate increased. This phenomenon suggested a diauxic-like behavior as was deduced from the successive appearance of two peaks in the hydrogen and carbon dioxide productivity. The observation could be described with a modified diauxic model including a second enzyme system with a higher affinity for glucose being expressed when pentose saccharides are consumed. This behavior was more pronounced when WSH was used as substrate. Conclusions The previously observed co-consumption of glucose and pentoses with a preference for the latter was herein confirmed. However, once all pentoses were consumed, C. saccharolyticus most probably expressed another uptake system to account for the observed increased glucose uptake rate. This phenomenon could be quantitatively captured in a kinetic model of the entire diauxic-like growth process. Moreover, the observation indicates a regulation system that has fundamental research relevance, since pentose and glucose uptake in C. saccharolyticus has only been described with ABC transporters, whereas previously reported diauxic growth phenomena have been correlated mainly to PTS systems for sugar uptake.


Background
The need for renewable energy is ever increasing to tackle the major challenges of global warming, energy demand, and limited resources. According to statistics published by the International Energy Agency [1], just over 86% of the Total Primary Energy Supply (TPES) in 2014 was produced from fossil resources, leaving a modest 14% originating from renewable energy sources. When putting these numbers in relation with the adopted Paris Agreement in 2015, targeting to keep the global average temperature increase below the 2 °C above pre-industrial levels [2], it is evident that actions need to be taken. There are, however, positive trends in that the supply of renewable energy sources has grown faster, with an average annual rate of 2.0% since 1990, compared to the growth of the world TPES of 1.8% [1].
Hydrogen has the potential of becoming an important renewable energy carrier. Currently, hydrogen is widely used as a reducing agent in the chemical and food industry. However, using hydrogen as an energy carrier in sustainable applications is of great interest due to its potentially high efficiency of conversion to usable power, its low emissions of pollutants and high energy density [3]. Up to 96% of the world's hydrogen production is fossil based, i.e., natural gas, oil, and coal [4]. A sustainable alternative to the conventional methods for producing hydrogen is by biological methods, i.e., biohydrogen. There are four major categories in which production of biological hydrogen can be classified, namely: photofermentation of organic compounds by photosynthetic bacteria, biophotolysis of water using algae and cyanobacteria, bioelectrohydrogenesis, and fermentative hydrogen production, so-called dark fermentation, from organic wastes or energy crops [5,6]. The latter is the focus of this study, where various sugars present in, e.g., agricultural waste like wheat straw, can be fermented by microorganisms for hydrogen production. This also addresses the challenge of converting lignocellulosic biomass to renewable energy.
The thermophilic microorganism Caldicellulosiruptor saccharolyticus is able to produce hydrogen from lignocellulosic biomass through dark fermentation and has previously shown the potential of producing hydrogen close to the maximum theoretical yield of 4 mol hydrogen per mol hexose [9][10][11]. C. saccharolyticus is cellulolytic and can utilize a broad range of di-and monosaccharides for hydrogen production [12]. Van de Werken et al. [13] showed that C. saccharolyticus coferments glucose and xylose as it lacks catabolite repression. VanFossen et al. [14] revealed that although C. saccharolyticus co-utilizes different sugars, it has a preference for some sugars over others. Xylose was discussed as a preferred sugar over glucose and is, therefore, utilized by the microorganism to a greater extent than glucose. However, the substrate uptake kinetics was not determined and a yeast extract (YE)-supplemented medium was used [13].
By developing a mathematical model for a biological process, it is possible to describe past and predict future performances as well as gaining a deeper understanding of the physiological mechanism behind the process. The aim of this study is to present a model that describes the growth of C. saccharolyticus on lignocellulosic sugar mixtures and how the uptake rate changes when the sugars are used simultaneously or individually. Similar kinds of models have been developed [15,16]; however, these models focus on single sugar uptake. The proposed model here builds on the one presented by Ljunggren et al. [15] by adding the consumption rates for each individual sugar in the sugar mixtures. Monod [17] first described the phenomenon of diauxic growth, where a microorganism is exposed to two substrates and first consumes the substrate that supports the most efficient growth rate. Several models have been developed in this area [18,19] describing how to capture the subsequent uptake of sugars when multiple sugars are present. This phenomenon can be modeled using a cybernetic approach to whether a particular enzyme, needed for a specific sugar to be metabolized, is upregulated or not. This paper describes the development of a substratebased uptake model using Monod-type kinetics including biomass growth, product formation, liquid-to-gas mass transfer, and enzyme synthesis with Hill kinetics, with C. saccharolyticus as model organism. The model presented in this paper takes into consideration the usage of different sugars, including hexoses, i.e., glucose, and pentoses, i.e., xylose and arabinose. The model describes the different sugar uptakes individually, exemplifying the rate at which each sugar is consumed when C. saccharolyticus grows on the sugar mixtures and on the individual sugars, respectively.

Strains and cultivation medium
Caldicellulosiruptor saccharolyticus DSM 8903 was obtained from the Deutsche Sammlung von Mikroorganismen und Zellkulturen (Braunschweig, Germany). Subcultivations were conducted in 250 mL serum flasks with 50 mL modified DSM 640 media [20]. The carbon source of each cultivation corresponded to that of the subsequent fermentor cultivation. The 1000× vitamin solution and modified SL-10 solution were prepared according to [20] and [21], respectively.
All bioreactor experiments used a modified DSM 640 medium with the exclusion of yeast extract according to Willquist and van Niel [20]. To quantify the kinetics of xylose and glucose uptake and the effect of when the sugars were mixed in pure and industrial medium, the growth and hydrogen production was monitored in four different cases, where the total sugar concentration in the medium was fixed to 10 g/L. Cultivations were performed using 10 g/L glucose (Case 1), 10 g/L xylose (Case 2), a sugar mixture (Case 3), and wheat straw hydrolysate (Case 4). In Case 4, a 9% solution of wheat straw hydrolysate was used corresponding to approximately 10 g/L total sugars. In Case 3, the sugar mixture contained pure sugars with the same concentration as the wheat straw hydrolysate (6.75 g/L glucose, 3.06 g/L xylose, and 0.173 g/L arabinose). The total sugar concentrations at the start of the fermentation included the sugar added as described above and the additional sugar added from the inoculum, which varied slightly in the different conditions. The starting sugar concentration was, therefore, as follows: Case 1, 12.11 ± 0.09 g/L glucose; Case 2, 10.96 ± 0.20 g/L xylose; Case 3, 8.69 ± 0.12 g/L glucose, 3.38 ± 0.19 g/L xylose, and 0.38 ± 0.01 g/L arabinose; Case 4, 7.31 ± 0.07 g/L glucose, 3.36 ± 0.06 g/L xylose, and 0.34 ± 0.00 g/L arabinose.

Fermentor setup
Batch cultivations were performed in a jacketed, 3-L fermentor equipped with ADI 1025 Bio-Console and ADI 1010 Bio-Controller (Applikon, Schiedam, The Netherlands). A working volume of 1 L was used for cultivations and the pH was maintained at optimal conditions 6.5 ± 0.1 at 70 °C by automatic titration with 4 M NaOH. The temperature was thermostatically kept at 70 ± 1 °C. Stirring was maintained at 250 rpm and nitrogen was sparged through the medium at a rate of 6 L/h. Sparging was initiated 4 h after inoculation and was continued throughout the cultivation. A condenser cooled with water at 4 °C was utilized to prevent evaporation of the medium. Samples were collected at regular time intervals for monitoring of the optical density. The supernatant from each culture was collected and stored at − 20 °C for further quantification of various sugars and organic acids. Gas samples were collected from the fermentor's headspace to quantify H 2 and CO 2 . The sugar mixture and wheat straw hydrolysate experiments were done in triplicate. The individual sugar fermentations were done in biological duplicate.
A defined medium was autoclaved in each fermentor, while anoxic solutions of cysteine HCl·H 2 O (1 g/L), MgCl 2 ·6H 2 O (0.4 g/L), and carbon source(s) were prepared separately and were added to the fermentor before inoculation. Just after inoculation, the fermentor was closed for 4 h to allow buildup of CO 2 as previously described [20] necessary to initiate growth.

Analytical methods
Optical density was determined using an Ultraspec 2100 pro spectrophotometer (Amersham Biosciences) at 620 nm. Sugars, organic acids, hydroxymethyl furfural (HMF), and furfural were detected using HPLC (Waters, Milford, MA, USA). For the quantification of organic acids, an HPLC equipped with an Aminex HPX-87H ion-exchange column (Bio-Rad, Hercules, USA) at 60 °C and 5 mM H 2 SO 4 as mobile phase was used at a flow rate of 0.6 mL/min. Glucose, xylose, and arabinose quantification was conducted using an HPLC with a Shodex SP-0810 Column (Shodex, Japan) with water as a mobile phase at a flow rate of 0.6 mL/min. CO 2 and H 2 were quantified with a dual channel Micro-GC (CP-4900; Varian, Micro-gas chromatography, Middelburg, The Netherlands), as previously described [21].

Mathematical model description
The model developed for C. saccharolyticus in this study takes into account the kinetics of biomass growth, consumption of glucose, xylose and arabinose, and formation of the products acetate, hydrogen, and carbon dioxide. Furthermore, the model includes liquid-to-gas mass transfer of hydrogen and carbon dioxide as well as the equilibrium between carbon dioxide, bicarbonate (HCO 3 − ) and carbonate (CO 3 2− ). The model is developed on a cmol basis. The formation of lactate was excluded to reduce the complexity of the model, as it constituted to less than 5% of the total product in the sugar mixture fermentations. In addition, inhibition due to high aqueous H 2 concentration and high osmolarity was not included in the model to reduce the number of unknown parameters. This is motivated by the fact that the focus of this study is mainly on the consumption behavior of C. saccharolyticus on the different sugars.
The model is constructed with a similar nomenclature and setup as in the anaerobic digestion model no 1 (ADM1) described by Batstone et al. [22] and was implemented in MATLAB R2015b (Mathworks, USA). The following biochemical degradation reactions are the basis for the model (Eqs. 1, 2).
Biomass formation from sugar [23]: Reaction 1 is not balanced, since there were elements in the fermentation medium that were not included in the model, i.e., cysteine. The value of the yield factor Y X is calculated from the data of the batch fermentations. It is assumed that nitrogen, sulfur, and phosphorus are in excess in the media and, therefore, are not included as separate entities in the mathematical model. Sugar degradation to product formation by C. saccharolyticus in cmol:

Model inputs and initial conditions
The model requires a range of input variables. The lag time was determined by calculating the intersection point between the lag phase and the exponential phase when taking the natural logarithm of the biomass concentration over time, as illustrated by Swinnen et al. [24]. Since the lag phase is dependent on the culture status before the fermentation, which was not addressed in this study, it was excluded from the experimental data when the latter were compared to model data and for initial input values for the model. The start values of the unknown state variables are listed in Table 1. The constants used in the model are presented in Table 2.

Mass balances for biomass growth, substrate consumption, and product formation in the liquid phase
The stoichiometric relationships and mass balances of the reactants and products present in the model are displayed in Table 3. The model is supplemented with an enzyme, E2, and cybernetic variables v and u as in [18], where the former controls the activity of the enzyme and the latter is the fractional allocation of a critical resource for the synthesis of the enzyme. We hypothesize that initially, there is a first enzyme system present aiding the subsequent uptake of both hexose and pentose sugars, but with a preference for the pentoses (phase I). This transporter is only available as long as pentoses are present. After depletion of the pentoses, a second enzyme system, E2, is synthesized allowing for uptake of the remaining hexose sugars by a second transporter (phase II). For the sake of convenience, we simplify the enzyme system, consisting of multiple proteins, using the word enzyme and using this abstraction also in the kinetic model. The mass balance for the biomass, X, is dependent on the rate of substrate consumption ρ, with Monodtype kinetics, and on the biomass decay rate, which is described with first-order kinetics, where r cd (h −1 ) is the cell death rate and Y x (cmol/cmol) is the yield of biomass from total sugar (Table 3). A second glucose rate equation ( ρ Glu, 2 ) is added to describe the diauxic-like growth appearance in the sugar mixture. The rate of the glucose consumption, when the pentose sugars are depleted, is dependent on enzyme E2. The rate of the enzyme synthesis, ρ E , is based on Hill kinetics, as in [19], the decay rate of the enzyme is first-order kinetics, and the third term, − 1·E2·ρ Glu, 2 , represents the dilution of the specific enzyme level as is described with kinetics similar to Hill, i.e., E2 2 . The parameters k m and k m,2 (h −1 ) are the maximal uptake rates in phase I and phase II, respectively, and K s,glu , K s,glu,2 , K s,xyl , K s,ara, and K s,E2 (cmol/L) are the affinity constants for the uptake of glucose, xylose, arabinose, and synthesis of enzyme E2, respectively. Finally, α is the enzyme synthesis rate (h −1 ) and β is the enzyme decay rate (h −1 ).
Acetate, hydrogen, and carbon dioxide are produced in the liquid phase. Y ac (cmol/cmol), Y H 2 (mol/cmol) and Y CO 2 (cmol/cmol) represent the conversion yields of acetate, hydrogen, and carbon dioxide, respectively, from both hexose and pentose sugars. The conversion yields were fitted with experimental data from the batch fermentations. Y X was determined by the slope of the curve: total sugar vs biomass; here, only phase I was considered. Y ac and Y CO 2 were determined by first taking the slope of the curves, total sugar vs acetate, and total sugar vs carbon dioxide, and then, the actual yields were calculated according to the following equation: When Y H 2 was calculated the same way as in Eq. 3, it gave a too high conversion yield. To obtain a more accurate yield, the effects of liquid-to-gas mass transport were considered and Y H 2 was instead determined as follows: H 2,end − H 2,start Tot sugar start − Tot sugar end .

Acid-base reactions
The acid-base reaction considered in the model is that of carbon dioxide, bicarbonate, and carbonate formation. ρ AB,CO 2 in Table 4 describes the rate of formation of bicarbonate and carbonate. CO 2,sol is the sum of the ionic species, HCO − 3 and CO 3 2− and Eq. 5 gives the differential equation for CO 2,sol :

Liquid-to-gas mass transfer and mass balances for product formation
Hydrogen and carbon dioxide are produced in the liquid phase and then transferred to the gas phase via liquid-togas mass transport. ρ t,H 2 describes the mass transfer rate of hydrogen and ρ t,CO 2 is the mass transfer rate of carbon dioxide (Table 5). p gas,H 2 and p gas,CO 2 (in atm then converted to Pa) are the partial pressures of H 2 and CO 2 , respectively.
The expression for the mass balances describing the gaseous products can be described as in Eqs. 6, 7, where q gas (L/h) is the total gas flow, and V liq and V gas (L) are the liquid and the gas volumes, respectively:

Sensitivity analysis
A sensitivity analysis can identify parameters that have great effect on the model output. The sensitivity analysis was done based on the OFAT approach, i.e., one-factorat-at-time [25]. The chosen parameter was altered with a factor δ, as described in [26], to see the effect on the different state variable output result, as in the following equation: where Γ i,j is the sensitivity of state variable i with respect to model parameter j in each timepoint of the Matlab simulation. Furthermore, y i (θ j ) is the value of state variable i in regard to parameter j and y i θ j + δ · θ j is the value of state variable i when parameter j has been altered with a factor δ. The parameters that were included in the sensitivity analysis were k m , k m,2 , K s,glu , K s,glu,2 , K s,xyl , K s,ara , K s,E2 , α, n, r cd, and k L a H2 and the state variables that were considered were Glu, Xyl, Ara, Ac, X, and H 2 . The presented sensitivity data of one parameter in regards to a specific state variable were calculated as the average of Γ i,j .

Model calibration
To get a better fit to the experimental data, the model parameters were calibrated using the knowledge that was revealed in the sensitivity analysis. This was done with the function lsqcurvefit in MATLAB which uses a least square method to find the right parameter value for a non-linear curve fitting by seeking to find coefficients x that solve the problem in the following equation: given the input data xdata and the observed output ydata, where xdata and ydata are matrices or vectors and CO 2 acid-base 1 − 1 ρ AB,CO 2 = k AB · (CO 2,aq · 10 −pK 1 10 −pH + 10 −pK1 · 10 −pK 2 (10 −pH ) 2 − CO 2,sol F(x,xdata) is a matrix-valued or vector-valued function of the same size as ydata.
The lsqcurvefit function starts at x0 and finds coefficient, i.e., parameter x, to best fit the non-linear function fun(x,xdata) to the data ydata: The uncertainties of the calibrated parameters were assessed by calculating the confidence interval. This was done with the function nlparci in MATLAB which computes the 95% confidence intervals for the non-linear least square parameters estimated.

Growth profiles on the various sugars
The growth profiles of the single sugar experiments (glucose; Case 1 and xylose; Case 2), sugar mixture experiments (Case 3) and wheat straw hydrolysate experiments (Case 4) are presented in Fig. 1a-d. Glucose is consumed approx. two times faster when used as sole substrate (10) x = lsqcurvefit(fun, x0, xdata, ydata).
(Case 1) than in the sugar mixtures (Cases 3 and 4). Xylose, on the other hand, is consumed approx. two times slower when used as sole substrate and is completely consumed after approx. 60 h compared to around 20 h when co-fermented with other sugars (Cases 3 and 4; Fig. 1c, d). The highest production rate of acetate and hydrogen occurred around 20 h both in the sugar mixture and in the wheat straw hydrolysate fermentations. Lactate was formed just after 20 (Case 3) and 30 h (Case 4) reaching in total 0.015 and 0.014 cmol/L, respectively.
The calculated lag phases differed for each experiment. The lag phases of the sugar mixture experiments ranged from 9 to 11 h, whereas the lag phase of the wheat straw hydrolysate experiment was 4 h. This observation could be correlated to the richer nutrient content of wheat straw than the defined sugar mixture medium. A similar observation was found by Pawar et al. [27]. The lag phase with glucose alone was 4 h, but there was no lag phase with xylose alone. It is worth noticing though that it took more effort to initiate  1, 3, and 4) failed. This is due to that precautions are needed to start a culture on xylose in the absence of yeast extract, such as no sparging for several hours. The profiles of the mixed sugars indicate a biphasic growth, where the uptake of glucose decreased after xylose was depleted, but then increased again (Fig. 1c,  d). The two-phased sugar uptake was more pronounced in the wheat straw hydrolysate fermentations. The behavior can be further illustrated by the hydrogen productivity and CO 2 productivity (Fig. 2a, b). This observation has, to our knowledge, not been reported for Caldicellulosiruptor previously, although the transcriptomics of multiple sugar uptake have been extensively studied [13,14]. One possible reason for this could be that many multi-sugar experimental studies on this genus have been performed on a yeast extractsupplemented medium [3]. Because yeast extract itself partly supports growth [20], it possibly masks biphasic behavior. Moreover, the initial ratio of pentose/hexose sugars was higher in those studies [14] than in the WSH used in this study. Thus, after xylose was consumed, the culture adapted to a hexose-only medium, which initiated a second phase of growth.
The emerging pattern resembles a diauxic growth behavior, which was first described by Monod [17], and is characterized by two growth phases often separated with a lag period. This normally occurs in the presence of two carbon sources, where the preferred one is consumed first by the microorganism followed by the second after a lag period [28][29][30]. However, in the case of C. saccharolyticus, both pentose and hexose sugars are consumed simultaneously, albeit with a slight preference for the former. When the pentose sugars are depleted hexose consumption continues, but in Case 4 that happened with an increased rate (Table 8).
To quantify this behavior and investigate whether the theory of diauxic growth could be used to explain the observations, a kinetic model was developed consisting of two phases. In the phase I, glucose was consumed simultaneously with xylose and arabinose. Van de Werken et al. [13] concluded that growth on glucose and xylose mixtures as well as growth on the individual sugars all trigger transcription of the genes encoding a xylose-specific ABC transport system. This supports our hypothesis that glucose, xylose, and arabinose were initially transported by the same uptake system. However, when xylose was depleted, phase II starts with a new uptake system being expressed that had a higher affinity for glucose, transporting glucose at an altered rate. It is relevant to observe, however, that diauxic growth behavior is generally considered to be related to PTS systems [31][32][33]. However, according to current knowledge, C. saccharolyticus only possesses ABC transport systems [13,14]. Still, it has been described that other transport systems can generate this diauxic growth profile. For example, in Streptomyces coelicolor and related species, the genes involved in carbon catabolite repression are PTS independent, and instead, glucose kinase is the main controlling enzyme [33].

Determination of conversion yields
The calculated conversion yields from the batch experiments differ from the stoichiometric yields (Table 6). To begin with, the single sugar fermentations the calculated yields are lower than the corresponding stoichiometric yields. This is in contrast to the yields calculated for the sugar mixture experiments, except for Y ac that was slightly lower. The lower yield for acetate could be due to that part of the acetate, or rather acetyl-CoA, which is used as a building block for cell mass production [34].
The carbon balances attained in the model were 90 and 102% with start data from the sugar mixture experiments  Fig. 2 a Hydrogen productivity and b CO 2 productivity in Cases 3 and 4, sugar mixture fermentation, and wheat straw hydrolysate fermentation, respectively and the WSH experiments, respectively, which are equal or close to the values calculated from the experimental data, 90 and 107%, respectively, Table 6. The higher values in the carbon balance, i.e., > 100%, for the WSH fermentations, could be due to that other carbon sources may be present, such as oligosaccharides, that are also converted to products giving a higher carbon and electron output.

Sensitivity analysis
Dynamic simulations using benchmark parameter values [15] showed discrepancies between the experimental results and the model predictions. To further improve the dynamic simulations, a sensitivity analysis was conducted to determine the most important parameters. This was done with start values both from the sugar mixture fermentations as well as from the wheat straw hydrolysate fermentations. The change, δ, in the parameter value was set to 1% as in [35]. The sensitivity analysis allowed ranking of the parameters, which was useful for the model calibration. The most sensitive parameters, i.e., with a sensitivity value of > 1%, in regard to each of the state variables are listed in Table 7. The state variables that were affected the most by a change in parameter value were Glu and Xyl. The sensitivities of the other parameters for the different state variables were less than 1%.

Parameter calibration
The sensitivity analysis served as a basis for the parameter calibration. The model was calibrated with data from the four different batch experiments, Cases 1-4. Start values of the state variables were taken from the experimental data (Table 1), and initial parameter values, i.e., benchmark values, were taken from the literature [15] or guesstimated, e.g., by manually fitting the curves of the data points. The calibrated parameters together with a confidence interval of 95% are given in Table 8. Some of the parameters were graphically calibrated and, therefore, are without a confidence interval. The simulations with start data from the single glucose and xylose fermentations were carried out without the diauxic-like growth additions; thus, only phase I was applied.
The k m values for Cases 3 and 4 describe the maximal simultaneous uptake rates of glucose, xylose, and arabinose (Table 8), and they are modeled with the same value for all the sugars in phase I. However, the K s values for glucose in phase I, K s,glu , are higher than the K s values for xylose, K s,xyl , which indicates a lower affinity for glucose in phase I, since xylose is present and preferred. Moreover, K s,glu in Case 4 is 18 times higher compared to K s,glu,2 and compared to K s,glu in Case 3. One explanation is the greater affinity for xylose in phase I and another possible explanation is that K s,glu in Case 4 also includes an inhibition term due to the characteristics of the wheat straw hydrolysate media, e.g., Eq. 11: (11) K s,glu = K s,glu, real · I,  Glu k m,2 , α, k m , r cd , k L a H2 , K s,glu,2 k L a H2 , α, k m,2 , r cd , k m , K s,glu,2 Xyl k L a H2 , k m , K s,ara , K s,xyl k L a H2 , k m , K s,xyl , K s,ara , K s,E2 Ara K s,ara , k m , k L a H2 K s,ara Ac -- where I represents a competitive inhibition, Eq. 12: with S I the concentration of the inhibitor and K I the inhibition parameter. This is possibly due to unknown inhibiting compounds in the wheat straw hydrolysate or other factors that inhibit glucose uptake in phase I in Case 4. The reason behind the competitive inhibition has not been identified, but we hypothesize the presence of oligosaccharides that might be preferably taken up instead of glucose. However, these sugars were not quantified in the HPLC analysis of WSH. The k m,2 value for Case 4 is 50% lower than the corresponding value for the glucose uptake rate in Case 1. One explanation for this is that the enzymes involved in the sugar uptake in Case 4 take some time to be synthesized making glucose consumption slower in the WSH compared to the single glucose fermentation. Again, the presence of inhibiting compounds or competitive oligosaccharides could further slow down the glucose uptake rate.
Furthermore, the results show that on single sugars and mineral medium, glucose uptake is approximately (12) I = 1 + S I K I 35% faster than xylose uptake (Table 8). Moreover, growth of C. saccharolyticus on glucose is approx. 40% faster than on xylose (Table 9). This outcome contradicts the previous results on these two sugars in media supplemented with yeast extract (YE), where growth is faster on xylose than on glucose [13,14]. An explanation for this observation could be that C. saccharolyticus needs other sugars (present in YE) to grow optimal on xylose. Indeed, when both sugars are present the growth on xylose is stimulated by the co-uptake of glucose. The stoichiometric relationship of glucose-toxylose uptake rate ρ(Glucose):ρ(Xylose) was affected by the media used and is approximately 0.7 and 0.3 in phase I for growth on defined sugar mixture and wheat  straw medium, respectively (data used from Fig. 1). Until xylose is depleted, the total glucose, xylose, and arabinose conversion rates, i.e., 0.54·3 h −1 , are similar to that of xylose conversion in the absence of glucose, i.e., 1.58 h −1 . This observation is supported by other studies with C. saccharolyticus using different sugar mixtures both with and without YE, e.g. in Willquist [36]. Xylose uptake increases if a small concentration of glucose is present or if either the fermentor is sparged with CO 2 instead of N 2 gas or closed, to allow buildup of HCO 3 − in the reactor.

Model prediction
Comparison between the model and experimental results for the combined sugars is depicted in  Table 10 shows the fitting between the experimental data and the model simulation displaying the regression analysis values. It is clear that the model is well able to describe the consumption of the different sugars as well as biomass growth, acetate formation, and accumulation of hydrogen in Cases 3 and 4. The model, without the diauxic-like additions, was better at describing the individual xylose fermentations (Case 2), rather than the individual glucose fermentations (Case 1) when it comes to biomass growth and hydrogen production ( Table 10).
The model only predicts a small second peak in hydrogen productivity compared to the data of the defined sugar mixture fermentations (Fig. 3g). However, the model succeeds in describing the diauxic-like behavior of the hydrogen productivity profile in the wheat straw hydrolysate fermentations (Fig. 4g). The uptake of the three sugars as well as the formation of acetate is well described by the model, both for Cases 3 and 4 (Figs. 3a-d, 4a-d).
According to the simulation, the enzyme (used to describe the diauxic behavior) concentration is very low, close to zero, in the beginning, and when phase I ends, the enzyme synthesis starts and the concentration increases up to a peak, where it begins decreasing just before t = 60 h in Case 3 and somewhat earlier in Case 4 (Figs. 3f, 4f ). The enzyme synthesis is dependent on the biomass concentration, which is why it follows the behavior of the latter. The two biomass growth phases are clearly displayed in Case 4 and expressed by the model (Fig. 4e), where a first growth phase takes place between 0 and 20 h and a second growth phase between 20 and 45 h. The phenomenon with two growth phases is characteristic for diauxic growth behavior as described in various literatures on the topic [18,28,37].
The hydrogen productivity profile, both in Cases 3 and 4, is a bit delayed in the model (Figs. 3g, 4g). This could be due to a slight underestimation of the k L a H 2 value. The benchmark k L a H 2 value used, from Ljunggren et al. [15], was later on calibrated against experimental data resulting in a higher value (Table 8). Still, the mass transfer seems to be less efficient in the model not being able to fully describe the experimental data.