The effects of alternate optimal solutions in constraint-based genome-scale metabolic models

A diel multi-tissue genome-scale metabolic model of Vitis vinifera

Vitis vinifera, also known as grapevine, is widely cultivated and commercialized, particularly to produce wine. As wine quality is directly linked to fruit quality, studying grapevine metabolism is important to understand the processes underlying grape composition. Genome-scale metabolic models (GSMMs) have been used for the study of plant metabolism and advances have been made, allowing the integration of omics datasets with GSMMs. On the other hand, Machine learning (ML) has been used to analyze and integrate omics data, and while the combination of ML with GSMMs has shown promising results, it is still scarcely used to study plants. Here, the first GSSM of V. vinifera was reconstructed and validated, comprising 7199 genes, 5399 reactions, and 5141 metabolites across 8 compartments. Tissue-specific models for the stem, leaf, and berry of the Cabernet Sauvignon cultivar were generated from the original model, through the integration of RNA-Seq data. These models have been merged into diel multi-tissue models to study the interactions between tissues at light and dark phases. The potential of combining ML with GSMMs was explored by using ML to analyze the fluxomics data generated by green and mature grape GSMMs and provide insights regarding the metabolism of grapes at different developmental stages. Therefore, the models developed in this work are useful tools to explore different aspects of grapevine metabolism and understand the factors influencing grape quality.

The marine environmental microbiome mediates physiological outcomes in host nematodes

BACKGROUND

Nematodes are the most abundant metazoans in marine sediments, many of which are bacterivores; however, how habitat bacteria affect physiological outcomes in marine nematodes remains largely unknown.  RESULTS: Here, we used a Litoditis marina inbred line to assess how native bacteria modulate host nematode physiology. We characterized seasonal dynamic bacterial compositions in L. marina habitats and examined the impacts of 448 habitat bacteria isolates on L. marina development, then focused on HQbiome with 73 native bacteria, of which we generated 72 whole genomes sequences. Unexpectedly, we found that the effects of marine native bacteria on the development of L. marina and its terrestrial relative Caenorhabditis elegans were significantly positively correlated. Next, we reconstructed bacterial metabolic networks and identified several bacterial metabolic pathways positively correlated with L. marina development (e.g., ubiquinol and heme b biosynthesis), while pyridoxal 5'-phosphate biosynthesis pathway was negatively associated. Through single metabolite supplementation, we verified CoQ10, heme b, acetyl-CoA, and acetaldehyde promoted L. marina development, while vitamin B6 attenuated growth. Notably, we found that only four development correlated metabolic pathways were shared between L. marina and C. elegans. Furthermore, we identified two bacterial metabolic pathways correlated with L. marina lifespan, while a distinct one in C. elegans. Strikingly, we found that glycerol supplementation significantly extended L. marina but not C. elegans longevity. Moreover, we comparatively demonstrated the distinct gut microbiota characteristics and their effects on L. marina and C. elegans physiology.

CONCLUSIONS

Given that both bacteria and marine nematodes are dominant taxa in sedimentary ecosystems, the resource presented here will provide novel insights to identify mechanisms underpinning how habitat bacteria affect nematode biology in a more natural context. Our integrative approach will provide a microbe-nematodes framework for microbiome mediated effects on host animal fitness.

Metabolic modelling identifies mitochondrial Pi uptake and pyruvate efflux as key aspects of daytime metabolism and proton homeostasis in crassulacean acid metabolism leaves

Crassulacean acid metabolism (CAM) leaves are characterized by nocturnal acidification and diurnal deacidification processes related with the timed actions of phosphoenolpyruvate carboxylase and Rubisco, respectively. How CAM leaves manage cytosolic proton homeostasis, particularly when facing massive diurnal proton effluxes from the vacuole, remains unclear. A 12-phase flux balance analysis (FBA) model was constructed for a mature malic enzyme-type CAM mesophyll cell in order to predict diel kinetics of intracellular proton fluxes. The charge- and proton-balanced FBA model identified the mitochondrial phosphate carrier (PiC, Pi/H+ symport), which provides Pi to the matrix to sustain ATP biosynthesis, as a major consumer of cytosolic protons during daytime (> 50%). The delivery of Pi to the mitochondrion, co-transported with protons, is required for oxidative phosphorylation and allows sufficient ATP to be synthesized to meet the high energy demand during CAM Phase III. Additionally, the model predicts that mitochondrial pyruvate originating from decarboxylation of malate is exclusively exported to the cytosol, probably via a pyruvate channel mechanism, to fuel gluconeogenesis. In this biochemical cycle, glyceraldehyde 3-phosphate dehydrogenase (GAPDH) acts as another important cytosolic proton consumer. Overall, our findings emphasize the importance of mitochondria in CAM and uncover a hitherto unappreciated role in metabolic proton homeostasis.

Insights into yeast response to chemotherapeutic agent through time series genome-scale metabolic models

Organism-specific genome-scale metabolic models (GSMMs) can unveil molecular mechanisms within cells and are commonly used in diverse applications, from synthetic biology, biotechnology, and systems biology to metabolic engineering. There are limited studies incorporating time-series transcriptomics in GSMM simulations. Yeast is an easy-to-manipulate model organism for tumor research. Here, a novel approach (TS-GSMM) was proposed to integrate time-series transcriptomics with GSMMs to narrow down the feasible solution space of all possible flux distributions and attain time-series flux samples. The flux samples were clustered using machine learning techniques, and the clusters' functional analysis was performed using reaction set enrichment analysis. A time series transcriptomics response of Yeast cells to a chemotherapeutic reagent-doxorubicin-was mapped onto a Yeast GSMM. Eleven flux clusters were obtained with our approach, and pathway dynamics were displayed. Induction of fluxes related to bicarbonate formation and transport, ergosterol and spermidine transport, and ATP production were captured. Integrating time-series transcriptomics data with GSMMs is a promising approach to reveal pathway dynamics without any kinetic modeling and detects pathways that cannot be identified through transcriptomics-only analysis. The codes are available at https://github.com/karabekmez/TS-GSMM.

Understanding flux switching in metabolic networks through an analysis of synthetic lethals

Biological systems are robust and redundant. The redundancy can manifest as alternative metabolic pathways. Synthetic double lethals are pairs of reactions that, when deleted simultaneously, abrogate cell growth. However, removing one reaction allows the rerouting of metabolites through alternative pathways. Little is known about these hidden linkages between pathways. Understanding them in the context of pathogens is useful for therapeutic innovations. We propose a constraint-based optimisation approach to identify inter-dependencies between metabolic pathways. It minimises rerouting between two reaction deletions, corresponding to a synthetic lethal pair, and outputs the set of reactions vital for metabolic rewiring, known as the synthetic lethal cluster. We depict the results for different pathogens and show that the reactions span across metabolic modules, illustrating the complexity of metabolism. Finally, we demonstrate how the two classes of synthetic lethals play a role in metabolic networks and influence the different properties of a synthetic lethal cluster.

Refining microbial community metabolic models derived from metagenomics using reference-based taxonomic profiling

Characterization of microbial community metabolic output is crucial to understanding their functions. Construction of genome-scale metabolic models from metagenome-assembled genomes (MAG) has enabled prediction of metabolite production by microbial communities, yet little is known about their accuracy. Here, we examined the performance of two approaches for metabolite prediction from metagenomes, one that is MAG-guided and another that is taxonomic reference-guided. We applied both on shotgun metagenomics data from human and environmental samples, and validated findings in the human samples using untargeted metabolomics. We found that in human samples, where taxonomic profiling is optimized and reference genomes are readily available, when number of input taxa was normalized, the reference-guided approach predicted more metabolites than the MAG-guided approach. The two approaches showed significant overlap but each identified metabolites not predicted in the other. Pathway enrichment analyses identified significant differences in inferences derived from data based on the approach, highlighting the need for caution in interpretation. In environmental samples, when the number of input taxa was normalized, the reference-guided approach predicted more metabolites than the MAG-guided approach for total metabolites in both sample types and non-redundant metabolites in seawater samples. Nonetheless, as was observed for the human samples, the approaches overlapped substantially but also predicted metabolites not observed in the other. Our findings report on utility of a complementary input to genome-scale metabolic model construction that is less computationally intensive forgoing MAG assembly and refinement, and that can be applied on shallow shotgun sequencing where MAGs cannot be generated.IMPORTANCELittle is known about the accuracy of genome-scale metabolic models (GEMs) of microbial communities despite their influence on inferring community metabolic outputs and culture conditions. The performance of GEMs for metabolite prediction from metagenomes was assessed by applying two approaches on shotgun metagenomics data from human and environmental samples, and validating findings in the human samples using untargeted metabolomics. The performance of the approach was found to be dependent on sample type, but collectively, the reference-guided approach predicted more metabolites than the MAG-guided approach. Despite the differences, the predictions from the approaches overlapped substantially but each identified metabolites not predicted in the other. We found significant differences in biological inferences based on the approach, with some examples of uniquely enriched pathways in one group being invalidated when using the alternative approach, highlighting the need for caution in interpretation of GEMs.

Comparing genome-scale metabolic models of the non-resistant Enterococcus faecalis ATCC 19433 and the multi-resistant Enterococcus faecalis V583

Enterococcus faecalis is a versatile lactic acid bacterium with a large variety of implications for humans. While some strains of this species are pathobionts being resistant against most of the common antibiotics, other strains are regarded as biological protectants or even probiotics. Accordingly, E. faecalis strains largely differ in the size and content of their accessory genome. In this study, we describe the genome-scale metabolic network reconstruction of E. faecalis ATCC 19433, a non-resistant human-associated strain. A comparison of the genome-scale metabolic model (GSM) of E. faecalis ATCC 19433 with a previously published GSM of the multi-resistant pathobiontic E. faecalis V583 reveals high similarities in the central metabolic abilities of these two human associated strains. This is reflected, e.g., in the identical amino acid auxotrophies. The ATCC 19433 strain, however, has a 14.1 % smaller genome than V583 and lacks the multiple antibiotic resistance genes and genes involved in capsule formation. Based on the measured metabolic fluxes at different growth rates, the energy demand at zero growth was calculated to be about 40 % lower for the ATCC 19433 strain compared to V583. Furthermore, the ATCC 19433 strain seems less prone to the depletion of amino acids utilizable for energy metabolism. This might hint at a lower overall energy demand of the ATCC 19433 strain as compared to V583.

Deciphering Antibiotic-Targeted Metabolic Pathways in Acinetobacter baumannii: Insights from Transcriptomics and Genome-Scale Metabolic Modeling

In the ongoing battle against antibiotic-resistant infections, Acinetobacter baumannii has emerged as a critical pathogen in healthcare settings. To understand its response to antibiotic-induced stress, we integrated transcriptomic data from various antibiotics (amikacin sulfate, ciprofloxacin, polymyxin-B, and meropenem) with metabolic modeling techniques. Key metabolic pathways, including arginine and proline metabolism, glycine-serine and threonine metabolism, glyoxylate and dicarboxylate metabolism, and propanoate metabolism, were significantly impacted by all four antibiotics across multiple strains. Specifically, biotin metabolism was consistently down-regulated under polymyxin-B treatment, while fatty acid metabolism was perturbed under amikacin sulfate. Ciprofloxacin induced up-regulation in glycerophospholipid metabolism. Validation with an independent dataset focusing on colistin treatment confirmed alterations in fatty acid degradation, elongation, and arginine metabolism. By harmonizing genetic data with metabolic modeling and a metabolite-centric approach, our findings offer insights into the intricate adaptations of A. baumannii under antibiotic pressure, suggesting more effective strategies to combat antibiotic-resistant infections.

Measuring and understanding information storage and transfer in a simulated human gut microbiome

Considering biological systems as information processing entities and analyzing their organizational structure via information-theoretic measures has become an established approach in life sciences. We transfer this framework to a field of broad general interest, the human gut microbiome. We use BacArena, a software combining agent-based modelling and flux-balance analysis, to simulate a simplified human intestinal microbiome (SIHUMI). In a first step, we derive information theoretic measures from the simulated abundance data, and, in a second step, relate them to the metabolic processes underlying the abundance data. Our study provides further evidence on the role of active information storage as an indicator of unexpected structural change in the observed system. Besides, we show that information transfer reflects coherent behavior in the microbial community, both as a reaction to environmental changes and as a result of direct effective interaction. In this sense, purely abundance-based information theoretic measures can provide meaningful insight on metabolic interactions within bacterial communities. Furthermore, we shed light on the important however little noticed technical aspect of distinguishing immediate and delayed effects in the interpretation of local information theoretical measures.

Multi-omic characterization of antibody-producing CHO cell lines elucidates metabolic reprogramming and nutrient uptake bottlenecks

Characterizing the phenotypic diversity and metabolic capabilities of industrially relevant manufacturing cell lines is critical to bioprocess optimization and cell line development. Metabolic capabilities of production hosts limit nutrient and resource channeling into desired cellular processes and can have a profound impact on productivity. These limitations cannot be directly inferred from measured data such as spent media concentrations or transcriptomics. Here, we present an integrated multi-omic analysis pipeline combining exo-metabolomics, transcriptomics, and genome-scale metabolic network analysis and apply it to three antibody-producing Chinese Hamster Ovary cell lines to identify reprogramming features associated with high-producing clones and metabolic bottlenecks limiting product formation in an industrial bioprocess. Analysis of individual datatypes revealed a decreased nitrogenous byproduct secretion in high-producing clones and the topological changes in peripheral metabolic pathway expression associated with phase shifts. An integrated omics analysis in the context of the genome-scale metabolic model elucidated the differences in central metabolism and identified amino acid utilization bottlenecks limiting cell growth and antibody production that were not evident from exo-metabolomics or transcriptomics alone. Thus, we demonstrate the utility of a multi-omics characterization in providing an in-depth understanding of cellular metabolism, which is critical to efforts in cell engineering and bioprocess optimization.

Mitochondrial ATP generation is more proteome efficient than glycolysis

Metabolic efficiency profoundly influences organismal fitness. Nonphotosynthetic organisms, from yeast to mammals, derive usable energy primarily through glycolysis and respiration. Although respiration is more energy efficient, some cells favor glycolysis even when oxygen is available (aerobic glycolysis, Warburg effect). A leading explanation is that glycolysis is more efficient in terms of ATP production per unit mass of protein (that is, faster). Through quantitative flux analysis and proteomics, we find, however, that mitochondrial respiration is actually more proteome efficient than aerobic glycolysis. This is shown across yeast strains, T cells, cancer cells, and tissues and tumors in vivo. Instead of aerobic glycolysis being valuable for fast ATP production, it correlates with high glycolytic protein expression, which promotes hypoxic growth. Aerobic glycolytic yeasts do not excel at aerobic growth but outgrow respiratory cells during oxygen limitation. We accordingly propose that aerobic glycolysis emerges from cells maintaining a proteome conducive to both aerobic and hypoxic growth.

A hierarchy of metabolite exchanges in metabolic models of microbial species and communities

The metabolic network of an organism can be analyzed as a constraint-based model. This analysis can be biased, optimizing an objective such as growth rate, or unbiased, aiming to describe the full feasible space of metabolic fluxes through pathway analysis or random flux sampling. In particular, pathway analysis can decompose the flux space into fundamental and formally defined metabolic pathways. Unbiased methods scale poorly with network size due to combinatorial explosion, but a promising approach to improve scalability is to focus on metabolic subnetworks, e.g., cells' metabolite exchanges with each other and the environment, rather than the full metabolic networks. Here, we applied pathway enumeration and flux sampling to metabolite exchanges in microbial species and a microbial community, using models ranging from central carbon metabolism to genome-scale and focusing on pathway definitions that allow direct targeting of subnetworks such as metabolite exchanges (elementary conversion modes, elementary flux patterns, and minimal pathways). Enumerating growth-supporting metabolite exchanges, we found that metabolite exchanges from different pathway definitions were related through a hierarchy, and we show that this hierarchical relationship between pathways holds for metabolic networks and subnetworks more generally. Metabolite exchange frequencies, defined as the fraction of pathways in which each metabolite was exchanged, were similar across pathway definitions, with a few specific exchanges explaining large differences in pathway counts. This indicates that biological interpretation of predicted metabolite exchanges is robust to the choice of pathway definition, and it suggests strategies for more scalable pathway analysis. Our results also signal wider biological implications, facilitating detailed and interpretable analysis of metabolite exchanges and other subnetworks in fields such as metabolic engineering and synthetic biology.

Machine learning for the advancement of genome-scale metabolic modeling

Constraint-based modeling (CBM) has evolved as the core systems biology tool to map the interrelations between genotype, phenotype, and external environment. The recent advancement of high-throughput experimental approaches and multi-omics strategies has generated a plethora of new and precise information from wide-ranging biological domains. On the other hand, the continuously growing field of machine learning (ML) and its specialized branch of deep learning (DL) provide essential computational architectures for decoding complex and heterogeneous biological data. In recent years, both multi-omics and ML have assisted in the escalation of CBM. Condition-specific omics data, such as transcriptomics and proteomics, helped contextualize the model prediction while analyzing a particular phenotypic signature. At the same time, the advanced ML tools have eased the model reconstruction and analysis to increase the accuracy and prediction power. However, the development of these multi-disciplinary methodological frameworks mainly occurs independently, which limits the concatenation of biological knowledge from different domains. Hence, we have reviewed the potential of integrating multi-disciplinary tools and strategies from various fields, such as synthetic biology, CBM, omics, and ML, to explore the biochemical phenomenon beyond the conventional biological dogma. How the integrative knowledge of these intersected domains has improved bioengineering and biomedical applications has also been highlighted. We categorically explained the conventional genome-scale metabolic model (GEM) reconstruction tools and their improvement strategies through ML paradigms. Further, the crucial role of ML and DL in omics data restructuring for GEM development has also been briefly discussed. Finally, the case-study-based assessment of the state-of-the-art method for improving biomedical and metabolic engineering strategies has been elaborated. Therefore, this review demonstrates how integrating experimental and in silico strategies can help map the ever-expanding knowledge of biological systems driven by condition-specific cellular information. This multiview approach will elevate the application of ML-based CBM in the biomedical and bioengineering fields for the betterment of society and the environment.

Elucidation of triacylglycerol catabolism in Yarrowia lipolytica: How cells balance acetyl-CoA and excess reducing equivalents

Yarrowia lipolytica is an industrial yeast that can convert waste oil to value-added products. However, it is unclear how this yeast metabolizes lipid feedstocks, specifically triacylglycerol (TAG) substrates. This study used 13C-metabolic flux analysis (13C-MFA), genome-scale modeling, and transcriptomics analyses to investigate Y. lipolytica W29 growth with oleic acid, glycerol, and glucose. Transcriptomics data were used to guide 13C-MFA model construction and to validate the 13C-MFA results. The 13C-MFA data were then used to constrain a genome-scale model (GSM), which predicted Y. lipolytica fluxes, cofactor balance, and theoretical yields of terpene products. The three data sources provided new insights into cellular regulation during catabolism of glycerol and fatty acid components of TAG substrates, and how their consumption routes differ from glucose catabolism. We found that (1) over 80% of acetyl-CoA from oleic acid is processed through the glyoxylate shunt, a pathway that generates less CO2 compared to the TCA cycle, (2) the carnitine shuttle is a key regulator of the cytosolic acetyl-CoA pool in oleic acid and glycerol cultures, (3) the oxidative pentose phosphate pathway and mannitol cycle are key routes for NADPH generation, (4) the mannitol cycle and alternative oxidase activity help balance excess NADH generated from β-oxidation of oleic acid, and (5) asymmetrical gene expressions and GSM simulations of enzyme usage suggest an increased metabolic burden for oleic acid catabolism.

Genome-scale community modelling elucidates the metabolic interaction in Indian type-2 diabetic gut microbiota

Type-2 diabetes (T2D) is a rapidly growing multifactorial metabolic disorder that induces the onset of various diseases in the human body. The compositional and metabolic shift of the gut microbiota is a crucial factor behind T2D. Hence, gaining insight into the metabolic profile of the gut microbiota is essential for revealing their role in regulating the metabolism of T2D patients. Here, we have focused on the genome-scale community metabolic model reconstruction of crucial T2D-associated gut microbes. The model-based analysis of biochemical flux in T2D and healthy gut conditions showed distinct biochemical signatures and diverse metabolic interactions in the microbial community. The metabolic interactions encompass cross-feeding of short-chain fatty acids, amino acids, and vitamins among individual microbes within the community. In T2D conditions, a reduction in the metabolic flux of acetate, butyrate, vitamin B5, and bicarbonate was observed in the microbial community model, which can impact carbohydrate metabolism. The decline in butyrate levels is correlated with both insulin resistance and diminished glucose metabolism in T2D patients. Compared to the healthy gut, an overall reduction in glucose consumption and SCFA production flux was estimated in the T2D gut environment. Moreover, the decreased consumption profiles of branch chain amino acids (BCAAs) and aromatic amino acids (AAAs) in the T2D gut microbiota can be a distinct biomarker for T2D. Hence, the flux-level analysis of the microbial community model can provide insights into the metabolic reprogramming in diabetic gut microbiomes, which may be helpful in personalized therapeutics and diet design against T2D.

Metabolic model guided CRISPRi identifies a central role for phosphoglycerate mutase in Chlamydia trachomatis persistence

Upon nutrient starvation, Chlamydia trachomatis serovar L2 (CTL) shifts from its normal growth to a non-replicating form, termed persistence. It is unclear if persistence reflects an adaptive response or a lack thereof. To understand this, transcriptomics data were collected for CTL grown under nutrient-replete and nutrient-starved conditions. Applying K-means clustering on transcriptomics data revealed a global transcriptomic rewiring of CTL under stress conditions in the absence of any canonical global stress regulator. This is consistent with previous data that suggested that CTL's stress response is due to a lack of an adaptive response mechanism. To investigate the impact of this on CTL metabolism, we reconstructed a genome-scale metabolic model of CTL (iCTL278) and contextualized it with the collected transcriptomics data. Using the metabolic bottleneck analysis on contextualized iCTL278, we observed that phosphoglycerate mutase (pgm) regulates the entry of CTL to the persistence state. Our data indicate that pgm has the highest thermodynamics driving force and lowest enzymatic cost. Furthermore, CRISPRi-driven knockdown of pgm in the presence or absence of tryptophan revealed the importance of this gene in modulating persistence. Hence, this work, for the first time, introduces thermodynamics and enzyme cost as tools to gain a deeper understanding on CTL persistence.

IMPORTANCE

This study uses a metabolic model to investigate factors that contribute to the persistence of Chlamydia trachomatis serovar L2 (CTL) under tryptophan and iron starvation conditions. As CTL lacks many canonical transcriptional regulators, the model was used to assess two prevailing hypotheses on persistence-that the chlamydial response to nutrient starvation represents a passive response due to the lack of regulators or that it is an active response by the bacterium. K-means clustering of stress-induced transcriptomics data revealed striking evidence in favor of the lack of adaptive (i.e., a passive) response. To find the metabolic signature of this, metabolic modeling pin-pointed pgm as a potential regulator of persistence. Thermodynamic driving force, enzyme cost, and CRISPRi knockdown of pgm supported this finding. Overall, this work introduces thermodynamic driving force and enzyme cost as a tool to understand chlamydial persistence, demonstrating how systems biology-guided CRISPRi can unravel complex bacterial phenomena.

Panera: An innovative framework for surmounting uncertainty in microbial community modeling using pan-genera metabolic models

Utilization of 16S rRNA data in constraint-based modeling to characterize microbial communities confronts a major hurdle of lack of species-level resolution, impeding the construction of community models. We introduce "Panera," an innovative framework designed to model communities under this uncertainty and yet perform metabolic inferences using pan-genus metabolic models (PGMMs). We demonstrated PGMMs' utility for comprehending the metabolic capabilities of a genus and in characterizing community models using amplicon data. The unique, adaptable nature of PGMMs unlocks their potential in building hybrid communities, combining genome-scale metabolic models (GSMMs) and PGMMs. Notably, these models provide predictions comparable to the standard GSMM-based community models, while achieving a nearly 46% reduction in error compared to the genus model-based communities. In essence, "Panera" presents a potent and effective approach to aid in metabolic modeling by enabling robust predictions of community metabolic potential when dealing with amplicon data, and offers insights into genus-level metabolic landscapes.

Bioenergetic stress potentiates antimicrobial resistance and persistence

Antimicrobial resistance (AMR) is a global health crisis and there is an urgent need to better understand AMR mechanisms. Antibiotic treatment alters several aspects of bacterial physiology, including increased ATP utilization, carbon metabolism, and reactive oxygen species (ROS) formation. However, how the "bioenergetic stress" induced by increased ATP utilization affects treatment outcomes is unknown. Here we utilized a synthetic biology approach to study the direct effects of bioenergetic stress on antibiotic efficacy. We engineered a genetic system that constitutively hydrolyzes ATP or NADH in Escherichia coli. We found that bioenergetic stress potentiates AMR evolution via enhanced ROS production, mutagenic break repair, and transcription-coupled repair. We also find that bioenergetic stress potentiates antimicrobial persistence via potentiated stringent response activation. We propose a unifying model that antibiotic-induced antimicrobial resistance and persistence is caused by antibiotic-induced. This has important implications for preventing or curbing the spread of AMR infections.

Semi-Automatic Detection of Errors in Genome-Scale Metabolic Models

Background

Genome-Scale Metabolic Models (GSMMs) are used for numerous tasks requiring computational estimates of metabolic fluxes, from predicting novel drug targets to engineering microbes to produce valuable compounds. A key limiting step in most applications of GSMMs is ensuring their representation of the target organism's metabolism is complete and accurate. Identifying and visualizing errors in GSMMs is complicated by the fact that they contain thousands of densely interconnected reactions. Furthermore, many errors in GSMMs only become apparent when considering pathways of connected reactions collectively, as opposed to examining reactions individually.

Results

We present Metabolic Accuracy Check and Analysis Workflow (MACAW), a collection of algorithms for detecting errors in GSMMs. The relative frequencies of errors we detect in manually curated GSMMs appear to reflect the different approaches used to curate them. Changing the method used to automatically create a GSMM from a particular organism's genome can have a larger impact on the kinds of errors in the resulting GSMM than using the same method with a different organism's genome. Our algorithms are particularly capable of identifying errors that are only apparent at the pathway level, including loops, and nontrivial cases of dead ends.

Conclusions

MACAW is capable of identifying inaccuracies of varying severity in a wide range of GSMMs. Correcting these errors can measurably improve the predictive capacity of a GSMM. The relative prevalence of each type of error we identify in a large collection of GSMMs could help shape future efforts for further automation of error correction and GSMM creation.

Development of a genome-scale metabolic model for the lager hybrid yeast S. pastorianus to understand the evolution of metabolic pathways in industrial settings

In silico tools such as genome-scale metabolic models have shown to be powerful for metabolic engineering of microorganisms. Saccharomyces pastorianus is a complex aneuploid hybrid between the mesophilic Saccharomyces cerevisiae and the cold-tolerant Saccharomyces eubayanus. This species is of biotechnological importance because it is the primary yeast used in lager beer fermentation and is also a key model for studying the evolution of hybrid genomes, including expression pattern of ortholog genes, composition of protein complexes, and phenotypic plasticity. Here, we created the iSP_1513 GSMM for S. pastorianus CBS1513 to allow top-down computational approaches to predict the evolution of metabolic pathways and to aid strain optimization in production processes. The iSP_1513 comprises 4,062 reactions, 1,808 alleles, and 2,747 metabolites, and takes into account the functional redundancy in the gene-protein-reaction rule caused by the presence of orthologous genes. Moreover, a universal algorithm to constrain GSMM reactions using transcriptome data was developed as a python library and enabled the integration of temperature as parameter. Essentiality data sets, growth data on various carbohydrates and volatile metabolites secretion were used to validate the model and showed the potential of media engineering to improve specific flavor compounds. The iSP_1513 also highlighted the different contributions of the parental sub-genomes to the oxidative and non-oxidative parts of the pentose phosphate pathway. Overall, the iSP_1513 GSMM represent an important step toward understanding the metabolic capabilities, evolutionary trajectories, and adaptation potential of S. pastorianus in different industrial settings.

IMPORTANCE

Genome-scale metabolic models (GSMM) have been successfully applied to predict cellular behavior and design cell factories in several model organisms, but no models to date are currently available for hybrid species due to their more complex genetics and general lack of molecular data. In this study, we generated a bespoke GSMM, iSP_1513, for this industrial aneuploid hybrid Saccharomyces pastorianus, which takes into account the aneuploidy and functional redundancy from orthologous parental alleles. This model will (i) help understand the metabolic capabilities and adaptive potential of S. pastorianus (domestication processes), (ii) aid top-down predictions for strain development (industrial biotechnology), and (iii) allow predictions of evolutionary trajectories of metabolic pathways in aneuploid hybrids (evolutionary genetics).

Genome-scale metabolic modeling of Thraustochytrium sp. RT2316-16: Effects of nutrients on metabolism

Marine thraustochytrids produce metabolically important lipids such as the long-chain omega-3 polyunsaturated fatty acids, carotenoids, and sterols. The growth and lipid production in thraustochytrids depends on the composition of the culture medium that often contains yeast extract as a source of amino acids. This work discusses the effects of individual amino acids provided in the culture medium as the only source of nitrogen, on the production of biomass and lipids by the thraustochytrid Thraustochytrium sp. RT2316-16. A reconstructed metabolic network based on the annotated genome of RT2316-16 in combination with flux balance analysis was used to explain the observed growth and consumption of the nutrients. The culture kinetic parameters estimated from the experimental data were used to constrain the flux via the nutrient consumption rates and the specific growth rate of the triacylglycerol-free biomass in the genome-scale metabolic model (GEM) to predict the specific rate of ATP production for cell maintenance. A relationship was identified between the specific rate of ATP production for maintenance and the specific rate of glucose consumption. The GEM and the derived relationship for the production of ATP for maintenance were used in linear optimization problems, to successfully predict the specific growth rate of RT2316-16 in different experimental conditions.

Accounting for NAD Concentrations in Genome-Scale Metabolic Models Captures Important Metabolic Alterations in NAD-Depleted Systems

Nicotinamide adenine dinucleotide (NAD) is a ubiquitous molecule found within all cells, acting as a crucial coenzyme in numerous metabolic reactions. It plays a vital role in energy metabolism, cellular signaling, and DNA repair. Notably, NAD levels decline naturally with age, and this decline is associated with the development of various age-related diseases. Despite this established link, current genome-scale metabolic models, which offer powerful tools for understanding cellular metabolism, do not account for the dynamic changes in NAD concentration. This impedes our understanding of a fluctuating NAD level's impact on cellular metabolism and its contribution to age-related pathologies. To bridge this gap in our knowledge, we have devised a novel method that integrates altered NAD concentration into genome-scale models of human metabolism. This approach allows us to accurately reflect the changes in fatty acid metabolism, glycolysis, and oxidative phosphorylation observed experimentally in an engineered human cell line with a compromised level of subcellular NAD.

Models and molecular mechanisms for trade-offs in the context of metabolism

Accumulating evidence for trade-offs involving metabolic traits has demonstrated their importance in the evolution of organisms. Metabolic models with different levels of complexity have already been considered when investigating mechanisms that explain various metabolic trade-offs. Here we provide a systematic review of modelling approaches that have been used to study and explain trade-offs between: (i) the kinetic properties of individual enzymes, (ii) rates of metabolic reactions, (iii) the rate and yield of metabolic pathways and networks, (iv) different metabolic objectives in single organisms and in metabolic communities, and (v) metabolic concentrations. In providing insights into the mechanisms underlying these five types of metabolic trade-offs obtained from constraint-based metabolic modelling, we emphasize the relationship of metabolic trade-offs to the classical black box Y-model that provides a conceptual explanation for resource acquisition-allocation trade-offs. In addition, we identify several pressing concerns and offer a perspective for future research in the identification and manipulation of metabolic trade-offs by relying on the toolbox provided by constraint-based metabolic modelling for single organisms and microbial communities.

A Multidisciplinary Hyper-Modeling Scheme in Personalized In Silico Oncology: Coupling Cell Kinetics with Metabolism, Signaling Networks, and Biomechanics as Plug-In Component Models of a Cancer Digital Twin

The massive amount of human biological, imaging, and clinical data produced by multiple and diverse sources necessitates integrative modeling approaches able to summarize all this information into answers to specific clinical questions. In this paper, we present a hypermodeling scheme able to combine models of diverse cancer aspects regardless of their underlying method or scale. Describing tissue-scale cancer cell proliferation, biomechanical tumor growth, nutrient transport, genomic-scale aberrant cancer cell metabolism, and cell-signaling pathways that regulate the cellular response to therapy, the hypermodel integrates mutation, miRNA expression, imaging, and clinical data. The constituting hypomodels, as well as their orchestration and links, are described. Two specific cancer types, Wilms tumor (nephroblastoma) and non-small cell lung cancer, are addressed as proof-of-concept study cases. Personalized simulations of the actual anatomy of a patient have been conducted. The hypermodel has also been applied to predict tumor control after radiotherapy and the relationship between tumor proliferative activity and response to neoadjuvant chemotherapy. Our innovative hypermodel holds promise as a digital twin-based clinical decision support system and as the core of future in silico trial platforms, although additional retrospective adaptation and validation are necessary.

Metabolic modeling of Halomonas campaniensis improves polyhydroxybutyrate production under nitrogen limitation

Poly-hydroxybutyrate (PHB) is an environmentally friendly alternative for conventional fossil fuel-based plastics that is produced by various microorganisms. Large-scale PHB production is challenging due to the comparatively higher biomanufacturing costs. A PHB overproducer is the haloalkaliphilic bacterium Halomonas campaniensis, which has low nutritional requirements and can grow in cultures with high salt concentrations, rendering it resistant to contamination. Despite its virtues, the metabolic capabilities of H. campaniensis as well as the limitations hindering higher PHB production remain poorly studied. To address this limitation, we present HaloGEM, the first high-quality genome-scale metabolic network reconstruction, which encompasses 888 genes, 1528 reactions (1257 gene-associated), and 1274 metabolites. HaloGEM not only displays excellent agreement with previous growth data and experiments from this study, but it also revealed nitrogen as a limiting nutrient when growing aerobically under high salt concentrations using glucose as carbon source. Among different nitrogen source mixtures for optimal growth, HaloGEM predicted glutamate and arginine as a promising mixture producing increases of 54.2% and 153.4% in the biomass yield and PHB titer, respectively. Furthermore, the model was used to predict genetic interventions for increasing PHB yield, which were consistent with the rationale of previously reported strategies. Overall, the presented reconstruction advances our understanding of the metabolic capabilities of H. campaniensis for rationally engineering this next-generation industrial biotechnology platform. KEY POINTS: A comprehensive genome-scale metabolic reconstruction of H. campaniensis was developed. Experiments and simulations predict N limitation in minimal media under aerobiosis. In silico media design increased experimental biomass yield and PHB titer.

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.

Prediction and integration of metabolite-protein interactions with genome-scale metabolic models

Metabolites, as small molecules, can act not only as substrates to enzymes, but also as effectors of activity of proteins with different functions, thereby affecting various cellular processes. While several experimental techniques have started to catalogue the metabolite-protein interactions (MPIs) present in different cellular contexts, characterizing the functional relevance of MPIs remains a challenging problem. Computational approaches from the constrained-based modeling framework allow for predicting MPIs and integrating their effects in the in silico analysis of metabolic and physiological phenotypes, like cell growth. Here, we provide a classification of all existing constraint-based approaches that predict and integrate MPIs using genome-scale metabolic networks as input. In addition, we benchmark the performance of the approaches to predict MPIs in a comparative study using different features extracted from the model structure and predicted metabolic phenotypes with the state-of-the-art metabolic networks of Escherichia coli and Saccharomyces cerevisiae. Lastly, we provide an outlook for future, feasible directions to expand the consideration of MPIs in constraint-based modeling approaches with wide biotechnological applications.

Genome scale metabolic network modelling for metabolic profile predictions

Metabolic profiling (metabolomics) aims at measuring small molecules (metabolites) in complex samples like blood or urine for human health studies. While biomarker-based assessment often relies on a single molecule, metabolic profiling combines several metabolites to create a more complex and more specific fingerprint of the disease. However, in contrast to genomics, there is no unique metabolomics setup able to measure the entire metabolome. This challenge leads to tedious and resource consuming preliminary studies to be able to design the right metabolomics experiment. In that context, computer assisted metabolic profiling can be of strong added value to design metabolomics studies more quickly and efficiently. We propose a constraint-based modelling approach which predicts in silico profiles of metabolites that are more likely to be differentially abundant under a given metabolic perturbation (e.g. due to a genetic disease), using flux simulation. In genome-scale metabolic networks, the fluxes of exchange reactions, also known as the flow of metabolites through their external transport reactions, can be simulated and compared between control and disease conditions in order to calculate changes in metabolite import and export. These import/export flux differences would be expected to induce changes in circulating biofluid levels of those metabolites, which can then be interpreted as potential biomarkers or metabolites of interest. In this study, we present SAMBA (SAMpling Biomarker Analysis), an approach which simulates fluxes in exchange reactions following a metabolic perturbation using random sampling, compares the simulated flux distributions between the baseline and modulated conditions, and ranks predicted differentially exchanged metabolites as potential biomarkers for the perturbation. We show that there is a good fit between simulated metabolic exchange profiles and experimental differential metabolites detected in plasma, such as patient data from the disease database OMIM, and metabolic trait-SNP associations found in mGWAS studies. These biomarker recommendations can provide insight into the underlying mechanism or metabolic pathway perturbation lying behind observed metabolite differential abundances, and suggest new metabolites as potential avenues for further experimental analyses.

Adjusting for false discoveries in constraint-based differential metabolic flux analysis

One of the critical steps to characterize metabolic alterations in multifactorial diseases, as well as their heterogeneity across different patients, is the identification of reactions that exhibit significantly different usage (or flux) between cohorts. However, since metabolic fluxes cannot be determined directly, researchers typically use constraint-based metabolic network models, customized on post-genomics datasets. The use of random sampling within the feasible region of metabolic networks is becoming more prevalent for comparing these networks. While many algorithms have been proposed and compared for efficiently and uniformly sampling the feasible region of metabolic networks, their impact on the risk of making false discoveries when comparing different samples has not been investigated yet, and no sampling strategy has been so far specifically designed to mitigate the problem. To be able to precisely assess the False Discovery Rate (FDR), in this work we compared different samples obtained from the very same metabolic model. We compared the FDR obtained for different model scales, sample sizes, parameters of the sampling algorithm, and strategies to filter out non-significant variations. To be able to compare the largely used hit-and-run strategy with the much less investigated corner-based strategy, we first assessed the intrinsic capability of current corner-based algorithms and of a newly proposed one to visit all vertices of a constraint-based region. We show that false discoveries can occur at high rates even for large samples of small-scale networks. However, we demonstrate that a statistical test based on the empirical null distribution of Kullback-Leibler divergence can effectively correct for false discoveries. We also show that our proposed corner-based algorithm is more efficient than state-of-the-art alternatives and much less prone to false discoveries than hit-and-run strategies. We report that the differences in the marginal distributions obtained with the two strategies are related to but not fully explained by differences in sample standard deviation, as previously thought. Overall, our study provides insights into the impact of sampling strategies on FDR in metabolic network analysis and offers new guidelines for more robust and reproducible analyses.

Targeting the Bottlenecks in Levan Biosynthesis Pathway in Bacillus subtilis and Strain Optimization by Computational Modeling and Omics Integration

Levan is a fructan polymer with many industrial applications such as the formulation of hydrogels, drug delivery, and wound healing, among others. To this end, metabolic systems engineering is a valuable method to improve the yield of a specific metabolite in a wide range of bacterial and eukaryotic organisms. In this study, we report a systems biology approach integrating genomics data for the Bacillus subtilis model, wherein the metabolic pathway for levan biosynthesis is unpacked. We analyzed a revised genome-scale enzyme-constrained metabolic model (ecGEM) and performed simulations to increase levan biopolymer production capacity in B. subtilis. We used the model ec_iYO844_lvn to (1) identify the essential genes and bottlenecks in levan production, and (2) specifically design an engineered B. subtilis strain capable of producing higher levan yields. The FBA and FVA analysis showed the maximal growth rate of the organism up to 0.624 hr-1 at 20 mmol gDw-1 hr-1 of sucrose intake. Gene knockout analyses were performed to identify gene knockout targets to increase the levan flux in B. subtilis. Importantly, we found that the pgk and ctaD genes are the two target genes for the knockout. The perturbation of these two genes has flux gains for levan production reactions with 1.3- and 1.4-fold the relative flux span in the mutant strains, respectively, compared to the wild type. In all, this work identifies the bottlenecks in the production of levan and possible ways to overcome them. Our results provide deeper insights on the bacterium's physiology and new avenues for strain engineering.

FastKnock: an efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed by making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies (with a predefined maximum number of reaction deletions) for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using various Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple- and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more beneficial and important practical solutions. The availability of all the solutions provides the opportunity to further characterize, rank and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock .

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Model validation and selection in metabolic flux analysis and flux balance analysis

13C-Metabolic Flux Analysis (13C-MFA) and Flux Balance Analysis (FBA) are widely used to investigate the operation of biochemical networks in both biological and biotechnological research. Both methods use metabolic reaction network models of metabolism operating at steady state so that reaction rates (fluxes) and the levels of metabolic intermediates are constrained to be invariant. They provide estimated (MFA) or predicted (FBA) values of the fluxes through the network in vivo, which cannot be measured directly. These fluxes can shed light on basic biology and have been successfully used to inform metabolic engineering strategies. Several approaches have been taken to test the reliability of estimates and predictions from constraint-based methods and to compare alternative model architectures. Despite advances in other areas of the statistical evaluation of metabolic models, such as the quantification of flux estimate uncertainty, validation and model selection methods have been underappreciated and underexplored. We review the history and state-of-the-art in constraint-based metabolic model validation and model selection. Applications and limitations of the χ2 -test of goodness-of-fit, the most widely used quantitative validation and selection approach in 13C-MFA, are discussed, and complementary and alternative forms of validation and selection are proposed. A combined model validation and selection framework for 13C-MFA incorporating metabolite pool size information that leverages new developments in the field is presented and advocated for. Finally, we discuss how adopting robust validation and selection procedures can enhance confidence in constraint-based modeling as a whole and ultimately facilitate more widespread use of FBA in biotechnology.

Machine Learning and Metabolic Model Guided CRISPRi Reveals a Central Role for Phosphoglycerate Mutase in Chlamydia trachomatis Persistence

Upon nutrient starvation, Chlamydia trachomatis serovar L2 (CTL) shifts from its normal growth to a non-replicating form, termed persistence. It is unclear if persistence is an adaptive response or lack of it. To understand that transcriptomics data were collected for nutrient-sufficient and nutrient-starved CTL. Applying machine learning approaches on transcriptomics data revealed a global transcriptomic rewiring of CTL under stress conditions without having any global stress regulator. This indicated that CTL's stress response is due to lack of an adaptive response mechanism. To investigate the impact of this on CTL metabolism, we reconstructed a genome-scale metabolic model of CTL (iCTL278) and contextualized it with the collected transcriptomics data. Using the metabolic bottleneck analysis on contextualized iCTL278, we observed phosphoglycerate mutase (pgm) regulates the entry of CTL to the persistence. Later, pgm was found to have the highest thermodynamics driving force and lowest enzymatic cost. Furthermore, CRISPRi-driven knockdown of pgm and tryptophan starvation experiments revealed the importance of this gene in inducing persistence. Hence, this work, for the first time, introduced thermodynamics and enzyme-cost as tools to gain deeper understanding on CTL persistence.

Applications of genome-scale metabolic models to investigate microbial metabolic adaptations in response to genetic or environmental perturbations

Environmental perturbations are encountered by microorganisms regularly and will require metabolic adaptations to ensure an organism can survive in the newly presenting conditions. In order to study the mechanisms of metabolic adaptation in such conditions, various experimental and computational approaches have been used. Genome-scale metabolic models (GEMs) are one of the most powerful approaches to study metabolism, providing a platform to study the systems level adaptations of an organism to different environments which could otherwise be infeasible experimentally. In this review, we are describing the application of GEMs in understanding how microbes reprogram their metabolic system as a result of environmental variation. In particular, we provide the details of metabolic model reconstruction approaches, various algorithms and tools for model simulation, consequences of genetic perturbations, integration of '-omics' datasets for creating context-specific models and their application in studying metabolic adaptation due to the change in environmental conditions.

Inferring mitochondrial and cytosolic metabolism by coupling isotope tracing and deconvolution

The inability to inspect metabolic activities within distinct subcellular compartments has been a major barrier to our understanding of eukaryotic cell metabolism. Previous work addressed this challenge by analyzing metabolism in isolated organelles, which grossly bias metabolic activity. Here, we describe a method for inferring physiological metabolic fluxes and metabolite concentrations in mitochondria and cytosol based on isotope tracing experiments performed with intact cells. This is made possible by computational deconvolution of metabolite isotopic labeling patterns and concentrations into cytosolic and mitochondrial counterparts, coupled with metabolic and thermodynamic modelling. Our approach lowers the uncertainty regarding compartmentalized fluxes and concentrations by one and three orders of magnitude compared to existing modelling approaches, respectively. We derive a quantitative view of mitochondrial and cytosolic metabolic activities in central carbon metabolism across cultured cell lines without performing cell fractionation, finding major variability in compartmentalized malate-aspartate shuttle fluxes. We expect our approach for inferring metabolism at a subcellular resolution to be instrumental for a variety of studies of metabolic dysfunction in human disease and for bioengineering.

Reconstruction of a Genome-Scale Metabolic Model of Scenedesmus obliquus and Its Application for Lipid Production under Three Trophic Modes

Green microalgae have emerged as beneficial feedstocks for biofuel production. A systems-level understanding of the biochemical network is needed to harness the microalgal metabolic capacity for bioproduction. Genome-scale metabolic modeling (GEM) showed immense potential in rational metabolic engineering, utilizing biochemical flux distribution analysis. Here, we report the first GEM for the green microalga, Scenedesmus obliquus (iAR632), a promising biodiesel feedstock with high lipid-storing capability. iAR632 comprises 1467 reactions, 734 metabolites, and 632 genes distributed among 7 compartments. The model was optimized under three different trophic modes of microalgal cultivation, i.e., autotrophy, mixotrophy, and heterotrophy. The robustness of the reconstructed network was confirmed by analyzing its sensitivity to the biomass components. Pathway-level flux profiles were analyzed, and significant flux space expansion was noticed majorly in reactions associated with lipid biosynthesis. In agreement with the experimental observation, iAR632 predicted about 3.8-fold increased biomass and almost 4-fold higher lipid under mixotrophy than the other trophic modes. Thus, the assessment of the condition-specific metabolic flux distribution of iAR632 suggested that mixotrophy is the preferred cultivation condition for improved microalgal growth and lipid production. Overall, the reconstructed GEM and subsequent analyses will provide a systematic framework for developing model-driven strategies to improve microalgal bioproduction.

A Markov constraint to uniquely identify elementary flux mode weights in unimolecular metabolic networks

Elementary flux modes (EFMs) are minimal, steady state pathways characterizing a flux network. Fundamentally, all steady state fluxes in a network are decomposable into a linear combination of EFMs. While there is typically no unique set of EFM weights that reconstructs these fluxes, several optimization-based methods have been proposed to constrain the solution space by enforcing some notion of parsimony. However, it has long been recognized that optimization-based approaches may fail to uniquely identify EFM weights and return different feasible solutions across objective functions and solvers. Here we show that, for flux networks only involving single molecule transformations, these problems can be avoided by imposing a Markovian constraint on EFM weights. Our Markovian constraint guarantees a unique solution to the flux decomposition problem, and that solution is arguably more biophysically plausible than other solutions. We describe an algorithm for computing Markovian EFM weights via steady state analysis of a certain discrete-time Markov chain, based on the flux network, which we call the cycle-history Markov chain. We demonstrate our method with a differential analysis of EFM activity in a lipid metabolic network comparing healthy and Alzheimer's disease patients. Our method is the first to uniquely decompose steady state fluxes into EFM weights for any unimolecular metabolic network.

A reactive species reactions module for integration into genome-scale metabolic models for improved insights: Application to cancer

Reactive species (RS) play significant roles in many disease contexts. Despite their crucial roles in diseases including cancer, the RS are not adequately modeled in the genome-scale metabolic (GSM) models, which are used to understand cell metabolism in disease contexts. We have developed a scalable RS reactions module that can be integrated with any Recon 3D-derived human metabolic model, or after fine-tuning, with any metabolic model. With RS-integration, the GSM models of three cancers (basal-like triple negative breast cancer (TNBC), high grade serous ovarian carcinoma (HGSOC) and colorectal cancer (CRC)) built from Recon 3D, precisely highlighted the increases/decreases in fluxes (dysregulation) occurring in important pathways of these cancers. These dysregulations were not prominent in the standard cancer models without the RS module. Further, the results from these RS-integrated cancer GSM models suggest the following decreasing order in the ease of ferroptosis-targeting to treat the cancers: TNBC > HGSOC > CRC.

Unveiling abundance-dependent metabolic phenotypes of microbial communities

IMPORTANCE

In nature, organisms live in communities and not as isolated species, and their interactions provide a source of resilience to environmental disturbances. Despite their importance in ecology, human health, and industry, understanding how organisms interact in different environments remains an open question. In this work, we provide a novel approach that, only using genomic information, studies the metabolic phenotype exhibited by communities, where the exploration of suboptimal growth flux distributions and the composition of a community allows to unveil its capacity to respond to environmental changes, shedding light of the degrees of metabolic plasticity inherent to the community.

Model-driven insights into the effects of temperature on metabolism

Temperature affects cellular processes at different spatiotemporal scales, and identifying the genetic and molecular mechanisms underlying temperature responses paves the way to develop approaches for mitigating the effects of future climate scenarios. A systems view of the effects of temperature on cellular physiology can be obtained by focusing on metabolism since: (i) its functions depend on transcription and translation and (ii) its outcomes support organisms' development, growth, and reproduction. Here we provide a systematic review of modelling efforts directed at investigating temperature effects on properties of single biochemical reactions, system-level traits, metabolic subsystems, and whole-cell metabolism across different prokaryotes and eukaryotes. We compare and contrast computational approaches and theories that facilitate modelling of temperature effects on key properties of enzymes and their consideration in constraint-based as well as kinetic models of metabolism. In addition, we provide a summary of insights from computational approaches, facilitating integration of omics data from temperature-modulated experiments with models of metabolic networks, and review the resulting biotechnological applications. Lastly, we provide a perspective on how different types of metabolic modelling can profit from developments in machine learning and models of different cellular layers to improve model-driven insights into the effects of temperature relevant for biotechnological applications.

A genome-scale metabolic reconstruction of soybean and Bradyrhizobium diazoefficiens reveals the cost-benefit of nitrogen fixation

Nitrogen-fixing symbioses allow legumes to thrive in nitrogen-poor soils at the cost of diverting some photoassimilate to their microsymbionts. Effort is being made to bioengineer nitrogen fixation into nonleguminous crops. This requires a quantitative understanding of its energetic costs and the links between metabolic variations and symbiotic efficiency. A whole-plant metabolic model for soybean (Glycine max) with its associated microsymbiont Bradyrhizobium diazoefficiens was developed and applied to predict the cost-benefit of nitrogen fixation with varying soil nitrogen availability. The model predicted a nitrogen-fixation cost of c. 4.13 g C g-1 N, which when implemented into a crop scale model, translated to a grain yield reduction of 27% compared with a non-nodulating plant receiving its nitrogen from the soil. Considering the lower nitrogen content of cereals, the yield cost to a hypothetical N-fixing cereal is predicted to be less than half that of soybean. Soybean growth was predicted to be c. 5% greater when the nodule nitrogen export products were amides versus ureides. This is the first metabolic reconstruction in a tropical crop species that simulates the entire plant and nodule metabolism. Going forward, this model will serve as a tool to investigate carbon use efficiency and key mechanisms within N-fixing symbiosis in a tropical species forming determinate nodules.

Maximizing microbial bioproduction from sustainable carbon sources using iterative systems engineering

Maximizing the production of heterologous biomolecules is a complex problem that can be addressed with a systems-level understanding of cellular metabolism and regulation. Specifically, growth-coupling approaches can increase product titers and yields and also enhance production rates. However, implementing these methods for non-canonical carbon streams is challenging due to gaps in metabolic models. Over four design-build-test-learn cycles, we rewire Pseudomonas putida KT2440 for growth-coupled production of indigoidine from para-coumarate. We explore 4,114 potential growth-coupling solutions and refine one design through laboratory evolution and ensemble data-driven methods. The final growth-coupled strain produces 7.3 g/L indigoidine at 77% maximum theoretical yield in para-coumarate minimal medium. The iterative use of growth-coupling designs and functional genomics with experimental validation was highly effective and agnostic to specific hosts, carbon streams, and final products and thus generalizable across many systems.

Optimal energy and redox metabolism in the cyanobacterium Synechocystis sp. PCC 6803

Understanding energy and redox homeostasis and carbon partitioning is crucial for systems metabolic engineering of cell factories. Carbon metabolism alone cannot achieve maximal accumulation of metabolites in production hosts, since an efficient production of target molecules requires energy and redox balance, in addition to carbon flow. The interplay between cofactor regeneration and heterologous production in photosynthetic microorganisms is not fully explored. To investigate the optimality of energy and redox metabolism, while overproducing alkenes-isobutene, isoprene, ethylene and 1-undecene, in the cyanobacterium Synechocystis sp. PCC 6803, we applied stoichiometric metabolic modelling. Our network-wide analysis indicates that the rate of NAD(P)H regeneration, rather than of ATP, controls ATP/NADPH ratio, and thereby bioproduction. The simulation also implies that energy and redox balance is interconnected with carbon and nitrogen metabolism. Furthermore, we show that an auxiliary pathway, composed of serine, one-carbon and glycine metabolism, supports cellular redox homeostasis and ATP cycling. The study revealed non-intuitive metabolic pathways required to enhance alkene production, which are mainly driven by a few key reactions carrying a high flux. We envision that the presented comparative in-silico metabolic analysis will guide the rational design of Synechocystis as a photobiological production platform of target chemicals.

Mathematical reconstruction of the metabolic network in an in-vitro multiple myeloma model

It is increasingly apparent that cancer cells, in addition to remodelling their metabolism to survive and proliferate, adapt and manipulate the metabolism of other cells. This property may be a telling sign that pre-clinical tumour metabolism studies exclusively utilising in-vitro mono-culture models could prove to be limited for uncovering novel metabolic targets able to translate into clinical therapies. Although this is increasingly recognised, and work towards addressing the issue is becoming routinary much remains poorly understood. For instance, knowledge regarding the biochemical mechanisms through which cancer cells manipulate non-cancerous cell metabolism, and the subsequent impact on their survival and proliferation remains limited. Additionally, the variations in these processes across different cancer types and progression stages, and their implications for therapy, also remain largely unexplored. This study employs an interdisciplinary approach that leverages the predictive power of mathematical modelling to enrich experimental findings. We develop a functional multicellular in-silico model that facilitates the qualitative and quantitative analysis of the metabolic network spawned by an in-vitro co-culture model of bone marrow mesenchymal stem- and myeloma cell lines. To procure this model, we devised a bespoke human genome constraint-based reconstruction workflow that combines aspects from the legacy mCADRE & Metabotools algorithms, the novel redHuman algorithm, along with 13C-metabolic flux analysis. Our workflow transforms the latest human metabolic network matrix (Recon3D) into two cell-specific models coupled with a metabolic network spanning a shared growth medium. When cross-validating our in-silico model against the in-vitro model, we found that the in-silico model successfully reproduces vital metabolic behaviours of its in-vitro counterpart; results include cell growth predictions, respiration rates, as well as support for observations which suggest cross-shuttling of redox-active metabolites between cells.

The first multi-tissue genome-scale metabolic model of a woody plant highlights suberin biosynthesis pathways in Quercus suber

Over the last decade, genome-scale metabolic models have been increasingly used to study plant metabolic behaviour at the tissue and multi-tissue level under different environmental conditions. Quercus suber, also known as the cork oak tree, is one of the most important forest communities of the Mediterranean/Iberian region. In this work, we present the genome-scale metabolic model of the Q. suber (iEC7871). The metabolic model comprises 7871 genes, 6231 reactions, and 6481 metabolites across eight compartments. Transcriptomics data was integrated into the model to obtain tissue-specific models for the leaf, inner bark, and phellogen, with specific biomass compositions. The tissue-specific models were merged into a diel multi-tissue metabolic model to predict interactions among the three tissues at the light and dark phases. The metabolic models were also used to analyse the pathways associated with the synthesis of suberin monomers, namely the acyl-lipids, phenylpropanoids, isoprenoids, and flavonoids production. The models developed in this work provide a systematic overview of the metabolism of Q. suber, including its secondary metabolism pathways and cork formation.

Maximizing multi-reaction dependencies provides more accurate and precise predictions of intracellular fluxes than the principle of parsimony

Intracellular fluxes represent a joint outcome of cellular transcription and translation and reflect the availability and usage of nutrients from the environment. While approaches from the constraint-based metabolic framework can accurately predict cellular phenotypes, such as growth and exchange rates with the environment, accurate prediction of intracellular fluxes remains a pressing problem. Parsimonious flux balance analysis (pFBA) has become an approach of choice to predict intracellular fluxes by employing the principle of efficient usage of protein resources. Nevertheless, comparative analyses of intracellular flux predictions from pFBA against fluxes estimated from labeling experiments remain scarce. Here, we posited that steady-state flux distributions derived from the principle of maximizing multi-reaction dependencies are of improved accuracy and precision than those resulting from pFBA. To this end, we designed a constraint-based approach, termed complex-balanced FBA (cbFBA), to predict steady-state flux distributions that support the given specific growth rate and exchange fluxes. We showed that the steady-state flux distributions resulting from cbFBA in comparison to pFBA show better agreement with experimentally measured fluxes from 17 Escherichia coli strains and are more precise, due to the smaller space of alternative solutions. We also showed that the same principle holds in eukaryotes by comparing the predictions of pFBA and cbFBA against experimentally derived steady-state flux distributions from 26 knock-out mutants of Saccharomyces cerevisiae. Furthermore, our results showed that intracellular fluxes predicted by cbFBA provide better support for the principle of minimizing metabolic adjustment between mutants and wild types. Together, our findings point that other principles that consider the dynamics and coordination of steady states may govern the distribution of intracellular fluxes.

Uses of Multi-Objective Flux Analysis for Optimization of Microbial Production of Secondary Metabolites

Secondary metabolites are not essential for the growth of microorganisms, but they play a critical role in how microbes interact with their surroundings. In addition to this important ecological role, secondary metabolites also have a variety of agricultural, medicinal, and industrial uses, and thus the examination of secondary metabolism of plants and microbes is a growing scientific field. While the chemical production of certain secondary metabolites is possible, industrial-scale microbial production is a green and economically attractive alternative. This is even more true, given the advances in bioengineering that allow us to alter the workings of microbes in order to increase their production of compounds of interest. This type of engineering requires detailed knowledge of the "chassis" organism's metabolism. Since the resources and the catalytic capacity of enzymes in microbes is finite, it is important to examine the tradeoffs between various bioprocesses in an engineered system and alter its working in a manner that minimally perturbs the robustness of the system while allowing for the maximum production of a product of interest. The in silico multi-objective analysis of metabolism using genome-scale models is an ideal method for such examinations.

Integrating proteomic data with metabolic modeling provides insight into key pathways of Bordetella pertussis biofilms

Pertussis, commonly known as whooping cough is a severe respiratory disease caused by the bacterium, Bordetella pertussis. Despite widespread vaccination, pertussis resurgence has been observed globally. The development of the current acellular vaccine (ACV) has been based on planktonic studies. However, recent studies have shown that B. pertussis readily forms biofilms. A better understanding of B. pertussis biofilms is important for developing novel vaccines that can target all aspects of B. pertussis infection. This study compared the proteomic expression of biofilm and planktonic B. pertussis cells to identify key changes between the conditions. Major differences were identified in virulence factors including an upregulation of toxins (adenylate cyclase toxin and dermonecrotic toxin) and downregulation of pertactin and type III secretion system proteins in biofilm cells. To further dissect metabolic pathways that are altered during the biofilm lifestyle, the proteomic data was then incorporated into a genome scale metabolic model using the Integrative Metabolic Analysis Tool (iMAT). The generated models predicted that planktonic cells utilised the glyoxylate shunt while biofilm cells completed the full tricarboxylic acid cycle. Differences in processing aspartate, arginine and alanine were identified as well as unique export of valine out of biofilm cells which may have a role in inter-bacterial communication and regulation. Finally, increased polyhydroxybutyrate accumulation and superoxide dismutase activity in biofilm cells may contribute to increased persistence during infection. Taken together, this study modeled major proteomic and metabolic changes that occur in biofilm cells which helps lay the groundwork for further understanding B. pertussis pathogenesis.