Genome-scale metabolic models reveal determinants of phenotypic differences in non-Saccharomyces yeasts.

BACKGROUND: Use of alternative non-Saccharomyces yeasts in wine and beer brewing has gained more attention the recent years. This is both due to the desire to obtain a wider variety of flavours in the product and to reduce the final alcohol content. Given the metabolic differences between the yeast species, we wanted to account for some of the differences by using in silico models. RESULTS: We created and studied genome-scale metabolic models of five different non-Saccharomyces species using an automated processes. These were: Metschnikowia pulcherrima, Lachancea thermotolerans, Hanseniaspora osmophila, Torulaspora delbrueckii and Kluyveromyces lactis. Using the models, we predicted that M. pulcherrima, when compared to the other species, conducts more respiration and thus produces less fermentation products, a finding which agrees with experimental data. Complex I of the electron transport chain was to be present in M. pulcherrima, but absent in the others. The predicted importance of Complex I was diminished when we incorporated constraints on the amount of enzymatic protein, as this shifts the metabolism towards fermentation. CONCLUSIONS: Our results suggest that Complex I in the electron transport chain is a key differentiator between Metschnikowia pulcherrima and the other yeasts considered. Yet, more annotations and experimental data have the potential to improve model quality in order to increase fidelity and confidence in these results. Further experiments should be conducted to confirm the in vivo effect of Complex I in M. pulcherrima and its respiratory metabolism.

Inducible Synthetic Growth Regulation Using the ClpXP Proteasome Enhances cis,cis-Muconic Acid and Glycolic Acid Yields in Saccharomyces cerevisiae.

Engineered microbial cells can produce sustainable chemistry, but the production competes for resources with growth. Inducible synthetic control over the resource use would enable fast accumulation of sufficient biomass and then divert the resources to production. We developed inducible synthetic resource-use control overSaccharomyces cerevisiae by expressing a bacterial ClpXP proteasome from an inducible promoter. By individually targeting growth-essential metabolic enzymes Aro1, Hom3, and Acc1 to the ClpXP proteasome, cell growth could be efficiently repressed during cultivation. The ClpXP proteasome was specific to the target proteins, and there was no reduction in the targets when ClpXP was not induced. The inducible growth repression improved product yields from glucose (cis,cis-muconic acid) and per biomass (cis,cis-muconic acid and glycolic acid). The inducible ClpXP proteasome tackles uncertainties in strain optimization by enabling model-guided repression of competing, growth-essential, and metabolic enzymes. Most importantly, it allows improving production without compromising biomass accumulation when uninduced; therefore, it is expected to mitigate strain stability and low productivity challenges.

Predictive evolution of metabolic phenotypes using model-designed environments.

Adaptive evolution under controlled laboratory conditions has been highly effective in selecting organisms with beneficial phenotypes such as stress tolerance. The evolution route is particularly attractive when the organisms are either difficult to engineer or the genetic basis of the phenotype is complex. However, many desired traits, like metabolite secretion, have been inaccessible to adaptive selection due to their trade-off with cell growth. Here, we utilize genome-scale metabolic models to design nutrient environments for selecting lineages with enhanced metabolite secretion. To overcome the growth-secretion trade-off, we identify environments wherein growth becomes correlated with a secondary trait termed tacking trait. The latter is selected to be coupled with the desired trait in the application environment where the trait manifestation is required. Thus, adaptive evolution in the model-designed selection environment and subsequent return to the application environment is predicted to enhance the desired trait. We experimentally validate this strategy by evolving Saccharomyces cerevisiae for increased secretion of aroma compounds, and confirm the predicted flux-rerouting using genomic, transcriptomic, and proteomic analyses. Overall, model-designed selection environments open new opportunities for predictive evolution.

Production of biopolymer precursors beta-alanine and L-lactic acid from CO2 with metabolically versatile Rhodococcus opacus DSM 43205.

Hydrogen oxidizing autotrophic bacteria are promising hosts for conversion of CO2 into chemicals. In this work, we engineered the metabolically versatile lithoautotrophic bacterium R. opacus strain DSM 43205 for synthesis of polymer precursors. Aspartate decarboxylase (panD) or lactate dehydrogenase (ldh) were expressed for beta-alanine or L-lactic acid production, respectively. The heterotrophic cultivations on glucose produced 25 mg L-1 beta-alanine and 742 mg L-1 L-lactic acid, while autotrophic cultivations with CO2, H2, and O2 resulted in the production of 1.8 mg L-1 beta-alanine and 146 mg L-1 L-lactic acid. Beta-alanine was also produced at 345 µg L-1 from CO2 in electrobioreactors, where H2 and O2 were provided by water electrolysis. This work demonstrates that R. opacus DSM 43205 can be engineered to produce chemicals from CO2 and provides a base for its further metabolic engineering.

Operation of a TCA cycle subnetwork in the mammalian nucleus.

Nucleic acid and histone modifications critically depend on the tricarboxylic acid (TCA) cycle for substrates and cofactors. Although a few TCA cycle enzymes have been reported in the nucleus, the corresponding pathways are considered to operate in mitochondria. Here, we show that a part of the TCA cycle is operational also in the nucleus. Using 13C-tracer analysis, we identified activity of glutamine-to-fumarate, citrate-to-succinate, and glutamine-to-aspartate routes in the nuclei of HeLa cells. Proximity labeling mass spectrometry revealed a spatial vicinity of the involved enzymes with core nuclear proteins. We further show nuclear localization of aconitase 2 and 2-oxoglutarate dehydrogenase in mouse embryonic stem cells. Nuclear localization of the latter enzyme, which produces succinyl-CoA, changed from pluripotency to a differentiated state with accompanying changes in the nuclear protein succinylation. Together, our results demonstrate operation of an extended metabolic pathway in the nucleus, warranting a revision of the canonical view on metabolic compartmentalization.

Strain design optimization using reinforcement learning.

Engineered microbial cells present a sustainable alternative to fossil-based synthesis of chemicals and fuels. Cellular synthesis routes are readily assembled and introduced into microbial strains using state-of-the-art synthetic biology tools. However, the optimization of the strains required to reach industrially feasible production levels is far less efficient. It typically relies on trial-and-error leading into high uncertainty in total duration and cost. New techniques that can cope with the complexity and limited mechanistic knowledge of the cellular regulation are called for guiding the strain optimization. In this paper, we put forward a multi-agent reinforcement learning (MARL) approach that learns from experiments to tune the metabolic enzyme levels so that the production is improved. Our method is model-free and does not assume prior knowledge of the microbe's metabolic network or its regulation. The multi-agent approach is well-suited to make use of parallel experiments such as multi-well plates commonly used for screening microbial strains. We demonstrate the method's capabilities using the genome-scale kinetic model of Escherichia coli, k-ecoli457, as a surrogate for an in vivo cell behaviour in cultivation experiments. We investigate the method's performance relevant for practical applicability in strain engineering i.e. the speed of convergence towards the optimum response, noise tolerance, and the statistical stability of the solutions found. We further evaluate the proposed MARL approach in improving L-tryptophan production by yeast Saccharomyces cerevisiae, using publicly available experimental data on the performance of a combinatorial strain library. Overall, our results show that multi-agent reinforcement learning is a promising approach for guiding the strain optimization beyond mechanistic knowledge, with the goal of faster and more reliably obtaining industrially attractive production levels.

PHB production from cellobiose with Saccharomyces cerevisiae.

Replacement of petrochemical-based materials with microbially produced biodegradable alternatives calls for industrially attractive fermentation processes. Lignocellulosic materials offer non-edible alternatives for cultivated sugars, but require often use of expensive sugar releasing enzymes, such as ß-glucosidases. These cellulose treatment costs could be reduced if microbial production hosts could use short cellodextrins such as cellobiose directly as their substrates. In this study, we demonstrate production of poly(hydroxybutyrate) (PHB) in yeast Saccharomyces cerevisiae using cellobiose as a sole carbon source. Yeast strains expressing PHB pathway genes from Cupriavidus necator and cellodextrin transporter gene CDT-1 from Neurospora crassa were complemented either with ß-glucosidase gene GH1-1 from N. crassa or with cellobiose phosphorylase gene cbp from Ruminococcus flavefaciens. These cellobiose utilization routes either with Gh1-1 or Cbp enzymes differ in energetics and dynamics. However, both routes enabled higher PHB production per consumed sugar and higher PHB accumulation % of cell dry weight (CDW) than use of glucose as a carbon source. As expected, the strains with Gh1-1 consumed cellobiose faster than the strains with Cbp, both in flask and bioreactor batch cultures. In shake flasks, higher final PHB accumulation % of CDW was reached with Cbp route (10.0 ± 0.3%) than with Gh1-1 route (8.1 ± 0.2%). However, a higher PHB accumulation was achieved in better aerated and pH-controlled bioreactors, in comparison to shake flasks, and the relative performance of strains switched. In bioreactors, notable PHB accumulation levels per CDW of 13.4 ± 0.9% and 18.5 ± 3.9% were achieved with Cbp and Gh1-1 routes, respectively. The average molecular weights of accumulated PHB were similar using both routes; approximately 500 kDa and 450 kDa for strains expressing either cbp or GH1-1 genes, respectively. The formation of PHB with high molecular weights, combined with efficient cellobiose conversion, demonstrates a highly potential solution for improving attractiveness of sustainable polymer production using microbial cells.

Adaptive laboratory evolution of microbial co-cultures for improved metabolite secretion.

Adaptive laboratory evolution has proven highly effective for obtaining microorganisms with enhanced capabilities. Yet, this method is inherently restricted to the traits that are positively linked to cell fitness, such as nutrient utilization. Here, we introduce coevolution of obligatory mutualistic communities for improving secretion of fitness-costly metabolites through natural selection. In this strategy, metabolic cross-feeding connects secretion of the target metabolite, despite its cost to the secretor, to the survival and proliferation of the entire community. We thus co-evolved wild-type lactic acid bacteria and engineered auxotrophic Saccharomyces cerevisiae in a synthetic growth medium leading to bacterial isolates with enhanced secretion of two B-group vitamins, viz., riboflavin and folate. The increased production was specific to the targeted vitamin, and evident also in milk, a more complex nutrient environment that naturally contains vitamins. Genomic, proteomic and metabolomic analyses of the evolved lactic acid bacteria, in combination with flux balance analysis, showed altered metabolic regulation towards increased supply of the vitamin precursors. Together, our findings demonstrate how microbial metabolism adapts to mutualistic lifestyle through enhanced metabolite exchange.

Model-guided development of an evolutionarily stable yeast chassis.

First-principle metabolic modelling holds potential for designing microbial chassis that are resilient against phenotype reversal due to adaptive mutations. Yet, the theory of model-based chassis design has rarely been put to rigorous experimental test. Here, we report the development of Saccharomyces cerevisiae chassis strains for dicarboxylic acid production using genome-scale metabolic modelling. The chassis strains, albeit geared for higher flux towards succinate, fumarate and malate, do not appreciably secrete these metabolites. As predicted by the model, introducing product-specific TCA cycle disruptions resulted in the secretion of the corresponding acid. Adaptive laboratory evolution further improved production of succinate and fumarate, demonstrating the evolutionary robustness of the engineered cells. In the case of malate, multi-omics analysis revealed a flux bypass at peroxisomal malate dehydrogenase that was missing in the yeast metabolic model. In all three cases, flux balance analysis integrating transcriptomics, proteomics and metabolomics data confirmed the flux re-routing predicted by the model. Taken together, our modelling and experimental results have implications for the computer-aided design of microbial cell factories.

Yeast Genome-Scale Metabolic Models for Simulating Genotype-Phenotype Relations.

Understanding genotype-phenotype dependency is a universal aim for all life sciences. While the complete genotype-phenotype relations remain challenging to resolve, metabolic phenotypes are moving within the reach through genome-scale metabolic model simulations. Genome-scale metabolic models are available for commonly investigated yeasts, such as model eukaryote and domesticated fermentation species Saccharomyces cerevisiae, and automatic reconstruction methods facilitate obtaining models for any sequenced species. The models allow for investigating genotype-phenotype relations through simulations simultaneously considering the effects of nutrient availability, and redox and energy homeostasis in cells. Genome-scale models also offer frameworks for omics data integration to help to uncover how the translation of genotypes to the apparent phenotypes is regulated at different levels. In this chapter, we provide an overview of the yeast genome-scale metabolic models and the simulation approaches for using these models to interrogate genotype-phenotype relations. We review the methodological approaches according to the underlying biological reasoning in order to inspire formulating novel questions and applications that the genome-scale metabolic models could contribute to. Finally, we discuss current challenges and opportunities in the genome-scale metabolic model simulations.

Enhanced Triacylglycerol Production With Genetically Modified Trichosporon oleaginosus.

Mitochondrial pyruvate dehydrogenase (PDH) is important in the production of lipids in oleaginous yeast, but other yeast may bypass the mitochondria (PDH bypass), converting pyruvate in the cytosol to acetaldehyde, then acetate and acetyl CoA which is further converted to lipids. Using a metabolic model based on the oleaginous yeast Yarrowia lipolytica, we found that introduction of this bypass to an oleaginous yeast should result in enhanced yield of triacylglycerol (TAG) on substrate. Trichosporon oleaginosus (formerly Cryptococcus curvatus) is an oleaginous yeast which can produce TAGs from both glucose and xylose. Based on the sequenced genome, it lacks at least one of the enzymes needed to complete the PDH bypass, acetaldehyde dehydrogenase (ALD), and may also be deficient in pyruvate decarboxylase and acetyl-CoA synthetase under production conditions. We introduced these genes to T. oleaginosus in various combinations and demonstrated that the yield of TAG on both glucose and xylose was improved, particularly at high C/N ratio. Expression of a phospholipid:diacyltransferase encoding gene in conjunction with the PDH bypass further enhanced lipid production. The yield of TAG on xylose (0.27 g/g) in the engineered strain approached the theoretical maximum yield of 0.289 g/g. Interestingly, TAG production was also enhanced compared to the control in some strains which were given only part of the bypass pathway, suggesting that these genes may contribute to alternative routes to cytoplasmic acetyl CoA. The metabolic model indicated that the improved yield of TAG on substrate in the PDH bypass was dependent on the production of NADPH by ALD. NADPH for lipid synthesis is otherwise primarily supplied by the pentose phosphate pathway (PPP). This would contribute to the greater improvement of TAG production from xylose compared to that observed from glucose when the PDH bypass was introduced, since xylose enters metabolism through the non-oxidative part of the PPP. Yield of TAG from xylose in the engineered strains (0.21-0.27 g/g) was comparable to that obtained from glucose and the highest so far reported for lipid or TAG production from xylose.

Nutritional preferences of human gut bacteria reveal their metabolic idiosyncrasies.

Bacterial metabolism plays a fundamental role in gut microbiota ecology and host-microbiome interactions. Yet the metabolic capabilities of most gut bacteria have remained unknown. Here we report growth characteristics of 96 phylogenetically diverse gut bacterial strains across 4 rich and 15 defined media. The vast majority of strains (76) grow in at least one defined medium, enabling accurate assessment of their biosynthetic capabilities. These do not necessarily match phylogenetic similarity, thus indicating a complex evolution of nutritional preferences. We identify mucin utilizers and species inhibited by amino acids and short-chain fatty acids. Our analysis also uncovers media for in vitro studies wherein growth capacity correlates well with in vivo abundance. Further value of the underlying resource is demonstrated by correcting pathway gaps in available genome-scale metabolic models of gut microorganisms. Together, the media resource and the extracted knowledge on growth abilities widen experimental and computational access to the gut microbiota.

Laboratory evolution reveals regulatory and metabolic trade-offs of glycerol utilization in Saccharomyces cerevisiae.

Most microbial species, including model eukaryote Saccharomyces cerevisiae, possess genetic capability to utilize many alternative nutrient sources. Yet, it remains an open question whether these manifest into assimilatory phenotypes. Despite possessing all necessary pathways, S. cerevisiae grows poorly or not at all when glycerol is the sole carbon source. Here we discover, through multiple evolved lineages, genetic determinants underlying glycerol catabolism and the associated fitness trade-offs. Most evolved lineages adapted through mutations in the HOG pathway, but showed hampered osmotolerance. In the other lineages, we find that only three mutations cause the improved phenotype. One of these contributes counter-intuitively by decoupling the TCA cycle from oxidative phosphorylation, and thereby hampers ethanol utilization. Transcriptomics, proteomics and metabolomics analysis of the re-engineered strains affirmed the causality of the three mutations at molecular level. Introduction of these mutations resulted in improved glycerol utilization also in industrial strains. Our findings not only have a direct relevance for improving glycerol-based bioprocesses, but also illustrate how a metabolic pathway can remain unexploited due to fitness trade-offs in other, ecologically important, traits.

Low Phenotypic Penetrance and Technological Impact of Yeast [GAR +] Prion-Like Elements on Winemaking.

[GAR +] prion-like elements partially relieve carbon catabolite repression in Saccharomyces cerevisiae. They have been hypothesized to contribute to wine yeast survival and alcohol level reduction, as well as communication with bacteria and stuck fermentation. In this work, we selected [GAR +] derivatives from several genetic backgrounds. They were characterized for phenotypic penetrance, heritability and confirmed as prion-like through curing by desiccation. In terms of fermentation kinetics, the impact of the prion on anaerobic wine fermentation (natural grape juice) was either neutral or negative, depending on the genetic background. Likewise, residual sugars were higher or similar for [GAR +] as compared to the cognate [gar -] strains. The prions had little or no impact on glycerol and ethanol yields; while acetic acid yields experienced the highest variations between [GAR +] and [gar -] strains. Strains analyzed under aerobic conditions followed the same pattern, with either little or no impact on fermentation kinetics, ethanol or glycerol yield; and a clearer influence on volatile acidity. Although no clear winemaking advantages were found for [GAR +] strains in this work, they might eventually show interest for some combinations of genetic background or winemaking conditions, e.g., for reducing acetic acid yield under aerated fermentation.

13 C metabolic flux profiling of Pichia pastoris grown in aerobic batch cultures on glucose revealed high relative anabolic use of TCA cycle and limited incorporation of provided precursors of branched-chain amino acids.

Carbon metabolism of Crabtree-negative yeast Pichia pastoris was profiled using 13 C nuclear magnetic resonance (NMR) to delineate regulation during exponential growth and to study the import of two precursors for branched-chain amino acid biosynthesis, α-ketoisovalerate and α-ketobutyrate. Cells were grown in aerobic batch cultures containing (a) only glucose, (b) glucose along with the precursors, or (c) glucose and Val. The study provided the following new insights. First, 13 C flux ratio analyses of central metabolism reveal an unexpectedly high anaplerotic supply of the tricarboxylic acid cycle for a Crabtree-negative yeast, and show that a substantial fraction of glucose catabolism proceeds through the pentose phosphate pathway. A comparison with previous flux ratio analyses for batch cultures of Crabtree-negative Pichia stipitis and Crabtree-positive Saccharomyces cerevisiae indicate that the overall regulation of central carbon metabolism in P. pastoris is intermediate in between P. stipitis and S. cerevisiae. Second, excess α-ketoisovalerate in the medium is not transported into the cytoplasm indicating that P. pastoris lacks a suitable transporter. In contrast, excess Val is efficiently taken up and largely fulfills demands for both Val and Leu for protein synthesis. Third, excess α-ketobutyrate is transported into the mitochondria for Ile biosynthesis. However, the import does not efficiently inhibit the synthesis of α-ketobutyrate from pyruvate indicating that P. pastoris has not been optimized evolutionarily to take full advantage of this carbon source. These findings have direct implications for preparing uniformly 2 H,13 C,15 N-labeled proteins containing protonated Ile, Val, and Leu methyl groups in P. pastoris for NMR-based structural biology. ENZYMES: Acetohydroxy acid isomeroreductase (EC 1.1.1.86), branched-chain amino acid aminotransferase (BCAT, EC 2.6.1.42), fumarase (EC 4.2.1.2), malic enzyme (EC 1.1.1.39/1.1.1.40), phosphoenolpyruvate carboxykinase (EC 4.1.1.49), pyruvate carboxylase (EC 6.4.1.1), pyruvate kinase (EC 2.7.1.40), l-serine hydroxymethyltransferase (EC 2.1.2.1), threonine aldolase (EC 4.1.2.5), threonine dehydratase (EC 4.3.1.19); transketolase (EC 2.2.1.1), transaldolase (EC 2.2.1.2).

Metabolic anchor reactions for robust biorefining.

Microbial cell factories based on renewable carbon sources are fundamental to a sustainable bio-economy. The economic feasibility of producer cells requires robust performance balancing growth and production. However, the inherent competition between these two objectives often leads to instability and reduces productivity. While algorithms exist to design metabolic network reduction strategies for aligning these objectives, the biochemical basis of the growth-product coupling has remained unresolved. Here, we reveal key reactions in the cellular biochemical repertoire as universal anchor reactions for aligning cell growth and production. A necessary condition for a reaction to be an anchor is that it splits a substrate into two or more molecules. By searching the currently known biochemical reaction space, we identify 62 C-C cleaving anchor reactions, such as isocitrate lyase (EC 4.1.3.1) and L-tryptophan indole-lyase (EC 4.1.99.1), which are relevant for biorefining. The here identified anchor reactions mark network nodes for basing growth-coupled metabolic engineering and novel pathway designs.

Saccharomyces cerevisiae metabolism in ecological context.

The architecture and regulation of Saccharomyces cerevisiae metabolic network are among the best studied owing to its widespread use in both basic research and industry. Yet, several recent studies have revealed notable limitations in explaining genotype-metabolic phenotype relations in this yeast, especially when concerning multiple genetic/environmental perturbations. Apparently unexpected genotype-phenotype relations may originate in the evolutionarily shaped cellular operating principles being hidden in common laboratory conditions. Predecessors of laboratory S. cerevisiae strains, the wild and the domesticated yeasts, have been evolutionarily shaped by highly variable environments, very distinct from laboratory conditions, and most interestingly by social life within microbial communities. Here we present a brief review of the genotypic and phenotypic peculiarities of S. cerevisiae in the context of its social lifestyle beyond laboratory environments. Accounting for this ecological context and the origin of the laboratory strains in experimental design and data analysis would be essential in improving the understanding of genotype-environment-phenotype relationships.

Yeast metabolic chassis designs for diverse biotechnological products.

The diversity of industrially important molecules for which microbial production routes have been experimentally demonstrated is rapidly increasing. The development of economically viable producer cells is, however, lagging behind, as it requires substantial engineering of the host metabolism. A chassis strain suitable for production of a range of molecules is therefore highly sought after but remains elusive. Here, we propose a genome-scale metabolic modeling approach to design chassis strains of Saccharomyces cerevisiae - a widely used microbial cell factory. For a group of 29 products covering a broad range of biochemistry and applications, we identified modular metabolic engineering strategies for re-routing carbon flux towards the desired product. We find distinct product families with shared targets forming the basis for the corresponding chassis cells. The design strategies include overexpression targets that group products by similarity in precursor and cofactor requirements, as well as gene deletion strategies for growth-product coupling that lead to non-intuitive product groups. Our results reveal the extent and the nature of flux re-routing necessary for producing a diverse range of products in a widely used cell factory and provide blueprints for constructing pre-optimized chassis strains.

Hypoxia and GABA shunt activation in the pathogenesis of Alzheimer's disease.

We have previously observed that the conversion of mild cognitive impairment to definitive Alzheimer's disease (AD) is associated with a significant increase in the serum level of 2,4-dihydroxybutyrate (2,4-DHBA). The metabolic generation of 2,4-DHBA is linked to the activation of the γ-aminobutyric acid (GABA) shunt, an alternative energy production pathway activated during cellular stress, when the function of Krebs cycle is compromised. The GABA shunt can be triggered by local hypoperfusion and subsequent hypoxia in AD brains caused by cerebral amyloid angiopathy. Succinic semialdehyde dehydrogenase (SSADH) is a key enzyme in the GABA shunt, converting succinic semialdehyde (SSA) into succinate, a Krebs cycle intermediate. A deficiency of SSADH activity stimulates the conversion of SSA into γ-hydroxybutyrate (GHB), an alternative route from the GABA shunt. GHB can exert not only acute neuroprotective activities but unfortunately also chronic detrimental effects which may lead to cognitive impairment. Subsequently, GHB can be metabolized to 2,4-DHBA and secreted from the brain. Thus, the activation of the GABA shunt and the generation of GHB and 2,4-DHBA can have an important role in the early phase of AD pathogenesis.

Xylose-induced dynamic effects on metabolism and gene expression in engineered Saccharomyces cerevisiae in anaerobic glucose-xylose cultures.

Xylose is present with glucose in lignocellulosic streams available for valorisation to biochemicals. Saccharomyces cerevisiae has excellent characteristics as a host for the bioconversion, except that it strongly prefers glucose to xylose, and the co-consumption remains a challenge. Further, since xylose is not a natural substrate of S. cerevisiae, the regulatory response it induces in an engineered strain cannot be expected to have evolved for its utilisation. Xylose-induced effects on metabolism and gene expression during anaerobic growth of an engineered strain of S. cerevisiae on medium containing both glucose and xylose medium were quantified. The gene expression of S. cerevisiae with an XR-XDH pathway for xylose utilisation was analysed throughout the cultivation: at early cultivation times when mainly glucose was metabolised, at times when xylose was co-consumed in the presence of low glucose concentrations, and when glucose had been depleted and only xylose was being consumed. Cultivations on glucose as a sole carbon source were used as a control. Genome-scale dynamic flux balance analysis models were simulated to analyse the metabolic dynamics of S. cerevisiae. The simulations quantitatively estimated xylose-dependent flux dynamics and challenged the utilisation of the metabolic network. A relative increase in xylose utilisation was predicted to induce the bi-directionality of glycolytic flux and a redox challenge even at low glucose concentrations. Remarkably, xylose was observed to specifically delay the glucose-dependent repression of particular genes in mixed glucose-xylose cultures compared to glucose cultures. The delay occurred at a cultivation time when the metabolic flux activities were similar in the both cultures.