×

The QSSA in chemical kinetics: as taught and as practiced. (English) Zbl 1290.92013

Jonoska, Nataša (ed.) et al., Discrete and topological models in molecular biology. Berlin: Springer (ISBN 978-3-642-40192-3/hbk; 978-3-642-40193-0/ebook). Natural Computing Series, 419-442 (2014).
Summary: Chemical mechanisms for even simple reaction networks involve many highly reactive and short-lived species (intermediates), present in small concentrations, in addition to the main reactants and products, present in larger concentrations. The chemical mechanism also often contains many rate constants whose values are unknown a priori and must be determined from experimental measurements of the large species concentrations. A classic model reduction method known as the quasi-steady-state assumption (QSSA) is often used to eliminate the highly reactive intermediate species and remove the large rate constants that cannot be determined from concentration measurements of the reactants and products. Mathematical analysis based on the QSSA is ubiquitous in modeling enzymatic reactions. In this chapter, we focus attention on the QSSA, how it is “taught” to students of chemistry, biology, and chemical and biological engineering, and how it is “practiced” when researchers confront realistic and complex examples. We describe the main types of difficulties that appear when trying to apply the standard ideas of the QSSA, and propose a new strategy for overcoming them, based on rescaling the reactive intermediate species.
First, we prove mathematically that the program taught to beginning students for applying the 100-year-old approach of classic QSSA model reduction cannot be carried out for many of the relevant kinetics problems, and perhaps even most of them. By using Galois theory, we prove that the required algebraic equations cannot be solved for as few as five bimolecular reactions between five species (with three intermediates). We expect that many practitioners have suspected this situation regarding nonsolvability to exist, but we have seen no statement or proof of this fact, especially when the kinetics are restricted to unimolecular and bimolecular reactions. We describe algorithms that can test any mechanism for solvability. We also show that an alternative to solving the QSSA equations, the Horiuti-Temkin theory, also does not work for many examples.
Of course, the reduced model (and the full model, for that matter) can be solved numerically, which is the standard approach in practice. The remaining difficulty, however, is how to obtain the values of the large kinetic parameters appearing in the model. These parameters cannot be estimated from measurements of the large-concentration reactants and products. We show here how the concept of rescaling the reactive intermediate species allows the large kinetic parameters to be removed from the parameter estimation problem. In general, the number of parameters that can be removed from the full model is less than or equal to the number of intermediate species. The outcome is a reduced model with a set of rescaled parameters that is often identifiable from routinely available measurements. New and freely available computational software (parest_dae) for estimating the reduced model’s kinetic parameters and confidence intervals is briefly described.
For the entire collection see [Zbl 1279.92005].

MSC:

92C45 Kinetics in biochemical problems (pharmacokinetics, enzyme kinetics, etc.)

Software:

parest_dae
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Abel, N. H., Mémoire sur les équations algébriques, où l’on démontre l’impossibilité de la résolution de l’équation générale du cinquième degré in Oeuvres complètes de Niels Henrik Abel, Édition de Christiana, 28-33 (1881), Oslo: Gröndahl and Son, Oslo
[2] Atkinson, A. C.; Donev, A. N., Optimum Experimental Designs (1992), New York: Oxford University Press, New York · Zbl 0829.62070
[3] Bard, Y., Nonlinear Parameter Estimation (1974), New York: Academic, New York · Zbl 0345.62045
[4] Bodenstein, M., Eine Theorie der photochemischen Reaktionsgeschwindigkeiten, Z. Phys. Chem., 85, 329-397 (1913) · JFM 44.0905.03
[5] Chapman, D. L.; Underhill, L. K., The interaction of chlorine and hydrogen, The influence of mass. J. Chem. Soc. Trans., 103, 496-508 (1913) · doi:10.1039/ct9130300496
[6] Ciliberto, Andrea; Capuani, Fabrizio; Tyson, John J., Modeling Networks of Coupled Enzymatic Reactions Using the Total Quasi-Steady State Approximation, PLoS Computational Biology, 3, 3, e45 (2007)
[7] Cox, D.; Little, J.; O’Shea, D., Ideals, Varieties, and Algorithms: An Introduction to Computational Algebraic Geometry and Commutative Algebra (2006), New York: Springer, New York
[8] Craciun, G.; Pantea, C., Identifiability of chemical reaction networks, J. Math. Chem., 44, 244-259 (2007) · Zbl 1145.92040 · doi:10.1007/s10910-007-9307-x
[9] Goldbeter, A.; Koshland, D. E., An amplified sensitivity arising from covalent modification in biological systems, Proc. Natl. Acad. Sci. U.S.A., 78, 11, 6840-6844 (1981) · doi:10.1073/pnas.78.11.6840
[10] Hanson, S.; Schnell, S., Reactant stationary approximation in enzyme kinetics, J. Phys. Chem. A, 112, 37, 8654-8658 (2008) · doi:10.1021/jp8026226
[11] Hill, A. V., The possible effects of the aggregation of the molecules of hemoglobin on its dissociation curves, J. Physiol., 40, 4-7 (1910)
[12] Hindmarsh, A. C.; Brown, P. N.; Grant, K. E.; Lee, S. L.; Serban, R.; Shumaker, D. E.; Woodward, C. S., SUNDIALS: suite of nonlinear and differential/algebraic equation solvers, ACM Trans. Math. Softw., 31, 3, 363-396 (2005) · Zbl 1136.65329 · doi:10.1145/1089014.1089020
[13] Horiuti, J.; Nakamura, T., Stoichiometric number and the theory of steady reaction, Z. Phys. Chem., 11, 358-365 (1957) · doi:10.1524/zpch.1957.11.5_6.358
[14] Isaacs, M., Algebra: A Graduate Course (1993), Pacific Grove: Brooks Cole, Pacific Grove · Zbl 0805.00001
[15] Khovanskii, A. G., On solvability and unsolvability of equations in explicit form, Russ. Math. Surv., 59, 4, 661-736 (2004) · Zbl 1069.34133 · doi:10.1070/RM2004v059n04ABEH000759
[16] Lam, S. H.; Goussis, D. A., The CSP method for simplifying kinetics, Int. J. Chem. Kinet., 26, 461-486 (1994) · doi:10.1002/kin.550260408
[17] Maple is a division of Waterloo Maple, Inc. (2007)
[18] J.D. Murray, Mathematical Biology: I. An Introduction (Springer, New York, 2002) · Zbl 1006.92001
[19] Rawlings, J. B.; Ekerdt, J. G., Chemical Reactor Analysis and Design Fundamentals (2004), Madison: Nob Hill Publishing, Madison
[20] Schnell, S.; Maini, P. K., A century of enzyme kinetics: reliability of the K_M and v_max estimates, Comments Theor. Biol., 8, 169-187 (2003) · doi:10.1080/08948550302453
[21] Segel, L. A.; Slemrod, M., The quasi-steady state assumption: a case study in perturbation, SIAM Rev., 31, 3, 446-477 (1989) · Zbl 0679.34066 · doi:10.1137/1031091
[22] Temkin, M. I., The kinetics of some industrial heterogeneous catalytic reactions, Adv. Catal., 28, 173-291 (1979) · doi:10.1016/S0360-0564(08)60135-2
[23] Turanyi, T.; Tomlin, A. S.; Pilling, M. J., On the error of the quasi-steady-state approximation, J. Phys. Chem., 97, 1, 163-172 (1993) · doi:10.1021/j100103a028
[24] Tyson, J. J.; Chen, K. C.; Novak, B., Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell, Curr. Opin. Cell Biol., 15, 221-231 (2003) · doi:10.1016/S0955-0674(03)00017-6
[25] Wright, M. R., An Introduction to Chemical Kinetics (2004), Chichester/Hoboken: Wiley, Chichester/Hoboken · doi:10.1002/047009060X
[26] Zagaris, A.; Kaper, H. G.; Kaper, T. J., Analysis of the computational singular perturbation reduction method for chemical kinetics, J. Nonlinear Sci., 14, 59-91 (2004) · Zbl 1053.92051 · doi:10.1007/s00332-003-0582-9
[27] Zagaris, A.; Kaper, H. G.; Kaper, T. J., Fast and slow dynamics for the computational singular perturbation method, Multiscale Model. Simul., 2, 4, 613-638 (2004) · Zbl 1065.34049 · doi:10.1137/040603577
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.