SUMOFLUX: A Generalized Method for Targeted 13C Metabolic Flux Ratio Analysis

Machine learning for the advancement of genome-scale metabolic modeling

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

Functional link hybrid artificial neural network for predicting continuous biohydrogen production in dynamic membrane bioreactor

Conventional machine learning approaches have shown limited predictive power when applied to continuous biohydrogen production due to nonlinearity and instability. This study was aimed at forecasting the dynamic membrane reactor performance in terms of the hydrogen production rate (HPR) and hydrogen yield (HY) using laboratory-based daily operation datapoints for twelve input variables. Hybrid algorithms were developed by integrating particle swarm optimized with functional link artificial neural network (PSO-FLN) which outperformed other hybrid algorithms for both HPR and HY, with determination coefficients (R2) of 0.97 and 0.80 and mean absolute percentage errors of 0.014 % and 0.023 %, respectively. Shapley additive explanations (SHAP) explained the two positive-influencing parameters, OLR_added (1.1-1.3 mol/L/d) and butyric acid (7.5-16.5 g COD/L) supports the highest HPR (40-60 L/L/d). This research indicates that PSO-FLN model are capable of handling complicated datasets with high precision in less computational timeat 9.8 sec for HPR and 10.0 sec for HY prediction.

PyMIDA: A Graphical User Interface for Mass Isotopomer Distribution Analysis

Mass isotopomer distribution analysis (MIDA) is an analytical technique that measures the synthesis rate of biological polymers using combinatorial probabilities and stable isotope labeling. Over the past few decades, this method has been developed and applied to a wide range of uses that have increased our understanding of metabolism and the etiology and monitoring of disease. There is currently no publicly available piece of software for performing MIDA calculations in a targeted manner without its functionality being limited to a specific use case. We present a cross-platform Python graphical user interface implementation for research to obtain kinetic parameters easily from stable-isotope labeling studies and provide the code and user manual on GitHub.

FreeFlux: A Python Package for Time-Efficient Isotopically Nonstationary Metabolic Flux Analysis

13C metabolic flux analysis is a powerful tool for metabolism characterization in metabolic engineering and synthetic biology. However, the widespread adoption of this tool is hindered by limited software availability and computational efficiency. Currently, the most widely accepted 13C-flux tools, such as INCA and 13CFLUX2, are developed in a closed-source environment. While several open-source packages or software are available, they are either computationally inefficient or only suitable for flux estimation at isotopic steady state. To address the need for a time-efficient computational tool for the more complicated flux analysis at an isotopically nonstationary state, especially for understanding the single-carbon substrate metabolism, we present FreeFlux. FreeFlux is an open-source Python package that performs labeling pattern simulation and flux analysis at both isotopic steady state and transient state, enabling a more comprehensive analysis of cellular metabolism. FreeFlux provides a set of interfaces to manipulate the objects abstracted from a labeling experiment and computational process, making it easy to integrate into other programs or pipelines. The flux estimation by FreeFlux is fast and reliable, and its validity has been confirmed by comparison with results from other computational tools using both synthetic and experimental data. FreeFlux is freely available at https://github.com/Chaowu88/freeflux with a detailed online tutorial and documentation provided at https://freeflux.readthedocs.io/en/latest/index.html.

Recent advances in machine learning applications in metabolic engineering

Metabolic engineering encompasses several widely-used strategies, which currently hold a high seat in the field of biotechnology when its potential is manifesting through a plethora of research and commercial products with a strong societal impact. The genomic revolution that occurred almost three decades ago has initiated the generation of large omics-datasets which has helped in gaining a better understanding of cellular behavior. The itinerary of metabolic engineering that has occurred based on these large datasets has allowed researchers to gain detailed insights and a reasonable understanding of the intricacies of biosystems. However, the existing trail-and-error approaches for metabolic engineering are laborious and time-intensive when it comes to the production of target compounds with high yields through genetic manipulations in host organisms. Machine learning (ML) coupled with the available metabolic engineering test instances and omics data brings a comprehensive and multidisciplinary approach that enables scientists to evaluate various parameters for effective strain design. This vast amount of biological data should be standardized through knowledge engineering to train different ML models for providing accurate predictions in gene circuits designing, modification of proteins, optimization of bioprocess parameters for scaling up, and screening of hyper-producing robust cell factories. This review briefs on the premise of ML, followed by mentioning various ML methods and algorithms alongside the numerous omics datasets available to train ML models for predicting metabolic outcomes with high-accuracy. The combinative interplay between the ML algorithms and biological datasets through knowledge engineering have guided the recent advancements in applications such as CRISPR/Cas systems, gene circuits, protein engineering, metabolic pathway reconstruction, and bioprocess engineering. Finally, this review addresses the probable challenges of applying ML in metabolic engineering which will guide the researchers toward novel techniques to overcome the limitations.

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.

Metabolic flux analysis: a comprehensive review on sample preparation, analytical techniques, data analysis, computational modelling, and main application areas

Metabolic flux analysis (MFA) quantitatively describes cellular fluxes to understand metabolic phenotypes and functional behaviour after environmental and/or genetic perturbations. In the last decade, the application of stable isotopes became extremely important to determine and integrate in vivo measurements of metabolic reactions in systems biology. 13C-MFA is one of the most informative methods used to study central metabolism of biological systems. This review aims to outline the current experimental procedure adopted in 13C-MFA, starting from the preparation of cell cultures and labelled tracers to the quenching and extraction of metabolites and their subsequent analysis performed with very powerful software. Here, the limitations and advantages of nuclear magnetic resonance spectroscopy and mass spectrometry techniques used in carbon labelled experiments are elucidated by reviewing the most recent published papers. Furthermore, we summarise the most successful approaches used for computational modelling in flux analysis and the main application areas with a particular focus in metabolic engineering.

Isotope Calculation Gadgets: A Series of Software for Isotope-Tracing Experiments in Garuda Platform

Isotope tracing is a powerful technique for elucidating intracellular metabolism. Experiments utilizing this technique involve various processes, such as the correction of natural isotopes. Although some previously developed software are available for these procedures, there are still time-consuming steps in isotope tracing including the creation of an isotope measurement method in mass spectrometry (MS) and the interpretation of obtained labeling data. Additionally, these multi-step tasks often require data format conversion, which is also time-consuming. In this study, the Isotope Calculation Gadgets, a series of software that supports an entire workflow of isotope-tracing experiments, was developed in the Garuda platform, an open community. Garuda is a graphical user interface-based platform that allows individual operations to be sequentially performed, without data format conversion, which significantly reduces the required time and effort. The developed software includes new features that construct channels for isotopomer measurements, as well as conventional functions such as natural isotope correction, the calculation of fractional labeling and split ratio, and data mapping, thus facilitating an overall workflow of isotope-tracing experiments through smooth functional integration.

Hot isopropanol quenching procedure for automated microtiter plate scale 13C-labeling experiments

Background

Currently, the generation of genetic diversity for microbial cell factories outpaces the screening of strain variants with omics-based phenotyping methods. Especially isotopic labeling experiments, which constitute techniques aimed at elucidating cellular phenotypes and supporting rational strain design by growing microorganisms on substrates enriched with heavy isotopes, suffer from comparably low throughput and the high cost of labeled substrates.

Results

We present a miniaturized, parallelized, and automated approach to ¹³C-isotopic labeling experiments by establishing and validating a hot isopropanol quenching method on a robotic platform coupled with a microbioreactor cultivation system. This allows for the first time to conduct automated labeling experiments at a microtiter plate scale in up to 48 parallel batches. A further innovation enabled by the automated quenching method is the analysis of free amino acids instead of proteinogenic ones on said microliter scale. Capitalizing on the latter point and as a proof of concept, we present an isotopically instationary labeling experiment in Corynebacterium glutamicum ATCC 13032, generating dynamic labeling data of free amino acids in the process.

Conclusions

Our results show that a robotic liquid handler is sufficiently fast to generate informative isotopically transient labeling data. Furthermore, the amount of biomass obtained from a sub-milliliter cultivation in a microbioreactor is adequate for the detection of labeling patterns of free amino acids. Combining the innovations presented in this study, isotopically stationary and instationary automated labeling experiments can be conducted, thus fulfilling the prerequisites for ¹³C-metabolic flux analyses in high-throughput.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01806-4.

Computational Framework for Machine-Learning-Enabled 13C Fluxomics

¹³C metabolic flux analysis (MFA) has emerged as a powerful tool for synthetic biology. This optimization-based approach suffers long computation time and unstable solutions depending on the initial guess. Here, we develop a machine-learning-based framework for ¹³C fluxomics. Specifically, training and test data sets are generated by metabolic network decomposition and flux sampling, in which flux ratios at metabolic nodes and simulated labeling patterns of metabolites are used as training targets and features, respectively. To improve prediction accuracy and simplify the model, automated processes are developed for flux ratio selection based on solvability and feature screening based on importance. We found that predictive performance can be significantly improved using both amino acids and central carbon metabolites in comparison with amino acids alone. Together with measured external fluxes, the predicted flux ratios determine the mass balance system, yielding global flux distributions. This approach is validated by flux estimation using both simulated and experimental data in comparison with canonical ¹³C MFA. The approach represents a reliable fluxomics method readily applicable to high-throughput metabolic phenotyping, which highlights the advances of intelligent learning algorithms in synthetic biology, specifically in the Test and Learn stage of the Design-Build-Test-Learn cycle.

Stable Isotopes for Tracing Cardiac Metabolism in Diseases

Although metabolic remodeling during cardiovascular diseases has been well-recognized for decades, the recent development of analytical platforms and mathematical tools has driven the emergence of assessing cardiac metabolism using tracers. Metabolism is a critical component of cellular functions and adaptation to stress. The pathogenesis of cardiovascular disease involves metabolic adaptation to maintain cardiac contractile function even in advanced disease stages. Stable-isotope tracer measurements are a powerful tool for measuring flux distributions at the whole organism level and assessing metabolic changes at a systems level in vivo. The goal of this review is to summarize techniques and concepts for in vivo or ex vivo stable isotope labeling in cardiovascular research, to highlight mathematical concepts and their limitations, to describe analytical methods at the tissue and single-cell level, and to discuss opportunities to leverage metabolic models to address important mechanistic questions relevant to all patients with cardiovascular disease.

Traditional Herbal Remedies Used for Managing Anxiety and Insomnia in Italy: An Ethnopharmacological Overview

Anxiety and insomnia are among the most common mental health disorders and are a major cause of disability around the world. Traditional herbal medicines are receiving significant attention in global health debates. Several Italian regions maintain rural traditions and are among the most extensively studied areas of Europe regarding medicinal plant uses. The present overview aims to highlight the use of wild and cultivated plants, specifically as sedatives and for insomnia treatment in Italy, and to collect, analyze, and summarize the available literature about their pharmacological activity as well as clinical and pre-clinical studies concerning the most cited plants. In total, 106 wild taxa are used in Italy for sedative purposes. The plant species belong to 76 genera and 32 families, of which the most cited are Asteraceae (24.2%) and Lamiaceae (21.1%). Leaves (29%) and flowers (27%) are the plant parts mostly used as infusion (70%) and decoction (25%). Out of 106 taxa documented, only the most cited are analyzed in this overview (A. arvensis L., C. nepeta L., C. monogyna Jacq., H. lupulus L., L. nobilis L., L. angustifolia Mill., M. sylvestris L., M. chamomilla L., M. officinalis L., O. basilicum L., P. rhoeas L., P. somniferum L., R. officinalis L., T. platyphyllus Scop., and V. officinalis L.). Among the fifteen species selected, only seven have been studied for their pharmacological activity as hypnotic-sedatives. Future pre-clinical and clinical studies are needed to better clarify the mechanism of action of bioactive compounds and confirm the potential of these alternative therapies.

Metabolomics in infectious diseases and drug discovery

Metabolomics has emerged as an invaluable tool that can be used along with genomics, transcriptomics and proteomics to understand host–pathogen interactions at small-molecule levels. Metabolomics has been used to study a variety of infectious diseases and applications. The most common application of metabolomics is for prognostic and diagnostic purposes, specifically the screening of disease-specific biomarkers by either NMR-based or mass spectrometry-based metabolomics. In addition, metabolomics is of great significance for the discovery of druggable metabolic enzymes and/or metabolic regulators through the use of state-of-the-art flux analysis, for example, via the elucidation of metabolic mechanisms. This review discusses the application of metabolomics technologies to biomarker screening, the discovery of drug targets in infectious diseases such as viral, bacterial and parasite infections and immunometabolomics, highlights the challenges associated with accessing metabolite compartmentalization and discusses the available tools for determining local metabolite concentrations.

Metabolomics has emerged as an invaluable tool that can be used along with genomics, transcriptomics and proteomics to understand host–pathogen interactions at small-molecule levels.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

The era of big data: Genome-scale modelling meets machine learning

With omics data being generated at an unprecedented rate, genome-scale modelling has become pivotal in its organisation and analysis. However, machine learning methods have been gaining ground in cases where knowledge is insufficient to represent the mechanisms underlying such data or as a means for data curation prior to attempting mechanistic modelling. We discuss the latest advances in genome-scale modelling and the development of optimisation algorithms for network and error reduction, intracellular constraining and applications to strain design. We further review applications of supervised and unsupervised machine learning methods to omics datasets from microbial and mammalian cell systems and present efforts to harness the potential of both modelling approaches through hybrid modelling.

IsoSearch: An Untargeted and Unbiased Metabolite and Lipid Isotopomer Tracing Strategy from HR-LC-MS/MS Datasets

Stable isotopic tracer analysis is a technique used to determine carbon or nitrogen atom incorporation into biological systems. A number of mass spectrometry based approaches have been developed for this purpose, including high-resolution tandem mass spectrometry (HR-LC-MS/MS), selected reaction monitoring (SRM) and parallel reaction monitoring (PRM). We have developed an approach for analyzing untargeted metabolomic and lipidomic datasets using high-resolution mass spectrometry with polarity switching and implemented our approach in the open-source R script IsoSearch and in Scaffold Elements software. Using our strategy, which requires an unlabeled reference dataset and isotope labeled datasets across various biological conditions, we traced metabolic isotopomer alterations in breast cancer cells (MCF-7) treated with the metabolic drugs 2-deoxy-glucose, 6-aminonicotinamide, compound 968, and rapamycin. Metabolites and lipids were first identified by the commercial software Scaffold Elements and LipidSearch, then IsoSearch successfully profiled the ¹³C-isotopomers extracted metabolites and lipids from ¹³C-glucose labeled MCF-7 cells. The results interpreted known models, such as glycolysis and pentose phosphate pathway inhibition, but also helped to discover new metabolic/lipid flux patterns, including a reactive oxygen species (ROS) defense mechanism induced by 6AN and triglyceride accumulation in rapamycin treated cells. The results suggest the IsoSearch/Scaffold Elements platform is effective for studying metabolic tracer analysis in diseases, drug metabolism, and metabolic engineering for both polar metabolites and non-polar lipids.

Enzymes as Parts in Need of Replacement - and How to Extend Their Working Life

Enzymes catalyze reactions in vivo at different rates and each enzyme molecule has a lifetime limit before it is degraded and replaced to enable catalysis to continue. Considering these rates together as a unitless ratio of catalytic cycles until replacement (CCR) provides a new quantitative tool to assess the replacement schedule of and energy investment into enzymes as they relate to function. Here, we outline the challenges of determining CCRs and new approaches to overcome them and then assess the CCRs of selected enzymes in bacteria and plants to reveal a range of seven orders of magnitude for this ratio. Modifying CCRs in plants holds promise to lower cellular costs, to tailor enzymes for particular environments, and to breed enzyme improvements for crop productivity.

Staring into the void: demystifying microbial metabolomics

Metabolites give us a window into the chemistry of microbes and are split into two subclasses: primary and secondary. Primary metabolites are required for life whereas secondary metabolites have historically been classified as those appearing after exponential growth and are not necessarily needed for survival. Many microbial species are estimated to produce hundreds of metabolites and can be affected by differing nutrients. Using various analytical techniques, metabolites can be directly detected in order to elucidate their biological significance. Currently, a single experiment can produce anywhere from megabytes to terabytes of data. This big data has motivated scientists to develop informatics tools to help target specific metabolites or sets of metabolites. Broadly, it is imperative to identify clear biological questions before embarking on a study of metabolites (metabolomics). For instance, studying the effect of a transposon insertion on phenazine biosynthesis in Pseudomonas is a very different from asking what molecules are present in a specific banana-derived strain of Pseudomonas. This review is meant to serve as a primer for a 'choose your own adventure' approach for microbiologists with limited mass spectrometry expertise, with a strong focus on liquid chromatography mass spectrometry based workflows developed or optimized within the past five years.

Machine Learning Applications for Mass Spectrometry-Based Metabolomics

The metabolome of an organism depends on environmental factors and intracellular regulation and provides information about the physiological conditions. Metabolomics helps to understand disease progression in clinical settings or estimate metabolite overproduction for metabolic engineering. The most popular analytical metabolomics platform is mass spectrometry (MS). However, MS metabolome data analysis is complicated, since metabolites interact nonlinearly, and the data structures themselves are complex. Machine learning methods have become immensely popular for statistical analysis due to the inherent nonlinear data representation and the ability to process large and heterogeneous data rapidly. In this review, we address recent developments in using machine learning for processing MS spectra and show how machine learning generates new biological insights. In particular, supervised machine learning has great potential in metabolomics research because of the ability to supply quantitative predictions. We review here commonly used tools, such as random forest, support vector machines, artificial neural networks, and genetic algorithms. During processing steps, the supervised machine learning methods help peak picking, normalization, and missing data imputation. For knowledge-driven analysis, machine learning contributes to biomarker detection, classification and regression, biochemical pathway identification, and carbon flux determination. Of important relevance is the combination of different omics data to identify the contributions of the various regulatory levels. Our overview of the recent publications also highlights that data quality determines analysis quality, but also adds to the challenge of choosing the right model for the data. Machine learning methods applied to MS-based metabolomics ease data analysis and can support clinical decisions, guide metabolic engineering, and stimulate fundamental biological discoveries.

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

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

Dynamic 13C Flux Analysis Captures the Reorganization of Adipocyte Glucose Metabolism in Response to Insulin

Cellular metabolism is dynamic, but quantifying non-steady metabolic fluxes by stable isotope tracers presents unique computational challenges. Here, we developed an efficient 13C-tracer dynamic metabolic flux analysis (13C-DMFA) framework for modeling central carbon fluxes that vary over time. We used B-splines to generalize the flux parameterization system and to improve the stability of the optimization algorithm. As proof of concept, we investigated how 3T3-L1 cultured adipocytes acutely metabolize glucose in response to insulin. Insulin rapidly stimulates glucose uptake, but intracellular pathways responded with differing speeds and magnitudes. Fluxes in lower glycolysis increased faster than those in upper glycolysis. Glycolysis fluxes rose disproportionally larger and faster than the tricarboxylic acid cycle, with lactate a primary glucose end product. The uncovered array of flux dynamics suggests that glucose catabolism is additionally regulated beyond uptake to help shunt glucose into appropriate pathways. This work demonstrates the value of using dynamic intracellular fluxes to understand metabolic function and pathway regulation.

Tracing metabolic flux through time and space with isotope labeling experiments

Metabolism is dynamic and must function in context-specific ways to adjust to changes in the surrounding cellular and ecological environment. When isotopic tracers are used, metabolite flow (i.e. metabolic flux) can be quantified through biochemical networks to assess metabolic pathway operation. The cellular activities considered across multiple tissues and organs result in the observed phenotype and can be analyzed to discover emergent, whole-system properties of biology and elucidate misconceptions about network operation. However, temporal and spatial challenges remain significant hurdles and require novel approaches and creative solutions. We survey current investigations in higher plant and animal systems focused on dynamic isotope labeling experiments, spatially resolved measurement strategies, and observations from re-analysis of our own studies that suggest prospects for future work. Related discoveries will be necessary to push the frontier of our understanding of metabolism to suggest novel solutions to cure disease and feed a growing future world population.

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.

The challenge and potential of photosynthesis: unique considerations for metabolic flux measurements in photosynthetic microorganisms

Photosynthetic microorganisms have the potential for sustainable production of chemical feedstocks and products but have had limited success due to a lack of tools and deeper understanding of metabolic pathway regulation. The application of instationary metabolic flux analysis (INST-MFA) to photosynthetic microorganisms has allowed researchers to quantify fluxes and identify bottlenecks and metabolic inefficiencies to improve strain performance or gain insight into cellular physiology. Additionally, flux measurements can also highlight deviations between measured and predicted fluxes, revealing weaknesses in metabolic models and highlighting areas where a lack of understanding still exists. In this review, we outline the experimental steps necessary to successfully perform photosynthetic flux experiments and analysis. We also discuss the challenges unique to photosynthetic microorganisms and how to account for them, including: light supply, quenching, concentration, extraction, analysis, and flux calculation. We hope that this will enable a larger number of researchers to successfully apply isotope assisted metabolic flux analysis (¹³C-MFA) to their favorite photosynthetic organism.

13C based proteinogenic amino acid (PAA) and metabolic flux ratio analysis of Lactococcus lactis reveals changes in pentose phosphate (PP) pathway in response to agitation and temperature related stresses

Lactococcus lactis subsp. cremoris MG1363 is an important starter culture for dairy fermentation. During industrial fermentations, L. lactis is constantly exposed to stresses that affect the growth and performance of the bacterium. Although the response of L. lactis to several stresses has been described, the adaptation mechanisms at the level of in vivo fluxes have seldom been described. To gain insights into cellular metabolism, ¹³C metabolic flux analysis and gas chromatography mass spectrometry (GC-MS) were used to measure the flux ratios of active pathways in the central metabolism of L. lactis when subjected to three conditions varying in temperature (30°C, 37°C) and agitation (with and without agitation at 150 rpm). Collectively, the concentrations of proteinogenic amino acids (PAAs) and free fatty acids (FAAs) were compared, and Pearson correlation analysis (r) was calculated to measure the pairwise relationship between PAAs. Branched chain and aromatic amino acids, threonine, serine, lysine and histidine were correlated strongly, suggesting changes in flux regulation in glycolysis, the pentose phosphate (PP) pathway, malic enzyme and anaplerotic reaction catalysed by pyruvate carboxylase (pycA). Flux ratio analysis revealed that glucose was mainly converted by glycolysis, highlighting the stability of L. lactis’ central carbon metabolism despite different conditions. Higher flux ratios through oxaloacetate (OAA) from pyruvate (PYR) reaction in all conditions suggested the activation of pyruvate carboxylate (pycA) in L. lactis, in response to acid stress during exponential phase. Subsequently, more significant flux ratio differences were seen through the oxidative and non-oxidative pentose phosphate (PP) pathways, malic enzyme, and serine and C1 metabolism, suggesting NADPH requirements in response to environmental stimuli. These reactions could play an important role in optimization strategies for metabolic engineering in L. lactis. Overall, the integration of systematic analysis of amino acids and flux ratio analysis provides a systems-level understanding of how L. lactis regulates central metabolism under various conditions.

Estimation of flux ratios without uptake or release data: Application to serine and methionine metabolism

Model-based metabolic flux analysis (MFA) using isotope-labeled substrates has provided great insight into intracellular metabolic activities across a host of organisms. One challenge with applying MFA in mammalian systems, however, is the need for absolute quantification of nutrient uptake, biomass composition, and byproduct release fluxes. Such measurements are often not feasible in complex culture systems or in vivo. One way to address this issue is to estimate flux ratios, the fractional contribution of a flux to a metabolite pool, which are independent of absolute measurements and yet informative for cellular metabolism. Prior work has focused on "local" estimation of a handful of flux ratios for specific metabolites and reactions. Here, we perform systematic, model-based estimation of all flux ratios in a metabolic network using isotope labeling data, in the absence of uptake/release data. In a series of examples, we investigate what flux ratios can be well estimated with reasonably tight confidence intervals, and contrast this with confidence intervals on normalized fluxes. We find that flux ratios can provide useful information on the metabolic state, and is complementary to normalized fluxes: for certain metabolic reactions, only flux ratios can be well estimated, while for others normalized fluxes can be obtained. Simulation studies of a large human metabolic network model suggest that estimation of flux ratios is technically feasible for complex networks, but additional studies on data from actual isotopomer labeling experiments are needed to validate these results. Finally, we experimentally study serine and methionine metabolism in cancer cells using flux ratios. We find that, in these cells, the methionine cycle is truncated with little remethylation from homocysteine, and polyamine synthesis in the absence of methionine salvage leads to loss of 5-methylthioadenosine, suggesting a new mode of overflow metabolism in cancer cells. This work highlights the potential for flux ratio analysis in the absence of absolute quantification, which we anticipate will be important for both in vitro and in vivo studies of cancer metabolism.

Recent advances in high-throughput 13C-fluxomics

The rise of high throughput (HT) strain engineering tools accompanying the area of synthetic biology is supporting the generation of a large number of microbial cell factories. A current bottleneck in process development is our limited capacity to rapidly analyze the metabolic state of the engineered strains, and in particular their intracellular fluxes. HT ¹³C-fluxomics workflows have not yet become commonplace, despite the existence of several HT tools at each of the required stages. This includes cultivation and sampling systems, analytics for isotopic analysis, and software for data processing and flux calculation. Here, we review recent advances in the field and highlight bottlenecks that must be overcome to allow the emergence of true HT ¹³C-fluxomics workflows.