Dynamic sensitivity analysis of biological systems

BackgroundA mathematical model to understand, predict, control, or even design a real biological system is a central theme in systems biology. A dynamic biological system is always modeled as a nonlinear ordinary differential equation (ODE) system. How to simulate the dynamic behavior and dynamic parameter sensitivities of systems described by ODEs efficiently and accurately is a critical job. In many practical applications, e.g., the fed-batch fermentation systems, the system admissible input (corresponding to independent variables of the system) can be time-dependent. The main difficulty for investigating the dynamic log gains of these systems is the infinite dimension due to the time-dependent input. The classical dynamic sensitivity analysis does not take into account this case for the dynamic log gains.ResultsWe present an algorithm with an adaptive step size control that can be used for computing the solution and dynamic sensitivities of an autonomous ODE system simultaneously. Although our algorithm is one of the decouple direct methods in computing dynamic sensitivities of an ODE system, the step size determined by model equations can be used on the computations of the time profile and dynamic sensitivities with moderate accuracy even when sensitivity equations are more stiff than model equations. To show this algorithm can perform the dynamic sensitivity analysis on very stiff ODE systems with moderate accuracy, it is implemented and applied to two sets of chemical reactions: pyrolysis of ethane and oxidation of formaldehyde. The accuracy of this algorithm is demonstrated by comparing the dynamic parameter sensitivities obtained from this new algorithm and from the direct method with Rosenbrock stiff integrator based on the indirect method. The same dynamic sensitivity analysis was performed on an ethanol fed-batch fermentation system with a time-varying feed rate to evaluate the applicability of the algorithm to realistic models with time-dependent admissible input.ConclusionBy combining the accuracy we show with the efficiency of being a decouple direct method, our algorithm is an excellent method for computing dynamic parameter sensitivities in stiff problems. We extend the scope of classical dynamic sensitivity analysis to the investigation of dynamic log gains of models with time-dependent admissible input.

[1]  W. V. Loscutoff,et al.  General sensitivity theory , 1972 .

[2]  Igor Lutoshkin,et al.  The parameterization method in optimal control problems and differential-algebraic equations , 2006 .

[3]  W. E. Stewart,et al.  Sensitivity analysis of initial value problems with mixed odes and algebraic equations , 1985 .

[4]  Feng-Sheng Wang,et al.  A modified collocation method for solving differential-algebraic equations , 2000, Appl. Math. Comput..

[5]  P. Albrecht.  The Runge-Kutta Theory in a Nutshell , 1996 .

[6]  Guang-Ren Duan,et al.  Control parameterization enhancing transform for optimal control of switched systems , 2006, Math. Comput. Model..

[7]  Robert J. Gelinas,et al.  Sensitivity analysis of ordinary differential equation systems—A direct method , 1976 .

[8]  Peter Albrecht,et al.  A new theoretical approach to Runge-Kutta methods , 1987 .

[9]  H. Kacser,et al.  The control of flux. , 1995, Biochemical Society transactions.

[10]  R. Heinrich,et al.  The Regulation of Cellular Systems , 1996, Springer US.

[11]  Eberhard O. Voit,et al.  Computational Analysis of Biochemical Systems: A Practical Guide for Biochemists and Molecular Biologists , 2000 .

[12]  Herbert M Sauro,et al.  Sensitivity analysis of stoichiometric networks: an extension of metabolic control analysis to non-steady state trajectories. , 2003, Journal of theoretical biology.

[13]  K. F. Tipton,et al.  Biochemical systems analysis: A study of function and design in molecular biology , 1978 .

[14]  Eugene P. Dougherty,et al.  Further developments and applications of the Green’s function method of sensitivity analysis in chemical kinetics , 1979 .

[15]  Fumihide Shiraishi,et al.  A Taylor-series solution in Cartesian space to GMA-system equations and its application to initial-value problems , 2002, Appl. Math. Comput..

[16]  Yuji Hatoh,et al.  An efficient method for calculation of dynamic logarithmic gains in biochemical systems theory. , 2005 .

[17]  Curtis F. Gerald.  Applied numerical analysis , 1970 .

[18]  M. Savageau.  Biochemical systems analysis. II. The steady-state solutions for an n-pool system using a power-law approximation. , 1969, Journal of theoretical biology.

[19]  M. Savageau,et al.  Parameter Sensitivity as a Criterion for Evaluating and Comparing the Performance of Biochemical Systems , 1971, Nature.

[20]  Mark A. Kramer,et al.  The simultaneous solution and sensitivity analysis of systems described by ordinary differential equations , 1988, TOMS.

[21]  J. V. Villadsen.  Solution of boundary-value problems by orthogonal collocation: J. V. Villadsen and W. E. Stewart, Chem. Engng Sci.22, 1483–1501, 1967 , 1995 .

[22]  Herschel Rabitz,et al.  The Green’s function method of sensitivity analysis in chemical kinetics , 1978 .

[23]  A. Ghosh,et al.  Automated sensitivity analysis of stiff biochemical systems using a fourth-order adaptive step size Rosenbrock integration method. , 2006, Systems biology.

[24]  A. M. Dunker.  The decoupled direct method for calculating sensitivity coefficients in chemical kinetics , 1984 .

[25]  Feng-Sheng Wang,et al.  Hybrid Differential Evolution for Problems of Kinetic Parameter Estimation and Dynamic Optimization of an Ethanol Fermentation Process , 2001 .

[26]  E. Kreyszig.  Introductory Functional Analysis With Applications , 1978 .

[27]  V. K. Gorbunov.  A method for the parametrization of optimal control problems , 1979 .

[28]  Matthias Reuss,et al.  Dynamic sensitivity analysis for metabolic systems , 1997 .

[29]  Andres Kriete,et al.  Integrated spatio-temporal model of cell signaling , 2005 .

[30]  D. Fell.  Metabolic control analysis: a survey of its theoretical and experimental development. , 1992, The Biochemical journal.

[31]  W. E. Stewart,et al.  Solution of boundary-value problems by orthogonal collocation , 1995 .

[32]  Hasan Bulut,et al.  A Comparison of Numerical ODE Solvers based on Euler Methods , 1998 .

[33]  Yuji Hatoh,et al.  Erratum: An efficient method for calculation of dynamic logarithmic gains in biochemical systems theory (Journal of Theoretical Biology (2005) 234 (79-85) DOI: 10.1016/j.jtbi.2004.11.015) , 2006 .

[34]  Ji Shen,et al.  A DIRECT METHOD OF CALCULATING SENSITIVITY COEFFICIENTS OF CHEMICAL KINETICS , 1999 .

[35]  K. Teo,et al.  The control parameterization enhancing transform for constrained optimal control problems , 1999, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.