Mesoscopic Biochemical Basis of Isogenetic Inheritance and Canalization: Stochasticity, Nonlinearity, and Emergent Landscape

Biochemical reaction systems in mesoscopic volume, under sustained environmental chemical gradient(s), can have multiple stochastic attractors. Two distinct mechanisms are known for their origins: ($a$) Stochastic single-molecule events, such as gene expression, with slow gene on-off dynamics; and ($b$) nonlinear networks with feedbacks. These two mechanisms yield different volume dependence for the sojourn time of an attractor. As in the classic Arrhenius theory for temperature dependent transition rates, a landscape perspective provides a natural framework for the system’s behavior. However, due to the nonequilibrium nature of the open chemical systems, the landscape, and the attractors it represents, are all themselves {\em emergent properties} of complex, mesoscopic dynamics. In terms of the landscape, we show a generalization of Kramers’ approach is possible to provide a rate theory. The emergence of attractors is a form of self-organization in the mesoscopic system; stochastic attractors in biochemical systems such as gene regulation and cellular signaling are naturally inheritable via cell division. Delbr"{u}ck-Gillespie’s mesoscopic reaction system theory, therefore, provides a biochemical basis for spontaneous isogenetic switching and canalization.

💡 Research Summary

The paper investigates how mesoscopic biochemical reaction systems—those operating in cellular volumes on the order of femtoliters—can exhibit multiple stochastic attractors when subjected to sustained environmental chemical gradients. Two distinct mechanisms for the emergence of these attractors are identified. The first mechanism (a) involves stochastic single‑molecule events, such as transcription factor binding/unbinding or gene on/off switching, which occur on a slow timescale relative to other reactions. In this regime the residence time of an attractor scales inversely with system volume; the transition rate follows an Arrhenius‑like expression k ≈ k₀ exp(−ΔG/RT) · V⁻¹, so that larger volumes dramatically suppress rare switching. The second mechanism (b) relies on nonlinear network feedbacks (positive or negative) that generate multistability. Here the dynamics are described by stochastic differential equations whose small‑noise limit yields a non‑equilibrium potential (or “landscape”) Φ(x) ≈ −ε ln P_ss(x). The transition rate between attractors follows a generalized Kramers formula k ≈ A exp(−ΔΦ/ε), where the effective noise strength ε scales as V⁻¹. Consequently, the volume dependence is non‑linear (often a fractional power) rather than a simple inverse proportionality.

Because the systems are open and driven far from equilibrium, the landscape itself is an emergent property: it is continuously reshaped by the external chemical fluxes and internal reaction currents, and it contains both gradient (potential) and rotational (non‑conservative) components. The authors show mathematically that, despite cell division halving the volume, the shape of the landscape and the positions of its wells remain essentially unchanged. Thus, daughter cells inherit the same stochastic attractor distribution as the mother cell—a phenomenon the authors term “isogenetic switching.” This inheritance is demonstrated through Gillespie simulations that reproduce the predicted volume‑dependent transition statistics before and after division.

The paper further connects the landscape framework to Waddington’s concept of canalization. A sustained external gradient can deepen particular wells in the landscape, making the system robust to fluctuations and guiding developmental trajectories along preferred paths. This provides a mechanistic explanation for how cellular differentiation can be both flexible (allowing stochastic switches) and robust (maintaining a canalized fate) within the same biochemical network.

Finally, the theoretical predictions are compared with experimental observations in systems such as λ‑phage lysogeny/lysis switching, synthetic gene toggle switches, and MAPK signaling cascades. Measured residence times, probability distributions of expression states, and the preservation of phenotypic states across cell divisions align with the model’s volume‑dependent scaling laws and landscape‑based rate theory. The authors conclude that Delbrück‑Gillespie mesoscopic reaction theory offers a unified biochemical basis for spontaneous isogenetic switching and developmental canalization, bridging stochastic molecular events, nonlinear network dynamics, and emergent non‑equilibrium landscapes.