COnstraint-Based Reconstruction and Analysis (COBRA) provides a molecular mechanistic framework for integrative analysis of experimental data and quantitative prediction of physicochemically and biochemically feasible phenotypic states. The COBRA Toolbox is a comprehensive software suite of interoperable COBRA methods. It has found widespread applications in biology, biomedicine, and biotechnology because its functions can be flexibly combined to implement tailored COBRA protocols for any biochemical network. Version 3.0 includes new methods for quality controlled reconstruction, modelling, topological analysis, strain and experimental design, network visualisation as well as network integration of chemoinformatic, metabolomic, transcriptomic, proteomic, and thermochemical data. New multi-lingual code integration also enables an expansion in COBRA application scope via high-precision, high-performance, and nonlinear numerical optimisation solvers for multi-scale, multi-cellular and reaction kinetic modelling, respectively. This protocol can be adapted for the generation and analysis of a constraint-based model in a wide variety of molecular systems biology scenarios. This protocol is an update to the COBRA Toolbox 1.0 and 2.0. The COBRA Toolbox 3.0 provides an unparalleled depth of constraint-based reconstruction and analysis methods. ]); 61 | The MUST sets are the sets of reactions that must increase or decrease their flux in order to achieve the desired phenotype in the mutant strain. As shown in Figure 6, the first order MUST sets are MustU and MustL while second order MUST sets are denoted as MustUU, MustLL, and MustUL. After parameters and constraints are defined, the functions findMustL and findMustU are run to determine the mustU and mustL sets, respectively. Define an ID of the run with:Each time the MUST sets are determined, folders are generated to read inputs and store outputs, i.e., reports. These folders are located in the directory defined by the uniquely defined runID.62 | In order to find the first order MUST sets, constraints should be defined: >> constrOpt = struct('rxnList', {{'EX_gluc', 'R75', 'EX_suc'}}, 'values', [-100; 0; 155.5]); 63 | The first order MUST set MustL is determined by running: >> [mustLSet, pos_mustL] = findMustL(model, minFluxesW, maxFluxesW, ... 'constrOpt', constrOpt, 'runID', runID);If runID is set to 'TestoptForceL', a folder TestoptForceL is created, in which two additional folders InputsMustL and OutputsMustL are created. The InputsMustL folder contains all the inputs required to run the function findMustL, while the OutputsMustL folder contains the mustL set found and a report that summarises all the inputs and outputs. In order to maintain a chronological order of computational experiments, the report is timestamped.64 | Display the reactions that belong to the mustL set using: >> disp(mustLSet) 65 | The first order MUST set MustU is determined by running: >> [mustUSet, pos_mustU] = findMustU(model, minFluxesW, maxFluxesW, ... 'constrOpt', constrOpt, 'runID', runID);...

1 r e s o u r c eChanges in the composition of the human gut microbiota have been associated with the development of chronic diseases including type 2 diabetes, obesity, and colorectal cancer 1 . Gut bacterial functions, such as synthesis of amino acids and vitamins 2 , breakdown of indigestible plant polysaccharides 3 , and production of metabolites involved in energy metabolism 4 , have been linked to human health. The use of 'omics approaches to study human microbiome communities has led to the generation of enormous data sets whose interpretations require systems biology tools to shed light on the functional capacity of gut microbiomes and their interactions with the human host 5 .In order to infer the metabolic repertoire of a gut metagenome data set, researchers usually map sequenced genes or organisms onto metabolic networks derived from the Kyoto Encyclopedia of Genes and Genomes (KEGG) 6 , and functional annotations from KEGG ontologies 7 . However, this approach cannot identify the contribution of each bacterial species to the metabolic repertoire of the whole gut microbiome, nor can it infer the effects of different gut microbial communities on host metabolism.A technique that can bridge this gap is constraint-based reconstruction and analysis (COBRA) 8 using genome-scale metabolic reconstructions (GENREs) of individual human gut microbes. GENREs are assembled using the genome sequence and experimental information 9 . These reconstructions form the basis for the development of condition-specific metabolic models whose functions are simulated and validated by comparison with experimental results. The models can be used to investigate genotype-phenotype relationships 10 , microbe-microbe interactions 11 , and host-microbe interactions 11 . Numerous tools can be used to automatically generate draft GENREs but such models contain errors 12 and are incomplete.Manual curation of draft reconstructions is time consuming because it involves an extensive literature review and experimental validation of metabolic functions 9 .To provide an extensive resource of GENREs for human gut microbes, we developed a comparative metabolic reconstruction method that enables any refinement to one metabolic reconstruction to be propagated to others. This accelerates reconstruction and improves model quality. We generated AGORA, which includes 773 gut microbes, comprising 205 genera and 605 species. All reconstructions were based on literature-derived experimental data and comparative genomics. The metabolic predictions of two AGORA reconstructions and their derived metabolic models were validated against experimental data. RESULTS Metabolic reconstruction pipelineWe devised a comparative metabolic reconstruction method (Fig. 1a,c), which is analogous to the comparative microbial genome annotation approach 13 that has enabled accelerated annotation by propagation of refinements to one genome to others. First, we downloaded draft GENREs using Model SEED 14 and KBase (US Department of Energy Systems Biology Knowledgebase, http:/...

Genome-scale network reconstructions have helped uncover the molecular basis of metabolism. Here we present Recon3D, a computational resource that includes three-dimensional (3D) metabolite and protein structure data and enables integrated analyses of metabolic functions in humans. We use Recon3D to functionally characterize mutations associated with disease, and identify metabolic response signatures that are caused by exposure to certain drugs. Recon3D represents the most comprehensive human metabolic network model to date, accounting for 3,288 open reading frames (representing 17% of functionally annotated human genes), 13,543 metabolic reactions involving 4,140 unique metabolites, and 12,890 protein structures. These data provide a unique resource for investigating molecular mechanisms of human metabolism. Recon3D is available at http://vmh.life.

A multitude of factors contribute to complex diseases and can be measured with ‘omics’ methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the Virtual Metabolic Human (VMH, www.vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources ‘Human metabolism’, ‘Gut microbiome’, ‘Disease’, ‘Nutrition’, and ‘ReconMaps’. The VMH captures 5180 unique metabolites, 17 730 unique reactions, 3695 human genes, 255 Mendelian diseases, 818 microbes, 632 685 microbial genes and 8790 food items. The VMH’s unique features are (i) the hosting of the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; (ii) seven human metabolic maps for data visualization; (iii) a nutrition designer; (iv) a user-friendly webpage and application-programming interface to access its content; (v) user feedback option for community engagement and (vi) the connection of its entities to 57 other web resources. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community.

A multitude of factors contribute to complex diseases and can be measured with "omics" methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the Virtual Metabolic Human (VMH, http://vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources "Human metabolism", "Gut microbiome", "Disease", "Nutrition", and "ReconMaps". The VMH captures 5,180 unique metabolites, 17,730 unique reactions, 3,288 human genes, 255 Mendelian diseases, 818 microbes, 632,685 microbial genes, and 8,790 food items. The VMH's unique features are i) the hosting the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; ii) seven human metabolic maps for data visualization; iii) a nutrition designer; iv) a user-friendly webpage and application-programming interface to access its content; and v) user feedback option for community engagement. We demonstrate with four examples the VMH's utility. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community.

No abstract

Constraint-based analysis of genomescale models (GEMs) arose shortly after the first genome sequences became available. As numerous reviews of the field show, this approach and methodology has proven to be successful in studying a wide range of biological phenomena (McCloskey et al, 2013;Bordbar et al, 2014). However, efforts to expand the user base are impeded by hurdles in correctly formulating these problems to obtain numerical solutions. In particular, in a study entitled "An exact arithmetic toolbox for a consistent and reproducible structural analysis of metabolic network models" (Chindelevitch et al, 2014), the authors apply an exact solver to 88 genome-scale constraint-based models of metabolism. The authors claim that COBRA calculations (Orth et al, 2010) are inconsistent with their results and that many published and actively

MotivationA genome-scale reconstruction of human metabolism, Recon 2, is available but no interface exists to interactively visualize its content integrated with omics data and simulation results.ResultsWe manually drew a comprehensive map, ReconMap 2.0, that is consistent with the content of Recon 2. We present it within a web interface that allows content query, visualization of custom datasets and submission of feedback to manual curators.Availability and ImplementationReconMap can be accessed via http://vmh.uni.lu, with network export in a Systems Biology Graphical Notation compliant format released under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License. A Constraint-Based Reconstruction and Analysis (COBRA) Toolbox extension to interact with ReconMap is available via https://github.com/opencobra/cobratoolbox.