The relationship between stochastic and deterministic quasi-steady state approximations

Coupled oscillator cooperativity as a control mechanism in chronobiology

Control of dynamical processes is vital for maintaining correct cell regulation and cell-fate decisions. Numerous regulatory networks show oscillatory behavior; however, our knowledge of how one oscillator behaves when stimulated by two or more external oscillatory signals is still missing. We explore this problem by constructing a synthetic oscillatory system in yeast and stimulate it with two external oscillatory signals. Letting model verification and prediction operate in a tight interplay with experimental observations, we find that stimulation with two external signals expands the plateau of entrainment and reduces the fluctuations of oscillations. Furthermore, by adjusting the phase differences of external signals, one can control the amplitude of oscillations, which is understood through the signal delay of the unperturbed oscillatory network. With this we reveal a direct amplitude dependency of downstream gene transcription. Taken together, these results suggest a new path to control oscillatory systems by coupled oscillator cooperativity.

Model reduction for the Chemical Master Equation: An information-theoretic approach

The complexity of mathematical models in biology has rendered model reduction an essential tool in the quantitative biologist's toolkit. For stochastic reaction networks described using the Chemical Master Equation, commonly used methods include time-scale separation, Linear Mapping Approximation, and state-space lumping. Despite the success of these techniques, they appear to be rather disparate, and at present, no general-purpose approach to model reduction for stochastic reaction networks is known. In this paper, we show that most common model reduction approaches for the Chemical Master Equation can be seen as minimizing a well-known information-theoretic quantity between the full model and its reduction, the Kullback-Leibler divergence defined on the space of trajectories. This allows us to recast the task of model reduction as a variational problem that can be tackled using standard numerical optimization approaches. In addition, we derive general expressions for propensities of a reduced system that generalize those found using classical methods. We show that the Kullback-Leibler divergence is a useful metric to assess model discrepancy and to compare different model reduction techniques using three examples from the literature: an autoregulatory feedback loop, the Michaelis-Menten enzyme system, and a genetic oscillator.

Enhanced DNA repair through droplet formation and p53 oscillations

Living organisms are constantly exposed to DNA damage, and optimal repair is therefore crucial. A characteristic hallmark of the response is the formation of sub-compartments around the site of damage, known as foci. Following multiple DNA breaks, the transcription factor p53 exhibits oscillations in its nuclear concentration, but how this dynamics can affect the repair remains unknown. Here, we formulate a theory for foci formation through droplet condensation and discover how oscillations in p53, with its specific periodicity and amplitude, optimize the repair process by preventing Ostwald ripening and distributing protein material in space and time. Based on the theory predictions, we reveal experimentally that the oscillatory dynamics of p53 does enhance the repair efficiency. These results connect the dynamical signaling of p53 with the microscopic repair process and create a new paradigm for the interplay of complex dynamics and phase transitions in biology.

Stochastic enzyme kinetics and the quasi-steady-state reductions: Application of the slow scale linear noise approximation à la Fenichel

The linear noise approximation models the random fluctuations from the mean-field model of a chemical reaction that unfolds near the thermodynamic limit. Specifically, the fluctuations obey a linear Langevin equation up to order [Formula: see text], where [Formula: see text] is the size of the chemical system (usually the volume). In the presence of disparate timescales, the linear noise approximation admits a quasi-steady-state reduction referred to as the slow scale linear noise approximation (ssLNA). Curiously, the ssLNAs reported in the literature are slightly different. The differences in the reported ssLNAs lie at the mathematical heart of the derivation. In this work, we derive the ssLNA directly from geometric singular perturbation theory and explain the origin of the different ssLNAs in the literature. Moreover, we discuss the loss of normal hyperbolicity and we extend the ssLNA derived from geometric singular perturbation theory to a non-classical singularly perturbed problem. In so doing, we disprove a commonly-accepted qualifier for the validity of stochastic quasi-steady-state approximation of the Michaelis -Menten reaction mechanism.

Nonmodular oscillator and switch based on RNA decay drive regeneration of multimodal gene expression

Periodic gene expression dynamics are key to cell and organism physiology. Studies of oscillatory expression have focused on networks with intuitive regulatory negative feedback loops, leaving unknown whether other common biochemical reactions can produce oscillations. Oscillation and noise have been proposed to support mammalian progenitor cells’ capacity to restore heterogenous, multimodal expression from extreme subpopulations, but underlying networks and specific roles of noise remained elusive. We use mass-action-based models to show that regulated RNA degradation involving as few as two RNA species—applicable to nearly half of human protein-coding genes—can generate sustained oscillations without explicit feedback. Diverging oscillation periods synergize with noise to robustly restore cell populations’ bimodal expression on timescales of days. The global bifurcation organizing this divergence relies on an oscillator and bistable switch which cannot be decomposed into two structural modules. Our work reveals surprisingly rich dynamics of post-transcriptional reactions and a potentially widespread mechanism underlying development, tissue regeneration, and cancer cell heterogeneity.

On the Validity of the Stochastic Quasi-Steady-State Approximation in Open Enzyme Catalyzed Reactions: Timescale Separation or Singular Perturbation?

The quasi-steady-state approximation is widely used to develop simplified deterministic or stochastic models of enzyme catalyzed reactions. In deterministic models, the quasi-steady-state approximation can be mathematically justified from singular perturbation theory. For several closed enzymatic reactions, the homologous extension of the quasi-steady-state approximation to the stochastic regime, known as the stochastic quasi-steady-state approximation, has been shown to be accurate under the analogous conditions that permit the quasi-steady-state reduction in the deterministic counterpart. However, it was recently demonstrated that the extension of the stochastic quasi-steady-state approximation to an open Michaelis-Menten reaction mechanism is only valid under a condition that is far more restrictive than the qualifier that ensures the validity of its corresponding deterministic quasi-steady-state approximation. In this paper, we suggest a possible explanation for this discrepancy from the lens of geometric singular perturbation theory. In so doing, we illustrate a misconception in the application of the quasi-steady-state approximation: timescale separation does not imply singular perturbation.

Universally valid reduction of multiscale stochastic biochemical systems using simple non-elementary propensities

Biochemical systems consist of numerous elementary reactions governed by the law of mass action. However, experimentally characterizing all the elementary reactions is nearly impossible. Thus, over a century, their deterministic models that typically contain rapid reversible bindings have been simplified with non-elementary reaction functions (e.g., Michaelis-Menten and Morrison equations). Although the non-elementary reaction functions are derived by applying the quasi-steady-state approximation (QSSA) to deterministic systems, they have also been widely used to derive propensities for stochastic simulations due to computational efficiency and simplicity. However, the validity condition for this heuristic approach has not been identified even for the reversible binding between molecules, such as protein-DNA, enzyme-substrate, and receptor-ligand, which is the basis for living cells. Here, we find that the non-elementary propensities based on the deterministic total QSSA can accurately capture the stochastic dynamics of the reversible binding in general. However, serious errors occur when reactant molecules with similar levels tightly bind, unlike deterministic systems. In that case, the non-elementary propensities distort the stochastic dynamics of a bistable switch in the cell cycle and an oscillator in the circadian clock. Accordingly, we derive alternative non-elementary propensities with the stochastic low-state QSSA, developed in this study. This provides a universally valid framework for simplifying multiscale stochastic biochemical systems with rapid reversible bindings, critical for efficient stochastic simulations of cell signaling and gene regulation. To facilitate the framework, we provide a user-friendly open-source computational package, ASSISTER, that automatically performs the present framework.

Three-dimensional stochastic simulation of chemoattractant-mediated excitability in cells

During the last decade, a consensus has emerged that the stochastic triggering of an excitable system drives pseudopod formation and subsequent migration of amoeboid cells. The presence of chemoattractant stimuli alters the threshold for triggering this activity and can bias the direction of migration. Though noise plays an important role in these behaviors, mathematical models have typically ignored its origin and merely introduced it as an external signal into a series of reaction-diffusion equations. Here we consider a more realistic description based on a reaction-diffusion master equation formalism to implement these networks. In this scheme, noise arises naturally from a stochastic description of the various reaction and diffusion terms. Working on a three-dimensional geometry in which separate compartments are divided into a tetrahedral mesh, we implement a modular description of the system, consisting of G-protein coupled receptor signaling (GPCR), a local excitation-global inhibition mechanism (LEGI), and signal transduction excitable network (STEN). Our models implement detailed biochemical descriptions whenever this information is available, such as in the GPCR and G-protein interactions. In contrast, where the biochemical entities are less certain, such as the LEGI mechanism, we consider various possible schemes and highlight the differences between them. Our simulations show that even when the LEGI mechanism displays perfect adaptation in terms of the mean level of proteins, the variance shows a dose-dependence. This differs between the various models considered, suggesting a possible means for determining experimentally among the various potential networks. Overall, our simulations recreate temporal and spatial patterns observed experimentally in both wild-type and perturbed cells, providing further evidence for the excitable system paradigm. Moreover, because of the overall importance and ubiquity of the modules we consider, including GPCR signaling and adaptation, our results will be of interest beyond the field of directed migration.

Though the term noise usually carries negative connotations, it can also contribute positively to the characteristic dynamics of a system. In biological systems, where noise arises from the stochastic interactions between molecules, its study is usually confined to genetic regulatory systems in which copy numbers are small and fluctuations large. However, noise can have important roles when the number of signaling molecules is large. The extension of pseudopods and the subsequent motion of amoeboid cells arises from the noise-induced trigger of an excitable system. Chemoattractant signals bias this triggering thereby directing cell motion. To date, this paradigm has not been tested by mathematical models that account accurately for the noise that arises in the corresponding reactions. In this study, we employ a reaction-diffusion master equation approach to investigate the effects of noise. Using a modular approach and a three-dimensional cell model with specific subdomains attributed to the cell membrane and cortex, we explore the spatiotemporal dynamics of the system. Our simulations recreate many experimentally-observed cell behaviors thereby supporting the biased-excitable network hypothesis.

Variance Reduction with Array-RQMC for Tau-Leaping Simulation of Stochastic Biological and Chemical Reaction Networks

We explore the use of Array-RQMC, a randomized quasi-Monte Carlo method designed for the simulation of Markov chains, to reduce the variance when simulating stochastic biological or chemical reaction networks with [Formula: see text]-leaping. The task is to estimate the expectation of a function of molecule copy numbers at a given future time T by the sample average over n sample paths, and the goal is to reduce the variance of this sample-average estimator. We find that when the method is properly applied, variance reductions by factors in the thousands can be obtained. These factors are much larger than those observed previously by other authors who tried RQMC methods for the same examples. Array-RQMC simulates an array of realizations of the Markov chain and requires a sorting function to reorder these chains according to their states, after each step. The choice of sorting function is a key ingredient for the efficiency of the method, although in our experiments, Array-RQMC was never worse than ordinary Monte Carlo, regardless of the sorting method. The expected number of reactions of each type per step also has an impact on the efficiency gain.

On rapid oscillations driving biological processes at disparate timescales

We consider a generic biological process described by a dynamical system, subject to an input signal with a high-frequency periodic component. The rapid oscillations of the input signal induce inherently multiscale dynamics, motivating order-reduction techniques. It is intuitive that the system behaviour is well approximated by its response to the averaged input signal. However, changes to the high-frequency component that preserve the average signal are beyond the reach of such intuitive reasoning. In this study, we explore system response under the influence of such an input signal by exploiting the timescale separation between high-frequency input variations and system response time. Employing the asymptotic method of multiple scales, we establish that, in some circumstances, the intuitive approach is simply the leading-order asymptotic contribution. We focus on higher-order corrections that capture the response to the details of the high-frequency component beyond its average. This approach achieves a reduction in system complexity while providing valuable insight into the structure of the response to the oscillations. We develop the general theory for nonlinear systems, while highlighting the important case of systems affine in the state and the input signal, presenting examples of both discrete and continuum state spaces. Importantly, this class of systems encompasses biochemical reaction networks described by the chemical master equation and its continuum approximations. Finally, we apply the framework to a nonlinear system describing mRNA translation and protein expression previously studied in the literature. The analysis shines new light on several aspects of the system quantification and both extends and simplifies results previously obtained.

Dynamic bistable switches enhance robustness and accuracy of cell cycle transitions

Bistability is a common mechanism to ensure robust and irreversible cell cycle transitions. Whenever biological parameters or external conditions change such that a threshold is crossed, the system abruptly switches between different cell cycle states. Experimental studies have uncovered mechanisms that can make the shape of the bistable response curve change dynamically in time. Here, we show how such a dynamically changing bistable switch can provide a cell with better control over the timing of cell cycle transitions. Moreover, cell cycle oscillations built on bistable switches are more robust when the bistability is modulated in time. Our results are not specific to cell cycle models and may apply to other bistable systems in which the bistable response curve is time-dependent.

Many systems in nature show bistability, which means they can evolve to one of two stable steady states under exactly the same conditions. Which state they evolve to depends on where the system comes from. Such bistability underlies the switching behavior that is essential for cells to progress in the cell division cycle. A quick switch happens when the cell jumps from one steady state to another steady state. Typical of this switching behavior is its robustness and irreversibility. In this paper, we expand this viewpoint of the dynamics of the cell cycle by considering bistable switches which themselves are changing in time. This gives the cell an extra layer of control over transitions both in time and in space, and can make those transitions more robust. Such dynamically changing bistability can appear very naturally. We show this in a model of mitotic entry, in which we include a nuclear and cytoplasmic compartment. The activity of a crucial cell cycle protein follows a bistable switch in each compartment, but the shape of its response is changing in time as proteins are imported into and exported from the nucleus.

Stochastic time-dependent enzyme kinetics: Closed-form solution and transient bimodality

We derive an approximate closed-form solution to the chemical master equation describing the Michaelis-Menten reaction mechanism of enzyme action. In particular, assuming that the probability of a complex dissociating into an enzyme and substrate is significantly larger than the probability of a product formation event, we obtain expressions for the time-dependent marginal probability distributions of the number of substrate and enzyme molecules. For delta function initial conditions, we show that the substrate distribution is either unimodal at all times or else becomes bimodal at intermediate times. This transient bimodality, which has no deterministic counterpart, manifests when the initial number of substrate molecules is much larger than the total number of enzyme molecules and if the frequency of enzyme-substrate binding events is large enough. Furthermore, we show that our closed-form solution is different from the solution of the chemical master equation reduced by means of the widely used discrete stochastic Michaelis-Menten approximation, where the propensity for substrate decay has a hyperbolic dependence on the number of substrate molecules. The differences arise because the latter does not take into account enzyme number fluctuations, while our approach includes them. We confirm by means of a stochastic simulation of all the elementary reaction steps in the Michaelis-Menten mechanism that our closed-form solution is accurate over a larger region of parameter space than that obtained using the discrete stochastic Michaelis-Menten approximation.

Misuse of the Michaelis-Menten rate law for protein interaction networks and its remedy

For over a century, the Michaelis-Menten (MM) rate law has been used to describe the rates of enzyme-catalyzed reactions and gene expression. Despite the ubiquity of the MM rate law, it accurately captures the dynamics of underlying biochemical reactions only so long as it is applied under the right condition, namely, that the substrate is in large excess over the enzyme-substrate complex. Unfortunately, in circumstances where its validity condition is not satisfied, especially so in protein interaction networks, the MM rate law has frequently been misused. In this review, we illustrate how inappropriate use of the MM rate law distorts the dynamics of the system, provides mistaken estimates of parameter values, and makes false predictions of dynamical features such as ultrasensitivity, bistability, and oscillations. We describe how these problems can be resolved with a slightly modified form of the MM rate law, based on the total quasi-steady state approximation (tQSSA). Furthermore, we show that the tQSSA can be used for accurate stochastic simulations at a lower computational cost than using the full set of mass-action rate laws. This review describes how to use quasi-steady state approximations in the right context, to prevent drawing erroneous conclusions from in silico simulations.

Basal leakage in oscillation: Coupled transcriptional and translational control using feed-forward loops

The circadian clock is a complex system that plays many important roles in most organisms. Previously, many mathematical models have been used to sharpen our understanding of the Arabidopsis clock, which brought to light the roles of each transcriptional and post-translational regulations. However, the presence of both regulations, instead of either transcription or post-translation, raised curiosity of whether the combination of these two regulations is important for the clock’s system. In this study, we built a series of simplified oscillators with different regulations to study the importance of post-translational regulation (specifically, 26S proteasome degradation) in the clock system. We found that a simple transcriptional-based oscillator can already generate sustained oscillation, but the oscillation can be easily destroyed in the presence of transcriptional leakage. Coupling post-translational control with transcriptional-based oscillator in a feed-forward loop will greatly improve the robustness of the oscillator in the presence of basal leakage. Using these general models, we were able to replicate the increased variability observed in the E3 ligase mutant for both plant and mammalian clocks. With this insight, we also predict a plausible regulator of several E3 ligase genes in the plant’s clock. Thus, our results provide insights into and the plausible importance in coupling transcription and post-translation controls in the clock system.

For circadian clocks, several current models had successfully captured the essential dynamic behavior of the clock system mainly with transcriptional regulation. Previous studies have shown that the 26S proteasome degradation controls are important in maintaining the stability of circadian rhythms. However, how the loss-of-function or over-expression mutant of this targeted degradations lead to unstable oscillation is still unclear. In this work, we investigate the importance of coupled transcriptional and post-translational feedback loop in the circadian oscillator. With general models our study indicate that the unstable behavior of degradation mutants could be caused by the increase in the basal level of the clock genes. We found that coupling a non-linear degradation control into this transcriptional based oscillator using feed-forward loop improves the robustness of the oscillator. Using this finding, we further predict some plausible regulators of Arabidopsis’s E3 ligase protein such as COP1 and SINAT5. Hence, our results provide insights on the importance of coupling transcription and post-translation controls in the clock system.

Nonlinear phenomena in models of the circadian clock

The mammalian circadian clock is well-known to be important for our sleep-wake cycles, as well as other daily rhythms such as temperature regulation, hormone release or feeding-fasting cycles. Under normal conditions, these daily cyclic events follow 24 h limit cycle oscillations, but under some circumstances, more complex nonlinear phenomena, such as the emergence of chaos, or the splitting of physiological dynamics into oscillations with two different periods, can be observed. These nonlinear events have been described at the organismic and tissue level, but whether they occur at the cellular level is still unknown. Our results show that period-doubling, chaos and splitting appear in different models of the mammalian circadian clock with interlocked feedback loops and in the absence of external forcing. We find that changes in the degradation of clock genes and proteins greatly alter the dynamics of the system and can induce complex nonlinear events. Our findings highlight the role of degradation rates in determining the oscillatory behaviour of clock components, and can contribute to the understanding of molecular mechanisms of circadian dysregulation.

Stochastic Modeling of Autoregulatory Genetic Feedback Loops: A Review and Comparative Study

Autoregulatory feedback loops are one of the most common network motifs. A wide variety of stochastic models have been constructed to understand how the fluctuations in protein numbers in these loops are influenced by the kinetic parameters of the main biochemical steps. These models differ according to 1) which subcellular processes are explicitly modeled, 2) the modeling methodology employed (discrete, continuous, or hybrid), and 3) whether they can be analytically solved for the steady-state distribution of protein numbers. We discuss the assumptions and properties of the main models in the literature, summarize our current understanding of the relationship between them, and highlight some of the insights gained through modeling.

Revisiting the Reduction of Stochastic Models of Genetic Feedback Loops with Fast Promoter Switching

Propensity functions of the Hill type are commonly used to model transcriptional regulation in stochastic models of gene expression. This leads to an effective reduced master equation for the mRNA and protein dynamics only. Based on deterministic considerations, it is often stated or tacitly assumed that such models are valid in the limit of rapid promoter switching. Here, starting from the chemical master equation describing promoter-protein interactions, mRNA transcription, protein translation, and decay, we prove that in the limit of fast promoter switching, the distribution of protein numbers is different than that given by standard stochastic models with Hill-type propensities. We show the differences are pronounced whenever the protein-DNA binding rate is much larger than the unbinding rate, a special case of fast promoter switching. Furthermore, we show using both theory and simulations that use of the standard stochastic models leads to drastically incorrect predictions for the switching properties of positive feedback loops and that these differences decrease with increasing mean protein burst size. Our results confirm that commonly used stochastic models of gene regulatory networks are only accurate in a subset of the parameter space consistent with rapid promoter switching.

A new Michaelis-Menten equation valid everywhere multi-scale dynamics prevails

The Michaelis-Menten reaction scheme is among the most influential models in the field of biochemistry, since it led to a very popular expression for the rate of an enzyme-catalysed reaction. After the realisation that this expression is valid in a limited region of the parameter space, two additional expressions were later introduced. The range of validity of these three expressions has been studied thoroughly, since the significance of a reliable rate is not based only on the accuracy of its predictive abilities but also on the physical insight that is acquired in the process of its construction. Here a new expression for the rate is introduced that is valid in practically the full parameter space and reduces to the expressions in the literature, when considering appropriate limits. The new expression is produced by employing algorithmic tools for asymptotic analysis, so that its construction is not hindered by the dimensional form and the complexity of the full model or by a wide parameter range of interest. These tools can be employed for the derivation of enzyme rate expressions of much more complex kinetics mechanisms.

A Robust Model for Circadian Redox Oscillations

The circadian clock is an endogenous oscillator that controls daily rhythms in metabolism, physiology, and behavior. Although the timekeeping components differ among species, a common design principle is a transcription-translation negative feedback loop. However, it is becoming clear that other mechanisms can contribute to the generation of 24 h rhythms. Peroxiredoxins (Prxs) exhibit 24 h rhythms in their redox state in all kingdoms of life. In mammalian adrenal gland, heart and brown adipose tissue, such rhythms are generated as a result of an inactivating hyperoxidation reaction that is reduced by coordinated import of sulfiredoxin (Srx) into the mitochondria. However, a quantitative description of the Prx/Srx oscillating system is still missing. We investigate the basic principles that generate mitochondrial Prx/Srx rhythms using computational modeling. We observe that the previously described delay in mitochondrial Srx import, in combination with an appropriate separation of fast and slow reactions, is sufficient to generate robust self-sustained relaxation-like oscillations. We find that our conceptual model can be regarded as a series of three consecutive phases and two temporal switches, highlighting the importance of delayed negative feedback and switches in the generation of oscillations.

Limits of noise for autoregulated gene expression

Gene expression is influenced by extrinsic noise (involving a fluctuating environment of cellular processes) and intrinsic noise (referring to fluctuations within a cell under constant environment). We study the standard model of gene expression including an (in-)active gene, mRNA and protein. Gene expression is regulated in the sense that the protein feeds back and either represses (negative feedback) or enhances (positive feedback) its production at the stage of transcription. While it is well-known that negative (positive) feedback reduces (increases) intrinsic noise, we give a precise result on the resulting fluctuations in protein numbers. The technique we use is an extension of the Langevin approximation and is an application of a central limit theorem under stochastic averaging for Markov jump processes (Kang et al. in Ann Appl Probab 24:721-759, 2014). We find that (under our scaling and in equilibrium), negative feedback leads to a reduction in the Fano factor of at most 2, while the noise under positive feedback is potentially unbounded. The fit with simulations is very good and improves on known approximations.

Ensemble methods for stochastic networks with special reference to the biological clock of Neurospora crassa

A major challenge in systems biology is to infer the parameters of regulatory networks that operate in a noisy environment, such as in a single cell. In a stochastic regime it is hard to distinguish noise from the real signal and to infer the noise contribution to the dynamical behavior. When the genetic network displays oscillatory dynamics, it is even harder to infer the parameters that produce the oscillations. To address this issue we introduce a new estimation method built on a combination of stochastic simulations, mass action kinetics and ensemble network simulations in which we match the average periodogram and phase of the model to that of the data. The method is relatively fast (compared to Metropolis-Hastings Monte Carlo Methods), easy to parallelize, applicable to large oscillatory networks and large (~2000 cells) single cell expression data sets, and it quantifies the noise impact on the observed dynamics. Standard errors of estimated rate coefficients are typically two orders of magnitude smaller than the mean from single cell experiments with on the order of ~1000 cells. We also provide a method to assess the goodness of fit of the stochastic network using the Hilbert phase of single cells. An analysis of phase departures from the null model with no communication between cells is consistent with a hypothesis of Stochastic Resonance describing single cell oscillators. Stochastic Resonance provides a physical mechanism whereby intracellular noise plays a positive role in establishing oscillatory behavior, but may require model parameters, such as rate coefficients, that differ substantially from those extracted at the macroscopic level from measurements on populations of millions of communicating, synchronized cells.

Beyond the Michaelis-Menten equation: Accurate and efficient estimation of enzyme kinetic parameters

Examining enzyme kinetics is critical for understanding cellular systems and for using enzymes in industry. The Michaelis-Menten equation has been widely used for over a century to estimate the enzyme kinetic parameters from reaction progress curves of substrates, which is known as the progress curve assay. However, this canonical approach works in limited conditions, such as when there is a large excess of substrate over enzyme. Even when this condition is satisfied, the identifiability of parameters is not always guaranteed, and often not verifiable in practice. To overcome such limitations of the canonical approach for the progress curve assay, here we propose a Bayesian approach based on an equation derived with the total quasi-steady-state approximation. In contrast to the canonical approach, estimates obtained with this proposed approach exhibit little bias for any combination of enzyme and substrate concentrations. Importantly, unlike the canonical approach, an optimal experiment to identify parameters with certainty can be easily designed without any prior information. Indeed, with this proposed design, the kinetic parameters of diverse enzymes with disparate catalytic efficiencies, such as chymotrypsin, fumarase, and urease, can be accurately and precisely estimated from a minimal amount of timecourse data. A publicly accessible computational package performing such accurate and efficient Bayesian inference for enzyme kinetics is provided.

Stability of Wake-Sleep Cycles Requires Robust Degradation of the PERIOD Protein

Robustness in biology is the stability of phenotype under diverse genetic and/or environmental perturbations. The circadian clock has remarkable stability of period and phase that-unlike other biological oscillators-is maintained over a wide range of conditions. Here, we show that the high fidelity of the circadian system stems from robust degradation of the clock protein PERIOD. We show that PERIOD degradation is regulated by a balance between ubiquitination and deubiquitination, and that disruption of this balance can destabilize the clock. In mice with a loss-of-function mutation of the E3 ligase gene β-Trcp2, the balance of PERIOD degradation is perturbed and the clock becomes dramatically unstable, presenting a unique behavioral phenotype unlike other circadian mutant animal models. We believe that our data provide a molecular explanation for how circadian phases, such as wake-sleep onset times, can become unstable in humans, and we present a unique mouse model to study human circadian disorders with unstable circadian rhythm phases.

Reduction of multiscale stochastic biochemical reaction networks using exact moment derivation

Biochemical reaction networks (BRNs) in a cell frequently consist of reactions with disparate timescales. The stochastic simulations of such multiscale BRNs are prohibitively slow due to high computational cost for the simulations of fast reactions. One way to resolve this problem uses the fact that fast species regulated by fast reactions quickly equilibrate to their stationary distribution while slow species are unlikely to be changed. Thus, on a slow timescale, fast species can be replaced by their quasi-steady state (QSS): their stationary conditional expectation values for given slow species. As the QSS are determined solely by the state of slow species, such replacement leads to a reduced model, where fast species are eliminated. However, it is challenging to derive the QSS in the presence of nonlinear reactions. While various approximation schemes for the QSS have been developed, they often lead to considerable errors. Here, we propose two classes of multiscale BRNs which can be reduced by deriving an exact QSS rather than approximations. Specifically, if fast species constitute either a feedforward network or a complex balanced network, the reduced model based on the exact QSS can be derived. Such BRNs are frequently observed in a cell as the feedforward network is one of fundamental motifs of gene or protein regulatory networks. Furthermore, complex balanced networks also include various types of fast reversible bindings such as bindings between transcriptional factors and gene regulatory sites. The reduced models based on exact QSS, which can be calculated by the computational packages provided in this work, accurately approximate the slow scale dynamics of the original full model with much lower computational cost.

Interrogating the topological robustness of gene regulatory circuits by randomization

One of the most important roles of cells is performing their cellular tasks properly for survival. Cells usually achieve robust functionality, for example, cell-fate decision-making and signal transduction, through multiple layers of regulation involving many genes. Despite the combinatorial complexity of gene regulation, its quantitative behavior has been typically studied on the basis of experimentally verified core gene regulatory circuitry, composed of a small set of important elements. It is still unclear how such a core circuit operates in the presence of many other regulatory molecules and in a crowded and noisy cellular environment. Here we report a new computational method, named random circuit perturbation (RACIPE), for interrogating the robust dynamical behavior of a gene regulatory circuit even without accurate measurements of circuit kinetic parameters. RACIPE generates an ensemble of random kinetic models corresponding to a fixed circuit topology, and utilizes statistical tools to identify generic properties of the circuit. By applying RACIPE to simple toggle-switch-like motifs, we observed that the stable states of all models converge to experimentally observed gene state clusters even when the parameters are strongly perturbed. RACIPE was further applied to a proposed 22-gene network of the Epithelial-to-Mesenchymal Transition (EMT), from which we identified four experimentally observed gene states, including the states that are associated with two different types of hybrid Epithelial/Mesenchymal phenotypes. Our results suggest that dynamics of a gene circuit is mainly determined by its topology, not by detailed circuit parameters. Our work provides a theoretical foundation for circuit-based systems biology modeling. We anticipate RACIPE to be a powerful tool to predict and decode circuit design principles in an unbiased manner, and to quantitatively evaluate the robustness and heterogeneity of gene expression.

Cells are able to robustly carry out their essential biological functions, possibly because of multiple layers of tight regulation via complex, yet well-designed, gene regulatory networks involving a substantial number of genes. State-of-the-art genomics technology has enabled the mapping of these large gene networks, yet it remains a tremendous challenge to elucidate their design principles and the regulatory mechanisms underlying their biological functions such as signal processing and decision-making. One of the key barriers is the absence of accurate kinetics for the regulatory interactions, especially from in vivo experiments. To this end, we have developed a new computational modeling method, Random Circuit Perturbation (RACIPE), to explore the dynamic behaviors of gene regulatory circuits without the requirement of detailed kinetic parameters. RACIPE takes a network topology as the input, and generates an unbiased ensemble of models with varying kinetic parameters. Each model is subjected to simulation, followed by statistical analysis for the ensemble. We tested RACIPE on several gene circuits, and found that the predicted gene expression patterns from all of the models converge to experimentally observed gene state clusters. We expect RACIPE to be a powerful method to identify the role of network topology in determining network operating principles.

A thermodynamically consistent model of the post-translational Kai circadian clock

The principal pacemaker of the circadian clock of the cyanobacterium S. elongatus is a protein phosphorylation cycle consisting of three proteins, KaiA, KaiB and KaiC. KaiC forms a homohexamer, with each monomer consisting of two domains, CI and CII. Both domains can bind and hydrolyze ATP, but only the CII domain can be phosphorylated, at two residues, in a well-defined sequence. While this system has been studied extensively, how the clock is driven thermodynamically has remained elusive. Inspired by recent experimental observations and building on ideas from previous mathematical models, we present a new, thermodynamically consistent, statistical-mechanical model of the clock. At its heart are two main ideas: i) ATP hydrolysis in the CI domain provides the thermodynamic driving force for the clock, switching KaiC between an active conformational state in which its phosphorylation level tends to rise and an inactive one in which it tends to fall; ii) phosphorylation of the CII domain provides the timer for the hydrolysis in the CI domain. The model also naturally explains how KaiA, by acting as a nucleotide exchange factor, can stimulate phosphorylation of KaiC, and how the differential affinity of KaiA for the different KaiC phosphoforms generates the characteristic temporal order of KaiC phosphorylation. As the phosphorylation level in the CII domain rises, the release of ADP from CI slows down, making the inactive conformational state of KaiC more stable. In the inactive state, KaiC binds KaiB, which not only stabilizes this state further, but also leads to the sequestration of KaiA, and hence to KaiC dephosphorylation. Using a dedicated kinetic Monte Carlo algorithm, which makes it possible to efficiently simulate this system consisting of more than a billion reactions, we show that the model can describe a wealth of experimental data.

Circadian clocks are biological timekeeping devices with a rhythm of 24 hours in living cells pertaining to all kingdoms of life. They help organisms to coordinate their behavior with the day-night cycle. The circadian clock of the cyanobacterium Synechococcus elongatus is one of the simplest and best characterized clocks in biology. The central clock component is the protein KaiC, which is phosphorylated and dephosphorylated in a cyclical manner with a 24 hr period. While we know from elementary thermodynamics that oscillations require a net turnover of fuel molecules, in this case ATP, how ATP hydrolysis drives the clock has remained elusive. Based on recent experimental observations and building on ideas from existing models, we construct the most detailed mathematical model of this system to date. KaiC consists of two domains, CI and CII, which each can bind ATP, yet only CII can be phosphorylated. Moreover, KaiC can exist in two conformational states, an active one in which the phosphorylation level tends to rise, and an inactive one in which it tends to fall. Our model predicts that ATP hydrolysis in the CI domain is the principal energetic driver of the clock, driving the switching between the two conformational states, while phosphorylation in the CII domain provides the timer for the conformational switch. The coupling between ATP hydrolysis in the CI domain and phosphorylation in the CII domain leads to novel testable predictions.