MUFINS: multi-formalism interaction network simulator.

Systems Biology has established numerous approaches for mechanistic modeling of molecular networks in the cell and a legacy of models. The current frontier is the integration of models expressed in different formalisms to address the multi-scale biological system organization challenge. We present MUFINS (MUlti-Formalism Interaction Network Simulator) software, implementing a unique set of approaches for multi-formalism simulation of interaction networks. We extend the constraint-based modeling (CBM) framework by incorporation of linear inhibition constraints, enabling for the first time linear modeling of networks simultaneously describing gene regulation, signaling and whole-cell metabolism at steady state. We present a use case where a logical hypergraph model of a regulatory network is expressed by linear constraints and integrated with a Genome-Scale Metabolic Network (GSMN) of mouse macrophage. We experimentally validate predictions, demonstrating application of our software in an iterative cycle of hypothesis generation, validation and model refinement. MUFINS incorporates an extended version of our Quasi-Steady State Petri Net approach to integrate dynamic models with CBM, which we demonstrate through a dynamic model of cortisol signaling integrated with the human Recon2 GSMN and a model of nutrient dynamics in physiological compartments. Finally, we implement a number of methods for deriving metabolic states from ~omics data, including our new variant of the iMAT congruency approach. We compare our approach with iMAT through the analysis of 262 individual tumor transcriptomes, recovering features of metabolic reprogramming in cancer. The software provides graphics user interface with network visualization, which facilitates use by researchers who are not experienced in coding and mathematical modeling environments.

MultiMetEval: comparative and multi-objective analysis of genome-scale metabolic models.

Comparative metabolic modelling is emerging as a novel field, supported by the development of reliable and standardized approaches for constructing genome-scale metabolic models in high throughput. New software solutions are needed to allow efficient comparative analysis of multiple models in the context of multiple cellular objectives. Here, we present the user-friendly software framework Multi-Metabolic Evaluator (MultiMetEval), built upon SurreyFBA, which allows the user to compose collections of metabolic models that together can be subjected to flux balance analysis. Additionally, MultiMetEval implements functionalities for multi-objective analysis by calculating the Pareto front between two cellular objectives. Using a previously generated dataset of 38 actinobacterial genome-scale metabolic models, we show how these approaches can lead to exciting novel insights. Firstly, after incorporating several pathways for the biosynthesis of natural products into each of these models, comparative flux balance analysis predicted that species like Streptomyces that harbour the highest diversity of secondary metabolite biosynthetic gene clusters in their genomes do not necessarily have the metabolic network topology most suitable for compound overproduction. Secondly, multi-objective analysis of biomass production and natural product biosynthesis in these actinobacteria shows that the well-studied occurrence of discrete metabolic switches during the change of cellular objectives is inherent to their metabolic network architecture. Comparative and multi-objective modelling can lead to insights that could not be obtained by normal flux balance analyses. MultiMetEval provides a powerful platform that makes these analyses straightforward for biologists. Sources and binaries of MultiMetEval are freely available from https://github.com/PiotrZakrzewski/MetEval/downloads.

SurreyFBA: a command line tool and graphics user interface for constraint-based modeling of genome-scale metabolic reaction networks.

UNLABELLED: Constraint-based modeling of genome-scale metabolic networks has been successfully used in numerous applications such as prediction of gene essentiality and metabolic engineering. We present SurreyFBA, which provides constraint-based simulations and network map visualization in a free, stand-alone software. In addition to basic simulation protocols, the tool also implements the analysis of minimal substrate and product sets, which is useful for metabolic engineering and prediction of nutritional requirements in complex in vivo environments, but not available in other commonly used programs. The SurreyFBA is based on a command line interface to the GLPK solver distributed as binary and source code for the three major operating systems. The command line tool, implemented in C++, is easily executed within scripting languages used in the bioinformatics community and provides efficient implementation of tasks requiring iterative calls to the linear programming solver. SurreyFBA includes JyMet, a graphics user interface allowing spreadsheet-based model presentation, visualization of numerical results on metabolic networks represented in the Petri net convention, as well as in charts and plots. AVAILABILITY: SurreyFBA is distributed under GNU GPL license and available from http://sysbio3.fhms.surrey.ac.uk/SurreyFBA.zip.

Building and analysing genome-scale metabolic models.

Reconstructing a model of the metabolic network of an organism from its annotated genome sequence would seem, at first sight, to be one of the most straightforward tasks in functional genomics, even if the various data sources required were never designed with this application in mind. The number of genome-scale metabolic models is, however, lagging far behind the number of sequenced genomes and is likely to continue to do so unless the model-building process can be accelerated. Two aspects that could usefully be improved are the ability to find the sources of error in a nascent model rapidly, and the generation of tenable hypotheses concerning solutions that would improve a model. We will illustrate these issues with approaches we have developed in the course of building metabolic models of Streptococcus agalactiae and Arabidopsis thaliana.

Integration of metabolic databases for the reconstruction of genome-scale metabolic networks.

BACKGROUND: Genome-scale metabolic reconstructions have been recognised as a valuable tool for a variety of applications ranging from metabolic engineering to evolutionary studies. However, the reconstruction of such networks remains an arduous process requiring a high level of human intervention. This process is further complicated by occurrences of missing or conflicting information and the absence of common annotation standards between different data sources. RESULTS: In this article, we report a semi-automated methodology aimed at streamlining the process of metabolic network reconstruction by enabling the integration of different genome-wide databases of metabolic reactions. We present results obtained by applying this methodology to the metabolic network of the plant Arabidopsis thaliana. A systematic comparison of compounds and reactions between two genome-wide databases allowed us to obtain a high-quality core consensus reconstruction, which was validated for stoichiometric consistency. A lower level of consensus led to a larger reconstruction, which has a lower quality standard but provides a baseline for further manual curation. CONCLUSION: This semi-automated methodology may be applied to other organisms and help to streamline the process of genome-scale network reconstruction in order to accelerate the transfer of such models to applications.

Detection of stoichiometric inconsistencies in biomolecular models.

MOTIVATION: Metabolic modelling provides a mathematically rigorous basis for system-level analysis of biochemical networks. However, the growing sizes of metabolic models can lead to serious problems in their construction and validation. In this work, we describe a relatively poorly investigated type of modelling error, called stoichiometric inconsistencies. These errors are caused by incorrect definitions of reaction stoichiometries and result in conflicts between two fundamental physical constraints to be satisfied by any valid metabolic model: positivity of molecular masses of all metabolites and mass conservation in all interconversions. RESULTS: We introduce formal definitions of stoichiometric inconsistencies, inconsistent net stoichiometries, elementary leakage modes and other important fundamental properties of incorrectly defined biomolecular networks. Algorithms are described for the verification of stoichiometric consistency of a model, detection of unconserved metabolites and inconsistent minimal net stoichiometries. The usefulness of these algorithms for effective resolving of inconsistencies and for detection of input errors is demonstrated on a published genome-scale metabolic model of Saccharomyces cerevisiae and one of Streptococcus agalactiae constructed using the KEGG database. AVAILABILITY: http://mudshark.brookes.ac.uk/index.php/Albert_Gevorgyan.