GPTomics/bio-metabolomics-isotope-tracing

Version Compatibility

Reference examples tested with: isocor 2.2+, numpy 1.26+

Before using code patterns, verify installed versions match. If versions differ:

• Python: pip show <package> then help(module.function) to check signatures
• R: packageVersion('<pkg>') then ?function_name to verify parameters

If code throws ImportError, AttributeError, or TypeError, introspect the installed package and adapt the example to match the actual API rather than retrying.

Isotope Tracing / Stable-Isotope-Resolved Metabolomics

"What is this pathway actually doing, not just how much metabolite is there?" -> Feed a labeled tracer, then measure how label propagates into downstream metabolites as a mass-isotopomer distribution (MID) over time.

• Python: isocor.mscorrectors.MetaboliteCorrectorFactory().correct() (IsoCor) for natural-abundance correction
• R: accucor::natural_abundance_correction() (AccuCor) for high-resolution correction
• Modeling layer (separate discipline): INCA, 13CFLUX2, OpenFLUX for 13C-MFA / INST-MFA flux fitting

The Single Most Important Insight -- Labeling Reports Flux; Pool Size Does Not

A metabolite's concentration is how much is there; its labeling pattern (MID) is where the carbon came from and how fast it got there. These are independent measurements and frequently move in OPPOSITE directions: block a downstream-consuming enzyme and the intermediate pool rises (it backs up) while the labeling of downstream products falls (flux through them dropped). Reading the pool alone reports the opposite of the biology. An isotope-tracing experiment therefore answers a fundamentally different question than untargeted or targeted abundance profiling, and its analysis is dominated by two mandatory corrections that fabricate flux if skipped: natural-abundance / tracer-purity correction (raw isotopologue areas are NOT the labeling), and the steady-state assumption (a single MID is a snapshot whose meaning depends on whether labeling has plateaued). Fractional enrichment is concentration-independent (it is a ratio within one pool), which is why it survives the recovery/matrix problems that plague absolute quant -- but it says nothing about amount.

Core Concepts

| Term | Meaning | Why it matters | |---|---|---| | Tracer / tracee | The labeled substrate fed (tracer, e.g. U-13C6-glucose) vs the unlabeled endogenous pool (tracee) | The experiment measures how tracer atoms replace tracee atoms over time | | Isotopologue | A molecule differing only in number of heavy atoms (M+0, M+1, M+2 ...) | Resolved by MASS; this is what MS measures and what an MID counts | | Isotopomer | Same number of heavy atoms but at different POSITIONS (e.g. 1-13C vs 6-13C lactate) | Resolved by POSITION; needs NMR or positional tracers, NOT mass spectra alone | | MID (mass-isotopomer distribution) | The fractional vector of M+0, M+1, ... for one metabolite | The primary readout; its shape encodes which route carbon took | | Fractional / mean enrichment | Weighted-mean labeled-atom fraction = sum(i * MID_i) / n_atoms | One-number summary of how labeled a pool is; concentration-independent | | Atom transitions | The map of which substrate carbons land on which product carbons per reaction | Defines the expected MID for each pathway; the basis of flux models | | Metabolic steady state | Pool sizes constant in time | Required for classical MFA; if pools drift, plateau MIDs do not give fluxes | | Isotopic steady state | Labeling has equilibrated to a stable plateau | Classical MFA reads fluxes from the plateau; sampling before it is invalid |

Decision Tree by Scenario

| Goal / situation | Do | Why | |---|---|---| | Want amount/concentration, units, biomarker level | Use abundance profiling -> metabolomics/targeted-analysis | Pool size is not flux; tracing cannot give a concentration | | Want pathway ACTIVITY/route, central carbon metabolism | 13C tracing (U-13C6-glucose, 13C5-glutamine); measure MIDs | Labeling reports flux through the route the carbon took | | Trace nitrogen handling (transamination, urea, nucleotides) | 15N tracer (e.g. 15N2-glutamine, 15N-ammonia) | N-flux is invisible to a 13C tracer | | Distinguish two carbon entry points into one pool | Positional / partially-labeled tracer (e.g. 1,2-13C2-glucose) | The M+1 vs M+2 split of products separates PPP from glycolysis | | Fast-labeling small pools, cultured cells, clear metabolic steady state | Steady-state 13C-MFA from plateau MIDs (INCA, 13CFLUX2) | Plateau labeling + atom transitions -> flux estimates | | Slow labeling, large pools, autotrophs, primary/quiescent cells | INST-MFA from the labeling TIME COURSE (INCA) | Drops the isotopic-steady-state assumption; fits transient + pool sizes | | Have raw isotopologue areas (low-res QqQ / high-res Orbitrap) | Natural-abundance + purity correction FIRST (IsoCor / AccuCor) | Uncorrected MID is wrong by construction; see below | | Want genome-scale predicted flux without a tracer | systems-biology/flux-balance-analysis | FBA is constraint-based prediction, NOT empirical label measurement |

Natural-Abundance + Tracer-Purity Correction (mandatory)

Goal: Turn raw measured isotopologue areas into a true MID that reflects only tracer-derived label.

Approach: Even a fully unlabeled molecule shows an M+1, M+2 ladder because ~1.07% of carbon is naturally 13C (plus 15N, 2H, 18O, 34S, and derivatization Si). Build the natural-abundance ladder from the molecular (and derivative) formula, deconvolve it out, then correct for the tracer not being 100% isotopically pure. Feeding uncorrected areas to a flux model is the equivalent of reporting an uncalibrated peak area as a concentration.

import isocor

# corrector knows the formula's natural-abundance ladder and the tracer
corrector = isocor.mscorrectors.MetaboliteCorrectorFactory(
    'C6H12O6', tracer='13C',
    correct_NA_tracer=True,           # also strip the labeled element's own natural abundance
    tracer_purity=[0.01, 0.99])       # [unlabeled, labeled] per-position purity of the tracer

# raw measured areas M+0..M+6 for a partially labeled glucose pool
corrected_area, iso_fraction, residuum, mean_enrichment = corrector.correct(
    [50000., 8000., 12000., 3000., 1500., 6000., 25000.])
# iso_fraction is the corrected MID; mean_enrichment is fractional enrichment

High-resolution Orbitrap data resolves 13C from 15N/2H by exact mass, enabling a different (often simpler) correction; AccuCor (R) is tuned for that case:

library(accucor)
# El-MAVEN / MAVEN isotopologue table; Resolution is the instrument resolving power
corrected <- natural_abundance_correction(path = 'elmaven_export.xlsx',
                                          resolution = 100000, purity = 0.99)

Pick the corrector by tracer count and resolution: IsoCor handles any tracer at any resolution; AccuCor (single tracer) and AccuCor2 (dual 13C-15N / 13C-2H) target high-res. Verify the chosen tool's current argument names before running -- both APIs drift across versions.

Computing and Plotting an MID / Fractional Enrichment

Goal: Summarize a corrected isotopologue vector as an MID and one fractional-enrichment number, comparably across conditions.

Approach: Normalize corrected areas to sum 1 (the MID), then take the atom-weighted mean over isotopologue index divided by the number of tracer atoms.

import numpy as np

corrected = np.array([26000., 2200., 5600., 1200., 500., 2300., 12500.])
mid = corrected / corrected.sum()                              # M+0..M+n fractions
fractional_enrichment = np.sum(np.arange(len(mid)) * mid) / (len(mid) - 1)
# stacked-bar MID per condition is the standard visualization; never plot raw (uncorrected) areas

Steady-State Check (does the labeling number mean anything yet?)

Goal: Decide whether a measured MID may be read as flux-informative or is still a kinetic transient.

Approach: Sample labeling at several timepoints; isotopic steady state is reached when the MID stops changing (plateau). Only plateau MIDs license classical-MFA flux inference; a rising MID is kinetic data requiring INST-MFA.

import numpy as np

# fractional enrichment per timepoint (minutes) for one metabolite
t = np.array([0, 5, 15, 30, 60, 120])
fe = np.array([0.00, 0.18, 0.31, 0.39, 0.42, 0.43])
reached_plateau = abs(fe[-1] - fe[-2]) < 0.02     # <2% change between last points = plateau
# if not reached_plateau: the pool is still labeling -> use the full time course (INST-MFA), not one point

Per-Method Failure Modes

Skipping natural-abundance correction

• Trigger: Reporting or modeling raw isotopologue areas straight from El-MAVEN / Skyline.
• Mechanism: ~1.07% natural 13C (plus 15N, 2H, derivatization Si) creates an M+1/M+2 ladder on every molecule independent of the tracer; raw M+1 is mostly natural abundance for short-chain metabolites.
• Symptom: Apparent labeling in unlabeled controls; inflated M+1; flux fits with tight CIs that are simply wrong.
• Fix: Always run IsoCor/AccuCor with the correct formula (and derivative formula for GC-MS), tracer element, and tracer purity before any interpretation.

Assuming steady state when it is not reached

• Trigger: Inferring flux from a single early-timepoint MID.
• Mechanism: Classical MFA assumes both metabolic AND isotopic steady state; a transient MID encodes kinetics, not the flux plateau.
• Symptom: Fluxes that change with sampling time; large residuals; biologically implausible splits.
• Fix: Verify plateau across a time course, or switch to INST-MFA (INCA) which fits the transient and estimates pool sizes too.

Tracer impurity ignored

• Trigger: Treating a "U-13C6" tracer as 100% labeled.
• Mechanism: Per-position purity is ~99%, so a fraction of tracer molecules carry a 12C, distorting the fully-labeled isotopologue; the error compounds with atom count.
• Symptom: Fully-labeled isotopologue (M+n) systematically under-counted; enrichment biased low.
• Fix: Supply the measured tracer purity to the corrector (tracer_purity / purity).

Pool-size-vs-labeling confound

• Trigger: Concluding "flux changed" from a changed pool concentration (or vice versa).
• Mechanism: Pool and labeling are independent and can anticorrelate; a rising intermediate can mean LESS downstream flux.
• Symptom: Pool-based and label-based conclusions disagree; "activation" that is actually a backup.
• Fix: Interpret MID/enrichment for flux and concentration for amount separately; report both, never substitute one for the other.

Quench/extraction continuing turnover

• Trigger: Slow quench between harvest and metabolism arrest.
• Mechanism: High-turnover metabolites keep reacting post-harvest, scrambling labeling before extraction.
• Symptom: Variable, sample-dependent MIDs; collapsed nucleotide/energy-charge metabolites.
• Fix: Fast cold quench (-40 to -80 C aqueous methanol/acetonitrile); standardize and minimize harvest-to-quench time.

Quantitative Thresholds

| Threshold | Source | Rationale | |---|---|---| | 13C natural abundance ~1.07% | IUPAC isotopic composition | Sets the natural-abundance ladder corrected out of every MID | | Tracer purity ~99% per position | Vendor U-13C specs | Must be supplied to correction; compounds with atom count | | Isotopic-steady-state = <~2% MID change between timepoints | Convention | Below this, plateau reached; classical MFA licensed | | Quench at -40 to -80 C aqueous organic | Quenching literature (convention) | Arrests metabolism fast enough for high-turnover pools | | INST-MFA when labeling is slow / pools large / autotrophic | Cheah & Young 2018 | Isotopic steady state is unreachable in time, so fit the transient |

Common Errors

| Error / symptom | Cause | Solution | |---|---|---| | correct() length mismatch in IsoCor | Measurement vector is not n_tracer_atoms + 1 long | Pass M+0..M+n with n = count of tracer-element atoms in the formula | | Labeling appears in unlabeled control | No natural-abundance correction | Run IsoCor/AccuCor before interpreting | | M+n isotopologue under-reported | Tracer purity left at 1.0 | Set tracer_purity / purity to the measured value | | GC-MS MID still wrong after correction | Derivatization atoms (TMS/TBDMS Si, extra C) omitted | Provide the derivative formula to the corrector | | Flux estimates shift with sampling time | Isotopic steady state not reached | Use a time course + INST-MFA, not a single MID | | ValueError half-defined resolution in IsoCor | Gave mz_of_resolution/charge without resolution | Provide all high-res parameters together or none |

References

• Millard P, Delepine B, Guionnet M, Heuillet M, Bellvert F, Letisse F. 2019. IsoCor: isotope correction for high-resolution MS labeling experiments. Bioinformatics 35:4484-4487.
• Su X, Lu W, Rabinowitz JD. 2017. Metabolite Spectral Accuracy on Orbitraps. Analytical Chemistry 89:5940-5948.
• Clasquin MF, Melamud E, Rabinowitz JD. 2012. LC-MS Data Processing with MAVEN: A Metabolomic Analysis and Visualization Engine. Current Protocols in Bioinformatics 37:14.11.1-14.11.23.
• Cheah YE, Young JD. 2018. Isotopically nonstationary metabolic flux analysis (INST-MFA): putting theory into practice. Current Opinion in Biotechnology 54:80-87.
• Young JD. 2014. INCA: a computational platform for isotopically non-stationary metabolic flux analysis. Bioinformatics 30:1333-1335.
• Antoniewicz MR. 2018. A guide to 13C metabolic flux analysis for the cancer biologist. Experimental & Molecular Medicine 50:1-13.

Related Skills

• metabolomics/targeted-analysis - Absolute pool quantification and MRM/SRM mechanics
• metabolomics/xcms-preprocessing - Upstream LC-MS feature detection
• metabolomics/pathway-mapping - Pathway enrichment that interprets pools, not flux
• systems-biology/flux-balance-analysis - Constraint-based predicted flux, distinct from empirical tracing