Computational quantification of metabolic fluxes from a single isotope snapshot: application to an animal biopsy

• Published in: Bioinformatics Vol. 26; no. 5; pp. 653 - 660
• Main Authors: Binsl, Thomas W., Alders, David J.C., Heringa, Jaap, Groeneveld, A.B. Johan, van Beek, Johannes H.G.M.
• Format: Journal Article
• Language: English
• Published: Oxford Oxford University Press 01.03.2010
• Subjects: Algorithms; Animals; Biological and medical sciences; Biopsy; Carbon Isotopes; Citric Acid Cycle - physiology; Computer Simulation; Fundamental and applied biological sciences. Psychology; General aspects; Isotope Labeling - methods; Magnetic Resonance Spectroscopy; Mathematics in biology. Statistical analysis. Models. Metrology. Data processing in biology (general aspects); Oxygen Consumption; Flux; Biopsy; Application; Animal; Quantitative analysis
• ISSN: 1367-4803; 1460-2059; 1367-4811
• DOI: 10.1093/bioinformatics/btq018

Motivation: Quantitative determination of metabolic fluxes in single tissue biopsies is difficult. We report a novel analysis approach and software package for in vivo flux quantification using stable isotope labeling. Results: We developed a protocol based on brief, timed infusion of 13C isotope-enriched substrates for the tricarboxylic acid (TCA) cycle followed by quick freezing of tissue biopsies. NMR measurements of tissue extracts were used for flux estimation based on a computational model of carbon transitions between TCA cycle metabolites and related amino acids. To this end, we developed a computational framework in which metabolic systems can be flexibly assembled, simulated and analyzed. Flux parameters were quantified from NMR multiplets by a partial grid search followed by repeated Nelder–Mead optimizations implemented on a computer grid. We implemented a model of the TCA cycle and showed by extensive simulations that the timed infusion protocol reliably quantitates multiple fluxes. Experimental validation of the method was done in vivo on hearts of anesthetized pigs under two different conditions: basal state (n = 7) and cardiac stress caused by infusion of dobutamine (n = 7). About nine tissue samples (40–200 mg dry-weight) were taken per heart. TCA cycle flux was 6.11 ± 0.28 (SEM) μmol/min · gdw at baseline versus 9.29 ± 1.03 μmol/min · gdw for dobutamine stress. Oxygen consumption calculated from the TCA cycle flux and from ‘gold standard’ blood gas-based measurements were close, correlating with r=0.88 (P < 10−4). Spatial heterogeneity in metabolic fluxes is detectable amongst the small samples. We propose that our novel isotope snapshot methodology is suitable for flux measurements in biopsies in vivo. Availability: Non-profit organizations will, upon request, be granted a non-exclusive license to use the software for internal research and teaching purposes at no charge. A web interface for using the software on our computer grid is available under http://www.ibi.vu.nl/programs/ Contact: hans.van.beek@falw.vu.nl Supplementary information: Supplementary data are available at Bioinformatics online.