Highly accurate computation of dynamic sensitivities in metabolic reaction systems by a Taylor series method

Fumihide Shiraishi, Masaaki Egashira, Michio Iwata

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

We have previously developed the software for calculation of dynamic sensitivities, SoftCADS, in which one can calculate dynamic sensitivities with high accuracy by just setting the differential equations for metabolite concentrations. However, SoftCADS did not always provide calculated values with the machine accuracy of a computer, although a Taylor series method was employed to numerically solve the differential equations. This is because numerical derivatives calculated from an approximate formula were directly used in the derivation of the differential equations for sensitivities from those for metabolite concentrations. The present work therefore attempts to further enhance the performance of SoftCADS, including not only the accuracies of the calculated values but also the calculation time. To overcome the problem, the approximate formula is expanded into a Taylor series in time and the first-term value of the series is replaced by the exact coefficient on the second term of the flux function expanded into a Taylor series in an independent or dependent variable. The result reveals that this replacement certainly provides not only numerical derivatives but also dynamic sensitivities with superhigh accuracies comparable to the machine accuracy, regardless of the degree of stiffness of the differential equations. Moreover, a comparison indicates that the improved SoftCADS shortens the calculation time of the dynamic sensitivities without reducing their accuracies, even when the simplest approximate derivative formula is used.

Original languageEnglish
Pages (from-to)59-67
Number of pages9
JournalMathematical Biosciences
Volume233
Issue number1
DOIs
Publication statusPublished - Sep 1 2011

Fingerprint

Taylor series
Differential equations
Numerical derivative
Differential equation
Metabolites
Derivatives
Software
methodology
metabolites
Term
Stiffness
Replacement
Fluxes
High Accuracy
time series analysis
Calculate
Derivative
Series
Dependent
Coefficient

All Science Journal Classification (ASJC) codes

  • Statistics and Probability
  • Modelling and Simulation
  • Biochemistry, Genetics and Molecular Biology(all)
  • Immunology and Microbiology(all)
  • Agricultural and Biological Sciences(all)
  • Applied Mathematics

Cite this

Highly accurate computation of dynamic sensitivities in metabolic reaction systems by a Taylor series method. / Shiraishi, Fumihide; Egashira, Masaaki; Iwata, Michio.

In: Mathematical Biosciences, Vol. 233, No. 1, 01.09.2011, p. 59-67.

Research output: Contribution to journalArticle

@article{988fe1d081194609aab17f1746740e68,
title = "Highly accurate computation of dynamic sensitivities in metabolic reaction systems by a Taylor series method",
abstract = "We have previously developed the software for calculation of dynamic sensitivities, SoftCADS, in which one can calculate dynamic sensitivities with high accuracy by just setting the differential equations for metabolite concentrations. However, SoftCADS did not always provide calculated values with the machine accuracy of a computer, although a Taylor series method was employed to numerically solve the differential equations. This is because numerical derivatives calculated from an approximate formula were directly used in the derivation of the differential equations for sensitivities from those for metabolite concentrations. The present work therefore attempts to further enhance the performance of SoftCADS, including not only the accuracies of the calculated values but also the calculation time. To overcome the problem, the approximate formula is expanded into a Taylor series in time and the first-term value of the series is replaced by the exact coefficient on the second term of the flux function expanded into a Taylor series in an independent or dependent variable. The result reveals that this replacement certainly provides not only numerical derivatives but also dynamic sensitivities with superhigh accuracies comparable to the machine accuracy, regardless of the degree of stiffness of the differential equations. Moreover, a comparison indicates that the improved SoftCADS shortens the calculation time of the dynamic sensitivities without reducing their accuracies, even when the simplest approximate derivative formula is used.",
author = "Fumihide Shiraishi and Masaaki Egashira and Michio Iwata",
year = "2011",
month = "9",
day = "1",
doi = "10.1016/j.mbs.2011.06.004",
language = "English",
volume = "233",
pages = "59--67",
journal = "Mathematical Biosciences",
issn = "0025-5564",
publisher = "Elsevier Inc.",
number = "1",

}

TY - JOUR

T1 - Highly accurate computation of dynamic sensitivities in metabolic reaction systems by a Taylor series method

AU - Shiraishi, Fumihide

AU - Egashira, Masaaki

AU - Iwata, Michio

PY - 2011/9/1

Y1 - 2011/9/1

N2 - We have previously developed the software for calculation of dynamic sensitivities, SoftCADS, in which one can calculate dynamic sensitivities with high accuracy by just setting the differential equations for metabolite concentrations. However, SoftCADS did not always provide calculated values with the machine accuracy of a computer, although a Taylor series method was employed to numerically solve the differential equations. This is because numerical derivatives calculated from an approximate formula were directly used in the derivation of the differential equations for sensitivities from those for metabolite concentrations. The present work therefore attempts to further enhance the performance of SoftCADS, including not only the accuracies of the calculated values but also the calculation time. To overcome the problem, the approximate formula is expanded into a Taylor series in time and the first-term value of the series is replaced by the exact coefficient on the second term of the flux function expanded into a Taylor series in an independent or dependent variable. The result reveals that this replacement certainly provides not only numerical derivatives but also dynamic sensitivities with superhigh accuracies comparable to the machine accuracy, regardless of the degree of stiffness of the differential equations. Moreover, a comparison indicates that the improved SoftCADS shortens the calculation time of the dynamic sensitivities without reducing their accuracies, even when the simplest approximate derivative formula is used.

AB - We have previously developed the software for calculation of dynamic sensitivities, SoftCADS, in which one can calculate dynamic sensitivities with high accuracy by just setting the differential equations for metabolite concentrations. However, SoftCADS did not always provide calculated values with the machine accuracy of a computer, although a Taylor series method was employed to numerically solve the differential equations. This is because numerical derivatives calculated from an approximate formula were directly used in the derivation of the differential equations for sensitivities from those for metabolite concentrations. The present work therefore attempts to further enhance the performance of SoftCADS, including not only the accuracies of the calculated values but also the calculation time. To overcome the problem, the approximate formula is expanded into a Taylor series in time and the first-term value of the series is replaced by the exact coefficient on the second term of the flux function expanded into a Taylor series in an independent or dependent variable. The result reveals that this replacement certainly provides not only numerical derivatives but also dynamic sensitivities with superhigh accuracies comparable to the machine accuracy, regardless of the degree of stiffness of the differential equations. Moreover, a comparison indicates that the improved SoftCADS shortens the calculation time of the dynamic sensitivities without reducing their accuracies, even when the simplest approximate derivative formula is used.

UR - http://www.scopus.com/inward/record.url?scp=79961126416&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=79961126416&partnerID=8YFLogxK

U2 - 10.1016/j.mbs.2011.06.004

DO - 10.1016/j.mbs.2011.06.004

M3 - Article

VL - 233

SP - 59

EP - 67

JO - Mathematical Biosciences

JF - Mathematical Biosciences

SN - 0025-5564

IS - 1

ER -