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

Amir Akbari,Zachary B. Haiman,Bernhard O. Palsson
DOI: https://doi.org/10.1128/msystems.01001-23
2024-01-24
mSystems
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.
microbiology
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to understand and analyze the dynamic characteristics exhibited by biological systems in a constantly changing environment. Specifically, the paper proposes a data - driven method for time - scale decomposition of biochemical reaction networks to distill fundamental patterns from complex dynamics. This method can identify time - scale hierarchies, concentration pools, and coherent structures, thereby providing a system - level description of the reaction network on physiologically important time - scales. The paper verifies the effectiveness of this method by analyzing models of hypothesized and actual biological pathways. In particular, by analyzing the time - scale hierarchy of the glycolytic pathway, the paper clarifies the connection between the stoichiometry and dissipative structure of the reaction network and the coherently structured time organization. The study found that glycolysis is a co - factor - driven pathway, and its slowest dynamic feature is described by the balance between high - energy phosphate bonds and redox transport. Overall, this method provides a more biologically interpretable characterization of network dynamics than large - scale kinetic models, which is helpful for the development of model simplification and personalized medical applications. The main contribution of the paper lies in the development of a computational framework - Dynamic Mode Analysis (DMA). This method combines the advantages of top - down and bottom - up techniques and can determine the time - scale hierarchy of the reaction network from experimentally measured or numerically generated time - series data. DMA can not only reliably identify concentration pools in complex networks with uncharacterized rate laws like top - down techniques, but also provide a mechanistic description of concentration pools and their relationship to the energetics and stoichiometric organization of the reaction network like bottom - up techniques.