Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2006, p. 2134-2140, Vol. 72, No. 3
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.3.2134-2140.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Sandrine Alfenore,
Xavier Cameleyre,
Carole Molina-Jouve,
Jean-Louis Uribelarrea, and
Stéphane E. Guillouet*
Biotechnology and Bioprocess Laboratory, UMR-CNRS 5504, UMR-INRA 792, Département de Génie Biochimique et Alimentaire, Institut National des Sciences Appliquées, 135 Avenue de Rangueil, 31077 Toulouse Cedex, France
Received 21 July 2005/ Accepted 29 November 2005
|
|
|---|
|
|
|---|
Complementing biochemical and metabolic engineering experimental approaches, mathematical models may also be useful for interpretation and prediction of biochemical processes. Process models consist in coupling mass balance equations at the reactor scale and at the cellular scale (29). Cells can be modeled at different levels: unstructured models on the population level (15, 27) or structured models on the cellular level. The latter can be a compartmental model (41) or a metabolic regulator model (28). Constraints on carbon, energy, and redox potential linked to the metabolic network can be taken into account by a metabolic model that contains a complete set of metabolic pathways involved in biomass synthesis, energy production, substrate degradation, and cometabolite production (30, 35).
Metabolic flux analysis is generally used to quantify intracellular fluxes in the central metabolism of microorganisms (7, 9, 11, 22, 24, 43), but it has had a limited impact due to the lack of regulatory information in model formulation (10). Kinetic information can be introduced by the description of enzymatic rates with mathematical equations whose parameters are estimated with experimental data (26, 38). Another way of accounting for the regulatory aspects of cell metabolism consists in applying cybernetic principles (37) by using neural networks (34) or fuzzy logic-based models (13) when mechanistic details are missing.
Metabolic models were used in a few studies to predict (36) and to optimize (32) growth and by-product secretion. Metabolic flux analysis, combined with a black box kinetic model, then permits the introduction of dynamic information keeping the constraints imposed by metabolism.
With the aim of increasing the production of ethanol for use as a biofuel, glycerol production has to be minimized in order to better monitor the fate of carbon during the fermentation process and to limit technological problems in distillation units and/or separation processes after the fermentation stage. Reducing by-product yields is indeed a general focus in developing a competitive bioprocess, since the price of raw materials is a very important parameter in determining the overall economy of the process.
In this paper, the adopted strategy of reducing the production of glycerol during ethanolic fermentation was based on the combination of "process" (simulation, prediction, and process control) and metabolic approaches. On the basis of knowledge of the biological role of glycerol in redox balances, a strategy of fermentation monitoring was defined to reduce the need of surplus NADH oxidation responsible for glycerol synthesis during ethanol fermentation. A metabolic model was used to predict the operating conditions that would reduce glycerol production during ethanol fermentation in fed-batch process. Experimental validation of the simulation results was done by monitoring the inlet substrate feeding in order to maintain the respiratory quotient (RQ) (defined as the CO2 production to O2 consumption ratio) value between 4 and 5. Fermentation data were then discussed and compared in terms of growth, cell viability, and the production of ethanol and glycerol with previous highly productive ethanolic fermentations.
|
|
|---|
Fermentations.
Fed-batch experiments were performed in a 20-liter fermentor using the Braun Biostat E fermenting system without oxygen limitation. The temperature was maintained at 30°C, and the pH was maintained at 4 with the addition of a 14% (vol/vol) NH3 solution. All fermentations were performed with an exponential feeding of vitamins on the basis of the growth profile (2). The fermentor was connected to a computer. A homemade software program enabled the online acquisition of the controlled parameters (stirring rate, pH, temperature, partial pressure of the dissolved oxygen, NH3 and antifoam additions) and the monitoring and regulation of these parameters. The pressure in the bioreactor was regulated at 0.2 bar (relative pressure). The fermentor was flushed continuously with air through the sparger placed at the bottom of the reactor. The airflow rate was 313 normal liters h1. The stirring rate was fixed at 400 rpm until the concentration of dissolved oxygen reached 20% of saturation, and then it was increased in order to avoid any oxygen limitation in the culture.
The bioreactor was also equipped with a detection system of foam and controlled addition of Struktol (Schill & Seilacher, Germany) antifoam. The maximum amount of added Struktol was 1.2 ml. The fermentor was fed with a sterile concentrated glucose solution (about 700 g liter1) using a peristaltic pump (Masterflex). The glucose concentration in the reactor was automatically calculated by the software program on the basis of all the inlet and outlet volumes (base addition, sampling, and antifoam addition), the added mass (estimated instantly by weighing), and density, and the glucose concentration in the feeding solution was measured by high-performance liquid chromatography. Glucose density was measured by weighing a known volume.
Both the instantaneous CO2 production and the O2 consumption rates, and thus, the respiratory quotient, were calculated by gas analysis. As a consequence of the model prediction, the RQ value had to be regulated between 4 and 5 by controlling glucose flux. Thus, a PID (proportional-integral-derivative) controller was built to command the glucose pump in order to maintain the RQ value.
Gas analysis.
Outlet gas analysis was performed by mass spectroscopy after the condenser. The mass spectrometer (PRIMA 600s; VG Gas, Manchester, United Kingdom) was used for its accuracy to measure CO2, O2, N2, and Ar. The gas flow was ionized by electronic bombardment, and the gases were then separated according to their mass/charge ratio (m/e) by the application of a magnetic field and then collected by a detector. Gas analysis was conducted on the outlet flow of the reactor every 20 seconds, and gas analysis was performed on the inlet air every hour by commutation of the valves. The O2 consumption rate and the CO2 production rate were calculated from the differences between the inlet and outlet gas compositions, taking into account the evolution of the liquid volume in the reactor, inlet airflow (measured by a mass flowmeter), temperature, and pressure.
Analytical methods.
Biomass and concentrations of substrate and product were measured by the previously described analytical methods (2).
Determination of cell viability.
The viability of cells was determined by the methylene blue technique (23). A 200-µl sterile solution of methylene blue (0.3 mM in 68 mM sodium citrate) was mixed with 200 µl of a yeast suspension diluted to reach an optical density at 620 nm of 0.4 to 0.7. The mixture was shaken and placed in a Thoma counting chamber after 5 min of incubation. The numbers of stained (nonactive) or unstained (active) cells and the number of budding cells were counted in five different fields with a total of at least 200 to 300 cells. The percentage of viable cells was the number of unstained cells (lively cells) divided by the total number of cells (stained and unstained). Under our conditions, the mean viability (v) is estimated with an accuracy of 10% (17), which means that the interval v ± 0.1v contains the true value of viability with a probability of 95%.
Chemicals.
Chemical products (glucose, salts, oligo-elements, orthophosphoric acid, and NH3) were from Prolabo, the vitamins were from Sigma, and the sodium glutamate was from Merck. All products were of the highest analytical grade available.
|
|
|---|
Simulation by the metabolic descriptor model.
The metabolic descriptor model was composed in two equation systems: a metabolic component describing the cell metabolism through stoichiometric equations and classical mass balance equations describing the interactions between the cells and the bioreactor environment.
A metabolic matrix was constructed on the basis of the mass balance and on the pseudo-steady-state hypothesis on the intracellular intermediate metabolites. The formulation resulted in a set of linear equations expressed as a stoichiometric matrix S (m x n) with vectors for reaction rates or fluxes r (m x 1) and vectors for consumption or secretion rates of metabolites R (1 x n). Thus, S x r = R. This linear system was undetermined, and calculation of a unique solution was possible by considering a set of rates whose number was equal to the degrees of freedom of the system. For the major reactions of catabolism, anabolism, and energetics that we considered, see the supplemental material. The metabolic matrix construction resulted in a system containing 110 metabolites (n = 110) and 99 intracellular rates (m = 99) and a degree of freedom equal to 6.
Starting from stoichiometric equations of the metabolic pathway of Saccharomyces cerevisiae, the metabolic descriptor model allowed us to obtain expression of intracellular fluxes and extracellular specific rates as a function of six specific rates, such as the specific growth rate (µ), specific rate for glucose consumption (qglucose), specific rate for ethanol production (qethanol), specific rate for acetate production (qacetate), specific rate for glycerol production (qglycerol), and specific rate for succinate production (qsuccinate). In order to verify the possibility of producing ethanol during the fermentation process without any other by-products (particularly glycerol), the sign of metabolic fluxes and extracellular specific rates was first studied when the specific rates of production of glycerol, acetate, and succinate were null. A set of inequality functions of µ, qglucose, and qethanol was obtained. The resolution of these inequalities allowed us to determine a variation domain for qglucose [qglucose-min; qglucose-max], each being a function of µ and qethanol in such a way that qglycerol, qacetate, and qsuccinate were null: qglucose-min(µ, qethanol
qglucose
qglucose-max (µ, qethanol).
The ordinary differential equations (ODE) of the bioreactor (as a function of the specific rates) were solved by the Runge-Kutta method using Matlab software. The specific rates were determined each time ODE were computed with the metabolic descriptor model by constrained optimization. Aiming to minimize by-product production and concomitantly to maximize the ethanol production, the following objective function had to be minimized:
![]() | (1) |
![]() | (2) |
![]() | (3) |
![]() | (4) |
(coupling constant for the ethanol production specific rate) and ß (uncoupling constant for the ethanol production specific rate) were determined from experimental data. (
= 0.0612 and ß = 0.0063) (2). The simulation results of the metabolic descriptor model showed that ethanol production without any by-product synthesis was stoichiometrically possible, provided that a specific glucose feeding profile was applied during the fed-batch culture. This simulated profile for glucose feeding led us to monitor the glucose-limited condition in the medium during the culture. As a result, calculation of the respiratory quotient showed that its value had to be lower than 6.2.
The estimated biomass and ethanol concentrations after 30 h of culture were 66 g · liter1 and 100 g · liter1, respectively.
Experimental validation. (i) Fermentation parameters.
In order to validate the results of the simulation carried out by the metabolic descriptor model, two fed-batch fermentations were performed in a 20-liter fermentor in a mineral medium. Fermentation conditions were identical to those used for very high-performance ethanolic fermentation (2) except for the glucose feeding. As a result of the simulation, the RQ value was regulated within the range 4 to 5 by monitoring the glucose feeding rate using a PID controller.
The mass of ethanol and glycerol and biomass produced during the fed-batch fermentation process are shown in Fig. 1. The RQ value was kept constant at 4.5. As a result, the residual glucose concentration remained limiting (as required by the predictive model) during the first 17 hours and then increased slightly up to 5 g · liter1. Concomitant growth and production of glycerol and ethanol was observed during the 30 h of fermentation. Final concentrations of 18 g liter1 biomass, 85 g liter1 ethanol, and 1.7 g liter1 glycerol were reached. The fermentation led to global yields on glucose of 0.071 g g1, 0.34 g g1, and 0.006 g g1, for biomass, ethanol, and glycerol, respectively (Table 1).
![]() View larger version (23K): [in a new window] |
FIG. 1. Evolution of the biomass and ethanol, glycerol, and residual glucose concentrations (A) and respiratory quotient (B) in the RQ-controlled aerobic batch-fed fermentation.
|
|
View this table: [in a new window] |
TABLE 1. Comparison of ethanolic fermentation under full aeration, microaeration, and RQ-controlled conditionsa
|
(ii) Fermentation kinetic parameters.
The specific growth rate, the specific ethanol production rate, and the specific glycerol production rate were drawn versus the ethanol concentration along the fermentation (Fig. 2). µ and qethanol were coupled and decreased gradually, likely due to inhibition by the ethanol accumulated during the process. Maximum values were obtained at the beginning of the fermentation and reached 0.42 h1 for growth and 0.62 g g1 h1 for ethanol production. The specific glycerol production rate showed a different pattern, as its specific production rate remained constant within the range of 20 to 70 g liter1 ethanol.
![]() View larger version (17K): [in a new window] |
FIG. 2. Specific growth rate, specific ethanol production rate ( p), and glycerol production rates versus ethanol (EtOH) concentration for the RQ-controlled batch-fed fermentation.
|
![]() View larger version (15K): [in a new window] |
FIG. 3. Comparison of the percentage of methylene blue-stained cells (% MB staining, a measure of viability) as a function of ethanol concentration between RQ-controlled fermentation (open symbols) and high-performance ethanol (EtOH) fermentation (black symbols).
|
|
|
|---|
The validation was experimentally tested in glucose-limited fed-batch cultivation monitored by glucose feeding in order to control the RQ value at 4.5. The fermentation led to a final glycerol concentration of 1.7 g liter1, biomass concentration of 18 g liter1, and ethanol concentration of 85 g liter1. These data were compared with a fully aerated fed-batch fermentation with a sequenced glucose feeding (3) (Table 1). This reference high-performance ethanolic fermentation allowed us to reach an ethanol concentration of 132 g liter1 in 30 h under nonlimiting aeration conditions.
Here, the experimental results showed that the inhibition of the specific growth rate by ethanol was obtained for a lower concentration than in the reference experiment used for the prediction (85 g liter1 versus 100 g liter1). It was also observed that the oxygen consumption specific rate and oxygen to biomass yield (Table 1) were greater in this case than in the reference fermentation, leading to a lower expected value of the yield of ATP (YATP).
A study of the sensitivity of all the symbolic variables and parameters of the prediction model (Pc, YATP, phosphorus versus oxygen atom ratio, biomass macromolecular composition, ...) showed that only ATP yield and Pc induced a significant change in the biomass production (Fig. 4). By integrating into the model the ATP yield calculated from experimental data and a Pc value of 85 g · liter1, the predicted final biomass concentration can be decreased from 66 g liter1 to 20 g liter1 concomitantly with a final ethanol concentration of 85 g liter1 in 30 h of fermentation.
![]() View larger version (39K): [in a new window] |
FIG. 4. Sensitivity analysis of variables and parameters of the prediction model (µmax, maximum specific growth rate; Pc, critical ethanol concentration; , coupling constant for the ethanol production specific rate; ß, uncoupling constant for the ethanol production specific rate; Y(X;ATP), global ATP/biomass yield; P/O, phosphorus versus oxygen atom ratio).
|
Maintaining cell viability is one of the most important parameters in order to reach high ethanol concentrations in ethanolic fermentation. Comparing the results between our high-performance fed-batch process (2) and this RQ-controlled fed-batch process, a drop in cell viability was observed when glycerol production was minimized. This observation is unlikely to be explained by some potential energetic or redox imbalance resulting from the reduction in glycerol synthesis. The lower glycerol formation could have led to a perturbation in the redox balance, since anabolic and catabolic NADH production increased due to the higher biomass yield and lower glycerol yield. However, this deficit of NADH oxidation was likely compensated for by the observed increase in the oxygen consumption (Table 1). Such an increase in oxygen consumption is in agreement with reported data where overexpressing the transcription factor HAP4 gene in Saccharomyces cerevisiae led to a 3.7-fold decrease in qglycerol concomitantly with a twofold increase in qO2, a 25% reduction of qethanol, and a 2.2-fold increase in qacetate during aerobic batch culture compared to the values for the wild-type strain (6).
From an energetic point of view, the reduction of glycerol production in the RQ-controlled fermentation could also have resulted in an imbalance for ATP, since ATP production decreased due to a higher biomass yield and a lower ethanol yield. As previously, the increasing oxygen consumption may have balanced the cell energetics. Our data could then suggest a protective role of glycerol against ethanolic stress during ethanolic fermentation. It is well-known that glycerol formation in yeast is a response to the osmotic pressure of the medium (4, 16). Increasing glycerol production was also reported after submitting yeast cells to heat shock (12, 20). A growth-coupled glycerol formation was also reported to be enhanced during aerobic ethanolic fermentation with increasing temperatures within the range of 27 to 33°C, whereas a decoupling phenomenon occurred above 33°C, pointing out a possible role of glycerol in yeast thermal protection (1).
The industrial ethanol production industry is increasingly interested in the reduction of fermentation by-products, such as glycerol, since the price of the sugar source is a very important criterion for evaluation of the economic process. The hope is to increase the ethanol yield by minimizing glycerol production during ethanol fermentation. The data here showed that reducing the glycerol yield by 50% by glucose feeding monitoring led to a 20% decrease in the ethanol yield to the benefit of biomass (30%), acetate, and succinate compared to the reference fermentation (Table 1). When the formation of glycerol was reduced in a microaerobic ethanolic fermentation in continuous culture by a carefully controlled oxygenation, a decrease in the ethanol yield was also observed (8). RQ-controlled fermentation led to a lower maximum specific ethanol production rate and average ethanol productivity; however, it was possible to reach high ethanol production (85 g · liter1 in 30 h). The obtained mass ratio of ethanol production to glycerol production of 50 was the highest among the highly productive ethanol production processes reported in the literature (within 10 to 20) (33, 39). Previous attempts at reducing the glycerol formation by construction of mutant strains were reported in the literature (18, 19). The authors succeeded in increasing the ethanol yield by 10% with a 38% lower glycerol yield in anaerobic cultivations by changing the cofactor requirement in amino acid synthesis (18). Mutant strains impaired in glycerol synthesis were also obtained by deleting the glycerol dehydrogenase isoenzymes GPD1 (glycerol 3-phosphate dehydrogenase isoform 1) and GPD2 (19). Under aerobic conditions, this deletion resulted in a 12% increase in the ethanol yield with no glycerol production and higher acetate and succinate yields, but it also resulted in a dramatic reduction in the maximum specific growth rate and biomass formation. However, these two types of mutant strains led to poor ethanol production, since 10 g liter1 ethanol were produced in 16 h with no glycerol under aerobic conditions (19) and in 30 h with 1.3 g liter1 glycerol in anaerobic cultivations (18).
We showed here that it was possible to reduce glycerol production in very high-performance ethanolic fermentation process by monitoring the glucose feeding and that metabolic models may be valuable tools for predicting optimal experimental conditions. However, the microbiological engineering and metabolic engineering approaches cited above showed that managing glycerol production during ethanol fermentation remains challenging. Strategies will have to deal with metabolic and stress management issues in highly productive systems, since viable biomass formation is critical in achieving high-performance ethanolic fermentation.
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
These two authors worked equally on this project. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»