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

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 (13)C 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) micromol/min x gdw at baseline versus 9.29 +/- 1.03 micromol/min x 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/

[1]  S. Vatner,et al.  Limited transfer of cytosolic NADH into mitochondria at high cardiac workload. , 2004, American journal of physiology. Heart and circulatory physiology.

[2]  A. Sherry,et al.  TCA cycle kinetics in the rat heart by analysis of (13)C isotopomers using indirect (1)H. , 2001, American journal of physiology. Heart and circulatory physiology.

[3]  L. Quek,et al.  OpenFLUX: efficient modelling software for 13C-based metabolic flux analysis , 2009, Microbial cell factories.

[4]  H. Brunengraber,et al.  Partitioning of pyruvate between oxidation and anaplerosis in swine hearts. , 2000, American journal of physiology. Heart and circulatory physiology.

[5]  Heinrich Taegtmeyer,et al.  Substrate–Enzyme Competition Attenuates Upregulated Anaplerotic Flux Through Malic Enzyme in Hypertrophied Rat Heart and Restores Triacylglyceride Content: Attenuating Upregulated Anaplerosis in Hypertrophy , 2009, Circulation research.

[6]  E M Chance,et al.  Mathematical analysis of isotope labeling in the citric acid cycle with applications to 13C NMR studies in perfused rat hearts. , 1983, The Journal of biological chemistry.

[7]  W. Wiechert 13C metabolic flux analysis. , 2001, Metabolic engineering.

[8]  M. Chandler,et al.  Acute hibernation decreases myocardial pyruvate carboxylation and citrate release. , 2001, American journal of physiology. Heart and circulatory physiology.

[9]  A. Sherry,et al.  Effects of aminooxyacetate on glutamate compartmentation and TCA cycle kinetics in rat hearts. , 1998, American journal of physiology. Heart and circulatory physiology.

[10]  James Kennedy,et al.  Particle swarm optimization , 2002, Proceedings of ICNN'95 - International Conference on Neural Networks.

[11]  Vanhamme,et al.  Improved method for accurate and efficient quantification of MRS data with use of prior knowledge , 1997, Journal of magnetic resonance.

[12]  Christopher R. Myers,et al.  Universally Sloppy Parameter Sensitivities in Systems Biology Models , 2007, PLoS Comput. Biol..

[13]  T Soddemann,et al.  散逸粒子動力学:平衡および非平衡分子動力学シミュレーションのための有用なサーモスタット(原標題は英語) , 2003 .

[14]  David B. Searls A View from the Dark Side , 2007, PLoS Comput. Biol..

[15]  Jaap Heringa,et al.  FluxSimulator: An R Package to Simulate Isotopomer Distributions in Metabolic Networks , 2007 .

[16]  F. Prinzen,et al.  Blood flow distributions by microsphere deposition methods. , 2000, Cardiovascular research.

[17]  G. Stephanopoulos,et al.  Elementary metabolite units (EMU): a novel framework for modeling isotopic distributions. , 2007, Metabolic engineering.

[18]  Craig R Malloy,et al.  Analytical solutions for (13)C isotopomer analysis of complex metabolic conditions: substrate oxidation, multiple pyruvate cycles, and gluconeogenesis. , 2004, Metabolic engineering.

[19]  J B Bassingthwaighte,et al.  Fractal 15O-labeled water washout from the heart. , 1995, Circulation research.

[20]  R. Denton,et al.  Control of the tricarboxylate cycle and its interactions with glycolysis during acetate utilization in rat heart. , 1970, The Biochemical journal.

[21]  Wolfgang Wiechert,et al.  Experimental design principles for isotopically instationary 13C labeling experiments , 2006, Biotechnology and bioengineering.

[22]  Peter Schattner Automated Querying of Genome Databases , 2007, PLoS Comput. Biol..

[23]  Edvard A. Falch Dynamic determination of the thermal characteristics of fermentation tanks , 1968 .

[24]  J M O'Donnell,et al.  Dynamic 13C NMR analysis of oxidative metabolism in the in vivo canine myocardium. , 1993, The Journal of biological chemistry.

[25]  K. S. Brown,et al.  Statistical mechanical approaches to models with many poorly known parameters. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  Johannes H G M van Beek,et al.  A 13C NMR double-labeling method to quantitate local myocardial O2 consumption using frozen tissue samples. , 1999, American journal of physiology. Heart and circulatory physiology.

[27]  Craig R. Malloy,et al.  TCA cycle kinetics in the rat heart by analysis of 13C isotopomers using indirect 1H[13C] detection , 2001 .

[28]  Henri Brunengraber,et al.  Quantitative assessment of anaplerosis from propionate in pig heart in vivo. , 2003, American journal of physiology. Endocrinology and metabolism.

[29]  J. V. van Beek,et al.  Simple model analysis of 13C NMR spectra to measure oxygen consumption using frozen tissue samples. , 1998, Advances in experimental medicine and biology.

[30]  J K Kelleher,et al.  Flux estimation using isotopic tracers: common ground for metabolic physiology and metabolic engineering. , 2001, Metabolic engineering.

[31]  Johannes H G M van Beek,et al.  Myocardial O2 consumption in porcine left ventricle is heterogeneously distributed in parallel to heterogeneous O2 delivery. , 2004, American journal of physiology. Heart and circulatory physiology.

[32]  Silvia Kuhlmann,et al.  Novel expression hosts for complex secondary metabolite megasynthetases: Production of myxochromide in the thermopilic isolate Corallococcus macrosporus GT-2 , 2009, Microbial cell factories.

[33]  E M Chance,et al.  Calculation of Absolute Metabolic Flux and the Elucidation of the Pathways of Glutamate Labeling in Perfused Rat Heart by 13C NMR Spectroscopy and Nonlinear Least Squares Analysis (*) , 1995, The Journal of Biological Chemistry.

[34]  D. Challoner,et al.  Respiration in Myocardium , 1968, Nature.

[35]  A. Sherry,et al.  Analysis of tricarboxylic acid cycle of the heart using 13C isotope isomers. , 1990, The American journal of physiology.