Yeast 5 - an expanded reconstruction of the Saccharomyces cerevisiae metabolic network

Network thermodynamic curation of human and yeast genome-scale metabolic models

Genome-scale models are used for an ever-widening range of applications. Although there has been much focus on specifying the stoichiometric matrix, the predictive power of genome-scale models equally depends on reaction directions. Two-thirds of reactions in the two eukaryotic reconstructions Homo sapiens Recon 1 and Yeast 5 are specified as irreversible. However, these specifications are mainly based on biochemical textbooks or on their similarity to other organisms and are rarely underpinned by detailed thermodynamic analysis. In this study, a to our knowledge new workflow combining network-embedded thermodynamic and flux variability analysis was used to evaluate existing irreversibility constraints in Recon 1 and Yeast 5 and to identify new ones. A total of 27 and 16 new irreversible reactions were identified in Recon 1 and Yeast 5, respectively, whereas only four reactions were found with directions incorrectly specified against thermodynamics (three in Yeast 5 and one in Recon 1). The workflow further identified for both models several isolated internal loops that require further curation. The framework also highlighted the need for substrate channeling (in human) and ATP hydrolysis (in yeast) for the essential reaction catalyzed by phosphoribosylaminoimidazole carboxylase in purine metabolism. Finally, the framework highlighted differences in proline metabolism between yeast (cytosolic anabolism and mitochondrial catabolism) and humans (exclusively mitochondrial metabolism). We conclude that network-embedded thermodynamics facilitates the specification and validation of irreversibility constraints in compartmentalized metabolic models, at the same time providing further insight into network properties.

Cancer cell metabolism as new targets for novel designed therapies

Metabolic processes are altered in cancer cells, which obtain advantages from this metabolic reprogramming in terms of energy production and synthesis of biomolecules that sustain their uncontrolled proliferation. Due to the conceptual progresses in the last decade, metabolic reprogramming was recently included as one of the new hallmarks of cancer. The advent of high-throughput technologies to amass an abundance of omic data, together with the development of new computational methods that allow the integration and analysis of omic data by using genome-scale reconstructions of human metabolism, have increased and accelerated the discovery and development of anticancer drugs and tumor-specific metabolic biomarkers. Here we review and discuss the latest advances in the context of metabolic reprogramming and the future in cancer research.

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.

The LASER database: Formalizing design rules for metabolic engineering

The ability of metabolic engineers to conceptualize, implement, and evaluate strain designs has dramatically increased in the last decade. Unlike other engineering fields, no centralized, open-access, and easily searched repository exists for cataloging these designs and the lessons learned from their construction and evaluation. To address this issue, we have developed a repository for metabolic engineering strain designs, known as LASER (Learning Assisted Strain EngineeRing, laser.colorado.edu) and a formal standard for disseminating designs to metabolic engineers. Curation of every available genetically-defined E. coli and S. cerevisiae strain from 310 metabolic engineering papers published over the last 21 years yields a total of 417 designs containing a total of 2661 genetic modifications. This collection has been deposited in LASER and represents the known bibliome of genetically defined and tested metabolic engineering designs in the academic literature. Properties of LASER designs and the analysis pipeline are examined to provide insight into LASER capabilities. Several future research directions utilizing LASER capabilities are discussed to highlight the potential of the LASER database for metabolic engineering.

Sequence-based Network Completion Reveals the Integrality of Missing Reactions in Metabolic Networks

Genome-scale metabolic models are central in connecting genotypes to metabolic phenotypes. However, even for well studied organisms, such as Escherichia coli, draft networks do not contain a complete biochemical network. Missing reactions are referred to as gaps. These gaps need to be filled to enable functional analysis, and gap-filling choices influence model predictions. To investigate whether functional networks existed where all gap-filling reactions were supported by sequence similarity to annotated enzymes, four draft networks were supplemented with all reactions from the Model SEED database for which minimal sequence similarity was found in their genomes. Quadratic programming revealed that the number of reactions that could partake in a gap-filling solution was vast: 3,270 in the case of E. coli, where 72% of the metabolites in the draft network could connect a gap-filling solution. Nonetheless, no network could be completed without the inclusion of orphaned enzymes, suggesting that parts of the biochemistry integral to biomass precursor formation are uncharacterized. However, many gap-filling reactions were well determined, and the resulting networks showed improved prediction of gene essentiality compared with networks generated through canonical gap filling. In addition, gene essentiality predictions that were sensitive to poorly determined gap-filling reactions were of poor quality, suggesting that damage to the network structure resulting from the inclusion of erroneous gap-filling reactions may be predictable.

RobOKoD: microbial strain design for (over)production of target compounds

Sustainable production of target compounds such as biofuels and high-value chemicals for pharmaceutical, agrochemical, and chemical industries is becoming an increasing priority given their current dependency upon diminishing petrochemical resources. Designing these strains is difficult, with current methods focusing primarily on knocking-out genes, dismissing other vital steps of strain design including the overexpression and dampening of genes. The design predictions from current methods also do not translate well-into successful strains in the laboratory. Here, we introduce RobOKoD (Robust, Overexpression, Knockout and Dampening), a method for predicting strain designs for overproduction of targets. The method uses flux variability analysis to profile each reaction within the system under differing production percentages of target-compound and biomass. Using these profiles, reactions are identified as potential knockout, overexpression, or dampening targets. The identified reactions are ranked according to their suitability, providing flexibility in strain design for users. The software was tested by designing a butanol-producing Escherichia coli strain, and was compared against the popular OptKnock and RobustKnock methods. RobOKoD shows favorable design predictions, when predictions from these methods are compared to a successful butanol-producing experimentally-validated strain. Overall RobOKoD provides users with rankings of predicted beneficial genetic interventions with which to support optimized strain design.

On correlated reaction sets and coupled reaction sets in metabolic networks

Two reactions are in the same "correlated reaction set" (or "Co-Set") if their fluxes are linearly correlated. On the other hand, two reactions are "coupled" if nonzero flux through one reaction implies nonzero flux through the other reaction. Flux correlation analysis has been previously used in the analysis of enzyme dysregulation and enzymopathy, while flux coupling analysis has been used to predict co-expression of genes and to model network evolution. The goal of this paper is to emphasize, through a few examples, that these two concepts are inherently different. In other words, except for the case of full coupling, which implies perfect correlation between two fluxes (R(2) = 1), there are no constraints on Pearson correlation coefficients (CC) in case of any other type of (un)coupling relations. In other words, Pearson CC can take any value between 0 and 1 in other cases. Furthermore, by analyzing genome-scale metabolic networks, we confirm that there are some examples in real networks of bacteria, yeast and human, which approve that flux coupling and flux correlation cannot be used interchangeably.

Applications of genome-scale metabolic network model in metabolic engineering

Genome-scale metabolic network model (GEM) is a fundamental framework in systems metabolic engineering. GEM is built upon extensive experimental data and literature information on gene annotation and function, metabolites and enzymes so that it contains all known metabolic reactions within an organism. Constraint-based analysis of GEM enables the identification of phenotypic properties of an organism and hypothesis-driven engineering of cellular functions to achieve objectives. Along with the advances in omics, high-throughput technology and computational algorithms, the scope and applications of GEM have substantially expanded. In particular, various computational algorithms have been developed to predict beneficial gene deletion and amplification targets and used to guide the strain development process for the efficient production of industrially important chemicals. Furthermore, an Escherichia coli GEM was integrated with a pathway prediction algorithm and used to evaluate all possible routes for the production of a list of commodity chemicals in E. coli. Combined with the wealth of experimental data produced by high-throughput techniques, much effort has been exerted to add more biological contexts into GEM through the integration of omics data and regulatory network information for the mechanistic understanding and improved prediction capabilities. In this paper, we review the recent developments and applications of GEM focusing on the GEM-based computational algorithms available for microbial metabolic engineering.

Critical assessment of genome-scale metabolic networks: the need for a unified standard

Genome-scale metabolic networks have been reconstructed for several organisms. These metabolic networks provide detailed information about the metabolism inside the cells, coupled with the genomic, proteomic and thermodynamic information. These networks are widely simulated using 'constraint-based' modelling techniques and find applications ranging from strain improvement for metabolic engineering to prediction of drug targets in pathogenic organisms. Components of these metabolic networks are represented in multiple file formats and also using different markup languages, with varying levels of annotations; this leads to inconsistencies and increases the complexities in comparing and analysing reconstructions on multiple platforms. In this work, we critically examine nearly 100 published genome-scale metabolic networks and their corresponding constraint-based models and discuss various issues with respect to model quality. One of the major concerns is the lack of annotations using standard identifiers that can uniquely describe several components such as metabolites, genes, proteins and reactions. We also find that many models do not have complete information regarding constraints on reactions fluxes and objective functions for carrying out simulations. Overall, our analysis highlights the need for a widely acceptable standard for representing constraint-based models. A rigorous standard can help in streamlining the process of reconstruction and improve the quality of reconstructed metabolic models.

BioPreDyn-bench: a suite of benchmark problems for dynamic modelling in systems biology

Background

Dynamic modelling is one of the cornerstones of systems biology. Many research efforts are currently being invested in the development and exploitation of large-scale kinetic models. The associated problems of parameter estimation (model calibration) and optimal experimental design are particularly challenging. The community has already developed many methods and software packages which aim to facilitate these tasks. However, there is a lack of suitable benchmark problems which allow a fair and systematic evaluation and comparison of these contributions.

Results

Here we present BioPreDyn-bench, a set of challenging parameter estimation problems which aspire to serve as reference test cases in this area. This set comprises six problems including medium and large-scale kinetic models of the bacterium E. coli, baker’s yeast S. cerevisiae, the vinegar fly D. melanogaster, Chinese Hamster Ovary cells, and a generic signal transduction network. The level of description includes metabolism, transcription, signal transduction, and development. For each problem we provide (i) a basic description and formulation, (ii) implementations ready-to-run in several formats, (iii) computational results obtained with specific solvers, (iv) a basic analysis and interpretation.

Conclusions

This suite of benchmark problems can be readily used to evaluate and compare parameter estimation methods. Further, it can also be used to build test problems for sensitivity and identifiability analysis, model reduction and optimal experimental design methods. The suite, including codes and documentation, can be freely downloaded from the BioPreDyn-bench website, https://sites.google.com/site/biopredynbenchmarks/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0144-4) contains supplementary material, which is available to authorized users.

Computational prediction of essential metabolic genes using constraint-based approaches

In this chapter, we describe the application of constraint-based modeling to predict the impact of gene deletions on a metabolic phenotype. The metabolic reactions taking place inside cells form large networks, which have been reconstructed at a genome-scale for several organisms at increasing levels of detail. By integrating mathematical modeling techniques with biochemical principles, constraint-based approaches enable predictions of metabolite fluxes and growth under specific environmental conditions or for genetically modified microorganisms. Similar to the experimental knockout of a gene, predicting the essentiality of a metabolic gene for a phenotype further allows to generate hypotheses on its biological function and design of genetic engineering strategies for biotechnological applications. Here, we summarize the principles of constraint-based approaches and provide a detailed description of the procedure to predict the essentiality of metabolic genes with respect to a specific metabolic function. We exemplify the approach by predicting the essentiality of reactions in the citric acid cycle for the production of glucose from fatty acids.

MetDFBA: incorporating time-resolved metabolomics measurements into dynamic flux balance analysis

This paper presents MetDFBA, a new approach incorporating experimental metabolomics time-series into constraint-based modeling. The method can be used for hypothesis testing and predicting dynamic flux profiles.

Model-guided identification of gene deletion targets for metabolic engineering in Saccharomyces cerevisiae

Identification of metabolic engineering strategies for rerouting intracellular fluxes towards a desired product is often a challenging task owing to the topological and regulatory complexity of metabolic networks. Genome-scale metabolic models help tackling this complexity through systematic consideration of mass balance and reaction directionality constraints over the entire network. Here, we describe how genome-scale metabolic models can be used for identifying gene deletion targets leading to increased production of the desired product. Vanillin production in Saccharomyces cerevisiae is used as a case study throughout this chapter.

Systems biology of monovalent cation homeostasis in yeast: the translucent contribution

Maintenance of monovalent cation homeostasis (mainly K(+) and Na(+)) is vital for cell survival, and cation toxicity is at the basis of a myriad of relevant phenomena, such as salt stress in crops and diverse human diseases. Full understanding of the importance of monovalent cations in the biology of the cell can only be achieved from a systemic perspective. Translucent is a multinational project developed within the context of the SysMO (System Biology of Microorganisms) initiative and focussed in the study of cation homeostasis using the well-known yeast Saccharomyces cerevisiae as a model. The present review summarize how the combination of biochemical, genetic, genomic and computational approaches has boosted our knowledge in this field, providing the basis for a more comprehensive and coherent vision of the role of monovalent cations in the biology of the cell.

Pathway connectivity and signaling coordination in the yeast stress-activated signaling network

Stressed cells coordinate a multi-faceted response spanning many levels of physiology. Yet knowledge of the complete stress-activated regulatory network as well as design principles for signal integration remains incomplete. We developed an experimental and computational approach to integrate available protein interaction data with gene fitness contributions, mutant transcriptome profiles, and phospho-proteome changes in cells responding to salt stress, to infer the salt-responsive signaling network in yeast. The inferred subnetwork presented many novel predictions by implicating new regulators, uncovering unrecognized crosstalk between known pathways, and pointing to previously unknown ‘hubs’ of signal integration. We exploited these predictions to show that Cdc14 phosphatase is a central hub in the network and that modification of RNA polymerase II coordinates induction of stress-defense genes with reduction of growth-related transcripts. We find that the orthologous human network is enriched for cancer-causing genes, underscoring the importance of the subnetwork's predictions in understanding stress biology.

Prediction of metabolic flux distribution from gene expression data based on the flux minimization principle

Prediction of possible flux distributions in a metabolic network provides detailed phenotypic information that links metabolism to cellular physiology. To estimate metabolic steady-state fluxes, the most common approach is to solve a set of macroscopic mass balance equations subjected to stoichiometric constraints while attempting to optimize an assumed optimal objective function. This assumption is justifiable in specific cases but may be invalid when tested across different conditions, cell populations, or other organisms. With an aim to providing a more consistent and reliable prediction of flux distributions over a wide range of conditions, in this article we propose a framework that uses the flux minimization principle to predict active metabolic pathways from mRNA expression data. The proposed algorithm minimizes a weighted sum of flux magnitudes, while biomass production can be bounded to fit an ample range from very low to very high values according to the analyzed context. We have formulated the flux weights as a function of the corresponding enzyme reaction's gene expression value, enabling the creation of context-specific fluxes based on a generic metabolic network. In case studies of wild-type Saccharomyces cerevisiae, and wild-type and mutant Escherichia coli strains, our method achieved high prediction accuracy, as gauged by correlation coefficients and sums of squared error, with respect to the experimentally measured values. In contrast to other approaches, our method was able to provide quantitative predictions for both model organisms under a variety of conditions. Our approach requires no prior knowledge or assumption of a context-specific metabolic functionality and does not require trial-and-error parameter adjustments. Thus, our framework is of general applicability for modeling the transcription-dependent metabolism of bacteria and yeasts.

Computational analysis of reciprocal association of metabolism and epigenetics in the budding yeast: a genome-scale metabolic model (GSMM) approach

Metaboloepigenetics is a newly coined term in biological sciences that investigates the crosstalk between epigenetic modifications and metabolism. The reciprocal relation between biochemical transformations and gene expression regulation has been experimentally demonstrated in cancers and metabolic syndromes. In this study, we explored the metabolism-histone modifications crosstalk by topological analysis and constraint-based modeling approaches in the budding yeast. We constructed nine models through the integration of gene expression data of four mutated histone tails into a genome-scale metabolic model of yeast. Accordingly, we defined the centrality indices of the lowly expressed enzymes in the undirected enzyme-centric network of yeast by CytoHubba plug-in in Cytoscape. To determine the global effects of histone modifications on the yeast metabolism, the growth rate and the range of possible flux values of reactions, we used constraint-based modeling approach. Centrality analysis shows that the lowly expressed enzymes could affect and control the yeast metabolic network. Besides, constraint-based modeling results are in a good agreement with the experimental findings, confirming that the mutations in histone tails lead to non-lethal alterations in the yeast, but have diverse effects on the growth rate and reveal the functional redundancy.

Assessing the metabolic impact of nitrogen availability using a compartmentalized maize leaf genome-scale model

Maize (Zea mays) is an important C4 plant due to its widespread use as a cereal and energy crop. A second-generation genome-scale metabolic model for the maize leaf was created to capture C4 carbon fixation and investigate nitrogen (N) assimilation by modeling the interactions between the bundle sheath and mesophyll cells. The model contains gene-protein-reaction relationships, elemental and charge-balanced reactions, and incorporates experimental evidence pertaining to the biomass composition, compartmentalization, and flux constraints. Condition-specific biomass descriptions were introduced that account for amino acids, fatty acids, soluble sugars, proteins, chlorophyll, lignocellulose, and nucleic acids as experimentally measured biomass constituents. Compartmentalization of the model is based on proteomic/transcriptomic data and literature evidence. With the incorporation of information from the MetaCrop and MaizeCyc databases, this updated model spans 5,824 genes, 8,525 reactions, and 9,153 metabolites, an increase of approximately 4 times the size of the earlier iRS1563 model. Transcriptomic and proteomic data have also been used to introduce regulatory constraints in the model to simulate an N-limited condition and mutants deficient in glutamine synthetase, gln1-3 and gln1-4. Model-predicted results achieved 90% accuracy when comparing the wild type grown under an N-complete condition with the wild type grown under an N-deficient condition.

A blueprint of the amino acid biosynthesis network of hemiascomycetes

The structure and regulation of biosynthesis pathways in Saccharomyces cerevisiae have been detailed extensively. For other hemiascomycetes, genomic sequences are primarily available, whereas biochemical information on them is scarce. The resulting biochemical networks that are used for research in basic science and biotechnology are often biased by data from S. cerevisiae, assuming that there are often implicitly conserved structures between species. We examined the structure of the amino acid biosynthesis network in nine hemiascomycetes, spanning the phylogenetic clade. Differences in the genetic inventory included the presence and absence of isoenzymes and compartmentation of the pathways. Notably, no two hemiascomycetes had identical genetic inventories. For example, the lack of the mitochondrial αIPMS isoenzyme and presence of only one copy of the BCAA aminotransferase in Pichia pastoris indicate a disparately compartmented leucine biosynthesis pathway. Our findings suggest that αIPMS and BCAA aminotransferase are solely located in the cytosol of P. pastoris, requiring correction of the leucine biosynthesis pathway layout in this species. Our results argue for careful use of information from S. cerevisiae and for joint efforts to fill the knowledge gaps in other species. Such analysis will lead to contributions in biotechnology disciplines, such as protein production and compartment engineering.

The yeast deletion collection: a decade of functional genomics

The yeast deletion collections comprise >21,000 mutant strains that carry precise start-to-stop deletions of ∼6000 open reading frames. This collection includes heterozygous and homozygous diploids, and haploids of both MATa and MATα mating types. The yeast deletion collection, or yeast knockout (YKO) set, represents the first and only complete, systematically constructed deletion collection available for any organism. Conceived during the Saccharomyces cerevisiae sequencing project, work on the project began in 1998 and was completed in 2002. The YKO strains have been used in numerous laboratories in >1000 genome-wide screens. This landmark genome project has inspired development of numerous genome-wide technologies in organisms from yeast to man. Notable spinoff technologies include synthetic genetic array and HIPHOP chemogenomics. In this retrospective, we briefly describe the yeast deletion project and some of its most noteworthy biological contributions and the impact that these collections have had on the yeast research community and on genomics in general.

A framework for application of metabolic modeling in yeast to predict the effects of nsSNV in human orthologs

BACKGROUND

We have previously suggested a method for proteome wide analysis of variation at functional residues wherein we identified the set of all human genes with nonsynonymous single nucleotide variation (nsSNV) in the active site residue of the corresponding proteins. 34 of these proteins were shown to have a 1:1:1 enzyme:pathway:reaction relationship, making these proteins ideal candidates for laboratory validation through creation and observation of specific yeast active site knock-outs and downstream targeted metabolomics experiments. Here we present the next step in the workflow toward using yeast metabolic modeling to predict human metabolic behavior resulting from nsSNV.

RESULTS

For the previously identified candidate proteins, we used the reciprocal best BLAST hits method followed by manual alignment and pathway comparison to identify 6 human proteins with yeast orthologs which were suitable for flux balance analysis (FBA). 5 of these proteins are known to be associated with diseases, including ribose 5-phosphate isomerase deficiency, myopathy with lactic acidosis and sideroblastic anaemia, anemia due to disorders of glutathione metabolism, and two porphyrias, and we suspect the sixth enzyme to have disease associations which are not yet classified or understood based on the work described herein.

CONCLUSIONS

Preliminary findings using the Yeast 7.0 FBA model show lack of growth for only one enzyme, but augmentation of the Yeast 7.0 biomass function to better simulate knockout of certain genes suggested physiological relevance of variations in three additional proteins. Thus, we suggest the following four proteins for laboratory validation: delta-aminolevulinic acid dehydratase, ferrochelatase, ribose-5 phosphate isomerase and mitochondrial tyrosyl-tRNA synthetase. This study indicates that the predictive ability of this method will improve as more advanced, comprehensive models are developed. Moreover, these findings will be useful in the development of simple downstream biochemical or mass-spectrometric assays to corroborate these predictions and detect presence of certain known nsSNVs with deleterious outcomes. Results may also be useful in predicting as yet unknown outcomes of active site nsSNVs for enzymes that are not yet well classified or annotated.

BioMet Toolbox 2.0: genome-wide analysis of metabolism and omics data

Analysis of large data sets using computational and mathematical tools have become a central part of biological sciences. Large amounts of data are being generated each year from different biological research fields leading to a constant development of software and algorithms aimed to deal with the increasing creation of information. The BioMet Toolbox 2.0 integrates a number of functionalities in a user-friendly environment enabling the user to work with biological data in a web interface. The unique and distinguishing feature of the BioMet Toolbox 2.0 is to provide a web user interface to tools for metabolic pathways and omics analysis developed under different platform-dependent environments enabling easy access to these computational tools.

Inferring host gene subnetworks involved in viral replication

Systematic, genome-wide loss-of-function experiments can be used to identify host factors that directly or indirectly facilitate or inhibit the replication of a virus in a host cell. We present an approach that combines an integer linear program and a diffusion kernel method to infer the pathways through which those host factors modulate viral replication. The inputs to the method are a set of viral phenotypes observed in single-host-gene mutants and a background network consisting of a variety of host intracellular interactions. The output is an ensemble of subnetworks that provides a consistent explanation for the measured phenotypes, predicts which unassayed host factors modulate the virus, and predicts which host factors are the most direct interfaces with the virus. We infer host-virus interaction subnetworks using data from experiments screening the yeast genome for genes modulating the replication of two RNA viruses. Because a gold-standard network is unavailable, we assess the predicted subnetworks using both computational and qualitative analyses. We conduct a cross-validation experiment in which we predict whether held-aside test genes have an effect on viral replication. Our approach is able to make high-confidence predictions more accurately than several baselines, and about as well as the best baseline, which does not infer mechanistic pathways. We also examine two kinds of predictions made by our method: which host factors are nearest to a direct interaction with a viral component, and which unassayed host genes are likely to be involved in viral replication. Multiple predictions are supported by recent independent experimental data, or are components or functional partners of confirmed relevant complexes or pathways. Integer program code, background network data, and inferred host-virus subnetworks are available at http://www.biostat.wisc.edu/~craven/chasman_host_virus/.

Direct calculation of elementary flux modes satisfying several biological constraints in genome-scale metabolic networks

MOTIVATION

The concept of Elementary Flux Mode (EFM) has been widely used for the past 20 years. However, its application to genome-scale metabolic networks (GSMNs) is still under development because of methodological limitations. Therefore, novel approaches are demanded to extend the application of EFMs. A novel family of methods based on optimization is emerging that provides us with a subset of EFMs. Because the calculation of the whole set of EFMs goes beyond our capacity, performing a selective search is a proper strategy.

RESULTS

Here, we present a novel mathematical approach calculating EFMs fulfilling additional linear constraints. We validated our approach based on two metabolic networks in which all the EFMs can be obtained. Finally, we analyzed the performance of our methodology in the GSMN of the yeast Saccharomyces cerevisiae by calculating EFMs producing ethanol with a given minimum carbon yield. Overall, this new approach opens new avenues for the calculation of EFMs in GSMNs.

AVAILABILITY AND IMPLEMENTATION

Matlab code is provided in the supplementary online materials

CONTACT

fplanes@ceit.es.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Broad metabolic sensitivity profiling of a prototrophic yeast deletion collection

Background

Genome-wide sensitivity screens in yeast have been immensely popular following the construction of a collection of deletion mutants of non-essential genes. However, the auxotrophic markers in this collection preclude experiments on minimal growth medium, one of the most informative metabolic environments. Here we present quantitative growth analysis for mutants in all 4,772 non-essential genes from our prototrophic deletion collection across a large set of metabolic conditions.

Results

The complete collection was grown in environments consisting of one of four possible carbon sources paired with one of seven nitrogen sources, for a total of 28 different well-defined metabolic environments. The relative contributions to mutants' fitness of each carbon and nitrogen source were determined using multivariate statistical methods. The mutant profiling recovered known and novel genes specific to the processing of nutrients and accurately predicted functional relationships, especially for metabolic functions. A benchmark of genome-scale metabolic network modeling is also given to demonstrate the level of agreement between current in silico predictions and hitherto unavailable experimental data.

Conclusions

These data address a fundamental deficiency in our understanding of the model eukaryote Saccharomyces cerevisiae and its response to the most basic of environments. While choice of carbon source has the greatest impact on cell growth, specific effects due to nitrogen source and interactions between the nutrients are frequent. We demonstrate utility in characterizing genes of unknown function and illustrate how these data can be integrated with other whole-genome screens to interpret similarities between seemingly diverse perturbation types.

Incorporation of flexible objectives and time-linked simulation with flux balance analysis

We present two modifications of the flux balance analysis (FBA) metabolic modeling framework which relax implicit assumptions of the biomass reaction. Our flexible flux balance analysis (flexFBA) objective removes the fixed proportion between reactants, and can therefore produce a subset of biomass reactants. Our time-linked flux balance analysis (tFBA) simulation removes the fixed proportion between reactants and byproducts, and can therefore describe transitions between metabolic steady states. Used together, flexFBA and tFBA model a time scale shorter than the regulatory and growth steady state encoded by the biomass reaction. This combined short-time FBA method is intended for integrated modeling applications to enable detailed and dynamic depictions of microbial physiology such as whole-cell modeling. For example, when modeling Escherichia coli, it avoids artifacts caused by low-copy-number enzymes in single-cell models with kinetic bounds. Even outside integrated modeling contexts, the detailed predictions of flexFBA and tFBA complement existing FBA techniques. We show detailed metabolite production of in silico knockouts used to identify when correct essentiality predictions are made for the wrong reason.

Transcriptomic insights into the molecular response of Saccharomyces cerevisiae to linoleic acid hydroperoxide

Eukaryotic microorganisms are constantly challenged by reactive oxygen species derived endogenously or encountered in their environment. Such adversity is particularly applied to Saccharomyces cerevisiae under harsh industrial conditions. One of the major oxidants to challenge S. cerevisiae is linoleic acid hydroperoxide (LoaOOH). This study, which used genome-wide microarray analysis in conjunction with deletion mutant screening, uncovered the molecular pathways of S. cerevisiae that were altered by an arresting concentration of LoaOOH (75 μM). The oxidative stress response, iron homeostasis, detoxification through PDR transport and direct lipid β-oxidation were evident through the induction of the genes encoding for peroxiredoxins (GPX2, TSA2), the NADPH:oxidoreductase (OYE3), iron uptake (FIT2, ARN2, FET3), PDR transporters (PDR5, PDR15, SNQ2) and β-oxidation machinery (FAA2, POX1). Further, we discovered that Gpx3p, the dual redox sensor and peroxidase, is required for protection against LoaOOH, indicated by the sensitivity of gpx3Δ to a mild dose of LoaOOH (37.5 μM). Deletion of GPX3 conferred a greater sensitivity to LoaOOH than the loss of its signalling partner YAP1. Deletion of either of the iron homeostasis regulators AFT1 or AFT2 also resulted in sensitivity to LoaOOH. These novel findings for Gpx3p, Aft1p and Aft2p point to their distinct roles in response to the lipid peroxide. Finally, the expression of 89 previously uncharacterised genes was significantly altered against LoaOOH, which will contribute to their eventual annotation.

Exploration and comparison of inborn capacity of aerobic and anaerobic metabolisms of Saccharomyces cerevisiae for microbial electrical current production

Saccharomyces cerevisiae possesses numerous advantageous biological features, such as being robust, easily handled, mostly non-pathogenic and having high catabolic rates, etc., which can be considered as merits for being used as a promising biocatalyst in microbial fuel cells (MFCs) for electricity generation. Previous studies have developed efficient MFC configurations to convert metabolic electron shuttles, such as cytoplasmic NADH, into usable electric current. However, no studies have elucidated the maximum potential of S. cerevisiae for current output and the underlying metabolic pathways, resulting from the interaction of thousands of reactions inside the cell during MFC operation. To address these two key issues, this study used in silico metabolic engineering techniques, flux balance analysis (FBA), and flux variability analysis with target flux minimization (FATMIN), to model the metabolic perturbation of S. cerevisiae under the MFC-energy extraction. The FBA results showed that, in the cytoplasmic NADH-dependent mediated electron transfer (MET) mode, S. cerevisiae had a potential to produce currents at up to 5.781 A/gDW for the anaerobic and 6.193 A/gDW for the aerobic environments. The FATMIN results showed that the aerobic and anaerobic metabolisms are resilient, relying on six and five contributing reactions respectively for high current production. Two reactions, catalyzed by glutamate dehydrogenase (NAD) (EC 1.4.1.3) and methylene tetrahydrofolate dehydrogenase (NAD) (EC 1.5.1.5), were shared in both current-production modes and contributed to over 80% of the identified maximum current outputs. It is also shown that the NADH regeneration was much less energy costly than biomass production rate. Taken together, our finding suggests that S. cerevisiae should receive more research effort for MFC electricity production.

Version 6 of the consensus yeast metabolic network refines biochemical coverage and improves model performance

Updates to maintain a state-of-the art reconstruction of the yeast metabolic network are essential to reflect our understanding of yeast metabolism and functional organization, to eliminate any inaccuracies identified in earlier iterations, to improve predictive accuracy and to continue to expand into novel subsystems to extend the comprehensiveness of the model. Here, we present version 6 of the consensus yeast metabolic network (Yeast 6) as an update to the community effort to computationally reconstruct the genome-scale metabolic network of Saccharomyces cerevisiae S288c. Yeast 6 comprises 1458 metabolites participating in 1888 reactions, which are annotated with 900 yeast genes encoding the catalyzing enzymes. Compared with Yeast 5, Yeast 6 demonstrates improved sensitivity, specificity and positive and negative predictive values for predicting gene essentiality in glucose-limited aerobic conditions when analyzed with flux balance analysis. Additionally, Yeast 6 improves the accuracy of predicting the likelihood that a mutation will cause auxotrophy. The network reconstruction is available as a Systems Biology Markup Language (SBML) file enriched with Minimium Information Requested in the Annotation of Biochemical Models (MIRIAM)-compliant annotations. Small- and macromolecules in the network are referenced to authoritative databases such as Uniprot or ChEBI. Molecules and reactions are also annotated with appropriate publications that contain supporting evidence. Yeast 6 is freely available at http://yeast.sf.net/ as three separate SBML files: a model using the SBML level 3 Flux Balance Constraint package, a model compatible with the MATLAB® COBRA Toolbox for backward compatibility and a reconstruction containing only reactions for which there is experimental evidence (without the non-biological reactions necessary for simulating growth). Database URL: http://yeast.sf.net/

An analysis of a 'community-driven' reconstruction of the human metabolic network

Following a strategy similar to that used in baker's yeast (Herrgård et al. Nat Biotechnol 26:1155-1160, 2008). A consensus yeast metabolic network obtained from a community approach to systems biology (Herrgård et al. 2008; Dobson et al. BMC Syst Biol 4:145, 2010). Further developments towards a genome-scale metabolic model of yeast (Dobson et al. 2010; Heavner et al. BMC Syst Biol 6:55, 2012). Yeast 5-an expanded reconstruction of the Saccharomyces cerevisiae metabolic network (Heavner et al. 2012) and in Salmonella typhimurium (Thiele et al. BMC Syst Biol 5:8, 2011). A community effort towards a knowledge-base and mathematical model of the human pathogen Salmonellatyphimurium LT2 (Thiele et al. 2011), a recent paper (Thiele et al. Nat Biotechnol 31:419-425, 2013). A community-driven global reconstruction of human metabolism (Thiele et al. 2013) described a much improved 'community consensus' reconstruction of the human metabolic network, called Recon 2, and the authors (that include the present ones) have made it freely available via a database at http://humanmetabolism.org/ and in SBML format at Biomodels (http://identifiers.org/biomodels.db/MODEL1109130000). This short analysis summarises the main findings, and suggests some approaches that will be able to exploit the availability of this model to advantage.

Revising the Representation of Fatty Acid, Glycerolipid, and Glycerophospholipid Metabolism in the Consensus Model of Yeast Metabolism

Genome-scale metabolic models are built using information from an organism's annotated genome and, correspondingly, information on reactions catalyzed by the set of metabolic enzymes encoded by the genome. These models have been successfully applied to guide metabolic engineering to increase production of metabolites of industrial interest. Congruity between simulated and experimental metabolic behavior is influenced by the accuracy of the representation of the metabolic network in the model. In the interest of applying the consensus model of Saccharomyces cerevisiae metabolism for increased productivity of triglycerides, we manually evaluated the representation of fatty acid, glycerophospholipid, and glycerolipid metabolism in the consensus model (Yeast v6.0). These areas of metabolism were chosen due to their tightly interconnected nature to triglyceride synthesis. Manual curation was facilitated by custom MATLAB functions that return information contained in the model for reactions associated with genes and metabolites within the stated areas of metabolism. Through manual curation, we have identified inconsistencies between information contained in the model and literature knowledge. These inconsistencies include incorrect gene-reaction associations, improper definition of substrates/products in reactions, inappropriate assignments of reaction directionality, nonfunctional β-oxidation pathways, and missing reactions relevant to the synthesis and degradation of triglycerides. Suggestions to amend these inconsistencies in the Yeast v6.0 model can be implemented through a MATLAB script provided in theSupplementary Materials, Supplementary Data S1(Supplementary Data are available online at www.liebertpub.com/ind).

Reconstruction and validation of a genome-scale metabolic model for the filamentous fungus Neurospora crassa using FARM

The filamentous fungus Neurospora crassa played a central role in the development of twentieth-century genetics, biochemistry and molecular biology, and continues to serve as a model organism for eukaryotic biology. Here, we have reconstructed a genome-scale model of its metabolism. This model consists of 836 metabolic genes, 257 pathways, 6 cellular compartments, and is supported by extensive manual curation of 491 literature citations. To aid our reconstruction, we developed three optimization-based algorithms, which together comprise Fast Automated Reconstruction of Metabolism (FARM). These algorithms are: LInear MEtabolite Dilution Flux Balance Analysis (limed-FBA), which predicts flux while linearly accounting for metabolite dilution; One-step functional Pruning (OnePrune), which removes blocked reactions with a single compact linear program; and Consistent Reproduction Of growth/no-growth Phenotype (CROP), which reconciles differences between in silico and experimental gene essentiality faster than previous approaches. Against an independent test set of more than 300 essential/non-essential genes that were not used to train the model, the model displays 93% sensitivity and specificity. We also used the model to simulate the biochemical genetics experiments originally performed on Neurospora by comprehensively predicting nutrient rescue of essential genes and synthetic lethal interactions, and we provide detailed pathway-based mechanistic explanations of our predictions. Our model provides a reliable computational framework for the integration and interpretation of ongoing experimental efforts in Neurospora, and we anticipate that our methods will substantially reduce the manual effort required to develop high-quality genome-scale metabolic models for other organisms.

Few organisms have been as foundational to the development of modern genetics and cellular metabolism as Neurospora crassa. Given the wealth of knowledge available for this filamentous fungus, the effort required to manually curate a high-quality genome-scale metabolic reconstruction would be daunting. To aid the reconstruction process, we developed three optimization-based algorithms. The first algorithm predicts flux while linearly accounting for metabolite dilution; the second algorithm removes blocked reactions with one compact linear program; and the third algorithm reconciles differences between in silico predictions and experimental observations of mutant viability. We have used these algorithms to develop the first genome-scale metabolic model for Neurospora. We have validated the accuracy of our model against an independent test set of more than 300 growth/no-growth phenotypes, and our model displays 93% sensitivity and specificity. Simulating the biochemical genetics experiments originally performed on Neurospora, we comprehensively predicted essential genes, nutrient rescues of auxotroph mutants and synthetic lethal interactions. With these predictions, we provide potential mechanistic insight into known mutant phenotypes, and testable hypotheses for novel mutant phenotypes. The model, the algorithms and the testable hypotheses provide a computational foundation for the study of Neurospora crassa metabolism.

A model of yeast glycolysis based on a consistent kinetic characterisation of all its enzymes

We present an experimental and computational pipeline for the generation of kinetic models of metabolism, and demonstrate its application to glycolysis in Saccharomyces cerevisiae. Starting from an approximate mathematical model, we employ a “cycle of knowledge” strategy, identifying the steps with most control over flux. Kinetic parameters of the individual isoenzymes within these steps are measured experimentally under a standardised set of conditions. Experimental strategies are applied to establish a set of in vivo concentrations for isoenzymes and metabolites. The data are integrated into a mathematical model that is used to predict a new set of metabolite concentrations and reevaluate the control properties of the system. This bottom-up modelling study reveals that control over the metabolic network most directly involved in yeast glycolysis is more widely distributed than previously thought.

Responses to light intensity in a genome-scale model of rice metabolism

We describe the construction and analysis of a genome-scale metabolic model representing a developing leaf cell of rice (Oryza sativa) primarily derived from the annotations in the RiceCyc database. We used flux balance analysis to determine that the model represents a network capable of producing biomass precursors (amino acids, nucleotides, lipid, starch, cellulose, and lignin) in experimentally reported proportions, using carbon dioxide as the sole carbon source. We then repeated the analysis over a range of photon flux values to examine responses in the solutions. The resulting flux distributions show that (1) redox shuttles between the chloroplast, cytosol, and mitochondrion may play a significant role at low light levels, (2) photorespiration can act to dissipate excess energy at high light levels, and (3) the role of mitochondrial metabolism is likely to vary considerably according to the balance between energy demand and availability. It is notable that these organelle interactions, consistent with many experimental observations, arise solely as a result of the need for mass and energy balancing without any explicit assumptions concerning kinetic or other regulatory mechanisms.

Mapping condition-dependent regulation of metabolism in yeast through genome-scale modeling

Background

The genome-scale metabolic model of Saccharomyces cerevisiae, first presented in 2003, was the first genome-scale network reconstruction for a eukaryotic organism. Since then continuous efforts have been made in order to improve and expand the yeast metabolic network.

Results

Here we present iTO977, a comprehensive genome-scale metabolic model that contains more reactions, metabolites and genes than previous models. The model was constructed based on two earlier reconstructions, namely iIN800 and the consensus network, and then improved and expanded using gap-filling methods and by introducing new reactions and pathways based on studies of the literature and databases. The model was shown to perform well both for growth simulations in different media and gene essentiality analysis for single and double knock-outs. Further, the model was used as a scaffold for integrating transcriptomics, and flux data from four different conditions in order to identify transcriptionally controlled reactions, i.e. reactions that change both in flux and transcription between the compared conditions.

Conclusion

We present a new yeast model that represents a comprehensive up-to-date collection of knowledge on yeast metabolism. The model was used for simulating the yeast metabolism under four different growth conditions and experimental data from these four conditions was integrated to the model. The model together with experimental data is a useful tool to identify condition-dependent changes of metabolism between different environmental conditions.

A community-driven global reconstruction of human metabolism

Multiple models of human metabolism have been reconstructed, but each represents only a subset of our knowledge. Here we describe Recon 2, a community-driven, consensus 'metabolic reconstruction', which is the most comprehensive representation of human metabolism that is applicable to computational modeling. Compared with its predecessors, the reconstruction has improved topological and functional features, including ∼2× more reactions and ∼1.7× more unique metabolites. Using Recon 2 we predicted changes in metabolite biomarkers for 49 inborn errors of metabolism with 77% accuracy when compared to experimental data. Mapping metabolomic data and drug information onto Recon 2 demonstrates its potential for integrating and analyzing diverse data types. Using protein expression data, we automatically generated a compendium of 65 cell type-specific models, providing a basis for manual curation or investigation of cell-specific metabolic properties. Recon 2 will facilitate many future biomedical studies and is freely available at http://humanmetabolism.org/.

Iron, glucose and intrinsic factors alter sphingolipid composition as yeast cells enter stationary phase

Survival of Saccharomyces cerevisiae cells, like most microorganisms, requires switching from a rapidly dividing to a non-dividing or stationary state. To further understand how cells navigate this switch, we examined sphingolipids since they are key structural elements of membranes and also regulate signaling pathways vital for survival. During and after the switch to a non-dividing state there is a large increase in total free and sphingolipid-bound long chain-bases and an even larger increase in free and bound C20-long-chain bases, which are nearly undetectable in dividing cells. These changes are due to intrinsic factors including Orm1 and Orm2, ceramide synthase, Lcb4 kinase and the Tsc3 subunit of serine palmitoyltransferase as well as extrinsic factors including glucose and iron. Lowering the concentration of glucose, a form of calorie restriction, decreases the level of LCBs, which is consistent with the idea that reducing the level of some sphingolipids enhances lifespan. In contrast, iron deprivation increases LCB levels and decreases long term survival; however, these phenomena may not be related because iron deprivation disrupts many metabolic pathways. The correlation between increased LCBs and shorter lifespan is unsupported at this time. The physiological rise in LCBs that we observe may serve to modulate nutrient transporters and possibly other membrane phenomena that contribute to enhanced stress resistance and survival in stationary phase.

What mRNA Abundances Can Tell us about Metabolism

Inferring decreased or increased metabolic functions from transcript profiles is at first sight a bold and speculative attempt because of the functional layers in between: proteins, enzymatic activities, and reaction fluxes. However, the growing interest in this field can easily be explained by two facts: the high quality of genome-scale metabolic network reconstructions and the highly developed technology to obtain genome-covering RNA profiles. Here, an overview of important algorithmic approaches is given by means of criteria by which published procedures can be classified. The frontiers of the methods are sketched and critical voices are being heard. Finally, an outlook for the prospects of the field is given.

Improving metabolic flux predictions using absolute gene expression data

Background

Constraint-based analysis of genome-scale metabolic models typically relies upon maximisation of a cellular objective function such as the rate or efficiency of biomass production. Whilst this assumption may be valid in the case of microorganisms growing under certain conditions, it is likely invalid in general, and especially for multicellular organisms, where cellular objectives differ greatly both between and within cell types. Moreover, for the purposes of biotechnological applications, it is normally the flux to a specific metabolite or product that is of interest rather than the rate of production of biomass per se.

Results

An alternative objective function is presented, that is based upon maximising the correlation between experimentally measured absolute gene expression data and predicted internal reaction fluxes. Using quantitative transcriptomics data acquired from Saccharomyces cerevisiae cultures under two growth conditions, the method outperforms traditional approaches for predicting experimentally measured exometabolic flux that are reliant upon maximisation of the rate of biomass production.

Conclusion

Due to its improved prediction of experimentally measured metabolic fluxes, and of its lack of a requirement for knowledge of the biomass composition of the organism under the conditions of interest, the approach is likely to be of rather general utility. The method has been shown to predict fluxes reliably in single cellular systems. Subsequent work will investigate the method’s ability to generate condition- and tissue-specific flux predictions in multicellular organisms.