Expanding a dynamic flux balance model of yeast fermentation to genome-scale

Developing a metabolic model-based fed-batch feeding strategy for Pichia pastoris fermentation through fine-tuning of the methanol utilization pathway

Pichia pastoris is a commonly used microbial host for recombinant protein production. It is mostly cultivated in fed-batch mode, in which the environment of the cell is continuously changing. Hence, it is vital to understand the influence of feeding strategy parameters on the intracellular reaction network to fine-tune bioreactor performance. This study used dynamic flux balance analysis (DFBA) integrated with transcriptomics data to simulate the recombinant P. pastoris (Mut^s ) growth during the induction phase for three fed-batch strategies, conducted at constant specific growth rates (μ-stat). The induction phase was split into equal time intervals, and the correlated reactions with protein yield were identified in the three fed-batch strategies using the Pearson correlation coefficient. Subsequently, principal component analysis (PCA) was applied to cluster induction phase time intervals and identify the role of correlated reactions on metabolic differentiation of time intervals. It was found that increasing fluxes through the methanol dissimilation pathway increased protein yield. By adding a methanol assimilation pathway inhibitor (HgCl₂ ) to the shake flask medium growing on glycerol: methanol mixture (10%: 90%, v/v), the protein titre increased by 60%. As per DFBA, the higher the methanol to biomass flux ratio (R_meoh/Δx ), the higher the protein yield. Finally, a novel feeding strategy was developed to increase the amount of R_meoh/Δx compared to the three feeding strategies. The concentration of recombinant human growth hormone (rhGH), used as the model protein, increased by 16% compared to the optimal culture result obtained previously (800 mg L^-1 to 928 mg L^-1 ), while production yield improved by 85% (24.8 mg g_DCW ^-1 to 46 mg g_DCW ^-1 ).

Consideration of Maintenance in Wine Fermentation Modeling

We show that a simple model with a maintenance term can satisfactorily reproduce the simulations of several existing models of wine fermentation from the literature, as well as experimental data. The maintenance describes a consumption of the nitrogen that is not entirely converted into biomass. We show also that considering a maintenance term in the model is equivalent to writing a model with a variable yield that can be estimated from data.

Metabolic Modeling of Wine Fermentation at Genome Scale

Wine fermentation is an ancient biotechnological process mediated by different microorganisms such as yeast and bacteria. Understanding of the metabolic and physiological phenomena taking place during this process can be now attained at a genome scale with the help of metabolic models. In this chapter, we present a detailed protocol for modeling wine fermentation using genome-scale metabolic models. In particular, we illustrate how metabolic fluxes can be computed, optimized and interpreted, for both yeast and bacteria under winemaking conditions. We also show how nutritional requirements can be determined and simulated using these models in relevant test cases. This chapter introduces fundamental concepts and practical steps for applying flux balance analysis in wine fermentation, and as such, it is intended for a broad microbiology audience as well as for practitioners in the metabolic modeling field.

Metabolic flux sampling predicts strain-dependent differences related to aroma production among commercial wine yeasts

BACKGROUND

Metabolomics coupled with genome-scale metabolic modeling approaches have been employed recently to quantitatively analyze the physiological states of various organisms, including Saccharomyces cerevisiae. Although yeast physiology in laboratory strains is well-studied, the metabolic states under industrially relevant scenarios such as winemaking are still not sufficiently understood, especially as there is considerable variation in metabolism between commercial strains. To study the potential causes of strain-dependent variation in the production of volatile compounds during enological conditions, random flux sampling and statistical methods were used, along with experimental extracellular metabolite flux data to characterize the differences in predicted intracellular metabolic states between strains.

RESULTS

It was observed that four selected commercial wine yeast strains (Elixir, Opale, R2, and Uvaferm) produced variable amounts of key volatile organic compounds (VOCs). Principal component analysis was performed on extracellular metabolite data from the strains at three time points of cell cultivation (24, 58, and 144 h). Separation of the strains was observed at all three time points. Furthermore, Uvaferm at 24 h, for instance, was most associated with propanol and ethyl hexanoate. R2 was found to be associated with ethyl acetate and Opale could be associated with isobutanol while Elixir was most associated with phenylethanol and phenylethyl acetate. Constraint-based modeling (CBM) was employed using the latest genome-scale metabolic model of yeast (Yeast8) and random flux sampling was performed with experimentally derived fluxes at various stages of growth as constraints for the model. The flux sampling simulations allowed us to characterize intracellular metabolic flux states and illustrate the key parts of metabolism that likely determine the observed strain differences. Flux sampling determined that Uvaferm and Elixir are similar while R2 and Opale exhibited the highest degree of differences in the Ehrlich pathway and carbon metabolism, thereby causing strain-specific variation in VOC production. The model predictions also established the top 20 fluxes that relate to phenotypic strain variation (e.g. at 24 h). These fluxes indicated that Opale had a higher median flux for pyruvate decarboxylase reactions compared with the other strains. Conversely, R2 which was lower in all VOCs, had higher median fluxes going toward central metabolism. For Elixir and Uvaferm, the differences in metabolism were most evident in fluxes pertaining to transaminase and hexokinase associated reactions. The applied analysis of metabolic divergence unveiled strain-specific differences in yeast metabolism linked to fusel alcohol and ester production.

CONCLUSIONS

Overall, this approach proved useful in elucidating key reactions in amino acid, carbon, and glycerophospholipid metabolism which suggest genetic divergence in activity in metabolic subsystems among these wine strains related to the observed differences in VOC formation. The findings in this study could steer more focused research endeavors in developing or selecting optimal aroma-producing yeast stains for winemaking and other types of alcoholic fermentations.

Isobutanol tolerance and production of Saccharomyces cerevisiae can be improved by engineering its TATA-binding protein Spt15

Low isobutanol tolerance of Saccharomyces cerevisiae limits its application in isobutanol fermentation. Here, we used global transcription machinery engineering to screen mutants with higher isobutanol tolerance and elevated isobutanol titres. TATA-binding protein Spt15 was used as the target of global transcription machinery engineering for improvement of such complex phenotypes. A random mutagenesis library of S. cerevisiae TATA-binding protein Spt15 was constructed and subjected to screening under isobutanol stress. A mutant strain (denoted as spt15-3) with improved isobutanol tolerance was identified. There were three mutations of Spt15 in strain spt15-3, including deletion of A at position -132 nt upstream of initiation codon, insertion of G at position -65 nt upstream of initiation codon and a synonymous mutation at position 315 nt (T → C) downstream of initiation codon. We then metabolically engineered isobutanol synthesis in strains harbouring plasmids YCplac22 containing these Spt15 mutations. Delta integration was used to overexpress ILV3 gene, and 2μ plasmids carrying PGK1p-ILV2 and PGK1p-ARO10 were used to overexpress ILV2 and ARO10 genes. After 24-h micro-aerobic fermentation, Engi-3 produced 0·556 g l^-1 isobutanol, which was 404% and 25·3% greater than isobutanol produced by control Engi-1 and engineered Engi-2, respectively. After 28 h, Engi-4 produced 0·459 g l^-1 isobutanol, which was 315% and 3·2% greater than isobutanol produced Engi-1 and Engi-2, respectively. RNA-Seq-based transcriptome analysis shows that mutations of Spt15 in strain spt15-3 increased the expression of SPT15. Meanwhile, compared with strain Engi-3, the spt15-3 mutation downregulated the expression of genes involved in the TCA cycle and glyoxylic acid cycle, but increased the expression of genes related to cell stability. This work demonstrates that isobutanol tolerance and production of S. cerevisiae can be improved by engineering its TATA-binding protein Spt15. This study clarified the molecular mechanisms regulating isobutanol production and tolerance in S. cerevisiae.

Monod model is insufficient to explain biomass growth in nitrogen-limited yeast fermentation

The yeast Saccharomyces cerevisiae is an essential microorganism in food biotechnology; particularly, in wine and beer making. During wine fermentation, yeasts transform sugars present in the grape juice into ethanol and carbon dioxide. The process occurs in batch conditions and is, for the most part, an anaerobic process. Previous studies linked limited-nitrogen conditions with problematic fermentations, with negative consequences for the performance of the process and the quality of the final product. It is, therefore, of the highest interest to anticipate such problems through mathematical models. Here we propose a model to explain fermentations under nitrogen-limited anaerobic conditions. We separated the biomass formation into two phases: growth and carbohydrate accumulation. Growth was modelled using the well-known Monod equation while carbohydrate accumulation was modelled by an empirical function, analogous to a proportional controller activated by the limitation of available nitrogen. We also proposed to formulate the fermentation rate as a function of the total protein content when relevant data are available. The final model was used to successfully explain experiments taken from the literature, performed under normal and nitrogen-limited conditions. Our results revealed that Monod model is insufficient to explain biomass formation kinetics in nitrogen-limited fermentations of S. cerevisiae. The goodness-of-fit of the herewith proposed model is superior to that of previously published models, offering the means to predict, and thus control fermentations. Importance: Problematic fermentations still occur in the winemaking industrial practise. Problems include sluggish rates of fermentation, which have been linked to insufficient levels of assimilable nitrogen. Data and relevant models can help anticipate poor fermentation performance. In this work, we proposed a model to predict biomass growth and fermentation rate under nitrogen-limited conditions and tested its performance with previously published experimental data. Our results show that the well-known Monod equation does not suffice to explain biomass formation.

A Multiphase Multiobjective Dynamic Genome-Scale Model Shows Different Redox Balancing among Yeast Species of the Saccharomyces Genus in Fermentation

Yeasts constitute over 1,500 species with great potential for biotechnology. Still, the yeast Saccharomyces cerevisiae dominates industrial applications, and many alternative physiological capabilities of lesser-known yeasts are not being fully exploited. While comparative genomics receives substantial attention, little is known about yeasts’ metabolic specificity in batch cultures. Here, we propose a multiphase multiobjective dynamic genome-scale model of yeast batch cultures that describes the uptake of carbon and nitrogen sources and the production of primary and secondary metabolites. The model integrates a specific metabolic reconstruction, based on the consensus Yeast8, and a kinetic model describing the time-varying culture environment. In addition, we proposed a multiphase multiobjective flux balance analysis to compute the dynamics of intracellular fluxes. We then compared the metabolism of S. cerevisiae and Saccharomyces uvarum strains in a rich medium fermentation. The model successfully explained the experimental data and brought novel insights into how cryotolerant strains achieve redox balance. The proposed model (along with the corresponding code) provides a comprehensive picture of the main steps occurring inside the cell during batch cultures and offers a systematic approach to prospect or metabolically engineering novel yeast cell factories.

IMPORTANCE Nonconventional yeast species hold the promise to provide novel metabolic routes to produce industrially relevant compounds and tolerate specific stressors, such as cold temperatures. This work validated the first multiphase multiobjective genome-scale dynamic model to describe carbon and nitrogen metabolism throughout batch fermentation. To test and illustrate its performance, we considered the comparative metabolism of three yeast strains of the Saccharomyces genus in rich medium fermentation. The study revealed that cryotolerant Saccharomyces species might use the γ-aminobutyric acid (GABA) shunt and the production of reducing equivalents as alternative routes to achieve redox balance, a novel biological insight worth being explored further. The proposed model (along with the provided code) can be applied to a wide range of batch processes started with different yeast species and media, offering a systematic and rational approach to prospect nonconventional yeast species metabolism and engineering novel cell factories.

Recent advances in systems and synthetic biology approaches for developing novel cell-factories in non-conventional yeasts

Microbial bioproduction of chemicals, proteins, and primary metabolites from cheap carbon sources is currently an advancing area in industrial research. The model yeast, Saccharomyces cerevisiae, is a well-established biorefinery host that has been used extensively for commercial manufacturing of bioethanol from myriad carbon sources. However, its Crabtree-positive nature often limits the use of this organism for the biosynthesis of commercial molecules that do not belong in the fermentative pathway. To avoid extensive strain engineering of S. cerevisiae for the production of metabolites other than ethanol, non-conventional yeasts can be selected as hosts based on their natural capacity to produce desired commodity chemicals. Non-conventional yeasts like Kluyveromyces marxianus, K. lactis, Yarrowia lipolytica, Pichia pastoris, Scheffersomyces stipitis, Hansenula polymorpha, and Rhodotorula toruloides have been considered as potential industrial eukaryotic hosts owing to their desirable phenotypes such as thermotolerance, assimilation of a wide range of carbon sources, as well as ability to secrete high titers of protein and lipid. However, the advanced metabolic engineering efforts in these organisms are still lacking due to the limited availability of systems and synthetic biology methods like in silico models, well-characterised genetic parts, and optimized genome engineering tools. This review provides an insight into the recent advances and challenges of systems and synthetic biology as well as metabolic engineering endeavours towards the commercial usage of non-conventional yeasts. Particularly, the approaches in emerging non-conventional yeasts for the production of enzymes, therapeutic proteins, lipids, and metabolites for commercial applications are extensively discussed here. Various attempts to address current limitations in designing novel cell factories have been highlighted that include the advances in the fields of genome-scale metabolic model reconstruction, flux balance analysis, 'omics'-data integration into models, genome-editing toolkit development, and rewiring of cellular metabolisms for desired chemical production. Additionally, the understanding of metabolic networks using ¹³C-labelling experiments as well as the utilization of metabolomics in deciphering intracellular fluxes and reactions have also been discussed here. Application of cutting-edge nuclease-based genome editing platforms like CRISPR/Cas9, and its optimization towards efficient strain engineering in non-conventional yeasts have also been described. Additionally, the impact of the advances in promising non-conventional yeasts for efficient commercial molecule synthesis has been meticulously reviewed. In the future, a cohesive approach involving systems and synthetic biology will help in widening the horizon of the use of unexplored non-conventional yeast species towards industrial biotechnology.

Curation and Analysis of a Saccharomyces cerevisiae Genome-Scale Metabolic Model for Predicting Production of Sensory Impact Molecules under Enological Conditions

One approach for elucidating strain-to-strain metabolic differences is the use of genome-scale metabolic models (GSMMs). To date GSMMs have not focused on the industrially important area of flavor production and, as such; do not cover all the pathways relevant to flavor formation in yeast. Moreover, current models for Saccharomyces cerevisiae generally focus on carbon-limited and/or aerobic systems, which is not pertinent to enological conditions. Here, we curate a GSMM (iWS902) to expand on the existing Ehrlich pathway and ester formation pathways central to aroma formation in industrial winemaking, in addition to the existing sulfur metabolism and medium-chain fatty acid (MCFA) pathways that also contribute to production of sensory impact molecules. After validating the model using experimental data, we predict key differences in metabolism for a strain (EC 1118) in two distinct growth conditions, including differences for aroma impact molecules such as acetic acid, tryptophol, and hydrogen sulfide. Additionally, we propose novel targets for metabolic engineering for aroma profile modifications employing flux variability analysis with the expanded GSMM. The model provides mechanistic insights into the key metabolic pathways underlying aroma formation during alcoholic fermentation and provides a potential framework to contribute to new strategies to optimize the aroma of wines.

QTL mapping of modelled metabolic fluxes reveals gene variants impacting yeast central carbon metabolism

The yeast Saccharomyces cerevisiae is an attractive industrial microorganism for the production of foods and beverages as well as for various bulk and fine chemicals, such as biofuels or fragrances. Building blocks for these biosyntheses are intermediates of yeast central carbon metabolism (CCM), whose intracellular availability depends on balanced single reactions that form metabolic fluxes. Therefore, efficient product biosynthesis is influenced by the distribution of these fluxes. We recently demonstrated great variations in CCM fluxes between yeast strains of different origins. However, we have limited understanding of flux modulation and the genetic basis of flux variations. In this study, we investigated the potential of quantitative trait locus (QTL) mapping to elucidate genetic variations responsible for differences in metabolic flux distributions (fQTL). Intracellular metabolic fluxes were estimated by constraint-based modelling and used as quantitative phenotypes, and differences in fluxes were linked to genomic variations. Using this approach, we detected four fQTLs that influence metabolic pathways. The molecular dissection of these QTLs revealed two allelic gene variants, PDB1 and VID30, contributing to flux distribution. The elucidation of genetic determinants influencing metabolic fluxes, as reported here for the first time, creates new opportunities for the development of strains with optimized metabolite profiles for various applications.

Engineering microbial chemical factories using metabolic models

Living organisms in analogy with chemical factories use simple molecules such as sugars to produce a variety of compounds which are necessary for sustaining life and some of which are also commercially valuable. The metabolisms of simple (such as bacteria) and higher organisms (such as plants) alike can be exploited to convert low value inputs into high value outputs. Unlike conventional chemical factories, microbial production chassis are not necessarily tuned for a single product overproduction. Despite the same end goal, metabolic and industrial engineers rely on different techniques for achieving productivity goals. Metabolic engineers cannot affect reaction rates by manipulating pressure and temperature, instead they have at their disposal a range of enzymes and transcriptional and translational processes to optimize accordingly. In this review, we first highlight how various analytical approaches used in metabolic engineering and synthetic biology are related to concepts developed in systems and control engineering. Specifically, how algorithmic concepts derived in operations research can help explain the structure and organization of metabolic networks. Finally, we consider the future directions and challenges faced by the field of metabolic network modeling and the possible contributions of concepts drawn from the classical fields of chemical and control engineering. The aim of the review is to offer a current perspective of metabolic engineering and all that it entails without requiring specialized knowledge of bioinformatics or systems biology.

Clostridium butyricum population balance model: Predicting dynamic metabolic flux distributions using an objective function related to extracellular glycerol content

BACKGROUND

Extensive experimentation has been conducted to increment 1,3-propanediol (PDO) production using Clostridium butyricum cultures in glycerol, but computational predictions are limited. Previously, we reconstructed the genome-scale metabolic (GSM) model iCbu641, the first such model of a PDO-producing Clostridium strain, which was validated at steady state using flux balance analysis (FBA). However, the prediction ability of FBA is limited for batch and fed-batch cultures, which are the most often employed industrial processes.

RESULTS

We used the iCbu641 GSM model to develop a dynamic flux balance analysis (DFBA) approach to predict the PDO production of the Colombian strain Clostridium sp IBUN 158B. First, we compared the predictions of the dynamic optimization approach (DOA), static optimization approach (SOA), and direct approach (DA). We found no differences between approaches, but the DOA simulation duration was nearly 5000 times that of the SOA and DA simulations. Experimental results at glycerol limitation and glycerol excess allowed for validating dynamic predictions of growth, glycerol consumption, and PDO formation. These results indicated a 4.4% error in PDO prediction and therefore validated the previously proposed objective functions. We performed two global sensitivity analyses, finding that the kinetic input parameters of glycerol uptake flux had the most significant effect on PDO predictions. The other input parameters evaluated during global sensitivity analysis were biomass composition (precursors and macromolecules), death constants, and the kinetic parameters of acetic acid secretion flux. These last input parameters, all obtained from other Clostridium butyricum cultures, were used to develop a population balance model (PBM). Finally, we simulated fed-batch cultures, predicting a final PDO production near to 66 g/L, almost three times the PDO predicted in the best batch culture.

CONCLUSIONS

We developed and validated a dynamic approach to predict PDO production using the iCbu641 GSM model and the previously proposed objective functions. This validated approach was used to propose a population model and then an increment in predictions of PDO production through fed-batch cultures. Therefore, this dynamic model could predict different scenarios, including its integration into downstream processes to predict technical-economic feasibilities and reducing the time and costs associated with experimentation.

Grape and Wine Metabolomics to Develop New Insights Using Untargeted and Targeted Approaches

Chemical analysis of grape juice and wine has been performed for over 50 years in a targeted manner to determine a limited number of compounds using Gas Chromatography, Mass-Spectrometry (GC-MS) and High Pressure Liquid Chromatography (HPLC). Therefore, it only allowed the determination of metabolites that are present in high concentration, including major sugars, amino acids and some important carboxylic acids. Thus, the roles of many significant but less concentrated metabolites during wine making process are still not known. This is where metabolomics shows its enormous potential, mainly because of its capability in analyzing over 1000 metabolites in a single run due to the recent advancements of high resolution and sensitive analytical instruments. Metabolomics has predominantly been adopted by many wine scientists as a hypothesis-generating tool in an unbiased and non-targeted way to address various issues, including characterization of geographical origin (terroir) and wine yeast metabolic traits, determination of biomarkers for aroma compounds, and the monitoring of growth developments of grape vines and grapes. The aim of this review is to explore the published literature that made use of both targeted and untargeted metabolomics to study grapes and wines and also the fermentation process. In addition, insights are also provided into many other possible avenues where metabolomics shows tremendous potential as a question-driven approach in grape and wine research.

Engineering TATA-binding protein Spt15 to improve ethanol tolerance and production in Kluyveromyces marxianus

BACKGROUND

Low ethanol tolerance of Kluyveromyces marxianus limits its application in high-temperature ethanol fermentation. As a complex phenotype, ethanol tolerance involves synergistic actions of many genes that are widely distributed throughout the genome, thereby being difficult to engineer. TATA-binding protein is the most common target of global transcription machinery engineering for improvement of complex phenotypes.

RESULTS

A random mutagenesis library of K. marxianus TATA-binding protein Spt15 was constructed and subjected to screening under ethanol stress. Two mutant strains with improved ethanol tolerance were identified, one of which (denoted as M2) exhibited increased ethanol productivity. The mutant of Spt15 in strain M2 (denoted as Spt15-M2) has a single amino acid substitution at position 31 (Lys → Glu). RNA-Seq-based transcriptomic analysis revealed cellular transcription profile changes resulting from Spt15-M2. Spt15-M2 caused changes in transcriptional level of most of the genes in the central carbon metabolism network. Compared with control strain, 444 differentially expressed genes (DEGs) were identified in strain M2 (fold change > 2, Padj < 0.05), including 48 up-regulated and 396 down-regulated. The up-regulated DEGs are involved in amino acid transport, long-chain fatty acid biosynthesis and MAPK signaling pathway, while the down-regulated DEGs are related to ribosome biogenesis, translation and protein synthesis. Five candidate genes (GAP1, GNP1, FAR1, STE2 and TEC1), which were found to be up-regulated in M2 strain, were overexpressed for a gain-of-function assay. However, the overexpression of no single gene helped improve ethanol tolerance as SPT15-M2 did.

CONCLUSIONS

This work demonstrates that ethanol tolerance of K. marxianus can be improved by engineering its TATA-binding protein. A single amino acid substitution (K31E) of TATA-binding protein Spt15 is able to bring differential expression of hundreds of genes that acted as an interconnected network for the phenotype of ethanol tolerance. Future perspectives of this technique in K. marxianus were discussed.

Saccharomyces cerevisiae and S. kudriavzevii Synthetic Wine Fermentation Performance Dissected by Predictive Modeling

Wineries face unprecedented challenges due to new market demands and climate change effects on wine quality. New yeast starters including non-conventional Saccharomyces species, such as S. kudriavzevii, may contribute to deal with some of these challenges. The design of new fermentations using non-conventional yeasts requires an improved understanding of the physiology and metabolism of these cells. Dynamic modeling brings the potential of exploring the most relevant mechanisms and designing optimal processes more systematically. In this work we explore mechanisms by means of a model selection, reduction and cross-validation pipeline which enables to dissect the most relevant fermentation features for the species under consideration, Saccharomyces cerevisiae T73 and Saccharomyces kudriavzevii CR85. The pipeline involved the comparison of a collection of models which incorporate several alternative mechanisms with emphasis on the inhibitory effects due to temperature and ethanol. We focused on defining a minimal model with the minimum number of parameters, to maximize the identifiability and the quality of cross-validation. The selected model was then used to highlight differences in behavior between species. The analysis of model parameters would indicate that the specific growth rate and the transport of hexoses at initial times are higher for S. cervisiae T73 while S. kudriavzevii CR85 diverts more flux for glycerol production and cellular maintenance. As a result, the fermentations with S. kudriavzevii CR85 are typically slower; produce less ethanol but higher glycerol. Finally, we also explored optimal initial inoculation and process temperature to find the best compromise between final product characteristics and fermentation duration. Results reveal that the production of glycerol is distinctive in S. kudriavzevii CR85, it was not possible to achieve the same production of glycerol with S. cervisiae T73 in any of the conditions tested. This result brings the idea that the optimal design of mixed cultures may have an enormous potential for the improvement of final wine quality.

Genome-Scale Reconstruction of the Metabolic Network in Oenococcus oeni to Assess Wine Malolactic Fermentation

Oenococcus oeni is the main responsible agent for malolactic fermentation in wine, an unpredictable and erratic process in winemaking. To address this, we have constructed and exhaustively curated the first genome-scale metabolic model of Oenococcus oeni, comprising 660 reactions, 536 metabolites and 454 genes. In silico experiments revealed that nutritional requirements are predicted with an accuracy of 93%, while 14 amino acids were found to be essential for the growth of this bacterial species. When the model was applied to determine the non-growth associated maintenance, results showed that O. oeni grown at 12% ethanol concentration spent 30 times more ATP to stay alive than in the absence of ethanol. Most of this ATP is employed for extruding protons outside of the cell. A positive relationship was also found between specific consumption rates of fructose, amino acids, oxygen, and malic acid and the specific production rates of erythritol, lactate, and acetate, according to the ethanol content of the medium. The metabolic model reconstructed here represents a unique tool to predict the successful completion of wine malolactic fermentation carried out either by different strains of Oenococcus oeni, as well as at any particular physico-chemical composition of wine. It will also allow the development of consortium metabolic models that could be applied to winemaking to simulate and understand the interactions between O. oeni and other microorganisms that share this ecological niche.

Dynamic genome-scale metabolic modeling of the yeast Pichia pastoris

BACKGROUND

Pichia pastoris shows physiological advantages in producing recombinant proteins, compared to other commonly used cell factories. This yeast is mostly grown in dynamic cultivation systems, where the cell's environment is continuously changing and many variables influence process productivity. In this context, a model capable of explaining and predicting cell behavior for the rational design of bioprocesses is highly desirable. Currently, there are five genome-scale metabolic reconstructions of P. pastoris which have been used to predict extracellular cell behavior in stationary conditions.

RESULTS

In this work, we assembled a dynamic genome-scale metabolic model for glucose-limited, aerobic cultivations of Pichia pastoris. Starting from an initial model structure for batch and fed-batch cultures, we performed pre/post regression diagnostics to ensure that model parameters were identifiable, significant and sensitive. Once identified, the non-relevant ones were iteratively fixed until a priori robust modeling structures were found for each type of cultivation. Next, the robustness of these reduced structures was confirmed by calibrating the model with new datasets, where no sensitivity, identifiability or significance problems appeared in their parameters. Afterwards, the model was validated for the prediction of batch and fed-batch dynamics in the studied conditions. Lastly, the model was employed as a case study to analyze the metabolic flux distribution of a fed-batch culture and to unravel genetic and process engineering strategies to improve the production of recombinant Human Serum Albumin (HSA). Simulation of single knock-outs indicated that deviation of carbon towards cysteine and tryptophan formation improves HSA production. The deletion of methylene tetrahydrofolate dehydrogenase could increase the HSA volumetric productivity by 630%. Moreover, given specific bioprocess limitations and strain characteristics, the model suggests that implementation of a decreasing specific growth rate during the feed phase of a fed-batch culture results in a 25% increase of the volumetric productivity of the protein.

CONCLUSION

In this work, we formulated a dynamic genome scale metabolic model of Pichia pastoris that yields realistic metabolic flux distributions throughout dynamic cultivations. The model can be calibrated with experimental data to rationally propose genetic and process engineering strategies to improve the performance of a P. pastoris strain of interest.

Management of Multiple Nitrogen Sources during Wine Fermentation by Saccharomyces cerevisiae

During fermentative growth in natural and industrial environments, Saccharomyces cerevisiae must redistribute the available nitrogen from multiple exogenous sources to amino acids in order to suitably fulfill anabolic requirements. To exhaustively explore the management of this complex resource, we developed an advanced strategy based on the reconciliation of data from a set of stable isotope tracer experiments with labeled nitrogen sources. Thus, quantifying the partitioning of the N compounds through the metabolism network during fermentation, we demonstrated that, contrary to the generally accepted view, only a limited fraction of most of the consumed amino acids is directly incorporated into proteins. Moreover, substantial catabolism of these molecules allows for efficient redistribution of nitrogen, supporting the operative de novo synthesis of proteinogenic amino acids. In contrast, catabolism of consumed amino acids plays a minor role in the formation of volatile compounds. Another important feature is that the α-keto acid precursors required for the de novo syntheses originate mainly from the catabolism of sugars, with a limited contribution from the anabolism of consumed amino acids. This work provides a comprehensive view of the intracellular fate of consumed nitrogen sources and the metabolic origin of proteinogenic amino acids, highlighting a strategy of distribution of metabolic fluxes implemented by yeast as a means of adapting to environments with changing and scarce nitrogen resources.IMPORTANCE A current challenge for the wine industry, in view of the extensive competition in the worldwide market, is to meet consumer expectations regarding the sensory profile of the product while ensuring an efficient fermentation process. Understanding the intracellular fate of the nitrogen sources available in grape juice is essential to the achievement of these objectives, since nitrogen utilization affects both the fermentative activity of yeasts and the formation of flavor compounds. However, little is known about how the metabolism operates when nitrogen is provided as a composite mixture, as in grape must. Here we quantitatively describe the distribution through the yeast metabolic network of the N moieties and C backbones of these nitrogen sources. Knowledge about the management of a complex resource, which is devoted to improvement of the use of the scarce N nutrient for growth, will be useful for better control of the fermentation process and the sensory quality of wines.

Effect of amino acid supplementation on titer and glycosylation distribution in hybridoma cell cultures-Systems biology-based interpretation using genome-scale metabolic flux balance model and multivariate data analysis

Genome-scale flux balance analysis (FBA) is a powerful systems biology tool to characterize intracellular reaction fluxes during cell cultures. FBA estimates intracellular reaction rates by optimizing an objective function, subject to the constraints of a metabolic model and media uptake/excretion rates. A dynamic extension to FBA, dynamic flux balance analysis (DFBA), can calculate intracellular reaction fluxes as they change during cell cultures. In a previous study by Read et al. (2013), a series of informed amino acid supplementation experiments were performed on twelve parallel murine hybridoma cell cultures, and this data was leveraged for further analysis (Read et al., Biotechnol Prog. 2013;29:745-753). In order to understand the effects of media changes on the model murine hybridoma cell line, a systems biology approach is applied in the current study. Dynamic flux balance analysis was performed using a genome-scale mouse metabolic model, and multivariate data analysis was used for interpretation. The calculated reaction fluxes were examined using partial least squares and partial least squares discriminant analysis. The results indicate media supplementation increases product yield because it raises nutrient levels extending the growth phase, and the increased cell density allows for greater culture performance. At the same time, the directed supplementation does not change the overall metabolism of the cells. This supports the conclusion that product quality, as measured by glycoform assays, remains unchanged because the metabolism remains in a similar state. Additionally, the DFBA shows that metabolic state varies more at the beginning of the culture but less by the middle of the growth phase, possibly due to stress on the cells during inoculation. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:1163-1173, 2016.

Genome scale models of yeast: towards standardized evaluation and consistent omic integration

Genome scale models (GEMs) have enabled remarkable advances in systems biology, acting as functional databases of metabolism, and as scaffolds for the contextualization of high-throughput data. In the case of Saccharomyces cerevisiae (budding yeast), several GEMs have been published and are currently used for metabolic engineering and elucidating biological interactions. Here we review the history of yeast's GEMs, focusing on recent developments. We study how these models are typically evaluated, using both descriptive and predictive metrics. Additionally, we analyze the different ways in which all levels of omics data (from gene expression to flux) have been integrated in yeast GEMs. Relevant conclusions and current challenges for both GEM evaluation and omic integration are highlighted.

Analyzing redox balance in a synthetic yeast platform to improve utilization of brown macroalgae as feedstock

Macroalgae have high potential to be an efficient, and sustainable feedstock for the production of biofuels and other more valuable chemicals. Attempts have been made to enable the co-fermentation of alginate and mannitol by Saccharomyces cerevisiae to unlock the full potential of this marine biomass. However, the efficient use of the sugars derived from macroalgae depends on the equilibrium of cofactors derived from the alginate and mannitol catabolic pathways. There are a number of strong metabolic limitations that have to be tackled before this bioconversion can be carried out efficiently by engineered yeast cells.

An analysis of the redox balance during ethanol fermentation from alginate and mannitol by Saccharomyces cerevisiae using metabolic engineering tools was carried out. To represent the strain designed for conversion of macroalgae carbohydrates to ethanol, a context-specific model was derived from the available yeast genome-scale metabolic reconstructions. Flux balance analysis and dynamic simulations were used to determine the flux distributions. The model indicates that ethanol production is determined by the activity of 4-deoxy-l-erythro-5-hexoseulose uronate (DEHU) reductase (DehR) and its preferences for NADH or NADPH which influences strongly the flow of cellular resources. Different scenarios were explored to determine the equilibrium between NAD(H) and NADP(H) that will lead to increased ethanol yields on mannitol and DEHU under anaerobic conditions. When rates of mannitol dehydrogenase and DehR_NADH tend to be close to a ratio in the range 1–1.6, high growth rates and ethanol yields were predicted. The analysis shows a number of metabolic limitations that are not easily identified through experimental procedures such as quantifying the impact of the cofactor preference by DEHU reductase in the system, the low flux into the alginate catabolic pathway, and a detailed analysis of the redox balance. These results show that production of ethanol and other chemicals can be optimized if a redox balance is achieved. A possible methodology to achieve this balance is presented. This paper shows how metabolic engineering tools are essential to comprehend and overcome this limitation.

•

We studied a strain designed for bioconversion of macroalgae sugars to ethanol.

•

A genome-scale model was used to simulate biomass and by-product formation.

•

The characterization of the current metabolic state of the strain was achieved.

•

Biofuel production depends on the redox balance derived from alginate and mannitol.

•

Flux split into DehR determines the redox balance, by-products and ethanol level.

Dynamic flux balancing elucidates NAD(P)H production as limiting response to furfural inhibition in Saccharomyces cerevisiae

Achieving efficient and economical lignocellulose-based bioprocess requires a robust organism tolerant to furfural, a major inhibitory compound present in lignocellulosic hydrolysate. The aim of this study was to develop a model that could generate quantitative descriptions of cell metabolism for elucidating the cell's adaptive response to furfural. Such a modelling tool could provide strategies for the design of more robust cells. A dynamic flux balance (dFBA) model of Saccharomyces cerevisiae was created by coupling a kinetic fermentation model with a previously published genome-scale stoichiometric model. The dFBA model was used for studying intracellular and extracellular flux responses to furfural perturbations under steady state and dynamic conditions. The predicted effects of furfural on dynamic flux profiles agreed well with previously published experimental results. The model showed that the yeast cell adjusts its metabolism in response to furfural challenge by increasing fluxes through the pentose phosphate pathway, TCA cycle, and proline and serine biosynthesis in order to meet the high demand of NAD(P)H cofactors. The model described here can be used to aid in systematic optimization of the yeast, as well as of the fermentation process, for efficient lignocellulosic ethanol production.

Predictive sulfur metabolism - a field in flux

The key role of sulfur metabolites in response to biotic and abiotic stress in plants, as well as their importance in diet and health has led to a significant interest and effort in trying to understand and manipulate the production of relevant compounds. Metabolic engineering utilizes a set of theoretical tools to help rationally design modifications that enhance the production of a desired metabolite. Such approaches have proven their value in bacterial systems, however, the paucity of success stories to date in plants, suggests that challenges remain. Here, we review the most commonly used methods for understanding metabolic flux, focusing on the sulfur assimilatory pathway. We highlight known issues with both experimental and theoretical approaches, as well as presenting recent methods for integrating different modeling strategies, and progress toward an understanding of flux at the whole plant level.

Use of chemostat cultures mimicking different phases of wine fermentations as a tool for quantitative physiological analysis

Background

Saccharomyces cerevisiae is the most relevant yeast species conducting the alcoholic fermentation that takes place during winemaking. Although the physiology of this model organism has been extensively studied, systematic quantitative physiology studies of this yeast under winemaking conditions are still scarce, thus limiting the understanding of fermentative metabolism of wine yeast strains and the systematic description, modelling and prediction of fermentation processes. In this study, we implemented and validated the use of chemostat cultures as a tool to simulate different stages of a standard wine fermentation, thereby allowing to implement metabolic flux analyses describing the sequence of metabolic states of S. cerevisae along the wine fermentation.

Results

Chemostat cultures mimicking the different stages of standard wine fermentations of S. cerevisiae EC1118 were performed using a synthetic must and strict anaerobic conditions. The simulated stages corresponded to the onset of the exponential growth phase, late exponential growth phase and cells just entering stationary phase, at dilution rates of 0.27, 0.04, 0.007 h⁻¹, respectively. Notably, measured substrate uptake and product formation rates at each steady state condition were generally within the range of corresponding conversion rates estimated during the different batch fermentation stages.

Moreover, chemostat data were further used for metabolic flux analysis, where biomass composition data for each condition was considered in the stoichiometric model. Metabolic flux distributions were coherent with previous analyses based on batch cultivations data and the pseudo-steady state assumption.

Conclusions

Steady state conditions obtained in chemostat cultures reflect the environmental conditions and physiological states of S. cerevisiae corresponding to the different growth stages of a typical batch wine fermentation, thereby showing the potential of this experimental approach to systematically study the effect of environmental relevant factors such as temperature, sugar concentration, C/N ratio or (micro) oxygenation on the fermentative metabolism of wine yeast strains.

Age-correlated stress resistance improves fitness of yeast: support from agent-based simulations

BACKGROUND

Resistance to stress is often heterogeneous among individuals within a population, which helps protect against intermittent stress (bet hedging). This is also the case for heat shock resistance in the budding yeast Saccharomyces cerevisiae. Interestingly, the resistance appears to be continuously distributed (vs. binary, switch-like) and correlated with replicative age (vs. random). Older, slower-growing cells are more resistant than younger, faster-growing ones. Is there a fitness benefit to age-correlated stress resistance?

RESULTS

Here this hypothesis is explored using a simple agent-based model, which simulates a population of individual cells that grow and replicate. Cells age by accumulating damage, which lowers their growth rate. They synthesize trehalose at a metabolic cost, which helps protect against heat shock. Proteins Tsl1 and Tps3 (trehalose synthase complex regulatory subunit TSL1 and TPS3) represent the trehalose synthesis complex and they are expressed using constant, age-dependent and stochastic terms. The model was constrained by calibration and comparison to data from the literature, including individual-based observations obtained using high-throughput microscopy and flow cytometry. A heterogeneity network was developed, which highlights the predominant sources and pathways of resistance heterogeneity. To determine the best trehalose synthesis strategy, model strains with different Tsl1/Tps3 expression parameters were placed in competition in an environment with intermittent heat shocks.

CONCLUSIONS

For high severities and low frequencies of heat shock, the winning strain used an age-dependent bet hedging strategy, which shows that there can be a benefit to age-correlated stress resistance. The study also illustrates the utility of combining individual-based observations and modeling to understand mechanisms underlying population heterogeneity, and the effect on fitness.

Systematic construction of kinetic models from genome-scale metabolic networks

The quantitative effects of environmental and genetic perturbations on metabolism can be studied in silico using kinetic models. We present a strategy for large-scale model construction based on a logical layering of data such as reaction fluxes, metabolite concentrations, and kinetic constants. The resulting models contain realistic standard rate laws and plausible parameters, adhere to the laws of thermodynamics, and reproduce a predefined steady state. These features have not been simultaneously achieved by previous workflows. We demonstrate the advantages and limitations of the workflow by translating the yeast consensus metabolic network into a kinetic model. Despite crudely selected data, the model shows realistic control behaviour, a stable dynamic, and realistic response to perturbations in extracellular glucose concentrations. The paper concludes by outlining how new data can continuously be fed into the workflow and how iterative model building can assist in directing experiments.

Metabolic flux analysis during the exponential growth phase of Saccharomyces cerevisiae in wine fermentations

As a consequence of the increase in global average temperature, grapes with the adequate phenolic and aromatic maturity tend to be overripe by the time of harvest, resulting in increased sugar concentrations and imbalanced C/N ratios in fermenting musts. This fact sets obvious additional hurdles in the challenge of obtaining wines with reduced alcohols levels, a new trend in consumer demands. It would therefore be interesting to understand Saccharomyces cerevisiae physiology during the fermentation of must with these altered characteristics. The present study aims to determine the distribution of metabolic fluxes during the yeast exponential growth phase, when both carbon and nitrogen sources are in excess, using continuous cultures. Two different sugar concentrations were studied under two different winemaking temperature conditions. Although consumption and production rates for key metabolites were severely affected by the different experimental conditions studied, the general distribution of fluxes in central carbon metabolism was basically conserved in all cases. It was also observed that temperature and sugar concentration exerted a higher effect on the pentose phosphate pathway and glycerol formation than on glycolysis and ethanol production. Additionally, nitrogen uptake, both quantitatively and qualitatively, was strongly influenced by environmental conditions. This work provides the most complete stoichiometric model used for Metabolic Flux Analysis of S. cerevisiae in wine fermentations employed so far, including the synthesis and release of relevant aroma compounds and could be used in the design of optimal nitrogen supplementation of wine fermentations.

Characterization of the metabolic requirements in yeast meiosis

The diploid yeast Saccharomyces cerevisiae undergoes mitosis in glucose-rich medium but enters meiosis in acetate sporulation medium. The transition from mitosis to meiosis involves a remarkable adaptation of the metabolic machinery to the changing environment to meet new energy and biosynthesis requirements. Biochemical studies indicate that five metabolic pathways are active at different stages of sporulation: glutamate formation, tricarboxylic acid cycle, glyoxylate cycle, gluconeogenesis, and glycogenolysis. A dynamic synthesis of macromolecules, including nucleotides, amino acids, and lipids, is also observed. However, the metabolic requirements of sporulating cells are poorly understood. In this study, we apply flux balance analyses to uncover optimal principles driving the operation of metabolic networks over the entire period of sporulation. A meiosis-specific metabolic network is constructed, and flux distribution is simulated using ten objective functions combined with time-course expression-based reaction constraints. By systematically evaluating the correlation between computational and experimental fluxes on pathways and macromolecule syntheses, the metabolic requirements of cells are determined: sporulation requires maximization of ATP production and macromolecule syntheses in the early phase followed by maximization of carbohydrate breakdown and minimization of ATP production in the middle and late stages. Our computational models are validated by in silico deletion of enzymes known to be essential for sporulation. Finally, the models are used to predict novel metabolic genes required for sporulation. This study indicates that yeast cells have distinct metabolic requirements at different phases of meiosis, which may reflect regulation that realizes the optimal outcome of sporulation. Our meiosis-specific network models provide a framework for an in-depth understanding of the roles of enzymes and reactions, and may open new avenues for engineering metabolic pathways to improve sporulation efficiency.

Metabolic modelling in the development of cell factories by synthetic biology

Cell factories are commonly microbial organisms utilized for bioconversion of renewable resources to bulk or high value chemicals. Introduction of novel production pathways in chassis strains is the core of the development of cell factories by synthetic biology. Synthetic biology aims to create novel biological functions and systems not found in nature by combining biology with engineering. The workflow of the development of novel cell factories with synthetic biology is ideally linear which will be attainable with the quantitative engineering approach, high-quality predictive models, and libraries of well-characterized parts. Different types of metabolic models, mathematical representations of metabolism and its components, enzymes and metabolites, are useful in particular phases of the synthetic biology workflow. In this minireview, the role of metabolic modelling in synthetic biology will be discussed with a review of current status of compatible methods and models for the in silico design and quantitative evaluation of a cell factory.

Genomic analysis of Saccharomyces cerevisiae isolates that grow optimally with glucose as the sole carbon source

A population of Saccharomyces cerevisiae was cultured for approximately 450 generations in the presence of high glucose to select for genetic variants that grew optimally under these conditions. Using the parental strain BY4741 as the starting population, an evolved culture was obtained after aerobic growth in a high glucose medium for approximately 450 generations. After the evolution period, three single colony isolates were selected for analysis. Next-generation Ion Torrent sequencing was used to evaluate genetic changes. Greater than 100 deletion/insertion changes were found with approximately half of these effecting genes. Additionally, over 180 SNPs were identified with more than one-quarter of these resulting in a nonsynonymous mutation. Affymetrix DNA microarrays and RNseq analysis were used to determine differences in gene expression in the evolved strains compared to the parental strain. It was established that approximately 900 genes demonstrated significantly altered expression in the evolved strains relative to the parental strain. Many of these genes showed similar alterations in their expression in all three evolved strains. Interestingly, genes with altered expression in the three evolved strains included genes with a role in oxidative metabolism. Overall these results are consistent with the physiological observations of optimal growth with glucose as the carbon source. Namely, the decreased ethanol production suggest that the underlying metabolism switched from fermentation to respiration during the selection for optimal growth on glucose.