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.

Study of the structure and dynamics of complex biological networks

Modeling biological networks: from single gene systems to large microbial communities

The relationship between structure and function in locally observed complex networks

Size Dependent Growth in Metabolic Networks

Modeling and Analyzing Biological Oscillations in Molecular Networks

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

Function Constrains Network Architecture and Dynamics: a Case Study on the Yeast Cell Cycle Boolean Network.

The large-scale organization of metabolic networks

The regulatory network of E. coli metabolism as a Boolean dynamical system exhibits both homeostasis and flexibility of response

Network analyses in systems biology: new strategies for dealing with biological complexity

An enzyme-centric approach for modelling non-linear biological complexity

A system-wide network reconstruction of gene regulation and metabolism in Escherichia coli

Discovering universal statistical laws of complex networks

Network dynamics and cell physiology

Metabolic Network Modelling: Including Stochastic Effects.

Network biology: understanding the cell's functional organization

On the basic computational structure of gene regulatory networks

Reconstruction of metabolic networks from genome data and analysis of their global structure for various organisms

Communication on structure of biological networks

The role of modularity in self-organisation dynamics in biological networks

Systems biology: From the cell to the brain