Paula Jouhten

A multi-level study of recombinant Pichia pastoris in different oxygen conditions

BackgroundYeasts are attractive expression platforms for many recombinant proteins, and there is evidence for an important interrelation between the protein secretion machinery and environmental stresses. While adaptive responses to such stresses are extensively studied in Saccharomyces cerevisiae, little is known about their impact on the physiology of Pichia pastoris. We have recently reported a beneficial effect of hypoxia on recombinant Fab secretion in P. pastoris chemostat cultivations. As a consequence, a systems biology approach was used to comprehensively identify cellular adaptations to low oxygen availability and the additional burden of protein production. Gene expression profiling was combined with proteomic analyses and the 13C isotope labelling based experimental determination of metabolic fluxes in the central carbon metabolism.ResultsThe physiological adaptation of P. pastoris to hypoxia showed distinct traits in relation to the model yeast S. cerevisiae. There was a positive correlation between the transcriptomic, proteomic and metabolic fluxes adaptation of P. pastoris core metabolism to hypoxia, yielding clear evidence of a strong transcriptional regulation component of key pathways such as glycolysis, pentose phosphate pathway and TCA cycle. In addition, the adaptation to reduced oxygen revealed important changes in lipid metabolism, stress responses, as well as protein folding and trafficking.ConclusionsThis systems level study helped to understand the physiological adaptations of cellular mechanisms to low oxygen availability in a recombinant P. pastoris strain. Remarkably, the integration of data from three different levels allowed for the identification of differences in the regulation of the core metabolism between P. pastoris and S. cerevisiae. Detailed comparative analysis of the transcriptomic data also led to new insights into the gene expression profiles of several cellular processes that are not only susceptible to low oxygen concentrations, but might also contribute to enhanced protein secretion.

Oxygen dependence of metabolic fluxes and energy generation of Saccharomyces cerevisiae CEN.PK113-1A

BackgroundThe yeast Saccharomyces cerevisiae is able to adjust to external oxygen availability by utilizing both respirative and fermentative metabolic modes. Adjusting the metabolic mode involves alteration of the intracellular metabolic fluxes that are determined by the cells multilevel regulatory network. Oxygen is a major determinant of the physiology of S. cerevisiae but understanding of the oxygen dependence of intracellular flux distributions is still scarce.ResultsMetabolic flux distributions of S. cerevisiae CEN.PK113-1A growing in glucose-limited chemostat cultures at a dilution rate of 0.1 h-1 with 20.9%, 2.8%, 1.0%, 0.5% or 0.0% O2 in the inlet gas were quantified by 13C-MFA. Metabolic flux ratios from fractional [U-13C]glucose labelling experiments were used to solve the underdetermined MFA system of central carbon metabolism of S. cerevisiae.While ethanol production was observed already in 2.8% oxygen, only minor differences in the flux distribution were observed, compared to fully aerobic conditions. However, in 1.0% and 0.5% oxygen the respiratory rate was severely restricted, resulting in progressively reduced fluxes through the TCA cycle and the direction of major fluxes to the fermentative pathway. A redistribution of fluxes was observed in all branching points of central carbon metabolism. Yet only when oxygen provision was reduced to 0.5%, was the biomass yield exceeded by the yields of ethanol and CO2. Respirative ATP generation provided 59% of the ATP demand in fully aerobic conditions and still a substantial 25% in 0.5% oxygenation. An extensive redistribution of fluxes was observed in anaerobic conditions compared to all the aerobic conditions. Positive correlation between the transcriptional levels of metabolic enzymes and the corresponding fluxes in the different oxygenation conditions was found only in the respirative pathway.Conclusion13C-constrained MFA enabled quantitative determination of intracellular fluxes in conditions of different redox challenges without including redox cofactors in metabolite mass balances. A redistribution of fluxes was observed not only for respirative, respiro-fermentative and fermentative metabolisms, but also for cells grown with 2.8%, 1.0% and 0.5% oxygen. Although the cellular metabolism was respiro-fermentative in each of these low oxygen conditions, the actual amount of oxygen available resulted in different contributions through respirative and fermentative pathways.

Metabolic flux profiling of recombinant protein secreting Pichia pastoris growing on glucose:methanol mixtures

BackgroundThe methylotrophic yeast Pichia pastoris has emerged as one of the most promising yeast hosts for the production of heterologous proteins. Mixed feeds of methanol and a multicarbon source instead of methanol as sole carbon source have been shown to improve product productivities and alleviate metabolic burden derived from protein production. Nevertheless, systematic quantitative studies on the relationships between the central metabolism and recombinant protein production in P. pastoris are still rather limited, particularly when growing this yeast on mixed carbon sources, thus hampering future metabolic network engineering strategies for improved protein production.ResultsThe metabolic flux distribution in the central metabolism of P. pastoris growing on a mixed feed of glucose and methanol was analyzed by Metabolic Flux Analysis (MFA) using 13C-NMR-derived constraints. For this purpose, we defined new flux ratios for methanol assimilation pathways in P. pastoris cells growing on glucose:methanol mixtures. By using this experimental approach, the metabolic burden caused by the overexpression and secretion of a Rhizopus oryzae lipase (Rol) in P. pastoris was further analyzed. This protein has been previously shown to trigger the unfolded protein response in P. pastoris. A series of 13C-tracer experiments were performed on aerobic chemostat cultivations with a control and two different Rol producing strains growing at a dilution rate of 0.09 h−1 using a glucose:methanol 80:20 (w/w) mix as carbon source.The MFA performed in this study reveals a significant redistristribution of carbon fluxes in the central carbon metabolism when comparing the two recombinant strains vs the control strain, reflected in increased glycolytic, TCA cycle and NADH regeneration fluxes, as well as higher methanol dissimilation rates.ConclusionsOverall, a further 13C-based MFA development to characterise the central metabolism of methylotrophic yeasts when growing on mixed methanol:multicarbon sources has been implemented, thus providing a new tool for the investigation of the relationships between central metabolism and protein production. Specifically, the study points at a limited but significant impact of the conformational stress associated to secretion of recombinant proteins on the central metabolism, occurring even at modest production levels.

Comparative Genome-Scale Reconstruction of Gapless Metabolic Networks for Present and Ancestral Species

We introduce a novel computational approach, CoReCo, for comparative metabolic reconstruction and provide genome-scale metabolic network models for 49 important fungal species. Leveraging on the exponential growth in sequenced genome availability, our method reconstructs genome-scale gapless metabolic networks simultaneously for a large number of species by integrating sequence data in a probabilistic framework. High reconstruction accuracy is demonstrated by comparisons to the well-curated Saccharomyces cerevisiae consensus model and large-scale knock-out experiments. Our comparative approach is particularly useful in scenarios where the quality of available sequence data is lacking, and when reconstructing evolutionary distant species. Moreover, the reconstructed networks are fully carbon mapped, allowing their use in 13C flux analysis. We demonstrate the functionality and usability of the reconstructed fungal models with computational steady-state biomass production experiment, as these fungi include some of the most important production organisms in industrial biotechnology. In contrast to many existing reconstruction techniques, only minimal manual effort is required before the reconstructed models are usable in flux balance experiments. CoReCo is available at http://esaskar.github.io/CoReCo/.

Correlation of gene expression and protein production rate - a system wide study

BackgroundGrowth rate is a major determinant of intracellular function. However its effects can only be properly dissected with technically demanding chemostat cultivations in which it can be controlled. Recent work on Saccharomyces cerevisiae chemostat cultivations provided the first analysis on genome wide effects of growth rate. In this work we study the filamentous fungus Trichoderma reesei (Hypocrea jecorina) that is an industrial protein production host known for its exceptional protein secretion capability. Interestingly, it exhibits a low growth rate protein production phenotype.ResultsWe have used transcriptomics and proteomics to study the effect of growth rate and cell density on protein production in chemostat cultivations of T. reesei. Use of chemostat allowed control of growth rate and exact estimation of the extracellular specific protein production rate (SPPR). We find that major biosynthetic activities are all negatively correlated with SPPR. We also find that expression of many genes of secreted proteins and secondary metabolism, as well as various lineage specific, mostly unknown genes are positively correlated with SPPR. Finally, we enumerate possible regulators and regulatory mechanisms, arising from the data, for this response.ConclusionsBased on these results it appears that in low growth rate protein production energy is very efficiently used primarly for protein production. Also, we propose that flux through early glycolysis or the TCA cycle is a more fundamental determining factor than growth rate for low growth rate protein production and we propose a novel eukaryotic response to this i.e. the lineage specific response (LSR).

An analytic and systematic framework for estimating metabolic flux ratios from 13C tracer experiments

BackgroundMetabolic fluxes provide invaluable insight on the integrated response of a cell to environmental stimuli or genetic modifications. Current computational methods for estimating the metabolic fluxes from 13C isotopomer measurement data rely either on manual derivation of analytic equations constraining the fluxes or on the numerical solution of a highly nonlinear system of isotopomer balance equations. In the first approach, analytic equations have to be tediously derived for each organism, substrate or labelling pattern, while in the second approach, the global nature of an optimum solution is difficult to prove and comprehensive measurements of external fluxes to augment the 13C isotopomer data are typically needed.ResultsWe present a novel analytic framework for estimating metabolic flux ratios in the cell from 13C isotopomer measurement data. In the presented framework, equation systems constraining the fluxes are derived automatically from the model of the metabolism of an organism. The framework is designed to be applicable with all metabolic network topologies, 13C isotopomer measurement techniques, substrates and substrate labelling patterns.By analyzing nuclear magnetic resonance (NMR) and mass spectrometry (MS) measurement data obtained from the experiments on glucose with the model micro-organisms Bacillus subtilis and Saccharomyces cerevisiae we show that our framework is able to automatically produce the flux ratios discovered so far by the domain experts with tedious manual analysis. Furthermore, we show by in silico calculability analysis that our framework can rapidly produce flux ratio equations – as well as predict when the flux ratios are unobtainable by linear means – also for substrates not related to glucose.ConclusionThe core of 13C metabolic flux analysis framework introduced in this article constitutes of flow and independence analysis of metabolic fragments and techniques for manipulating isotopomer measurements with vector space techniques. These methods facilitate efficient, analytic computation of the ratios between the fluxes of pathways that converge to a common junction metabolite. The framework can been seen as a generalization and formalization of existing tradition for computing metabolic flux ratios where equations constraining flux ratios are manually derived, usually without explicitly showing the formal proofs of the validity of the equations.

L-galactonate dehydratase is part of the fungal path for D-galacturonic acid catabolism.

An l‐galactonate dehydratase and the corresponding gene were identified from the mould Hypocrea jecorina (Trichoderma reesei). This novel enzyme converts l‐galactonate to l‐threo‐3‐deoxy‐hexulosonate (2‐keto‐3‐deoxy‐l‐galactonate). The enzyme is part of the fungal pathway for d‐galacturonic acid catabolism, a pathway which is only partly known. It is the second enzyme of this pathway after the d‐galacturonic acid reductase. l‐galactonate dehydratase activity is present in H. jecorina cells grown on d‐galacturonic acid but absent when other carbon sources are used for growth. A deletion of the l‐galactonate dehydratase gene in H. jecorina results in a strain with no growth on d‐galacturonic acid. The active enzyme was produced in the heterologous host Saccharomyces cerevisiae and characterized. It exhibited activity with l‐galactonate and d‐arabonate where the hydroxyl group of the C2 is in l‐ and the hydroxyl group of the C3 is in d‐configuration in the Fischer projection. However, it did not exhibit activity with d‐galactonate, d‐gluconate, l‐gulonate or d‐xylonate where the hydroxyl groups of the C2 and C3 are in different configuration.

Transcriptional responses of Saccharomyces cerevisiae to shift from respiratory and respirofermentative to fully fermentative metabolism.

In industrial fermentations of Saccharomyces cerevisiae, transient changes in oxygen concentration commonly occur and it is important to understand the behavior of cells during these changes. Glucose-limited chemostat cultivations were used to study the time-dependent effect of sudden oxygen depletion on the transcriptome of S. cerevisiae cells initially in fully aerobic or oxygen-limited conditions. The overall responses to anaerobic conditions of cells initially in different conditions were very similar. Independent of initial culture conditions, transient downregulation of genes related to growth and cell proliferation, mitochondrial translation and protein import, and sulphate assimilation was seen. In addition, transient or permanent upregulation of genes related to protein degradation, and phosphate and amino acid uptake was observed in all cultures. However, only in the initially oxygen-limited cultures was a transient upregulation of genes related to fatty acid oxidation, peroxisomal biogenesis, oxidative phosphorylation, TCA cycle, response to oxidative stress, and pentose phosphate pathway observed. Furthermore, from the initially oxygen-limited conditions, a rapid response around the metabolites of upper glycolysis and the pentose phosphate pathway was seen, while from the initially fully aerobic conditions, a slower response around the pathways for utilization of respiratory carbon sources was observed.

Dynamic flux balance analysis of the metabolism of Saccharomyces cerevisiae during the shift from fully respirative or respirofermentative metabolic states to anaerobiosis

Dynamic flux balance analysis was utilized to simulate the metabolic behaviour of initially fully respirative and respirofermentative steady‐state cultures of Saccharomyces cerevisiae during sudden oxygen depletion. The hybrid model for the dynamic flux balance analysis included a stoichiometric genome‐scale metabolic model as a static part and dynamic equations for the uptake of glucose and the cessation of respirative metabolism. The yeast consensus genome‐scale metabolic model [Herrgård MJ et al. (2008) Nat Biotechnol26, 1155–1160; Dobson PD et al. (2010) BMC Syst Biol4, 145] was refined with respect to oxygen‐dependent energy metabolism and further modified to reflect S. cerevisiae anabolism in the absence of oxygen. Dynamic flux balance analysis captured well the essential features of the dynamic metabolic behaviour of S. cerevisiae during adaptation to anaerobiosis. Modelling and simulation enabled the identification of short time‐scale flux distribution dynamics under the transition to anaerobic metabolism, during which the specific growth rate was reduced, as well as longer time‐scale process dynamics when the specific growth rate recovered. Expression of the metabolic genes was set into the context of the identified dynamics. Metabolic gene expression responses associated with the specific growth rate and with the cessation of respirative metabolism were distinguished.

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.