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.

Computational translation framework identifies biochemical reaction networks with special topologies and their long-term dynamics

Translated Chemical Reaction Networks

A Scalable Computational Framework for Establishing Long-Term Behavior of Stochastic Reaction Networks

Network transformation-based analysis of biochemical systems

Network translation and steady state properties of chemical reaction systems

A Stochastic Framework for the Design of Transient and Steady State Behavior of Biochemical Reaction Networks

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

Reduction of dynamical biochemical reaction networks in computational biology

A portable structural analysis library for reaction networks

Identification of Parameters and Initial Values for Reaction-Diffusion Systems in Protein Networks (Extended Version)

An Ingenious Characterization of Reaction Network Using Sub-Network Reconstruction.

Characterization of Biologically Relevant Network Structures form Time-series Data

A unique transformation from ordinary differential equations to reaction networks

A Computational Approach to Steady State Correspondence of Regular and Generalized Mass Action Systems

A Decomposition-Based Approach for Deriving Positive Steady States of a Class of Chemical Reaction Networks with Non-Mass-Action Kinetics

Reduction of Chemical Reaction Networks with Approximate Conservation Laws

React-OT: Optimal Transport for Generating Transition State in Chemical Reactions

Multi-species consensus network of DNA strand displacement for concentration-to-strand translation

Analytical Equilibrium Solutions of Biochemical Systems with Synthesis and Degradation

Optimality-preserving Reduction of Chemical Reaction Networks