OpenFLUX2: 13C-MFA modeling software package adjusted for the comprehensive analysis of single and parallel labeling experiments

被引:37
作者
Shupletsov, Mikhail S. [1 ,2 ]
Golubeva, Lyubov I. [1 ]
Rubina, Svetlana S. [1 ]
Podvyaznikov, Dmitry A. [1 ,3 ]
Iwatani, Shintaro [1 ]
Mashko, Sergey V. [1 ,3 ,4 ]
机构
[1] Ajinomoto Genet Res Inst, Moscow 117545, Russia
[2] Moscow MV Lomonosov State Univ, Computat Math & Cybernet Dept, Moscow 119991, Russia
[3] Tech Univ, Moscow Phys Engn Inst, Dept Theoret & Expt Phys, Moscow 115409, Russia
[4] Moscow MV Lomonosov State Univ, Dept Biol, Moscow 119991, Russia
关键词
Non-linear least-squares minimization problem; Normalized flux precision function; Partial optimization of experimental design; Convergence control of flux confidence interval bounds; METABOLIC FLUX ANALYSIS; BIDIRECTIONAL REACTION STEPS; CORYNEBACTERIUM-GLUTAMICUM; ISOTOPOMER DISTRIBUTIONS; LYSINE PRODUCTION; NETWORK ANALYSIS; UNITS EMU; DESIGN; QUANTIFICATION; SYSTEMS;
D O I
10.1186/s12934-014-0152-x
中图分类号
Q81 [生物工程学(生物技术)]; Q93 [微生物学];
学科分类号
071005 ; 0836 ; 090102 ; 100705 ;
摘要
Background: Steady-state C-13-based metabolic flux analysis (C-13-MFA) is the most powerful method available for the quantification of intracellular fluxes. These analyses include concertedly linked experimental and computational stages: (i) assuming the metabolic model and optimizing the experimental design; (ii) feeding the investigated organism using a chosen C-13-labeled substrate (tracer); (iii) measuring the extracellular effluxes and detecting the C-13-patterns of intracellular metabolites; and (iv) computing flux parameters that minimize the differences between observed and simulated measurements, followed by evaluating flux statistics. In its early stages, C-13-MFA was performed on the basis of data obtained in a single labeling experiment (SLE) followed by exploiting the developed high-performance computational software. Recently, the advantages of parallel labeling experiments (PLEs), where several LEs are conducted under the conditions differing only by the tracer(s) choice, were demonstrated, particularly with regard to improving flux precision due to the synergy of complementary information. The availability of an open-source software adjusted for PLE-based C-13-MFA is an important factor for PLE implementation. Results: The open-source software OpenFLUX, initially developed for the analysis of SLEs, was extended for the computation of PLE data. Using the OpenFLUX2, in silico simulation confirmed that flux precision is improved when C-13-MFA is implemented by fitting PLE data to the common model compared with SLE-based analysis. Efficient flux resolution could be achieved in the PLE-mediated analysis when the choice of tracer was based on an experimental design computed to minimize the flux variances from different parts of the metabolic network. The analysis provided by OpenFLUX2 mainly includes (i) the optimization of the experimental design, (ii) the computation of the flux parameters from LEs data, (iii) goodness-of-fit testing of the model's adequacy, (iv) drawing conclusions concerning the identifiability of fluxes and construction of a contribution matrix reflecting the relative contribution of the measurement variances to the flux variances, and (v) precise determination of flux confidence intervals using a fine-tunable and convergence-controlled Monte Carlo-based method. Conclusions: The developed open-source OpenFLUX2 provides a friendly software environment that facilitates beginners and existing OpenFLUX users to implement LEs for steady-state C-13-MFA including experimental design, quantitative evaluation of flux parameters and statistics.
引用
收藏
页数:25
相关论文
共 76 条
[41]   OpenFLUX: efficient modelling software for 13C-based metabolic flux analysis [J].
Quek, Lake-Ee ;
Wittmann, Christoph ;
Nielsen, Lars K. ;
Kromer, Jens O. .
MICROBIAL CELL FACTORIES, 2009, 8
[42]   Planning optimal measurements of isotopomer distributions for estimation of metabolic fluxes [J].
Rantanen, A ;
Mielikäinen, T ;
Rousu, J ;
Maaheimo, H ;
Ukkonen, E .
BIOINFORMATICS, 2006, 22 (10) :1198-1206
[43]   Construction of an E. Coli Genome-Scale Atom Mapping Model for MFA Calculations [J].
Ravikirthi, Prabhasa ;
Suthers, Patrick F. ;
Maranas, Costas D. .
BIOTECHNOLOGY AND BIOENGINEERING, 2011, 108 (06) :1372-1382
[44]  
Sauer U, 2001, METAB ENG, V3, P195
[45]  
Schellenberger Jan, 2012, BMC Syst Biol, V6, P9, DOI 10.1186/1752-0509-6-9
[46]  
Schmidt K, 1997, BIOTECHNOL BIOENG, V55, P831, DOI 10.1002/(SICI)1097-0290(19970920)55:6<831::AID-BIT2>3.0.CO
[47]  
2-H
[48]   Identification of glyA (encoding serine hydroxymethyltransferase) and its use together with the exporter ThrE to increase L-threonine accumulation by Corynebacterium glutamicum [J].
Simic, P ;
Willuhn, J ;
Sahm, H ;
Eggeling, L .
APPLIED AND ENVIRONMENTAL MICROBIOLOGY, 2002, 68 (07) :3321-3327
[49]   influx_s: increasing numerical stability and precision for metabolic flux analysis in isotope labelling experiments [J].
Sokol, Serguei ;
Millard, Pierre ;
Portais, Jean-Charles .
BIOINFORMATICS, 2012, 28 (05) :687-693
[50]   FLUX PARTITIONING IN THE SPLIT PATHWAY OF LYSINE SYNTHESIS IN CORYNEBACTERIUM-GLUTAMICUM QUANTIFICATION BY C-13-NMR AND H-1-NMR SPECTROSCOPY [J].
SONNTAG, K ;
EGGELING, L ;
DEGRAAF, AA ;
SAHM, H .
EUROPEAN JOURNAL OF BIOCHEMISTRY, 1993, 213 (03) :1325-1331