SUMMARY: COnstraint-Based Reconstruction and Analysis of genome-scale metabolic models has become a widely used tool to understand metabolic network behavior at a large scale. However, existing reconstruction tools lack functionalities to address modellers' common objective to study metabolic networks on the pathway level. Thus, we developed CobraMod-a Python package for pathway-centric modification and extension of genome-scale metabolic networks. CobraMod can integrate data from various metabolic pathway databases as well as user-curated information. Our tool tests newly added metabolites, reactions and pathways against multiple curation criteria, suggests manual curation steps and provides the user with records of changes to ensure high quality metabolic reconstructions. CobraMod uses the visualization tool Escher for pathway representation and offers simple customization options for comparison of pathways and flux distributions. Our package enables coherent and reproducible workflows as it can be seamlessly integrated with COBRApy and Escher. AVAILABILITY AND IMPLEMENTATION: The source code can be found at https://github.com/Toepfer-Lab/cobramod/ and can be installed with pip. The documentation including tutorials is available at https://cobramod.readthedocs.io/.

MOTIVATION: Understanding the rerouting of metabolic reaction fluxes upon perturbations has the potential to link changes in molecular state of a cellular system to alteration of growth. Yet, differential flux profiling on a genome-scale level remains one of the biggest challenges in systems biology. This is particularly relevant in plants, for which fluxes in autotrophic growth necessitate time-consuming instationary labeling experiments and costly computations, feasible for small-scale networks. RESULTS: Here we present a computationally and experimentally facile approach, termed iReMet-Flux, which integrates relative metabolomics data in a metabolic model to predict differential fluxes at a genome-scale level. Our approach and its variants complement the flux estimation methods based on radioactive tracer labeling. We employ iReMet-Flux with publically available metabolic profiles to predict reactions and pathways with altered fluxes in photo-autotrophically grown Arabidopsis and four photorespiratory mutants undergoing high-to-low CO2 acclimation. We also provide predictions about reactions and pathways which are most strongly regulated in the investigated experiments. The robustness and variability analyses, tailored to the formulation of iReMet-Flux, demonstrate that the findings provide biologically relevant information that is validated with external measurements of net CO2 exchange and biomass production. Therefore, iReMet-Flux paves the wave for mechanistic dissection of the interplay between pathways of primary and secondary metabolisms at a genome-scale. AVAILABILITY AND IMPLEMENTATION: The source code is available from the authors upon request. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Constraint-based approaches have been used for integrating data in large-scale metabolic networks to obtain insights into metabolism of various organisms. Due to the underlying steady-state assumption, these approaches are usually not suited for making predictions about metabolite levels. Here, we ask whether we can make inferences about the variability of metabolite levels from a constraint-based analysis based on the integration of transcriptomics data. To this end, we analyze time-resolved transcriptomics and metabolomics data from Arabidopsis thaliana under a set of eight different light and temperature conditions. In a previous study, the gene expression data have already been integrated in a genome-scale metabolic network to predict pathways, termed modulators and sustainers, which are differentially regulated with respect to a biochemically meaningful data-driven null model. Here, we present a follow-up analysis which bridges the gap between flux- and metabolite-centric methods. One of our main findings demonstrates that under certain environmental conditions, the levels of metabolites acting as substrates in modulators or sustainers show significantly lower temporal variations with respect to the remaining measured metabolites. This observation is discussed within the context of a systems-view of plasticity and robustness of metabolite contents and pathway fluxes. Our study paves the way for investigating the existence of similar principles in other species for which both genome-scale networks and high-throughput metabolomics data of high quality are becoming increasingly available.