A reliable Taylor series-based computational method for the calculation of dynamic sensitivities in large-scale metabolic reaction systems: algorithm and software evaluation. (English) Zbl 1178.92018
Summary: Dynamic sensitivity analysis has become an important tool to successfully characterize all sorts of biological systems. However, when the analysis is carried out on large scale systems, it becomes imperative to employ a highly accurate computational method in order to obtain reliable values. Furthermore, the preliminary laborious mathematical operations required by current software before the computation of dynamic sensitivities makes it inconvenient for a significant number of unacquainted users. To satisfy these needs, the present work investigates a newly developed algorithm consisting of a combination of Taylor series method that can directly execute Taylor expansions for simultaneous nonlinear differential equations and a simple but highly-accurate numerical differentiation method based on finite-difference formulas.
Applications to three examples of biochemical systems indicate that the proposed method makes it possible to compute the dynamic sensitivity values with highly-reliable accuracies and also allows to readily compute them by setting up only the differential equations for metabolite concentrations in the computer program. Also, it is found that the Padé approximation introduced in the Taylor series method shortens the computation time greatly because it stabilizes the computation so that it allows us to use larger step sizes in the numerical integration. Consequently, the calculated results suggest that the proposed computational method, in addition to being user-friendly, makes it possible to perform dynamic sensitivity analysis in large-scale metabolic reaction systems both efficiently and reliably.
Applications to three examples of biochemical systems indicate that the proposed method makes it possible to compute the dynamic sensitivity values with highly-reliable accuracies and also allows to readily compute them by setting up only the differential equations for metabolite concentrations in the computer program. Also, it is found that the Padé approximation introduced in the Taylor series method shortens the computation time greatly because it stabilizes the computation so that it allows us to use larger step sizes in the numerical integration. Consequently, the calculated results suggest that the proposed computational method, in addition to being user-friendly, makes it possible to perform dynamic sensitivity analysis in large-scale metabolic reaction systems both efficiently and reliably.
MSC:
| 92C40 | Biochemistry, molecular biology |
| 65D25 | Numerical differentiation |
| 41A21 | Padé approximation |
| 65D30 | Numerical integration |
Keywords:
dynamic sensitivity; high accuracy and reliability; reaction engineering; simulation; non-linear dynamics; metabolismSoftware:
ESSYNSReferences:
| [1] | Varma, A.; Morbidelli, M.; Wu, H., Parametric Sensitivity in Chemical Systems (1999), Cambridge University: Cambridge University Cambridge |
| [2] | Rabitz, H.; Kramer, M.; Dacol, D., Sensitivity analysis in chemical kinetics, Annu. Rev. Phys. Chem., 34, 419 (1983) |
| [3] | Seigneur, C.; Stephanopoulos, G.; Carr, R. W., Dynamic sensitivity analysis of chemical reaction systems, Chem. Eng. Sci., 37, 845 (1982) |
| [4] | Shiraishi, F.; Hatoh, Y., Dynamic sensitivities in chaotic dynamical systems, Appl. Math. Comput., 186, 1347 (2007) · Zbl 1117.65167 |
| [5] | Ingalls, B. P.; Sauro, H. M., Sensitivity analysis of stoichiometric networks: an extension of metabolic control analysis to non-steady state trajectories, J. Theor. Biol., 222, 23 (2003) · Zbl 1464.92110 |
| [6] | Schwacke, J. H.; Voit, E. O., Computation and analysis of time-dependent sensitivities in generalized mass action systems, J. Theor. Biol., 236, 21 (2005) · Zbl 1442.92061 |
| [7] | Shiraishi, F.; Hatoh, Y.; Irie, T., An efficient method for calculation of dynamic logarithmic gains in biochemical systems theory, J. Theor. Biol., 234, 79 (2005) · Zbl 1445.92126 |
| [8] | Mauch, K.; Arnold, S.; Reuss, M., Dynamic sensitivity analysis for metabolic systems, Chem. Eng. Sci., 52, 2589 (1997) |
| [9] | Wu, S. H.; Wang, F. S.; Chang, M. S., Dynamic sensitivity analysis of biological systems, BMC Bioinf., 9, S1 (2008) |
| [10] | Irvine, D. H.; Savageau, M. A., Efficient solution of nonlinear ordinary differential equations expressed in S-system canonical form, SIAM J. Numer. Anal., 27, 704 (1990) · Zbl 0702.65068 |
| [11] | Savageau, M. A.; Voit, E. O., Recasting nonlinear differential equations as S-systems: a canonical nonlinear form, Math. Biosci., 87, 83 (1987) · Zbl 0631.34015 |
| [12] | Shiraishi, F.; Takeuchi, H.; Hasegawa, T.; Nagasue, H., A Taylor-series solution in Cartesian space to GMA-system equations and its application to initial-value problems, Appl. Math. Comput., 127, 103 (2002) · Zbl 1023.65064 |
| [13] | Shiraishi, F., Highly accurate solution of the axial dispersion model expressed in S-system canonical form by Taylor-series method, Chem. Eng. J., 83, 175 (2001) |
| [14] | Shiraishi, F.; Fujiwara, S., An efficient method for solving two-point boundary value problems with extremely high accuracy, J. Chem. Eng. Jpn., 29, 88 (1996) |
| [15] | Shiraishi, F.; Hasegawa, T.; Nagasue, H., Numerical solution of the two-point boundary value problem by the combined Taylor series method with a technique for rapidly selecting suitable stepsizes, J. Chem. Eng. Jpn., 28, 306 (1995) |
| [16] | Miyakawa, H.; Nagasue, H.; Shiraishi, F., A highly accurate numerical method for calculating apparent kinetic parameters of immobilized enzyme reactions: 1. Theory, Biochem. Eng. J., 3, 91 (1999) |
| [17] | Miyakawa, H.; Nagasue, H.; Shiraishi, F., A highly accurate numerical method for calculating apparent kinetic parameters of immobilized enzyme reactions: 2. Accuracies of calculated values, Biochem. Eng. J., 3, 103 (1999) |
| [18] | H. Nagasue, Dissertation, Kyushu University, 1994.; H. Nagasue, Dissertation, Kyushu University, 1994. |
| [19] | Ascher, U. M.; Petzold, L. R., Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations (1998), SIAM: SIAM Philadelphia, PA · Zbl 0908.65055 |
| [20] | Shiraishi, F.; Furuta, S.; Ishimatsu, T.; Akhter, J., A simple and highly accurate numerical differentiation method for sensitivity analysis of large-scale metabolic reaction systems, Math. Biosci., 208, 590 (2007) · Zbl 1119.92035 |
| [21] | Hirayama, H.; Komiya, S.; Sato, S., Solving ordinary differential equations by Taylor series, JSIAM, 12, 1 (2002) |
| [22] | Baker, G. A.; Graves-Morris, P., Padé Approximants (1996), Cambridge University: Cambridge University New York · Zbl 0923.41001 |
| [23] | Curto, R.; Sorribas, A.; Cascante, M., Comparative characterization of the fermentation pathway of Saccharomyces cerevisiae using biochemical systems theory and metabolic control analysis: model definition and nomenclature, Math. Biosci., 130, 25 (1995) · Zbl 0835.92015 |
| [24] | Albe, K. R.; Wright, B. E., Systems analysis of the tricarboxylic acid cycle in Dictyostelium discoideum. II. Control analysis, J. Biol. Chem., 267, 3106 (1992) |
| [25] | Wright, B. E.; Butler, M. H.; Albe, K. R., Systems analysis of the tricarboxylic acid cycle in Dictyostelium discoideum. I. The basis for model construction, J. Biol. Chem., 267, 3101 (1992) |
| [26] | Shiraishi, F.; Savageau, M. A., The tricarboxylic acid cycle in Dictyostelium discoideum. I. Formulation of alternative kinetic representations, J. Biol. Chem., 267, 22912 (1992) |
| [27] | Shiraishi, F.; Savageau, M. A., The tricarboxylic acid cycle in Dictyostelium discoideum. II. Evaluation of model consistency and robustness, J. Biol. Chem., 267, 22919 (1992) |
| [28] | Shiraishi, F.; Savageau, M. A., The tricarboxylic acid cycle in Dictyostelium discoideum. III. Analysis of steady-state and dynamic behavior, J. Biol. Chem., 267, 22926 (1992) |
| [29] | Shiraishi, F.; Savageau, M. A., The tricarboxylic acid cycle in Dictyostelium discoideum. IV. Resolution of discrepancies between alternative methods of analysis, J. Biol. Chem., 267, 22934 (1992) |
| [30] | Shiraishi, F.; Savageau, M. A., The tricarboxylic acid cycle in Dictyostelium discoideum: systemic effects of including protein turnover in the current model, J. Biol. Chem., 268, 16917 (1993) |
| [31] | Dunker, A. M., The decoupled direct method for calculating sensitivities coefficients in chemical kinetics, J. Chem. Phys., 81, 2385 (1984) |
| [32] | Kramer, M. S.; Calo, J. M.; Rabits, H., An improved computational method for sensitivity analysis: Green’s function method with AIM, Appl. Math. Model., 5, 432 (1981) · Zbl 0469.65049 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
