A design-build-test cycle using modeling and experiments reveals interdependencies between upper glycolysis and xylose uptake in recombinant S. cerevisiae and improves predictive capabilities of large-scale kinetic models
2017

Ljubisa Miskovic
Ecole Polytechnique Federale de Lausanne (EPFL)
1015LausanneSwitzerland

Susanne Alff-Tuomala
Ecole Polytechnique Federale de Lausanne (EPFL)
1015LausanneSwitzerland

Keng Cher Soh
Ecole Polytechnique Federale de Lausanne (EPFL)
1015LausanneSwitzerland

Dorothee Barth
Ecole Polytechnique Federale de Lausanne (EPFL)
1015LausanneSwitzerland

Laura Salusjärvi
Ecole Polytechnique Federale de Lausanne (EPFL)
1015LausanneSwitzerland

Juha-Pekka Pitkänen
Ecole Polytechnique Federale de Lausanne (EPFL)
1015LausanneSwitzerland

Laura Ruohonen
Ecole Polytechnique Federale de Lausanne (EPFL)
1015LausanneSwitzerland

Merja Penttilä
Ecole Polytechnique Federale de Lausanne (EPFL)
1015LausanneSwitzerland

Vassily Hatzimanikatis
Ecole Polytechnique Federale de Lausanne (EPFL)
1015LausanneSwitzerland

A design-build-test cycle using modeling and experiments reveals interdependencies between upper glycolysis and xylose uptake in recombinant S. cerevisiae and improves predictive capabilities of large-scale kinetic models

Biotechnol Biofuels
10166201710.1186/s13068-017-0838-5RESEARCH
Background: Recent advancements in omics measurement technologies have led to an ever-increasing amount of available experimental data that necessitate systems-oriented methodologies for efficient and systematic integration of data into consistent large-scale kinetic models. These models can help us to uncover new insights into cellular physiology and also to assist in the rational design of bioreactor or fermentation processes. Optimization and Risk Analysis of Complex Living Entities (ORACLE) framework for the construction of large-scale kinetic models can be used as guidance for formulating alternative metabolic engineering strategies.Results:We used ORACLE in a metabolic engineering problem: improvement of the xylose uptake rate during mixed glucose-xylose consumption in a recombinant Saccharomyces cerevisiae strain. Using the data from bioreactor fermentations, we characterized network flux and concentration profiles representing possible physiological states of the analyzed strain. We then identified enzymes that could lead to improved flux through xylose transporters (XTR). For some of the identified enzymes, including hexokinase (HXK), we could not deduce if their control over XTR was positive or negative. We thus performed a follow-up experiment, and we found out that HXK2 deletion improves xylose uptake rate. The data from the performed experiments were then used to prune the kinetic models, and the predictions of the pruned population of kinetic models were in agreement with the experimental data collected on the HXK2-deficient S. cerevisiae strain.Conclusions:We present a design-build-test cycle composed of modeling efforts and experiments with a glucosexylose co-utilizing recombinant S. cerevisiae and its HXK2-deficient mutant that allowed us to uncover interdependencies between upper glycolysis and xylose uptake pathway. Through this cycle, we also obtained kinetic models with improved prediction capabilities. The present study demonstrates the potential of integrated "modeling and experiments" systems biology approaches that can be applied for diverse applications ranging from biotechnology to drug discovery.

Background

Nature provides abundant plant biomass that can be potentially converted into biofuels as a promising sustainable source of energy [1,2]. The main sugar constituents of biomass are stored in the form of hexoses, primarily glucose, and pentoses, primarily xylose. Saccharomyces cerevisiae, which is widely used in industry for production of a variety of chemicals, including bioethanol from glucose, is a promising candidate for bioethanol production from sugars. However, the native strain of this organism is unable to utilize xylose because of the lack of a mechanism converting xylose into metabolic intermediates [3,4].

In efforts to enable xylose utilization in S. cerevisiae, the wild-type strains have been modified in two ways: (1) by introducing heterologous XYL1 and XYL2 genes from Scheffersomyces stipitis (previously Pichia stipitis) [5] or XYL1 from Candida tenuis and XYL2 from Galactocandida mastotermitis [6,7] that encode for xylose reductase (XR) and xylitol dehydrogenase (XDH), respectively, which enable the two-step transformation from xylose to xylulose; (2) by expressing a heterologous xylose isomerase (XI) from fungi such as Piromyces [8] and Orpinomyces [9], or from Clostridium phytofermentans [10] which converts xylose to xylulose in one step. Comparison of the two engineered pathways in the same background strain revealed that strains with XR/XDH pathway show higher xylose uptake rates under anaerobic conditions and better aerobic growth [11]. In contrast, it appears that strains with XI have higher ethanol yields due to better cofactor balancing as compared to XR/XDH strains.

While the genetically engineered XR/XDH yeast strains successfully acquired the capability of utilizing xylose as a carbon source, a number of factors cause the suboptimal performance of bioconversion from xylose to bioethanol: low transport efficiency of xylose across the cellular membrane through hexose transporters; different cofactor specificities of XR and XDH reactions; high accumulation of xylitol; and excessive activity of the oxidative pentose phosphate pathway [12].

In an attempt to improve xylose uptake rate and obtain a high ethanol yield from xylose fermentation using recombinant S. cerevisiae strains, a number of metabolic engineering efforts have been carried out targeting different enzymes of the metabolic network [13]. Some of the enzymes that have been investigated include hexose transporters (GTR and XTR) [14][15][16][17], xylose reductase (XR) [18,19], xylulose kinase (XK) [20], xylitol dehydrogenase (XDH) [21], and glucose-6-phosphate-1-dehydrogenase (ZWF) [22]. More recently, alteration of the cofactor specificity of XR and XDH has become an important part of efforts to alleviate the cofactor imbalance and to reduce production of xylitol [7, 22,23]. Despite their sophisticated genetic design, these approaches have not achieved yet commercially viable improvements in the utilization of xylose sugar and its subsequent conversion to ethanol. Therefore, it appears that metabolic engineering strategies stemming from a systems-level analysis of cellular pathways can play an important role in creating potential breakthroughs and advancing the research in this area.

Mathematical methods are well suited to address complex metabolic engineering problems that lack extensive experimental data. These methods allow quantification of the impact of individual enzymes on the overall performance of the metabolic network through response analysis of metabolic fluxes with respect to various enzymatic modifications [24]. A mathematical framework named metabolic control analysis (MCA) has been introduced for this purpose [25]. As an outcome of this analysis, the enzymes whose activity has greater impact on desired metabolic flux rates (e.g., xylose uptake) could be identified and they can be considered as high-priority targets for metabolic engineering.

Kinetic models of xylose metabolism have been proposed by different groups previously [26][27][28]. Eliasson et al. built a four-reaction model (XRI, XRII, XDH, and XK) to find a XR/XDH/XK ratio that minimizes xylitol formation during xylose utilization [26]. Parachin et al. constructed models for XR/XDH and XI xylose uptake pathways starting with xylose transporters and ending with XK, and found that in both cases increased XK activity led to improved xylose utilization [27]. Recently, Trausinger et al. integrated data from enzyme activity analyses and quantitative metabolite profiling into an ad hoc reduced model that included XR/XDH pathway, pentose phosphate pathway (PPP), and two lumped reactions describing glycolysis, and used this model for MCA [28].

In this work, we used the Optimization and Risk Analysis of Complex Living Entities (ORACLE) framework [29][30][31][32][33][34][35][36] to construct a population of large-scale kinetic models of glucose-xylose co-utilizing S. cerevisiae that includes XR/XDH pathway, glycolysis, PPP, tricarboxylic cycle (TCA), and electron transport chain (ETC). ORACLE accounts explicitly for mechanistic properties of enzymes and integrates available experimental data, Keywords: Bioethanol, Metabolic control analysis, Large-scale kinetic models, Hexokinase, HXK2 deletion, Xylose utilization, S. cerevisiae network thermodynamics, and physico-chemical constraints of metabolic networks. ORACLE employs Monte Carlo sampling methods to explore the kinetic space of metabolic networks whenever information about kinetic properties is incomplete or missing, and it generates populations of models that are consistent with the experimental information. One of the steps in the ORACLE workflow involves pruning, where the populations of the generated models are further classified into subpopulations with distinct characteristics based on existing or follow-up experiments. The basic principles of ORACLE have been introduced in [32,33,37], and the method was developed and extended in [24, 29-31, 34, 36].

The main goals of this study were to: (1) analyze the impact of the network enzymes on pentose sugar utilization in a genetically recombinant XR/XDH S. cerevisiae strain caused by network-wide couplings and limitations rather than by a substrate competition for hexose transporters; (2) identify the sources of interdependencies between upper glycolysis and xylose uptake rate; (3) engineer genetically a recombinant S. cerevisiae strain with improved xylose uptake capabilities based on hypotheses generated in (2); and (4) use the experimental data acquired on the engineered strain to further improve predictive capabilities of the kinetic models.

For this purpose, we opted for a rational metabolic engineering approach involving a loop between experimental lab work and intensive computational studies. We constructed large-scale kinetic models of S. cerevisiae central metabolism with the XR/XDH xylose uptake pathway wherein we integrated the data collected from a recombinant S. cerevisiae strain growing anaerobically in a medium with a mixture of glucose and xylose. We performed an ORACLE analysis and we identified the enzymes that control the flux through XTR. However, for some of these enzymes such as HXK we could not discern with confidence whether an increase of their activities would lead to a positive or negative effect on the flux through XTR. In order to address this issue, we designed a follow-up computational pruning experiment that helped us classify models into two subpopulations of models with HXK having a positive and a negative control over XTR. We next tested experimentally how changes in hexokinase (HXK) activity would affect xylose uptake rates. The performed experiment with an HXK2 deletion strain showed an increase in the maximum specific consumption rates of xylose. We used this result to further prune our kinetic models, and the predictions of the refined models were consistent with the cultivation data of the engineered HXK-negative recombinant S. cerevisiae strain.

Results

Cultivation and flux analysis of the glucose-xylose utilizing base strain

The strain VTT C-10880 was anaerobically cultivated in a batch reactor using minimal mineral medium with 20 g/l glucose and 50 g/l xylose as the main carbon source for 170 h (for more details see "Methods"). The concentrations (g/l) of glucose, xylose, xylitol, glycerol, acetate, and ethanol were measured (Fig. 1). We observed that the organism started to consume perceptible quantities of xylose between 8 and 10 h, whereas at this time period glucose was at 33% of its starting concentration. Glucose was rapidly consumed, and from T = 15 h onwards the yeast consumed xylose only. At T = 140 h, xylose was not completely consumed yet and xylitol accumulated to approximately 22 g/l.

We formulated a stoichiometric model that describes glucose-xylose co-metabolism under the studied conditions ( Fig. 2; Additional file 1), and we performed the flux analysis using the specific accumulation rates calculated as the output of the kinetic model based on the cultivation data ("Methods"). We used in the flux analysis the calculated rates at time point T = 12 h when glucose and xylose were consumed simultaneously. We then performed thermodynamics-based flux analysis (details in "Methods") to find a flux profile that was consistent both thermodynamically and with the available experimental data. We computed the displacement of reactions from thermodynamic equilibrium, and we found that a majority of reactions were far from equilibrium and no reactions were near equilibrium ( Fig. 2; Additional file 2). This implies that for all the enzymes in the network their control over fluxes and concentrations would depend on their kinetic properties and saturation state [33,35,38]. We assumed that the reaction catalyzed by malic enzyme (MAE) was operating in the direction of decarboxylation because it has been observed that pyruvate carboxylasenegative strains of S. cerevisiae are not capable of growing on glucose indicating that MAE cannot assume the anaplerotic role of pyruvate carboxylase [39]. The upper limit on oxygen uptake rate was set to 2 mmol/gDW h for the following reasons: (1) the medium was saturated with air oxygen at the beginning of the fermentation and its concentration gradually reduced over time; we set the upper bound on oxygen utilization rate to reflect the possibility that there was a residual oxygen activity at the analyzed time point; (2) molecular oxygen appears as a substrate of the used biomass reaction (Additional files 1, 3); at the analyzed time point, the organism was still growing and the model required a small amount of oxygen to describe the growth.

Factors affecting xylose uptake during glucose-xylose co-metabolism

It has been proposed that the reason for slow consumption of xylose in the beginning of the glucose-xylose co-utilization phase is due to substrate competition for high-affinity hexose transporters [40]. In this study, we investigated if the other factors that emerge from the properties of the metabolic network and that are unrelated to the competition for transporters could also impact the xylose uptake. For example, it has been discussed in the literature that availability of the cofactors NAD+ and NADPH in the metabolic network can be one of the limiting factors in xylose utilization [12].

We used the ORACLE methodology ("Methods") to build a population of consistent models. More precisely, we generated ~467,000 system-wide profiles of flux and concentration control coefficients. 1 We then focused on 1 Matlab files containing sets of computed flux and concentration control coefficients together with sampled metabolite concentrations, steady-state reference flux profile and stoichiometric matrix are available at https://doi. org/10.6084/m9.figshare.5028770. the flux control coefficients of the xylose uptake (XTR) and identified the top 18 enzymes with respect to the absolute value of their control coefficients (Fig. 3). Among the 18 enzymes, 7 were from the pentose phosphate pathway (ZWF, XTR, XRI, XRII, XDH, XK, and XLT), 3 were from the upper glycolysis (HXK, PGI, and TPI), 3 were related to energy balance through oxidative phosphorylation and its consumption through maintenance (ASN, NDR, and ATPM), 3 were enzymes controlling the formation of ethanol (ADH1), acetate (ALD), and glycerol-3-phosphate (GPD1), and 2 were enzymes catalyzing exchanges with the extracellular environment (ETOHt and PIt). Downstream enzymes relating to lower glycolysis, citric acid cycle, and glyoxylate shunt were found to have little to almost no control over xylose uptake (Additional file 4).

The positive control over the xylose uptake was distributed over ten enzymes, i.e., ZWF, XTR, XRI, XRII, XDH, XK, XLT, TPI, ALD, and PIt, whereas 5 enzymes, i.e., PGI, ASN, ADH1, GPD1, and ETOHt, had negative control (Fig. 3). XRI had the highest positive control over XTR with C XTR XRI ~ 0.29 indicating that a twofold increase (100%) of XRI activity would result in a 1.29-fold increase  i.e., 0.1 ≤ Γ ≤ 0.9; and (iii) strictly near equilibrium (light red), i.e., 0.9 ≤ Γ < 1. Reactions whose displacements belonged to more than one of these ranges were denoted with (iv) dark red, for 0 < Γ ≤ 0.9; (v) dark green, for 0.1 ≤ Γ < 1; and (vi) gray, for 0 < Γ < 1. The numerical values next to reactions denote flux values (29%) in the xylose uptake rate. ADH1 had the highest negative control over this rate with C XTR ADH1 ~ −0.07 indicating that a twofold increase of ADH1 activity would result in an ~7% decrease in the xylose uptake rate.

For three enzymes, i.e., HXK, ATPM, and NDR, the distributions of the control coefficients, C XTR HXK , C XTR ATPM , and C XTR NDR , were spread around the zero value (Fig. 3), and with the available data we were unable to predict with certainty if the manipulation of these enzymes would affect the xylose uptake rate in a positive or negative manner. This implied that we might have two distinct subpopulations of models with one population predicting a negative control of HXK over XTR (C XTR HXK < 0), and the other predicting a positive control (C XTR HXK > 0). In addition, we observed a rather strong negative correlation between the values of C XTR HXK and the values of both C XTR ATPM and C XTR NDR . We computed the correlation between these control coefficients and the Pearson coefficient for the correlation of C XTR HXK and C XTR ATPM was −0.68, the one of C XTR HXK and C XTR NDR was −0.67, whereas the one for the correlation of C XTR ATPM and C XTR NDR was +0.76.

Factors affecting ATP consumption and generation

Interestingly, the enzymes participating in the generation and consumption of ATP appeared to have considerable control over xylose uptake, even though they have been traditionally overlooked when the system was studied from the perspective of carbon flow.

Therefore, we further investigated how changes in the activities of the enzymes would affect the levels of cytosolic ATP (Fig. 4). We used the population of ~467,000 kinetic models to identify the top 16 enzymes controlling the levels of this metabolite. Out of the top 16 enzymes, three were from the upper glycolysis (HXT, HXK, and TPI), two were from the pentose phosphate pathway (RKI and XDH), two were directly related to ATP consumption and production (ADK and ATPM), two were from acetate metabolism (ACS and ACE), five were exchange reactions with the extracellular environment (ETOHt, NH 4 t, SO 4 t, PIt, and CO 2 t), and one was COH which is catalyzing the reversible conversion of carbonic acid to CO 2 and water. HXK had the primary negative control on the levels of ATP with C Considering that ATP impacts a large number of reactions in the network including HXK and that HXK had a strong control over ATP implied that there was a strong coupling through ATP in the network. This together with the hypothesis that an increase in the activity of HXK would decrease the xylose uptake rate was intriguing, and we decided to verify this hypothesis by performing an experiment on a strain with attenuated HXK activity.

HXK2 deletion improves xylose uptake

Saccharomyces cerevisiae has three enzymes phosphorylating glucose, HXK1, HXK2, and GLK1, of which HXK2 is the predominant enzyme under consumption of glucose.

HXK2 is expressed when growing on glucose, but when growing on non-fermentable carbon sources HXK2 expression is repressed and the expression of HXK1 and GLK1 de-repressed [41]. As HXK2 contributes to the predominant hexokinase activity on glucose and is the main enzyme conveying catabolite repression of these two hexokinases, we constructed a recombinant S. cerevisiae HXK2-deficient strain (see Methods). The engineered strain indeed displayed a clear improvement of the xylose uptake rate with values of specific rate of xylose consumption consistently higher than the ones of the base strain with the maximum values showing improvement by approximately 60% (Fig. 5a). The increase in the xylose consumption led also to an increase in the production of xylitol with a 60% increase in the maximum specific xylitol production rate (Fig. 5c). Conversely, deletion of the HXK2 gene resulted in (a) reduced maximal glucose specific consumption rate by 26% with the depletion of glucose occurring at approximately T = 36 h compared to T = 16 h of the base strain (Fig. 5b); and (b) reduction of the maximal specific accumulation rate of ethanol by 13% (Fig. 5d).

The computed control coefficients for specific glucose uptake, C HXT HXK , and for specific productivity of ethanol, C ADH1 HXK , were in qualitative agreement with the experiments, i.e., the trends of the expected metabolic responses and the relative magnitudes of these control coefficients were in agreement with the experiments. The computed mean values of C HXT HXK and C ADH1 HXK were, respectively, 0.18 and 0.095 (Fig. 6b, d), i.e., a twofold decrease in the activity of HXK would lead to 9 and 4.75% decrease in the corresponding specific productivities.

Pruning the models

Based on the experimentally obtained information, we rejected all sampled models with C XTR HXK > 0 and ended up with a population of ~204,000 models. The mean value of C XTR HXK shifted to a value of −0.11 after the pruning compared to 0.02 before the pruning (Fig. 6a). The value of −0.11 indicated that a twofold increase in the activity of HXK would result in the reduction of the specific xylose uptake by 11%. Interestingly, the pruning based on C XTR HXK did not have a significant effect on the distribution of most of the other control coefficients with the exception of C XTR ATPM and C XTR NDR . Before the pruning, these two control coefficients had distributions spread around the zero value with the means of −0.02 and −0.05, respectively (Fig. 6a) Fig. 6 Control coefficients of the original (purple) and refined (green) models for xylose uptake rate, XTR (a), glucose uptake rate, HXT (b), xylitol production, XLT (c), and alcohol dehydrogenase, ADH1 (d). The bars represent the mean values of the control coefficients, and the error bars denote the 1st and the 3rd quartile of the control coefficients with respect to their mean value, i.e., 50% of the samples closest to the mean value are within the error bars. With # are denoted the enzymes whose distributions of control coefficients were spread around zero in the original (purple) models that the negative correlation between the values of C XTR HXK and the values of both C XTR ATPM and C XTR NDR that we discovered on the whole population of 467,000 models still existed on the pruned population. We computed the correlations for the pruned population of control coefficients and the Pearson coefficient for the correlation of C XTR HXK and C XTR ATPM was −0.56, the one of C XTR HXK and C XTR NDR was −0.58, whereas the one for the correlation of C XTR ATPM and C XTR NDR was +0. 73.

In silico predictions of the pruned control coefficients were consistent with all the experimentally observed effects of HXK2 deletion on the behavior of the engineered strain. More precisely, in qualitative agreement with the performed experiments the refined population of control coefficients predicted that HXK2 deletion would (a) decrease the glucose uptake, i.e., the mean of computed C HXT HXK was ~0.14 ( Fig. 6b); (b) reduce the ethanol production, i.e., the mean C ADH1 HXK was ~0.05 (Fig. 6d); and (c) increase the xylitol production, i.e., the mean C XLT HXK was ~−0.14 ( Fig. 6c). The refined models also consistently predicted that the HXK2 deletion would increase the flux through xylulose kinase (XK) with the mean C XK HXK of ~−0.08 (Additional file 5). Indeed, computed as the difference between the experimental values for xylose consumption and xylitol production fluxes (Fig. 4a, c), the estimated flux through XK was higher in the HXK2 deletion strain than in the reference strain.

Network-wide effects of altered HXK activity

Metabolic control analysis indicated that HXK2 deletion would cause significant rewiring of the metabolic network. Reactions from the glycolysis, pentose phosphate pathway with the exception of xylose uptake pathway, oxidative phosphorylation with the exception of lactate-related reactions, glycerol metabolism, acetate and ethanol production, and exchange with the extracellular environment would reduce their fluxes compared to the parental strain ( Fig. 7; Additional file 4). For instance, a twofold decrease of the activity of HXK would reduce the fluxes in the glycolysis from 5.6% (PYK) to 14% (HXK) and the fluxes in the pentose phosphate pathway from 4% (TKL2) to 141% (RPE). In contrast, decreasing the activity of HXK would result in an increased activity of the xylose uptake pathway, glyoxylate shunt, tricarboxylic acid cycle, and ATPM which manages a series of energy-requiring events inside the cell without leading to net formation of biomass including futile cycles and maintenance of the proton gradient and electrical potential (Fig. 7). A twofold decrease of activity of HXK would increase fluxes through XTR (by 5%), XRI (by 7%), XLT (by 6.6%), XDH (by 4%), and XK (by 4%) in the xylose uptake pathway, fluxes through MDHc (by 30%) and ICL, MLS, and ACOc in the glyoxylate shunt (by 43%), fluxes through KGD, LSC, SCD, and FUM in the TCA cycle (by 6%), and ATPM (by 18%).

Furthermore, a twofold reduction in the HXK activity would reduce ATP turnover (by ~3.2%) where all ATPproducing reactions (PGK, PYK, and AAC) and most of ATP-consuming reactions (HXK, PFK, PYC, ADK, ACS, PEPCK, and biomass) would have reduced fluxes (Fig. 7; Additional file 4). The only exceptions would be XK and ATPM, whose fluxes would increase (by 4.2 and 18%, respectively).

Network-wide correlations of HXK, NDR, and ATPM

Interestingly enough, an increased activity of NDR and ATPM, whose control over XTR negatively correlated with the one of HXK ("Results"), would have similar effects on the metabolic network as the HXK2 deletion (Additional files 4, 6, 7). Notable differences in the metabolic responses would be in (1) a reduced flux through XK and XDH (by 15%) with a twofold ATPM increase, and by 21% with a twofold increase of NDR activity; this implies that metabolic engineering strategies involving NDR and ATPM would improve xylose uptake but the carbon taken up would be used almost exclusively for xylitol production; (2) reduction in the TCA fluxes (ACO, IDH, KGD, LSC, SCD, FUM, and MDH) from 33 to 94%, and reduced ATPM by 11%, with a twofold increase of the NDR activity. The twofold increase in ATPM and in NDR activity would also result in reduced turnover of cytosolic ATP by 9 and 17%, respectively (versus 3.2% of the one in the case of HXK) and reduced turnovers of cytosolic NADH (by 11% for ATPM and 14% for NDR versus 3.4% for HXK) and NADPH (by 41% for ATPM and 31% for NDR versus 15% for HXK).

The negative correlation between HXK and ATPM was expected as the reduced HXK activity results in reduced glycolytic fluxes and as a consequence in reduced ATP production and availability. A similar effect on ATP can be reproduced by increasing non-growth-associated ATP maintenance (ATPM). In contrast, the correlation between NDR and both HXK and ATPM was more surprising. Our analysis revealed that this correlation is due to the stoichiometric coupling between NDR and the growth in the analyzed physiological condition (Additional file 8). Specifically, we performed the flux variability analysis (FVA), sampled the space of steady-state fluxes, and found that the growth is perfectly coupled with NDR. For the smallest feasible values of the flux through NDR (0.016 mmol/gDW/h), the growth was small (0.002 1/h). In contrast, for the highest feasible values of NDR flux (0.48 mmol/gDW/h), the growth increased almost 100-fold to 0.19 1/h (Additional file 8). The growth increase resulted in an increase of the growth-associated ATP requirements, and the

Discussion

HXK2 is the gene that has a main role in glucose phosphorylation, and it could be expected for a strain that lacks this gene to exhibit decreased glucose consumption levels. In the performed experiment, HXK2 deletion decreased the maximal glucose consumption rate by 26% compared to the parental strain. However, for predicting the responses of the other fluxes in the network like xylose uptake rate, a systems-oriented tool such as MCA is needed. HXK2 deletion would have been an unlikely target for improving xylose consumption without the MCA analysis. Indeed, the pruned MCA models successfully predicted an increase in xylose consumption and xylitol production and a decrease in glucose consumption and ethanol production of the HXK2-deficient strain.

Considering that we modeled the transporters of glucose and xylose as enzymes without competitive inhibition and that we neglected the HXK2 catabolic repression, MCA results implied that there are factors other than competition for the transporters or glucose repression that could negatively affect the xylose uptake. MCA suggested that decreasing HXK activity would increase the concentration of cytosolic ATP (Fig. 4) and the rate of xylulose kinase (XK) (Additional file 5), thus allowing improved xylose uptake (Fig. 7).

Previous experimental studies investigated some of the enzymes identified here for improved xylose uptake rate (Fig. 3). Jeppsson et al. [22] reported that a reduced flux through the oxidative branch of pentose phosphate pathway (PPP) causes reduced xylose uptake rate in a xylose-only fed recombinant S. cerevisiae strain. Consistent with [22] where the ZWF1 gene was deleted to reduce the NADPH-producing flux, our predictions indicate that lowering the activity of glucose-6-phosphate-1-dehydrogenase (ZWF) would result in reduced xylose uptake rate. Our predictions are also in agreement with two experimental studies regarding the effects of xylose reductases (XRI and XRII) [18] and xylitol dehydrogenase (XDH) [21] on xylose uptake rate.

Our predictions additionally suggest that increasing the PGI activity would result in a reduced xylose uptake rate. During glucose-xylose fermentations studied here, PGI converts glucose 6-phosphate to fructose 6-phosphate, and therefore a higher PGI activity results in a reduced flux through the oxidative branch of PPP. This is consistent with the experimental study [22], since in batch fermentations with xylose as the sole carbon source performed in [22] PGI converts fructose 6-phosphate to glucose 6-phosphate, and a lower PGI activity reduces the flux through the oxidative branch of PPP.

Toivari et al. [20] obtained improved xylose uptake with a CEN.PK2 strain with an increased xylulose kinase (XK) activity both in xylose-only and glucose-xylose anaerobic batch fermentations. These experimentally observed effects of XKS1 overexpression are in agreement with the quantitative predictions (Fig. 3). In contrast, XKS1 overexpression in a different strain background (including overexpressed non-oxidative PPP genes, gre3 deletion, and xylB originating from Escherichia coli) and in an anaerobic xylose-only batch fermentation resulted in slightly reduced xylose uptake [27]. These findings are in line with the results from other studies on XK role in xylose uptake for S. cerevisiae which indicate that only a fine-tuned XK overexpression leads to improved xylose fermentation [42][43][44][45].

Furthermore, in a recent study, Trausinger et al. [28] performed MCA of xylose assimilation of S. cerevisiae strain grown under anaerobic conditions on xylose as a sole carbon source. The computed flux control coefficients for xylose uptake with respect to XR, XDH, XK, and ZWF in [28] are consistent with the ones found in the current study.

While the refined models successfully predicted the trends of the metabolic responses of the HXK2 deletion, the magnitudes of the responses predicted by the refined control coefficients did not accurately match the experimentally observed ones. These discrepancies can originate from any of the following reasons or their combination: (a) control coefficients represent a log-linear approximation of the metabolic responses upon small changes around the studied steady state in an enzyme activity, and they typically provide a good ranking of enzymes that control the studied fluxes or concentrations. In this study, we performed a gene deletion that is a major perturbation of a metabolic network; (b) it has been shown that in a HXK2-deficient S. cerevisiae strain the transporters HXT2, HXT4, HXT6, HXT7, HXT8, and HXT16 have higher gene expression levels [46]. Of these hexose transporters, especially the high-affinity transporters HXT2, HXT4, and HXT7 are active in xylose transport [14,47]. The absence of HXK2 repressive effect on the hexose transporters involved in the transport of xylose might be the reason for the higher than predicted consumption rate of xylose in the HXK2-deficient strain; and (c) the population of kinetic models was computed for the time point T = 12 h when glucose and xylose were co-utilized, whereas the maximal consumption and production rates compared in two strains were measured at T = 3.5 h when glucose was predominantly utilized.

Finally, although the estimated values for specific rates of xylose and glucose conversion and ethanol and xylitol production may be affected by variations in experimental measurements and by a bias in data smoothing, there was a clear and significant improvement in specific xylose conversion rate upon HXK2 deletion and it was predicted by the model. ORACLE is designed to handle the uncertainty in data, and the populations of sampled models capture the trends in the metabolic responses upon HXK2 deletion.

Conclusions

We described here a design-build-test cycle including the steps of modeling, identifying targets for engineering, testing the gene modification in a engineered strain, and further refining the models. More specifically, we used MCA to explore potential improvements of xylose utilization in glucose-xylose co-fermenting S. cerevisiae. From this analysis, we postulated a hypothesis that the HXK enzyme has negative control over specific xylose uptake rate and we performed an experiment to verify it. Although the effects of HXK2 deletions on glucose repression have been studied before [48][49][50], a deletion of HXK2 has not been studied as a means to improve xylose consumption in glucose-xylose co-fermenting S. cerevisiae strains. The influence of carbon catabolite repression on xylose metabolism in S. cerevisiae strains through deletion of MIG1 and MIG2 genes was studied previously; however, no improvement of xylose consumption during batch cultivations was observed [51].

The performed in silico analysis demonstrated that modifying HXK activity would substantially rebalance the network fluxes, and that these effects strongly correlate with the ones of NDR and ATPM. We discovered that the correlation of effects of HXK, ATPM, and NDR on the whole network is an emerging property of the stoichiometric coupling through ATP.

We engineered an HXK2-deficient strain and the cultivation results validated the postulated hypothesis. We used the experimental information to prune the population of our kinetic models and the refined models were in qualitative agreement not only with the response of XTR to HXK2 deletion but also with the responses of HXT, ADH1, XLT, and XK. Pruning is one of the important features of ORACLE as it allows a straightforward refinement of models based on experimental information about the metabolic responses to gene modifications and it improves the efficiency of the design-build-test cycle. This is an advantage with respect to the conventional parameter estimation methods where the integration of this kind of information is challenging.

The present study aims to improve the xylose uptake rate during a mixed glucose-xylose co-utilization and to elucidate a complex interplay between upper glycolysis and xylose uptake. It is a step toward obtaining strains with economically viable yields and specific productivity of ethanol from both glucose and xylose as carbon sources. The next step in this direction is to identify targets that will allow for recovered glucose uptake and ethanol specific productivity while retaining the improved xylose uptake of the HXK2 deletion strain. Our models suggest that one of such targets might be TPI which can also be used together with XDH to improve the yield of ethanol from both glucose and xylose.

Methods

Strains and plasmids

The yeast strain used in the study was VTT C-10880, i.e., [CEN.PK113-1A (URA3. HIS3. LEU2. TRP1. MAL2-8 c .SUC2), ura3::XYL1 XYL2, xks1::XKS1]. XYL1 and XYL2 genes of S. stipitis were chromosomally integrated into URA3 locus under PGK1 and ADH1 promoters, respectively. The integration cassette was constructed as described earlier [20]. XKS1 of S. cerevisiae was integrated into XKS1 locus as described below.

The XKS1 expression cassette pADH1m-XKS1-tADH1 was released as a BamHI fragment from plasmid pMLV14 (BluescribeM13) and ligated into the bacterial plasmid pMLV24. The plasmid pMLV24 contains loxP-S. cerevisiae MEL5 (encoding α-galactosidase)-loxP marker cassette with 60-bp flanking regions for targeting to the XKS1 locus. The XKS1 flanking regions were from nucleotides −250 to −192 and from nucleotides 1801 to 1861, where numbers are relative to the nucleotide A in the XKS1 ATG start codon. The BamHI cloning site was included in the XKS1 5′ flanking sequence. The resulting plasmid, containing XKS1, was named pMLV46. The expression cassette for XKS1, together with the loxP-MEL5-loxP marker, was released from pMLV46 with PstI and SpeI and introduced into yeast cells by transformation. Blue-colored, MEL5 (α-galactosidase)expressing yeast colonies were collected from agar-solidified YP containing 2% w/v d-galactose, supplemented with 5-bromo-4-chloro-3-indolyl-α-d-galactopyranoside (X-α-Gal, 40 μg/ml). To remove the MEL5 marker cassette from the yeast chromosome, the transformant was retransformed with a plasmid pSH47 [52], expressing the Cre recombinase. Integration of the expression cassette pADH1m-XKS1-tADH1 into the XKS1 locus in the genome was verified with PCR. The yeast strain generated was named VTT C-10880.

Plasmid pKK27-1 was constructed by releasing KanMX4 from pFA6a-KanMX4 [53] by BamHI + EcoRI digestion. KanMX4 was further ligated into YEplac195 [54] digested with BamHI + EcoRI resulting in pKK27-1.

HXK2 of VTT C-10880 was deleted by integration of KanMX4 into HXK2 locus. KanMX4 deletion cassette was amplified by PCR from pKK27-1 plasmid using forward 5′ TCTGTTTAGCTTGCCTCGTC and reverse 5′ CACTGGATGGCGGCGTTAGTA primers with the homologous flanking sequences to HXK2 locus (−27 to −33 for the forward primer and +1 to +60 for the reverse primer). The resulting PCR product was integrated into the genome of VTT C-10880 by homologous recombination by following the transformation procedure of Woods and Gietz [55]. Transformants were selected on YPD medium containing 200 µg/ml G418. The deletion of HXK2 was verified by PCR using primer pairs 5′GGTTG-TAGGAATATAATTCTC3′ and 5′CAACGCTACCTTT-GCCATGT3′ and 5′ATCCTGATATGAATAAATTG3′ and 5′CATGTTCACATAAGTAGAAAAAGGGCACC3′.

Bioreactor cultivations

The strain VTT C-10880 and the HXK2 deletion strain were studied for 170 h in 1.5-l batch cultures in Biostat B bioreactor (Sartorius, Göttingen, Germany) at pH 5 and 30 °C using agitation at 200 rpm, under anaerobic conditions with 0.5 l/min nitrogen flushed to reactor headspace. The HXK2 deletion strain was cultivated in duplicate experiments. The inocula were prepared by transferring the cells from an YPD plate to 25 ml mineral medium in 100-ml Erlenmeyer flasks and incubated on a plate shaker (150 rpm, 30 °C) overnight. The inocula were further transferred into 75 ml mineral medium in 250-ml Erlenmeyer flasks and incubated (150 rpm, 30 °C) for 4-6 h. The inocula were centrifuged (2000 rpm, 4 °C, 5 min) and the cells were resuspended into the growth medium lacking the source of carbon. Initial cell dry weight in bioreactor was 0.15 g/l. The minimal mineral medium [56] contained 20 g glucose/l, 50 g xylose/l, 22 mg uracil/l, 5 g (NH 4 ) 2 SO 4 /l, 3 g KH 2 PO 4 /l, 0.5 g MgSO 4 ·7H 2 O/l, 15 mg C 10 H 14 N 2 Na 2 O 8 ·2H 2 O (EDTA)/l, 4.5 mg ZnSO 4 ·7H 2 O/l, 0.84 mg MnCl 2 ·2H 2 O/l, 0.3 mg CoCl 2 ·6H 2 O/l, 0.3 mg CuSO 4 ·5H 2 O/l, 0.4 mg Na 2 MoO 4 ·2H 2 O/l, 4.5 mg Na 2 MoO 4 ·2H 2 O/l, 3.0 mg FeSO 4 ·7H 2 O/l, 1.0 mg H 3 BO 3 /l, 0.1 mg KI/l, 0.05 mg biotin/l, 1.0 mg Ca-Pantothenate/l, 5 mg nicotinic acid/l, 25 mg myoinositol/l, 1.0 mg thiamine HCl/l, 1.0 mg pyridoxol HCl/l, and 0.2 mg p-aminobenzoic acid/l.

Sampling and analyses

Samples were taken frequently with a sampling robot (Medicel, Helsinki, Finland) into a cold bath where the samples froze at the temperature of −30 °C. The frozen samples were thawed, OD600 and cell dry weight (CDW) measurements were performed, and supernatants were collected for HPLC and CE measurements. Concentrations of xylose, glucose, glycerol, xylitol, ethanol, and acetate were measured with a Waters Alliance 2690 HPLC system with a Waters 2410 RI detector (Waters Corporation, Milford, USA). Aminex HPX-87H column and Aminex fast acid analysis column (BioRad, USA) were used in series with 5 mM H 2 SO 4 (Merck Titrisol) with a flow rate of 0.6 ml/min to separate the analytes at 55 °C. Concentrations of lactic acid, succinic acid, and malic acid were measured using a P/ACE MDQ capillary electrophoresis system (Beckmann Coulter Inc., Fullerton, USA) equipped with a PDA detector [57]. Concentrations of CO 2 , O 2 , and ethanol in the exhaust gas of the bioreactor were measured using an Innova gas analyzer (Innova Air Tech Instruments A/S, Ballerup, Denmark).

Quantitative analysis of fermentation measurements

A kinetic model was applied to smooth the measurement data to reduce the noise for better estimation of the calculated continuous variables such as specific accumulation rates. Performing these calculations directly from the measured data would provide data too noisy for interpretations. The dynamic system of anaerobic yeast cultivation was described by a set of non-linear ordinary differential equations (ODEs), with the system components comprising concentrations of biomass, the substrates glucose and xylose, the products ethanol, glycerol, xylitol, acetate, lactate, succinate, and malate, and the offgas CO 2 .

The consumption of both substrates, glucose and xylose, was described by Monod kinetics with an initial lag phase. Concentrations are written as c, the subscript s stands for either of the substrates, glucose and xylose, respectively, and the subscript x for biomass. The maximum specific substrate uptake rate is denoted by Y s , and the monod dynamics are characterized by the parameter k m,s . The lag phase is determined by the time constant τ.

Growth of the biomass is determined by the rate of substrate consumption, and thus with the specific maximum growth rates µ glu and µ xyl on glucose and xylose, respectively.

The products, i.e., ethanol, xylitol, glycerol, acetate, and minor acids, as well as the off-gas CO 2 , denoted by the subscript p, are released proportionally to substrate consumption with a yield parameter Y s,p specific to the substrate and the respective product.

Ethanol is evaporating from the broth during the cultivation in non-negligible amounts. In order to fit the The parameters of the model were fit to the available measurement data by means of a maximum likelihood estimator in two subsequent steps. In a first step, the parameters for the coupled system containing the ODEs for glucose, xylose, and biomass concentrations were identified, and in a following step all the remaining parameters were estimated (Additional file 9). Parameter estimations and simulations were performed with Matlab, version 7.10.0 (R2010a) (Mathworks, Natick, MA, USA).
dc s dt = −Y s c x c s k m, s + c s 1 − e − t τ dc x dt = −µ glu dc glu dt − µ xyl dc xyl dt dc p dt = −Y glu, p dc glu dt − Y xyl

Mathematical models of xylose-glucose utilizing metabolic network Stoichiometric model

The model of the S. cerevisiae recombinant xylose-glucose co-utilizing metabolic network consists of 102 reactions and 96 intracellular metabolites distributed over cytosolic and mitochondrial compartment and extracellular environment (Fig. 2). The model contains XR/XDH xylose assimilation pathway [12]. The metabolites have been categorized as cytosolic or mitochondrial according to the physiological information on their cellular compartmentalization. All reactions in the model are completely balanced, i.e., reactions are atomically balanced with respect to carbon, protons, nitrogen, phosphorus, etc.

Although a family of hexose transporters, Hxt1p-Hxt17p, Gal2p, and Snf3p [40,58,59], facilitates the xylose assimilation in S. cerevisiae, in this study the xylose and glucose transport were modeled as independent reactions catalyzed by separate enzymes. This way, we are able to investigate the relations between upper glycolysis and xylose uptake pathway that are not caused by the competition for transporters of these two pathways.

Kinetic model

For each of the reactions within the studied metabolic network, a kinetic mechanism has been assigned based on the available literature data. The kinetic mechanisms used include reversible Michaelis-Menten kinetics, ordered Bi-Bi, Bi-Ter, and Ter-Bi [60]. Whenever the exact kinetic mechanism was not known, generalized approximations of enzymatic mechanisms such as generalized reversible Hill [61] or convenience kinetics [62] dc eth, mess dt = dc eth dt − ec eth e = a − bv.

have been used. The kinetic model includes also a model of allosteric regulation for the phosphofructokinase (PFK) reaction. The mechanism for this reaction is modeled as Hill kinetics with the Hill coefficient h = 4, where adenosine monophosphate (AMP), acting as an external activator, and adenosine triphosphate (ATP), acting as an inhibitor, bind to the same site. Detailed account about the used kinetic mechanisms is given in the Additional files 3 and 10.

Displacement from thermodynamic equilibrium

For any metabolic reaction, the displacement from the thermodynamic equilibrium can be expressed as a function of the equilibrium constant, K eq , defined as the ratio between the products of the forward and backward rate constants, and the concentration of the involved substrates, S i , and products, P j , i.e.

The equilibrium displacement, Γ , is related to the Gibbs free energy difference, G, through the following thermodynamic equation [38,63]:

with R being the ideal gas constant and T the standard temperature. We assume that the reactions produce spontaneously the products P j , i.e., �G < 0 and consequently Γ < 1. When a reaction gets close to its thermodynamic equilibrium, Γ is approaching 1. At equilibrium, Γ = 1 and where G • denotes the standard Gibbs energy of reaction.

Efficient building of large-scale kinetic models

The optimization and risk analysis of complex living entities (ORACLE) methodology used in this contribution is based on the MCA paradigm. ORACLE is composed of several successive computational procedures, each of them bringing a new level of information about metabolic networks (Fig. 8). Conceptually, ORACLE is organized as follows:

Level 1

We first determine the stoichiometry of the metabolic network using biochemical data or genome reconstructions [64][65][66][67][68][69][70]. For organisms with incompletely sequenced genomes, the part of the stoichiometry concerning missing pathways or compounds can
Γ = 1 K eq P j S i . �G = RT ln Γ P j S i eq = −RT ln K eq = G • ,
be hypothesized. We proceed further by integrating the available information from metabolomics and fluxomics analyses [71][72][73], and using the estimates of the standard free energy of reactions [74][75][76][77] we perform the thermodynamics-based flux balance analysis (TFA) [78,79].

Integration of thermodynamics at this stage allows us to eliminate thermodynamically unfeasible reactions, to establish reaction directionality, and to compute thermodynamically feasible flux profiles.

Level 2

Accurate and consistent determination of the space of metabolite concentrations plays a significant role in appraising all viable kinetic models sharing the same steady-state flux profile for two reasons: (1) the metabolite concentration levels intrinsically affect the Gibbs free energy difference; consequently, these have to be consistent with the reaction feasibility and directionality of computed flux profiles;

(2) the metabolite concentration levels are one of the factors determining local stability of the kinetic models [32]. Uniform random sampling techniques are commonly used to investigate the space of metabolite concentrations [80]. However, the volume reduction of the thermodynamically feasible metabolite concentration space becomes so important in higher dimensional metabolite networks that acceptance-rejection methods are prohibitively inefficient in generating feasible random samples [81]. Furthermore, the thermodynamically feasible metabolite concentration space is not in the form of a parallelotope (nor it is in the form of a simplex or a hypersphere), and consequently the transformation techniques for the generation of random samples are not applicable [82,83]. Therefore, we recur to Monte Carlo Markov Chain methods to generate random samples with distributions approaching the uniform distribution asymptotically. When experimental measurements or estimates from experimental data under similar physiological conditions are available for some metabolites in the network [84][85][86], we integrate these data as bounds on metabolite levels for the sampling.

We use the levels of allowable metabolite concentrations along with the standard Gibbs energy of reaction, G • , to compute the equilibrium displacements of thermodynamically feasible reactions, Γ . Knowledge about equilibrium allow us, before even going to kinetics, to discern the reactions that are near thermodynamic equilibrium, i.e., the reactions that are not potential targets for metabolic engineering. The generated pairs of the metabolite concentration levels and the equilibrium displacements will subsequently be called 'equilibrium displacement profiles' as they are intrinsically inseparable.

Level 3

We integrate the kinetic properties of enzymes at this level. Available kinetic data from literature are incorporated. Whenever kinetic information about enzymes is incomplete or missing, we sample to recover the missing kinetic data. This is performed either through sampling of enzyme states [31], or through sampling of the degree of saturation of the enzyme active site [32].

At this level, we also perform consistency checks and pruning. In these tests, we evaluate the local stability of the resulting kinetic models and the consistency with the experimental information. As a result of the computational procedures performed at this level, we obtain the populations of elasticities consistent with thermodynamically feasible metabolite concentrations.

Level 4

Using the information gathered in previous levels, we compute populations of control coefficients and store them. Only the control coefficients stemming from stable systems are retained. If at a later time experimental evidence about the control coefficients becomes available, e.g., the response of a metabolic flux to an increase of activity of an enzyme becomes known, we prune the samples of the population with inconsistent control coefficients.

Level 5

We perform statistical analysis and data mining on the populations of control coefficients in order to quantify the importance of the enzymes in possible metabolic engineering strategies [87][88][89]. Hypotheses about possible couplings within the metabolic network are postulated at this level as well.

All the above-mentioned levels include visualization of the obtained knowledge about the metabolic network.

Turnover control coefficients

The turnover of a metabolite, tMet, can be defined as the sum of fluxes that produce, or alternatively as the sum of fluxes that consume, the metabolite: tMet = n i=1 v i . Here, n represents the number of producing (or consuming) fluxes v i . The turnover control coefficients can then be calculated as
C tMet q = q tMet dtMet dq = n i=1 v i tMet C v i q ,
where q denotes the set of system parameters, such as the enzyme activity and extracellular metabolite concentrations, and C v i q is the flux control coefficient of the flux v i .  Additional file 2. Thermodynamic displacement of reactions, Γ. Reactions can operate: (i) strictly far from thermodynamic equilibrium, i.e. 0 < Γ ≤ 0.1; (ii) with the middle displacements, i.e. 0.1 ≤ Γ ≤ 0.9; and (iii) strictly near equilibrium, i.e. 0.9 ≤ Γ < 1. Reactions whose displacements spanned more than one of these ranges were denoted with: (iv) I + II, for 0 < Γ ≤ 0.9; (v) II + III, for 0.1 ≤ Γ < 1; and (vi) I + II + III, for 0 < Γ < 1.

Additional file 3. List of reactions, their substrates and products, together with the used kinetic mechanisms (Additional file 10).

Additional file 4. Sheet "Flux Control Coefficients": (i) Control coefficients of the xylose uptake (XTR) for a set of most important enzymes; (ii) Correlation of the distributions of the XTR control coefficients for HXK, ATPM and ZWF; (iii) The effects of manipulations of HXK, ATPM and NDR on metabolic fluxes. Sheet "Concentration Control Coeffs": Control coefficients of the cytosolic ATP, cytosolic NADH and cytosolic NADPH during glucose-xylose co-utilization. Sheet "Turnover Control Coefficients": (i) Turnover control coefficients of cytosolic ATP, mitochondrial ATP, cytosolic NADH, cytosolic NADPH, cytosolic protons and mitochondrial protons for HXK, ATPM and NDR; (ii) Flux control coefficients of producing and consuming reactions of cytosolic ATP, mitochondrial ATP, cytosolic NADH, cytosolic NADPH, cytosolic protons and mitochondrial protons for HXK, ATPM and NDR.

Additional file 5. Control coefficients of the flux through xylulose kinase (XK) during glucose-xylose co-utilisation for the original (purple) and refined (green) kinetic models. The bars represent the mean values of the control coefficients, and the error bars denote the 1st and the 3rd quartile of the control coefficients with respect to their mean value, i.e. 50% of the samples closest to the mean value are within the error bars.

Fig. 1
1Substrates and products time evolutions of glucose-xylose co-utilization obtained for the VTT C-10880 strain growing in a bioreactor under anaerobic conditions: glucose (blue squares), xylose (red diamonds), xylitol (green triangles), glycerol (cyan pluses), acetate (magenta crosses), and ethanol (orange asterisks) Fig. 2 Schematic representation of the VTT C-10880 metabolism together with the displacement of the reactions from thermodynamic equilibrium. Reactions can operate (i) strictly far from thermodynamic equilibrium (light green), i.e., 0 < Γ ≤ 0.1; (ii) with the middle displacements (blue),

HXK
~ −1.07 indicating an important increase in the levels of cytosolic ATP upon reduction in the activity or deletion of HXK2. HXT also had a strong negative control over the levels of cytosolic ATP (C atp c HXT ~ −0.93). The strongest positive control over this metabolite was surprisingly in CO 2

Fig. 3 Fig. 4
34Control coefficients of the xylose uptake during glucosexylose co-utilization. The bars represent the mean values of the control coefficients through xylose transporters (XTR). The error bars denote the 1st and the 3rd quartile of the control coefficients with respect to their mean value, i.e., 50% of the samples closest to the mean value are within the error bars. The enzymes whose distributions of control coefficients were spread around zero, i.e., whose values did not allow us to predict with certainty the responses of metabolic network, are marked in gray Control coefficients of the cytosolic ATP concentration during glucose-xylose co-utilization. The bars represent the mean values of the control coefficients, and the error bars denote the 1st and the 3rd quartile of the control coefficients with respect to their mean value, i.e., 50% of the samples closest to the mean value are within the error bars

Fig. 5
5Fermentation profiles of the C-10880 strain (blue) and of the engineered HXK2-deficient strain (red): specific xylose consumption rates (a), specific glucose consumption rates (b), specific xylitol production rates (c), and specific ethanol production rates (d)

Fig. 7
7Predicted effects of HXK manipulations on the metabolic fluxes. An increased activity of HXK would result in reactions with: increased flux (red), decreased flux (green), or negligible changes in flux (gray). Higher intensity of red or green indicates larger changes in fluxes. The numerical values shown above reactions denote the mean control coefficients upon changes in HXK, i.e., C * HXK where asterisk denotes any metabolic flux in the network maximal feasible flux value of ATPM for low growth of 18.28 mmol/gDW/h reduced to the one of 12.55 mmol/ gDW/h for high growth.

, p dc xyl dt model to the measurement data, ethanol evaporation is included in the model by first-order dynamics. The measured ethanol concentration is then described by with the evaporation parameter e being a linear function of the fermentation broth volume v

Fig. 8
8Flowchart of the ORACLE framework. The successive application of computational procedures integrates biological information from different sources, thereby refining kinetic models and providing guidance for metabolic engineeringAdditional filesAdditional file 1. Stoichiometry of used reactions and the corresponding mass balances.

Additional file 6 .
6Predicted effects of changes in ATPM on the metabolic fluxes. The increased ATPM would result in reactions with: increased flux (red), decreased flux (green) or negligible changes in flux (grey). Higher intensity of red or green indicates larger changes in fluxes. The numerical values shown above reactions denote the mean control coefficients upon changes in ATPM, i.e. C * ATPM where * denotes any metabolic flux in the network. Additional file 7. Predicted effects of NDR manipulations on the metabolic fluxes. The increased activity of NDR would result in reactions with: increased flux (red), decreased flux (green) or negligible changes in flux (grey). Higher intensity of red or green indicates larger changes in fluxes. The numerical values shown above reactions denote the mean control coefficients upon activity changes in NDR, i.e. C * NDR where * denotes any metabolic flux in the network. Additional file 8. Flux Variability Analysis (FVA): (i) obtained ranges of the fluxes of reactions involved in the cytosolic ATP and NADPH metabolism

. After the pruning, the distributions of both C XTR ATPM and C XTR NDR were shifted toward positive values with their respective means of 0.10 and 0.13. The pruning results suggestedNDR

HXK
PGI
TPI
ZWF

ATPM
ASN

GPD1
ADH1

NDR

ALD
ETOHt
PIt

XTR
XRI
XRII
XDH
XK
XLT

−0.2
0
0.2
0 .4

C

*

XTR

HXK
PGI
PFK
TPI
TDH
PYK

ATPM

ADH1

NDR

PDC
CO2t
O2t
ETOHt
NH4t
SO4t
PIt
COH

XDH

−0.4
−0.2
0
0 .2
0.4
0 .6

C

*

ADH1

HXT
HXK
PGI
PFK
FBA
TPI
TDH
PYK
ATPM

ADH1
PDC
CO2t
O2t
NH4t
SO4t
PIt
COH

−0.4
−0.2
0
0.2
0 .4
0.6

C

*

HXT

HXK
TPI

ATPM
NDH
ASN

ADH1

NDR

ALD
ACE
CO2t
ETOHt
COH

XTR
XRI
XRII
XDH
XK
XLT

−0.5
0
0.5

C

*

XLT

a
b

c
d

#

#

#

Populations of concentration and flux control coefficientsMetabolic engineering strategies identifiedIncreasing knowledge about metabolic networkStoichiometry
Metabolomics
Fluxomics

Sample
elasticities

Sample
metabolites

MCA

Consistency
check
Reject

Analyze results
& design strategies

∆G ' o

No

Yes

Enzyme
kinetics

Thermodynamically
feasible flux profiles

Thermodynamically feasible
metabolite concentrations
and equilibrium displacements

Elasticities consistent with
thermodynamically feasible
metabolite concentrations

AcknowledgementsHeidi Turkia is thanked for the analysis of minor acids in the cultivation samples with CE.l-malate glyoxylate-lyase (CoA-acetylating); MDHc: (S)-malate:NAD+ oxidoreductase; CITc: citrate oxaloacetate-lyase cytosolic; ACOc: citrate hydro-lyase cytosolic; LACm2r: d-lactate transport, mitochondrial; CITt2m: citrate transport, mitochondrial; LDH: (R)-lactate:ferricytochrome-c 2-oxidoreductase; O 2 m: O 2 transport (diffusion); CO 2 m: CO 2 transport (diffusion), mitochondrial; PIm: phosphate transporter, mitochondrial; CO 2 t: CO 2 transport via diffusion; GLYCt: glycerol transport in/out via diffusion reversible; PYRst: pyruvate transport via proton symport; SO 4 t: sulfate transport via proton symport; Pit: phosphate transport via proton symport; O 2 t: O 2 transport via diffusion; NH 4 t: ammonia transport via diffusion; LACt2r: d-lactate transport via proton symport; SUCCt2r: succinate transporter in/out via proton symport; MALt2r: l-malate transport in via proton symport; GND1: 6-phospho-d-glucono-1,5-lactone lactonohydrolase; GND2: 6-phospho-d-gluconate:NADP+ 2-oxidoreductase (decarboxylating); GPD1: glycerol-3-phosphate:NAD+ 2-oxidoreductase; GPD2: glycerol-3-phosphate phosphohydrolase; XL: xylose; XLT: xylitol; XYLL: xylulose; GLC: glucose; G6P: glucose-6-phosphate; F6P: fructose-6-phosphate;Competing interestsThe authors declare that they have no competing interests.Availability of data and materialsAuthors can confirm that all relevant data are included in the article and/or its Additional files 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10.FundingThe research leading to these results has received funding from the European Community's Seventh Framework Programme for research and technological development (FP7) under the project NEMO for bioethanol (Grant Agreement No. 222699).Publisher's NoteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Additional file 9. Estimated parameter values of the kinetic model used to smooth the measurement data together with the raw data (dots) and fitted model outputs (lines) for the cultivation of the VTT C-10880 and HXK2 deletion strains. Correlation parameters for the measurement data and simulated values are also provided.Additional file 10. Rate expressions for used kinetic mechanisms and the expressions for the corresponding metabolite elasticities.• We accept pre-submission inquiries • Our selector tool helps you to find the most relevant journal Submit your next manuscript to BioMed Central and we will help you at every step:
Functional expression of a bacterial xylose isomerase in Saccharomyces cerevisiae. D Brat, E Boles, B Wiedemann, 10.1128/Aem.02522-08Appl Environ Microbiol. 758Brat D, Boles E, Wiedemann B. Functional expression of a bacterial xylose isomerase in Saccharomyces cerevisiae. Appl Environ Microbiol. 2009;75(8):2304-11. doi:10.1128/Aem.02522-08.

Comparison of the xylose reductase-xylitol dehydrogenase and the xylose isomerase pathways for xylose fermentation by recombinant Saccharomyces cerevisiae. Microbial Cell Factories. K Karhumaa, R G Sanchez, B Hahn-Hagerdal, M F Gorwa-Grauslund, 10.1186/1475-2859-6-5Karhumaa K, Sanchez RG, Hahn-Hagerdal B, Gorwa-Grauslund MF. Comparison of the xylose reductase-xylitol dehydrogenase and the xylose isomerase pathways for xylose fermentation by recom- binant Saccharomyces cerevisiae. Microbial Cell Factories. 2007. doi:10.1186/1475-2859-6-5.

Metabolic flux analysis of xylose metabolism in recombinant Saccharomyces cerevisiae using continuous culture. J P Pitkanen, A Aristidou, L Salusjarvi, L Ruohonen, M Penttila, 10.1016/S1096-7176(02)00012-5doi:10.1016/ S1096-7176(02Metab Eng. 51Pitkanen JP, Aristidou A, Salusjarvi L, Ruohonen L, Penttila M. Metabolic flux analysis of xylose metabolism in recombinant Saccharomyces cerevi- siae using continuous culture. Metab Eng. 2003;5(1):16-31. doi:10.1016/ S1096-7176(02)00012-5.

Metabolic engineering for improved fermentation of pentoses by yeasts. T W Jeffries, Y S Jin, 10.1007/S00253-003-1450-0Appl Microbiol Biotechnol. 635Jeffries TW, Jin YS. Metabolic engineering for improved fermentation of pentoses by yeasts. Appl Microbiol Biotechnol. 2004;63(5):495-509. doi:10.1007/S00253-003-1450-0.

Characterization of the xylose-transporting properties of yeast hexose transporters and their influence on xylose utilization. T Hamacher, J Becker, M Gardonyi, B Hahn-Hagerdal, E Boles, Microbiology. 148Hamacher T, Becker J, Gardonyi M, Hahn-Hagerdal B, Boles E. Characteri- zation of the xylose-transporting properties of yeast hexose transporters and their influence on xylose utilization. Microbiology. 2002;148:2783-8.

Control of xylose consumption by xylose transport in recombinant Saccharomyces cerevisiae. M Gardonyi, M Jeppsson, G Liden, M F Gorwa-Grausland, B Hahn-Hagerdal, 10.1002/Bit.10631doi:10.1002/ Bit.10631Biotechnol Bioeng. 827Gardonyi M, Jeppsson M, Liden G, Gorwa-Grausland MF, Hahn-Hagerdal B. Control of xylose consumption by xylose transport in recombinant Sac- charomyces cerevisiae. Biotechnol Bioeng. 2003;82(7):818-24. doi:10.1002/ Bit.10631.

Xylose fermentation by Saccharomyces-cerevisiae. P Kotter, M Ciriacy, Appl Microbiol Biotechnol. 386Kotter P, Ciriacy M. Xylose fermentation by Saccharomyces-cerevisiae. Appl Microbiol Biotechnol. 1993;38(6):776-83.

Evolved hexose transporter enhances xylose uptake and glucose/xylose coutilization in Saccharomyces cerevisiae. A R Apel, M Ouellet, H Szmidt-Middleton, J D Keasling, A Mukhopadhyay, 10.1038/srep19512Sci Rep. Apel AR, Ouellet M, Szmidt-Middleton H, Keasling JD, Mukhopadhyay A. Evolved hexose transporter enhances xylose uptake and glucose/xylose co- utilization in Saccharomyces cerevisiae. Sci Rep. 2016. doi:10.1038/srep19512.

Effect of enhanced xylose reductase activity on xylose consumption and product distribution in xylose-fermenting recombinant Saccharomyces cerevisiae. M Jeppsson, K Traff, B Johansson, B Hahn-Hagerdal, M F Gorwa-Grauslund, FEMS Yeast Res. 32Jeppsson M, Traff K, Johansson B, Hahn-Hagerdal B, Gorwa-Grauslund MF. Effect of enhanced xylose reductase activity on xylose consumption and product distribution in xylose-fermenting recombinant Saccharomyces cerevisiae. FEMS Yeast Res. 2003;3(2):167-75.

Endogenous NADPH-dependent aldose reductase activity influences product formation during xylose consumption in recombinant Saccharomyces cerevisiae. K L Traff-Bjerre, M Jeppsson, B Hahn-Hagerdal, M F Gorwa-Grauslund, Yeast. 212Traff-Bjerre KL, Jeppsson M, Hahn-Hagerdal B, Gorwa-Grauslund MF. Endogenous NADPH-dependent aldose reductase activity influences product formation during xylose consumption in recombinant Saccharo- myces cerevisiae. Yeast. 2004;21(2):141-50.

Conversion of xylose to ethanol by recombinant Saccharomyces cerevisiae: importance of xylulokinase (XKS1) and oxygen availability. M H Toivari, A Aristidou, L Ruohonen, M Penttila, 10.1006/Mben.2000.0191Metab Eng. 33Toivari MH, Aristidou A, Ruohonen L, Penttila M. Conversion of xylose to ethanol by recombinant Saccharomyces cerevisiae: importance of xylulokinase (XKS1) and oxygen availability. Metab Eng. 2001;3(3):236-49. doi:10.1006/Mben.2000.0191.

High expression of XYL2 coding for xylitol dehydrogenase is necessary for efficient xylose fermentation by engineered Saccharomyces cerevisiae. S R Kim, S J Ha, I I Kong, Y S Jin, 10.1016/j.ymben.2012.04.001Metab Eng. 144Kim SR, Ha SJ, Kong II, Jin YS. High expression of XYL2 coding for xylitol dehydrogenase is necessary for efficient xylose fermentation by engineered Saccharomyces cerevisiae. Metab Eng. 2012;14(4):336-43. doi:10.1016/j.ymben.2012.04.001.

Reduced oxidative pentose phosphate pathway flux in recombinant xylose-utilizing Saccharomyces cerevisiae strains improves the ethanol yield from xylose. M Jeppsson, B Johansson, B Hahn-Hagerdal, M F Gorwa-Grauslund, Appl Environ Microbiol. 684Jeppsson M, Johansson B, Hahn-Hagerdal B, Gorwa-Grauslund MF. Reduced oxidative pentose phosphate pathway flux in recombinant xylose-utilizing Saccharomyces cerevisiae strains improves the ethanol yield from xylose. Appl Environ Microbiol. 2002;68(4):1604-9.

Systems metabolic engineering, industrial biotechnology and microbial cell factories. Microbial Cell Factories. S Y Lee, D Mattanovich, A Villaverde, 10.1186/1475-2859-11-15611156Lee SY, Mattanovich D, Villaverde A. Systems metabolic engineering, industrial biotechnology and microbial cell factories. Microbial Cell Facto- ries. 2012;11:156. doi:10.1186/1475-2859-11-156.

Rites of passage: requirements and standards for building kinetic models of metabolic phenotypes. L Miskovic, M Tokic, G Fengos, V Hatzimanikatis, Curr Opin Biotechnol. 36Miskovic L, Tokic M, Fengos G, Hatzimanikatis V. Rites of passage: require- ments and standards for building kinetic models of metabolic pheno- types. Curr Opin Biotechnol. 2015;36:1-8.

The control of flux. H Kacser, J A Burns, Symp Soc Exp Biol. 27Kacser H, Burns JA. The control of flux. Symp Soc Exp Biol. 1973;27:65-104.

The xylose reductase/xylitol dehydrogenase/xylulokinase ratio affects product formation in recombinant xylose-utilising Saccharomyces cerevisiae. A Eliasson, Jhs Hofmeyr, S Pedler, B Hahn-Hagerdal, 10.1016/S0141-0229(01)00386-6Enzyme Microb Technol. 294-5Eliasson A, Hofmeyr JHS, Pedler S, Hahn-Hagerdal B. The xylose reduc- tase/xylitol dehydrogenase/xylulokinase ratio affects product formation in recombinant xylose-utilising Saccharomyces cerevisiae. Enzyme Microb Technol. 2001;29(4-5):288-97. doi:10.1016/S0141-0229(01)00386-6.

Kinetic modelling reveals current limitations in the production of ethanol from xylose by recombinant Saccharomyces cerevisiae. N S Parachin, B Bergdahl, E W Van Niel, M F Gorwa-Grauslund, Metab Eng. 135Parachin NS, Bergdahl B, van Niel EW, Gorwa-Grauslund MF. Kinetic mod- elling reveals current limitations in the production of ethanol from xylose by recombinant Saccharomyces cerevisiae. Metab Eng. 2011;13(5):508-17.

Identification of novel metabolic interactions controlling carbon flux from xylose to ethanol in natural and recombinant yeasts. G Trausinger, C Gruber, S Krahulec, C Magnes, B Nidetzky, M Klimacek, 10.1186/s13068-015-0340-xBiotechnol Biofuels. 8157Trausinger G, Gruber C, Krahulec S, Magnes C, Nidetzky B, Klimacek M. Identification of novel metabolic interactions controlling carbon flux from xylose to ethanol in natural and recombinant yeasts. Biotechnol Biofuels. 2015;8:157. doi:10.1186/s13068-015-0340-x.

Production of biofuels and biochemicals: in need of an ORACLE. L Miskovic, V Hatzimanikatis, Trends Biotechnol. 288Miskovic L, Hatzimanikatis V. Production of biofuels and biochemicals: in need of an ORACLE. Trends Biotechnol. 2010;28(8):391-7.

Towards kinetic modeling of genome-scale metabolic networks without sacrificing stoichiometric, thermodynamic and physiological constraints. A Chakrabarti, L Miskovic, K C Soh, V Hatzimanikatis, 10.1002/biot.201300091Biotechnol J. 89Chakrabarti A, Miskovic L, Soh KC, Hatzimanikatis V. Towards kinetic modeling of genome-scale metabolic networks without sacrificing stoichiometric, thermodynamic and physiological constraints. Biotechnol J. 2013;8(9):1043-57. doi:10.1002/biot.201300091.

Modeling of uncertainties in biochemical reactions. L Miskovic, V Hatzimanikatis, 10.1002/bit.22932doi:10.1002/ bit.22932Biotechnol Bioeng. 1082Miskovic L, Hatzimanikatis V. Modeling of uncertainties in biochemi- cal reactions. Biotechnol Bioeng. 2011;108(2):413-23. doi:10.1002/ bit.22932.

Metabolic control analysis under uncertainty: framework development and case studies. L Wang, I Birol, V Hatzimanikatis, Biophys J. 87Wang L, Birol I, Hatzimanikatis V. Metabolic control analysis under uncertainty: framework development and case studies. Biophys J. 2004;87:3750-63.

Metabolic engineering under uncertainty-II: analysis of yeast metabolism. L Q Wang, V Hatzimanikatis, 10.1016/J.Yinben.2005.11.002doi:10.1016/J. Yinben.2005.11.002Metab Eng. 82Wang LQ, Hatzimanikatis V. Metabolic engineering under uncertainty-II: analysis of yeast metabolism. Metab Eng. 2006;8(2):142-59. doi:10.1016/J. Yinben.2005.11.002.

From network models to network responses: integration of thermodynamic and kinetic properties of yeast genome-scale metabolic networks. K S Soh, L Miskovic, V Hatzimanikatis, FEMS Yeast Res. 12Soh KS, Miskovic L, Hatzimanikatis V. From network models to network responses: integration of thermodynamic and kinetic properties of yeast genome-scale metabolic networks. FEMS Yeast Res. 2012;12:129-43.

Identification of metabolic engineering targets for the enhancement of 1,4-butanediol production in recombinant E. coli using large-scale kinetic models. S Andreozzi, A Chakrabarti, K C Soh, A Burgard, T H Yang, S Van Dien, 10.1016/j.ymben.2016.01.009Metab Eng. 35Andreozzi S, Chakrabarti A, Soh KC, Burgard A, Yang TH, Van Dien S, et al. Identification of metabolic engineering targets for the enhancement of 1,4-butanediol production in recombinant E. coli using large-scale kinetic models. Metab Eng. 2016;35:148-59. doi:10.1016/j.ymben.2016.01.009.

iSCHRUNK-in silico approach to characterization and reduction of uncertainty in the kinetic models of genome-scale metabolic networks. S Andreozzi, L Miskovic, V Hatzimanikatis, 10.1016/j.ymben.2015.10.002Metab Eng. 33Andreozzi S, Miskovic L, Hatzimanikatis V. iSCHRUNK-in silico approach to characterization and reduction of uncertainty in the kinetic models of genome-scale metabolic networks. Metab Eng. 2016;33:158-68. doi:10.1016/j.ymben.2015.10.002.

Metabolic engineering under uncertainty. I: framework development. L Q Wang, V Hatzimanikatis, 10.1016/J.Ymben.2005.11.003doi:10.1016/J. Ymben.2005.11.003Metab Eng. 82Wang LQ, Hatzimanikatis V. Metabolic engineering under uncertainty. I: framework development. Metab Eng. 2006;8(2):133-41. doi:10.1016/J. Ymben.2005.11.003.

The regulation of cellular systems. R Heinrich, S Schuster, SpringerBerlinHeinrich R, Schuster S. The regulation of cellular systems. Berlin: Springer; 1996.

DNA sequences in chromosomes II and VII code for pyruvate carboxylase isoenzymes in Saccharomyces cerevisiae: analysis of pyruvate carboxylase-deficient strains. R Stucka, S Dequin, J M Salmon, C Gancedo, Mol Gen Genet: MGG. 2292Stucka R, Dequin S, Salmon JM, Gancedo C. DNA sequences in chromo- somes II and VII code for pyruvate carboxylase isoenzymes in Saccharo- myces cerevisiae: analysis of pyruvate carboxylase-deficient strains. Mol Gen Genet: MGG. 1991;229(2):307-15.

Hexose and pentose transport in ascomycetous yeasts: an overview. M J Leandro, C Fonseca, P Goncalves, 10.1111/j.1567-1364.2009.00509.xFEMS Yeast Res. 94Leandro MJ, Fonseca C, Goncalves P. Hexose and pentose transport in ascomycetous yeasts: an overview. FEMS Yeast Res. 2009;9(4):511-25. doi:10.1111/j.1567-1364.2009.00509.x.

The hexokinase 2 protein regulates the expression of the GLK1, HXK1 and HXK2 genes of Saccharomyces cerevisiae. A Rodriguez, T De La Cera, P Herrero, F Moreno, Biochem J. 355Rodriguez A, de la Cera T, Herrero P, Moreno F. The hexokinase 2 protein regulates the expression of the GLK1, HXK1 and HXK2 genes of Saccharo- myces cerevisiae. Biochem J. 2001;355:625-31.

Ethanol production from xylose in engineered Saccharomyces cerevisiae strains: current state and perspectives. A Matsushika, H Inoue, T Kodaki, S Sawayama, 10.1007/s00253-009-2101-xdoi:10.1007/ s00253-009-2101-xAppl Microbiol Biotechnol. 841Matsushika 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(1):37-53. doi:10.1007/ s00253-009-2101-x.

The YGR194c (XKS1) gene encodes the xylulokinase from the budding yeast Saccharomyces cerevisiae. J M Rodriguez-Pena, V J Cid, J Arroyo, C Nombela, 10.1111/j.1574-6968.1998.tb12993.xFEMS Microbiol Lett. 1621Rodriguez-Pena JM, Cid VJ, Arroyo J, Nombela C. The YGR194c (XKS1) gene encodes the xylulokinase from the budding yeast Saccharomyces cerevisiae. FEMS Microbiol Lett. 1998;162(1):155-60. doi:10.1111/j.1574-6968.1998.tb12993.x.

Optimal growth and ethanol production from xylose by recombinant Saccharomyces cerevisiae require moderate D-xylulokinase activity. Y S Jin, H Y Ni, J M Laplaza, T W Jeffries, 10.1128/Aem.69.1.495-503.2003Appl Environ Microbiol. 691Jin YS, Ni HY, Laplaza JM, Jeffries TW. Optimal growth and ethanol production from xylose by recombinant Saccharomyces cerevisiae require moderate D-xylulokinase activity. Appl Environ Microbiol. 2003;69(1):495- 503. doi:10.1128/Aem.69.1.495-503.2003.

Efficient bioethanol production from xylose by recombinant Saccharomyces cerevisiae requires high activity of xylose reductase and moderate xylulokinase activity. A Matsushika, S Sawayama, 10.1263/jbb.106.306J Biosci Bioeng. 1063Matsushika A, Sawayama S. Efficient bioethanol production from xylose by recombinant Saccharomyces cerevisiae requires high activity of xylose reductase and moderate xylulokinase activity. J Biosci Bioeng. 2008;106(3):306-9. doi:10.1263/jbb.106.306.

A systems biology approach to study glucose repression in the yeast Saccharomyces cerevisiae. S L Westergaard, A P Oliveira, C Bro, L Olsson, J Nielsen, 10.1002/bit.21135Biotechnol Bioeng. 961Westergaard SL, Oliveira AP, Bro C, Olsson L, Nielsen J. A systems biology approach to study glucose repression in the yeast Saccharomyces cerevi- siae. Biotechnol Bioeng. 2007;96(1):134-45. doi:10.1002/bit.21135.

A systems-level approach for metabolic engineering of yeast cell factories. I K Kim, A Roldao, V Siewers, J Nielsen, 10.1111/j.1567-1364.2011.00779.xFEMS Yeast Res. 122Kim IK, Roldao A, Siewers V, Nielsen J. A systems-level approach for meta- bolic engineering of yeast cell factories. FEMS Yeast Res. 2012;12(2):228- 48. doi:10.1111/j.1567-1364.2011.00779.x.

Yeast carbon catabolite repression. J M Gancedo, Microbiol Mol Biol Rev. 622Gancedo JM. Yeast carbon catabolite repression. Microbiol Mol Biol Rev. 1998;62(2):334-61.

Carbon source-dependent phosphorylation of hexokinase PII and its role in the glucose-signaling response in yeast. F Randez-Gil, P Sanz, K D Entian, J A Prieto, Mol Cell Biol. 185Randez-Gil F, Sanz P, Entian KD, Prieto JA. Carbon source-dependent phosphorylation of hexokinase PII and its role in the glucose-signaling response in yeast. Mol Cell Biol. 1998;18(5):2940-8.

Differential glucose repression in common yeast strains in response to HXK2 deletion. A Kummel, J C Ewald, S M Fendt, S J Jol, P Picotti, R Aebersold, 10.1111/j.1567-1364.2010.00609.xFEMS Yeast Res. 103Kummel A, Ewald JC, Fendt SM, Jol SJ, Picotti P, Aebersold R, et al. Differential glucose repression in common yeast strains in response to HXK2 deletion. FEMS Yeast Res. 2010;10(3):322-32. doi:10.1111/j.1567-1364.2010.00609.x.

Engineering of carbon catabolite repression in recombinant xylose fermenting Saccharomyces cerevisiae. C Roca, M B Haack, L Olsson, 10.1007/s00253-003-1408-2Appl Microbiol Biotechnol. 635Roca C, Haack MB, Olsson L. Engineering of carbon catabolite repression in recombinant xylose fermenting Saccharomyces cerevisiae. Appl Micro- biol Biotechnol. 2004;63(5):578-83. doi:10.1007/s00253-003-1408-2.

A new efficient gene disruption cassette for repeated use in budding yeast. U Guldener, S Heck, T Fiedler, J Beinhauer, J H Hegemann, 10.1093/Nar/24.13.2519Nucleic Acids Res. 2413Guldener U, Heck S, Fiedler T, Beinhauer J, Hegemann JH. A new efficient gene disruption cassette for repeated use in budding yeast. Nucleic Acids Res. 1996;24(13):2519-24. doi:10.1093/Nar/24.13.2519.

New heterologous modules for classical or pcr-based gene disruptions in Saccharomyces-cerevisiae. A Wach, A Brachat, R Pohlmann, P Philippsen, 10.1002/Yea.320101310Yeast. 1013Wach A, Brachat A, Pohlmann R, Philippsen P. New heterologous modules for classical or pcr-based gene disruptions in Saccharomyces-cerevisiae. Yeast. 1994;10(13):1793-808. doi:10.1002/Yea.320101310.

New yeast-Escherichia-coli shuttle vectors constructed with invitro mutagenized yeast genes lacking 6-base pair restriction sites. R D Gietz, A Sugino, 10.1016/0378-1119(88)90185-0Gene. 742Gietz RD, Sugino A. New yeast-Escherichia-coli shuttle vectors constructed with invitro mutagenized yeast genes lacking 6-base pair restriction sites. Gene. 1988;74(2):527-34. doi:10.1016/0378-1119(88)90185-0.

High-efficiency transformation of plasmid DNA into yeast. R A Woods, R D Gietz, 10.1385/1-59259-210-4:085Methods Mol Biol. 177Woods RA, Gietz RD. High-efficiency transformation of plasmid DNA into yeast. Methods Mol Biol. 2001;177:85-97. doi:10.1385/1-59259-210-4:085.

Effect of benzoic acid on metabolic fluxes in yeasts: a continuous-culture study on the regulation of respiration and alcoholic fermentation. C Verduyn, E Postma, W A Scheffers, J P Van Dijken, 10.1002/yea.320080703Yeast. 87Verduyn C, Postma E, Scheffers WA, Van Dijken JP. Effect of benzoic acid on metabolic fluxes in yeasts: a continuous-culture study on the regula- tion of respiration and alcoholic fermentation. Yeast. 1992;8(7):501-17. doi:10.1002/yea.320080703.

Capillary electrophoresis for the monitoring of carboxylic acid production by Gluconobacter oxydans. H Turkia, H Siren, J P Pitkanen, M Wiebe, M Penttila, 10.1016/j.chroma.2009.12.075doi:10.1016/j. chroma.2009.12.075J Chromatogr A. 12179Turkia H, Siren H, Pitkanen JP, Wiebe M, Penttila M. Capillary electro- phoresis for the monitoring of carboxylic acid production by Glucono- bacter oxydans. J Chromatogr A. 2010;1217(9):1537-42. doi:10.1016/j. chroma.2009.12.075.

Characterisation of glucose transport in Saccharomyces cerevisiae with plasma membrane vesicles (countertransport) and intact cells (initial uptake) with single Hxt1, Hxt2, Hxt3, Hxt4, Hxt6, Hxt7 or Gal2 transporters. A Maier, B Volker, E Boles, G F Fuhrmann, FEMS Yeast Res. 24Maier A, Volker B, Boles E, Fuhrmann GF. Characterisation of glucose trans- port in Saccharomyces cerevisiae with plasma membrane vesicles (coun- tertransport) and intact cells (initial uptake) with single Hxt1, Hxt2, Hxt3, Hxt4, Hxt6, Hxt7 or Gal2 transporters. FEMS Yeast Res. 2002;2(4):539-50.

Characterization of the effectiveness of hexose transporters for transporting xylose during glucose and xylose co-fermentation by a recombinant Saccharomyces yeast. M Sedlak, N W Ho, 10.1002/yea.1060Yeast. 218Sedlak M, Ho NW. Characterization of the effectiveness of hexose trans- porters for transporting xylose during glucose and xylose co-fermen- tation by a recombinant Saccharomyces yeast. Yeast. 2004;21(8):671-84. doi:10.1002/yea.1060.

Enzyme kinetics. I H Segel, WileyNew YorkSegel IH. Enzyme kinetics. New York: Wiley; 1975.

The reversible Hill equation: how to incorporate cooperative enzymes into metabolic models. J Hofmeyr, A Cornish-Bowden, Comp Appl Biosci. 13Hofmeyr J, Cornish-Bowden A. The reversible Hill equation: how to incor- porate cooperative enzymes into metabolic models. Comp Appl Biosci. 1997;13:377-85.

Bringing metabolic networks to life: convenience rate law and thermodynamic constraints. W Liebermeister, E Klipp, 10.1186/1742-4682-3-41Theor Biol Med Model. 341Liebermeister W, Klipp E. Bringing metabolic networks to life: conveni- ence rate law and thermodynamic constraints. Theor Biol Med Model. 2006;3:41. doi:10.1186/1742-4682-3-41.

Relationship between thermodynamic driving force and one-way fluxes in reversible processes. D Beard, H Qian, 10.1371/journal.pone.0000144PLoS ONE. 21Beard D, Qian H. Relationship between thermodynamic driving force and one-way fluxes in reversible processes. PLoS ONE. 2007;2(1):e144. doi:10.1371/journal.pone.0000144.

Toward the automated generation of genome-scale metabolic networks in the SEED. M Dejongh, K Formsma, P Boillot, J Gould, M Rycenga, A Best, 10.1186/1471-2105-8-139BMC Bioinform. 8139DeJongh M, Formsma K, Boillot P, Gould J, Rycenga M, Best A. Toward the automated generation of genome-scale metabolic networks in the SEED. BMC Bioinform. 2007;8:139. doi:10.1186/1471-2105-8-139.

A protocol for generating a high-quality genomescale metabolic reconstruction. I Thiele, B O Palsson, 10.1038/nprot.2009.203Nat Protoc. 51Thiele I, Palsson BO. A protocol for generating a high-quality genome- scale metabolic reconstruction. Nat Protoc. 2010;5(1):93-121. doi:10.1038/nprot.2009.203.

Global reconstruction of the human metabolic network based on genomic and bibliomic data. N C Duarte, S Becker, N Jamshidi, I Thiele, M Mo, T D Vo, 10.1073/pnas.0610772104doi:10.1073/ pnas.0610772104Proc Natl Acad Sci USA. 1046Duarte NC, Becker S, Jamshidi N, Thiele I, Mo M, Vo TD, et al. Global reconstruction of the human metabolic network based on genomic and bibliomic data. Proc Natl Acad Sci USA. 2007;104(6):1777-82. doi:10.1073/ pnas.0610772104.

GrowMatch: an automated method for reconciling in silico/in vivo growth predictions. V S Kumar, C D Maranas, 10.1371/journal.pcbi.1000308PLoS Comput Biol. 531000308Kumar VS, Maranas CD. GrowMatch: an automated method for reconciling in silico/in vivo growth predictions. PLoS Comput Biol. 2009;5(3):e1000308. doi:10.1371/journal.pcbi.1000308.

Reconstruction of biochemical networks in microorganisms. A M Feist, M J Herrgard, I Thiele, J L Reed, B O Palsson, 10.1038/nrmicro1949Nat Rev Microbiol. 72Feist AM, Herrgard MJ, Thiele I, Reed JL, Palsson BO. Reconstruction of bio- chemical networks in microorganisms. Nat Rev Microbiol. 2009;7(2):129- 43. doi:10.1038/nrmicro1949.

iBsu1103: a new genomescale metabolic model of Bacillus subtilis based on SEED annotations. C S Henry, J F Zinner, M P Cohoon, R L Stevens, Henry CS, Zinner JF, Cohoon MP, Stevens RL. iBsu1103: a new genome- scale metabolic model of Bacillus subtilis based on SEED annotations.

. 10.1186/gb-2009-10-6-r69Genome Biol. 10669Genome Biol. 2009;10(6):R69. doi:10.1186/gb-2009-10-6-r69.

A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. A M Feist, C S Henry, J L Reed, M Krummenacker, A R Joyce, P D Karp, 10.1038/Msb4100155Mol Syst Biol. 3121Feist AM, Henry CS, Reed JL, Krummenacker M, Joyce AR, Karp PD, et al. A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Mol Syst Biol. 2007;3:121. doi:10.1038/Msb4100155.

Novel biological insights through metabolomics and C-13-flux analysis. N Zamboni, U Sauer, 10.1016/j.mib.2009.08.003Curr Opin Microbiol. 125Zamboni N, Sauer U. Novel biological insights through metabo- lomics and C-13-flux analysis. Curr Opin Microbiol. 2009;12(5):553-8. doi:10.1016/j.mib.2009.08.003.

Industrial systems biology. J M Otero, J Nielsen, 10.1002/bit.22592Biotechnol Bioeng. 1053Otero JM, Nielsen J. Industrial systems biology. Biotechnol Bioeng. 2010;105(3):439-60. doi:10.1002/bit.22592.

C-13 metabolic flux analysis. W Wiechert, 10.1006/mben.2001.0187Metab Eng. 33Wiechert W. C-13 metabolic flux analysis. Metab Eng. 2001;3(3):195-206. doi:10.1006/mben.2001.0187.

Group contribution method for thermodynamic analysis of complex metabolic networks. M Jankowski, C Henry, L Broadbelt, V Hatzimanikatis, Biophys J. 95Jankowski M, Henry C, Broadbelt L, Hatzimanikatis V. Group contribution method for thermodynamic analysis of complex metabolic networks. Biophys J. 2008;95:1487-99.

Thermodynamics of enzyme-catalyzed reactions-a database for quantitative biochemistry. R N Goldberg, Y B Tewari, T N Bhat, 10.1093/bioinformatics/bth314Bioinformatics. 2016Goldberg RN, Tewari YB, Bhat TN. Thermodynamics of enzyme-catalyzed reactions-a database for quantitative biochemistry. Bioinformatics. 2004;20(16):2874-7. doi:10.1093/bioinformatics/bth314.

Thermodynamics of enzyme-catalyzed reactions: part 7-2007 update. R N Goldberg, Y B Tewari, T N Bhat, 10.1063/1.2789450J Phys Chem Ref Data. 364Goldberg RN, Tewari YB, Bhat TN. Thermodynamics of enzyme-catalyzed reactions: part 7-2007 update. J Phys Chem Ref Data. 2007;36(4):1347- 97. doi:10.1063/1.2789450.

Estimation of standard Gibbs energy changes of biotransformations. M L Mavrovouniotis, J Biol Chem. 26622Mavrovouniotis ML. Estimation of standard Gibbs energy changes of biotransformations. J Biol Chem. 1991;266(22):14440-5.

Thermodynamics-based metabolic flux analysis. C Henry, L Broadbelt, V Hatzimanikatis, Biophys J. 92Henry C, Broadbelt L, Hatzimanikatis V. Thermodynamics-based meta- bolic flux analysis. Biophys J. 2007;92:1792-805.

Constraining the flux space using thermodynamics and integration of metabolomics data. K S Soh, V Hatzimanikatis, Methods Mol Biol. 1191Soh KS, Hatzimanikatis V. Constraining the flux space using thermo- dynamics and integration of metabolomics data. Methods Mol Biol. 2014;1191:49-63.

Random number generation and Monte Carlo methods. J Gentle, SpringerNew York, NY2nd edGentle J. Random number generation and Monte Carlo methods. 2nd ed. New York, NY: Springer; 2003.

Various techniques used in connection with random digits. Von Neumann, J , John von Neumann. Taub AH5Pergamon PressVon Neumann J. Various techniques used in connection with random digits. In: Taub AH, editor. The collected works of John von Neumann. Vol. 5. Oxford: Pergamon Press; 1963. P. 768-770.

Simulation and the Monte Carlo method. R Rubinstein, WileyNew YorkRubinstein R. Simulation and the Monte Carlo method. New York: Wiley; 1981.

Generating random vectors uniformly distributed inside and on the surface of different regions. R Rubinstein, Europ J Oper Res. 10Rubinstein R. Generating random vectors uniformly distributed inside and on the surface of different regions. Europ J Oper Res. 1982;10:205-9.

Growthlimiting intracellular metabolites in yeast growing under diverse nutrient limitations. V M Boer, C A Crutchfield, P H Bradley, D Botstein, J D Rabinowitz, 10.1091/mbc.E09-07-0597Mol Biol Cell. 211Boer VM, Crutchfield CA, Bradley PH, Botstein D, Rabinowitz JD. Growth- limiting intracellular metabolites in yeast growing under diverse nutrient limitations. Mol Biol Cell. 2010;21(1):198-211. doi:10.1091/mbc. E09-07-0597.

Absolute metabolite concentrations and implied enzyme active site occupancy in Escherichia coli. B D Bennett, E H Kimball, M Gao, R Osterhout, S J Van Dien, J D Rabinowitz, 10.1038/nchembio.186Nat Chem Biol. 58Bennett BD, Kimball EH, Gao M, Osterhout R, Van Dien SJ, Rabinowitz JD. Absolute metabolite concentrations and implied enzyme active site occupancy in Escherichia coli. Nat Chem Biol. 2009;5(8):593-9. doi:10.1038/nchembio.186.

Cellular metabolomics of Escherichia coli. J D Rabinowitz, 10.1586/147894188.8.131.52Expert Rev Proteom. 42Rabinowitz JD. Cellular metabolomics of Escherichia coli. Expert Rev Proteom. 2007;4(2):187-98. doi:10.1586/14789422.214.171.124.

Practical nonparametric statistics. W J Conover, WileyNew YorkConover WJ. Practical nonparametric statistics. New York: Wiley; 1999.

Biological data mining. J Chen, S Lonardi, Chapman & Hall/ CRCBoca RatonChen J, Lonardi S. Biological data mining. Boca Raton: Chapman & Hall/ CRC; 2010.

Density estimation for statistics and data analysis. B W Silverman, CRC PressBoca RatonSilverman BW. Density estimation for statistics and data analysis. Boca Raton: CRC Press; 1986.