Generation and analysis of context-specific genome-scale metabolic models derived from single-cell RNA-Seq data.

Single-cell RNA sequencing combined with genome-scale metabolic models (GEMs) has the potential to unravel the differences in metabolism across both cell types and cell states but requires new computational methods. Here, we present a method for generating cell-type-specific genome-scale models from clusters of single-cell RNA-Seq profiles. Specifically, we developed a method to estimate the minimum number of cells required to pool to obtain stable models, a bootstrapping strategy for estimating statistical inference, and a faster version of the task-driven integrative network inference for tissues algorithm for generating context-specific GEMs. In addition, we evaluated the effect of different RNA-Seq normalization methods on model topology and differences in models generated from single-cell and bulk RNA-Seq data. We applied our methods on data from mouse cortex neurons and cells from the tumor microenvironment of lung cancer and in both cases found that almost every cell subtype had a unique metabolic profile. In addition, our approach was able to detect cancer-associated metabolic differences between cancer cells and healthy cells, showcasing its utility. We also contextualized models from 202 single-cell clusters across 19 human organs using data from Human Protein Atlas and made these available in the web portal Metabolic Atlas, thereby providing a valuable resource to the scientific community. With the ever-increasing availability of single-cell RNA-Seq datasets and continuously improved GEMs, their combination holds promise to become an important approach in the study of human metabolism.

HGTphyloDetect: facilitating the identification and phylogenetic analysis of horizontal gene transfer.

BACKGROUND: Horizontal gene transfer (HGT) is an important driver in genome evolution, gain-of-function, and metabolic adaptation to environmental niches. Genome-wide identification of putative HGT events has become increasingly practical, given the rapid growth of genomic data. However, existing HGT analysis toolboxes are not widely used, limited by their inability to perform phylogenetic reconstruction to explore potential donors, and the detection of HGT from both evolutionarily distant and closely related species. RESULTS: In this study, we have developed HGTphyloDetect, which is a versatile computational toolbox that combines high-throughput analysis with phylogenetic inference, to facilitate comprehensive investigation of HGT events. Two case studies with Saccharomyces cerevisiae and Candida versatilis demonstrate the ability of HGTphyloDetect to identify horizontally acquired genes with high accuracy. In addition, HGTphyloDetect enables phylogenetic analysis to illustrate a likely path of gene transmission among the evolutionarily distant or closely related species. CONCLUSIONS: The HGTphyloDetect computational toolbox is designed for ease of use and can accurately find HGT events with a very low false discovery rate in a high-throughput manner. The HGTphyloDetect toolbox and its related user tutorial are freely available at https://github.com/SysBioChalmers/HGTphyloDetect.

Engineering of Saccharomyces cerevisiae for enhanced metabolic robustness and L-lactic acid production from lignocellulosic biomass.

Metabolic engineering for high productivity and increased robustness is needed to enable sustainable biomanufacturing of lactic acid from lignocellulosic biomass. Lactic acid is an important commodity chemical used for instance as a monomer for production of polylactic acid, a biodegradable polymer. Here, rational and model-based optimization was used to engineer a diploid, xylose fermenting Saccharomyces cerevisiae strain to produce L-lactic acid. The metabolic flux was steered towards lactic acid through the introduction of multiple lactate dehydrogenase encoding genes while deleting ERF2, GPD1, and CYB2. A production of 93 g/L of lactic acid with a yield of 0.84 g/g was achieved using xylose as the carbon source. To increase xylose utilization and reduce acetic acid synthesis, PHO13 and ALD6 were also deleted from the strain. Finally, CDC19 encoding a pyruvate kinase was overexpressed, resulting in a yield of 0.75 g lactic acid/g sugars consumed, when the substrate used was a synthetic lignocellulosic hydrolysate medium, containing hexoses, pentoses and inhibitors such as acetate and furfural. Notably, modeling also provided leads for understanding the influence of oxygen in lactic acid production. High lactic acid production from xylose, at oxygen-limitation could be explained by a reduced flux through the oxidative phosphorylation pathway. On the contrast, higher oxygen levels were beneficial for lactic acid production with the synthetic hydrolysate medium, likely as higher ATP concentrations are needed for tolerating the inhibitors therein. The work highlights the potential of S. cerevisiae for industrial production of lactic acid from lignocellulosic biomass.

Genome-scale metabolic modeling reveals metabolic trade-offs associated with lipid production in Rhodotorula toruloides.

Rhodotorula toruloides is a non-conventional, oleaginous yeast able to naturally accumulate high amounts of microbial lipids. Constraint-based modeling of R. toruloides has been mainly focused on the comparison of experimentally measured and model predicted growth rates, while the intracellular flux patterns have been analyzed on a rather general level. Hence, the intrinsic metabolic properties of R. toruloides that make lipid synthesis possible are not thoroughly understood. At the same time, the lack of diverse physiological data sets has often been the bottleneck to predict accurate fluxes. In this study, we collected detailed physiology data sets of R. toruloides while growing on glucose, xylose, and acetate as the sole carbon source in chemically defined medium. Regardless of the carbon source, the growth was divided into two phases from which proteomic and lipidomic data were collected. Complemental physiological parameters were collected in these two phases and altogether implemented into metabolic models. Simulated intracellular flux patterns demonstrated the role of phosphoketolase in the generation of acetyl-CoA, one of the main precursors during lipid biosynthesis, while the role of ATP citrate lyase was not confirmed. Metabolic modeling on xylose as a carbon substrate was greatly improved by the detection of chirality of D-arabinitol, which together with D-ribulose were involved in an alternative xylose assimilation pathway. Further, flux patterns pointed to metabolic trade-offs associated with NADPH allocation between nitrogen assimilation and lipid biosynthetic pathways, which was linked to large-scale differences in protein and lipid content. This work includes the first extensive multi-condition analysis of R. toruloides using enzyme-constrained models and quantitative proteomics. Further, more precise kcat values should extend the application of the newly developed enzyme-constrained models that are publicly available for future studies.

Adaptations in metabolism and protein translation give rise to the Crabtree effect in yeast.

Aerobic fermentation, also referred to as the Crabtree effect in yeast, is a well-studied phenomenon that allows many eukaryal cells to attain higher growth rates at high glucose availability. Not all yeasts exhibit the Crabtree effect, and it is not known why Crabtree-negative yeasts can grow at rates comparable to Crabtree-positive yeasts. Here, we quantitatively compared two Crabtree-positive yeasts, Saccharomyces cerevisiae and Schizosaccharomyces pombe, and two Crabtree-negative yeasts, Kluyveromyces marxianus and Scheffersomyces stipitis, cultivated under glucose excess conditions. Combining physiological and proteome quantification with genome-scale metabolic modeling, we found that the two groups differ in energy metabolism and translation efficiency. In Crabtree-positive yeasts, the central carbon metabolism flux and proteome allocation favor a glucose utilization strategy minimizing proteome cost as proteins translation parameters, including ribosomal content and/or efficiency, are lower. Crabtree-negative yeasts, however, use a strategy of maximizing ATP yield, accompanied by higher protein translation parameters. Our analyses provide insight into the underlying reasons for the Crabtree effect, demonstrating a coupling to adaptations in both metabolism and protein translation.

Biosynthesis of Hesperetin, Homoeriodictyol, and Homohesperetin in a Transcriptomics-Driven Engineered Strain of Streptomyces albidoflavus.

Flavonoids exhibit various bioactivities including anti-oxidant, anti-tumor, anti-inflammatory, and anti-viral properties. Methylated flavonoids are particularly significant due to their enhanced oral bioavailability, improved intestinal absorption, and greater stability. The heterologous production of plant flavonoids in bacterial factories involves the need for enough biosynthetic precursors to allow for high production levels. These biosynthetic precursors are malonyl-CoA and l-tyrosine. In this work, to enhance flavonoid biosynthesis in Streptomyces albidoflavus, we conducted a transcriptomics study for the identification of candidate genes involved in l-tyrosine catabolism. The hypothesis was that the bacterial metabolic machinery would detect an excess of this amino acid if supplemented with the conventional culture medium and would activate the genes involved in its catabolism towards energy production. Then, by inactivating those overexpressed genes (under an excess of l-tyrosine), it would be possible to increase the intracellular pools of this precursor amino acid and eventually the final flavonoid titers in this bacterial factory. The RNAseq data analysis in the S. albidoflavus wild-type strain highlighted the hppD gene encoding 4-hydroxyphenylpyruvate dioxygenase as a promising target for knock-out, exhibiting a 23.2-fold change (FC) in expression upon l-tyrosine supplementation in comparison to control cultivation conditions. The subsequent knock-out of the hppD gene in S. albidoflavus resulted in a 1.66-fold increase in the naringenin titer, indicating enhanced flavonoid biosynthesis. Leveraging the improved strain of S. albidoflavus, we successfully synthesized the methylated flavanones hesperetin, homoeriodictyol, and homohesperetin, achieving titers of 2.52 mg/L, 1.34 mg/L, and 0.43 mg/L, respectively. In addition, the dimethoxy flavanone homohesperetin was produced as a byproduct of the endogenous metabolism of S. albidoflavus. To our knowledge, this is the first time that hppD deletion was utilized as a strategy to augment the biosynthesis of flavonoids. Furthermore, this is the first report where hesperetin and homoeriodictyol have been synthesized from l-tyrosine as a precursor. Therefore, transcriptomics is, in this case, a successful approach for the identification of catabolism reactions affecting key precursors during flavonoid biosynthesis, allowing the generation of enhanced production strains.

Improved production of Taxol® precursors in S. cerevisiae using combinatorial in silico design and metabolic engineering.

BACKGROUND: Integrated metabolic engineering approaches that combine system and synthetic biology tools enable the efficient design of microbial cell factories for synthesizing high-value products. In this study, we utilized in silico design algorithms on the yeast genome-scale model to predict genomic modifications that could enhance the production of early-step Taxol® in engineered Saccharomyces cerevisiae cells. RESULTS: Using constraint-based reconstruction and analysis (COBRA) methods, we narrowed down the solution set of genomic modification candidates. We screened 17 genomic modifications, including nine gene deletions and eight gene overexpressions, through wet-lab studies to determine their impact on taxadiene production, the first metabolite in the Taxol® biosynthetic pathway. Under different cultivation conditions, most single genomic modifications resulted in increased taxadiene production. The strain named KM32, which contained four overexpressed genes (ILV2, TRR1, ADE13, and ECM31) involved in branched-chain amino acid biosynthesis, the thioredoxin system, de novo purine synthesis, and the pantothenate pathway, respectively, exhibited the best performance. KM32 achieved a 50% increase in taxadiene production, reaching 215 mg/L. Furthermore, KM32 produced the highest reported yields of taxa-4(20),11-dien-5α-ol (T5α-ol) at 43.65 mg/L and taxa-4(20),11-dien-5-α-yl acetate (T5αAc) at 26.2 mg/L among early-step Taxol® metabolites in S. cerevisiae. CONCLUSIONS: This study highlights the effectiveness of computational and integrated approaches in identifying promising genomic modifications that can enhance the performance of yeast cell factories. By employing in silico design algorithms and wet-lab screening, we successfully improved taxadiene production in engineered S. cerevisiae strains. The best-performing strain, KM32, achieved substantial increases in taxadiene as well as production of T5α-ol and T5αAc. These findings emphasize the importance of using systematic and integrated strategies to develop efficient yeast cell factories, providing potential implications for the industrial production of high-value isoprenoids like Taxol®.

Abolishing storage lipids induces protein misfolding and stress responses in Yarrowia lipolytica.

Yarrowia lipolytica naturally saves excess carbon as storage lipids. Engineering efforts allow redirecting the high precursor flux required for lipid synthesis toward added-value chemicals such as polyketides, flavonoids, and terpenoids. To redirect precursor flux from storage lipids to other products, four genes involved in triacylglycerol and sterol ester synthesis (DGA1, DGA2, LRO1, and ARE1) can be deleted. To elucidate the effect of the deletions on cell physiology and regulation, we performed chemostat cultivations under carbon and nitrogen limitations, followed by transcriptome analysis. We found that storage lipid-free cells show an enrichment of the unfolded protein response, and several biological processes related to protein refolding and degradation are enriched. Additionally, storage lipid-free cells show an altered lipid class distribution with an abundance of potentially cytotoxic free fatty acids under nitrogen limitation. Our findings not only highlight the importance of lipid metabolism on cell physiology and proteostasis, but can also aid the development of improved chassy strains of Y. lipolytica for commodity chemical production.

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection.

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome-scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence-level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.

Validated Growth Rate-Dependent Regulation of Lipid Metabolism in Yarrowia lipolytica.

Given the strong potential of Yarrowia lipolytica to produce lipids for use as renewable fuels and oleochemicals, it is important to gain in-depth understanding of the molecular mechanism underlying its lipid accumulation. As cellular growth rate affects biomass lipid content, we performed a comparative proteomic analysis of Y. lipolytica grown in nitrogen-limited chemostat cultures at different dilution rates. After confirming the correlation between growth rate and lipid accumulation, we were able to identify various cellular functions and biological mechanisms involved in oleaginousness. Inspection of significantly up- and downregulated proteins revealed nonintuitive processes associated with lipid accumulation in this yeast. This included proteins related to endoplasmic reticulum (ER) stress, ER-plasma membrane tether proteins, and arginase. Genetic engineering of selected targets validated that some genes indeed affected lipid accumulation. They were able to increase lipid content and were complementary to other genetic engineering strategies to optimize lipid yield.

Evaluating accessibility, usability and interoperability of genome-scale metabolic models for diverse yeasts species.

Metabolic network reconstructions have become an important tool for probing cellular metabolism in the field of systems biology. They are used as tools for quantitative prediction but also as scaffolds for further knowledge contextualization. The yeast Saccharomyces cerevisiae was one of the first organisms for which a genome-scale metabolic model (GEM) was reconstructed, in 2003, and since then 45 metabolic models have been developed for a wide variety of relevant yeasts species. A systematic evaluation of these models revealed that-despite this long modeling history-the sequential process of tracing model files, setting them up for basic simulation purposes and comparing them across species and even different versions, is still not a generalizable task. These findings call the yeast modeling community to comply to standard practices on model development and sharing in order to make GEMs accessible and useful for a wider public.

Systems-level approaches for understanding and engineering of the oleaginous cell factory Yarrowia lipolytica.

Concerns about climate change and the search for renewable energy sources together with the goal of attaining sustainable product manufacturing have boosted the use of microbial platforms to produce fuels and high-value chemicals. In this regard, Yarrowia lipolytica has been known as a promising yeast with potentials in diverse array of biotechnological applications such as being a host for different oleochemicals, organic acid, and recombinant protein production. Having a rapidly increasing number of molecular and genetic tools available, Y. lipolytica has been well studied amongst oleaginous yeasts and metabolic engineering has been used to explore its potentials. More recently, with the advancement in systems biotechnology and the implementation of mathematical modeling and high throughput omics data-driven approaches, in-depth understanding of cellular mechanisms of cell factories have been made possible resulting in enhanced rational strain design. In case of Y. lipolytica, these systems-level studies and the related cutting-edge technologies have recently been initiated which is expected to result in enabling the biotechnology sector to rationally engineer Y. lipolytica-based cell factories with favorable production metrics. In this regard, here, we highlight the current status of systems metabolic engineering research and assess the potential of this yeast for future cell factory design development.

The yeastGemMap: A process diagram to assist yeast systems-metabolic studies.

Visualization is a key aspect of the analysis of omics data. Although many tools can generate pathway maps for visualization of yeast metabolism, they fail in reconstructing genome-scale metabolic diagrams of compartmentalized metabolism. Here we report on the yeastGemMap, a process diagram of whole yeast metabolism created to assist data analysis in systems-metabolic studies. The map was manually reconstructed with reactions from a compartmentalized genome-scale metabolic model, based on biochemical process diagrams typically found in educational and specialized literature. The yeastGemMap consists of 3815 reactions representing 1150 genes, 2742 metabolites, and 14 compartments. Computational functions for adapting the graphical representation of the map are also reported. These functions modify the visual representation of the map to assist in three systems-metabolic tasks: illustrating reaction networks, interpreting metabolic flux data, and visualizing omics data. The versatility of the yeastGemMap and algorithms to assist visualization of systems-metabolic data was demonstrated in various tasks, including for single lethal reaction evaluation, flux balance analysis, and transcriptomic data analysis. For instance, visual interpretation of metabolic transcriptomes of thermally evolved and parental yeast strains allowed to demonstrate that evolved strains activate a preadaptation response at 30°C, which enabled thermotolerance. A quick interpretation of systems-metabolic data is promoted with yeastGemMap visualizations.

Genome-scale model of Rhodotorula toruloides metabolism.

The basidiomycete red yeast Rhodotorula toruloides is a promising platform organism for production of biooils. We present rhto-GEM, the first genome-scale model (GEM) of R. toruloides metabolism, that was largely reconstructed using RAVEN toolbox. The model includes 852 genes, 2,731 reactions, and 2,277 metabolites, while lipid metabolism is described using the SLIMEr formalism allowing direct integration of lipid class and acyl chain experimental distribution data. The simulation results confirmed that the R. toruloides model provides valid growth predictions on glucose, xylose, and glycerol, while prediction of genetic engineering targets to increase production of linolenic acid, triacylglycerols, and carotenoids identified genes-some of which have previously been engineered to successfully increase production. This renders rtho-GEM valuable for future studies to improve the production of other oleochemicals of industrial relevance including value-added fatty acids and carotenoids, in addition to facilitate system-wide omics-data analysis in R. toruloides. Expanding the portfolio of GEMs for lipid-accumulating fungi contributes to both understanding of metabolic mechanisms of the oleaginous phenotype but also uncover particularities of the lipid production machinery in R. toruloides.

RAVEN 2.0: A versatile toolbox for metabolic network reconstruction and a case study on Streptomyces coelicolor.

RAVEN is a commonly used MATLAB toolbox for genome-scale metabolic model (GEM) reconstruction, curation and constraint-based modelling and simulation. Here we present RAVEN Toolbox 2.0 with major enhancements, including: (i) de novo reconstruction of GEMs based on the MetaCyc pathway database; (ii) a redesigned KEGG-based reconstruction pipeline; (iii) convergence of reconstructions from various sources; (iv) improved performance, usability, and compatibility with the COBRA Toolbox. Capabilities of RAVEN 2.0 are here illustrated through de novo reconstruction of GEMs for the antibiotic-producing bacterium Streptomyces coelicolor. Comparison of the automated de novo reconstructions with the iMK1208 model, a previously published high-quality S. coelicolor GEM, exemplifies that RAVEN 2.0 can capture most of the manually curated model. The generated de novo reconstruction is subsequently used to curate iMK1208 resulting in Sco4, the most comprehensive GEM of S. coelicolor, with increased coverage of both primary and secondary metabolism. This increased coverage allows the use of Sco4 to predict novel genome editing targets for optimized secondary metabolites production. As such, we demonstrate that RAVEN 2.0 can be used not only for de novo GEM reconstruction, but also for curating existing models based on up-to-date databases. Both RAVEN 2.0 and Sco4 are distributed through GitHub to facilitate usage and further development by the community (https://github.com/SysBioChalmers/RAVEN and https://github.com/SysBioChalmers/Streptomyces_coelicolor-GEM).

Improving the phenotype predictions of a yeast genome-scale metabolic model by incorporating enzymatic constraints.

Genome-scale metabolic models (GEMs) are widely used to calculate metabolic phenotypes. They rely on defining a set of constraints, the most common of which is that the production of metabolites and/or growth are limited by the carbon source uptake rate. However, enzyme abundances and kinetics, which act as limitations on metabolic fluxes, are not taken into account. Here, we present GECKO, a method that enhances a GEM to account for enzymes as part of reactions, thereby ensuring that each metabolic flux does not exceed its maximum capacity, equal to the product of the enzyme's abundance and turnover number. We applied GECKO to a Saccharomyces cerevisiae GEM and demonstrated that the new model could correctly describe phenotypes that the previous model could not, particularly under high enzymatic pressure conditions, such as yeast growing on different carbon sources in excess, coping with stress, or overexpressing a specific pathway. GECKO also allows to directly integrate quantitative proteomics data; by doing so, we significantly reduced flux variability of the model, in over 60% of metabolic reactions. Additionally, the model gives insight into the distribution of enzyme usage between and within metabolic pathways. The developed method and model are expected to increase the use of model-based design in metabolic engineering.

Advancing metabolic engineering of Yarrowia lipolytica using the CRISPR/Cas system.

The oleaginous yeast Yarrowia lipolytica is widely used for the production of both bulk and fine chemicals, including organic acids, fatty acid-derived biofuels and chemicals, polyunsaturated fatty acids, single-cell proteins, terpenoids, and other valuable products. Consequently, it is becoming increasingly popular for metabolic engineering applications. Multiple gene manipulation tools including URA blast, Cre/LoxP, and transcription activator-like effector nucleases (TALENs) have been developed for metabolic engineering in Y. lipolytica. However, the low efficiency and time-consuming procedures involved in these methods hamper further research. The emergence of the CRISPR/Cas system offers a potential solution for these problems due to its high efficiency, ease of operation, and time savings, which can significantly accelerate the genomic engineering of Y. lipolytica. In this review, we summarize the research progress on the development of CRISPR/Cas systems for Y. lipolytica, including Cas9 proteins and sgRNA expression strategies, as well as gene knock-out/knock-in and repression/activation applications. Finally, the most promising and tantalizing future prospects in this area are highlighted.

Probing the metabolic network in bloodstream-form Trypanosoma brucei using untargeted metabolomics with stable isotope labelled glucose.

Metabolomics coupled with heavy-atom isotope-labelled glucose has been used to probe the metabolic pathways active in cultured bloodstream form trypomastigotes of Trypanosoma brucei, a parasite responsible for human African trypanosomiasis. Glucose enters many branches of metabolism beyond glycolysis, which has been widely held to be the sole route of glucose metabolism. Whilst pyruvate is the major end-product of glucose catabolism, its transamination product, alanine, is also produced in significant quantities. The oxidative branch of the pentose phosphate pathway is operative, although the non-oxidative branch is not. Ribose 5-phosphate generated through this pathway distributes widely into nucleotide synthesis and other branches of metabolism. Acetate, derived from glucose, is found associated with a range of acetylated amino acids and, to a lesser extent, fatty acids; while labelled glycerol is found in many glycerophospholipids. Glucose also enters inositol and several sugar nucleotides that serve as precursors to macromolecule biosynthesis. Although a Krebs cycle is not operative, malate, fumarate and succinate, primarily labelled in three carbons, were present, indicating an origin from phosphoenolpyruvate via oxaloacetate. Interestingly, the enzyme responsible for conversion of phosphoenolpyruvate to oxaloacetate, phosphoenolpyruvate carboxykinase, was shown to be essential to the bloodstream form trypanosomes, as demonstrated by the lethal phenotype induced by RNAi-mediated downregulation of its expression. In addition, glucose derivatives enter pyrimidine biosynthesis via oxaloacetate as a precursor to aspartate and orotate.

TrypanoCyc: a community-led biochemical pathways database for Trypanosoma brucei.

The metabolic network of a cell represents the catabolic and anabolic reactions that interconvert small molecules (metabolites) through the activity of enzymes, transporters and non-catalyzed chemical reactions. Our understanding of individual metabolic networks is increasing as we learn more about the enzymes that are active in particular cells under particular conditions and as technologies advance to allow detailed measurements of the cellular metabolome. Metabolic network databases are of increasing importance in allowing us to contextualise data sets emerging from transcriptomic, proteomic and metabolomic experiments. Here we present a dynamic database, TrypanoCyc (http://www.metexplore.fr/trypanocyc/), which describes the generic and condition-specific metabolic network of Trypanosoma brucei, a parasitic protozoan responsible for human and animal African trypanosomiasis. In addition to enabling navigation through the BioCyc-based TrypanoCyc interface, we have also implemented a network-based representation of the information through MetExplore, yielding a novel environment in which to visualise the metabolism of this important parasite.

Genome-scale metabolic model of Pichia pastoris with native and humanized glycosylation of recombinant proteins.

Pichia pastoris is used for commercial production of human therapeutic proteins, and genome-scale models of P. pastoris metabolism have been generated in the past to study the metabolism and associated protein production by this yeast. A major challenge with clinical usage of recombinant proteins produced by P. pastoris is the difference in N-glycosylation of proteins produced by humans and this yeast. However, through metabolic engineering, a P. pastoris strain capable of producing humanized N-glycosylated proteins was constructed. The current genome-scale models of P. pastoris do not address native nor humanized N-glycosylation, and we therefore developed ihGlycopastoris, an extension to the iLC915 model with both native and humanized N-glycosylation for recombinant protein production, but also an estimation of N-glycosylation of P. pastoris native proteins. This new model gives a better prediction of protein yield, demonstrates the effect of the different types of N-glycosylation of protein yield, and can be used to predict potential targets for strain improvement. The model represents a step towards a more complete description of protein production in P. pastoris, which is required for using these models to understand and optimize protein production processes.