Metabolic Modelling as a Framework for Metabolomics Data Integration and Analysis

Personalized metabolic whole-body models for newborns and infants predict growth and biomarkers of inherited metabolic diseases

Comprehensive whole-body models (WBMs) accounting for organ-specific dynamics have been developed to simulate adult metabolism, but such models do not exist for infants. Here, we present a resource of 360 organ-resolved, sex-specific models of newborn and infant metabolism (infant-WBMs) spanning the first 180 days of life. These infant-WBMs were parameterized to represent the distinct metabolic characteristics of newborns and infants, including nutrition, energy requirements, and thermoregulation. We demonstrate that the predicted infant growth was consistent with the recommendation by the World Health Organization. We assessed the infant-WBMs' reliability and capabilities for personalization by simulating 10,000 newborns based on their blood metabolome and birth weight. Furthermore, the infant-WBMs accurately predicted changes in known biomarkers over time and metabolic responses to treatment strategies for inherited metabolic diseases. The infant-WBM resource holds promise for personalized medicine, as the infant-WBMs could be a first step to digital metabolic twins for newborn and infant metabolism.

Reconstruction of cell-specific models capturing the influence of metabolism on DNA methylation in cancer

The imbalance of epigenetic regulatory mechanisms such as DNA methylation, which can promote aberrant gene expression profiles without affecting the DNA sequence, may cause the deregulation of signaling, regulatory, and metabolic processes, contributing to a cancerous phenotype. Since some metabolites are substrates and cofactors of epigenetic regulators, their availability can be affected by characteristic cancer cell metabolic shifts, feeding cancer onset and progression through epigenetic deregulation. Hence, there is a need to study the influence of cancer metabolic reprogramming in DNA methylation to design new effective treatments. In this study, a generic Genome-Scale Metabolic Model (GSMM) of a human cell, integrating DNA methylation or demethylation reactions, was obtained and used for the reconstruction of Genome-Scale Metabolic Models enhanced with Enzymatic Constraints using Kinetic and Omics data (GECKOs) of 31 cancer cell lines. Furthermore, cell-line-specific DNA methylation levels were included in the models, as coefficients of a DNA composition pseudo-reaction, to depict the influence of metabolism over global DNA methylation in each of the cancer cell lines. Flux simulations demonstrated the ability of these models to provide simulated fluxes of exchange reactions similar to the equivalent experimentally measured uptake/secretion rates and to make good functional predictions. In addition, simulations found metabolic pathways, reactions and enzymes directly or inversely associated with the gene promoter methylation. Two potential candidates for targeted cancer epigenetic therapy were identified.

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.

Monitoring and modelling the glutamine metabolic pathway: a review and future perspectives

BACKGROUND

Analysis of the glutamine metabolic pathway has taken a special place in metabolomics research in recent years, given its important role in cell biosynthesis and bioenergetics across several disorders, especially in cancer cell survival. The science of metabolomics addresses the intricate intracellular metabolic network by exploring and understanding how cells function and respond to external or internal perturbations to identify potential therapeutic targets. However, despite recent advances in metabolomics, monitoring the kinetics of a metabolic pathway in a living cell in situ, real-time and holistically remains a significant challenge.

AIM

This review paper explores the range of analytical approaches for monitoring metabolic pathways, as well as physicochemical modeling techniques, with a focus on glutamine metabolism. We discuss the advantages and disadvantages of each method and explore the potential of label-free Raman microspectroscopy, in conjunction with kinetic modeling, to enable real-time and in situ monitoring of the cellular kinetics of the glutamine metabolic pathway.

KEY SCIENTIFIC CONCEPTS

Given its important role in cell metabolism, the ability to monitor and model the glutamine metabolic pathways are highlighted. Novel, label free approaches have the potential to revolutionise metabolic biosensing, laying the foundation for a new paradigm in metabolomics research and addressing the challenges in monitoring metabolic pathways in living cells.

Towards advanced bioprocess optimization: A multiscale modelling approach

Mammalian cells produce up to 80 % of the commercially available therapeutic proteins, with Chinese Hamster Ovary (CHO) cells being the primary production host. Manufacturing involves a train of reactors, the last of which is typically run in fed-batch mode, where cells grow and produce the required protein. The feeding strategy is decided a priori, from either past operations or the design of experiments and rarely considers the current state of the process. This work proposes a Model Predictive Control (MPC) formulation based on a hybrid kinetic-stoichiometric reactor model to provide optimal feeding policies in real-time, which is agnostic to the culture, hence transferable across CHO cell culture systems. The benefits of the proposed controller formulation are demonstrated through a comparison between an open-loop simulation and closed-loop optimization, using a digital twin as an emulator of the process.

Recent advances in mass spectrometry-based computational metabolomics

The computational metabolomics field brings together computer scientists, bioinformaticians, chemists, clinicians, and biologists to maximize the impact of metabolomics across a wide array of scientific and medical disciplines. The field continues to expand as modern instrumentation produces datasets with increasing complexity, resolution, and sensitivity. These datasets must be processed, annotated, modeled, and interpreted to enable biological insight. Techniques for visualization, integration (within or between omics), and interpretation of metabolomics data have evolved along with innovation in the databases and knowledge resources required to aid understanding. In this review, we highlight recent advances in the field and reflect on opportunities and innovations in response to the most pressing challenges. This review was compiled from discussions from the 2022 Dagstuhl seminar entitled "Computational Metabolomics: From Spectra to Knowledge".

Phenotype-specific estimation of metabolic fluxes using gene expression data

A cell's genome influences its metabolism via the expression of enzyme-related genes, but transcriptome and fluxome are not perfectly correlated as post-transcriptional mechanisms also regulate reaction's kinetics. Here, we addressed the question: given a transcriptome, how unobserved mechanisms of reaction kinetics should be systematically accounted for when inferring the fluxome? To infer the most likely and least biased fluxome, we present Pheflux, a constraint-based model maximizing Shannon's entropy of fluxes per mRNA. Benchmarked against ¹³C fluxes of yeast and bacteria, Pheflux accurately estimates the carbon core metabolism. We applied Pheflux to thousands of normal and tumor cell transcriptomes obtained from The Cancer Genome Atlas. Pheflux showed statistically significantly higher glucose yields on lactate in breast, kidney, and bronchus-lung tumoral cells than their normal counterparts. Results are consistent with the Warburg effect, a hallmark of cancer metabolism, suggesting that Pheflux can be efficiently used to study the metabolism of eukaryotic cells.

Towards inferring absolute concentrations from relative abundance in time-course GC-MS metabolomics data

Metabolomics, the large-scale study of metabolites, has significant appeal as a source of information for metabolic modeling and other scientific applications. One common approach for measuring metabolomics data is gas chromatography-mass spectrometry (GC-MS). However, GC-MS metabolomics data are typically reported as relative abundances, precluding their use with approaches and tools where absolute concentrations are necessary. While chemical standards can be used to help provide quantification, their use is time-consuming, expensive, or even impossible due to their limited availability. The ability to infer absolute concentrations from GC-MS metabolomics data without chemical standards would have significant value. We hypothesized that when analyzing time-course metabolomics datasets, the mass balances of metabolism and other biological information could provide sufficient information towards inference of absolute concentrations. To demonstrate this, we developed and characterized MetaboPAC, a computational framework that uses two approaches-one based on kinetic equations and another using biological heuristics-to predict the most likely response factors that allow translation between relative abundances and absolute concentrations. When used to analyze noiseless synthetic data generated from multiple types of kinetic rate laws, MetaboPAC performs significantly better than negative control approaches when 20% of kinetic terms are known a priori. Under conditions of lower sampling frequency and high noise, MetaboPAC is still able to provide significant inference of concentrations in 3 of 4 models studied. This provides a starting point for leveraging biological knowledge to extract concentration information from time-course intracellular GC-MS metabolomics datasets, particularly for systems that are well-studied and have partially known kinetic structures.

Specifics of Metabolite-Protein Interactions and Their Computational Analysis and Prediction

Computational approaches to the characterization and prediction of compound-protein interactions have a long research history and are well established, driven primarily by the needs of drug development. While, in principle, many of the computational methods developed in the context of drug development can also be applied directly to the investigation of metabolite-protein interactions, the interactions of metabolites with proteins (enzymes) are characterized by a number of particularities that result from their natural evolutionary origin and their biological and biochemical roles, as well as from a different problem setting when investigating them. In this review, these special aspects will be highlighted and recent research on them and developed computational approaches presented, along with available resources. They concern, among others, binding promiscuity, allostery, the role of posttranslational modifications, molecular steering and crowding effects, and metabolic conversion rate predictions. Recent breakthroughs in the field of protein structure prediction and newly developed machine learning techniques are being discussed as a tremendous opportunity for developing a more detailed molecular understanding of metabolism.

Metabolomics during the spontaneous fermentation in cocoa (Theobroma cacao L.): An exploraty review

Spontaneous fermentation is a process that depends on substrates' physical characteristics, crop variety, and postharvest practices; it induces variations in the metabolites that are responsible for the taste, aroma, and quality. Metabolomics makes it possible to detect key metabolites using chemometrics and makes it possible to establish patterns or identify biomarker behaviors under certain conditions at a given time. Therefore, sensitive and highly efficient analytical techniques allow for studying the metabolomic fingerprint changes during fermentation; which identify and quantify metabolites related to taste and aroma formation of an adequate processing time. This review shows that studying metabolomics in spontaneous fermentation permits the characterization of spontaneous fermentation in different stages. Also, it demonstrates the possibility of modulating the quality of cocoa by improving the spontaneous fermentation time (because of volatile aromatic compounds formation), thus standardizing the process to obtain attributes and quality that will later impact the chocolate quality.

The metabolomics of human aging: Advances, challenges, and opportunities

As the global population becomes older, understanding the impact of aging on health and disease becomes paramount. Recent advancements in multiomic technology have allowed for the high-throughput molecular characterization of aging at the population level. Metabolomics studies that analyze the small molecules in the body can provide biological information across a diversity of aging processes. Here, we review the growing body of population-scale metabolomics research on aging in humans, identifying the major trends in the field, implicated biological pathways, and how these pathways relate to health and aging. We conclude by assessing the main challenges in the research to date, opportunities for advancing the field, and the outlook for precision health applications.

Kinetic and data-driven modeling of pancreatic β-cell central carbon metabolism and insulin secretion

Pancreatic β-cells respond to increased extracellular glucose levels by initiating a metabolic shift. That change in metabolism is part of the process of glucose-stimulated insulin secretion and is of particular interest in the context of diabetes. However, we do not fully understand how the coordinated changes in metabolic pathways and metabolite products influence insulin secretion. In this work, we apply systems biology approaches to develop a detailed kinetic model of the intracellular central carbon metabolic pathways in pancreatic β-cells upon stimulation with high levels of glucose. The model is calibrated to published metabolomics datasets for the INS1 823/13 cell line, accurately capturing the measured metabolite fold-changes. We first employed the calibrated mechanistic model to estimate the stimulated cell's fluxome. We then used the predicted network fluxes in a data-driven approach to build a partial least squares regression model. By developing the combined kinetic and data-driven modeling framework, we gain insights into the link between β-cell metabolism and glucose-stimulated insulin secretion. The combined modeling framework was used to predict the effects of common anti-diabetic pharmacological interventions on metabolite levels, flux through the metabolic network, and insulin secretion. Our simulations reveal targets that can be modulated to enhance insulin secretion. The model is a promising tool to contextualize and extend the usefulness of metabolomics data and to predict dynamics and metabolite levels that are difficult to measure in vitro. In addition, the modeling framework can be applied to identify, explain, and assess novel and clinically-relevant interventions that may be particularly valuable in diabetes treatment.

Metabolomics and modelling approaches for systems metabolic engineering

Metabolic engineering involves the manipulation of microbes to produce desirable compounds through genetic engineering or synthetic biology approaches. Metabolomics involves the quantitation of intracellular and extracellular metabolites, where mass spectrometry and nuclear magnetic resonance based analytical instrumentation are often used. Here, the experimental designs, sample preparations, metabolite quenching and extraction are essential to the quantitative metabolomics workflow. The resultant metabolomics data can then be used with computational modelling approaches, such as kinetic and constraint-based modelling, to better understand underlying mechanisms and bottlenecks in the synthesis of desired compounds, thereby accelerating research through systems metabolic engineering. Constraint-based models, such as genome scale models, have been used successfully to enhance the yield of desired compounds from engineered microbes, however, unlike kinetic or dynamic models, constraint-based models do not incorporate regulatory effects. Nevertheless, the lack of time-series metabolomic data generation has hindered the usefulness of dynamic models till today. In this review, we show that improvements in automation, dynamic real-time analysis and high throughput workflows can drive the generation of more quality data for dynamic models through time-series metabolomics data generation. Spatial metabolomics also has the potential to be used as a complementary approach to conventional metabolomics, as it provides information on the localization of metabolites. However, more effort must be undertaken to identify metabolites from spatial metabolomics data derived through imaging mass spectrometry, where machine learning approaches could prove useful. On the other hand, single-cell metabolomics has also seen rapid growth, where understanding cell-cell heterogeneity can provide more insights into efficient metabolic engineering of microbes. Moving forward, with potential improvements in automation, dynamic real-time analysis, high throughput workflows, and spatial metabolomics, more data can be produced and studied using machine learning algorithms, in conjunction with dynamic models, to generate qualitative and quantitative predictions to advance metabolic engineering efforts.

Metabolic flux between organs measured by arteriovenous metabolite gradients

Mammalian organs convert dietary nutrients into circulating metabolites and share them to maintain whole-body metabolic homeostasis. While the concentrations of circulating metabolites have been frequently measured in a variety of pathophysiological conditions, the exchange flux of circulating metabolites between organs is not easily measurable due to technical difficulties. Isotope tracing is useful for measuring such fluxes for a metabolite of interest, but the shuffling of isotopic atoms between metabolites requires mathematical modeling. Arteriovenous metabolite gradient measurements can complement isotope tracing to infer organ-specific net fluxes of many metabolites simultaneously. Here, we review the historical development of arteriovenous measurements and discuss their advantages and limitations with key example studies that have revealed metabolite exchange flux between organs in diverse pathophysiological contexts.

Staphylococcus epidermidis RP62A’s Metabolic Network: Validation and Intervention Strategies

Increasingly, systems biology is gaining relevance in basic and applied research. The combination of computational biology with wet laboratory methods produces synergy that results in an exponential increase in knowledge of biological systems. The study of microorganisms such as Staphylococcus epidermidis RP62A enables the researcher to understand better their metabolic networks, which allows the design of effective strategies to treat infections caused by this species or others. S. epidermidis is the second most commoncause of infection in patients with joint implants, so treating its proliferation seems vital for public health. There are different approaches to the analysis of metabolic networks. Flux balance analysis (FBA) is one of the most widespread streams of research. It allows the study of large metabolic networks, the study their structural properties, the optimization of metabolic flux, and the search for intervention strategies to modify the state of the metabolic network. This work presents the validation of the Staphylococcus epidermidis RP62A metabolic network model elaborated by Díaz Calvo et al. Then, we elaborate further on the network analysis’s essential reactions. The full set of essential reactions (including a previously unobserved one) was computed, and we classified them into equivalence classes. Some proposals to intervene in the network and design knock-outs by studying minimal cut sets of small length are also introduced. In particular, minimal cut sets related to the medium (including exchange reactions associated with medium metabolites) have been computed. In this sense, the unique external MCS (composed of cysteine and sulfate ion) has been found, and all hybrid MCS (based on knocking out both internal and exchange reactions) of length two have also been computed. The paper also points out the possible importance of these new intervention strategies.

Metabolic systems analysis identifies a novel mechanism contributing to shock in patients with endotheliopathy of trauma (EoT) involving thromboxane A2 and LTC4

Purpose

Endotheliopathy of trauma (EoT), as defined by circulating levels of syndecan-1 ≥ 40 ng/mL, has been reported to be associated with significantly increased transfusion requirements and a doubled 30-day mortality. Increased shedding of the glycocalyx points toward the endothelial cell membrane composition as important for the clinical outcome being the rationale for this study.

Results

The plasma metabolome of 95 severely injured trauma patients was investigated by mass spectrometry, and patients with EoT vs. non-EoT were compared by partial least square-discriminant analysis, identifying succinic acid as the top metabolite to differentiate EoT and non-EoT patients (VIP score = 3). EoT and non-EoT patients' metabolic flux profile was inferred by integrating the corresponding plasma metabolome data into a genome-scale metabolic network reconstruction analysis and performing a functional study of the metabolic capabilities of each group. Model predictions showed a decrease in cholesterol metabolism secondary to impaired mevalonate synthesis in EoT compared to non-EoT patients. Intracellular task analysis indicated decreased synthesis of thromboxanA2 and leukotrienes, as well as a lower carnitine palmitoyltransferase I activity in EoT compared to non-EoT patients. Sensitivity analysis also showed a significantly high dependence of eicosanoid-associated metabolic tasks on alpha-linolenic acid as unique to EoT patients.

Conclusions

Model-driven analysis of the endothelial cells' metabolism identified potential novel targets as impaired thromboxane A2 and leukotriene synthesis in EoT patients when compared to non-EoT patients. Reduced thromboxane A2 and leukotriene availability in the microvasculature impairs vasoconstriction ability and may thus contribute to shock in EoT patients. These findings are supported by extensive scientific literature; however, further investigations are required on these findings.

Targeted metabolomics suggests a probable role of the FTO gene in the kynurenine pathway in prediabetes

Background

Genome-wide association studies have identified the alpha-ketoglutarate dependent dioxygenase gene (FTO) as the first susceptibility gene of obesity. In the present study, we utilized targeted metabolomics in an attempt to further elucidate mechanisms underlying the action of the FTO gene.

Methods

This study was part of a health survey of employees of the Electricity Generating Authority of Thailand (n = 79, 10 female and 69 male). Targeted metabolomics was performed by liquid chromatography-mass spectrometry using Biocrates AbsoluteIDQ-p180 kit. Genotyping of FTO rs9939609 was performed by real-time PCR (TaqMan™ MGB probes).

Results

Using OPLS-DA variable importance in projection (VIP), tryptophan was found to be among the metabolites with the 10 highest VIP scores. Pearson's correlation analysis showed that kynurenine and tryptophan were positively correlated only in subjects with the rs9939609 A allele (n = 32, r = 0.56, p < 0.001) and the correlation coefficients were significantly higher in subjects having the A allele than in those without the A allele (p < 0.05). Moreover, the kynurenine/tryptophan ratio was significantly associated with the presence of the A allele, independently of body mass index and sex.

Conclusions

The FTO gene is likely to influences the conversion of tryptophan to kynurenine.

Multi-omics approaches and genetic engineering of metabolism for improved biorefinery and wastewater treatment in microalgae

Microalgae, a group of photosynthetic microorganisms rich in diverse and novel bioactive metabolites, have been explored for the production of biofuels, high value-added compounds as food and feeds, and pharmaceutical chemicals as agents with therapeutic benefits. This article reviews the development of omics resources and genetic engineering techniques including gene transformation methodologies, mutagenesis, and genome-editing tools in microalgae biorefinery and wastewater treatment. The introduction of these enlisted techniques has simplified the understanding of complex metabolic pathways undergoing microalgal cells. The multiomics approach of the integrated omics datasets, big data analysis, and machine learning for the discovery of objective traits and genes responsible for metabolic pathways was reviewed. Recent advances and limitations of multiomics analysis and genetic bioengineering technology to facilitate the improvement of microalgae as the dual role of wastewater treatment and biorefinery feedstock production are discussed. This article is protected by copyright. All rights reserved.

Burden Imposed by Heterologous Protein Production in Two Major Industrial Yeast Cell Factories: Identifying Sources and Mitigation Strategies

Production of heterologous proteins, especially biopharmaceuticals and industrial enzymes, in living cell factories consumes cellular resources. Such resources are reallocated from normal cellular processes toward production of the heterologous protein that is often of no benefit to the host cell. This competition for resources is a burden to host cells, has a negative impact on cell fitness, and may consequently trigger stress responses. Importantly, this often causes a reduction in final protein titers. Engineering strategies to generate more burden resilient production strains offer sustainable opportunities to increase production and profitability for this growing billion-dollar global industry. We review recently reported impacts of burden derived from resource competition in two commonly used protein-producing yeast cell factories: Saccharomyces cerevisiae and Komagataella phaffii (syn. Pichia pastoris). We dissect possible sources of burden in these organisms, from aspects related to genetic engineering to protein translation and export of soluble protein. We also summarize advances as well as challenges for cell factory design to mitigate burden and increase overall heterologous protein production from metabolic engineering, systems biology, and synthetic biology perspectives. Lastly, future profiling and engineering strategies are highlighted that may lead to constructing robust burden-resistant cell factories. This includes incorporation of systems-level data into mathematical models for rational design and engineering dynamical regulation circuits in production strains.

Metabolic systems biology and multi-omics of cyanobacteria: Perspectives and future directions

Cyanobacteria are oxygenic photoautotrophs whose metabolism contains key biochemical pathways to fix atmospheric CO₂ and synthesize various metabolites. The development of bioengineering tools has enabled the manipulation of cyanobacterial chassis to produce various valuable bioproducts photosynthetically. However, effective utilization of cyanobacteria as photosynthetic cell factories needs a detailed understanding of their metabolism and its interaction with other cellular processes. Implementing systems and synthetic biology tools has generated a wealth of information on various metabolic pathways. However, to design effective engineering strategies for further improvement in growth, photosynthetic efficiency, and enhanced production of target biochemicals, in-depth knowledge of their carbon/nitrogen metabolism, pathway fluxe distribution, genetic regulation and integrative analyses are necessary. In this review, we discuss the recent advances in the development of genome-scale metabolic models (GSMMs), omics analyses (metabolomics, transcriptomics, proteomics, fluxomics), and integrative modeling approaches to showcase the current understanding of cyanobacterial metabolism.

Protocol for hybrid flux balance, statistical, and machine learning analysis of multi-omic data from the cyanobacterium Synechococcus sp. PCC 7002

Combining a computational framework for flux balance analysis with machine learning improves the accuracy of predicting metabolic activity across conditions, while enabling mechanistic interpretation. This protocol presents a guide to condition-specific metabolic modeling that integrates regularized flux balance analysis with machine learning approaches to extract key features from transcriptomic and fluxomic data. We demonstrate the protocol as applied to Synechococcus sp. PCC 7002; we also outline how it can be adapted to any species or community with available multi-omic data. For complete details on the use and execution of this protocol, please refer to Vijayakumar et al. (2020).

Toward modeling metabolic state from single-cell transcriptomics

Background

Single-cell metabolic studies bring new insights into cellular function, which can often not be captured on other omics layers. Metabolic information has wide applicability, such as for the study of cellular heterogeneity or for the understanding of drug mechanisms and biomarker development. However, metabolic measurements on single-cell level are limited by insufficient scalability and sensitivity, as well as resource intensiveness, and are currently not possible in parallel with measuring transcript state, commonly used to identify cell types. Nevertheless, because omics layers are strongly intertwined, it is possible to make metabolic predictions based on measured data of more easily measurable omics layers together with prior metabolic network knowledge.

Scope of Review

We summarize the current state of single-cell metabolic measurement and modeling approaches, motivating the use of computational techniques. We review three main classes of computational methods used for prediction of single-cell metabolism: pathway-level analysis, constraint-based modeling, and kinetic modeling. We describe the unique challenges arising when transitioning from bulk to single-cell modeling. Finally, we propose potential model extensions and computational methods that could be leveraged to achieve these goals.

Major Conclusions

Single-cell metabolic modeling is a rising field that provides a new perspective for understanding cellular functions. The presented modeling approaches vary in terms of input requirements and assumptions, scalability, modeled metabolic layers, and newly gained insights. We believe that the use of prior metabolic knowledge will lead to more robust predictions and will pave the way for mechanistic and interpretable machine-learning models.

•

Single-cell RNA sequencing and prior metabolic knowledge enable metabolic predictions.

•

When compared to bulk, single-cell modeling is linked to unique insights and challenges.

•

Computational modelling approaches differ in applicability and newly provided insights.

•

The use of prior metabolic knowledge paves the way for mechanistic machine-learning.

Topological Analysis of the Carbon-Concentrating CETCH Cycle and a Photorespiratory Bypass Reveals Boosted CO2-Sequestration by Plants

Synthetically designed alternative photorespiratory pathways increase the biomass of tobacco and rice plants. Likewise, some in planta-tested synthetic carbon-concentrating cycles (CCCs) hold promise to increase plant biomass while diminishing atmospheric carbon dioxide burden. Taking these individual contributions into account, we hypothesize that the integration of bypasses and CCCs will further increase plant productivity. To test this in silico, we reconstructed a metabolic model by integrating photorespiration and photosynthesis with the synthetically designed alternative pathway 3 (AP3) enzymes and transporters. We calculated fluxes of the native plant system and those of AP3 combined with the inhibition of the glycolate/glycerate transporter by using the YANAsquare package. The activity values corresponding to each enzyme in photosynthesis, photorespiration, and for synthetically designed alternative pathways were estimated. Next, we modeled the effect of the crotonyl-CoA/ethylmalonyl-CoA/hydroxybutyryl-CoA cycle (CETCH), which is a set of natural and synthetically designed enzymes that fix CO₂ manifold more than the native Calvin-Benson-Bassham (CBB) cycle. We compared estimated fluxes across various pathways in the native model and under an introduced CETCH cycle. Moreover, we combined CETCH and AP3-w/plgg1RNAi, and calculated the fluxes. We anticipate higher carbon dioxide-harvesting potential in plants with an AP3 bypass and CETCH-AP3 combination. We discuss the in vivo implementation of these strategies for the improvement of C3 plants and in natural high carbon harvesters.

Genome Scale Modeling to Study the Metabolic Competition between Cells in the Tumor Microenvironment

The tumor's physiology emerges from the dynamic interplay of numerous cell types, such as cancer cells, immune cells and stromal cells, within the tumor microenvironment. Immune and cancer cells compete for nutrients within the tumor microenvironment, leading to a metabolic battle between these cell populations. Tumor cells can reprogram their metabolism to meet the high demand of building blocks and ATP for proliferation, and to gain an advantage over the action of immune cells. The study of the metabolic reprogramming mechanisms underlying cancer requires the quantification of metabolic fluxes which can be estimated at the genome-scale with constraint-based or kinetic modeling. Constraint-based models use a set of linear constraints to simulate steady-state metabolic fluxes, whereas kinetic models can simulate both the transient behavior and steady-state values of cellular fluxes and concentrations. The integration of cell- or tissue-specific data enables the construction of context-specific models that reflect cell-type- or tissue-specific metabolic properties. While the available modeling frameworks enable limited modeling of the metabolic crosstalk between tumor and immune cells in the tumor stroma, future developments will likely involve new hybrid kinetic/stoichiometric formulations.

Thermodynamic Genome-Scale Metabolic Modeling of Metallodrug Resistance in Colorectal Cancer

Simple Summary

Cancer, but also its treatment, can lead to a reprogramming of cellular metabolism. These changes are observable in metabolite abundances, which can be unbiasedly measured via mass spectrometry metabolomics. However, even when the metabolome changes strongly, a (mechanistic) interpretation is difficult as metabolite levels do not necessarily directly correspond to pathway activities. Here we measure the changes of the cellular metabolome in colorectal cancer cell lines sensitive and resistant to the ruthenium-based drug BOLD-100/KP1339 and the platinum-based drug oxaliplatin. We map these changes onto a cancer-specific genome-scale metabolic model, which allows us not only to compute intracellular flux distributions, but also to disentangle drug-specific effects from growth differences from differences in metabolic adaptations due to resistance. Specifically, we find that resistance to BOLD-100/KP1339 induces more extensive reprogramming than oxaliplatin, especially with respect to fatty acid and amino acid metabolism.

Abstract

Background: Mass spectrometry-based metabolomics approaches provide an immense opportunity to enhance our understanding of the mechanisms that underpin the cellular reprogramming of cancers. Accurate comparative metabolic profiling of heterogeneous conditions, however, is still a challenge. Methods: Measuring both intracellular and extracellular metabolite concentrations, we constrain four instances of a thermodynamic genome-scale metabolic model of the HCT116 colorectal carcinoma cell line to compare the metabolic flux profiles of cells that are either sensitive or resistant to ruthenium- or platinum-based treatments with BOLD-100/KP1339 and oxaliplatin, respectively. Results: Normalizing according to growth rate and normalizing resistant cells according to their respective sensitive controls, we are able to dissect metabolic responses specific to the drug and to the resistance states. We find the normalization steps to be crucial in the interpretation of the metabolomics data and show that the metabolic reprogramming in resistant cells is limited to a select number of pathways. Conclusions: Here, we elucidate the key importance of normalization steps in the interpretation of metabolomics data, allowing us to uncover drug-specific metabolic reprogramming during acquired metal-drug resistance.

Model-guided dynamic control of essential metabolic nodes boosts acetyl-coenzyme A-dependent bioproduction in rewired Pseudomonas putida

Pseudomonas putida is evolutionarily endowed with features relevant for bioproduction, especially under harsh operating conditions. The rich metabolic versatility of this species, however, comes at the price of limited formation of acetyl-coenzyme A (CoA) from sugar substrates. Since acetyl-CoA is a key metabolic precursor for a number of added-value products, in this work we deployed an in silico-guided rewiring program of central carbon metabolism for upgrading P. putida as a host for acetyl-CoA-dependent bioproduction. An updated kinetic model, integrating fluxomics and metabolomics datasets in addition to manually-curated information of enzyme mechanisms, identified targets that would lead to increased acetyl-CoA levels. Based on these predictions, a set of plasmids based on clustered regularly interspaced short palindromic repeats (CRISPR) and dead CRISPR-associated protein 9 (dCas9) was constructed to silence genes by CRISPR interference (CRISPRi). Dynamic reduction of gene expression of two key targets (gltA, encoding citrate synthase, and the essential accA gene, encoding subunit A of the acetyl-CoA carboxylase complex) mediated an 8-fold increase in the acetyl-CoA content of rewired P. putida. Poly(3-hydroxybutyrate) (PHB) was adopted as a proxy of acetyl-CoA availability, and two synthetic pathways were engineered for biopolymer accumulation. By including cell morphology as an extra target for the CRISPRi approach, fully rewired P. putida strains programmed for PHB accumulation had a 5-fold increase in PHB titers in bioreactor cultures using glucose. Thus, the strategy described herein allowed for rationally redirecting metabolic fluxes in P. putida from central metabolism towards product biosynthesis-especially relevant when deletion of essential pathways is not an option.

Integration of transcriptional and metabolic control in macrophage activation

Macrophages react to microbial and endogenous danger signals by activating a broad panel of effector and homeostatic responses. Such responses entail rapid and stimulus-specific changes in gene expression programs accompanied by extensive rewiring of metabolism, with alterations in chromatin modifications providing one layer of integration of transcriptional and metabolic regulation. A systematic and mechanistic understanding of the mutual influences between signal-induced metabolic changes and gene expression is still lacking. Here, we discuss current evidence, controversies, knowledge gaps, and future areas of investigation on how metabolic and transcriptional changes are dynamically integrated during macrophage activation. The cross-talk between metabolism and inflammatory gene expression is in part accounted for by alterations in the production, usage, and availability of metabolic intermediates that impact the macrophage epigenome. In addition, stimulus-inducible gene expression changes alter the production of inflammatory mediators, such as nitric oxide, that in turn modulate the activity of metabolic enzymes thus determining complex regulatory loops. Critical issues remain to be understood, notably whether and how metabolic rewiring can bring about gene-specific (as opposed to global) expression changes.

Metabolomics-Guided Elucidation of Plant Abiotic Stress Responses in the 4IR Era: An Overview

Plants are constantly challenged by changing environmental conditions that include abiotic stresses. These are limiting their development and productivity and are subsequently threatening our food security, especially when considering the pressure of the increasing global population. Thus, there is an urgent need for the next generation of crops with high productivity and resilience to climate change. The dawn of a new era characterized by the emergence of fourth industrial revolution (4IR) technologies has redefined the ideological boundaries of research and applications in plant sciences. Recent technological advances and machine learning (ML)-based computational tools and omics data analysis approaches are allowing scientists to derive comprehensive metabolic descriptions and models for the target plant species under specific conditions. Such accurate metabolic descriptions are imperatively essential for devising a roadmap for the next generation of crops that are resilient to environmental deterioration. By synthesizing the recent literature and collating data on metabolomics studies on plant responses to abiotic stresses, in the context of the 4IR era, we point out the opportunities and challenges offered by omics science, analytical intelligence, computational tools and big data analytics. Specifically, we highlight technological advancements in (plant) metabolomics workflows and the use of machine learning and computational tools to decipher the dynamics in the chemical space that define plant responses to abiotic stress conditions.

Guided extraction of genome-scale metabolic models for the integration and analysis of omics data

Omics data can be integrated into a reference model using various model extraction methods (MEMs) to yield context-specific genome-scale metabolic models (GEMs). How to chose the appropriate MEM, thresholding rule and threshold remains a challenge. We integrated mouse transcriptomic data from a Cyp51 knockout mice diet experiment (GSE58271) using five MEMs (GIMME, iMAT, FASTCORE, INIT an tINIT) in a combination with a recently published mouse GEM iMM1865. Except for INIT and tINIT, the size of extracted models varied with the MEM used (t-test: p-value < 0.001). The Jaccard index of iMAT models ranged from 0.27 to 1.0. Out of the three factors under study in the experiment (diet, gender and genotype), gender explained most of the variability ( > 90%) in PC1 for FASTCORE. In iMAT, each of the three factors explained less than 40% of the variability within PC1, PC2 and PC3. Among all the MEMs, FASTCORE captured the most of the true variability in the data by clustering samples by gender. Our results show that for the efficient use of MEMs in the context of omics data integration and analysis, one should apply various MEMs, thresholding rules, and thresholding values to select the MEM and its configuration that best captures the true variability in the data. This selection can be guided by the methodology as proposed and used in this paper. Moreover, we describe certain approaches that can be used to analyse the results obtained with the selected MEM and to put these results in a biological context.

Challenges and Perspectives of Polyhydroxyalkanoate Production From Microalgae/Cyanobacteria and Bacteria as Microbial Factories: An Assessment of Hybrid Biological System

Polyhydroxyalkanoates (PHAs) are the biopolymer of choice if we look for a substitute of petroleum-based non-biodegradable plastics. Microbial production of PHAs as carbon reserves has been studied for decades and PHAs are gaining attention for a wide range of applications in various fields. Still, their uneconomical production is the major concern largely attributed to high cost of organic substrates for PHA producing heterotrophic bacteria. Therefore, microalgae/cyanobacteria, being photoautotrophic, prove to have an edge over heterotrophic bacteria. They have minimal metabolic requirements, such as inorganic nutrients (CO2, N, P, etc.) and light, and they can survive under adverse environmental conditions. PHA production under photoautotrophic conditions has been reported from cyanobacteria, the only candidate among prokaryotes, and few of the eukaryotic microalgae. However, an efficient cultivation system is still required for photoautotrophic PHA production to overcome the limitations associated with (1) stringent management of closed photobioreactors and (2) optimization of monoculture in open pond culture. Thus, a hybrid system is a necessity, involving the participation of microalgae/cyanobacteria and bacteria, i.e., both photoautotrophic and heterotrophic components having mutual interactive benefits for each other under different cultivation regime, e.g., mixotrophic, successive two modules, consortium based, etc. Along with this, further strategies like optimization of culture conditions (N, P, light exposure, CO2 dynamics, etc.), bioengineering, efficient downstream processes, and the application of mathematical/network modeling of metabolic pathways to improve PHA production are the key areas discussed here. Conclusively, this review aims to critically analyze cyanobacteria as PHA producers and proposes economically sustainable production of PHA from microbial autotrophs and heterotrophs in "hybrid biological system."

Systematic Review of Multi-Omics Approaches to Investigate Toxicological Effects in Macrophages

Insights into the modes of action (MoAs) of xenobiotics are of utmost importance for the definition of adverse outcome pathways (AOPs), which are essential for a mechanism-based risk assessment. A well-established strategy to reveal MoAs of xenobiotics is the use of omics. However, often an even more comprehensive approach is needed, which can be achieved using multi-omics. Since the immune system plays a central role in the defense against foreign substances and pathogens, with the innate immune system building a first barrier, we systematically reviewed multi-omics studies investigating the effects of xenobiotics on macrophages. Surprisingly, only nine publications were identified, combining proteomics with transcriptomics or metabolomics. We summarized pathways and single proteins, transcripts, or metabolites, which were described to be affected upon treatment with xenobiotics in the reviewed studies, thus revealing a broad range of effects. In summary, we show that macrophages are a relevant model system to investigate the toxicological effects induced by xenobiotics. Furthermore, the multi-omics approaches led to a more comprehensive overview compared to only one omics layer with slight advantages for combinations that complement each other directly, e.g., proteome and metabolome.

Seeing the forest for the trees: Retrieving plant secondary biochemical pathways from metabolome networks

Over the last decade, a giant leap forward has been made in resolving the main bottleneck in metabolomics, i.e., the structural characterization of the many unknowns. This has led to the next challenge in this research field: retrieving biochemical pathway information from the various types of networks that can be constructed from metabolome data. Searching putative biochemical pathways, referred to as biotransformation paths, is complicated because several flaws occur during the construction of metabolome networks. Multiple network analysis tools have been developed to deal with these flaws, while in silico retrosynthesis is appearing as an alternative approach. In this review, the different types of metabolome networks, their flaws, and the various tools to trace these biotransformation paths are discussed.