Abstract:Metabolic flux analysis (MFA) is a computational approach to deciphering labelling patterns based on machine learning principles. Differing from typical machine learning algorithms that train a model from known datasets to make predictions, the commonly used MFA algorithm trains a metabolic network with data from isotope tracing experiments and directly outputs the learned information — that is, all fluxes in the network that best fit data 3,5 (Fig. 1b). However, as a machine learning algorithm, current MFA methods often lack systematic evaluation and benchmarking, a standard practice in broader machine learning and artificial intelligence applications 6 . Issues such as algorithmic convergence, flux estimation accuracy and result robustness in MFA studies have been raised but remain largely unaddressed 3 , limiting the effectiveness and broader adoption of these automated tools in metabolic research. To advance the capabilities of MFA for complex metabolic networks and extensive isotope tracing datasets, we developed an automated analysis methodology alongside a large-scale metabolic network model. This model comprises over 100 fluxes across key pathways, including glycolysis, the tricarboxylic acid (TCA) cycle, the pentose phosphate pathway (PPP), one-carbon metabolism, and several amino acid (AA) biosynthetic pathways (Fig. 1c, Supplementary Methods). Compared to contemporary MFA tools 7,8,9 , a notable feature of our methodology is the incorporation of organelle compartmentalization, facilitating accurate quantification of exchange fluxes between mitochondria and cytosol in eukaryotic cells (Fig. 1c). While other tools typically require tens of minutes to obtain a solution 7,9 , our methodology can generate an optimized solution, with fluxes that accurately explain the labelling pattern from a 13 C tracing experiment on cultured cell lines, within 2 s on a desktop computer 10 (Supplementary Fig. 1a–d). Nonetheless, a challenge arose from the observation that these optimized solutions could diverge significantly, showing considerable variability in certain net fluxes even with similar loss values (Fig. 1d, Supplementary Fig. 1e,f).To address this problem, we developed an optimization-averaging algorithm that refines the computation process by selecting a subset of solutions with minimal loss (selected solutions) from the pool of optimized solutions and averaging them to produce a new, more stable solution set (averaged solutions) (Fig. 1e, Supplementary Methods). These solutions, along with those generated using the typical strategy used in contemporary software (Supplementary Methods, Supplementary Fig. 1c), were benchmarked using simulated 13 C tracing datasets generated from a known flux vector (Supplementary Fig. 2a). The results demonstrated that, relative to the benchmark, the optimization-averaging algorithm effectively reduced flux variability and improved the accuracy of the results in approximating the known flux, even with varying levels of data availability (Fig. 1f, Supplementary Figs. 2b–e and 3a–e).

Quantitative differences in anaerobic growth and metabolism of three different strains using C-Metabolic Flux Analysis

Comparative study on central metabolic fluxes of Bacillus megaterium strains in continuous culture using 13 C labelled substrates

Revealing differences in metabolic flux distributions between a mutant strain and its parent strain Gluconacetobacter xylinus CGMCC 2955.

High-resolution 13C metabolic flux analysis

Flux balance analysis reveals acetate metabolism modulates cyclic electron flow and alternative glycolytic pathways in Chlamydomonas reinhardtii

GC/MS-based 13C metabolic flux analysis resolves the parallel and cyclic photomixotrophic metabolism of Synechocystis sp. PCC 6803 and selected deletion mutants including the Entner-Doudoroff and phosphoketolase pathways

Metabolic mechanism of Saccharomyces cerevisiae under different physiological conditions based on non-stationary 13C metabolic flux analysis

Mutants lacking global regulators, fis and arcA, in Escherichia coli enhanced growth fitness under acetate metabolism by pathway reprogramming

Glycolysis/gluconeogenesis specialization in microbes is driven by biochemical constraints of flux sensing

A quantitative lens on anaerobic life: leveraging the state-of-the-art fluxomics approach to explore clostridial metabolism

Comparative Study of Two Saccharomyces Cerevisiae Strains with Kinetic Models at Genome-Scale.

Contextual Flexibility in Pseudomonas aeruginosa Central Carbon Metabolism during Growth in Single Carbon Sources

13C Metabolic Flux Analysis for Larger Scale Cultivation Using Gas Chromatography-Combustion-isotope Ratio Mass Spectrometry.

Metabolic flux analysis of the mixotrophic metabolisms in the green sulfur bacterium Chlorobaculum tepidum

Metabolic flux analysis of carbon balance in Lactobacillus strains

Investigating Glucose and Xylose Metabolism in Saccharomyces Cerevisiae and Scheffersomyces Stipitis Via 13C Metabolic Flux Analysis

Analysis of carbon metabolism in Escherichia coli strains with an inactive phosphotransferase system by (13)C labeling and NMR spectroscopy

Quantitation of metabolic activity from isotope tracing data using automated methodology

Metabolic flux analysis of pykF gene knockout Escherichia coli based on 13C-labeling experiments together with measurements of enzyme activities and intracellular metabolite concentrations

Different transcriptomic and metabolomic analysis of Saccharomyces cerevisiae BY4742 and CEN.PK2-1C strains

Catabolic Regulation Analysis of Escherichia Coli and Its Crp, Mlc, Mgsa, Pgi and Ptsg Mutants