Skip to content

Advertisement

You're viewing the new version of our site. Please leave us feedback.

Learn more

BMC Systems Biology

Open Access

MetaNET - a web-accessible interactive platform for biological metabolic network analysis

  • Pankaj Narang1, 2,
  • Shawez Khan1, 2,
  • Anmol Jaywant Hemrom1,
  • Andrew Michael Lynn1, 2Email author and
  • Open Source Drug Discovery Consortium2
BMC Systems Biology20148:130

https://doi.org/10.1186/s12918-014-0130-2

Received: 16 July 2014

Accepted: 12 November 2014

Published: 5 December 2014

Abstract

Background

Metabolic reactions have been extensively studied and compiled over the last century. These have provided a theoretical base to implement models, simulations of which are used to identify drug targets and optimize metabolic throughput at a systemic level. While tools for the perturbation of metabolic networks are available, their applications are limited and restricted as they require varied dependencies and often a commercial platform for full functionality. We have developed MetaNET, an open source user-friendly platform-independent and web-accessible resource consisting of several pre-defined workflows for metabolic network analysis.

Result

MetaNET is a web-accessible platform that incorporates a range of functions which can be combined to produce different simulations related to metabolic networks. These include (i) optimization of an objective function for wild type strain, gene/catalyst/reaction knock-out/knock-down analysis using flux balance analysis. (ii) flux variability analysis (iii) chemical species participation (iv) cycles and extreme paths identification and (v) choke point reaction analysis to facilitate identification of potential drug targets. The platform is built using custom scripts along with the open-source Galaxy workflow and Systems Biology Research Tool as components. Pre-defined workflows are available for common processes, and an exhaustive list of over 50 functions are provided for user defined workflows.

Conclusion

MetaNET, available at http://metanet.osdd.net, provides a user-friendly rich interface allowing the analysis of genome-scale metabolic networks under various genetic and environmental conditions. The framework permits the storage of previous results, the ability to repeat analysis and share results with other users over the internet as well as run different tools simultaneously using pre-defined workflows, and user-created custom workflows.

Keywords

Flux balance analysisMetabolic networkSystems biologyin silico gene knock-outPerturbation analysis

Background

Systems biology investigates the components of complex biological networks and can pinpoint drug targets through a combination of experimental and computational research [1]. Over the past few years, various approaches have been actively developed which attempt to provide a systems level analysis of these networks. Most of these approaches use system dynamics while others are based on a static representation of the networks which together form the core of systems biology [2]. Flux balance analysis (FBA) is a commonly used method in the field of systems biology for the quantitative simulation of metabolic networks using steady state stoichiometric models [3]. It is based on a constraint-based mathematical approach to calculate systemic phenotypes in the form of reaction fluxes. These fluxes are then used to interpret the metabolic capabilities of the system. The method’s widespread use is largely due to its independence from kinetic parameters, requiring only the stoichiometry of metabolic reactions.

Fundamentally, FBA requires four steps: (i) reconstruction of the metabolic network (ii) obtaining the stoichiometry matrix of the metabolic network (iii) defining the objective function and biochemical relevant constraints and (iv) optimization [4]. Reconstruction involves generating a network by identifying, compartmentalizing and interconnecting the components i.e. genes, proteins, reactions and metabolites involved in various activities of the network. The reconstructed network is then converted into a matrix Sm×n where m rows represent metabolites, n columns represent reactions and each element of the matrix represents the stoichiometric coefficient of the metabolite in the corresponding reaction. The dynamic mass balance of the metabolic system is
d x d t = S . v
where ν 1 = [v1, v2, ..... v n ] is the vector of the unsolved fluxes and × 1 is the vector representing metabolite concentration. At steady state, the flux distribution is defined by a system of linear equations as
d x d t = S . v = 0
In general, for a biological metabolic system the number of reactions is more than the number of metabolites leading to an under-determined solution of steady state. Therefore, additional biologically relevant constraints are imposed on fluxes of reactions in order to find substantive solutions of the system.
0 v i
b i
where ν i represents internal fluxes and bi signifies the exchange fluxes in the system. Next, a reaction of interest is defined as the objective function which formulate FBA as the following linear programming problem.
maximize C T v s . t . S . v = 0

where c represents a vector with the reaction of interest set to one, while all other reactions are set to zero [5],[6] . Examples of objective functions are maximization or minimization of biomass, ATP production or metabolite production. FBA has been successfully applied for predicting growth and metabolic by-product secretion in E.coli [7].

Over the last decade, FBA and other methods for metabolic network analysis have been implemented through various resources [8],[9]. The most prominent resources include the COBRA toolbox [10], Metatool [11], FluxAnalyzer [12], MetaFluxNet [13], SNA [14] and CellNetAnalyzer [15], a successor of FluxAnalyzer. COBRA, Metatool and CellNetAnalyzer require the proprietary MATLAB environment while SNA requires MATHEMATICA for their functionality. In addition to these commercial resources, many open source applications were also developed such as YANAsquare implemented in Java [16], Linear Inverse Model [17], sybil [18] and BiGGR [19] that run within the R environment [20]. Another open source package is the Java-based System Biology Research Tool (SBRT) that has only a basic graphical user interface [21]. All of these packages are of great value for researchers, but are stand-alone applications which often require installation of external dependencies and regular updates with time. On the other hand, web-based applications are platform-independent and require only internet access to construct and analyze networks. In the context of web-based applications, FAME [22] and MicrobesFlux [23] provide facilities to create, edit and analyze genome scale model for different micro-organism obtained from the KEGG database [24]. Flux-P is another efficient web-based tool for automating and standardizing 13C-based metabolic flux analysis, using the Bio-jETI workflow framework [25]. Webcoli [26] provides a framework for reconstructing a genome-scale metabolic model of E.coli. MetaNetX [27] is a valuable tool which offers analysis of several metabolic networks from BiGG [28] and MetRxn [29] in its repository. Existing web-based tools contain only pre-defined applications to simulate network properties. We have introduced a web accessible platform MetaNET, which is open source, user friendly and promises to simulate complex large metabolic networks without the requirement of typical hardware or software configuration at the user’s site. In addition to being web-accessible and platform independent, MetaNET offers a workflow editor to connect various tools together, along with a data-library to enable the sharing of data-sets. The framework contains a number of primitive tools which can be used stand-alone or combined to create higher level solutions by the user. This extensibility is showcased by the development of sample workflows that cover common metabolic simulations. A comparison of basic functionality of MetaNET with other existing stand-alone and web-based applications is provided as Table 1.
Table 1

Features comparison of MetaNET with similar software packages

Features

Installation requirement

External dependency

GUI

Programming knowledge

Workflow support

Functions

Tools

FBA

FVA

Choke points analysis

Chemical species participation

MetaNET

N

N

Y

N

Y

Y

Y

Y

Y

COBRA

Y

Y

N

N

N

Y

Y

N

N

FAME

N

N

Y

N

N

Y

Y

N

N

PathwayAnalyzer

Y

N

N

N

N

Y

N

N

N

Webcoli

N

N

Y

N

N

Y

N

N

N

CellNetAnalyzer

Y

Y

Y

N

N

Y

N

N

N

SBRT

N

N

N

Y

N

Y

Y

N

N

YANA Square

Y

N

Y

N

N

Y

N

N

N

Microbes flux

N

N

Y

N

N

Y

N

N

N

Metatool

Y

Y

N

Y

N

Y

N

N

N

MetanetX

N

N

Y

N

N

Y

Y

Y

N

Implementation

MetaNET has been implemented with different functionalities of the SBRT under a GNU/Linux operating system. Although, many processes to analyze biological networks are available in SBRT, it was designed with a basic non-interactive GUI which hinders its widespread use within the research community. More importantly, users of SBRT have to create various input files manually or through programs, which are time-consuming and prone to error while handling large metabolic networks. To overcome these limitations, we extended SBRT with a user friendly platform for metabolic network studies by integrating SBRT functions with the open source platform Galaxy, available over the internet [30]. Galaxy provides a web accessible platform to integrate different command line tools to make it more interactive. Along with processes of SBRT, we have also implemented some new tools using Perl and Java languages to provide a wide range of functionalities within a single platform (Figure 1). The R Language for statistical computing is also integrated with the system, and used to provide graphical output.
Figure 1

A schematic diagram representing the different components of MetaNET. The client–server architecture of MetaNET was developed using Galaxy framework and System Biology Research Tool (SBRT) as its main components. The functions of SBRT were called using Java wrappers and integrated with Galaxy using XML files. Additionally, the functions not present in SBRT were implemented using Perl and integrated using XML files.

(i) Installation of tools

All the three applications: Galaxy, SBRT and R were installed on a development machine running CentOS 6.5 with the prerequisites– python version 2.6.6 & Java Development kit. Galaxy was downloaded from http://wiki.galaxyproject.org/Admin/Get%20Galaxy & installed intentionally as a normal user – “galaxy”, to permit easy migration and for reasons of security. The application runs using HTTP with a user-specified port, and does not need root privileges. This results in the easy migration by porting the entire user space to another system and as root privileges are not required, it also secures the base operating system. SBRT was downloaded from and installed with the same user i.e. “galaxy”. Sybil [18] was used to find exchange reactions from SBML file. R is available as an installable package on most standard GNU/Linux operating systems, and was added by updating the package on the base operating system.

(ii) Application development

The application was implemented using XML files, Java and Perl wrappers at the following levels, as shown in Figure 2.
  1. (a)

    XML files were written for each integrated SBRT and in-house tool which include all the information required to execute the tool. Using these files, Galaxy automatically generates the user interface for user to select input files & parameters required by the tool. On execution, Galaxy calls the respective Perl/Java wrapper and runs the tool. Table 2 contains a list of available tools and their descriptions.

     
Table 2

Description of all tools available with MetaNET

Tool sub-category

Tool name

Tool description

Source

Upload Data

Upload file

To upload file on MetaNET.

Galaxy

Tool category: Format Converter

Import to MetaNET

BiGGSBML Reader

Converts BIGG SBML Format to reactions file.

SBRT

PalssonSBML Reader

Converts Palsson SBML Format to reactions file.

SBRT

Metatool Input File Reader

Converts input files of Metatool into reactions File.

SBRT

Export to MetaNET

Metatool Input File Writer

Converts reaction Files into input files of Metatool.

SBRT

Tool category: Network Information Retriever

Fetch IDs

Fetch Objective function

Retrieves reaction ID of objective function from SBML file.

In-house

Fetch Gene IDs

Retrieves list of all gene Ids from SBML file.

In-house

Fetch Reaction IDs

Retrieves list of reaction Ids from reactions file.

SBRT

Fetch Catalyst IDs

Retrieves list of all catalysts from SBML file.

In-house

Fetch Chemical Species (Metabolites)

Retrieves list of chemical species Ids from reactions file.

SBRT

Fetch Reaction attributes

Fetch Flux Constraints

Retrieves flux bounds from SBML file

In-house

Fetch Reactions-Genes Associations

Retrieves reaction-gene reactions associations from SBML file.

SBRT

Fetch Reactions-Catalysts Associations

Retrieves reaction-catalyst associations from SBML file.

SBRT

Fetch Other Info

Generate Stochiometric Matrix

Generates stoichiometry matrix.

SBRT

Generate Linear Equations

Generates system of linear equations.

SBRT

Fetch Reactions

Retrieves reactions without Ids.

In-house

Fetch Reactants

Retrieves list of all reactants from the given reactions file.

In-house

Fetch Products

Retrieves list of all products from the given reactions file.

In-house

Fetch Equivalent Reactions

Retrieves names of all stoichiometric equivalent reactions.

SBRT

Fetch Flux Vectors

Retrieves flux vectors from SBML file.

SBRT

Tool category: Flux Balance Analysis

Single Objective Optimizer

Single Objective Optimizer

Computes optimal value of fluxes in a stoichiometric network using multiple constraints.

SBRT

Multiple Objectives Optimizer

Multiple Objectives Optimizer

Computes optimal value of multiple objectives in a stoichiometric network using multiple constraints.

SBRT

Tool category: Perturbation Analysis

Single Entity Knock-outs

Reaction knock-out Optimizer

Computes effects of deleting sets of reactions from a stoichiometric network.

SBRT

Catalyst knock-out Optimizer

Computes effects of deleting sets of catalysts from a stoichiometric network.

SBRT

Gene Knock-out Optimizer

Computes effects of deleting sets of genes from a stoichiometric network.

SBRT

Pairwise Entity Knock-outs

Pairwise Genes Knock-out

Computes effects of deleting pairs of genes from a chemical reaction network for multi-target drug identification.

SBRT

Pairwise Catalysts knock-out

Computes effects of deleting pairs of catalysts from a chemical reaction network for multi-target drug identification.

SBRT

Pairwise Reactions knock-out

Computes effects of deleting pairs of reactions from a chemical reaction network for multi-target drug identification.

SBRT

Tool category: Essentiality Test

Find Essential Reactions

Essential Reactions Reporter

Finds essential reactions of a network.

In-house

Find Essential Catalysts

Essential Catalysts Reporter

Finds essential catalysts of a network.

In-house

Find Essential Genes

Essential Genes Reporter

Finds essential genes of a network.

In-house

Tool category: Flux Variability Analysis

Constraint Variations

Create Constraint Variation file

Creates constraint variation file for Flux Variability Analysis.

In-house

Flux Variability Optimizer

Computes optimal value of a single objective function for multiple set of flux constraints.

SBRT

Find Flux caps

Flux Cap Identification

Creates caps for each unbounded flux in a stoichiometric network.

SBRT

Tool category: Other Utilities

File Operations

Remove beginning

Removes constraints or objectives from a file.

Galaxy

Select First

Selects first n objectives or constraints from a file.

Galaxy

Select Last

Selects last n constraints or objectives from a file.

Galaxy

Add constraint

Adds constraint in a constraint file.

Galaxy

Replace constraint

Replaces constraint of particular reaction.

In-house

Make pairs

Makes all possible pairs of genes or catalysts or reactions IDs.

In-house

Tool category: Flux Comparisons

Flux Distribution Comparison

Flux Distribution Comparison

Compares flux distributions for equality within a given tolerance.

SBRT

Flux Intervals Comparison

Flux Interval Comparison

compares intervals for equality within a given tolerance.

SBRT

Tool category: Random Generator

Random Generator

Random Objectives Generator

Generates random objective functions.

SBRT

Random Constraints Generator

Generates random constraints.

SBRT

Tool category: Flux Distribution Plots

Flux Optimizer Plot

Flux Optimization Plotter

Plots result of FBA optimization.

In-house

Essentiality Plots

Reactions Essentiality Plotter

Plots result of deleting sets of reactions versus objective function.

In-house

Genes Essentiality Plotter

Plots result of deleting sets of genes versus objective function.

In-house

Catalysts Essentiality Plotter

Plots result of deleting sets of catalysts versus objective function.

In-house

Knock-Out Plots

Reactions knock-out Plotter

Plots the results of deleting sets of reactions or genes in a stoichiometric network.

In-house

Genes knock-out Plotter

Plots the results of deleting sets of genes in a stoichiometric network.

In-house

Catalysts knock-out Plotter

Plots the results of deleting sets of catalysts in a stoichiometric network.

In-house

Tool category: Chemical Species participation

Grouping Chemical Species

Chemical Species Participation as Reactants

Groups chemical reactions based on the given chemical species as reactants.

SBRT

Chemical Species Participation as Product

Groups chemical reactions based on the given chemical species as products.

SBRT

Single Species participation

Single Chemical species Participation

Finds chemical reactions containing given chemical species.

In-house

Chokepoint Analysis

Find Chokepoint Reactions

finds chokepoint reactions of biological network.

In-house

Tool category: Network Reducer

Reversible Reactions Breaker

Reversible Reactions Breaker

Breaks reversible reactions into pairs of irreversible “forward” and “reverse” reactions.

SBRT

Redundant Reactions Remover

Redundant Reactions Remover

Removes redundant reactions from stoichiometric networks.

SBRT

WW Network Reducer

WW Network Reducer

Reduces size of stoichiometric networks for the purpose of identifying the cycles they contain.

SBRT

MS Network Reducer

MS Network Reducer

Reduces size of stoichiometric networks for the purpose of identifying the cycles they contain.

SBRT

Tool category: Cycles/Paths Identifier

SLP Cycles

SLP Cycle Identifier

Identifies cycles in stoichiometric networks.

SBRT

Extreme Paths

Extreme Path Identifier

Identifies extreme currents in stoichiometric networks.

SBRT

  1. (b)

    For SBRT based methods, wrappers were written in Java to call functions of SBRT. These wrappers use SBRT as an application programming interface (API) by calling its respective functions.

     
  2. (c)

    For in-house tools, a second set of wrappers were written in Perl that provide missing functionality in SBRT. Additional Perl scripts were also written that take input from Galaxy and generate plots in MetaNET using R.

     
Figure 2

Data flow within MetaNET. Galaxy creates interfaces based on tool configuaration XML files for input of various parameters. User provides these parameters, input files and executes the tool through client machine. In turn, Galaxy passes these files to Java or Perl wrappers to run SBRT function or R scripts and returns results back to the client machine.

(iii) Creation and sharing of custom workflows

Custom workflows were designed by combining different tools using workflow editor. The framework can filter the coupling of tools for permitting only those cases where the output from one process can be used as input of the other. Sample workflows, created to perform common metabolic simulation tasks, are listed along with a description of their functionality in Table 3.
Table 3

Description of workflows implemented in MetaNET

MetaNET published workflows

Description

No. of steps involved

Automated Flux Balance Analysis Workflow

Used to perform FBA of metabolic network using default objective function

6

Customized Flux Balance Analysis Workflow

Used to perform FBA of metabolic network using user-defined objective function

5

Gene Essentiality Workflow

Used to investigate lethal genes

8

Reaction Essentiality Workflow

Used to investigate lethal reactions

7

Catalyst Essentiality workflow

Used to investigate lethal catalysts

8

Pairwise Genes knock-out Workflow

Used to knock-out genes pairwise

8

Pairwise Reactions knock-out Workflow

Used to knock-out reactions pairwise

8

Pairwise Catalysts knock-out workflow

Used to knock-out catalysts pairwise

8

Results and discussion

MetaNET currently has over 50 tools grouped into different categories (Table 2; MetaNET tool panel). Besides basic features for uploading data, format conversion and data editing, tools are available for FBA optimization, perturbation analysis and others that help to predict the cell phenotype under various genetic or environmental conditions. In addition, we have implemented a wide range of novel functionalities using in-house scripts, not available with existing web-based applications for metabolic simulations. For example, a systematic gene knock-out to identify lethal genes is a common function during analysis of biological networks, iteratively replacing the flux constraint of each reaction to zero and calculating the effective biomass. With MetaNET, users can additionally perform a pairwise knock-out of genes/catalysts/reactions for multi-target drug identification, which is otherwise a laborious task experimentally [31]. This involves switching off the reaction constraints by setting them to zero before calculating the effective biomass. The reactions, constraints and gene identifiers all extracted from the SBML file by other MetaNET tools (listed under category “Network Information Retriever”) or with the user supplied list of genes provides the input file for this function “Pairwise Entity Knock-outs” listed as a separate category on the tool panel. The results of various simulations can also be visualized graphically for better interpretation using different tools listed under category “Flux Distribution Plots”.

Model uploading and editing

MetaNET allows users to import the model network file in different formats (SBML [32], BiGGSBML [28], PalssonSBML, Metatool Format [11], or flat file) via three methods: uploading the user’s reconstructed model from a local machine, using existing models in the MetaNET data repositories or by specifying a URL to a network file from a standard database [26], ensuring the model is available for constraint based analysis. Once the model is uploaded, tools are available under “File operations” for adding, removing or replacing flux constraints. A more detailed description for various analysis tools with required input files and their formats, is available in the user manual with MetaNET.

Although MetaNET can be used in various types of analysis, it can be broadly classified into two major categories: Simulation studies and Topological analysis (Figure 3).
Figure 3

Tool categorization within MetaNET.

(i) Simulation studies

MetaNET provides various tools for simulation using FBA by applying different genetic or environmental constraints. To emphasize the core features of MetaNET, three case studies were performed with the publicly available genome-scale metabolic in silico model of Escherichia coli, (iJR904 GSM/GPR), from Reed et. al. [33]. The system consists of 761 metabolites (including external metabolites) and 1075 biochemical reactions (including transport reactions). The functionality of MetaNET has been validated in following case studies:

Simulation of Escherichia colimodel under aerobic conditions

To predict optimal biomass growth under aerobic conditions, the network was optimized through the biomass reaction by flux balance analysis using the “Single Objective optimizer” tool for the wild type strain (Figure 4A). The growth rate of biomass under aerobic conditions was found to be 0.92 gDW/gDW/hr while keeping glucose uptake rate fixed at 10 mmol/gDW/hr (Figure 4B). This result is in agreement with earlier published reports [34]. The flux distribution of all reactions after optimization was plotted using the “FBA optimization plotter” (Figure 4C).
Figure 4

FBA simulation of in silico model of Escherichia coli. (A) View of “Single Objective Optimizer” for FBA of in silico Escherichia coli model; (B) Flux distribution using Biomass reaction as objective function for wildtype; (C) Graphical representation of flux distribution using “FBA optimization plotter”.

Simulation of Escherichia colimodel under anaerobic conditions shows reduced growth

To predict optimal biomass growth under anaerobic conditions, the flux of the reaction which corresponds to exchange of oxygen from the extracellular to cytosolic compartment was constrained to zero using the “Replace constraint” tool and the network was again optimized using the “Single Objective optimizer” tool (Figure 5A). The growth rate of biomass under anaerobic conditions was found to be 0.22 gDW/gDW/hr (Figure 5B), reduction of 77% as compared to aerobic condition and similar to experimental determined values [35].
Figure 5

FBA simulation of in silico model of Escherichia coli. (A) Anaerobic conditions were created by reducing the constraint of reaction corresponding to exchange of oxygen to zero using “Replace constraint tool”; (B) FBA result under anaerobic conditions; (C) Results of reaction essentiality analysis.

Identification of essential reactions through perturbation analysis

To identify the critical reactions of the network, all 1075 reactions were iteratively knocked-out and its effect on biomass was calculated using the “Reaction knock-out optimizer” tool (listed under the category “perturbation analysis”). Using “Essential Reactions Reporter”, we found 135 reactions to be essential as their deletion reduce the flux of biomass significantly (Figure 5C).

(ii) Topological analysis

This category consists of tools designed to determine robustness and modularity based on the topology of the metabolic network. “SLP cyclic identifier” and “Extreme path identifier” are Java-based functions available within SBRT that have been have been integrated in MetaNET. Additionally, many tools are available in MetaNET which have been implemented using in-house Perl scripts. For instance, choke points are reactions present in the metabolic network which uniquely consume a particular substrate or produce a particular product [36],[37]. Any enzyme which catalyzes such a reaction can be a used as drug target because blocking such enzyme results in either accumulation of particular substrate or removal of a particular product essential for survival of cell. If an enzyme catalyzes at least one choke point reaction, it can be classified as potential drug target. MetaNET provides a tool “Choke point analysis” for identification of choke point reactions in a given metabolic network for proposing drug targets. Further, MetaNET also provides a tool “Chemical species participation” to identify reactions involving metabolites of interest as reactants or products.

Workflows support

One of the most powerful feature of MetaNET is to utilize the workflow support provided by the galaxy framework to couple tools into higher order solutions. Keeping this in mind, we designed different tools in such a way that these can be joined to create workflows in the MetaNET. Users can create custom workflows for analysis by linking appropriate tools of MetaNET using the workflow editor. The tools can only be joined if the output of one is compatible with the input of other. We have implemented eight different workflows which are generally used during metabolic network simulation (Table 3). These workflows have been publicly shared with all the users of MetaNET. Users can import these workflows to their own user-space from the “published workflows” section and extend it according to their requirements. As examples of workflows, we describe here the “Automated FBA”, “Gene Essentiality” and “Pairwise Reaction Knock-outs” workflows.

Automated flux balance analysis workflow

This is the simplest workflow which reads the SBML file and plots the flux distribution after optimization (Figure 6). It calls the tools “PalssonSBML reader”, “Fetch flux constraints” and ”Fetch Objective function” which retrieve the reactions, corresponding constraints and objective function from the SBML file. The resulting output is passed to the “Single objective optimizer” tool which optimizes the objective function and provides the output as fluxes of all reactions. The final tool which processes the output from the earlier workflow modules is “FBA optimization plotter” that plots the flux distribution. The “PalssonSBML reader” and “Single objective optimizer” are functions of SBRT while other tools of this workflow were developed in-house.
Figure 6

Flux balance analysis workflow. (A) stepwise representation (B) workflow representation.

Gene essentiality workflow

Genes are said to be lethal if their deletion stops the growth of the organism. This workflow finds lethal genes in the given metabolic network by iteratively knocking down each gene and optimizing the metabolic network of the mutant strain (Figure 7). This workflow was designed by joining “Fetch Gene IDs”, “Fetch flux constraints”, “Fetch objective function”, “PalssonSBML Reader” and “Fetch Reaction-Gene associations” tools which parse gene identifiers, reactions constraints, objective function, reactions and gene associated with each reaction respectively from the input SBML file. The resulting output is used by the tool “Gene Knock-out optimizer” which deletes all the genes one by one and optimizes the objective function. The final tool which processes the output from the earlier workflow module is the “Essential Gene Reporter” that selects genes which reduce biomass by a user provided percentage, and is a Perl script written in-house. As an example, if a user wants a list of genes, which when knocked out, reduce the biomass of the organism by more than 90%, one can either use this shared workflow or run the tools sequentially.
Figure 7

Gene essentiality workflow. Gene Essentiality workflow was designed by connecting various tools (shown in boxes) using workflow editor. Arrows represent the direction of flow of data from one tool to another tool.

Pairwise reactions knock-outs workflow

Pairwise reactions deletion experiments require considerable time and effort to perform on a large scale. This workflow is designed to rapidly predict the biomass flux by deleting a pair of reactions simultaneously. This workflow may predict some interesting pairs of reactions whose simultaneous deletion reduce the biomass, which otherwise are found to be non-essential during single-reaction deletion analysis. It calls the tools “Fetch Reaction IDs”, “Fetch objective function”, ”“Fetch flux constraints” and “PalssonSBML Reader” which parse the reaction identifiers, objective function, reactions constraints and reactions respectively from input SBML file. As the numbers of possible pairs is very large, it first finds the choke point reactions using “Find Choke Point Reactions” tool. The resulting output i.e. reaction Identifiers of choke point reactions is used by the tool “Make Pairs” tool implemented in Perl which makes all the possible pairs of reactions. Finally, another tool “Reaction knock-out optimizer” deletes the pairs sequentially and optimizes the objective function.

Conclusions

MetaNET is a user-friendly, extensible and platform independent framework for metabolic network analysis freely available at http://metanet.osdd.net. The framework is built with a set of tools for data management including data upload/download, format conversion, file operations and data extraction capabilities from SBML files, optimizing network using FBA, flux variability analysis [38], perturbation analysis via gene/reaction/catalyst knock-out/knock-down (single or pairwise) and visualizing results. The platform is developed by integrating the Systems Biology Research Tool with Galaxy as a proof of concept along with in-house scripts for various other functions. The rich list of tools can be interconnected through workflows to perform customized higher order simulations. End-users are encouraged to use published workflows with curated data models provided through the data libraries. Developers of data models and workflows are invited to join the mailing list set up for the purpose. The functionality of MetaNET can be further extended by including other simulation engines, and provisioning additional functions on demand from a user community, while keeping a consistent interface for models, data-sharing and functional use through shared workflows.

Availability and requirements

Project Name: MetaNET

Project home page: http://metanet.osdd.net/

Operating system: Platform Independent

Programming Language: Language independent

Other requirements: Internet

Any restrictions to use by non-academics: None

Abbreviations

FBA: 

Flux balance analysis

FVA: 

Flux variability analysis

SBML: 

System biology markup language

SBRT: 

System biology research tool

Declarations

Acknowledgements

PN acknowledges the Department of Biotechnology, Govt. of India for providing BINC Junior Research Fellowship. SK was supported by an Open Source Drug Discovery project Fellowship, and AJH a research fellowship from the Council of Scientific and Industrial Research, Govt. of India.

Authors’ Affiliations

(1)
School of Computational & Integrative Sciences, Jawaharlal Nehru University
(2)
The Open Source Drug Discovery (OSDD) Consortium, Council of Scientific and Industrial Research

References

  1. Kitano H: Computational systems biology. Nature. 2002, 420: 206-210. 10.1038/nature01254.View ArticlePubMedGoogle Scholar
  2. Feist AM, Palsson BO: The growing scope of applications of genome-scale metabolic reconstructions using Escherichia coli. Nat Biotechnol. 2008, 26: 659-667. 10.1038/nbt1401.PubMed CentralView ArticlePubMedGoogle Scholar
  3. Kauffman KJ, Prakash P, Edwards JS: Advances in flux balance analysis. Curr Opin Biotechnol. 2003, 14: 491-496. 10.1016/j.copbio.2003.08.001.View ArticlePubMedGoogle Scholar
  4. Lee JM, Gianchandani EP, Papin JA: Flux balance analysis in the era of metabolomics. Brief Bioinform. 2006, 7: 140-150. 10.1093/bib/bbl007.View ArticlePubMedGoogle Scholar
  5. Orth JD, Thiele I, Palsson BØ: What is flux balance analysis?. Nat Biotechnol. 2010, 28: 245-248. 10.1038/nbt.1614.PubMed CentralView ArticlePubMedGoogle Scholar
  6. Raman K, Chandra N: Flux balance analysis of biological systems: applications and challenges. Brief Bioinform. 2009, 10 (4): 435-449. 10.1093/bib/bbp011.View ArticlePubMedGoogle Scholar
  7. Varma A, Palsson BO: Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110. Appl Environ Microbiol. 1994, 60 (10): 3724-3731.PubMed CentralPubMedGoogle Scholar
  8. Dandekar T, Fieselmann A, Majeed S, Ahmed Z: Software applications toward quantitative metabolic flux analysis and modeling. Brief Bioinform 2012, doi:10.1093/bib/bbs065.,Google Scholar
  9. Lakshmanan M, Koh G, Chung BK, Lee DY: Software applications for flux balance analysis. Brief Bioinform. 2012, 15: 108-122. 10.1093/bib/bbs069.View ArticlePubMedGoogle Scholar
  10. Schellenberger J, Que R, Fleming RM, Thiele I, Orth JD, Feist AM, Zielinski DC, Bordbar A, Lewis NE, Rahmanian S, Kang J, Hyduke DR, Palsson BO: Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc. 2011, 6: 1290-1307. 10.1038/nprot.2011.308.PubMed CentralView ArticlePubMedGoogle Scholar
  11. Von KA, Schuster S: Metatool 5.0: fast and flexible elementary modes analysis. Bioinformatics. 2006, 22: 1930-1931. 10.1093/bioinformatics/btl267.View ArticleGoogle Scholar
  12. Klamt S, Stelling J, Ginkel M, Gilles ED: FluxAnalyzer: exploring structure, pathways, and flux distributions in metabolic networks on interactive flux maps. Bioinformatics. 2003, 19: 261-269. 10.1093/bioinformatics/19.2.261.View ArticlePubMedGoogle Scholar
  13. Lee DY, Yun H, Park S, Lee SY: MetaFluxNet: the management of metabolic reaction information and quantitative metabolic flux analysis. Bioinformatics. 2003, 19: 2144-2146. 10.1093/bioinformatics/btg271.View ArticlePubMedGoogle Scholar
  14. Urbanczik R: SNA–a toolbox for the stoichiometric analysis of metabolic networks. BMC Bioinformatics. 2006, 7: 129-10.1186/1471-2105-7-129.PubMed CentralView ArticlePubMedGoogle Scholar
  15. Klamt S, Saez-Rodriguez J, Gilles ED: Structural and functional analysis of cellular networks with Cell NetAnalyzer. BMC Syst Biol. 2007, 1: 2-10.1186/1752-0509-1-2.PubMed CentralView ArticlePubMedGoogle Scholar
  16. Schwarz R, Liang C, Kaleta C, Kuhnel M, Hoffmann E, Kuznetsov S, Hecker M, Griffiths G, Schuster S, Dandekar T: Integrated network reconstruction, visualization and analysis using YANAsquare. BMC Bioinformatics. 2007, 8: 313-10.1186/1471-2105-8-313.PubMed CentralView ArticlePubMedGoogle Scholar
  17. van Oevelen D, Van den Meersche K, Meysman FJR, Soetaert K, Middelburg JJ, Vezina AF: Quantifying food web flows using linear inverse models. Ecosystems. 2010, 13: 32-45. 10.1007/s10021-009-9297-6.View ArticleGoogle Scholar
  18. Gelius DG, Amer DA, Jonathan FC, Lercher M: sybil: efficient constrained based modelling in R. BMC Syst Biol. 2013, 7: 125-10.1186/1752-0509-7-125.View ArticleGoogle Scholar
  19. Gavai AK, Hettling H: BiGGR. Constraint based modeling in R using metabolic reconstruction databases. R package version1.2.0. 2013, [http://www.bioconductor.org/], [http://www.bioconductor.org/]Google Scholar
  20. R Core Team: R: a language and environment for statistical computing. In Vienna, Austria: R Foundation for Statistical Computing; 2014. URL ., [http://www.R-project.org/]
  21. Wright J, Wagner A: The systems biology research tool: evolvable open-source software. BMC Syst Biol. 2008, 2: 55-10.1186/1752-0509-2-55.PubMed CentralView ArticlePubMedGoogle Scholar
  22. Boele J, Olivier BG, Teusink B: FAME, the flux analysis and modeling environment. BMC Syst Biol. 2012, 6: 8-10.1186/1752-0509-6-8.PubMed CentralView ArticlePubMedGoogle Scholar
  23. Feng X, Xu Y, Chen Y, Tang YJ: MicrobesFlux: a web platform for drafting metabolic models from the KEGG database. BMC Syst Biol. 2012, 6: 94-10.1186/1752-0509-6-94.PubMed CentralView ArticlePubMedGoogle Scholar
  24. Kanehisa M, Goto S: KEGG: Kyoto encyclopedia of genes and genomes. MNucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.View ArticleGoogle Scholar
  25. Ebert BE, Lamprecht AL, Steffen B, Blank LM: Flux-P: automating metabolic flux analysis. Metabolites. 2012, 2: 872-890. 10.3390/metabo2040872.PubMed CentralView ArticlePubMedGoogle Scholar
  26. Jung TS, Yeo HC, Reddy SG, Cho WS, Lee DY: WEbcoli: an interactive and asynchronous web application for in silico design and analysis of genome-scale E.coli model. Bioinformatics. 2009, 25: 2850-2852. 10.1093/bioinformatics/btp496.View ArticlePubMedGoogle Scholar
  27. Ganter M, Bernard T, Moretti S, Stelling J, Pagni M: MetaNETX.org: a website and repository for accessing, analysing and manipulating metabolic networks. Bioinformatics. 2013, 29: 815-816. 10.1093/bioinformatics/btt036.PubMed CentralView ArticlePubMedGoogle Scholar
  28. Schellenberger J, Park JO, Conrad TM, Palsson BO: BiGG: a Biochemical Genetic and Genomic knowledgebase of large scale metabolic reconstructions. BMC Bioinformatics. 2010, 11: 213-10.1186/1471-2105-11-213.PubMed CentralView ArticlePubMedGoogle Scholar
  29. Kumar A, Suthers PF, Maranas CD: MetRxn: a knowledgebase of metabolites and reactions spanning metabolic models and databases. BMC Bioinformatics. 2012, 13: 6-10.1186/1471-2105-13-6.PubMed CentralView ArticlePubMedGoogle Scholar
  30. Giardine B, Riemer C, Hardison RC, Burhans R, Elnitski L, Shah P, Zhang Y, Blankenberg D, Albert I, Taylor J, Miller W, Kent WJ, Nekrutenko A: Galaxy: a platform for interactive large-scale genome analysis. Genome Res. 2005, 15: 1451-1455. 10.1101/gr.4086505.PubMed CentralView ArticlePubMedGoogle Scholar
  31. Thiele I, Vo TD, Price ND, Palsson BO: Expanded metabolic reconstruction of Helicobacter pylori (iIT341 GSM/GPR): an in silico genome-scale characterization of single- and double-deletion mutants. J Bacteriol. 2005, 187: 5818-5830. 10.1128/JB.187.16.5818-5830.2005.PubMed CentralView ArticlePubMedGoogle Scholar
  32. Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, Arkin AP, Bornstein BJ, Bray D, Cornish-Bowden A, Cuellar AA, Dronov S, Gilles ED, Ginkel M, Gor V, Goryanin II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, Juty NS, Kasberger JL, Kremling A, Kummer U, Le NN, Loew LM, Lucio D, Mendes P, Minch E, Mjolsness ED, et al: The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models. Bioinformatics. 2003, 19: 524-531. 10.1093/bioinformatics/btg015.View ArticlePubMedGoogle Scholar
  33. Reed JL, Vo TD, Schilling CH, Palsson BO: An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR). Genome Biol. 2003, 4: R54-10.1186/gb-2003-4-9-r54.PubMed CentralView ArticlePubMedGoogle Scholar
  34. Rocha I, Maia P, Evangelista P, Vilaça P, Soares S, Pinto JP, Nielsen J, Patil KR, Ferreira EC, Rocha M: OptFlux: an open-source software platform for in silico metabolic engineering. BMC Syst Biol. 2010, 4: 45-10.1186/1752-0509-4-45.PubMed CentralView ArticlePubMedGoogle Scholar
  35. Hasona A, Kim Y, Healy FG, Ingram LO, Shanmugam KT: Pyruvate formate and acetate kinase are essential for anaerobic growth of Escherichia coli on Xylose. J Bacteriol. 2004, 186: 7593-7600. 10.1128/JB.186.22.7593-7600.2004.PubMed CentralView ArticlePubMedGoogle Scholar
  36. Yeh I, Hanekamp T, Tsoka S, Karp PD, Altman RB: Computational analysis of Plasmodium falciparum metabolism: organizing genomic information to facilitate drug discovery. Genome Res. 2004, 14: 917-924. 10.1101/gr.2050304.PubMed CentralView ArticlePubMedGoogle Scholar
  37. Rahman SA, Schomburg D: Observing local and global properties of metabolic pathways: ‘load points’ and ‘choke points’ in the metabolic networks. Bioinformatics. 2006, 22: 1767-1774. 10.1093/bioinformatics/btl181.View ArticlePubMedGoogle Scholar
  38. Mahadevan R, Schilling CH: The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab Eng. 2003, 5: 264-276. 10.1016/j.ymben.2003.09.002.View ArticlePubMedGoogle Scholar

Copyright

© Narang et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Advertisement