Dynamic metabolic modeling of a microaerobic yeast co-culture: predicting and optimizing ethanol production from glucose/xylose mixtures
Biotechnology for Biofuels volume 6, Article number: 44 (2013)
A key step in any process that converts lignocellulose to biofuels is the efficient fermentation of both hexose and pentose sugars. The co-culture of respiratory-deficient Saccharomyces cerevisiae and wild-type Scheffersomyces stipitis has been identified as a promising system for microaerobic ethanol production because S. cerevisiae only consumes glucose while S. stipitis efficiently converts xylose to ethanol.
To better predict how these two yeasts behave in batch co-culture and to optimize system performance, a dynamic flux balance model describing co-culture metabolism was developed from genome-scale metabolic reconstructions of the individual organisms. First a dynamic model was developed for each organism by estimating substrate uptake kinetic parameters from batch pure culture data and evaluating model extensibility to different microaerobic growth conditions. The co-culture model was constructed by combining the two individual models assuming a cellular objective of total growth rate maximization. To obtain accurate predictions of batch co-culture data collected at different microaerobic conditions, the S. cerevisiae maximum glucose uptake rate was reduced from its pure culture value to account for more efficient S. stipitis glucose uptake in co-culture. The dynamic co-culture model was used to predict the inoculum concentration and aeration level that maximized batch ethanol productivity. The model predictions were validated with batch co-culture experiments performed at the optimal conditions. Furthermore, the dynamic model was used to predict how engineered improvements to the S. stipitis xylose transport system could improve co-culture ethanol production.
These results demonstrate the utility of the dynamic co-culture metabolic model for guiding process and metabolic engineering efforts aimed at increasing microaerobic ethanol production from glucose/xylose mixtures.
An essential component of the quest for energy independence is to develop renewable sources of energy via the conversion of plant biomass to liquid transportation fuels. Lignocellulosic biomass is a heterogeneous collection of polymers that compose plant cell walls, namely celluloses, hemicelluloses, pectins, lignins and proteoglycans . The production of liquid fuels from biomass currently occurs in four major steps: pretreatment to make the feedstock more amenable to enzymatic degradation, hydrolysis of the cellulose and hemicellulose fractions to release sugar monomers and oligomers, fermentation of the released hexose and pentose sugars to produce fuels, and recovery of the fuels from the reactor bulk using separation technologies such as distillation [2, 3].
The fermentation of sugar mixtures that result from biomass hydrolysis is a significant bottleneck in the overall process. Few fermentative microbes that efficiently convert both hexose and pentose sugars to liquid fuels such as ethanol have been identified [4, 5]. The majority of current research efforts are focused on engineering multiple metabolic functionalities, such as the introduction of exogenous pentose consumption pathways, into a single organism . However, this approach often results in conversion inefficiencies due to bottlenecks in metabolic pathways and may place a heavy metabolic burden on the organism [7, 8]. In recent years, research has increasingly focused on the use of defined microbial consortia for biotechnology applications . Microbial communities perform the task of biomass degradation in nature, albeit at a rate much slower than required for an industrial process . Mixed cultures allow for the selection of microbes that are best suited for performing one task of the overall conversion process  and moves the engineering focus from introducing new functionalities to improving existing metabolic pathways. Other benefits of mixed culture systems include tunability and increased resistance to environmental stress [12, 13].
Saccharomyces cerevisiae is a robust, budding yeast that has been widely used for fermentation of refined corn starch to fuel ethanol . Because this yeast is Crabtree-positive, excess sugar that would otherwise overload its limited respiratory capacity overflows into the fermentative pathway. Thus, the yeast produces significant titers of ethanol when grown aerobically in batch culture . However, S. cerevisiae is unable to utilize pentose sugars, such as xylose and arabinose, that result from the hydrolysis of hemicellulose. The engineering of pentose metabolism into the S. cerevisiae genome has been achieved, but problems with co-factor imbalances and gene expression have hindered the efficiency of these mutants .
Another species of yeast, Scheffersomyces stipitis (formerly known as Pichia stipitis), can natively ferment xylose to ethanol but it retains a preference for glucose as the carbon source. The growth and metabolite profile of this yeast is highly sensitive to the oxygenation level. Unlike S. cerevisiae, S. stipitis is a Crabtree-negative yeast that only produces ethanol under oxygen-limited conditions . S. stipitis is among the most efficient native fermenters of xylose when grown under microaerobic culture conditions. In certain aeration regimes, however, this yeast can reassimilate ethanol while often simultaneously consuming other growth substrates . Unlike most fermentative microorganisms, wild-type S. stipitis is unable to grow anaerobically. Insertion of the URA1 gene from S. cerevisiae has been shown to enable S. stipitis to grow anaerobically on glucose . However, there are no known gene insertions that allow for S. stipitis anaerobic growth on xylose [20, 21]. Thus, aeration level is a critical operating variable that must be tightly regulated to maximize xylose conversion to ethanol by wild-type S. stipitis.
S. cerevisiae and S. stipitis have been co-cultured for the production of ethanol from glucose and xylose mixtures [22–24]. In these studies, a respiratory-deficient strain of S. cerevisiae was used so the dissolved oxygen concentration could be more easily controlled at a level that was favorable for ethanol production by S. stipitis. Because it lacks the respiratory capability of a wild-type strain, respiratory-deficient S. cerevisiae could not utilize non-fermentable carbon sources such as ethanol once glucose had been exhausted . The inability of these strains to grow on the ethanol produced under certain microaerobic conditions resulted in high ethanol titers. In fact, co-culturing these two yeasts on a mixture of glucose and xylose has been shown to yield more ethanol than can be produced by either yeast alone .
Constraint-based analysis using genome-scale metabolic reconstructions is a widely used computational tool for predicting how fluxes through microbial metabolic pathways will respond to changes in the culture environment or gene deletions/insertions [26, 27]. With the addition of substrate uptake kinetics and extracellular mass balances on growth-limiting substrates and metabolic byproducts, these steady-state models can be adapted to predict culture dynamics that are critical in batch and fed-batch fermentations [28–30]. Recently, constraints-based modeling has been applied to mixed-culture systems [31–33]. In addition to describing individual species metabolism, mixed-cultures models must account for possible interactions between the species as well as postulate a community objective that captures the combined metabolic behavior. While the assumption that each species attempts to maximize its own growth rate is most common, other community objectives that capture more complex behavior such as altruism can be employed .
In this study, we used dynamic flux balance analysis (DFBA) to drive the experimental optimization of ethanol production from a respiratory-deficient S. cerevisiae and wild-type S. stipitis co-culture growing microaerobically on a mixture of glucose and xylose. A dynamic model of S. stipitis metabolism was developed from a recently published genome-scale reconstruction  by estimating glucose and xylose uptake parameters from batch pure culture data. Because the respiratory-deficient S. cerevisiae mutant had non-specific genetic alterations, we considered several plausible modifications of the wild-type metabolic network  to develop a S. cerevisiae dynamic model consistent with batch pure culture data. Co-culture experiments demonstrated that S. cerevisiae competed less successfully for glucose than expected from combining the pure culture dynamic models under the community objective of total biomass maximization. A revised co-culture model with reduced S. cerevisiae glucose uptake was shown to provide accurate predictions of batch co-culture data over a range of microaerobic growth conditions. The experimentally validated model was used to optimize the batch ethanol productivity by adjusting the inoculum concentrations and the aeration level. Finally, the model was used to examine what modifications to the S. stipitis xylose transport system would yield the largest improvements in co-culture ethanol production.
Results and discussion
S. stipitis pure cultures
Our initial dynamic model of S. stipitis metabolism did not include a balance on dissolved oxygen in the culture media. We were unable to satisfactorily fit this model to measured biomass and ethanol concentration profiles through adjustment of the sugar uptake rate parameters (see below). Because the dissolved oxygen (DO) concentration was above the assumed microaerobic value during the initial portion of the batch, we added the DO balance (8) to account for the apparent gas–liquid mass transfer limitations. The mass transfer coefficient (k L a) was determined from the gas sparge rate using a linear correlation (see Materials and Methods). Starting with literatures values when available , the uptake rate parameters for glucose (9), xylose (10) and oxygen (11) were determined by minimizing the least-squares difference between the experimental and predicted concentration profiles. Parameter adjustments were made by trial-and-error using glucose, xylose, biomass, and ethanol concentration profile measurements collected from two microaerobic batch fermentations performed at air sparging rates of 25 (k L a = 5.5 h-1) and 50 cc/min (k L a = 10.1 h-1). The resulting uptake parameters are compiled in Table 1.
The S. stipitis dynamic models with and without the DO balance (8) are compared in Figure 1 at an air sparging rate of 50 cc/min using initial conditions of 0.40 g/L biomass, 15.7 g/L glucose and 8.4 g/L xylose. A k L a value of 10.1 h-1 was used in the DO balance, while the model without the balance used a constant DO value of 0.0072 mM, the measured level at the end of the batch fermentation. The constant DO model could not reproduce the lag in ethanol production that resulted from the high level of dissolved oxygen present during the early stages of the fermentation. Additionally, the constant DO model produced lower biomass yields and higher ethanol yields than observed in experiment. Once the DO balance was added, the model could more accurately predict the time at which S. stipitis began microaerobic substrate consumption and ethanol production started to outpace biomass growth. The other S. stipitis pure culture experiment used for uptake parameter fitting of the variable DO model was performed with a k L a of 5.5 h-1 and initial conditions of 0.25 g/L biomass, 16.8 g/L glucose and 8.4 g/L xylose (Figure 2A). As before, the model produced very accurate predictions of the glucose, xylose, biomass and ethanol concentration profiles with the uptake parameter values in Table 1.
To examine extensibility of the S. stipitis dynamic model to different microaerobic conditions, batch experiments were performed at two additional k L a values and model predictions were generated with the same uptake parameter values listed in Table 1. Figure 2B shows a comparison of measured and predicted concentration profiles for a batch culture performed with a k L a of 7.6 h-1 and initial conditions of 0.25 g/L biomass, 16.8 g/L glucose and 8.0 g/L xylose. The simulation results were generally satisfactory, although the model predicted slightly faster consumption of both substrates and small overprediction of biomass production throughout most of the batch. A second validation experiment was performed with a larger k L a of 12.6 h-1 and initial conditions of 0.2 g/L biomass, 15.6 g/L glucose and 8.6 g/L xylose. As shown in Figure 2C, the model adequately captured the fermentation dynamics throughout the first 15 hours but it was unable to predict the subsequent reassimilation of ethanol that is evident in the data.
Although the S. stipitis dynamic metabolic model generated satisfactory predictions over a range of microaerobic conditions, the model failed to capture the ethanol and biomass concentration profiles at higher aeration levels once glucose had been exhausted. Ethanol can be simultaneously produced and reassimilated at these higher k L a values , resulting in overprediction of ethanol secretion and underprediction of biomass formation by our model. This phenomenon is likely due to the relative ease at which ethanol diffuses across the plasma membrane , while the uptake of xylose has been shown to be facilitated by oxygen . Since xylose was more energetically favorable, the LP solver returned a flux distribution with simultaneous xylose and ethanol uptake only if there was excess oxygen available beyond that required for xylose metabolism. Even in absence of xylose, the S. stipitis metabolic reconstruction predicted that k L a values greater than 34 h-1 would be required just to meet the ATP maintenance demand for growth on ethanol. By contrast, we observed ethanol assimilation in our experiments for k L a values as low as 12.6 h-1. Without altering the objective function or the ATP maintenance coefficient of the iBB814 model, the observed disparity between model and experiment cannot be resolved. Fortunately, the highest ethanol yields and titers observed experimentally were obtained at aeration levels under which ethanol assimilation did not occur. As a result, we do not consider k L a values above ~10 h-1 in the remainder of the paper.
S. cerevisiae 311 pure cultures
Because S. cerevisiae 311 was created by applying ethidium bromide to a wild-type strain, the genetic alterations that produced the respiratory-deficient mutant were non-specific. As a result, modifications to the wild-type iMM904 metabolic network necessary to describe S. cerevisiae 311 metabolism were unknown. We used the metabolic network model to identify putative gene knockouts consistent with the known alteration of mitochondrial DNA by ethidium bromide treatment  and to obtain predictions in agreement with our data. All the gene knockouts considered prevented growth on non-fermentable substrates and reduced oxygen demand. To determine the reactions to be deleted, model fluxes in the mitochondrial compartment were systematically set to zero until FBA simulations matched the S. cerevisiae 311 phenotype of lower biomass yields, higher ethanol yields and smaller oxygen demands than the wild-type strain. This process identified knockouts of genes encoding two enzymes: ubiquinol-6 cytochrome c reductase and mitochondrial cytochrome c oxidase. Modifications to these enzymes may well have occurred in S. cerevisiae 311, as other respiratory-deficient S. cerevisiae mutants have been created through the direct removal of the encoding genes .
The S. cerevisiae 311 dynamic model was built on the modified iMM904 metabolic network with glucose uptake parameters estimated from anaerobic culture data and oxygen uptake parameters estimated from aerobic culture data with the glucose parameters fixed . The resulting parameter values are listed in Table 1. Figure 3A shows a comparison of measured concentration profiles and model predictions for the anaerobic batch fermentation initiated with 0.28 g/L biomass, 18.0 g/L glucose and 1.8 g/L residual ethanol from the preculture shake flask. Figure 3B shows results for the aerobic fermentation initiated with 0.20 g/L biomass, 15.5 g/L glucose and 1.8 g/L ethanol. In each case, the dynamic model accurately predicted the entire ethanol profile and the glucose and biomass profiles in the second half of the batch. However, the model overpredicted glucose uptake and ethanol synthesis during the first half of the batch. While this discrepancy could have been attributable to the preculture cells used for inoculation being in stationary rather than exponential phase , we were unable to verify this hypothesis through additional experiments in which the preculture cells were harvested earlier.
To examine extensibility of the S. cerevisiae dynamic model to different aeration levels, a microaerobic batch experiment was performed at a k L a of 5.5 h-1 and model predictions were generated with the same uptake parameter values listed in Table 1. Figure 3C shows the results obtained for the fermentation initiated with 0.26 g/L biomass, 18.7 g/L glucose and 1.4 g/L residual ethanol. The model produced excellent agreement with data, and only small deviations in the biomass concentration during the initial portion of the batch and in the ethanol concentration following glucose exhaustion were observed.
A preliminary dynamic co-culture model was developed by direct combination of the two individual yeast models assuming no species interactions other than competition for glucose (see Materials and Methods). This approach resulted in faster glucose consumption and a higher final concentration of S. cerevisiae 311 than were observed in our microaerobic batch fermentations (not shown). We found that these discrepancies could be partially rectified by reducing the S. cerevisiae maximum glucose uptake rate from the pure culture value of 21.5 mmol/gdw/h to 18.5 mmol/gdw/h. Otherwise, the substrate uptake parameters remained fixed at the pure culture values listed in Table 1.
Figure 4A shows a comparison of the resulting model predictions and experimental data collected for three batch fermentations at a k L a of 5.5 h-1 with an equal inoculum of 0.15 g/L of each yeast species grown on 19.1 g/L glucose and 9.7 g/L xylose. Figure 4B shows corresponding results for a k L a of 10.1 h-1, 0.10 g/L of each yeast species, 16.5 g/L glucose and 7.9 g/L xylose. With the implemented change in the S. cerevisiae maximum glucose uptake rate, the dynamic co-culture model proved to be as accurate as the individual species models. The glucose consumption rate was slightly overpredicted during the initial portion of the batch at both aeration levels. Small discrepancies in the predicted biomass and ethanol concentrations were observed during the glucose consumption phase for k L a = 5.5 h-1, while small but longer lasting deviations were apparent in the biomass and xylose concentrations at k L a = 10.1 h-1.
To examine prediction accuracy of dynamic co-culture model at aeration levels not used for parameter adjustment, we performed three additional co-culture fermentations at a k L a of 7.6 h-1 with 0.145 g/L of each yeast as the inoculum and initial sugar concentrations of 16.5 g/L glucose and 8.5 g/L xylose. The model and data comparisons shown in Figure 4C are qualitatively similar to those obtained at the other two aeration levels. The results in Figure 4 collectively demonstrate that the co-culture model provided satisfactory predictions over a range of microaerobic conditions corresponding to k L a < ~10 h-1.
Following experimental validation, we utilized the dynamic co-culture model to computationally determine optimal bioreactor operating conditions for maximization of ethanol productivity. The productivity was defined as the final ethanol concentration divided by the batch time, which was chosen as the time at which the xylose concentration dropped below 0.5 g/L. Our simulations indicated that the aeration level had little effect on the ethanol titer but strongly affected the batch time due to the xylose consumption rate . Our simulations also demonstrated that the relative amount of each yeast in the inoculum strongly affected the productivity through both the ethanol titer and the batch time. Therefore, ethanol productivity was optimized by adjusting the aeration level and inoculum concentrations. Rather than perform rigorous optimization , dynamic simulations were run with different combinations of the k L a value and the initial yeast concentrations, and the case that yielded the largest productivity was deemed the optimal solution. All simulations were performed with a mixture of 16.0 g/L glucose and 8.0 g/L xylose, while the total inoculum concentration was fixed at 1.0 g/L to allow direct comparison of results for different inocula.
The k L a value was constrained to be below 10.1 h-1 since larger values had the potential to result in undesirable ethanol reassimilation by S. stipitis that was not captured by the co-culture model. Regardless of the inoculum, we found that increasing k L a improved productivity due to a decrease in batch time that resulted from enhanced xylose consumption. Consequently, the optimal solution was achieved at k L a = 10.1 h-1, and the optimization problem was reduced to determining the inoculum concentrations. This result demonstrated the importance of developing improved metabolic reconstructions that more accurately predict the relationship between aeration level and the onset of ethanol reassimilation. Figure 5 shows the effect of the initial S. cerevisiae 311 concentration on the ethanol titer, batch time and ethanol productivity. The initial S. stipitis concentration is not shown since the total inoculum concentration was constrained to be 1.0 g/L. The ethanol titer increased with increasing S. cerevisiae concentration since this yeast converts glucose to ethanol at higher yields than S. stipitis. Conversely, the batch time decreased with increasing S. cerevisiae concentration because S. stipitis could better compete for glucose during the initial growth phase and the larger S. stipitis biomass concentration that resulted allowed for more rapid xylose consumption. These two competing effects produced a clear optimum in the ethanol productivity at 0.1 g/L S. cerevisiae and 0.9 g/L S. stipitis.
To validate the model predictions, we performed three batch fermentations at the optimal conditions identified in silico. Figure 6 shows the averaged result of these experiments alongside the predicted growth and metabolite concentration curves. Although we were unable to obtain measurements during the first six hours, the optimized co-culture model produced excellent agreement with the measured concentration profiles after that time. The model predicted a final ethanol titer of 9.43 g/L and a productivity of 0.55 g/L/h, while the three fermentations were averaged to produce an ethanol titer of 9.07 g/L and a productivity of 0.521 g/L/h. According to our in silico analysis at k L a = 10.1 h-1, the optimized co-culture would outperform both a pure S. stipitis culture (9.43 g/L ethanol at a productivity of 0.505 g/L/h) and a co-culture with equal inoculum (9.07 g/L ethanol at a productivity of 0.469 g/L/h).
In silico transporter engineering
Our simulations demonstrated that S. stipitis xylose metabolism was the rate limiting process that limited co-culture conversion efficiency. Moreover, xylose transport has been identified as the main bottleneck in pentose sugar metabolism with pure S. stipitis cultures . Therefore, we used the dynamic co-culture model to predict the effects of modifying S. stipitis xylose transport parameters on ethanol productivity. To model engineering of the associated transport proteins, the nominal values listed in Table 1 for the maximum xylose uptake rate (v z,max ), the xylose uptake saturation constant (K z ), and the glucose inhibition constant for xylose uptake (K igz ) were perturbed both upward and downward to determine their impact. The ethanol inhibition constant for xylose uptake (K iez ) was excluded from detailed analysis because the ethanol concentrations achieved in silico were too small to cause significant xylose uptake inhibition. For the remaining three parameters, scaled sensitivity coefficients were calculated as:
where is the nominal parameter value, Δ p is the parameter change, is the optimal ethanol productivity obtained with the nominal parameter value, and Δ p is the predicted change in the ethanol productivity that results from the parameter change. All simulations were initialized with 16 g/L glucose and 8 g/L xylose at the optimum k L a = 10 h-1 and inoculum of 0.1 g/L S. cerevisiae 311 and 0.9 g/L S. stipitis.
Figure 7 shows the ethanol productivities and scaled sensitivity coefficients obtained when v z,max was changed upward and K z and K igz were changed downward. Increasing the maximum xylose uptake rate offered the largest improvement in ethanol productivity as indicated by the relatively large sensitivity coefficients, with a productivity enhancement of 14% achieved with a 25% v z,max increase (Figure 7A). As indicated by the small sensitivity coefficients, more modest gains in ethanol productivity were predicted for decreases in the xylose uptake saturation constant (Figure 7B) and the glucose inhibition constant for xylose uptake (Figure 7C). However, complete elimination of S. stipitis diauxic growth by removal of the glucose inhibition term from the xylose uptake expression resulted in a large productivity enhancement of 85% (not shown). The simulated parameter changes could potentially be realized through protein engineering of the sugar transporters [44, 45]. Because many S. stipitis transporters uptake both glucose and xylose , a complementary approach could be overexpression of specific transporters with more favorable xylose uptake characteristics and downregulation or deletion of less favorable transporters .
Due to the ability of Scheffersomyces stipitis to efficiently convert pentose sugars to ethanol, co-cultures of S. stipitis and respiratory-deficient Saccharomyces cerevisiae are promising for renewable ethanol production from glucose/xylose mixtures [23, 47]. Starting from genome-scale metabolic reconstructions of the individual yeasts, we developed a dynamic metabolic model of S. cerevisiae/S. stipitis co-cultures that was capable of predicting microaerobic metabolism in batch culture. As a first step towards co-culture modeling, dynamic models were developed for the individual yeasts using data from pure culture experiments. The Crabtree-negative yeast S. stipitis was shown to be very sensitive to the aeration level such that the incorporation of a dissolved oxygen balance was necessary to capture non-microaerobic oxygen levels during the initial batch phase and the subsequent switch from respiratory and to fermentative growth. The S. stipitis dynamic model was able to accurately reproduce measured glucose, xylose, biomass and ethanol concentration profiles over a range of microaerobic growth conditions characterized by oxygen-liquid mass transfer coefficient (k L a) values below ~10 h-1. However, the S. stipitis model was unable to capture ethanol reassimilation observed at higher oxygenation levels due to model energetics that required a k L a greater than 34 h-1 just to meet the ATP maintenance demand for growth on ethanol. Therefore, subsequent co-culture experiments and simulations were restricted to microaerobic conditions with k L a values less than ~10 h-1.
The procedure used to create the respiratory-deficient strain S. cerevisiae 311 introduced unknown genetic alterations that could not be directly implemented in the wild-type metabolic network. We used flux balance analysis to screen putative gene knockouts in the mitochondrial compartment that would reproduce the S. cerevisiae 311 phenotype of lower biomass yields, higher ethanol yields and smaller oxygen demands than the wild-type strain. The best agreement with S. cerevisiae 311 pure culture data was obtained for deletions of two genes that encode the ubiquinol-6 cytochrome c reductase and mitochondrial cytochrome c oxidase enzymes involved in the electron transport chain. With the fluxes through the two reactions catalyzed by these enzymes constrained to zero, the oxygen demand was greatly reduced because oxygen was not needed to serve as a final electron acceptor for ATP synthesis. These results were consistent with more directed genetic engineering efforts to create respiratory-deficient S. cerevisiae mutants through direct removal of these two genes .
A preliminary version of the dynamic co-culture model was developed by directly combining the dynamic models of the two yeast species under the assumption that each species attempted to maximize its individual growth rate and the only interspecies interaction was competition for glucose substrate. By comparison of model predictions to batch co-culture data collected at several microaerobic conditions, we found that the preliminary model overpredicted the glucose consumption and S. cerevisiae biomass formation rates. Much improved predictions were obtained by reducing the S. cerevisiae maximum glucose uptake rate from its pure culture value of 21.5 mmol/gdw/h to 18.5 mmol/gdw/h. This parameter change suggested the presence of unmodeled species interactions that resulted in S. stipitis having an antagonistic effect on S. cerevisiae growth. Crabtree-negative yeasts like S. stipitis have an established advantage when competing against a Crabtree-positive yeast such as S. cerevisiae for the same growth-limiting substrate . This effect is further magnified when the Crabtree-positive species is a respiratory deficient mutant . One possible cause for this effect was the competition for nutrients other than glucose. Growth under nitrogen limited conditions has been shown to slow the uptake of glucose and other hexoses in wild-type S. cerevisiae. Respiratory-deficient S. cerevisiae could be at a disadvantage in competing for nitrogen sources such as ammonium when grown with respiratory competent S. stipitis. Because sugar uptake parameters reflect the action of numerous transport systems, each having distinct mechanisms and affinities, the decrease in this parameter may indicate differences in transporter expression between pure and co-cultures of the two microbes.
Following experimental validation, the dynamic co-culture model was used to compute the aeration level and initial cell concentrations that maximized batch ethanol productivity. Our in silico analysis suggested that the co-culture should be operated at the largest possible k L a value before the onset of ethanol reassimilation by S. stipitis because this aeration level maximized the rate of xylose conversion to ethanol. Because the S. stipitis metabolic network model did not accurately predict this transition point, the optimization was constrained by experimental data such that the optimal solution was achieved at k L a = 10.1 h-1. Optimization of the initial cell concentrations produced an inoculum with 90% S. stipitis and only 10% S. cerevisiae to form sufficient S. stipitis biomass following glucose exhaustion for efficient conversion of the remaining xylose. Experimental validation of the optimal solution showed that the co-culture model provided excellent agreement with measured concentration profiles. Compared to an unoptimized co-culture with equal initial cell concentrations, the optimized co-culture was shown experimentally to produce a 11% improvement in ethanol productivity. Additional in silico analysis indicated that the co-culture would yield higher ethanol productivities than a S. stipitis pure culture due to the higher glucose to ethanol conversion efficiency of S. cerevisiae.
The relatively slow conversion of xylose by S. stipitis is the main bottleneck that limits overall co-culture performance. The dynamic co-culture model was used to explore the potential impact of transporter engineering efforts aimed at enhancing xylose uptake by increasing the maximum xylose uptake rate, reducing xylose uptake saturation and reducing glucose inhibition of xylose uptake. Based on sensitivity analysis for moderate parameter changes, the maximum xylose uptake rate was predicted to yield the largest improvement in ethanol productivity. However, substantially larger improvements were predicted for complete elimination of glucose catabolite repression in S. stipitis because this modification eliminated the diauxic growth pattern and allowed xylose consumption to commence at the beginning of the batch. Taken together, our computational results suggest that an engineered S. stipitis strain which rapidly consumes xylose in the presence of glucose and does not reassimilate ethanol under microaerobic conditions would be an ideal candidate for enhancing co-culture performance. Future work should focus on the creation of such S. stipitis strains as well as the model-based characterization of co-culture performance for growth on actual biomass hydrolysates with inhibitory compounds.
The wild-type S. stipitis strain NRRL Y-7124 (ATCC 58376) was used in this study. S. cerevisiae 311 (ATCC 42511), a mutant that was created by treating a wild-type strain with ethidium bromide , was chosen as the respiratory-deficient S. cerevisiae strain. Stocks of the two yeasts were stored at 4°C on YM agar slants.
All pure and mixed cultures were performed in a synthetic yeast minimal medium . The composition per liter of water was 1.00 g MgSO4 ·7 H2O, 1.10 g/L KCl, 0.15 g CaCl2 · 2 H2O, 1.00 g (NH4)2HPO4, 8.75 g/L (NH4)2SO4, 60.3 mg myo-inositol, 30.0 mg Ca-panthothenate, 6.0 mg thiamine-HCl, 1.5 mg pyridoxine-HCl, 0.03 mg biotin, 10.6 mg MnSO4 · H2O, 9.0 mg ZnSO4 · 7 H2O, 5.0 mg FeSO4 · 7 H2O, and 2.4 mg CuSO4 · 2 H2O. Pre-cultures in media containing 20 g/L glucose and 20 g/L xylose for S. cerevisiae and S. stipitis, respectively, were grown at 30°C for 36 hours on a shake table set at 175 RPM. The inoculum concentration for each experiment was determined by calculating the volume of preculture required to obtain the target initial concentration of each cell type using the measured biomass concentration in the shake flask media.
All fermentations were performed in a HEL BioX array of 4 250 mL vessels situated in a shared block that provided both electric heat and independent magnetic agitation (HEL Group Ltd., Barnet, UK). Electrochemical probes monitored the dissolved oxygen and pH in each vessel, while individual thermocouples recorded the media temperatures. Bioreactor cultivations were performed at a constant temperature of 30°C and pH of 5, the optimal growth conditions for each yeast species . The pH in each vessel was controlled by the automatic addition of 1 N sulfuric acid or 2 N NaOH. Glucose and xylose were autoclaved separately and added to the growth media in the amounts indicated for each experiment. Antifoam A was added to the reactors as necessary to prevent foaming.
Aeration of culture media was found to be a crucial operating variable. The agitation speed was held constant at 500 RPM for both pure and mixed culture fermentations. The gas flow rate into each reactor was altered according to the aeration level required for each experiment. A linear relationship between the gas sparge rate and the gas–liquid oxygen mass transfer coefficient (k L a) was determined using the static gassing out method . Purely aerobic cultures were aerated with pure oxygen, while microaerobic fermentations were aerated with house air passed through a HEPA-VENT filter (Whatman Ltd., Kent, UK).
Total cell weight was measured using a correlation between OD595 measured on a WPA UV1101 Biotech Photometer (Biochrom Ltd., Cambridge, UK) and dry cell weight. Cell counts of S. cerevisiae and S. stipitis in co-culture were performed on a hemacytometer in triplicate and averaged. A typical cell count considered approximately 50 S. cerevisiae cells and 200 S. stipitis cells. Conversion factors between dry cell weight and number of cells were found by drying pure culture samples of each yeast after cell counts had been performed. These factors were found to be 0.006943 gdw/L for S. cerevisiae and 0.001944 gdw/L for S. stipitis. Ethanol, glucose and xylose concentrations were measured by YSI 2700 SELECT biochemistry analyzers (YSI Inc., Yellow Springs, OH) configured with the enzyme-bound membranes specified for each metabolite. Raw readings were interpreted by 2700 Xylose PC Software (YSI Inc., Yellow Springs, OH) to resolve cross-talk between the xylose and glucose specific membranes.
The most comprehensive S. cerevisiae metabolic reconstruction currently available, iMM904 , was used for pure and mixed culture simulations. The fully compartmentalized network was reconstructed from 904 genes and accounts for 1228 metabolites and 1412 reactions. S. stipitis metabolism was simulated with iBB814 , the first published genome-scale reconstruction for this organism. This model accounts for 814 genes, 971 metabolites and 1371 reactions that are compartmentalized in the cytoplasm, mitochondria, and extracellular space. Following the publication of iBB814, a slightly more detailed S. stipitis metabolic reconstruction was developed . We do not anticipate that the use of this alternative reconstruction would significantly alter the results reported in this paper.
The S. cerevisiae/S. stipitis co-culture model was constructed by combining the iMM904 and iBB814 stoichiometric matrices into a single matrix . Flux distributions for S. cerevisiae (v c ) and S. stipitis (v s ) were calculated by solving the following linear program based on the assumption that the two species attempted to maximize their individual growth rates:
where the subscript i represents the species, A i is the matrix of stoichiometric coefficients, v i is the vector of reaction fluxes including exchange fluxes, vi,min and vi,max are vectors of lower and upper flux bounds, μ i is the growth rate, and w c and w s are vectors of experimentally determined weights that represent the contribution of each flux to biomass formation in S. cerevisiae and S. stipitis, respectively. Other than competing for the common substrate glucose, the two yeasts were assumed to grow independently without species interactions. Therefore the co-culture objective function μ was assumed to be the sum of the individual species growth rates, and the inclusion of multi-level objective functions  was deemed unnecessary. The co-culture model was also used to simulate pure cultures of S. cerevisiae and S. stipitis by constraining all fluxes of the unmodeled organism to zero.
The steady-state flux balance model (2) was extended to a dynamic model through the addition of the following extracellular mass balance equations:
where X c and X s are the biomass concentrations of S. cerevisiae and S. stipitis, respectively, G, Z, and E are the concentrations of glucose, xylose, and ethanol, respectively, ve,c and ve,s are ethanol exchange fluxes, v g,c is the glucose uptake rate for S. cerevisiae, and v g,s and v z,s are the glucose and xylose uptake rates, respectively, for S. stipitis. An equation for the dissolved oxygen concentration (O) was necessary to accurately describe microaerobic growth of S. stipitis (see results). In this equation (8), v o,c and v o,s are oxygen exchange fluxes, k L a is the volumetric mass transfer coefficient of oxygen from sparged gas to the culture medium, and O* is the saturation concentration of oxygen. For all simulations, O* was taken to be 0.24 mM, the saturation concentration for water at 30°C and 1 atm.
The following substrate uptake expressions were used to calculate upper bounds on the actual sugar and oxygen uptake rates:
where vg,max, vz,max and vo,max are the maximum uptake rates of each substrate, K g , K z and K o are corresponding saturation constants, K ieg and K iez are ethanol inhibition constants, and K igz is a glucose inhibition constant. The glucose (9) and xylose (10) uptake rates were assumed to follow Michaelis-Menten kinetics with an additional inhibitory term that reflects growth rate suppression at high ethanol concentrations . The glucose inhibition term added to the xylose uptake kinetics accounted for diauxic growth where S. stipitis favors glucose over xylose as the carbon source. The oxygen uptake rate was calculated from a Michaelis-Menten expression based on the dissolved oxygen content of the medium .
Pure and mixed culture dynamic flux balance models were solved using the Mosek optimization toolbox (Mosek ApS, Denmark) to resolve the linear program for intracellular metabolism within Matlab (Mathworks, Natick, MA) . Because S. cerevisiae could not meet the non-growth associated ATP maintenance demand during the xylose-only consumption phase, the maintenance flux was constrained to zero after glucose depletion to prevent the LP solver from returning zero fluxes for the S. stipitis network. Due to time-scale differences between the sugar and oxygen consumption rates, the differential equation system (3)--(8) exhibited a high degree of stiffness. To reduce the time required to generate large numbers of DFBA simulations for parameter fitting and in silico culture optimization, Matlab stiff ODE solvers ode15s and ode23tb were used to obtain approximate solutions. An ODE solver with greater accuracy, ode23, was used to generate model predictions once parameters had been estimated or an optimum had been determined. A typical co-culture batch simulation that was solved in two minutes with ode15s required five hours with ode23.
Growth rate (1/h)
Ethanol concentration (g/L)
Glucose concentration (g/L)
Glucose uptake saturation constant (g/L)
Ethanol inhibition constant (g/L)
Glucose inhibition constant (g/L)
Oxygen mass transfer coefficient (h-1)
Oxygen uptake saturation constant (mM)
Xylose uptake saturation constant (g/L)
Dissolved oxygen concentration (mM)
Dissolved oxygen saturation concentration (mM)
S. cerevisiae Flux distribution (mmol/gdw/h)
c: S. cerevisiae ethanol flux (mmol/gdw/h)
s: S. stipitis ethanol flux (mmol/gdw/h)
c: S. cerevisiae glucose uptake rate (mmol/gdw/h)
max: Maximum glucose uptake rate (mmol/gdw/h)
s: S. stipitis glucose uptake flux (mmol/gdw/h)
max: Vector of upper bound flux constraints (mmol/gdw/h)
min: Vector of lower bound flux constraints (mmol/gdw/h)
c: S. cerevisiae oxygen uptake rate (mmol/gdw/h)
max: Maximum oxygen uptake rate (mmol/gdw/h)
s: S. stipitis oxygen uptake rate (mmol/gdw/h)
S. stipitis Flux distribution (mmol/gdw/h)
max: Maximum xylose uptake rate (mmol/gdw/h)
s: S. stipitis xylose uptake flux (mmol/gdw/h)
S. cerevisiae biomass reaction
S. stipitis biomass reaction
S. cerevisiae concentration (g/L)
S. stipitis concentration (g/L)
Xylose concentration (g/L).
Lynd LR, Weimer PJ, van Zyl WH, Pretorius IS: Microbial Cellulose Utilization: Fundamentals and Biotechnology. Microbiol Mol Biol Rev 2002, 66: 506-577.
Wyman CE: What is (and is not) vital to advancing cellulosic ethanol. Trends Biotechnol 2007, 25: 153-157.
Lynd LR, Laser MS, Bransby D, Dale BE, Davison B, Hamilton R, Himmel M, Keller M, McMillan JD, Sheehan J, Wyman CE: How biotech can transform biofuels. Nat Biotech 2008, 26: 169-172.
Gowen CM, Fong SS: Exploring Biodiversity for Cellulosic Biofuel Production. Chem Biodivers 2010, 7: 1086-1097.
Gírio FM, Fonseca C, Carvalheiro F, Duarte LC, Marques S, Bogel-Łukasik R: Hemicelluloses for fuel ethanol: A review. Bioresour Technol 2010, 101: 4775-4800.
Ha S-J, Galazka JM, Rin Kim S, Choi J-H, Yang X, Seo J-H, Louise Glass N, Cate JHD, Jin Y-S: Engineered Saccharomyces cerevisiae capable of simultaneous cellobiose and xylose fermentation. Proc Natl Acad Sci 2011, 108: 504-509.
Pickens LB, Tang Y, Chooi Y-H: Metabolic Engineering for the Production of Natural Products. Ann Rev Chem Biomol Eng 2011, 2: 211-236.
Almeida JRM, Runquist D, Sànchez Nogué V, Lidén G, Gorwa-Grauslund MF: Stress-related challenges in pentose fermentation to ethanol by the yeast Saccharomyces cerevisiae. Biotechnol J 2011, 6: 286-299.
Brenner K, You L, Arnold FH: Engineering microbial consortia: a new frontier in synthetic biology. Trends Biotechnol 2008, 26: 483-489.
Wei H, Xu Q, Taylor Ii LE, Baker JO, Tucker MP, Ding S-Y: Natural paradigms of plant cell wall degradation. Curr Opin Biotechnol 2009, 20: 330-338.
Alper H, Stephanopoulos G: Engineering for biofuels: exploiting innate microbial capacity or importing biosynthetic potential? Nat Rev Micro 2009, 7: 715-723.
Eiteman MA, Lee SA, Altman R, Altman E: A substrate-selective co-fermentation strategy with Escherichia coli produces lactate by simultaneously consuming xylose and glucose. Biotechnol Bioeng 2009, 102: 822-827.
Fazzini RAB, Preto MJ, Quintas ACP, Bielecka A, Timmis KN, dos Santos VAPM: Consortia modulation of the stress response: proteomic analysis of single strain versus mixed culture. Environ Microbiol 2010, 12: 2436-2449.
Chemier JA, Fowler ZL, Koffas MAG: Trends In Microbial Synthesis of Natural Products and Biofuels. In Advances in Enzymology and Related Areas of Molecular Biology. John Wiley & Sons, Inc; 2009:151-217.
Sonnleitner B, Käppeli O: Growth of Saccharomyces cerevisiae is controlled by its limited respiratory capacity: Formulation and verification of a hypothesis. Biotechnol Bioeng 1986, 28: 927-937.
Matsushika A, Inoue H, Kodaki T, Sawayama S: Ethanol production from xylose in engineered Saccharomyces cerevisiae strains: current state and perspectives. Appl Microbiol Biotechnol 2009, 84: 37-53.
Passoth V, Zimmermann M, Klinner U: Peculiarities of the regulation of fermentation and respiration in the crabtree-negative, xylose-fermenting yeast Pichia stipitis . Appl Biochem Biotechnol 1996, 57–58: 201-212.
Silva JPA, Mussatto SI, Roberto IC, Teixeira JA: Fermentation medium and oxygen transfer conditions that maximize the xylose conversion to ethanol by Pichia stipitis. Renew Energy 2012, 37: 259-265.
Shi NQ, Jeffries TW: Anaerobic growth and improved fermentation of Pichia stipitis bearing a URA1 gene from Saccharomyces cerevisiae . Appl Microbiol Biotechnol 1998, 50: 339-345.
Balagurunathan B, Jonnalagadda S, Tan L, Srinivasan R: Reconstruction and analysis of a genome-scale metabolic model for Scheffersomyces stipitis. Microb Cell Fact 2012, 11: 27.
Caspeta L, Shoaie S, Agren R, Nookaew I, Nielsen J: Genome-scale metabolic reconstructions of Pichia stipitis and Pichia pastoris and in silico evaluation of their potentials. BMC Syst Biol 2012, 6: 24.
Delgenes JP, Escare MC, Laplace JM, Moletta R, Navarro JM: Biological production of industrial chemicals, i.e. xylitol and ethanol, from lignocelluloses by controlled mixed culture systems. Ind Crops Products 1998, 7: 101-111.
Taniguchi M, Tohma T, Itaya T, Fujii M: Ethanol production from a mixture of glucose and xylose by co-culture of Pichia stipitis and a respiratory-deficient mutant of Saccharomyces cerevisiae . J Fermen Bioeng 1997, 83: 364-370.
Srilekha Yadav K, Naseeruddin S, Sai Prashanthi G, Sateesh L, Venkateswar Rao L: Bioethanol fermentation of concentrated rice straw hydrolysate using co-culture of Saccharomyces cerevisiae and Pichia stipitis. Bioresour Technol 2011, 102: 6473-6478.
Goldring ES, Grossman LI, Marmur J: Petite Mutation in Yeast. J Bacteriol 1971, 107: 377-381.
Feist AM, Herrgard MJ, Thiele I, Reed JL, Palsson BO: Reconstruction of biochemical networks in microorganisms. Nat Rev Micro 2009, 7: 129-143.
Price ND, Papin JA, Schilling CH, Palsson BO: Genome-scale microbial in silico models: the constraints-based approach. Trends Biotechnol 2003, 21: 162-169.
Varma A, Palsson BO: Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110. Appl Environ Microbiol 1994, 60: 3724-3731.
Mahadevan R, Edwards JS, Doyle Iii FJ: Dynamic Flux Balance Analysis of Diauxic Growth in Escherichia coli . Biophys J 2002, 83: 1331-1340.
Hjersted JL, Henson MA: Optimization of fed-batch Saccharomyces cerevisiae fermentation using dynamic flux balance models. Biotechnol Progr 2006, 22: 1239-1248.
Stolyar S, Van Dien S, Hillesland KL, Pinel N, Lie TJ, Leigh JA, Stahl DA: Metabolic modeling of a mutualistic microbial community. Mol Syst Biol 2007, 3: 92.
Zhuang K, Izallalen M, Mouser P, Richter H, Risso C, Mahadevan R, Lovley DR: Genome-scale dynamic modeling of the competition between Rhodoferax and Geobacter in anoxic subsurface environments. ISME J 2011, 5: 305-316.
Salimi F, Zhuang K, Mahadevan R: Genome-scale metabolic modeling of a clostridial co-culture for consolidated bioprocessing. Biotech J 2010, 5: 726-738.
Zomorrodi AR, Maranas CD: OptCom: A Multi-Level Optimization Framework for the Metabolic Modeling and Analysis of Microbial Communities. PLoS Comput Biol 2012, 8: e1002363.
Feist A: Bernhard: The growing scope of applications of genome-scale metabolic reconstructions using Escherichia coli. Nat Biotechnol 2008, 26: 659-667.
Weierstall T, Hollenberg CP, Boles E: Cloning and characterization of three genes (SUT1–3) encoding glucose transporters of the yeast Pichia stipitis. Mol Microbiol 1999, 31: 871-883.
Guijarro JM, Lagunas R: Saccharomyces cerevisiae does not accumulate ethanol against a concentration gradient. J Bacteriol 1984, 160: 874-878.
Skoog K, Hahn-Hägerdal B: Effect of Oxygenation on Xylose Fermentation by Pichia stipitis. Appl Environ Microbiol 1990, 56: 3389-3394.
Küenzi MT, Tingle MA, Halvorson HO: Sporulation of Saccharomyces cerevisiae in the Absence of a Functional Mitochondrial Genome. J Bacteriol 1974, 117: 80-88.
Tzagoloff A, Dieckmann CL: PET genes of Saccharomyces cerevisiae. Microbiol Rev 1990, 54: 211-225.
Hanly T, Urello M, Henson M: Dynamic flux balance modeling of S. cerevisiae and E. coli co-cultures for efficient consumption of glucose/xylose mixtures. Appl Microbiol Biotechnol 2012, 93: 2529-2541.
Neves AA, Vieira LM, Menezes JC: Effects of preculture variability on clavulanic acid fermentation. Biotechnol Bioeng 2001, 72: 628-633.
Ligthelm ME, Prior BA, Preez JC, Brandt V: An investigation of d-xylose metabolism in Pichia stipitis under aerobic and anaerobic conditions. Appl Microbiol Biotechnol 1988, 28: 293-296.
Soberón X: Saier Jr MH: Engineering Transport Protein Function: Theoretical and Technical Considerations Using the Sugar-Transporting Phosphotransferase System of Escherichia coli as a Model System. J Mol Microbiol Biotechnol 2006, 11: 302-307.
Kasahara T, Maeda M, Ishiguro M, Kasahara M: Identification by Comprehensive Chimeric Analysis of a Key Residue Responsible for High Affinity Glucose Transport by Yeast HXT2. J Biol Chem 2007, 282: 13146-13150.
Bertilsson M, Andersson J, Lidén G: Modeling simultaneous glucose and xylose uptake in Saccharomyces cerevisiae from kinetics and gene expression of sugar transporters. Bioprocess Biosyst Eng 2008, 31: 369-377.
Taniguchi M, Tanaka T: Clarification of Interactions among Microorganisms and Development of Co-culture System for Production of Useful Substances. 2004.
Postma E, Kuiper A, Tomasouw WF, Scheffers WA, van Dijken JP: Competition for glucose between the yeasts Saccharomyces cerevisiae and Candida utilis. Appl Environ Microbiol 1989, 55: 3214-3220.
Laplace JM, Delgenes JP, Moletta R, Navarro JM: Effects of culture conditions on the co-fermentation of a glucose and xylose mixture to ethanol by a mutant of Saccharomyces diastaticus associated with Pichia stipitis. Appl Microbiol Biotechnol 1993, 39: 760-763.
Albers E, Larsson C, Andlid T, Walsh MC, Gustafsson L: Effect of Nutrient Starvation on the Cellular Composition and Metabolic Capacity of Saccharomyces cerevisiae. Appl Environ Microbiol 2007, 73: 4839-4848.
Heerde E, Radler F: Metabolism of the anaerobic formation of succinic acid by Saccharomyces cerevisiae . Arch Microbiol 1978, 117: 269-276.
Rieger M, Kappeli O, Fiechter A: The role of limited respiration in the incomplete oxidation of glucose by Saccharomyces cerevisiae . J Gen Microbiol 1983, 129: 653-661.
Laplace JM, Delgenes JP, Moletta R, Navarro JM: Alcoholic fermentation of glucose and xylose by Pichia stipitis , Candida shehatae , Saccharomyces cerevisiae ; and Zymomonas mobilis : oxygen requirement as a key factor. Appl Microbiol Biotechnol 1991, 36: 158-162.
Lamping SR, Zhang H, Allen B: Ayazi Shamlou P: Design of a prototype miniature bioreactor for high throughput automated bioprocessing. Chem Eng Sci 2003, 58: 747-758.
Mo M, Palsson B, Herrgard M: Connecting extracellular metabolomic measurements to intracellular flux states in yeast. BMC Syst Biol 2009, 3: 37.
Hanly TJ, Henson MA: Dynamic flux balance modeling of microbial co-cultures for efficient batch fermentation of glucose and xylose mixtures. Biotechnol Bioeng 2011, 108: 376-385.
Hjersted JL, Henson MA, Mahadevan R: Genome-scale analysis of Saccharomyces cerevisiae metabolism and ethanol production in fed-batch culture. Biotechnol Bioeng 2007, 97: 1190-1204.
This work was partially supported by the NSF-sponsored UMass Institute for Cellular Engineering IGERT program (Grant number DGE-0654128) and ReCommunity Recycling.
The authors declare that they have no competing interests.
TH carried out the fermentation and modeling work and drafted the manuscript. MH conceived of the study, and participated in its design and coordination and helped to draft the manuscript. Both authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Hanly, T.J., Henson, M.A. Dynamic metabolic modeling of a microaerobic yeast co-culture: predicting and optimizing ethanol production from glucose/xylose mixtures. Biotechnol Biofuels 6, 44 (2013). https://doi.org/10.1186/1754-6834-6-44