Landscape and flux quantify the stochastic transition dynamics for p53 cell fate decision

Extrinsic fluctuations in the p53 cycle

Fluctuations are inherent to biological systems, arising from the stochastic nature of molecular interactions, and influence various aspects of system behavior, stability, and robustness. These fluctuations can be categorized as intrinsic, stemming from the system's inherent structure and dynamics, and extrinsic, arising from external factors, such as temperature variations. Understanding the interplay between these fluctuations is crucial for obtaining a comprehensive understanding of biological phenomena. However, studying these effects poses significant computational challenges. In this study, we used an underexplored methodology to analyze the effect of extrinsic fluctuations in stochastic systems using ordinary differential equations instead of solving the master equation with stochastic parameters. By incorporating temperature fluctuations into reaction rates, we explored the impact of extrinsic factors on system dynamics. We constructed a master equation and calculated the equations for the dynamics of the first two moments, offering computational efficiency compared with directly solving the chemical master equation. We applied this approach to analyze a biological oscillator, focusing on the p53 model and its response to temperature-induced extrinsic fluctuations. Our findings underscore the impact of extrinsic fluctuations on the nature of oscillations in biological systems, with alterations in oscillatory behavior depending on the characteristics of extrinsic fluctuations. We observed an increased oscillation amplitude and frequency of the p53 concentration cycle. This study provides valuable insights into the effects of extrinsic fluctuations on biological oscillations and highlights the importance of considering them in more complex systems to prevent unwanted scenarios related to health issues.

Inertial effect of cell state velocity on the quiescence-proliferation fate decision

Energy landscapes can provide intuitive depictions of population heterogeneity and dynamics. However, it is unclear whether individual cell behavior, hypothesized to be determined by initial position and noise, is faithfully recapitulated. Using the p21-/Cdk2-dependent quiescence-proliferation decision in breast cancer dormancy as a testbed, we examined single-cell dynamics on the landscape when perturbed by hypoxia, a dormancy-inducing stress. Combining trajectory-based energy landscape generation with single-cell time-lapse microscopy, we found that a combination of initial position and velocity on a p21/Cdk2 landscape, but not position alone, was required to explain the observed cell fate heterogeneity under hypoxia. This is likely due to additional cell state information such as epigenetic features and/or other species encoded in velocity but missing in instantaneous position determined by p21 and Cdk2 levels alone. Here, velocity dependence manifested as inertia: cells with higher cell cycle velocities prior to hypoxia continued progressing along the cell cycle under hypoxia, resisting the change in landscape towards cell cycle exit. Such inertial effects may markedly influence cell fate trajectories in tumors and other dynamically changing microenvironments where cell state transitions are governed by coordination across several biochemical species.

Macrophage phenotype transitions in a stochastic gene-regulatory network model

Polarization is the process by which a macrophage cell commits to a phenotype based on external signal stimulation. To know how this process is affected by random fluctuations and events within a cell is of utmost importance to better understand the underlying dynamics and predict possible phenotype transitions. For this purpose, we develop a stochastic modeling approach for the macrophage polarization process. We classify phenotype states using the Robust Perron Cluster Analysis and quantify transition pathways and probabilities by applying Transition Path Theory. Depending on the model parameters, we identify four bistable and one tristable phenotype configuration. We find that bistable transitions are fast but their states less robust. In contrast, phenotype transitions in the tristable situation have a comparatively long time duration, which reflects the robustness of the states. The results indicate parallels in the overall transition behavior of macrophage cells with other heterogeneous and plastic cell types, such as cancer cells. Our approach allows for a probabilistic interpretation of macrophage phenotype transitions and biological inference on phenotype robustness. In general, the methodology can easily be adapted to other systems where random state switches are known to occur.

Inertial effect of cell state velocity on the quiescence-proliferation fate decision in breast cancer

Energy landscapes can provide intuitive depictions of population heterogeneity and dynamics. However, it is unclear whether individual cell behavior, hypothesized to be determined by initial position and noise, is faithfully recapitulated. Using the p21-/Cdk2-dependent quiescence-proliferation decision in breast cancer dormancy as a testbed, we examined single-cell dynamics on the landscape when perturbed by hypoxia, a dormancy-inducing stress. Combining trajectory-based energy landscape generation with single-cell time-lapse microscopy, we found that initial position on a p21/Cdk2 landscape did not fully explain the observed cell-fate heterogeneity under hypoxia. Instead, cells with higher cell state velocities prior to hypoxia, influenced by epigenetic parameters, tended to remain proliferative under hypoxia. Thus, the fate decision on this landscape is significantly influenced by "inertia", a velocity-dependent ability to resist directional changes despite reshaping of the underlying landscape, superseding positional effects. Such inertial effects may markedly influence cell-fate trajectories in tumors and other dynamically changing microenvironments.

An improved approach for calculating energy landscape of gene networks from moment equations

The energy landscape theory has widely been applied to study the stochastic dynamics of biological systems. Different methods have been developed to quantify the energy landscape for gene networks, e.g., using Gaussian approximation (GA) approach to calculate the landscape by solving the diffusion equation approximately from the first two moments. However, how high-order moments influence the landscape construction remains to be elucidated. Also, multistability exists extensively in biological networks. So, how to quantify the landscape for a multistable dynamical system accurately, is a paramount problem. In this work, we prove that the weighted summation from GA (WSGA), provides an effective way to calculate the landscape for multistable systems and limit cycle systems. Meanwhile, we proposed an extended Gaussian approximation (EGA) approach by considering the effects of the third moments, which provides a more accurate way to obtain probability distribution and corresponding landscape. By applying our generalized EGA approach to two specific biological systems: multistable genetic circuit and synthetic oscillatory network, we compared EGA with WSGA by calculating the KL divergence of the probability distribution between these two approaches and simulations, which demonstrated that the EGA provides a more accurate approach to calculate the energy landscape.

Sampling rare trajectories using stochastic bridges

The numerical quantification of the statistics of rare events in stochastic processes is a challenging computational problem. We present a sampling method that constructs an ensemble of stochastic trajectories that are constrained to have fixed start and end points (so-called stochastic bridges). We then show that by carefully choosing a set of such bridges and assigning an appropriate statistical weight to each bridge, one can focus more processing power on the rare events of a target stochastic process while faithfully preserving the statistics of these rare trajectories. Further, we also compare the stochastic bridges we produce to the Wentzel-Kramers-Brillouin (WKB) optimal paths of the target process, derived in the limit of low noise. We see that the generated paths, encoding the full statistics of the process, collapse onto the WKB optimal path as the level of noise is reduced. We propose that the method can also be used to judge the accuracy of the WKB approximation at finite levels of noise.

Unraveling the stochastic transition mechanism between oscillation states by landscape and minimum action path theory

Cell fate transitions have been studied from various perspectives, such as the transition between stable states, or the transition between stable states and oscillation states. However, there is a lack...

Quantifying the Landscape of Decision Making From Spiking Neural Networks

The decision making function is governed by the complex coupled neural circuit in the brain. The underlying energy landscape provides a global picture for the dynamics of the neural decision making system and has been described extensively in the literature, but often as illustrations. In this work, we explicitly quantified the landscape for perceptual decision making based on biophysically-realistic cortical network with spiking neurons to mimic a two-alternative visual motion discrimination task. Under certain parameter regions, the underlying landscape displays bistable or tristable attractor states, which quantify the transition dynamics between different decision states. We identified two intermediate states: the spontaneous state which increases the plasticity and robustness of changes of minds and the "double-up" state which facilitates the state transitions. The irreversibility of the bistable and tristable switches due to the probabilistic curl flux demonstrates the inherent non-equilibrium characteristics of the neural decision system. The results of global stability of decision-making quantified by barrier height inferred from landscape topography and mean first passage time are in line with experimental observations. These results advance our understanding of the stochastic and dynamical transition mechanism of decision-making function, and the landscape and kinetic path approach can be applied to other cognitive function related problems (such as working memory) in brain networks.