Comparison of metabolite production capability indices generated by network analysis methods

A framework of constraint-based reconstruction and analysis (COBRA) is used for modeling large-scale metabolic networks. In COBRA, extreme pathway and optimization analyses are commonly used to study the properties of networks. While the results of both methods are completely consistent, extreme pathway analysis is considered to be better because of its wider representational ability. In this study, we assessed these two methods by computational knockout experiments. We examined a simple pathway model and found that the extreme pathway method led to misguided conclusions in specific cases, while optimization analysis calculated the correct knockout effects. We also investigated the Escherichia coli metabolic pathway model, and found that these methods result in inconsistent interpretations of the network properties. IN CONCLUSION: it has been claimed that these two methods result in the same producible metabolites, but we found a difference in individual results for a biological pathway. Our results could provide helpful guidance for when to use the methods, particularly extreme pathway analysis.     

Paint4Net: COBRA Toolbox extension for visualization of stoichiometric models of metabolism

A visual analysis of reconstructions and large stoichiometric models with elastic change of the visualization scope and representation methods becomes increasingly important due to the rapidly growing size and number of available reconstructions. The Paint4Net is a novel COBRA Toolbox extension for automatic generation of a hypergraph layout of defined scope with the steady state rates of reaction fluxes of stoichiometric models. Directionalities and fluxes of reactions are constantly represented in the visualization while detailed information about reaction (ID, name and synonyms, and formula) and metabolite (ID, name and synonyms, and charged formula) appears placing the cursor on the item of interest. Additionally Paint4Net functionality can be used to: (1) get lists of involved metabolites and dead end metabolites of the visualized part of the network, (2) exclude (filter) particular metabolites from representation, (3) find isolated parts of a network and (4) find running cycles when all the substrates are cut down. Layout pictures can be saved in various formats and easily distributed. The Paint4Net is open source software under the GPL v3 license. Relevant documentation and sample data is available at http://www.biosystems.lv/paint4net. The Paint4Net works on MATLAB starting from version of 2009.     

基于约束建模法的结核菌H37Rv代谢网络分析

基于结核分枝杆菌国际标准强毒株H37Rv菌株的基因组尺度代谢网络模型iNJ661进行分析,以寻找代谢网络中培养基的关键成分和必要基因.该研究在Matlab平台上利用COBRA工具箱,采用基于约束的建模方法进行动态生长模拟、解空间抽样在酶活性水平上的具体化和基因删除模拟实验.结果发现培养基成分中铵盐、三价铁盐、磷酸盐、硫酸盐、甘油等可影响H37Rv的生长;培养基中去除磷酸盐后十种酶均在不同程度上受到抑制,其中丙糖磷酸异构酶、3-磷酸甘油醛脱氢酶、磷酸甘油酸变位酶、烯醇酶受限明显.通过基因删除得出188个必要基因以及非必要基因中的16个致死基因对.基于约束建模分析可初步了解结核杆菌H37Rv菌株代谢网络的性质,可为后续相关研究提供参考和借鉴.     

Elimination of thermodynamically infeasible loops in steady-state metabolic models

The constraint-based reconstruction and analysis (COBRA) framework has been widely used to study steady-state flux solutions in genome-scale metabolic networks. One shortcoming of current COBRA methods is the possible violation of the loop law in the computed steady-state flux solutions. The loop law is analogous to Kirchhoff's second law for electric circuits, and states that at steady state there can be no net flux around a closed network cycle. Although the consequences of the loop law have been known for years, it has been computationally difficult to work with. Therefore, the resulting loop-law constraints have been overlooked. Here, we present a general mixed integer programming approach called loopless COBRA (ll-COBRA), which can be used to eliminate all steady-state flux solutions that are incompatible with the loop law. We apply this approach to improve flux predictions on three common COBRA methods: flux balance analysis, flux variability analysis, and Monte Carlo sampling of the flux space. Moreover, we demonstrate that the imposition of loop-law constraints with ll-COBRA improves the consistency of simulation results with experimental data. This method provides an additional constraint for many COBRA methods, enabling the acquisition of more realistic simulation results.     

Systems Biology Toolbox for MATLAB: a computational platform for research in systems biology

We present a Systems Biology Toolbox for the widely used general purpose mathematical software MATLAB. The toolbox offers systems biologists an open and extensible environment, in which to explore ideas, prototype and share new algorithms, and build applications for the analysis and simulation of biological and biochemical systems. Additionally it is well suited for educational purposes. The toolbox supports the Systems Biology Markup Language (SBML) by providing an interface for import and export of SBML models. In this way the toolbox connects nicely to other SBML-enabled modelling packages. Models are represented in an internal model format and can be described either by entering ordinary differential equations or, more intuitively, by entering biochemical reaction equations. The toolbox contains a large number of analysis methods, such as deterministic and stochastic simulation, parameter estimation, network identification, parameter sensitivity analysis and bifurcation analysis.     

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.     

BiGG: a Biochemical Genetic and Genomic knowledgebase of large scale metabolic reconstructions

Background  

Genome-scale metabolic reconstructions under the Constraint Based Reconstruction and Analysis (COBRA) framework are valuable tools for analyzing the metabolic capabilities of organisms and interpreting experimental data. As the number of such reconstructions and analysis methods increases, there is a greater need for data uniformity and ease of distribution and use.     

An application programming interface for CellNetAnalyzer

CellNetAnalyzer (CNA) is a MATLAB toolbox providing computational methods for studying structure and function of metabolic and cellular signaling networks. In order to allow non-experts to use these methods easily, CNA provides GUI-based interactive network maps as a means of parameter input and result visualization. However, with the availability of high-throughput data, there is a need to make CNA's functionality also accessible in batch mode for automatic data processing. Furthermore, as some algorithms of CNA are of general relevance for network analysis it would be desirable if they could be called as sub-routines by other applications. For this purpose, we developed an API (application programming interface) for CNA allowing users (i) to access the content of network models in CNA, (ii) to use CNA's network analysis capabilities independent of the GUI, and (iii) to interact with the GUI to facilitate the development of graphical plugins.Here we describe the organization of network projects in CNA and the application of the new API functions to these projects. This includes the creation of network projects from scratch, loading and saving of projects and scenarios, and the application of the actual analysis methods. Furthermore, API functions for the import/export of metabolic models in SBML format and for accessing the GUI are described. Lastly, two example applications demonstrate the use and versatile applicability of CNA's API. CNA is freely available for academic use and can be downloaded from http://www.mpi-magdeburg.mpg.de/projects/cna/cna.html.     

MetGEMs Toolbox: Metagenome-scale models as integrative toolbox for uncovering metabolic functions and routes of human gut microbiome

Investigating metabolic functional capability of a human gut microbiome enables the quantification of microbiome changes, which can cause a phenotypic change of host physiology and disease. One possible way to estimate the functional capability of a microbial community is through inferring metagenomic content from 16S rRNA gene sequences. Genome-scale models (GEMs) can be used as scaffold for functional estimation analysis at a systematic level, however up to date, there is no integrative toolbox based on GEMs for uncovering metabolic functions. Here, we developed the MetGEMs (metagenome-scale models) toolbox, an open-source application for inferring metabolic functions from 16S rRNA gene sequences to facilitate the study of the human gut microbiome by the wider scientific community. The developed toolbox was validated using shotgun metagenomic data and shown to be superior in predicting functional composition in human clinical samples compared to existing state-of-the-art tools. Therefore, the MetGEMs toolbox was subsequently applied for annotating putative enzyme functions and metabolic routes related in human disease using atopic dermatitis as a case study.     

MapMaker and PathTracer for tracking carbon in genome‐scale metabolic models

Constraint‐based reconstruction and analysis (COBRA) modeling results can be difficult to interpret given the large numbers of reactions in genome‐scale models. While paths in metabolic networks can be found, existing methods are not easily combined with constraint‐based approaches. To address this limitation, two tools (MapMaker and PathTracer) were developed to find paths (including cycles) between metabolites, where each step transfers carbon from reactant to product. MapMaker predicts carbon transfer maps (CTMs) between metabolites using only information on molecular formulae and reaction stoichiometry, effectively determining which reactants and products share carbon atoms. MapMaker correctly assigned CTMs for over 97% of the 2,251 reactions in an Escherichia coli metabolic model (iJO1366). Using CTMs as inputs, PathTracer finds paths between two metabolites. PathTracer was applied to iJO1366 to investigate the importance of using CTMs and COBRA constraints when enumerating paths, to find active and high flux paths in flux balance analysis (FBA) solutions, to identify paths for putrescine utilization, and to elucidate a potential CO₂ fixation pathway in E. coli. These results illustrate how MapMaker and PathTracer can be used in combination with constraint‐based models to identify feasible, active, and high flux paths between metabolites.     

Quantitative analysis of Escherichia coli metabolic phenotypes within the context of phenotypic phase planes

In silico models of Escherichia coli metabolism have been developed to predict metabolic behavior and propose experimentally testable hypotheses. However, a thorough assessment of the metabolic phenotype requires well-designed experimentation and reproducible experimental techniques. A method for the quantitative analysis of E. coli metabolism in vivo within the framework of in silico phenotypic phase plane analysis is presented. Using this approach, we have quantitatively studied E. coli metabolism in various environmental conditions and nutritional media. Our experimental methodology, in combination with steady-state metabolic models, can be used to study biological properties and evaluate the metabolic capabilities of microbes.     

AutoClickChem: click chemistry in silico

Academic researchers and many in industry often lack the financial resources available to scientists working in "big pharma." High costs include those associated with high-throughput screening and chemical synthesis. In order to address these challenges, many researchers have in part turned to alternate methodologies. Virtual screening, for example, often substitutes for high-throughput screening, and click chemistry ensures that chemical synthesis is fast, cheap, and comparatively easy. Though both in silico screening and click chemistry seek to make drug discovery more feasible, it is not yet routine to couple these two methodologies. We here present a novel computer algorithm, called AutoClickChem, capable of performing many click-chemistry reactions in silico. AutoClickChem can be used to produce large combinatorial libraries of compound models for use in virtual screens. As the compounds of these libraries are constructed according to the reactions of click chemistry, they can be easily synthesized for subsequent testing in biochemical assays. Additionally, in silico modeling of click-chemistry products may prove useful in rational drug design and drug optimization. AutoClickChem is based on the pymolecule toolbox, a framework that may facilitate the development of future python-based programs that require the manipulation of molecular models. Both the pymolecule toolbox and AutoClickChem are released under the GNU General Public License version 3 and are available for download from http://autoclickchem.ucsd.edu.     

Metabolic modeling of microbial strains in silico

The large volume of genome-scale data that is being produced and made available in databases on the World Wide Web is demanding the development of integrated mathematical models of cellular processes. The analysis of reconstructed metabolic networks as systems leads to the development of an in silico or computer representation of collections of cellular metabolic constituents, their interactions and their integrated function as a whole. The use of quantitative analysis methods to generate testable hypotheses and drive experimentation at a whole-genome level signals the advent of a systemic modeling approach to cellular and molecular biology.     

Genome-scale models of microbial cells: evaluating the consequences of constraints

Microbial cells operate under governing constraints that limit their range of possible functions. With the availability of annotated genome sequences, it has become possible to reconstruct genome-scale biochemical reaction networks for microorganisms. The imposition of governing constraints on a reconstructed biochemical network leads to the definition of achievable cellular functions. In recent years, a substantial and growing toolbox of computational analysis methods has been developed to study the characteristics and capabilities of microorganisms using a constraint-based reconstruction and analysis (COBRA) approach. This approach provides a biochemically and genetically consistent framework for the generation of hypotheses and the testing of functions of microbial cells.     

BicAT: a biclustering analysis toolbox

SUMMARY: Besides classical clustering methods such as hierarchical clustering, in recent years biclustering has become a popular approach to analyze biological data sets, e.g. gene expression data. The Biclustering Analysis Toolbox (BicAT) is a software platform for clustering-based data analysis that integrates various biclustering and clustering techniques in terms of a common graphical user interface. Furthermore, BicAT provides different facilities for data preparation, inspection and postprocessing such as discretization, filtering of biclusters according to specific criteria or gene pair analysis for constructing gene interconnection graphs. The possibility to use different biclustering algorithms inside a single graphical tool allows the user to compare clustering results and choose the algorithm that best fits a specific biological scenario. The toolbox is described in the context of gene expression analysis, but is also applicable to other types of data, e.g. data from proteomics or synthetic lethal experiments. AVAILABILITY: The BicAT toolbox is freely available at http://www.tik.ee.ethz.ch/sop/bicat and runs on all operating systems. The Java source code of the program and a developer's guide is provided on the website as well. Therefore, users may modify the program and add further algorithms or extensions.     

Metabolic network modeling and simulation for drug targeting and discovery

Systems biology has greatly contributed toward the analysis and understanding of biological systems under various genotypic and environmental conditions on a much larger scale than ever before. One of the applications of systems biology can be seen in unraveling and understanding complicated human diseases where the primary causes for a disease are often not clear. The in silico genome-scale metabolic network models can be employed for the analysis of diseases and for the discovery of novel drug targets suitable for treating the disease. Also, new antimicrobial targets can be discovered by analyzing, at the systems level, the genome-scale metabolic network of pathogenic microorganisms. Such applications are possible as these genome-scale metabolic network models contain extensive stoichiometric relationships among the metabolites constituting the organism's metabolism and information on the associated biophysical constraints. In this review, we highlight applications of genome-scale metabolic network modeling and simulations in predicting drug targets and designing potential strategies in combating pathogenic infection. Also, the use of metabolic network models in the systematic analysis of several human diseases is examined. Other computational and experimental approaches are discussed to complement the use of metabolic network models in the analysis of biological systems and to facilitate the drug discovery pipeline.     

Investigation of metabolic objectives in cultured hepatocytes

Using optimization based methods to predict fluxes in metabolic flux balance models has been a successful approach for some microorganisms, enabling construction of in silico models and even inference of some regulatory motifs. However, this success has not been translated to mammalian cells. The lack of knowledge about metabolic objectives in mammalian cells is a major obstacle that prevents utilization of various metabolic engineering tools and methods for tissue engineering and biomedical purposes. In this work, we investigate and identify possible metabolic objectives for hepatocytes cultured in vitro. To achieve this goal, we present a special data-mining procedure for identifying metabolic objective functions in mammalian cells. This multi-level optimization based algorithm enables identifying the major fluxes in the metabolic objective from MFA data in the absence of information about critical active constraints of the system. Further, once the objective is determined, active flux constraints can also be identified and analyzed. This information can be potentially used in a predictive manner to improve cell culture results or clinical metabolic outcomes. As a result of the application of this method, it was found that in vitro cultured hepatocytes maximize oxygen uptake, coupling of urea and TCA cycles, and synthesis of serine and urea. Selection of these fluxes as the metabolic objective enables accurate prediction of the flux distribution in the system given a limited amount of flux data; thus presenting a workable in silico model for cultured hepatocytes. It is observed that an overall homeostasis picture is also emergent in the findings.