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 Analysis of the Physiological Contributions of Glucose to the TCA Cycle.

The Quantitative Relationship Between Isotopic and Net Contributions of Lactate and Glucose to the Tricarboxylic Acid (TCA) Cycle

Glucose feeds the TCA cycle via circulating lactate

The Source of Glycolytic Intermediates in Mammalian Tissues

Quantitative Fluxomics of Circulating Metabolites

Pathway analysis of oxygen utilization and tricarboxylic acid cycle activity in Saccharomyces cerevisiae growing on glucose

Respiratory Metabolism in Heterotrophic Plant Cells as Revealed by Isotopic Labeling and Metabolic Flux Analysis

Quantitative determinants of aerobic glycolysis identify flux through the enzyme GAPDH as a limiting step

Glucose feeds the TCA cycle via environmental ethanol in fermenting yeast

Correlation between TCA cycle flux and glucose uptake rate during respiro-fermentative growth of Saccharomyces cerevisiae

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

Determiners of cell fates: the tricarboxylic acid cycle versus the citrate-malate shuttle.

Quantification of lactate from various metabolic pathways and quantification issues of lactate isotopologues and isotopmers

Energy metabolism and glutamate-glutamine cycle in the brain: a stoichiometric modeling perspective

Global determination of reaction rates and lipid turnover kinetics in Mus musculus

Quantitation of metabolic activity from isotope tracing data using automated methodology

Direct mitochondrial import of lactate supports resilient carbohydrate oxidation

Glucose sparing by glycogenolysis (GSG) determines the relationship between brain metabolism and neurotransmission

Quantification of nutrient fluxes during acute exercise in mice

A Flux Balance of Glucose Metabolism Clarifies the Requirements of the Warburg Effect.

[Lactate in the brain--without turning sour]