A Numerical Aggregation Algorithm for the Enzyme-Catalyzed Substrate Conversion

摘要

Computational models of biochemical systems are usually very large, and moreover, if reaction frequencies of different reaction types differ in orders of magnitude, models possess the mathematical property of stiffness, which renders system analysis difficult and often even impossible with traditional methods. Recently, an accelerated stochastic simulation technique based on a system partitioning, the slow-scale stochastic simulation algorithm, has been applied to the enzyme-catalyzed substrate conversion to circumvent the inefficiency of standard stochastic simulation in the presence of stiffness. We propose a numerical algorithm based on a similar partitioning but without resorting to simulation. The algorithm exploits the connection to continuous-time Markov chains and decomposes the overall problem to significantly smaller sub-problems that become tractable. Numerical results show enormous efficiency improvements relative to accelerated stochastic simulation.