Dynamic estimation of specific fluxes in metabolic networks using non-linear dynamic optimization

BMC Systems Biology

Jan Van Impe¹  Filip Logist¹  Dominique Vercammen¹ 

[1] KU Leuven, BioTeC - Chemical and Biochemical Process Technology and Control & OPTEC - Center of Excellence: Optimization in Engineering, Department of Chemical Engineering, Willem de Croylaan 46/2423, Leuven 3001, Belgium

关键词: Parameter estimation;    Non-linear optimization;    B-spline parameterizations;    Dynamic metabolic flux analysis;   

Others  :  1091377
DOI  :  10.1186/s12918-014-0132-0

 received in 2014-03-31, accepted in 2014-08-13,  发布年份 2014

【 摘 要 】

Background

Metabolic network models describing the biochemical reaction network and material fluxes inside microorganisms open interesting routes for the model-based optimization of bioprocesses. Dynamic metabolic flux analysis (dMFA) has lately been studied as an extension of regular metabolic flux analysis (MFA), rendering a dynamic view of the fluxes, also in non-stationary conditions. Recent dMFA implementations suffer from some drawbacks, though. More specifically, the fluxes are not estimated as specific fluxes, which are more biologically relevant. Also, the flux profiles are not smooth, and additional constraints like, e.g., irreversibility constraints on the fluxes, cannot be taken into account. Finally, in all previous methods, a basis for the null space of the stoichiometric matrix, i.e., which set of free fluxes is used, needs to be chosen. This choice is not trivial, and has a large influence on the resulting estimates.

Results

In this work, a new methodology based on a B-spline parameterization of the fluxes is presented. Because of the high degree of non-linearity due to this parameterization, an incremental knot insertion strategy has been devised, resulting in a sequence of non-linear dynamic optimization problems. These are solved using state-of-the-art dynamic optimization methods and tools, i.e., orthogonal collocation, an interior-point optimizer and automatic differentiation. Also, a procedure to choose an optimal basis for the null space of the stoichiometric matrix is described, discarding the need to make a choice beforehand. The proposed methodology is validated on two simulated case studies: (i) a small-scale network with 7 fluxes, to illustrate the operation of the algorithm, and (ii) a medium-scale network with 68 fluxes, to show the algorithm’s capabilities for a realistic network. The results show an accurate correspondence to the reference fluxes used to simulate the measurements, both in a theoretically ideal setting with no experimental noise, and in a realistic noise setting.

Conclusions

Because, apart from a metabolic reaction network and the measurements, no extra input needs to be given, the resulting algorithm is a systematic, integrated and accurate methodology for dynamic metabolic flux analysis that can be run online in real-time if necessary.

【 授权许可】

   
2014 Vercammen et al.; licensee BioMed Central.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

Figure 6.

Figure 7.

Figure 8.

Figure 9.

Figure 10.

Figure 11.

Figure 12.

Figure 13.

Figure 14.

Figure 15.

Figure 16.

Figure 17.

Figure 18.

【 参考文献 】

• [1]Carinhas N, Bernal V, Teixeira AP, Carrondo MJT, Alves PM, Oliveira R: Hybrid metabolic flux analysis: combining stoichiometric and statistical constraints to model the formation of complex recombinant products. BMC Syst Biol 2011, 5:34. BioMed Central Full Text
• [2]Noor E, Lewis NE, Milo R: A proof for loop-law constraints in stoichiometric metabolic networks. BMC Syst Biol 2012, 6:140. BioMed Central Full Text
• [3]Zamorano F, Vande Wouwer A, Bastin G: A detailed metabolic flux analysis of an underdetermined network of CHO cells. J Biotechnol 2010, 150(4):497-508.
• [4]Fischer E, Zamboni N, Sauer U: High-throughput metabolic flux analysis based on gas chromatography–mass spectrometry derived 13C constraints. Biotechnol Bioeng 2001, 76(2):144-56.
• [5]Leighty RW, Antoniewicz MR: Dynamic metabolic flux analysis (DMFA): a framework for determining fluxes at metabolic non-steady state. Metab Eng 2011, 13(6):745-55.
• [6]Van Impe JF, Vercammen D, Van Derlinden E: Toward a next generation of predictive models: a systems biology primer. Food Control 2013, 29(2):336-42.
• [7]Llaneras F, Picȯ J: Stoichiometric modelling of cell metabolism. J Biosci Bioeng 2008, 105(1):1-11.
• [8]Stephanopoulos GN, Aristidou AA, Nielsen J: Metabolic Engineering: Principles and Methodologies . Academic Press, San Diego, USA; 1998.
• [9]Antoniewicz MR: Dynamic metabolic flux analysis - tools for probing transient states of metabolic networks. Curr Opin Biotech 2013, 24:973-8.
• [10]Lequeux G, Beauprez J, Maertens J, Van Horen E, Soetaert W, Vandamme E, Vanrolleghem PA: Dynamic metabolic flux analysis demonstrated on cultures where the limiting substrate is changed from carbon to nitrogen and vice versa. J Biomed Biotechnol 2010, 621645:1-19.
• [11]Niklas J, Schrȧder E, Sandig V, Noll T, Heinzle E: Quantitative characterization of metabolism and metabolic shifts during growth of the new human cell line AGE1.HN using time resolved metabolic flux analysis. Bioproc Biosyst Eng 2011, 34(5):533-45.
• [12]Llaneras F, Sala A, Picȯ J: Dynamic estimations of metabolic fluxes with constraint-based models and possibility theory. J Process Contr 2012, 22(10):1946-55.
• [13]Chou I, Voit EO: Estimation of dynamic flux profiles from metabolic time series data. BMC Syst Biol 2012, 6:84. BioMed Central Full Text
• [14]Zhan C, Yeung LF: Parameter estimation in systems biology models using spline approximation. BMC Syst Biol 2011, 5:14. BioMed Central Full Text
• [15]Dierckx P: Curve and Surface Fitting with Splines. New York . Oxford University Press, USA; 1993.
• [16]Van Loock W, Pipeleers G, De Schutter J, Swevers J: A convex optimization approach to curve fitting with, B-splines. In Proceedings of the 18th IFAC World Congress: 28 August - 2 September 2011; Milano, Italy . Edited by Bittanti S, Cenedese A, Zampieri S. International Federation of Automatic Control, Milano, Italy; 2011:2290-5.
• [17]de Boor C: A Practical Guide to Splines . Springer, New York, USA; 2001.
• [18]Biegler LT: Nonlinear Programming: Concepts, Algorithms, and Applications to Chemical Processes . MPS-SIAM, Philadelphia, USA; 2010.
• [19]Biegler LT: Solution of dynamic optimization problems by successive quadratic programming and orthogonal collocation. Comput Chem Eng 1984, 8:243-248.
• [20]Wächter A, Biegler LT: On the implementation of a primal-dual interior point filter line search algorithm for large-scale nonlinear programming. Math Program 2006, 106(1):25-57.
• [21]Andersson J, Åkesson J, Diehl M: CasADi – A symbolic package for automatic differentiation and optimal control. In Recent Advances in Algorithmic Differentiation. Lecture Notes in Computational Science and Engineering. Volume 87 . Edited by Forth S, Hovland P, Phipps E, Utke J, Walther A. Springer, Berlin Heidelberg; 2012:297-307.
• [22]Burnham KP, Anderson DR: Multimodel inference: understanding AIC and BIC in model selection. Social Method Res 2004, 33:261-304.
• [23]Motulsky H, Christopoulos A: Fitting Models to Biological Data Using Linear and Nonlinear Regression: A Practical Guide to Curve Fitting . Oxford University Press, New York, USA; 2004.
• [24]Snedecor GW, Cochran WG: Statistical Methods, Eighth Edition . Iowa State University Press, Ames, USA; 1989.
• [25]Efron B, Tibshirani R: An Introduction to the Bootstrap . Chapman and Hall/CRC, London, UK; 1994.
• [26]Schellenberger J, Zielinski DC, Choi W, Madireddi S, Portnoy V, Scott DA, Reed JL, Osterman AL, Palsson BØ: Predicting outcomes of steady-state13C isotope tracing experiments using Monte Carlo sampling. BMC Syst Biol 2012, 6:9. BioMed Central Full Text