Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOn the validity of the stochastic quasi-steady-state approximation in open enzyme catalyzed reactions: Timescale separation or singular perturbation?
70
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The quasi-steady-state approximation is widely used to develop simplified deterministic or stochastic models of enzyme catalyzed reactions. In deterministic models, the quasi-steady-state approximation can be mathematically justified from singular perturbation theory. For several closed enzymatic reactions, the homologous extension of the quasi-steady-state approximation to the stochastic regime, known as the stochastic quasi-steady-state approximation, has been shown to be accurate under the analogous conditions that permit the quasi-steady-state reduction of the deterministic counterpart. However, it was recently demonstrated that the extension of the stochastic quasi-steady-state approximation to an open Michaelis--Menten reaction mechanism is only valid under a condition that is far more restrictive than the qualifier that ensures the validity of its corresponding deterministic quasi-steady-state approximation. In this paper, we suggest a possible explanation for this discrepancy from the lens of geometric singular perturbation theory. In so doing, we illustrate a misconception in the application of the quasi-steady-state approximation: timescale separation does not imply singular perturbation.
rate research
Read More
The estimation of the kinetic parameters requires the careful design of experiments under a constrained set of conditions. Many estimates reported in the literature incorporate protocols that leverage simplified mathematical models known as quasi-steady-state reductions. Such reductions often - but not always - emerge as the result of a singular perturbation scenario. However, the utilization of the singular perturbation reduction method requires knowledge of a dimensionless parameter, $varepsilon$, that is proportional to the ratio of the reactions fast and slow timescales. Using techniques from differential equations, Fenichel theory, and center manifold theory, we derive the appropriate $varepsilon$ whose magnitude regulates the validity of the quasi-steady-state reduction employed in the reported experimental procedures for intermolecular autocatalytic zymogen activation reaction. Although the model equations are two-dimensional, the fast/slow dynamics are rich. The phase plane exhibits a dynamic transcritical bifurcation point in a particular singular limit. The existence of such a bifurcation is relevant, because the critical manifold losses normal hyperbolicity and classical Fenichel theory is inapplicable. Furthermore, we show that in some cases chemical reversibility can be interpreted dynamically as an imperfection, since the presence of reversibility can destroy the bifurcation structure present in the singular limit. We show that the reduction method by which QSS reductions are justified can depend on the path taken in parameter space. Specifically, we show that the standard quasi-steady-state reduction for this reaction is justifiable by center manifold theory in one limit, and via Fenichel theory in a different limit.
The linear noise approximation models the random fluctuations from the mean field model of a chemical reaction that unfolds near the thermodynamic limit. Specifically, the fluctuations obey a linear Langevin equation up to order $Omega^{-1/2}$, where $Omega$ is the size of the chemical system (usually the volume). Under the presence of disparate timescales, the linear noise approximation admits a quasi-steady-state reduction referred to as the slow scale linear noise approximation. However, the slow scale linear approximation has only been derived for fast/slow systems that are in Tikhonov standard form. In this work, we derive the slow scale linear noise approximation directly from Fenichel theory, without the need for a priori scaling and dimensional analysis. In so doing, we can apply for the first time the slow scale linear noise approximation to fast/slow systems that are not of standard form. This is important, because often times algorithms are only computationally expensive in parameter ranges where the system is singularly perturbed, but not in standard form. We also comment on the breakdown of the slow scale linear noise approximation near dynamic bifurcation points -- a topic that has remained absent in the chemical kinetics literature, despite the presence of bifurcations in simple biochemical reactions, such the Michaelis--Menten reaction mechanism.
The conditions for the validity of the standard quasi-steady-state approximation in the Michaelis--Menten mechanism in a closed reaction vessel have been well studied, but much less so the conditions for the validity of this approximation for the system with substrate inflow. We analyze quasi-steady-state scenarios for the open system attributable to singular perturbations, as well as less restrictive conditions. For both settings we obtain distinguished invariant slow manifolds and time scale estimates, and we highlight the special role of singular perturbation parameters in higher order approximations of slow manifolds. We close the paper with a discussion of distinguished invariant manifolds in the global phase portrait.
In biochemical networks, reactions often occur on disparate timescales and can be characterized as either fast or slow. The quasi-steady state approximation (QSSA) utilizes timescale separation to project models of biochemical networks onto lower-dimensional slow manifolds. As a result, fast elementary reactions are not modeled explicitly, and their effect is captured by non-elementary reaction rate functions (e.g. Hill functions). The accuracy of the QSSA applied to deterministic systems depends on how well timescales are separated. Recently, it has been proposed to use the non-elementary rate functions obtained via the deterministic QSSA to define propensity functions in stochastic simulations of biochemical networks. In this approach, termed the stochastic QSSA, fast reactions that are part of non-elementary reactions are not simulated, greatly reducing computation time. However, it is unclear when the stochastic QSSA provides an accurate approximation of the original stochastic simulation. We show that, unlike the deterministic QSSA, the validity of the stochastic QSSA does not follow from timescale separation alone, but also depends on the sensitivity of the non-elementary reaction rate functions to changes in the slow species. The stochastic QSSA becomes more accurate when this sensitivity is small. Different types of QSSAs result in non-elementary functions with different sensitivities, and the total QSSA results in less sensitive functions than the standard or the pre-factor QSSA. We prove that, as a result, the stochastic QSSA becomes more accurate when non-elementary reaction functions are obtained using the total QSSA. Our work provides a novel condition for the validity of the QSSA in stochastic simulations of biochemical reaction networks with disparate timescales.
In this work, we revisit the scaling analysis and commonly accepted conditions for the validity of the standard, reverse and total quasi-steady-state approximations through the lens of dimensional Tikhonov-Fenichel parameters and their respective critical manifolds. By combining Tikhonov-Fenichel parameters with scaling analysis and energy methods, we derive improved upper bounds on the approximation error for the standard, reverse and total quasi-steady-state approximations. Furthermore, previous analyses suggest that the reverse quasi-steady-state approximation is only valid when initial enzyme concentrations greatly exceed initial substrate concentrations. However, our results indicate that this approximation can be valid when initial enzyme and substrate concentrations are of equal magnitude. Using energy methods, we find that the condition for the validity of the reverse quasi-steady-state approximation is far less restrictive than was previously assumed, and we derive a new small parameter that determines the validity of this approximation. In doing so, we extend the established domain of validity for the reverse quasi-steady-state approximation. Consequently, this opens up the possibility of utilizing the reverse quasi-steady-state approximation to model enzyme catalyzed reactions and estimate kinetic parameters in enzymatic assays at much lower enzyme to substrate ratios than was previously thought. Moreover, we show for the first time that the critical manifold of the reverse quasi-steady-state approximation contains a singular point where normal hyperbolicity is lost. Associated with this singularity is a transcritical bifurcation, and the corresponding normal form of this bifurcation is recovered through scaling analysis.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
