Pathway identification using parallel optimization for a nonlinear hybrid system in batch culture

Abstract The complex bio-process for the bioconversion of glycerol to 1,3-propanediol (1,3-PD) can be characterized by a nonlinear hybrid dynamical system of enzyme-catalytic kinetics and genetic regulation. In this paper, in consideration of the possible ways that 3-hydroxypropionaldehyde (3-HPA) inhibits the cell growth, various possible transports and inhibition mechanisms, we consider 576 possible metabolic pathways and establish a fourteen-dimensional nonlinear hybrid dynamical system with uncertain system parameters and pathway parameters for describing the process of batch culture. Some important properties of the hybrid system are discussed. Taking into account the difficulty in accurately measuring the concentration of intracellular substances and the absence of equilibrium points for the hybrid system, we quantitatively define biological robustness of the intracellular substance concentrations for the overall process of batch culture. Our goal is to determine the most possible metabolic pathway and corresponding system parameters. To this end, taking the relative error between the experimental data and the computational values of the extracellular substances together with the proposed biological robustness of the intracellular substances as a cost function, we establish an identification model, in which 17 280 system parameters and 5760 pathway parameters are involved. Furthermore, the identification model is subject to the hybrid system, continuous state constraints and parameter constraints. As such, it is a very complicated task to solve the identification model by a serial program. With this in mind, we propose a parallel improved differential evolution algorithm (P-IDE), based on the constraint transcription and smoothing approximation techniques, for solving the identification model. An illustrative numerical example shows the appropriateness of the most possible metabolic pathway and the corresponding system parameters as well as the effectiveness of the parallel algorithm.

[1]  A. Zeng,et al.  Microbial production of 1,3-propanediol , 1999, Applied Microbiology and Biotechnology.

[2]  Kok Lay Teo,et al.  A Unified Computational Approach to Optimal Control Problems , 1991 .

[3]  Hiroaki Kitano,et al.  Biological robustness , 2008, Nature Reviews Genetics.

[4]  Enmin Feng,et al.  An improved model for multistage simulation of glycerol fermentation in batch culture and its parameter identification , 2009 .

[5]  Elijah Polak,et al.  Optimization: Algorithms and Consistent Approximations , 1997 .

[6]  Jibin Sun,et al.  Comparative Genomic Analysis of dha Regulon and Related Genes for Anaerobic Glycerol Metabolism in Bacteria , 2003, Biotechnology progress.

[7]  Chongyang Liu,et al.  From the SelectedWorks of Chongyang Liu 2013 Modelling and parameter identification for a nonlinear time-delay system in microbial batch fermentation , 2017 .

[8]  A. Zeng,et al.  Multiple product inhibition and growth modeling of clostridium butyricum and klebsiella pneumoniae in glycerol fermentation , 1994, Biotechnology and bioengineering.

[9]  José Carlos Pinto,et al.  Sequential experimental design for model discrimination: Taking into account the posterior covariance matrix of differences between model predictions , 2008 .

[10]  E. Biscaia,et al.  Nonlinear parameter estimation through particle swarm optimization , 2008 .

[11]  B. V. Babu,et al.  Modified differential evolution (MDE) for optimization of non-linear chemical processes , 2006, Comput. Chem. Eng..

[12]  Jiang-She Zhang,et al.  Global optimization by an improved differential evolutionary algorithm , 2007, Appl. Math. Comput..

[13]  Enmin Feng,et al.  Modelling and pathway identification involving the transport mechanism of a complex metabolic system in batch culture , 2014, Commun. Nonlinear Sci. Numer. Simul..

[14]  A. Zeng,et al.  A kinetic model for product formation of microbial and mammalian cells , 1995, Biotechnology and bioengineering.

[15]  Matjaž Perc,et al.  Noise enhances robustness of intracellular Ca2+ oscillations , 2003 .

[16]  Enmin Feng,et al.  Modeling and identification of a nonlinear hybrid dynamical system in batch fermentation of glycerol , 2011, Math. Comput. Model..

[17]  C. Nakamura,et al.  Metabolic engineering for the microbial production of 1,3-propanediol. , 2003, Current opinion in biotechnology.

[18]  H. Kitano.  Towards a theory of biological robustness , 2007, Molecular systems biology.

[19]  Ville Tirronen,et al.  Shuffle or update parallel differential evolution for large-scale optimization , 2011, Soft Comput..

[20]  Fabrice Heitz,et al.  Parallel Differential Evolution: Application to 3-D Medical Image Registration , 2005 .

[21]  Enmin Feng,et al.  Modeling nonlinear stochastic kinetic system and stochastic optimal control of microbial bioconversion process in batch culture , 2013 .

[22]  Heitor S. Lopes,et al.  An improved parallel differential evolution approach for protein structure prediction using both 2D and 3D off-lattice models , 2013, 2013 IEEE Symposium on Differential Evolution (SDE).

[23]  Enmin Feng,et al.  Nonlinear Mathematical Simulation and Analysis of Dha Regulon for Glycerol Metabolism in Klebsiella pneumoniae , 2012 .

[24]  Panagiotis D. Christofides,et al.  Stabilization of nonlinear systems with state and control constraints using Lyapunov-based predictive control , 2005, Proceedings of the 2005, American Control Conference, 2005..

[25]  Wolf-Dieter Deckwer,et al.  Glycerol conversion to 1,3-propanediol by newly isolated clostridia , 1992, Applied Microbiology and Biotechnology.

[26]  Dimitris K. Tasoulis,et al.  Parallel differential evolution , 2004, Proceedings of the 2004 Congress on Evolutionary Computation (IEEE Cat. No.04TH8753).

[27]  Amaro G. Barreto,et al.  A new approach for sequential experimental design for model discrimination , 2006 .

[28]  Enmin Feng,et al.  Complex metabolic network of glycerol fermentation by Klebsiella pneumoniae and its system identification via biological robustness , 2011 .

[29]  Enmin Feng,et al.  Pathway identification using parallel optimization for a complex metabolic system in microbial continuous culture , 2011 .

[30]  Enmin Feng,et al.  Identification and robustness analysis of nonlinear multi-stage enzyme-catalytic dynamical system in batch culture , 2015 .

[31]  Enmin Feng,et al.  Optimal Control for Multistage Nonlinear Dynamic System of Microbial Bioconversion in Batch Culture , 2011, J. Appl. Math..

[32]  Panagiotis D. Christofides,et al.  Predictive control of switched nonlinear systems with scheduled mode transitions , 2005, IEEE Transactions on Automatic Control.

[33]  Enmin Feng,et al.  Robust identification of enzymatic nonlinear dynamical systems for 1, 3-propanediol transport mechanisms in microbial batch culture , 2014, Appl. Math. Comput..

[34]  P. N. Suganthan,et al.  Differential Evolution Algorithm With Strategy Adaptation for Global Numerical Optimization , 2009, IEEE Transactions on Evolutionary Computation.

[35]  Martin Buss,et al.  Nonlinear Hybrid Dynamical Systems: Modeling, Optimal Control, and Applications , 2002 .

[36]  Zhijian Wu,et al.  Parallel differential evolution with self-adapting control parameters and generalized opposition-based learning for solving high-dimensional optimization problems , 2013, J. Parallel Distributed Comput..

[37]  Rainer Storn,et al.  Differential Evolution – A Simple and Efficient Heuristic for global Optimization over Continuous Spaces , 1997, J. Glob. Optim..

[38]  Matjaz Perc,et al.  Sensitivity and flexibility of regular and chaotic calcium oscillations. , 2003, Biophysical chemistry.

[39]  Lei Wang,et al.  Determining the transport mechanism of an enzyme-catalytic complex metabolic network based on biological robustness , 2013, Bioprocess and Biosystems Engineering.

[40]  J. Stelling,et al.  Robustness of Cellular Functions , 2004, Cell.

[41]  Enmin Feng,et al.  Sensitivity analysis and identification of kinetic parameters in batch fermentation of glycerol , 2012, J. Comput. Appl. Math..

[42]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[43]  D. C. Cameron,et al.  Enhancement of 1,3-Propanediol production by cofermentation inEscherichia coli expressingKlebsiella pneumoniae dha regulon genes , 1992, Applied biochemistry and biotechnology.

[44]  D. C. Cameron,et al.  1,3-Propanediol production by Escherichia coli expressing genes from the Klebsiella pneumoniae dha regulon , 1991, Applied and environmental microbiology.

[45]  Enmin Feng,et al.  Robust identification and its properties of nonlinear bilevel multi-stage dynamic system , 2013, Appl. Math. Comput..

[46]  Weihang Zhu,et al.  Massively parallel differential evolution—pattern search optimization with graphics hardware acceleration: an investigation on bound constrained optimization problems , 2011, J. Glob. Optim..

[47]  Wolf-Dieter Deckwer,et al.  Synthesis, properties and biodegradability of polyesters based on 1,3‐propanediol , 1994 .

[48]  Zhaohua Gong,et al.  Modelling and optimal control of a time-delayed switched system in fed-batch process , 2014, J. Frankl. Inst..

[49]  Enmin Feng,et al.  Identification and robustness analysis of nonlinear hybrid dynamical system concerning glycerol transport mechanism , 2012, Comput. Chem. Eng..

[50]  S. Leibler,et al.  Robustness in simple biochemical networks , 1997, Nature.

[51]  Enmin Feng,et al.  Joint estimation in batch culture by using unscented kalman filter , 2012, Biotechnology and Bioprocess Engineering.

[52]  Zhilong Xiu,et al.  Mathematical modeling of glycerol fermentation by Klebsiella pneumoniae: Concerning enzyme-catalytic reductive pathway and transport of glycerol and 1,3-propanediol across cell membrane , 2008 .

[53]  Behzad Moshiri,et al.  Optimal control of a nonlinear fed-batch fermentation process using model predictive approach , 2009 .