Abstract:ABSTRACT Understanding the dynamics of biological systems in evolving environments is a challenge due to their scale and complexity. Here, we present a computational framework for the timescale decomposition of biochemical reaction networks to distill essential patterns from their intricate dynamics. This approach identifies timescale hierarchies, concentration pools, and coherent structures from time-series data, providing a system-level description of reaction networks at physiologically important timescales. We apply this technique to kinetic models of hypothetical and biological pathways, validating it by reproducing analytically characterized or previously known concentration pools of these pathways. Moreover, by analyzing the timescale hierarchy of the glycolytic pathway, we elucidate the connections between the stoichiometric and dissipative structures of reaction networks and the temporal organization of coherent structures. Specifically, we show that glycolysis is a cofactor-driven pathway, the slowest dynamics of which are described by a balance between high-energy phosphate bond and redox trafficking. Overall, this approach provides more biologically interpretable characterizations of network dynamics than large-scale kinetic models, thus facilitating model reduction and personalized medicine applications. IMPORTANCE Complex interactions within interconnected biochemical reaction networks enable cellular responses to a wide range of unpredictable environmental perturbations. Understanding how biological functions arise from these intricate interactions has been a long-standing problem in biology. Here, we introduce a computational approach to dissect complex biological systems' dynamics in evolving environments. This approach characterizes the timescale hierarchies of complex reaction networks, offering a system-level understanding at physiologically relevant timescales. Analyzing various hypothetical and biological pathways, we show how stoichiometric properties shape the way energy is dissipated throughout reaction networks. Notably, we establish that glycolysis operates as a cofactor-driven pathway, where the slowest dynamics are governed by a balance between high-energy phosphate bonds and redox trafficking. This approach enhances our understanding of network dynamics and facilitates the development of reduced-order kinetic models with biologically interpretable components.

Simplicity of Completion Time Distributions for Common Complex Biochemical Processes

Characteristic, completion or matching timescales? An analysis of temporary boundaries in enzyme kinetics

Short relaxation times but long transient times in both simple and complex reaction networks

A data-driven approach for timescale decomposition of biochemical reaction networks

Patterns of stochastic behavior in dynamically unstable high-dimensional biochemical networks

Coarse-graining stochastic biochemical networks: quasi-stationary approximation and fast simulations using a stochastic path integral technique

Modelling non-Markovian dynamics in biochemical reactions

A path method for non-exponential ergodicity of Markov chains and its application for chemical reaction systems

Strong defocusing of molecular reaction times results from an interplay of geometry and reaction control

Statistical Analyses And Theoretical Models Of Single-Molecule Enzymatic Dynamics

Single-molecule Michaelis-Menten Equations.

The structure of inter-reaction times in reaction-diffusion processes and consequences for counting statistics

Analytical Refractory Period Distribution for a Class of Time-Variant Biochemical Systems with Second-Order Reactions.

Robust Stochastic Chemical Reaction Networks and Bounded Tau-Leaping

Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art

Multiscale Analysis of Reaction Networks

Coarse-grained entropy production with multiple reservoirs: unraveling the role of time-scales and detailed balance in biology-inspired systems

Noise and critical phenomena in biochemical signaling cycles at small molecule numbers

Single Molecule Michaelis-Menten Equation beyond Quasi-Static Disorder

Diffusion Controlled Reactions, Fluctuation Dominated Kinetics, and Living Cell Biochemistry

Complex balanced distributions for chemical reaction networks