How to measure metabolic fluxes: a taxonomic guide for (13)C fluxomics

Imbalance in Glucose Metabolism Regulates the Transition of Microglia from Homeostasis to Disease-Associated Microglia Stage 1

Microglia undergo two-stage activation in neurodegenerative diseases, known as disease-associated microglia (DAM). TREM2 mediates the DAM2 stage transition, but what regulates the first DAM1 stage transition is unknown. We report that glucose dyshomeostasis inhibits DAM1 activation and PKM2 plays a role. As in tumors, PKM2 was aberrantly elevated in both male and female human AD brains, but unlike in tumors, it is expressed as active tetramers, as well as among TREM2+ microglia surrounding plaques in 5XFAD male and female mice. snRNAseq analyses of microglia without Pkm2 in 5XFAD mice revealed significant increases in DAM1 markers in a distinct metabolic cluster, which is enriched in genes for glucose metabolism, DAM1, and AD risk. 5XFAD mice incidentally exhibited a significant reduction in amyloid pathology without microglial Pkm2 Surprisingly, microglia in 5XFAD without Pkm2 exhibited increases in glycolysis and spare respiratory capacity, which correlated with restoration of mitochondrial cristae alterations. In addition, in situ spatial metabolomics of plaque-bearing microglia revealed an increase in respiratory activity. These results together suggest that it is not only glycolytic but also respiratory inputs that are critical to the development of DAM signatures in 5XFAD mice.

Rethinking 13C-metabolic flux analysis - The Bayesian way of flux inference

Metabolic reaction rates (fluxes) play a crucial role in comprehending cellular phenotypes and are essential in areas such as metabolic engineering, biotechnology, and biomedical research. The state-of-the-art technique for estimating fluxes is metabolic flux analysis using isotopic labelling (13C-MFA), which uses a dataset-model combination to determine the fluxes. Bayesian statistical methods are gaining popularity in the field of life sciences, but the use of 13C-MFA is still dominated by conventional best-fit approaches. The slow take-up of Bayesian approaches is, at least partly, due to the unfamiliarity of Bayesian methods to metabolic engineering researchers. To address this unfamiliarity, we here outline similarities and differences between the two approaches and highlight particular advantages of the Bayesian way of flux analysis. With a real-life example, re-analysing a moderately informative labelling dataset of E. coli, we identify situations in which Bayesian methods are advantageous and more informative, pointing to potential pitfalls of current 13C-MFA evaluation approaches. We propose the use of Bayesian model averaging (BMA) for flux inference as a means of overcoming the problem of model uncertainty through its tendency to assign low probabilities to both, models that are unsupported by data, and models that are overly complex. In this capacity, BMA resembles a tempered Ockham's razor. With the tempered razor as a guide, BMA-based 13C-MFA alleviates the problem of model selection uncertainty and is thereby capable of becoming a game changer for metabolic engineering by uncovering new insights and inspiring novel approaches.

13CO2-labelling and Sampling in Algae for Flux Analysis of Photosynthetic and Central Carbon Metabolism

The flux in photosynthesis can be studied by performing 13CO2 pulse labelling and analysing the temporal labelling kinetics of metabolic intermediates using gas or liquid chromatography linked to mass spectrometry. Metabolic flux analysis (MFA) is the primary approach for analysing metabolic network function and quantifying intracellular metabolic fluxes. Different MFA approaches differ based on the metabolic state (steady vs. non-steady state) and the use of stable isotope tracers. The main methodology used to investigate metabolic systems is metabolite steady state associated with stable isotope labelling experiments. Specifically, in biological systems like photoautotrophic organisms, isotopic non-stationary 113C metabolic flux analysis at metabolic steady state with transient isotopic labelling (13C-INST-MFA) is required. The common requirement for metabolic steady state, alongside its very short half-timed reactions, complicates robust MFA of photosynthetic metabolism. While custom gas chambers design has addressed these challenges in various model plants, no similar tools were developed for liquid photosynthetic cultures (e.g., algae, cyanobacteria), where diffusion and equilibration of inorganic carbon species in the medium entails a new dimension of complexity. Recently, a novel tailor-made microfluidics labelling system has been introduced, supplying short 13CO2 pulses at steady state, and resolving fluxes across most photosynthetic metabolic pathways in algae. The system involves injecting algal cultures and medium containing pre-equilibrated inorganic 13C into a microfluidic mixer, followed by rapid metabolic quenching, enabling precise seconds-level label pulses. This was complemented by a 13CO2-bubbling-based open labelling system (photobioreactor), allowing long pulses (minutes-hours) required for investigating fluxes into central C metabolism and major products. This combined labelling procedure provides a comprehensive fluxome cover for most algal photosynthetic and central C metabolism pathways, thus allowing comparative flux analyses across algae and plants.

Time-resolved, deuterium-based fluxomics uncovers the hierarchy and dynamics of sugar processing by Pseudomonas putida

Pseudomonas putida, a microbial host widely adopted for metabolic engineering, processes glucose through convergent peripheral pathways that ultimately yield 6-phosphogluconate. The periplasmic gluconate shunt (PGS), composed by glucose and gluconate dehydrogenases, sequentially transforms glucose into gluconate and 2-ketogluconate. Although the secretion of these organic acids by P. putida has been extensively recognized, the mechanism and spatiotemporal regulation of the PGS remained elusive thus far. To address this challenge, we adopted a dynamic 13C- and 2H-metabolic flux analysis strategy, termed D-fluxomics. D-fluxomics demonstrated that the PGS underscores a highly dynamic metabolic architecture in glucose-dependent batch cultures of P. putida, characterized by hierarchical carbon uptake by the PGS throughout the cultivation. Additionally, we show that gluconate and 2-ketogluconate accumulation and consumption can be solely explained as a result of the interplay between growth rate-coupled and decoupled metabolic fluxes. As a consequence, the formation of these acids in the PGS is inversely correlated to the bacterial growth rate-unlike the widely studied overflow metabolism of Escherichia coli and yeast. Our findings, which underline survival strategies of soil bacteria thriving in their natural environments, open new avenues for engineering P. putida towards efficient, sugar-based bioprocesses.

A neural-mechanistic hybrid approach improving the predictive power of genome-scale metabolic models

Constraint-based metabolic models have been used for decades to predict the phenotype of microorganisms in different environments. However, quantitative predictions are limited unless labor-intensive measurements of media uptake fluxes are performed. We show how hybrid neural-mechanistic models can serve as an architecture for machine learning providing a way to improve phenotype predictions. We illustrate our hybrid models with growth rate predictions of Escherichia coli and Pseudomonas putida grown in different media and with phenotype predictions of gene knocked-out Escherichia coli mutants. Our neural-mechanistic models systematically outperform constraint-based models and require training set sizes orders of magnitude smaller than classical machine learning methods. Our hybrid approach opens a doorway to enhancing constraint-based modeling: instead of constraining mechanistic models with additional experimental measurements, our hybrid models grasp the power of machine learning while fulfilling mechanistic constrains, thus saving time and resources in typical systems biology or biological engineering projects.

A pathogen-specific isotope tracing approach reveals metabolic activities and fluxes of intracellular Salmonella

Pathogenic bacteria proliferating inside mammalian host cells need to rapidly adapt to the intracellular environment. How they achieve this and scavenge essential nutrients from the host has been an open question due to the difficulties in distinguishing between bacterial and host metabolites in situ. Here, we capitalized on the inability of mammalian cells to metabolize mannitol to develop a stable isotopic labeling approach to track Salmonella enterica metabolites during intracellular proliferation in host macrophage and epithelial cells. By measuring label incorporation into Salmonella metabolites with liquid chromatography-mass spectrometry (LC-MS), and combining it with metabolic modeling, we identify relevant carbon sources used by Salmonella, uncover routes of their metabolization, and quantify relative reaction rates in central carbon metabolism. Our results underline the importance of the Entner-Doudoroff pathway (EDP) and the phosphoenolpyruvate carboxylase for intracellularly proliferating Salmonella. More broadly, our metabolic labeling strategy opens novel avenues for understanding the metabolism of pathogens inside host cells.

Ammonium sulfate supplementation enhances erythromycin biosynthesis by augmenting intracellular metabolism and precursor supply in Saccharopolyspora erythraea

In this study, the cellular metabolic mechanisms regarding ammonium sulfate supplementation on erythromycin production were investigated by employing targeted metabolomics and metabolic flux analysis. The results suggested that the addition of ammonium sulfate stimulates erythromycin biosynthesis. Targeted metabolomics analysis uncovered that the addition of ammonium sulfate during the late stage of fermentation resulted in an augmented intracellular amino acid metabolism pool, guaranteeing an ample supply of precursors for organic acids and coenzyme A-related compounds. Therefore, adequate precursors facilitated cellular maintenance and erythromycin biosynthesis. Subsequently, an optimal supplementation rate of 0.02 g/L/h was determined. The results exhibited that erythromycin titer (1311.1 μg/mL) and specific production rate (0.008 mmol/gDCW/h) were 101.3% and 41.0% higher than those of the process without ammonium sulfate supplementation, respectively. Moreover, the erythromycin A component proportion increased from 83.2% to 99.5%. Metabolic flux analysis revealed increased metabolic fluxes with the supplementation of three ammonium sulfate rates.

Real-Time Assessment of Intracellular Metabolites in Single Cells through RNA-Based Sensors

Quantification of the concentration of particular cellular metabolites reports on the actual utilization of metabolic pathways in physiological and pathological conditions. Metabolite concentration also constitutes the readout for screening cell factories in metabolic engineering. However, there are no direct approaches that allow for real-time assessment of the levels of intracellular metabolites in single cells. In recent years, the modular architecture of natural bacterial RNA riboswitches has inspired the design of genetically encoded synthetic RNA devices that convert the intracellular concentration of a metabolite into a quantitative fluorescent signal. These so-called RNA-based sensors are composed of a metabolite-binding RNA aptamer as the sensor domain, connected through an actuator segment to a signal-generating reporter domain. However, at present, the variety of available RNA-based sensors for intracellular metabolites is still very limited. Here, we go through natural mechanisms for metabolite sensing and regulation in cells across all kingdoms, focusing on those mediated by riboswitches. We review the design principles underlying currently developed RNA-based sensors and discuss the challenges that hindered the development of novel sensors and recent strategies to address them. We finish by introducing the current and potential applicability of synthetic RNA-based sensors for intracellular metabolites.

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.

One-shot 13 C15 N-metabolic flux analysis for simultaneous quantification of carbon and nitrogen flux

Metabolic flux is the final output of cellular regulation and has been extensively studied for carbon but much less is known about nitrogen, which is another important building block for living organisms. For the tuberculosis pathogen, this is particularly important in informing the development of effective drugs targeting the pathogen's metabolism. Here we performed ¹³ C¹⁵ N dual isotopic labeling of Mycobacterium bovis BCG steady state cultures, quantified intracellular carbon and nitrogen fluxes and inferred reaction bidirectionalities. This was achieved by model scope extension and refinement, implemented in a multi-atom transition model, within the statistical framework of Bayesian model averaging (BMA). Using BMA-based ¹³ C¹⁵ N-metabolic flux analysis, we jointly resolve carbon and nitrogen fluxes quantitatively. We provide the first nitrogen flux distributions for amino acid and nucleotide biosynthesis in mycobacteria and establish glutamate as the central node for nitrogen metabolism. We improved resolution of the notoriously elusive anaplerotic node in central carbon metabolism and revealed possible operation modes. Our study provides a powerful and statistically rigorous platform to simultaneously infer carbon and nitrogen metabolism in any biological system.

Stochastic simulation algorithm for isotope-based dynamic flux analysis

Carbon isotope labeling method is a standard metabolic engineering tool for flux quantification in living cells. To cope with the high dimensionality of isotope labeling systems, diverse algorithms have been developed to reduce the number of variables or operations in metabolic flux analysis (MFA), but lacks generalizability to non-stationary metabolic conditions. In this study, we present a stochastic simulation algorithm (SSA) derived from the chemical master equation of the isotope labeling system. This algorithm allows to compute the time evolution of isotopomer concentrations in non-stationary conditions, with the valuable property that computational time does not scale with the number of isotopomers. The efficiency and limitations of the algorithm is benchmarked for the forward and inverse problems of 13C-DMFA in the pentose phosphate pathways, and is compared with EMU-based methods for NMFA and MFA including the central carbon metabolism. Overall, SSA constitutes an alternative class to deterministic approaches for metabolic flux analysis that is well adapted to comprehensive dataset including parallel labeling experiments, and whose limitations associated to the sampling size can be overcome by using Monte Carlo sampling approaches.

Metabolic Engineering: Methodologies and Applications

Metabolic engineering aims to improve the production of economically valuable molecules through the genetic manipulation of microbial metabolism. While the discipline is a little over 30 years old, advancements in metabolic engineering have given way to industrial-level molecule production benefitting multiple industries such as chemical, agriculture, food, pharmaceutical, and energy industries. This review describes the design, build, test, and learn steps necessary for leading a successful metabolic engineering campaign. Moreover, we highlight major applications of metabolic engineering, including synthesizing chemicals and fuels, broadening substrate utilization, and improving host robustness with a focus on specific case studies. Finally, we conclude with a discussion on perspectives and future challenges related to metabolic engineering.

Clinical stratification improves the diagnostic accuracy of small omics datasets within machine learning and genome-scale metabolic modelling methods

BACKGROUND

Recently, multi-omic machine learning architectures have been proposed for the early detection of cancer. However, for rare cancers and their associated small datasets, it is still unclear how to use the available multi-omics data to achieve a mechanistic prediction of cancer onset and progression, due to the limited data available. Hepatoblastoma is the most frequent liver cancer in infancy and childhood, and whose incidence has been lately increasing in several developed countries. Even though some studies have been conducted to understand the causes of its onset and discover potential biomarkers, the role of metabolic rewiring has not been investigated in depth so far.

METHODS

Here, we propose and implement an interpretable multi-omics pipeline that combines mechanistic knowledge from genome-scale metabolic models with machine learning algorithms, and we use it to characterise the underlying mechanisms controlling hepatoblastoma.

RESULTS AND CONCLUSIONS

While the obtained machine learning models generally present a high diagnostic classification accuracy, our results show that the type of omics combinations used as input to the machine learning models strongly affects the detection of important genes, reactions and metabolic pathways linked to hepatoblastoma. Our method also suggests that, in the context of computer-aided diagnosis of cancer, optimal diagnostic accuracy can be achieved by adopting a combination of omics that depends on the patient's clinical characteristics.

Isotopic Tracer for Absolute Quantification of Metabolites of the Pentose Phosphate Pathway in Bacteria

The pentose phosphate pathway (PPP) plays a key role in many metabolic functions, including the generation of NADPH, biosynthesis of nucleotides, and carbon homeostasis. In particular, the intermediates of PPP have been found to be significantly perturbed in bacterial metabolomic studies. Nonetheless, detailed analysis to gain mechanistic information of PPP metabolism remains limited as most studies are unable to report on the absolute levels of the metabolites. Absolute quantification of metabolites is a prerequisite to study the details of fluxes and its regulations. Isotope tracer or labeling studies are conducted in vivo and in vitro and have significantly improved the analysis and understanding of PPP. Due to the laborious procedure and limitations in the in vivo method, an in vitro approach known as Group Specific Internal Standard Technology (GSIST) has been successfully developed to measure the absolute levels of central carbon metabolism, including PPP. The technique adopts derivatization of an experimental sample and a corresponding internal standard with isotope-coded reagents to provide better precision for accurate identification and absolute quantification. In this review, we highlight bacterial studies that employed isotopic tracers as the tagging agents used for the absolute quantification analysis of PPP metabolites.

Isotope-Assisted Metabolic Flux Analysis: A Powerful Technique to Gain New Insights into the Human Metabolome in Health and Disease

Cell metabolism represents the coordinated changes in genes, proteins, and metabolites that occur in health and disease. The metabolic fluxome, which includes both intracellular and extracellular metabolic reaction rates (fluxes), therefore provides a powerful, integrated description of cellular phenotype. However, intracellular fluxes cannot be directly measured. Instead, flux quantification requires sophisticated mathematical and computational analysis of data from isotope labeling experiments. In this review, we describe isotope-assisted metabolic flux analysis (iMFA), a rigorous computational approach to fluxome quantification that integrates metabolic network models and experimental data to generate quantitative metabolic flux maps. We highlight practical considerations for implementing iMFA in mammalian models, as well as iMFA applications in in vitro and in vivo studies of physiology and disease. Finally, we identify promising new frontiers in iMFA which may enable us to fully unlock the potential of iMFA in biomedical research.

Single cell metabolism: current and future trends

Single cell metabolomics is an emerging and rapidly developing field that complements developments in single cell analysis by genomics and proteomics. Major goals include mapping and quantifying the metabolome in sufficient detail to provide useful information about cellular function in highly heterogeneous systems such as tissue, ultimately with spatial resolution at the individual cell level. The chemical diversity and dynamic range of metabolites poses particular challenges for detection, identification and quantification. In this review we discuss both significant technical issues of measurement and interpretation, and progress toward addressing them, with recent examples from diverse biological systems. We provide a framework for further directions aimed at improving workflow and robustness so that such analyses may become commonly applied, especially in combination with metabolic imaging and single cell transcriptomics and proteomics.

In-cell NMR: Why and how?

NMR spectroscopy has been applied to cells and tissues analysis since its beginnings, as early as 1950. We have attempted to gather here in a didactic fashion the broad diversity of data and ideas that emerged from NMR investigations on living cells. Covering a large proportion of the periodic table, NMR spectroscopy permits scrutiny of a great variety of atomic nuclei in all living organisms non-invasively. It has thus provided quantitative information on cellular atoms and their chemical environment, dynamics, or interactions. We will show that NMR studies have generated valuable knowledge on a vast array of cellular molecules and events, from water, salts, metabolites, cell walls, proteins, nucleic acids, drugs and drug targets, to pH, redox equilibria and chemical reactions. The characterization of such a multitude of objects at the atomic scale has thus shaped our mental representation of cellular life at multiple levels, together with major techniques like mass-spectrometry or microscopies. NMR studies on cells has accompanied the developments of MRI and metabolomics, and various subfields have flourished, coined with appealing names: fluxomics, foodomics, MRI and MRS (i.e. imaging and localized spectroscopy of living tissues, respectively), whole-cell NMR, on-cell ligand-based NMR, systems NMR, cellular structural biology, in-cell NMR… All these have not grown separately, but rather by reinforcing each other like a braided trunk. Hence, we try here to provide an analytical account of a large ensemble of intricately linked approaches, whose integration has been and will be key to their success. We present extensive overviews, firstly on the various types of information provided by NMR in a cellular environment (the "why", oriented towards a broad readership), and secondly on the employed NMR techniques and setups (the "how", where we discuss the past, current and future methods). Each subsection is constructed as a historical anthology, showing how the intrinsic properties of NMR spectroscopy and its developments structured the accessible knowledge on cellular phenomena. Using this systematic approach, we sought i) to make this review accessible to the broadest audience and ii) to highlight some early techniques that may find renewed interest. Finally, we present a brief discussion on what may be potential and desirable developments in the context of integrative studies in biology.

13C metabolic flux analysis: Classification and characterization from the perspective of mathematical modeling and application in physiological research of neural cell

¹³C metabolic flux analysis (¹³C-MFA) has emerged as a forceful tool for quantifying in vivo metabolic pathway activity of different biological systems. This technology plays an important role in understanding intracellular metabolism and revealing patho-physiology mechanism. Recently, it has evolved into a method family with great diversity in experiments, analytics, and mathematics. In this review, we classify and characterize the various branch of ¹³C-MFA from a unified perspective of mathematical modeling. By linking different parts in the model to each step of its workflow, the specific technologies of ¹³C-MFA are put into discussion, including the isotope labeling model (ILM), isotope pattern measuring technique, optimization algorithm and statistical method. Its application in physiological research in neural cell has also been reviewed.

Fluxomics - New Metabolomics Approaches to Monitor Metabolic Pathways

Fluxomics is an innovative -omics research field that measures the rates of all intracellular fluxes in the central metabolism of biological systems. Fluxomics gathers data from multiple different -omics fields, portraying the whole picture of molecular interactions. Recently, fluxomics has become one of the most relevant approaches to investigate metabolic phenotypes. Metabolic flux using ¹³C-labeled molecules is increasingly used to monitor metabolic pathways, to probe the corresponding gene-RNA and protein-metabolite interaction networks in actual time. Thus, fluxomics reveals the functioning of multi-molecular metabolic pathways and is increasingly applied in biotechnology and pharmacology. Here, we describe the main fluxomics approaches and experimental platforms. Moreover, we summarize recent fluxomic results in different biological systems.

Quantitative metabolic fluxes regulated by trans-omic networks

Cells change their metabolism in response to internal and external conditions by regulating the trans-omic network, which is a global biochemical network with multiple omic layers. Metabolic flux is a direct measure of the activity of a metabolic reaction that provides valuable information for understanding complex trans-omic networks. Over the past decades, techniques to determine metabolic fluxes, including 13C-metabolic flux analysis (13C-MFA), flux balance analysis (FBA), and kinetic modeling, have been developed. Recent studies that acquire quantitative metabolic flux and multi-omic data have greatly advanced the quantitative understanding and prediction of metabolism-centric trans-omic networks. In this review, we present an overview of 13C-MFA, FBA, and kinetic modeling as the main techniques to determine quantitative metabolic fluxes, and discuss their advantages and disadvantages. We also introduce case studies with the aim of understanding complex metabolism-centric trans-omic networks based on the determination of metabolic fluxes.

Metabolomics as a promising tool for improving understanding of Multiple Sclerosis: a review of recent advances

Multiple sclerosis (MS) is an inflammatory demyelinating disease of the central nervous system that usually affects young adults. The development of MS is closely related to the changes in the metabolome. Metabolomics studies have been performed using biofluids or tissue samples from rodent models and human patients to reveal metabolic alterations associated with MS progression. This review aims to provide an overview of the applications of metabolomics that for the investigations of the perturbed metabolic pathways in MS and to reveal the potential of metabolomics in personalizing treatments. In conclusion, informative variations of metabolites can be potential biomarkers in advancing our understanding of MS pathogenesis for MS diagnosis, predicting the progression of the disease, and estimating drug effects. Metabolomics will be a promising technique for improving clinical care in MS.

Probing the Biogenesis of Polysaccharide Granules in Algal Cells at Sub-Organellar Resolution via Raman Microscopy with Stable Isotope Labeling

Phototrophs assimilate CO₂ into organic compounds that accumulate in storage organelles. Elucidation of the carbon dynamics of storage organelles could enhance the production efficiency of valuable compounds and facilitate the screening of strains with high photosynthetic activity. To comprehensively elucidate the carbon dynamics of these organelles, the intraorganellar distribution of the carbon atoms that accumulate at specific time periods should be probed. In this study, the biosynthesis of polysaccharides in storage organelles was spatiotemporally probed via stimulated Raman scattering (SRS) microscopy using a stable isotope (¹³C) as the tracking probe. Paramylon granules (a storage organelle of β-1,3-glucan) accumulated in a unicellular photosynthetic alga, Euglena gracilis, were investigated as a model organelle. The carbon source of the culture medium was switched from NaH¹²CO₃ to NaH¹³CO₃ during the production of the paramylon granules; this resulted in the distribution of the ¹²C and ¹³C constituents in the granules, so that the biosynthetic process could be tracked. Taking advantage of high-resolution SRS imaging and label switching, the localization of the ¹²C and ¹³C constituents inside a single paramylon granule could be visualized in three dimensions, thus revealing the growth process of paramylon granules. We propose that this method can be used for comprehensive elucidation of the dynamic activities of storage organelles.

Nuclear Magnetic Resonance Spectroscopy in Clinical Metabolomics and Personalized Medicine: Current Challenges and Perspectives

Personalized medicine is probably the most promising area being developed in modern medicine. This approach attempts to optimize the therapies and the patient care based on the individual patient characteristics. Its success highly depends on the way the characterization of the disease and its evolution, the patient’s classification, its follow-up and the treatment could be optimized. Thus, personalized medicine must combine innovative tools to measure, integrate and model data. Towards this goal, clinical metabolomics appears as ideally suited to obtain relevant information. Indeed, the metabolomics signature brings crucial insight to stratify patients according to their responses to a pathology and/or a treatment, to provide prognostic and diagnostic biomarkers, and to improve therapeutic outcomes. However, the translation of metabolomics from laboratory studies to clinical practice remains a subsequent challenge. Nuclear magnetic resonance spectroscopy (NMR) and mass spectrometry (MS) are the two key platforms for the measurement of the metabolome. NMR has several advantages and features that are essential in clinical metabolomics. Indeed, NMR spectroscopy is inherently very robust, reproducible, unbiased, quantitative, informative at the structural molecular level, requires little sample preparation and reduced data processing. NMR is also well adapted to the measurement of large cohorts, to multi-sites and to longitudinal studies. This review focus on the potential of NMR in the context of clinical metabolomics and personalized medicine. Starting with the current status of NMR-based metabolomics at the clinical level and highlighting its strengths, weaknesses and challenges, this article also explores how, far from the initial “opposition” or “competition”, NMR and MS have been integrated and have demonstrated a great complementarity, in terms of sample classification and biomarker identification. Finally, a perspective discussion provides insight into the current methodological developments that could significantly raise NMR as a more resolutive, sensitive and accessible tool for clinical applications and point-of-care diagnosis. Thanks to these advances, NMR has a strong potential to join the other analytical tools currently used in clinical settings.

A guide to interrogating immunometabolism

The metabolic charts memorized in early biochemistry courses, and then later forgotten, have come back to haunt many immunologists with new recognition of the importance of these pathways. Metabolites and the activity of metabolic pathways drive energy production, macromolecule synthesis, intracellular signalling, post-translational modifications and cell survival. Immunologists who identify a metabolic phenotype in their system are often left wondering where to begin and what does it mean? Here, we provide a framework for navigating and selecting the appropriate biochemical techniques to explore immunometabolism. We offer recommendations for initial approaches to develop and test metabolic hypotheses and how to avoid common mistakes. We then discuss how to take things to the next level with metabolomic approaches, such as isotope tracing and genetic approaches. By proposing strategies and evaluating the strengths and weaknesses of different methodologies, we aim to provide insight, note important considerations and discuss ways to avoid common misconceptions. Furthermore, we highlight recent studies demonstrating the power of these metabolic approaches to uncover the role of metabolism in immunology. By following the framework in this Review, neophytes and seasoned investigators alike can venture into the emerging realm of cellular metabolism and immunity with confidence and rigour.

Integrating systemic and molecular levels to infer key drivers sustaining metabolic adaptations

Metabolic adaptations to complex perturbations, like the response to pharmacological treatments in multifactorial diseases such as cancer, can be described through measurements of part of the fluxes and concentrations at the systemic level and individual transporter and enzyme activities at the molecular level. In the framework of Metabolic Control Analysis (MCA), ensembles of linear constraints can be built integrating these measurements at both systemic and molecular levels, which are expressed as relative differences or changes produced in the metabolic adaptation. Here, combining MCA with Linear Programming, an efficient computational strategy is developed to infer additional non-measured changes at the molecular level that are required to satisfy these constraints. An application of this strategy is illustrated by using a set of fluxes, concentrations, and differentially expressed genes that characterize the response to cyclin-dependent kinases 4 and 6 inhibition in colon cancer cells. Decreases and increases in transporter and enzyme individual activities required to reprogram the measured changes in fluxes and concentrations are compared with down-regulated and up-regulated metabolic genes to unveil those that are key molecular drivers of the metabolic response.

Deciphering the essential events in the reprogramming of metabolic networks subjected to complex perturbations, including the response to pharmacological treatments in multifactorial diseases like cancer, is crucial for the design of efficient therapies. Yet, tools to infer the molecular drivers sustaining such metabolic responses remain elusive for large metabolic networks. Here we develop an efficient computational strategy that integrates measured changes at systemic and molecular levels and combines metabolic control analysis with linear programming tools to infer key molecular drivers sustaining the metabolic adaptations to complex perturbations, such as an antitumoral drug therapy. The collective behavior is approximated using linear expressions where the adaptation of systemic concentrations and fluxes to a perturbation is described as a function of the molecular reprogramming of transport and enzyme activities. Starting from measured changes in fluxes and concentrations, we identify changes in the reprogramming of transporter and enzyme activities that are required to orchestrate the metabolic adaptation of colon cancer cells to a cell cycle inhibitor.

Charting the Metabolic Landscape of the Facultative Methylotroph Bacillus methanolicus

Methanol is inexpensive, is easy to transport, and can be produced both from renewable and from fossil resources without mobilizing arable lands. As such, it is regarded as a potential carbon source to transition toward a greener industrial chemistry. Metabolic engineering of bacteria and yeast able to efficiently consume methanol is expected to provide cell factories that will transform methanol into higher-value chemicals in the so-called methanol economy. Toward that goal, the study of natural methylotrophs such as Bacillus methanolicus is critical to understand the origin of their efficient methylotrophy. This knowledge will then be leveraged to transform such natural strains into new cell factories or to design methylotrophic capability in other strains already used by the industry.

Bacillus methanolicus MGA3 is a thermotolerant and relatively fast-growing methylotroph able to secrete large quantities of glutamate and lysine. These natural characteristics make B. methanolicus a good candidate to become a new industrial chassis organism, especially in a methanol-based economy. Intriguingly, the only substrates known to support B. methanolicus growth as sole sources of carbon and energy are methanol, mannitol, and, to a lesser extent, glucose and arabitol. Because fluxomics provides the most direct readout of the cellular phenotype, we hypothesized that comparing methylotrophic and nonmethylotrophic metabolic states at the flux level would yield new insights into MGA3 metabolism. In this study, we designed and performed a ¹³C metabolic flux analysis (¹³C-MFA) of the facultative methylotroph B. methanolicus MGA3 growing on methanol, mannitol, and arabitol to compare the associated metabolic states. On methanol, results showed a greater flux in the ribulose monophosphate (RuMP) pathway than in the tricarboxylic acid (TCA) cycle, thus validating previous findings on the methylotrophy of B. methanolicus. New insights related to the utilization of cyclic RuMP versus linear dissimilation pathways and between the RuMP variants were generated. Importantly, we demonstrated that the linear detoxification pathways and the malic enzyme shared with the pentose phosphate pathway have an important role in cofactor regeneration. Finally, we identified, for the first time, the metabolic pathway used to assimilate arabitol. Overall, those data provide a better understanding of this strain under various environmental conditions.

IMPORTANCE Methanol is inexpensive, is easy to transport, and can be produced both from renewable and from fossil resources without mobilizing arable lands. As such, it is regarded as a potential carbon source to transition toward a greener industrial chemistry. Metabolic engineering of bacteria and yeast able to efficiently consume methanol is expected to provide cell factories that will transform methanol into higher-value chemicals in the so-called methanol economy. Toward that goal, the study of natural methylotrophs such as Bacillus methanolicus is critical to understand the origin of their efficient methylotrophy. This knowledge will then be leveraged to transform such natural strains into new cell factories or to design methylotrophic capability in other strains already used by the industry.

Metabolomic Approaches to Study Chemical Exposure-Related Metabolism Alterations in Mammalian Cell Cultures

Biological organisms are constantly exposed to an immense repertoire of molecules that cover environmental or food-derived molecules and drugs, triggering a continuous flow of stimuli-dependent adaptations. The diversity of these chemicals as well as their concentrations contribute to the multiplicity of induced effects, including activation, stimulation, or inhibition of physiological processes and toxicity. Metabolism, as the foremost phenotype and manifestation of life, has proven to be immensely sensitive and highly adaptive to chemical stimuli. Therefore, studying the effect of endo- or xenobiotics over cellular metabolism delivers valuable knowledge to apprehend potential cellular activity of individual molecules and evaluate their acute or chronic benefits and toxicity. The development of modern metabolomics technologies such as mass spectrometry or nuclear magnetic resonance spectroscopy now offers unprecedented solutions for the rapid and efficient determination of metabolic profiles of cells and more complex biological systems. Combined with the availability of well-established cell culture techniques, these analytical methods appear perfectly suited to determine the biological activity and estimate the positive and negative effects of chemicals in a variety of cell types and models, even at hardly detectable concentrations. Metabolic phenotypes can be estimated from studying intracellular metabolites at homeostasis in vivo, while in vitro cell cultures provide additional access to metabolites exchanged with growth media. This article discusses analytical solutions available for metabolic phenotyping of cell culture metabolism as well as the general metabolomics workflow suitable for testing the biological activity of molecular compounds. We emphasize how metabolic profiling of cell supernatants and intracellular extracts can deliver valuable and complementary insights for evaluating the effects of xenobiotics on cellular metabolism. We note that the concepts and methods discussed primarily for xenobiotics exposure are widely applicable to drug testing in general, including endobiotics that cover active metabolites, nutrients, peptides and proteins, cytokines, hormones, vitamins, etc.

A mechanism-aware and multiomic machine-learning pipeline characterizes yeast cell growth

Metabolic modeling and machine learning are key components in the emerging next generation of systems and synthetic biology tools, targeting the genotype-phenotype-environment relationship. Rather than being used in isolation, it is becoming clear that their value is maximized when they are combined. However, the potential of integrating these two frameworks for omic data augmentation and integration is largely unexplored. We propose, rigorously assess, and compare machine-learning-based data integration techniques, combining gene expression profiles with computationally generated metabolic flux data to predict yeast cell growth. To this end, we create strain-specific metabolic models for 1,143 Saccharomyces cerevisiae mutants and we test 27 machine-learning methods, incorporating state-of-the-art feature selection and multiview learning approaches. We propose a multiview neural network using fluxomic and transcriptomic data, showing that the former increases the predictive accuracy of the latter and reveals functional patterns that are not directly deducible from gene expression alone. We test the proposed neural network on a further 86 strains generated in a different experiment, therefore verifying its robustness to an additional independent dataset. Finally, we show that introducing mechanistic flux features improves the predictions also for knockout strains whose genes were not modeled in the metabolic reconstruction. Our results thus demonstrate that fusing experimental cues with in silico models, based on known biochemistry, can contribute with disjoint information toward biologically informed and interpretable machine learning. Overall, this study provides tools for understanding and manipulating complex phenotypes, increasing both the prediction accuracy and the extent of discernible mechanistic biological insights.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

Resolving Metabolic Heterogeneity in Experimental Models of the Tumor Microenvironment from a Stable Isotope Resolved Metabolomics Perspective

The tumor microenvironment (TME) comprises complex interactions of multiple cell types that determines cell behavior and metabolism such as nutrient competition and immune suppression. We discuss the various types of heterogeneity that exist in solid tumors, and the complications this invokes for studies of TME. As human subjects and in vivo model systems are complex and difficult to manipulate, simpler 3D model systems that are compatible with flexible experimental control are necessary for studying metabolic regulation in TME. Stable Isotope Resolved Metabolomics (SIRM) is a valuable tool for tracing metabolic networks in complex systems, but at present does not directly address heterogeneous metabolism at the individual cell level. We compare the advantages and disadvantages of different model systems for SIRM experiments, with a focus on lung cancer cells, their interactions with macrophages and T cells, and their response to modulators in the immune microenvironment. We describe the experimental set up, illustrate results from 3D cultures and co-cultures of lung cancer cells with human macrophages, and outline strategies to address the heterogeneous TME.

Parallel isotope differential modeling for instationary 13C fluxomics at the genome scale

BACKGROUND

A precise map of the metabolic fluxome, the closest surrogate to the physiological phenotype, is becoming progressively more important in the metabolic engineering of photosynthetic organisms for biofuel and biomass production. For photosynthetic organisms, the state-of-the-art method for this purpose is instationary 13C fluxomics, which has arisen as a sibling of transcriptomics or proteomics. Instationary 13C data processing requires solving high-dimensional nonlinear differential equations and leads to large computational and time costs when its scope is expanded to a genome-scale metabolic network.

RESULT

Here, we present a parallelized method to model instationary 13C labeling data. The elementary metabolite unit (EMU) framework is reorganized to allow treating individual mass isotopomers and breaking up of their networks into strongly connected components (SCCs). A variable domain parallel algorithm is introduced to process ordinary differential equations in a parallel way. 15-fold acceleration is achieved for constant-step-size modeling and ~ fivefold acceleration for adaptive-step-size modeling.

CONCLUSION

This algorithm is universally applicable to isotope granules such as EMUs and cumomers and can substantially accelerate instationary 13C fluxomics modeling. It thus has great potential to be widely adopted in any instationary 13C fluxomics modeling.

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For ¹³C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

Stable Isotopes for Tracing Mammalian-Cell Metabolism In Vivo

Metabolism is at the cornerstone of all cellular functions and mounting evidence of its deregulation in different diseases emphasizes the importance of a comprehensive understanding of metabolic regulation at the whole-organism level. Stable-isotope measurements are a powerful tool for probing cellular metabolism and, as a result, are increasingly used to study metabolism in in vivo settings. The additional complexity of in vivo metabolic measurements requires paying special attention to experimental design and data interpretation. Here, we review recent work where in vivo stable-isotope measurements have been used to address relevant biological questions within an in vivo context, summarize different experimental and data interpretation approaches and their limitations, and discuss future opportunities in the field.

NMR-based metabolomics and fluxomics: developments and future prospects

Recent NMR developments are acting as game changers for metabolomics and fluxomics – a critical and perspective review.

Metabolic flux ratio analysis by parallel 13C labeling of isoprenoid biosynthesis in Rhodobacter sphaeroides

Metabolic engineering for increased isoprenoid production often benefits from the simultaneous expression of the two naturally available isoprenoid metabolic routes, namely the 2-methyl-D-erythritol 4-phosphate (MEP) pathway and the mevalonate (MVA) pathway. Quantification of the contribution of these pathways to the overall isoprenoid production can help to obtain a better understanding of the metabolism within a microbial cell factory. Such type of investigation can benefit from ¹³C metabolic flux ratio studies. Here, we designed a method based on parallel labeling experiments (PLEs), using [1-¹³C]- and [4-¹³C]glucose as tracers to quantify the metabolic flux ratios in the glycolytic and isoprenoid pathways. By just analyzing a reporter isoprenoid molecule and employing only four equations, we could describe the metabolism involved from substrate catabolism to product formation. These equations infer ¹³C atom incorporation into the universal isoprenoid building blocks, isopentenyl-pyrophosphate (IPP) and dimethylallyl-pyrophosphate (DMAPP). Therefore, this renders the method applicable to the study of any of isoprenoid of interest. As proof of principle, we applied it to study amorpha-4,11-diene biosynthesis in the bacterium Rhodobacter sphaeroides. We confirmed that in this species the Entner-Doudoroff pathway is the major pathway for glucose catabolism, while the Embden-Meyerhof-Parnas pathway contributes to a lesser extent. Additionally, we demonstrated that co-expression of the MEP and MVA pathways caused a mutual enhancement of their metabolic flux capacity. Surprisingly, we also observed that the isoprenoid flux ratio remains constant under exponential growth conditions, independently from the expression level of the MVA pathway. Apart from proposing and applying a tool for studying isoprenoid biosynthesis within a microbial cell factory, our work reveals important insights from the co-expression of MEP and MVA pathways, including the existence of a yet unclear interaction between them.

Separation of 13C- and 15N-Isotopologues of Amino Acids with a Primary Amine without Mass Resolution by Means of O-Phthalaldehyde Derivatization and Collision Induced Dissociation

Computational and experimental advances of recent years have culminated in establishing ¹³C-Metabolic Flux Analysis (¹³C-MFA) as a routine methodology to unravel the fluxome. As the acronym suggests, ¹³C-MFA has relied on the relative abundance of ¹³C-isotopes in metabolites for flux inference, most commonly measured by mass spectrometry. In this manuscript we expand the scope of labeling measurements to the case of simultaneous ¹³C- and ¹⁵N-labeling of amino acids. Analytically, the separation of isotopologues of this metabolite class can only be achieved at resolving power beyond 65,000. In this manuscript we harvest an overlooked property of the collision induced dissociation of amino acid adducts to discern ¹³C- and ¹⁵N- isotopologues of amino acids with a primary amine without separating them in the m/z domain.

Tear Down the Fluorescent Curtain: A New Fluorescence Suppression Method for Raman Microspectroscopic Analyses

The near exponential proliferation of published Raman microspectroscopic applications over the last decade bears witness to the strengths and versatility of this technology. However, laser-induced fluorescence often severely impedes its application to biological samples. Here we report a new approach for near complete elimination of laser-induced background fluorescence in highly pigmented biological specimens (e.g., microalgae) enabling interrogation by Raman microspectroscopy. Our simple chemiphotobleaching method combines mild hydrogen peroxide oxidation with broad spectrum visible light irradiation of the entire specimen. This treatment permits observing intracellular distributions of macromolecular pools, isotopic tracers, and even viral propagation within cells previously not amenable to Raman microspectroscopic examination. Our approach demonstrates the potential for confocal Raman microspectroscopy becoming an indispensable tool to obtain spatially-resolved data on the chemical composition of highly fluorescent biological samples from individual cells to environmental samples.

Interoperable and scalable data analysis with microservices: applications in metabolomics

MOTIVATION

Developing a robust and performant data analysis workflow that integrates all necessary components whilst still being able to scale over multiple compute nodes is a challenging task. We introduce a generic method based on the microservice architecture, where software tools are encapsulated as Docker containers that can be connected into scientific workflows and executed using the Kubernetes container orchestrator.

RESULTS

We developed a Virtual Research Environment (VRE) which facilitates rapid integration of new tools and developing scalable and interoperable workflows for performing metabolomics data analysis. The environment can be launched on-demand on cloud resources and desktop computers. IT-expertise requirements on the user side are kept to a minimum, and workflows can be re-used effortlessly by any novice user. We validate our method in the field of metabolomics on two mass spectrometry, one nuclear magnetic resonance spectroscopy and one fluxomics study. We showed that the method scales dynamically with increasing availability of computational resources. We demonstrated that the method facilitates interoperability using integration of the major software suites resulting in a turn-key workflow encompassing all steps for mass-spectrometry-based metabolomics including preprocessing, statistics and identification. Microservices is a generic methodology that can serve any scientific discipline and opens up for new types of large-scale integrative science.

AVAILABILITY AND IMPLEMENTATION

The PhenoMeNal consortium maintains a web portal (https://portal.phenomenal-h2020.eu) providing a GUI for launching the Virtual Research Environment. The GitHub repository https://github.com/phnmnl/ hosts the source code of all projects.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

p13CMFA: Parsimonious 13C metabolic flux analysis

Deciphering the mechanisms of regulation of metabolic networks subjected to perturbations, including disease states and drug-induced stress, relies on tracing metabolic fluxes. One of the most informative data to predict metabolic fluxes are ¹³C based metabolomics, which provide information about how carbons are redistributed along central carbon metabolism. Such data can be integrated using ¹³C Metabolic Flux Analysis (¹³C MFA) to provide quantitative metabolic maps of flux distributions. However, ¹³C MFA might be unable to reduce the solution space towards a unique solution either in large metabolic networks or when small sets of measurements are integrated. Here we present parsimonious ¹³C MFA (p¹³CMFA), an approach that runs a secondary optimization in the ¹³C MFA solution space to identify the solution that minimizes the total reaction flux. Furthermore, flux minimization can be weighted by gene expression measurements allowing seamless integration of gene expression data with ¹³C data. As proof of concept, we demonstrate how p¹³CMFA can be used to estimate intracellular flux distributions from ¹³C measurements and transcriptomics data. We have implemented p¹³CMFA in Iso2Flux, our in-house developed isotopic steady-state ¹³C MFA software. The source code is freely available on GitHub (https://github.com/cfoguet/iso2flux/releases/tag/0.7.2).

¹³C Metabolic Flux Analysis (¹³C MFA) is a well-established technique that has proven to be a valuable tool in quantifying the metabolic flux profile of central carbon metabolism. When a biological system is incubated with a ¹³C-labeled substrate, ¹³C propagates to metabolites throughout the metabolic network in a flux and pathway-dependent manner. ¹³C MFA integrates measurements of ¹³C enrichment in metabolites to identify the flux distributions consistent with the measured ¹³C propagation. However, there is often a range of flux values that can lead to the observed ¹³C distribution. Indeed, either when the metabolic network is large or a small set of measurements are integrated, the range of valid solutions can be too wide to accurately estimate part of the underlying flux distribution. Here we propose to use flux minimization to select the best flux solution in the¹³C MFA solution space. Furthermore, this approach can integrate gene expression data to give greater weight to the minimization of fluxes through enzymes with low gene expression evidence in order to ensure that the selected solution is biologically relevant. The concept of using flux minimization to select the best solution is widely used in flux balance analysis, but it had never been applied in the framework of ¹³C MFA. We have termed this new approach parsimonious ¹³C MFA (p¹³CMFA).

Machine and deep learning meet genome-scale metabolic modeling

Omic data analysis is steadily growing as a driver of basic and applied molecular biology research. Core to the interpretation of complex and heterogeneous biological phenotypes are computational approaches in the fields of statistics and machine learning. In parallel, constraint-based metabolic modeling has established itself as the main tool to investigate large-scale relationships between genotype, phenotype, and environment. The development and application of these methodological frameworks have occurred independently for the most part, whereas the potential of their integration for biological, biomedical, and biotechnological research is less known. Here, we describe how machine learning and constraint-based modeling can be combined, reviewing recent works at the intersection of both domains and discussing the mathematical and practical aspects involved. We overlap systematic classifications from both frameworks, making them accessible to nonexperts. Finally, we delineate potential future scenarios, propose new joint theoretical frameworks, and suggest concrete points of investigation for this joint subfield. A multiview approach merging experimental and knowledge-driven omic data through machine learning methods can incorporate key mechanistic information in an otherwise biologically-agnostic learning process.

The Design of FluxML: A Universal Modeling Language for 13C Metabolic Flux Analysis

¹³C metabolic flux analysis (MFA) is the method of choice when a detailed inference of intracellular metabolic fluxes in living organisms under metabolic quasi-steady state conditions is desired. Being continuously developed since two decades, the technology made major contributions to the quantitative characterization of organisms in all fields of biotechnology and health-related research. ¹³C MFA, however, stands out from other "-omics sciences," in that it requires not only experimental-analytical data, but also mathematical models and a computational toolset to infer the quantities of interest, i.e., the metabolic fluxes. At present, these models cannot be conveniently exchanged between different labs. Here, we present the implementation-independent model description language FluxML for specifying ¹³C MFA models. The core of FluxML captures the metabolic reaction network together with atom mappings, constraints on the model parameters, and the wealth of data configurations. In particular, we describe the governing design processes that shaped the FluxML language. We demonstrate the utility of FluxML to represent many contemporary experimental-analytical requirements in the field of ¹³C MFA. The major aim of FluxML is to offer a sound, open, and future-proof language to unambiguously express and conserve all the necessary information for model re-use, exchange, and comparison. Along with FluxML, several powerful computational tools are supplied for easy handling, but also to maintain a maximum of flexibility. Altogether, the FluxML collection is an "all-around carefree package" for ¹³C MFA modelers. We believe that FluxML improves scientific productivity as well as transparency and therewith contributes to the efficiency and reproducibility of computational modeling efforts in the field of ¹³C MFA.

Systems Metabolic Engineering Meets Machine Learning: A New Era for Data-Driven Metabolic Engineering

The recent increase in high-throughput capacity of 'omics datasets combined with advances and interest in machine learning (ML) have created great opportunities for systems metabolic engineering. In this regard, data-driven modeling methods have become increasingly valuable to metabolic strain design. In this review, the nature of 'omics is discussed and a broad introduction to the ML algorithms combining these datasets into predictive models of metabolism and metabolic rewiring is provided. Next, this review highlights recent work in the literature that utilizes such data-driven methods to inform various metabolic engineering efforts for different classes of application including product maximization, understanding and profiling phenotypes, de novo metabolic pathway design, and creation of robust system-scale models for biotechnology. Overall, this review aims to highlight the potential and promise of using ML algorithms with metabolic engineering and systems biology related datasets.

PhenoMeNal: processing and analysis of metabolomics data in the cloud

BACKGROUND

Metabolomics is the comprehensive study of a multitude of small molecules to gain insight into an organism's metabolism. The research field is dynamic and expanding with applications across biomedical, biotechnological, and many other applied biological domains. Its computationally intensive nature has driven requirements for open data formats, data repositories, and data analysis tools. However, the rapid progress has resulted in a mosaic of independent, and sometimes incompatible, analysis methods that are difficult to connect into a useful and complete data analysis solution.

FINDINGS

PhenoMeNal (Phenome and Metabolome aNalysis) is an advanced and complete solution to set up Infrastructure-as-a-Service (IaaS) that brings workflow-oriented, interoperable metabolomics data analysis platforms into the cloud. PhenoMeNal seamlessly integrates a wide array of existing open-source tools that are tested and packaged as Docker containers through the project's continuous integration process and deployed based on a kubernetes orchestration framework. It also provides a number of standardized, automated, and published analysis workflows in the user interfaces Galaxy, Jupyter, Luigi, and Pachyderm.

CONCLUSIONS

PhenoMeNal constitutes a keystone solution in cloud e-infrastructures available for metabolomics. PhenoMeNal is a unique and complete solution for setting up cloud e-infrastructures through easy-to-use web interfaces that can be scaled to any custom public and private cloud environment. By harmonizing and automating software installation and configuration and through ready-to-use scientific workflow user interfaces, PhenoMeNal has succeeded in providing scientists with workflow-driven, reproducible, and shareable metabolomics data analysis platforms that are interfaced through standard data formats, representative datasets, versioned, and have been tested for reproducibility and interoperability. The elastic implementation of PhenoMeNal further allows easy adaptation of the infrastructure to other application areas and 'omics research domains.

Cellular determinants of metabolite concentration ranges

Cellular functions are shaped by reaction networks whose dynamics are determined by the concentrations of underlying components. However, cellular mechanisms ensuring that a component’s concentration resides in a given range remain elusive. We present network properties which suffice to identify components whose concentration ranges can be efficiently computed in mass-action metabolic networks. We show that the derived ranges are in excellent agreement with simulations from a detailed kinetic metabolic model of Escherichia coli. We demonstrate that the approach can be used with genome-scale metabolic models to arrive at predictions concordant with measurements from Escherichia coli under different growth scenarios. By application to 14 genome-scale metabolic models from diverse species, our approach specifies the cellular determinants of concentration ranges that can be effectively employed to make predictions for a variety of biotechnological and medical applications.

We present a computational approach for inferring concentration ranges from genome-scale metabolic models. The approach specifies a determinant and molecular mechanism underling facile control of concentration ranges for components in large-scale cellular networks. Most importantly, the predictions about concentration ranges do not require knowledge of kinetic parameters (which are difficult to specify at a genome scale), provided measurements of concentrations in a reference state. The approach assumes that reaction rates follow the mass action law used in the derivations of other types of kinetics. We apply the approach with large-scale kinetic and stoichiometric metabolic models of organisms from different kingdoms of life to show that we can identify a proportion of metabolites to which our approach is applicable. By challenging the predictions of concentration ranges in the genome-scale metabolic network of E. coli with real-world data sets, we further demonstrate the prediction power and limitations of the approach.

Assessing the Impact of the Nutrient Microenvironment on the Metabolism of Effector CD8+ T Cells

Immune cell function is tightly regulated by cellular metabolism, which in turn is strongly linked to the nutrient availability in the microenvironment surrounding the cells. This link is critical for effector CD8⁺ T cells which, after activation, must migrate from nutrient-rich environments into nutrient-scarce regions such as the tumor microenvironment. Assessing how nutrient availability modulates the metabolism of effector CD8⁺ T cells is thus key for understanding how harsh environments may impair their proliferation and effector function. Here, we describe an approach to systematically study the impact of the nutrient microenvironment on the metabolism of effector CD8⁺ T cells, based on performing stable ¹³C isotope labeling measurements on in vitro-differentiated murine effector CD8⁺ T cells.

A metabolome-wide characterization of the diabetic phenotype in ZDF rats and its reversal by pioglitazone

Type 2 diabetes (T2D) is a complex metabolic disease associated with alterations in glucose, lipid and protein metabolism. In order to characterize the biochemical phenotype of the Zucker diabetic fatty (ZDF) rat, the most common animal model for the study of T2D, and the impact of the insulin sensitizer pioglitazone, a global, mass spectrometry-based analysis of the metabolome was conducted. Overall, 420 metabolites in serum, 443 in the liver and 603 in the intestine were identified at study end. In comparison to two control groups, obese diabetic ZDF rats showed characteristic metabolic signatures that included hyperglycemia, elevated β-oxidation, dyslipidemia—featured by an increase in saturated and monounsaturated fatty acids and a decrease of medium chain and of polyunsaturated fatty acids in serum–and decreased amino acid levels, consistent with their utilization in hepatic gluconeogenesis. A 13-week treatment with the PPARγ agonist pioglitazone reversed most of these signatures: Pioglitazone improved glycemic control and the fatty acid profile, elevated amino acid levels in the liver, but decreased branched chain amino acids in serum. The hitherto most comprehensive metabolic profiling study identified a biochemical blueprint for the ZDF diabetic model and captured the impact of genetic, nutritional and pharmacological perturbations. The in-depth characterization on the molecular level deepens the understanding and further validates the ZDF rat as a suitable preclinical model of diabetes in humans.