BioMet Toolbox: genome-wide analysis of metabolism

mfapy: An open-source Python package for 13C-based metabolic flux analysis

¹³C-based metabolic flux analysis (¹³C-MFA) is an essential tool for estimating intracellular metabolic flux levels in metabolic engineering and biology. In ¹³C-MFA, a metabolic flux distribution that explains the observed isotope labeling data was computationally estimated using a non-linear optimization method. Herein, we report the development of mfapy, an open-source Python package developed for more flexibility and extensibility for ¹³C-MFA. mfapy compels users to write a customized Python code by describing each step in the data analysis procedures of the isotope labeling experiments. The flexibility and extensibility provided by mfapy can support trial-and-error performance in the routine estimation of metabolic flux distributions, experimental design by computer simulations of ¹³C-MFA experiments, and development of new data analysis techniques for stable isotope labeling experiments. mfapy is available to the public from the Github repository (https://github.com/fumiomatsuda/mfapy).

•

An open-source Python package, mfapy, is developed for ¹³C-MFA.

•

mfapy enables users to write Python codes for data analysis procedures of ¹³C-MFA.

•

mfapy has a flexibility and extensibility to support various data analysis procedures.

•

Computer simulations of ¹³C-MFA experiments is supported for experimental design.

AutoAnalyze in Systems Biology

AutoAnalyze is a highly customizable framework for the visualization and analysis of large-scale model graphs. Originally developed for use in the automotive domain, it also supports efficient computation within molecular networks represented by reaction equations. A static analysis approach is used for efficient treatment-condition-specific simulation. The chosen method relies on the computation of a global network data-flow resulting from the evaluation of individual genetic data. The approach facilitates complex analyses of biological components from a molecular network under specific therapeutic perturbations, as demonstrated in a case study. In addition to simulating the complex networks in a stable and reproducible way, kinetic constants can also be fine-tuned using a genetic algorithm and built-in statistical tools.

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

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

Understanding metabolism with flux analysis: From theory to application

Quantitative and qualitative knowledge of metabolic rates (i.e. fluxes) over a metabolic network and in specific cellular compartments gives insights into the regulation of metabolism and helps to understand the contribution of metabolic alterations to pathology. In this review we introduce methodology to resolve metabolic fluxes from stable isotope labeling and relevant techniques in model development, model simplification, flux uncertainty analysis and experimental design that together is termed metabolic flux analysis. Finally we discuss applications using metabolic flux analysis to elucidate mechanisms pertinent to tumor cell metabolism. We hope that this review gives the readers a brief introduction of how flux analysis is conducted, how technical issues related to it are addressed, and how its application has contributed to our knowledge of tumor cell metabolism.

Synergizing 13C Metabolic Flux Analysis and Metabolic Engineering for Biochemical Production

Metabolic engineering of industrial microorganisms to produce chemicals, fuels, and drugs has attracted increasing interest as it provides an environment-friendly and renewable route that does not depend on depleting petroleum sources. However, the microbial metabolism is so complex that metabolic engineering efforts often have difficulty in achieving a satisfactory yield, titer, or productivity of the target chemical. To overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been developed to investigate rigorously the cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, ¹³C-MFA has been widely used in academic labs and the biotechnology industry to pinpoint the key issues related to microbial-based chemical production and to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this chapter we introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied to synergize with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production.

An integrative analysis of transcriptomic response of ethanol tolerant strains to ethanol in Saccharomyces cerevisiae

The accumulation of ethanol is one of the main environmental stresses that Saccharomyces cerevisiae cells are exposed to in industrial alcoholic beverage and bioethanol production processes. Despite the known impacts of ethanol, the molecular mechanisms underlying ethanol tolerance are still not fully understood. Novel gene targets leading to ethanol tolerance were previously identified via a network approach and the investigations of the deletions of these genes resulted in the improved ethanol tolerance of pmt7Δ/pmt7Δ and yhl042wΔ/yhl042wΔ strains. In the present study, an integrative system based approach was used to investigate the global transcriptional changes in these two ethanol tolerant strains in response to ethanol and hence to elucidate the mechanisms leading to the observed tolerant phenotypes. In addition to strain specific biological processes, a number of common and already reported biological processes were found to be affected in the reference and both ethanol tolerant strains. However, the integrative analysis of the transcriptome with the transcriptional regulatory network and the ethanol tolerance network revealed that each ethanol tolerant strain had a specific organization of the transcriptomic response. Transcription factors around which most important changes occur were determined and active subnetworks in response to ethanol and functional clusters were identified in all strains.

Investigate the Metabolic Reprogramming of Saccharomyces cerevisiae for Enhanced Resistance to Mixed Fermentation Inhibitors via 13C Metabolic Flux Analysis

The fermentation inhibitors from the pretreatment of lignocellulosic materials, e.g., acetic acid and furfural, are notorious due to their negative effects on the cell growth and chemical production. However, the metabolic reprogramming of the cells under these stress conditions, especially metabolic response for resistance to mixed inhibitors, has not been systematically investigated and remains mysterious. Therefore, in this study, 13C metabolic flux analysis (13C-MFA), a powerful tool to elucidate the intracellular carbon flux distributions, has been applied to two Saccharomyces cerevisiae strains with different tolerances to the inhibitors under acetic acid, furfural, and mixed (i.e., acetic acid and furfural) stress conditions to unravel the key metabolic responses. By analyzing the intracellular carbon fluxes as well as the energy and cofactor utilization under different conditions, we uncovered varied metabolic responses to different inhibitors. Under acetate stress, ATP and NADH production was slightly impaired, while NADPH tended towards overproduction. Under furfural stress, ATP and cofactors (including both NADH and NADPH) tended to be overproduced. However, under dual-stress condition, production of ATP and cofactors was severely impaired due to synergistic stress caused by the simultaneous addition of two fermentation inhibitors. Such phenomenon indicated the pivotal role of the energy and cofactor utilization in resisting the mixed inhibitors of acetic acid and furfural. Based on the discoveries, valuable insights are provided to improve the tolerance of S. cerevisiae strain and further enhance lignocellulosic fermentation.

Genetic regulation and manipulation for natural product discovery

Natural products are an important source of modern medical development, e.g., antibiotics, anticancers, immune modulators, etc. and will continue to be a powerful driving force for the discovery of novel potential drugs. In the heterologous hosts, natural products are biosynthesized using dedicated metabolic networks. By gene engineering, pathway reconstructing, and enzyme engineering, metabolic networks can be modified to synthesize novel compounds containing enhanced structural feature or produce a large quantity of known valuable bioactive compounds. The review introduces some important technical platforms and relevant examples of genetic regulation and manipulation to improve natural product titers or drive novel secondary metabolite discoveries.

13C-Metabolic Flux Analysis: An Accurate Approach to Demystify Microbial Metabolism for Biochemical Production

Metabolic engineering of various industrial microorganisms to produce chemicals, fuels, and drugs has raised interest since it is environmentally friendly, sustainable, and independent of nonrenewable resources. However, microbial metabolism is so complex that only a few metabolic engineering efforts have been able to achieve a satisfactory yield, titer or productivity of the target chemicals for industrial commercialization. In order to overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been continuously developed and widely applied to rigorously investigate cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, many ¹³C-MFA studies have been performed in academic labs and biotechnology industries to pinpoint key issues related to microbe-based chemical production. Insightful information about the metabolic rewiring has been provided to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this review, we will introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied via integration with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production for various host microorganisms

Comparative Analysis of Yeast Metabolic Network Models Highlights Progress, Opportunities for Metabolic Reconstruction

We have compared 12 genome-scale models of the Saccharomyces cerevisiae metabolic network published since 2003 to evaluate progress in reconstruction of the yeast metabolic network. We compared the genomic coverage, overlap of annotated metabolites, predictive ability for single gene essentiality with a selection of model parameters, and biomass production predictions in simulated nutrient-limited conditions. We have also compared pairwise gene knockout essentiality predictions for 10 of these models. We found that varying approaches to model scope and annotation reflected the involvement of multiple research groups in model development; that single-gene essentiality predictions were affected by simulated medium, objective function, and the reference list of essential genes; and that predictive ability for single-gene essentiality did not correlate well with predictive ability for our reference list of synthetic lethal gene interactions (R = 0.159). We conclude that the reconstruction of the yeast metabolic network is indeed gradually improving through the iterative process of model development, and there remains great opportunity for advancing our understanding of biology through continued efforts to reconstruct the full biochemical reaction network that constitutes yeast metabolism. Additionally, we suggest that there is opportunity for refining the process of deriving a metabolic model from a metabolic network reconstruction to facilitate mechanistic investigation and discovery. This comparative study lays the groundwork for developing improved tools and formalized methods to quantitatively assess metabolic network reconstructions independently of any particular model application, which will facilitate ongoing efforts to advance our understanding of the relationship between genotype and cellular phenotype.

Reporter pathway analysis from transcriptome data: Metabolite-centric versus Reaction-centric approach

A systems-based investigation of the effect of perturbations on metabolic machinery is crucial to elucidate the mechanism behind perturbations. One way to investigate the perturbation-induced changes within the cell metabolism is to focus on pathway-level effects. In this study, three different perturbation types (genetic, environmental and disease-based) are analyzed to compute a list of reporter pathways, metabolic pathways which are significantly affected from a perturbation. The most common omics data type, transcriptome, is used as an input to the bioinformatic analysis. The pathways are scored by two alternative approaches: by averaging the changes in the expression levels of the genes controlling the associated reactions (reaction-centric), and by averaging the changes in the associated metabolites which were scored based on the associated genes (metabolite-centric). The analysis reveals the superiority of the novel metabolite-centric approach over the commonly used reaction-centric approach since it is based on metabolites which better represent the cross-talk among different pathways, enabling a more global and realistic cataloguing of network-wide perturbation effects.

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.

Computational approaches for microalgal biofuel optimization: a review

The increased demand and consumption of fossil fuels have raised interest in finding renewable energy sources throughout the globe. Much focus has been placed on optimizing microorganisms and primarily microalgae, to efficiently produce compounds that can substitute for fossil fuels. However, the path to achieving economic feasibility is likely to require strain optimization through using available tools and technologies in the fields of systems and synthetic biology. Such approaches invoke a deep understanding of the metabolic networks of the organisms and their genomic and proteomic profiles. The advent of next generation sequencing and other high throughput methods has led to a major increase in availability of biological data. Integration of such disparate data can help define the emergent metabolic system properties, which is of crucial importance in addressing biofuel production optimization. Herein, we review major computational tools and approaches developed and used in order to potentially identify target genes, pathways, and reactions of particular interest to biofuel production in algae. As the use of these tools and approaches has not been fully implemented in algal biofuel research, the aim of this review is to highlight the potential utility of these resources toward their future implementation in algal research.

Generation and Evaluation of a Genome-Scale Metabolic Network Model of Synechococcus elongatus PCC7942

The reconstruction of genome-scale metabolic models and their applications represent a great advantage of systems biology. Through their use as metabolic flux simulation models, production of industrially-interesting metabolites can be predicted. Due to the growing number of studies of metabolic models driven by the increasing genomic sequencing projects, it is important to conceptualize steps of reconstruction and analysis. We have focused our work in the cyanobacterium Synechococcus elongatus PCC7942, for which several analyses and insights are unveiled. A comprehensive approach has been used, which can be of interest to lead the process of manual curation and genome-scale metabolic analysis. The final model, iSyf715 includes 851 reactions and 838 metabolites. A biomass equation, which encompasses elementary building blocks to allow cell growth, is also included. The applicability of the model is finally demonstrated by simulating autotrophic growth conditions of Synechococcus elongatus PCC7942.

A system based network approach to ethanol tolerance in Saccharomyces cerevisiae

Background

Saccharomyces cerevisiae has been widely used for bio-ethanol production and development of rational genetic engineering strategies leading both to the improvement of productivity and ethanol tolerance is very important for cost-effective bio-ethanol production. Studies on the identification of the genes that are up- or down-regulated in the presence of ethanol indicated that the genes may be involved to protect the cells against ethanol stress, but not necessarily required for ethanol tolerance.

Results

In the present study, a novel network based approach was developed to identify candidate genes involved in ethanol tolerance. Protein-protein interaction (PPI) network associated with ethanol tolerance (tETN) was reconstructed by integrating PPI data with Gene Ontology (GO) terms. Modular analysis of the constructed networks revealed genes with no previously reported experimental evidence related to ethanol tolerance and resulted in the identification of 17 genes with previously unknown biological functions. We have randomly selected four of these genes and deletion strains of two genes (YDR307W and YHL042W) were found to exhibit improved tolerance to ethanol when compared to wild type strain.

The genome-wide transcriptomic response of yeast cells to the deletions of YDR307W and YHL042W in the absence of ethanol revealed that the deletion of YDR307W and YHL042W genes resulted in the transcriptional re-programming of the metabolism resulting from a mis-perception of the nutritional environment. Yeast cells perceived an excess amount of glucose and a deficiency of methionine or sulfur in the absence of YDR307W and YHL042W, respectively, possibly resulting from a defect in the nutritional sensing and signaling or transport mechanisms. Mutations leading to an increase in ribosome biogenesis were found to be important for the improvement of ethanol tolerance. Modulations of chronological life span were also identified to contribute to ethanol tolerance in yeast.

Conclusions

The system based network approach developed allows the identification of novel gene targets for improved ethanol tolerance and supports the highly complex nature of ethanol tolerance in yeast.

(Post-)genomics approaches in fungal research

To date, hundreds of fungal genomes have been sequenced and many more are in progress. This wealth of genomic information has provided new directions to study fungal biodiversity. However, to further dissect and understand the complicated biological mechanisms involved in fungal life styles, functional studies beyond genomes are required. Thanks to the developments of current -omics techniques, it is possible to produce large amounts of fungal functional data in a high-throughput fashion (e.g. transcriptome, proteome, etc.). The increasing ease of creating -omics data has also created a major challenge for downstream data handling and analysis. Numerous databases, tools and software have been created to meet this challenge. Facing such a richness of techniques and information, hereby we provide a brief roadmap on current wet-lab and bioinformatics approaches to study functional genomics in fungi.

OpenMebius: an open source software for isotopically nonstationary 13C-based metabolic flux analysis

The in vivo measurement of metabolic flux by ¹³C-based metabolic flux analysis (¹³C-MFA) provides valuable information regarding cell physiology. Bioinformatics tools have been developed to estimate metabolic flux distributions from the results of tracer isotopic labeling experiments using a ¹³C-labeled carbon source. Metabolic flux is determined by nonlinear fitting of a metabolic model to the isotopic labeling enrichment of intracellular metabolites measured by mass spectrometry. Whereas ¹³C-MFA is conventionally performed under isotopically constant conditions, isotopically nonstationary ¹³C metabolic flux analysis (INST-¹³C-MFA) has recently been developed for flux analysis of cells with photosynthetic activity and cells at a quasi-steady metabolic state (e.g., primary cells or microorganisms under stationary phase). Here, the development of a novel open source software for INST-¹³C-MFA on the Windows platform is reported. OpenMebius (Open source software for Metabolic flux analysis) provides the function of autogenerating metabolic models for simulating isotopic labeling enrichment from a user-defined configuration worksheet. Analysis using simulated data demonstrated the applicability of OpenMebius for INST-¹³C-MFA. Confidence intervals determined by INST-¹³C-MFA were less than those determined by conventional methods, indicating the potential of INST-¹³C-MFA for precise metabolic flux analysis. OpenMebius is the open source software for the general application of INST-¹³C-MFA.

Systematic analysis of transcription-level effects of neurodegenerative diseases on human brain metabolism by a newly reconstructed brain-specific metabolic network

Network-oriented analysis is essential to identify those parts of a cell affected by a given perturbation. The effect of neurodegenerative perturbations in the form of diseases of brain metabolism was investigated by using a newly reconstructed brain-specific metabolic network. The developed stoichiometric model correctly represents healthy brain metabolism, and includes 630 metabolic reactions in and between astrocytes and neurons, which are controlled by 570 genes. The integration of transcriptome data of six neurodegenerative diseases (Alzheimer's disease, Parkinson's disease, Huntington's disease, amyotrophic lateral sclerosis, multiple sclerosis, schizophrenia) with the model was performed to identify reporter features specific and common for these diseases, which revealed metabolites and pathways around which the most significant changes occur. The identified metabolites are potential biomarkers for the pathology of the related diseases. Our model indicated perturbations in oxidative stress, energy metabolism including TCA cycle and lipid metabolism as well as several amino acid related pathways, in agreement with the role of these pathways in the studied diseases. The computational prediction of transcription factors that commonly regulate the reporter metabolites was achieved through binding-site analysis. Literature support for the identified transcription factors such as USF1, SP1 and those from FOX families are known from the literature to have regulatory roles in the identified reporter metabolic pathways as well as in the neurodegenerative diseases. In essence, the reconstructed brain model enables the elucidation of effects of a perturbation on brain metabolism and the illumination of possible machineries in which a specific metabolite or pathway acts as a regulatory spot for cellular reorganization.

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.

Computer-aided design for metabolic engineering

The development and application of biotechnology-based strategies has had a great socio-economical impact and is likely to play a crucial role in the foundation of more sustainable and efficient industrial processes. Within biotechnology, metabolic engineering aims at the directed improvement of cellular properties, often with the goal of synthesizing a target chemical compound. The use of computer-aided design (CAD) tools, along with the continuously emerging advanced genetic engineering techniques have allowed metabolic engineering to broaden and streamline the process of heterologous compound-production. In this work, we review the CAD tools available for metabolic engineering with an emphasis, on retrosynthesis methodologies. Recent advances in genetic engineering strategies for pathway implementation and optimization are also reviewed as well as a range of bionalytical tools to validate in silico predictions. A case study applying retrosynthesis is presented as an experimental verification of the output from Retropath, the first complete automated computational pipeline applicable to metabolic engineering. Applying this CAD pipeline, together with genetic reassembly and optimization of culture conditions led to improved production of the plant flavonoid pinocembrin. Coupling CAD tools with advanced genetic engineering strategies and bioprocess optimization is crucial for enhanced product yields and will be of great value for the development of non-natural products through sustainable biotechnological processes.

Developing sustainable software solutions for bioinformatics by the " Butterfly" paradigm

Software design and sustainable software engineering are essential for the long-term development of bioinformatics software. Typical challenges in an academic environment are short-term contracts, island solutions, pragmatic approaches and loose documentation. Upcoming new challenges are big data, complex data sets, software compatibility and rapid changes in data representation. Our approach to cope with these challenges consists of iterative intertwined cycles of development (“ Butterfly” paradigm) for key steps in scientific software engineering. User feedback is valued as well as software planning in a sustainable and interoperable way. Tool usage should be easy and intuitive. A middleware supports a user-friendly Graphical User Interface (GUI) as well as a database/tool development independently. We validated the approach of our own software development and compared the different design paradigms in various software solutions.

In Pichia pastoris, growth rate regulates protein synthesis and secretion, mating and stress response

Protein production in yeasts is related to the specific growth rate μ. To elucidate on this correlation, we studied the transcriptome of Pichia pastoris at different specific growth rates by cultivating a strain secreting human serum albumin at μ = 0.015 to 0.15 h(-1) in glucose-limited chemostats. Genome-wide regulation revealed that translation-related as well as mitochondrial genes were upregulated with increasing μ, while autophagy and other proteolytic processes, carbon source-responsive genes and other targets of the TOR pathway as well as many transcriptional regulators were downregulated at higher μ. Mating and sporulation genes were most active at intermediate μ of 0.05 and 0.075 h(-1) . At very slow growth (μ = 0.015 h(-1) ) gene regulation differs significantly, affecting many transporters and glucose sensing. Analysis of a subset of genes related to protein folding and secretion reveals that unfolded protein response targets such as translocation, endoplasmic reticulum genes, and cytosolic chaperones are upregulated with increasing growth rate while proteolytic degradation of secretory proteins is downregulated. We conclude that a high μ positively affects specific protein secretion rates by acting on multiple cellular processes.

Investigating host-pathogen behavior and their interaction using genome-scale metabolic network models

Genome Scale Metabolic Modeling methods represent one way to compute whole cell function starting from the genome sequence of an organism and contribute towards understanding and predicting the genotype-phenotype relationship. About 80 models spanning all the kingdoms of life from archaea to eukaryotes have been built till date and used to interrogate cell phenotype under varying conditions. These models have been used to not only understand the flux distribution in evolutionary conserved pathways like glycolysis and the Krebs cycle but also in applications ranging from value added product formation in Escherichia coli to predicting inborn errors of Homo sapiens metabolism. This chapter describes a protocol that delineates the process of genome scale metabolic modeling for analysing host-pathogen behavior and interaction using flux balance analysis (FBA). The steps discussed in the process include (1) reconstruction of a metabolic network from the genome sequence, (2) its representation in a precise mathematical framework, (3) its translation to a model, and (4) the analysis using linear algebra and optimization. The methods for biological interpretations of computed cell phenotypes in the context of individual host and pathogen models and their integration are also discussed.

Software platforms to facilitate reconstructing genome-scale metabolic networks

System-level analyses of microbial metabolism are facilitated by genome-scale reconstructions of microbial biochemical networks. A reconstruction provides a structured representation of the biochemical transformations occurring within an organism, as well as the genes necessary to carry out these transformations, as determined by the annotated genome sequence and experimental data. Network reconstructions also serve as platforms for constraint-based computational techniques, which facilitate biological studies in a variety of applications, including evaluation of network properties, metabolic engineering and drug discovery. Bottom-up metabolic network reconstructions have been developed for dozens of organisms, but until recently, the pace of reconstruction has failed to keep up with advances in genome sequencing. To address this problem, a number of software platforms have been developed to automate parts of the reconstruction process, thereby alleviating much of the manual effort previously required. Here, we review four such platforms in the context of established guidelines for network reconstruction. While many steps of the reconstruction process have been successfully automated, some manual evaluation of the results is still required to ensure a high-quality reconstruction. Widespread adoption of these platforms by the scientific community is underway and will be further enabled by exchangeable formats across platforms.

Mapping condition-dependent regulation of lipid metabolism in Saccharomyces cerevisiae

Lipids play a central role in cellular function as constituents of membranes, as signaling molecules, and as storage materials. Although much is known about the role of lipids in regulating specific steps of metabolism, comprehensive studies integrating genome-wide expression data, metabolite levels, and lipid levels are currently lacking. Here, we map condition-dependent regulation controlling lipid metabolism in Saccharomyces cerevisiae by measuring 5636 mRNAs, 50 metabolites, 97 lipids, and 57 (13)C-reaction fluxes in yeast using a three-factor full-factorial design. Correlation analysis across eight environmental conditions revealed 2279 gene expression level-metabolite/lipid relationships that characterize the extent of transcriptional regulation in lipid metabolism relative to major metabolic hubs within the cell. To query this network, we developed integrative methods for correlation of multi-omics datasets that elucidate global regulatory signatures. Our data highlight many characterized regulators of lipid metabolism and reveal that sterols are regulated more at the transcriptional level than are amino acids. Beyond providing insights into the systems-level organization of lipid metabolism, we anticipate that our dataset and approach can join an emerging number of studies to be widely used for interrogating cellular systems through the combination of mathematical modeling and experimental biology.

Investigating the physiological response of Pichia (Komagataella) pastoris GS115 to the heterologous expression of misfolded proteins using chemostat cultures

Pichia pastoris is widely used as a host system for heterologous protein expression in both academia and industry. Production is typically accomplished by a fed-batch induction process that is known to have negative impacts on cell physiology that impose limits on both protein yields and quality. We have analysed recombinant protein production in chemostat cultures to understand the physiological responses associated with methanol-induced production of two human lysozyme variants with different degrees of misfolding by P. pastoris. Confounding variables associated with nutrient stress or growth-rate are minimised during steady-state growth in chemostats. Comparison of transcriptome-level data obtained during the non-inducing and inducing steady states identified changes in expression of only about 1 % of the genome during production of either an amyloidogenic human lysozyme variant prone to intracellular aggregation (I56T) or a misfolded but secretable variant (T70N), indicating near-complete acclimation to their production. A marked, but temporary, stress response involving both the unfolded protein response (UPR) and ER-associated degradation pathway was observed during the transient between steady states, particularly following induction of the T70N variant synthesis, and was accompanied by changes in expression of around 50 antisense transcripts. The results suggest that optimal heterologous protein production could best be achieved by a continuous process that minimises the number of methanol-induced transients experienced by the cultures. The processing of HAC1 mRNA required for the UPR was found to be constitutive in the culture conditions used, even in the absence of recombinant protein induction.

Synechocystis sp. PCC6803 metabolic models for the enhanced production of hydrogen

In the present economy, difficulties to access energy sources are real drawbacks to maintain our current lifestyle. In fact, increasing interests have been gathered around efficient strategies to use energy sources that do not generate high CO2 titers. Thus, science-funding agencies have invested more resources into research on hydrogen among other biofuels as interesting energy vectors. This article reviews present energy challenges and frames it into the present fuel usage landscape. Different strategies for hydrogen production are explained and evaluated. Focus is on biological hydrogen production; fermentation and photon-fuelled hydrogen production are compared. Mathematical models in biology can be used to assess, explore and design production strategies for industrially relevant metabolites, such as biofuels. We assess the diverse construction and uses of genome-scale metabolic models of cyanobacterium Synechocystis sp. PCC6803 to efficiently obtain biofuels. This organism has been studied as a potential photon-fuelled production platform for its ability to grow from carbon dioxide, water and photons, on simple culture media. Finally, we review studies that propose production strategies to weigh this organism's viability as a biofuel production platform. Overall, the work presented in this review unveils the industrial capabilities of cyanobacterium Synechocystis sp. PCC6803 to evolve interesting metabolites as a clean biofuel production platform.

Investigation of malic acid production in Aspergillus oryzae under nitrogen starvation conditions

Malic acid has great potential for replacing petrochemical building blocks in the future. For this application, high yields, rates, and titers are essential in order to sustain a viable biotechnological production process. Natural high-capacity malic acid producers like the malic acid producer Aspergillus flavus have so far been disqualified because of special growth requirements or the production of mycotoxins. As A. oryzae is a very close relative or even an ecotype of A. flavus, it is likely that its high malic acid production capabilities with a generally regarded as safe (GRAS) status may be combined with already existing large-scale fermentation experience. In order to verify the malic acid production potential, two wild-type strains, NRRL3485 and NRRL3488, were compared in shake flasks. As NRRL3488 showed a volumetric production rate twice as high as that of NRRL3485, this strain was selected for further investigation of the influence of two different nitrogen sources on malic acid secretion. The cultivation in lab-scale fermentors resulted in a higher final titer, 30.27 ± 1.05 g liter(-1), using peptone than the one of 22.27 ± 0.46 g liter(-1) obtained when ammonium was used. Through transcriptome analysis, a binding site similar to the one of the Saccharomyces cerevisiae yeast transcription factor Msn2/4 was identified in the upstream regions of glycolytic genes and the cytosolic malic acid production pathway from pyruvate via oxaloacetate to malate, which suggests that malic acid production is a stress response. Furthermore, the pyruvate carboxylase reaction was identified as a target for metabolic engineering, after it was confirmed to be transcriptionally regulated through the correlation of intracellular fluxes and transcriptional changes.

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.

Advances in metabolic pathway and strain engineering paving the way for sustainable production of chemical building blocks

Bio-based production of chemical building blocks from renewable resources is an attractive alternative to petroleum-based platform chemicals. Metabolic pathway and strain engineering is the key element in constructing robust microbial chemical factories within the constraints of cost effective production. Here we discuss how the development of computational algorithms, novel modules and methods, omics-based techniques combined with modeling refinement are enabling reduction in development time and thus advance the field of industrial biotechnology. We further discuss how recent technological developments contribute to the development of novel cell factories for the production of the building block chemicals: adipic acid, succinic acid and 3-hydroxypropionic acid.

The RAVEN toolbox and its use for generating a genome-scale metabolic model for Penicillium chrysogenum

We present the RAVEN (Reconstruction, Analysis and Visualization of Metabolic Networks) Toolbox: a software suite that allows for semi-automated reconstruction of genome-scale models. It makes use of published models and/or the KEGG database, coupled with extensive gap-filling and quality control features. The software suite also contains methods for visualizing simulation results and omics data, as well as a range of methods for performing simulations and analyzing the results. The software is a useful tool for system-wide data analysis in a metabolic context and for streamlined reconstruction of metabolic networks based on protein homology. The RAVEN Toolbox workflow was applied in order to reconstruct a genome-scale metabolic model for the important microbial cell factory Penicillium chrysogenum Wisconsin54-1255. The model was validated in a bibliomic study of in total 440 references, and it comprises 1471 unique biochemical reactions and 1006 ORFs. It was then used to study the roles of ATP and NADPH in the biosynthesis of penicillin, and to identify potential metabolic engineering targets for maximization of penicillin production.

Genome-scale models (GEMs) are large stoichiometric models of cell metabolism, where the goal is to incorporate every metabolic transformation that an organism can perform. Such models have been extensively used for the study of bacterial metabolism, in particular for metabolic engineering purposes. More recently, the use of GEMs for eukaryotic organisms has become increasingly widespread. Since these models typically involve thousands of metabolic reactions, the reconstruction and validation of them can be a very complex task. We have developed a software suite, RAVEN Toolbox, which aims at automating parts of the reconstruction process in order to allow for faster reconstruction of high-quality GEMs. The software is particularly well suited for reconstruction of models for eukaryotic organisms, due to how it deals with sub-cellular localization of reactions. We used the software for reconstructing a model of the filamentous fungi Penicillium chrysogenum, the organism used in penicillin production and an important microbial cell factory. The resulting model was validated through an extensive literature survey and by comparison with published fermentation data. The model was used for the identification of transcriptionally regulated metabolic bottlenecks in order to increase the yield in penicillin fermentations. In this paper we present the RAVEN Toolbox and the GEM for P. chrysogenum.

Systems biology for molecular life sciences and its impact in biomedicine

Modern systems biology is already contributing to a radical transformation of molecular life sciences and biomedicine, and it is expected to have a real impact in the clinical setting in the next years. In this review, the emergence of systems biology is contextualized with a historic overview, and its present state is depicted. The present and expected future contribution of systems biology to the development of molecular medicine is underscored. Concerning the present situation, this review includes a reflection on the "inflation" of biological data and the urgent need for tools and procedures to make hidden information emerge. Descriptions of the impact of networks and models and the available resources and tools for applying them in systems biology approaches to molecular medicine are provided as well. The actual current impact of systems biology in molecular medicine is illustrated, reviewing two cases, namely, those of systems pharmacology and cancer systems biology. Finally, some of the expected contributions of systems biology to the immediate future of molecular medicine are commented.

Enriching the gene set analysis of genome-wide data by incorporating directionality of gene expression and combining statistical hypotheses and methods

Gene set analysis (GSA) is used to elucidate genome-wide data, in particular transcriptome data. A multitude of methods have been proposed for this step of the analysis, and many of them have been compared and evaluated. Unfortunately, there is no consolidated opinion regarding what methods should be preferred, and the variety of available GSA software and implementations pose a difficulty for the end-user who wants to try out different methods. To address this, we have developed the R package Piano that collects a range of GSA methods into the same system, for the benefit of the end-user. Further on we refine the GSA workflow by using modifications of the gene-level statistics. This enables us to divide the resulting gene set P-values into three classes, describing different aspects of gene expression directionality at gene set level. We use our fully implemented workflow to investigate the impact of the individual components of GSA by using microarray and RNA-seq data. The results show that the evaluated methods are globally similar and the major separation correlates well with our defined directionality classes. As a consequence of this, we suggest to use a consensus scoring approach, based on multiple GSA runs. In combination with the directionality classes, this constitutes a more thorough basis for an enriched biological interpretation.

Genome-scale model management and comparison

Spurred by recent innovations in genome sequencing, the reconstruction of genome-scale models has increased in recent years. Genome-scale models are now available for a wide range of organisms, and models have been successfully applied to a number of research topics including metabolic engineering, genome annotation, biofuel production, and interpretation of omics data sets. The challenge is how to manage the large amount of data in genome-scale models and perform comparative analysis to gain new biological insights. In this chapter, important standards for genome-scale modeling are outlined. Furthermore, management strategies as well as existing repository and construction tools are discussed. As the comparison of models is an important aspect during the development and analysis stages, available methods are presented and existing software solutions are reviewed.

13CFLUX2--high-performance software suite for (13)C-metabolic flux analysis

Summary:¹³C-based metabolic flux analysis (¹³C-MFA) is the state-of-the-art method to quantitatively determine in vivo metabolic reaction rates in microorganisms. 13CFLUX2 contains all tools for composing flexible computational ¹³C-MFA workflows to design and evaluate carbon labeling experiments. A specially developed XML language, FluxML, highly efficient data structures and simulation algorithms achieve a maximum of performance and effectiveness. Support of multicore CPUs, as well as compute clusters, enables scalable investigations. 13CFLUX2 outperforms existing tools in terms of universality, flexibility and built-in features. Therewith, 13CFLUX2 paves the way for next-generation high-resolution ¹³C-MFA applications on the large scale.

Availability and implementation: 13CFLUX2 is implemented in C++ (ISO/IEC 14882 standard) with Java and Python add-ons to run under Linux/Unix. A demo version and binaries are available at www.13cflux.net.

Contact:info@13cflux.net or k.noeh@fz-juelich.de

Supplementary information:Supplementary data are available at Bioinformatics online.

MicrobesFlux: a web platform for drafting metabolic models from the KEGG database

Background

Concurrent with the efforts currently underway in mapping microbial genomes using high-throughput sequencing methods, systems biologists are building metabolic models to characterize and predict cell metabolisms. One of the key steps in building a metabolic model is using multiple databases to collect and assemble essential information about genome-annotations and the architecture of the metabolic network for a specific organism. To speed up metabolic model development for a large number of microorganisms, we need a user-friendly platform to construct metabolic networks and to perform constraint-based flux balance analysis based on genome databases and experimental results.

Results

We have developed a semi-automatic, web-based platform (MicrobesFlux) for generating and reconstructing metabolic models for annotated microorganisms. MicrobesFlux is able to automatically download the metabolic network (including enzymatic reactions and metabolites) of ~1,200 species from the KEGG database (Kyoto Encyclopedia of Genes and Genomes) and then convert it to a metabolic model draft. The platform also provides diverse customized tools, such as gene knockouts and the introduction of heterologous pathways, for users to reconstruct the model network. The reconstructed metabolic network can be formulated to a constraint-based flux model to predict and analyze the carbon fluxes in microbial metabolisms. The simulation results can be exported in the SBML format (The Systems Biology Markup Language). Furthermore, we also demonstrated the platform functionalities by developing an FBA model (including 229 reactions) for a recent annotated bioethanol producer, Thermoanaerobacter sp. strain X514, to predict its biomass growth and ethanol production.

Conclusion

MicrobesFlux is an installation-free and open-source platform that enables biologists without prior programming knowledge to develop metabolic models for annotated microorganisms in the KEGG database. Our system facilitates users to reconstruct metabolic networks of organisms based on experimental information. Through human-computer interaction, MicrobesFlux provides users with reasonable predictions of microbial metabolism via flux balance analysis. This prototype platform can be a springboard for advanced and broad-scope modeling of complex biological systems by integrating other “omics” data or ¹³ C- metabolic flux analysis results. MicrobesFlux is available at http://tanglab.engineering.wustl.edu/static/MicrobesFlux.html and will be continuously improved based on feedback from users.

GEMSiRV: a software platform for GEnome-scale metabolic model simulation, reconstruction and visualization

MOTIVATION

Genome-scale metabolic network models have become an indispensable part of the increasingly important field of systems biology. Metabolic systems biology studies usually include three major components-network model construction, objective- and experiment-guided model editing and visualization, and simulation studies based mainly on flux balance analyses. Bioinformatics tools are required to facilitate these complicated analyses. Although some of the required functions have been served separately by existing tools, a free software resource that simultaneously serves the needs of the three major components is not yet available.

RESULTS

Here we present a software platform, GEMSiRV (GEnome-scale Metabolic model Simulation, Reconstruction and Visualization), to provide functionalities of easy metabolic network drafting and editing, amenable network visualization for experimental data integration and flux balance analysis tools for simulation studies. GEMSiRV comes with downloadable, ready-to-use public-domain metabolic models, reference metabolite/reaction databases and metabolic network maps, all of which can be input into GEMSiRV as the starting materials for network construction or simulation analyses. Furthermore, all of the GEMSiRV-generated metabolic models and analysis results, including projects in progress, can be easily exchanged in the research community. GEMSiRV is a powerful integrative resource that may facilitate the development of systems biology studies.

AVAILABILITY

The software is freely available on the web at http://sb.nhri.org.tw/GEMSiRV.