DFBAlab: a fast and reliable MATLAB code for dynamic flux balance analysis

BackgroundDynamic Flux Balance Analysis (DFBA) is a dynamic simulation framework for biochemical processes. DFBA can be performed using different approaches such as static optimization (SOA), dynamic optimization (DOA), and direct approaches (DA). Few existing simulators address the theoretical and practical challenges of nonunique exchange fluxes or infeasible linear programs (LPs). Both are common sources of failure and inefficiencies for these simulators.ResultsDFBAlab, a MATLAB-based simulator that uses the LP feasibility problem to obtain an extended system and lexicographic optimization to yield unique exchange fluxes, is presented. DFBAlab is able to simulate complex dynamic cultures with multiple species rapidly and reliably, including differential-algebraic equation (DAE) systems. In addition, DFBAlab’s running time scales linearly with the number of species models. Three examples are presented where the performance of COBRA, DyMMM and DFBAlab are compared.ConclusionsLexicographic optimization is used to determine unique exchange fluxes which are necessary for a well-defined dynamic system. DFBAlab does not fail during numerical integration due to infeasible LPs. The extended system obtained through the LP feasibility problem in DFBAlab provides a penalty function that can be used in optimization algorithms.

[1]  H. O. Buhr,et al.  A dynamic model of the high-rate algal-bacterial wastewater treatment pond , 1983 .

[2]  P I Barton,et al.  A reliable simulator for dynamic flux balance analysis , 2013, Biotechnology and bioengineering.

[3]  B. Palsson,et al.  An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR) , 2003, Genome Biology.

[4]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.

[5]  Paul I. Barton,et al.  Design of Microbial Consortia for Industrial Biotechnology , 2014 .

[6]  Timothy J. Hanly,et al.  Dynamic flux balance modeling of microbial co‐cultures for efficient batch fermentation of glucose and xylose mixtures , 2011, Biotechnology and bioengineering.

[7]  Aidong Yang,et al.  Modeling and Evaluation of CO2 Supply and Utilization in Algal Ponds , 2011 .

[8]  John N. Tsitsiklis,et al.  Introduction to linear optimization , 1997, Athena scientific optimization and computation series.

[9]  R. Mahadevan,et al.  The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. , 2003, Metabolic engineering.

[10]  Radhakrishnan Mahadevan,et al.  Genome-scale dynamic modeling of the competition between Rhodoferax and Geobacter in anoxic subsurface environments , 2011, The ISME Journal.

[11]  Lake-Ee Quek,et al.  AlgaGEM – a genome-scale metabolic reconstruction of algae based on the Chlamydomonas reinhardtii genome , 2011, BMC Genomics.

[12]  Jeffrey D Orth,et al.  What is flux balance analysis? , 2010, Nature Biotechnology.

[13]  Markus J. Herrgård,et al.  Reconstruction and validation of Saccharomyces cerevisiae iND750, a fully compartmentalized genome-scale metabolic model. , 2004, Genome research.

[14]  Jason A. Papin,et al.  Metabolic network reconstruction of Chlamydomonas offers insight into light-driven algal metabolism , 2011, Molecular systems biology.

[15]  Radhakrishnan Mahadevan,et al.  The design of long‐term effective uranium bioremediation strategy using a community metabolic model , 2012, Biotechnology and bioengineering.

[16]  Michael R. Johns,et al.  Kinetic models for heterotrophic growth of Chlamydomonas reinhardtii in batch and fed-batch cultures , 1999 .

[17]  Wynand S. Verwoerd,et al.  ORCA: a COBRA toolbox extension for model-driven discovery and analysis , 2014, Bioinform..

[18]  F. Doyle,et al.  Dynamic flux balance analysis of diauxic growth in Escherichia coli. , 2002, Biophysical journal.

[19]  B. Colman,et al.  Mechanism of CO2 acquisition in an acid-tolerant Chlamydomonas. , 2007, Plant, cell & environment.