A simple and highly accurate numerical differentiation method for sensitivity analysis of large-scale metabolic reaction systems
Related Experiment Video
Updated: Jul 16, 2026

The Use of Chemostats in Microbial Systems Biology
Published on: October 14, 2013
Fumihide Shiraishi1, Shingo Furuta, Takaaki Ishimatsu
1Department of Bio-System Design, Bio-Architecture Center, Kyushu University, 3-1-1, Maidashi, Fukuoka 812-8582, Japan. fumishira@brs.kyushu-u.ac.jp
Abstract:
Numerical differentiation is known to be one of the most difficult numerical calculation methods to obtain reliable calculated values at all times. A simple numerical differentiation method using a combination of finite-difference formulas, derived by approximation of Taylor-series equations, is investigated in order to efficiently perform the sensitivity analysis of large-scale metabolic reaction systems. A result of the application to four basic mathematical functions reveals that the use of the eight-point differentiation formula with a non-dimensionalized stepsize close to 0.01 mostly provides more than 14 digits of accuracy in double precision for the numerical derivatives. Moreover, a result of the application to the modified TCA cycle model indicates that the numerical differentiation method gives the calculated values of steady-state metabolite concentrations within a range of round-off error and also makes it possible to transform the Michaelis-Menten equations into the S-system equations having the kinetic orders whose accuracies are mostly more than 14 significant digits. Because of the simple structure of the numerical differentiation formula and its promising high accuracy, it is evident that the present numerical differentiation method is useful for the analysis of large-scale metabolic reaction systems according to the systematic procedure of BST.
Related Concept Videos
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
Reaction Mechanisms: The Steady-State Approximation
Linear Differential Equations
Chemical Equilibria: Systematic Approach to Equilibrium Calculations
The first step is to identify all the chemical reactions involved, The...
Reaction Mechanisms: Rate-limiting Step Approximation
Determination of Multiple Dosing Parameters: Steady-State, Minimum and Maximum Concentrations

