Nonequilibrium work relation beyond the Boltzmann-Gibbs distribution

Discrete and continuous models of probability flux of switching dynamics: Uncovering stochastic oscillations in a toggle-switch system

The probability flux and velocity in stochastic reaction networks can help in characterizing dynamic changes in probability landscapes of these networks. Here, we study the behavior of three different models of probability flux, namely, the discrete flux model, the Fokker-Planck model, and a new continuum model of the Liouville flux. We compare these fluxes that are formulated based on, respectively, the chemical master equation, the stochastic differential equation, and the ordinary differential equation. We examine similarities and differences among these models at the nonequilibrium steady state for the toggle switch network under different binding and unbinding conditions. Our results show that at a strong stochastic condition of weak promoter binding, continuum models of Fokker-Planck and Liouville fluxes deviate significantly from the discrete flux model. Furthermore, we report the discovery of stochastic oscillation in the toggle-switch system occurring at weak binding conditions, a phenomenon captured only by the discrete flux model.

Discrete flux and velocity fields of probability and their global maps in reaction systems

Stochasticity plays important roles in reaction systems. Vector fields of probability flux and velocity characterize time-varying and steady-state properties of these systems, including high probability paths, barriers, checkpoints among different stable regions, as well as mechanisms of dynamic switching among them. However, conventional fluxes on continuous space are ill-defined and are problematic when at the boundaries of the state space or when copy numbers are small. By re-defining the derivative and divergence operators based on the discrete nature of reactions, we introduce new formulations of discrete fluxes. Our flux model fully accounts for the discreetness of both the state space and the jump processes of reactions. The reactional discrete flux satisfies the continuity equation and describes the behavior of the system evolving along directions of reactions. The species discrete flux directly describes the dynamic behavior in the state space of the reactants such as the transfer of probability mass. With the relationship between these two fluxes specified, we show how to construct time-evolving and steady-state global flow-maps of probability flux and velocity in the directions of every species at every microstate and how they are related to the outflow and inflow of probability fluxes when tracing out reaction trajectories. We also describe how to impose proper conditions enabling exact quantification of flux and velocity in the boundary regions, without the difficulty of enforcing artificial reflecting conditions. We illustrate the computation of probability flux and velocity using three model systems, namely, the birth-death process, the bistable Schlögl model, and the oscillating Schnakenberg model.

Escape rate for nonequilibrium processes dominated by strong non-detailed balance force

Quantifying the escape rate from a meta-stable state is essential to understand a wide range of dynamical processes. Kramers' classical rate formula is the product of an exponential function of the potential barrier height and a pre-factor related to the friction coefficient. Although many applications of the rate formula focused on the exponential term, the prefactor can have a significant effect on the escape rate in certain parameter regions, such as the overdamped limit and the underdamped limit. There have been continuous interests to understand the effect of non-detailed balance on the escape rate; however, how the prefactor behaves under strong non-detailed balance force remains elusive. In this work, we find that the escape rate formula has a vanishing prefactor with decreasing friction strength under the strong non-detailed balance limit. We both obtain analytical solutions in specific examples and provide a derivation for more general cases. We further verify the result by simulations and propose a testable experimental system of a charged Brownian particle in electromagnetic field. Our study demonstrates that a special care is required to estimate the effect of prefactor on the escape rate when non-detailed balance force dominates.

Endogenous Molecular-Cellular Network Cancer Theory: A Systems Biology Approach

In light of ever apparent limitation of the current dominant cancer mutation theory, a quantitative hypothesis for cancer genesis and progression, endogenous molecular-cellular network hypothesis has been proposed from the systems biology perspective, now for more than 10 years. It was intended to include both the genetic and epigenetic causes to understand cancer. Its development enters the stage of meaningful interaction with experimental and clinical data and the limitation of the traditional cancer mutation theory becomes more evident. Under this endogenous network hypothesis, we established a core working network of hepatocellular carcinoma (HCC) according to the hypothesis and quantified the working network by a nonlinear dynamical system. We showed that the two stable states of the working network reproduce the main known features of normal liver and HCC at both the modular and molecular levels. Using endogenous network hypothesis and validated working network, we explored genetic mutation pattern in cancer and potential strategies to cure or relieve HCC from a totally new perspective. Patterns of genetic mutations have been traditionally analyzed by posteriori statistical association approaches in light of traditional cancer mutation theory. One may wonder the possibility of a priori determination of any mutation regularity. Here, we found that based on the endogenous network theory the features of genetic mutations in cancers may be predicted without any prior knowledge of mutation propensities. Normal hepatocyte and cancerous hepatocyte stable states, specified by distinct patterns of expressions or activities of proteins in the network, provide means to directly identify a set of most probable genetic mutations and their effects in HCC. As the key proteins and main interactions in the network are conserved through cell types in an organism, similar mutational features may also be found in other cancers. This analysis yielded straightforward and testable predictions on an accumulated and preferred mutation spectrum in normal tissue. The validation of predicted cancer state mutation patterns demonstrates the usefulness and potential of a causal dynamical framework to understand and predict genetic mutations in cancer. We also obtained the following implication related to HCC therapy, (1) specific positive feedback loops are responsible for the maintenance of normal liver and HCC; (2) inhibiting proliferation and inflammation-related positive feedback loops, and simultaneously inducing liver-specific positive feedback loop is predicated as the potential strategy to cure or relieve HCC; (3) the genesis and regression of HCC is asymmetric. In light of the characteristic property of the nonlinear dynamical system, we demonstrate that positive feedback loops must be existed as a simple and general molecular basis for the maintenance of phenotypes such as normal liver and HCC, and regulating the positive feedback loops directly or indirectly provides potential strategies to cure or relieve HCC.

Potential landscape of high dimensional nonlinear stochastic dynamics with large noise

Quantifying stochastic processes is essential to understand many natural phenomena, particularly in biology, including the cell-fate decision in developmental processes as well as the genesis and progression of cancers. While various attempts have been made to construct potential landscape in high dimensional systems and to estimate transition rates, they are practically limited to the cases where either noise is small or detailed balance condition holds. A general and practical approach to investigate real-world nonequilibrium systems, which are typically high-dimensional and subject to large multiplicative noise and the breakdown of detailed balance, remains elusive. Here, we formulate a computational framework that can directly compute the relative probabilities between locally stable states of such systems based on a least action method, without the necessity of simulating the steady-state distribution. The method can be applied to systems with arbitrary noise intensities through A-type stochastic integration, which preserves the dynamical structure of the deterministic counterpart dynamics. We demonstrate our approach in a numerically accurate manner through solvable examples. We further apply the method to investigate the role of noise on tumor heterogeneity in a 38-dimensional network model for prostate cancer, and provide a new strategy on controlling cell populations by manipulating noise strength.

Decoding early myelopoiesis from dynamics of core endogenous network

A decade ago mainstream molecular biologists regarded it impossible or biologically ill-motivated to understand the dynamics of complex biological phenomena, such as cancer genesis and progression, from a network perspective. Indeed, there are numerical difficulties even for those who were determined to explore along this direction. Undeterred, seven years ago a group of Chinese scientists started a program aiming to obtain quantitative connections between tumors and network dynamics. Many interesting results have been obtained. In this paper we wish to test such idea from a different angle: the connection between a normal biological process and the network dynamics. We have taken early myelopoiesis as our biological model. A standard roadmap for the cell-fate diversification during hematopoiesis has already been well established experimentally, yet little was known for its underpinning dynamical mechanisms. Compounding this difficulty there were additional experimental challenges, such as the seemingly conflicting hematopoietic roadmaps and the cell-fate inter-conversion events. With early myeloid cell-fate determination in mind, we constructed a core molecular endogenous network from well-documented gene regulation and signal transduction knowledge. Turning the network into a set of dynamical equations, we found computationally several structurally robust states. Those states nicely correspond to known cell phenotypes. We also found the states connecting those stable states. They reveal the developmental routes-how one stable state would most likely turn into another stable state. Such interconnected network among stable states enabled a natural organization of cell-fates into a multi-stable state landscape. Accordingly, both the myeloid cell phenotypes and the standard roadmap were explained mechanistically in a straightforward manner. Furthermore, recent challenging observations were also explained naturally. Moreover, the landscape visually enables a prediction of a pool of additional cell states and developmental routes, including the non-sequential and cross-branch transitions, which are testable by future experiments. In summary, the endogenous network dynamics provide an integrated quantitative framework to understand the heterogeneity and lineage commitment in myeloid progenitors.

Anomalous free energy changes induced by topology

We report that nontrivial topology of a driven Brownian particle restricted on a ring leads to anomalous behaviors on free energy change. Starting from steady states with identical distribution and current on the ring, free energy changes are distinct and nonperiodic after the system is driven by the same periodic force protocol. We demonstrate our observation in examples through both exact solutions and numerical simulations. The free energy calculated here can be measured in recent experimental systems.

Summing over trajectories of stochastic dynamics with multiplicative noise

We demonstrate that previous path integral formulations for the general stochastic interpretation generate incomplete results exemplified by the geometric Brownian motion. We thus develop a novel path integral formulation for the overdamped Langevin equation with multiplicative noise. The present path integral leads to the corresponding Fokker-Planck equation, and naturally generates a normalized transition probability in examples. Our result solves the inconsistency of the previous path integral formulations for the general stochastic interpretation, and can have wide applications in chemical and physical stochastic processes.

Controlling symmetry-breaking states by a hidden quantity in multiplicative noise

The inhomogeneity of multiplicative white noise leads to various coupling modes between deterministic and stochastic forces. We investigate the phase transition induced by the variation of the coupling mode through manipulating its characteristic parameter continuously. Even when the noise strength is fixed, an increase of this parameter can enhance or inhibit the symmetry-breaking state. We also propose a scheme to implement these phase transitions experimentally. Our result demonstrates that the coupling mode previously considered to be a mathematical convention serves as an additional quantity leading to physically observable phase transitions. This observation provides a mechanism to control the effect of noise without regulating the noise strength.