Growth and depletion in linear stochastic reaction networks

Biological systems often involve transformations among substances through reactions. These processes can be modeled as continuous-time Markov jump processes. Prior research has shown that such models can represent a wide range of biological networks. However, the behavior of these systems under growth assumptions remains unclear. The mean-field behavior of these systems is typically described by ordinary differential equations. Yet, the connection between stochastic dynamics and deterministic approximations is not fully understood. The role of enzymes in reaction rates introduces another layer of complexity. No prior work had resolved how enzyme-dependent reaction rates affect long-term system behavior. This gap motivated the need for a rigorous analysis of such networks under growth conditions.

The goal of this study is to analyze the dynamics of stochastic reaction networks with enzyme-dependent rates. The researchers aim to understand how these systems behave under exponential growth assumptions. They explore whether mean-field ODE solutions can approximate the stochastic trajectories. The study also seeks to characterize the ω-limit sets of these systems. A second objective is to investigate depletion dynamics after a substance becomes unavailable. The researchers propose a mean-field description for this phase. They aim to provide a foundation for understanding bifurcations in reaction networks. This approach may clarify how intermittent production affects long-term behavior. The study addresses a key uncertainty in modeling biological systems with stochastic dynamics.

The researchers use continuous-time Markov jump processes to model the reaction networks. They assume exponential growth in the network size to simplify analysis. The mean-field behavior is described using ordinary differential equations. The team derives conditions under which stochastic trajectories can be approximated by ODE solutions. They apply rescaling techniques to analyze network behavior over time. The study also treats the model as a slow-fast system to study depletion dynamics. The researchers identify ω-limit sets as points or random tori. They examine how the number of enzymes affects dimension reduction in the system.

The study shows that ODE solutions can approximate stochastic trajectories under exponential growth. The researchers identify specific conditions for this approximation to hold with high probability. They characterize the ω-limit sets as either points or random tori. The results depend on the number of enzymes in the system. Depletion dynamics are analyzed after a substance becomes unavailable. The model reveals that the depleted substance can be produced intermittently. The researchers propose a mean-field description for this phase. Their analysis suggests a natural bifurcation in reaction networks.

The authors propose that ODE solutions can approximate stochastic trajectories under exponential growth. They suggest that ω-limit sets are either points or random tori. The study notes that dimension reduction depends on the number of enzymes. The researchers propose a mean-field description for depletion dynamics. They suggest that this phase involves intermittent production of the depleted substance. The analysis provides a first step toward understanding natural bifurcations in reaction networks. The results may help clarify how enzyme-dependent rates affect long-term behavior. The study highlights the importance of growth assumptions in modeling biological systems.