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);...

Genome-scale metabolic modeling has become widespread for analyzing microbial metabolism. Extending this established paradigm to more complex microbial communities is emerging as a promising way to unravel the interactions and biochemical repertoire of these omnipresent systems. While several modeling techniques have been developed for microbial communities, little emphasis has been placed on the need to impose a time-averaged constant growth rate across all members for a community to ensure co-existence and stability. In the absence of this constraint, the faster growing organism will ultimately displace all other microbes in the community. This is particularly important for predicting steady-state microbiota composition as it imposes significant restrictions on the allowable community membership, composition and phenotypes. In this study, we introduce the SteadyCom optimization framework for predicting metabolic flux distributions consistent with the steady-state requirement. SteadyCom can be rapidly converged by iteratively solving linear programming (LP) problem and the number of iterations is independent of the number of organisms. A significant advantage of SteadyCom is compatibility with flux variability analysis. SteadyCom is first demonstrated for a community of four E. coli double auxotrophic mutants and is then applied to a gut microbiota model consisting of nine species, with representatives from the phyla Bacteroidetes, Firmicutes, Actinobacteria and Proteobacteria. In contrast to the direct use of FBA, SteadyCom is able to predict the change in species abundance in response to changes in diets with minimal additional imposed constraints on the model. By randomizing the uptake rates of microbes, an abundance profile with a good agreement to experimental gut microbiota is inferred. SteadyCom provides an important step towards the cross-cutting task of predicting the composition of a microbial community in a given environment.

Supplementary data are available at Bioinformatics online.

The farnesoid X receptor (FXR) plays an important role in mediating the dialog between the host and gut microbiota, particularly through modulation of enterohepatic circulation of bile acids. Mounting evidence suggests that genetic ablation of Fxr in the gut or gut-restricted chemical antagonism of the FXR promotes beneficial health effects, including the prevention of nonalcoholic fatty liver disease in rodent models. However, questions remain unanswered, including whether modulation of FXR activity plays a role in shaping the gut microbiota community structure and function and what metabolic pathways of the gut microbiota contribute in an FXR-dependent manner to the host phenotype. In this report, new insights are gained into the metabolic contribution of the gut microbiota to the metabolic phenotypes, including establishing a link between FXR antagonism, bacterial bile salt hydrolase activity, and fermentation. Multiple approaches, including unique mouse models as well as metabolomics and genome-scale metabolic models, were employed to confirm these results.

Introductory Paragraph Initial microbial colonization and later succession in the gut of human infants are linked to health and disease later in life. The timing of the appearance of the first gut microbiome, and the consequences for the early life metabolome, are just starting to be defined. Here we evaluated the gut microbiome, proteome, and metabolome in 88 African American newborns using fecal samples collected in the first few days of life. Gut bacteria became detectable using molecular methods by 16 hours after birth. Detailed analysis of the three most common species, Escherichia coli , Enterococcus faecalis , and Bacteroides vulgatus , did not suggest a genomic signature for neonatal gut colonization. The appearance of bacteria was associated with reduced abundance of approximately 50 human proteins, decreased levels of free amino acids, and an increase in products of bacterial fermentation, including acetate and succinate. Using flux balance modeling and in vitro experiments, we provide evidence that fermentation of amino acids provides a mechanism for the initial growth of Escherichia coli , the most common early colonizer, under anaerobic conditions. These results provide a deep characterization of the first microbes in the human gut and show how the biochemical environment is altered by their appearance.

No abstract

Rhodosporidium toruloides is a red, basidiomycetes yeast that can accumulate a large amount of lipids and produce carotenoids. To better assess this non-model yeast’s metabolic capabilities, we reconstructed a genome-scale model of R. toruloides IFO0880’s metabolic network (iRhto1108) accounting for 2204 reactions, 1985 metabolites and 1108 genes. In this work, we integrated and supplemented the current knowledge with in-house generated biomass composition and experimental measurements pertaining to the organism’s metabolic capabilities. Predictions of genotype-phenotype relations were improved through manual curation of gene-protein-reaction rules for 543 reactions leading to correct recapitulations of 84.5% of gene essentiality data (sensitivity of 94.3% and specificity of 53.8%). Organism-specific macromolecular composition and ATP maintenance requirements were experimentally measured for two separate growth conditions: (i) carbon and (ii) nitrogen limitations. Overall, iRhto1108 reproduced R. toruloides’s utilization capabilities for 18 alternate substrates, matched measured wild-type growth yield, and recapitulated the viability of 772 out of 819 deletion mutants. As a demonstration to the model’s fidelity in guiding engineering interventions, the OptForce procedure was applied on iRhto1108 for triacylglycerol overproduction. Suggested interventions recapitulated many of the previous successful implementations of genetic modifications and put forth a few new ones.

Motivation Microbial metabolic interactions impact ecosystems, human health and biotechnology profoundly. However, their determination remains elusive, invoking an urgent need for predictive models seamlessly integrating metabolism with evolutionary principles that shape community interactions. Results Inspired by the evolutionary game theory, we formulated a bi-level optimization framework termed NECom for which any feasible solutions are Nash equilibria of microbial community metabolic models with/without an outer-level (community) objective function. Distinct from discrete matrix games, NECom models the continuous interdependent strategy space of metabolic fluxes. We showed that NECom successfully predicted several classical games in the context of metabolic interactions that were falsely or incompletely predicted by existing methods, including prisoner’s dilemma, snowdrift and cooperation. The improved capability originates from the novel formulation to prevent ‘forced altruism’ hidden in previous static algorithms while allowing for sensing all potential metabolite exchanges to determine evolutionarily favorable interactions between members, a feature missing in dynamic methods. The results provided insights into why mutualism is favorable despite seemingly costly cross-feeding metabolites and demonstrated similarities and differences between games in the continuous metabolic flux space and matrix games. NECom was then applied to a reported algae-yeast co-culture system that shares typical cross-feeding features of lichen, a model system of mutualism. 488 growth conditions corresponding to 3,221 experimental data points were simulated. Without training any parameters using the data, NECom is more predictive of species’ growth rates given uptake rates compared with flux balance analysis with an overall 63.5% and 81.7% reduction in root-mean-square error for the two species. Availability Simulation code and data are available at https://github.com/Jingyi-Cai/NECom.git Supplementary information Supplementary data are available at Bioinformatics online.