Abstract:We established a theoretical framework for studying nonequilibrium networks with two distinct natures essential for characterizing the global probabilistic dynamics: the underlying potential landscape and the corresponding curl flux. We applied the idea to two biological examples: biochemical oscillation network and stem-cell fate decision network for differentiation and development.We found that the underlying potential landscape for the oscillation limit cycle has a distinct closed ring valley (Mexican hat-like) shape when the fluctuations are small. This global landscape structure leads to attractions of the system to the ring valley. On the ring, we found that the nonequilibrium flux is the driving force for oscillations. Therefore, both structured landscape and flux are needed to guarantee a robust oscillating network. The barrier height separating the oscillation ring and other areas derived from the landscape topography is shown to be correlated with the escaping time from the limit cycle attractor and provides a quantitative measure of the robustness for the network. The landscape becomes shallower and the closed ring valley shape structure becomes weaker (lower barrier height) with larger fluctuations. We observe that the period and the amplitude of the oscillations are more dispersed and oscillations become less coherent when the fluctuations increase. We also found that the entropy production of the whole network, characterizing the dissipation costs from the combined effects of both landscapes and fluxes, decreases when the fluctuations decrease. Therefore, less dissipation leads to more robust networks. Our approach is quite general and applicable to other networks, dynamical systems, and biological evolution. It can help designing robust networks.Differentiation from a multipotent stem or progenitor state to a mature cell is an essential irreversible process. The associated changes in gene expression patterns exhibit time directionality. This "arrow of time" in the collective change of gene expression across multiple stable gene expression patterns (attractors) is neither explained by the regulated activation, or suppression of individual genes, which are bidirectional molecular processes, nor by standard dynamical models of the underlying gene circuit which only account for local stability of attractors. To capture the global dynamics of this nonequilibrium system and gain insight in the time asymmetry of state transitions, we computed the quasipotential landscape of the stochastic dynamics of a canonical gene circuit that governs branching cell fate commitment. The potential landscape reveals the global dynamics and permits the calculation of "potential barriers" between cell phenotypes imposed by the circuit architecture. The generic asymmetry of barrier heights indicates that the transition from the uncommitted multipotent state to differentiated states is inherently unidirectional. The model agrees with observations and predicts the extreme conditions for reprogramming cells back to the undifferentiated state.

Approximating the Energy Landscape of a Two-Dimensional Bistable Gene Autoregulation Model by Separating Slow and Fast Dynamics.

An Improved Approach for Calculating Energy Landscape of Gene Networks from Moment Equations

Bifurcation Analysis and Probabilistic Energy Landscapes of Two-Component Genetic Network

Quasi-potential landscape in complex multi-stable systems

Constructing the Energy Landscape for Genetic Switching System Driven by Intrinsic Noise

Free Energy Landscape for Biological Systems

Quantifying energy landscape of oscillatory systems: Explosion, pre-solution, and diffusion decomposition

Transition Path, Quasi-potential Energy Landscape and Stability of Genetic Switches

Estimation of Non-Equilibrium Transition Rate from Gene Expression Data

Potential Landscape And Flux Framework Of Nonequilibrium Biological Networks

Eddy Current and Coupled Landscapes for Nonadiabatic and Nonequilibrium Complex System Dynamics

Energy Landscape Decomposition for Cell Differentiation with Proliferation Effect.

The Potential and Flux Landscape Theory of Evolution

Kinetic Paths, Time Scale, and Underlying Landscapes: a Path Integral Framework to Study Global Natures of Nonequilibrium Systems and Networks.

Multiple Coupled Landscapes and Non-Adiabatic Dynamics with Applications to Self-Activating Genes

Potential landscape of high dimensional nonlinear stochastic dynamics and rare transitions with large noise

Construction of a Pathway Map on a Complicated Energy Landscape.

Towards Stable Kinetics of Large Metabolic Networks: Nonequilibrium Potential Function Approach

Dynamic landscapes and statistical limits on growth during cell fate specification

The Energy Pump and the Origin of the Non-Equilibrium Flux of the Dynamical Systems and the Networks

Landscape Framework and Global Stability for Stochastic Reaction Diffusion and General Spatially Extended Systems with Intrinsic Fluctuations