Fast-SNP: a fast matrix pre-processing algorithm for efficient loopless flux optimization of metabolic models

Genome-scale metabolic modeling of the human milk oligosaccharide utilization by Bifidobacterium longum subsp. infantis

Bifidobacterium longum subsp. infantis is a representative and dominant species in the infant gut and is considered a beneficial microbe. This organism displays multiple adaptations to thrive in the infant gut, regarded as a model for human milk oligosaccharides (HMOs) utilization. These carbohydrates are abundant in breast milk and include different molecules based on lactose. They contain fucose, sialic acid, and N-acetylglucosamine. Bifidobacterium metabolism is complex, and a systems view of relevant metabolic pathways and exchange metabolites during HMO consumption is missing. To address this limitation, a refined genome-scale network reconstruction of this bacterium is presented using a previous reconstruction of B. infantis ATCC 15967 as a template. The latter was expanded based on an extensive revision of genome annotations, current literature, and transcriptomic data integration. The metabolic reconstruction (iLR578) accounted for 578 genes, 1,047 reactions, and 924 metabolites. Starting from this reconstruction, we built context-specific genome-scale metabolic models using RNA-seq data from cultures growing in lactose and three HMOs. The models revealed notable differences in HMO metabolism depending on the functional characteristics of the substrates. Particularly, fucosyl-lactose showed a divergent metabolism due to a fucose moiety. High yields of lactate and acetate were predicted under growth rate maximization in all conditions, whereas formate, ethanol, and 1,2-propanediol were substantially lower. Similar results were also obtained under near-optimal growth on each substrate when varying the empirically observed acetate-to-lactate production ratio. Model predictions displayed reasonable agreement between central carbon metabolism fluxes and expression data across all conditions. Flux coupling analysis revealed additional connections between succinate exchange and arginine and sulfate metabolism and a strong coupling between central carbon reactions and adenine metabolism. More importantly, specific networks of coupled reactions under each carbon source were derived and analyzed. Overall, the presented network reconstruction constitutes a valuable platform for probing the metabolism of this prominent infant gut bifidobacteria.IMPORTANCEThis work presents a detailed reconstruction of the metabolism of Bifidobacterium longum subsp. infantis, a prominent member of the infant gut microbiome, providing a systems view of its metabolism of human milk oligosaccharides.

LooplessFluxSampler: an efficient toolbox for sampling the loopless flux solution space of metabolic models

BACKGROUND

Uniform random sampling of mass-balanced flux solutions offers an unbiased appraisal of the capabilities of metabolic networks. Unfortunately, it is impossible to avoid thermodynamically infeasible loops in flux samples when using convex samplers on large metabolic models. Current strategies for randomly sampling the non-convex loopless flux space display limited efficiency and lack theoretical guarantees.

RESULTS

Here, we present LooplessFluxSampler, an efficient algorithm for exploring the loopless mass-balanced flux solution space of metabolic models, based on an Adaptive Directions Sampling on a Box (ADSB) algorithm. ADSB is rooted in the general Adaptive Direction Sampling (ADS) framework, specifically the Parallel ADS, for which theoretical convergence and irreducibility results are available for sampling from arbitrary distributions. By sampling directions that adapt to the target distribution, ADSB traverses more efficiently the sample space achieving faster mixing than other methods. Importantly, the presented algorithm is guaranteed to target the uniform distribution over convex regions, and it provably converges on the latter distribution over more general (non-convex) regions provided the sample can have full support.

CONCLUSIONS

LooplessFluxSampler enables scalable statistical inference of the loopless mass-balanced solution space of large metabolic models. Grounded in a theoretically sound framework, this toolbox provides not only efficient but also reliable results for exploring the properties of the almost surely non-convex loopless flux space. Finally, LooplessFluxSampler includes a Markov Chain diagnostics suite for assessing the quality of the final sample and the performance of the algorithm.

Local flux coordination and global gene expression regulation in metabolic modeling

Genome-scale metabolic networks (GSMs) are fundamental systems biology representations of a cell's entire set of stoichiometrically balanced reactions. However, such static GSMs do not incorporate the functional organization of metabolic genes and their dynamic regulation (e.g., operons and regulons). Specifically, there are numerous topologically coupled local reactions through which fluxes are coordinated; the global growth state often dynamically regulates many gene expression of metabolic reactions via global transcription factor regulators. Here, we develop a GSM reconstruction method, Decrem, by integrating locally coupled reactions and global transcriptional regulation of metabolism by cell state. Decrem produces predictions of flux and growth rates, which are highly correlated with those experimentally measured in both wild-type and mutants of three model microorganisms Escherichia coli, Saccharomyces cerevisiae, and Bacillus subtilis under various conditions. More importantly, Decrem can also explain the observed growth rates by capturing the experimentally measured flux changes between wild-types and mutants. Overall, by identifying and incorporating locally organized and regulated functional modules into GSMs, Decrem achieves accurate predictions of phenotypes and has broad applications in bioengineering, synthetic biology, and microbial pathology.

Metabolic Modeling of Wine Fermentation at Genome Scale

Wine fermentation is an ancient biotechnological process mediated by different microorganisms such as yeast and bacteria. Understanding of the metabolic and physiological phenomena taking place during this process can be now attained at a genome scale with the help of metabolic models. In this chapter, we present a detailed protocol for modeling wine fermentation using genome-scale metabolic models. In particular, we illustrate how metabolic fluxes can be computed, optimized and interpreted, for both yeast and bacteria under winemaking conditions. We also show how nutritional requirements can be determined and simulated using these models in relevant test cases. This chapter introduces fundamental concepts and practical steps for applying flux balance analysis in wine fermentation, and as such, it is intended for a broad microbiology audience as well as for practitioners in the metabolic modeling field.

ErrorTracer: an algorithm for identifying the origins of inconsistencies in genome-scale metabolic models

MOTIVATION

The number and complexity of genome-scale metabolic models is steadily increasing, empowered by automated model-generation algorithms. The quality control of the models, however, has always remained a significant challenge, the most fundamental being reactions incapable of carrying flux. Numerous automated gap-filling algorithms try to address this problem, but can rarely resolve all of a model's inconsistencies. The need for fast inconsistency checking algorithms has also been emphasized with the recent community push for automated model-validation before model publication. Previously, we wrote a graphical software to allow the modeller to solve the remaining errors manually. Nevertheless, model size and complexity remained a hindrance to efficiently tracking origins of inconsistency.

RESULTS

We developed the ErrorTracer algorithm in order to address the shortcomings of existing approaches: ErrorTracer searches for inconsistencies, classifies them and identifies their origins. The algorithm is ∼2 orders of magnitude faster than current community standard methods, using only seconds even for large-scale models. This allows for interactive exploration in direct combination with model visualization, markedly simplifying the whole error-identification and correction work flow.

AVAILABILITY AND IMPLEMENTATION

Windows and Linux executables and source code are available under the EPL 2.0 Licence at https://github.com/TheAngryFox/ModelExplorer and https://www.ntnu.edu/almaaslab/downloads.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Resources and tools for the high-throughput, multi-omic study of intestinal microbiota

The human gut microbiome impacts several aspects of human health and disease, including digestion, drug metabolism and the propensity to develop various inflammatory, autoimmune and metabolic diseases. Many of the molecular processes that play a role in the activity and dynamics of the microbiota go beyond species and genic composition and thus, their understanding requires advanced bioinformatics support. This article aims to provide an up-to-date view of the resources and software tools that are being developed and used in human gut microbiome research, in particular data integration and systems-level analysis efforts. These efforts demonstrate the power of standardized and reproducible computational workflows for integrating and analysing varied omics data and gaining deeper insights into microbe community structure and function as well as host-microbe interactions.

OptFill: A Tool for Infeasible Cycle-Free Gapfilling of Stoichiometric Metabolic Models

Stoichiometric metabolic modeling, particularly genome-scale models (GSMs), is now an indispensable tool for systems biology. The model reconstruction process typically involves collecting information from public databases; however, incomplete systems knowledge leaves gaps in any reconstruction. Current tools for addressing gaps use databases of biochemical functionalities to address gaps on a per-metabolite basis and can provide multiple solutions but cannot avoid thermodynamically infeasible cycles (TICs), invariably requiring lengthy manual curation. To address these limitations, this work introduces an optimization-based multi-step method named OptFill, which performs TIC-avoiding whole-model gapfilling. We applied OptFill to three fictional prokaryotic models of increasing sizes and to a published GSM of Escherichia coli, iJR904. This application resulted in holistic and infeasible cycle-free gapfilling solutions. In addition, OptFill can be adapted to automate inherent TICs identification in any GSM. Overall, OptFill can address critical issues in automated development of high-quality GSMs.

•

This work presents an alternative to state-of-the-art methods for gapfilling

•

Unlike current methods, this method is holistic and infeasible cycle free

•

This method is applied to three tests and one published model

•

This method might also be used to address infeasible cycling

Contextualized genome-scale model unveils high-order metabolic effects of the specific growth rate and oxygenation level in recombinant Pichia pastoris

Pichia pastoris is recognized as a biotechnological workhorse for recombinant protein expression. The metabolic performance of this microorganism depends on genetic makeup and culture conditions, amongst which the specific growth rate and oxygenation level are critical. Despite their importance, only their individual effects have been assessed so far, and thus their combined effects and metabolic consequences still remain to be elucidated. In this work, we present a comprehensive framework for revealing high-order (i.e., individual and combined) metabolic effects of the above parameters in glucose-limited continuous cultures of P. pastoris, using thaumatin production as a case study. Specifically, we employed a rational experimental design to calculate statistically significant metabolic effects from multiple chemostat data, which were later contextualized using a refined and highly predictive genome-scale metabolic model of this yeast under the simulated conditions. Our results revealed a negative effect of the oxygenation on the specific product formation rate (thaumatin), and a positive effect on the biomass yield. Notably, we identified a novel positive combined effect of both the specific growth rate and oxygenation level on the specific product formation rate. Finally, model predictions indicated an opposite relationship between the oxygenation level and the growth-associated maintenance energy (GAME) requirement, suggesting a linear GAME decrease of 0.56 mmol ATP/g_DCW per each 1% increase in oxygenation level, which translated into a 44% higher metabolic cost under low oxygenation compared to high oxygenation. Overall, this work provides a systematic framework for mapping high-order metabolic effects of different culture parameters on the performance of a microbial cell factory. Particularly in this case, it provided valuable insights about optimal operational conditions for protein production in P. pastoris.

In silico model-guided identification of transcriptional regulator targets for efficient strain design

BACKGROUND

Cellular metabolism is tightly regulated by hard-wired multiple layers of biological processes to achieve robust and homeostatic states given the limited resources. As a result, even the most intuitive enzyme-centric metabolic engineering endeavours through the up-/down-regulation of multiple genes in biochemical pathways often deliver insignificant improvements in the product yield. In this regard, targeted engineering of transcriptional regulators (TRs) that control several metabolic functions in modular patterns is an interesting strategy. However, only a handful of in silico model-added techniques are available for identifying the TR manipulation candidates, thus limiting its strain design application.

RESULTS

We developed hierarchical-Beneficial Regulatory Targeting (h-BeReTa) which employs a genome-scale metabolic model and transcriptional regulatory network (TRN) to identify the relevant TR targets suitable for strain improvement. We then applied this method to industrially relevant metabolites and cell factory hosts, Escherichia coli and Corynebacterium glutamicum. h-BeReTa suggested several promising TR targets, many of which have been validated through literature evidences. h-BeReTa considers the hierarchy of TRs in the TRN and also accounts for alternative metabolic pathways which may divert flux away from the product while identifying suitable metabolic fluxes, thereby performing superior in terms of global TR target identification.

CONCLUSIONS

In silico model-guided strain design framework, h-BeReTa, was presented for identifying transcriptional regulator targets. Its efficacy and applicability to microbial cell factories were successfully demonstrated via case studies involving two cell factory hosts, as such suggesting several intuitive targets for overproducing various value-added compounds.

From reconstruction to C4 metabolic engineering: A case study for overproduction of polyhydroxybutyrate in bioenergy grasses

The compartmentalization of C₄ plants increases photosynthetic efficiency, while constraining how material and energy must flow in leaf tissues. To capture this metabolic phenomenon, a generic plant metabolic reconstruction was replicated into four connected spatiotemporal compartments, namely bundle sheath (B) and mesophyll (M) across the day and night cycle. The C₄ leaf model was used to explore how amenable polyhydroxybutyrate (PHB) production is with these four compartments working cooperatively. A strategic pattern of metabolite conversion and exchange emerged from a systems-level network that has very few constraints imposed; mainly the sequential two-step carbon capture in mesophyll, then bundle sheath and photosynthesis during the day only. The building of starch reserves during the day and their mobilization during the night connects day and night metabolism. Flux simulations revealed that PHB production did not require rerouting of metabolic pathways beyond what is already utilised for growth. PHB yield was sensitive to photoassimilation capacity, availability of carbon reserves, ATP maintenance, relative photosynthetic activity of B and M, and type of metabolites exchanged in the plasmodesmata, but not sensitive towards compartmentalization. Hence, the compartmentalization issues currently encountered are likely to be kinetic or thermodynamic limitations rather than stoichiometric.