Stochastic oscillations induced by intrinsic fluctuations in a self-repressing gene

Single-cell transcriptional dynamics in a living vertebrate

The ability to quantify transcriptional dynamics in individual cells via live imaging has revolutionized our understanding of gene regulation. However, such measurements are lacking in the context of vertebrate embryos. We addressed this deficit by applying MS2-MCP mRNA labeling to the quantification of transcription in zebrafish, a model vertebrate. We developed a platform of transgenic organisms, light sheet fluorescence microscopy, and optimized image analysis that enables visualization and quantification of MS2 reporters. We used these tools to obtain the first single-cell, real-time measurements of transcriptional dynamics of the segmentation clock. Our measurements challenge the traditional view of smooth clock oscillations and instead suggest a model of discrete transcriptional bursts that are organized in space and time. Together, these results highlight how measuring single-cell transcriptional activity can reveal unexpected features of gene regulation and how this data can fuel the dialogue between theory and experiment.

Quantitative control of noise in mammalian gene expression by dynamic histone regulation

Fluctuation ('noise') in gene expression is critical for mammalian cellular processes. Numerous mechanisms contribute to its origins, yet the mechanisms behind large fluctuations that are induced by single transcriptional activators remain elusive. Here, we probed putative mechanisms by studying the dynamic regulation of transcriptional activator binding, histone regulator inhibitors, chromatin accessibility, and levels of mRNAs and proteins in single cells. Using a light-induced expression system, we showed that the transcriptional activator could form an interplay with dual functional co-activator/histone acetyltransferases CBP/p300. This interplay resulted in substantial heterogeneity in H3K27ac, chromatin accessibility, and transcription. Simultaneous attenuation of CBP/p300 and HDAC4/5 reduced heterogeneity in the expression of endogenous genes, suggesting that this mechanism is universal. We further found that the noise was reduced by pulse-wide modulation of transcriptional activator binding possibly as a result of alternating the epigenetic states. Our findings suggest a mechanism for the modulation of noise in synthetic and endogenous gene expression systems.

Endogenous miRNA sponges mediate the generation of oscillatory dynamics for a non-coding RNA network

Oscillations are crucial to the normal function of living organisms, across a wide variety of biological processes. In eukaryotes, oscillatory dynamics are thought to arise from interactions at the protein and RNA levels; however, the role of non-coding RNA in regulating these dynamics remains understudied. In this work, we show how non-coding RNA acting as microRNA (miRNA) sponges in a conserved miRNA - transcription factor feedback motif, can give rise to oscillatory behaviour, and how to test for this experimentally. Control of these non-coding RNA can dynamically create oscillations or stability, and we show how this behaviour predisposes to oscillations in the stochastic limit. These results, supported by emerging evidence for the role of miRNA sponges in development, point towards key roles of different species of miRNA sponges, such as circular RNA, potentially in the maintenance of yet unexplained oscillatory behaviour. These results help to provide a paradigm for understanding functional differences between the many redundant, but distinct RNA species thought to act as miRNA sponges in nature, such as long non-coding RNA, pseudogenes, competing mRNA, circular RNA, and3' UTRs.

Oscillations in well-mixed, deterministic feedback systems: Beyond ring oscillators

•

I present a way of breaking down regulatory networks to find Hopf bifurcations. This helps find optimal conditions for oscillations in dynamical systems models of these networks.

•

In a model of negative auto-regulation of a gene by its dimeric protein, it is optimal for the monomer to degrade faster than the mRNA and the mRNA to degrade faster than the dimer.

•

Adding a weak positive feedback loop to a repressilator increases the probability of oscillations.

•

The optimal degradation rate of species in the sub-loop is higher than that of species outside it.

•

The opposite is true for a negative feedback sub-loop or a very strong positive feedback sub-loop.

A ring oscillator is a system in which one species regulates the next, which regulates the next and so on until the last species regulates the first. In addition, the number of the regulations which are negative, and so result in a reduction in the regulated species, is odd, making the overall feedback in the loop negative. In ring oscillators, the probability of oscillations is maximised if the degradation rates of the species are equal. When there is more than one loop in the regulatory network, the dynamics can be more complicated. Here, a systematic way of organising the characteristic equation of ODE models of regulatory networks is provided. This facilitates the identification of Hopf bifurcations. It is shown that the probability of oscillations in non-ring systems is maximised for unequal degradation rates. For example, when there is a ring and a second ring employing a subset of the genes in the first ring, then the probability of oscillations is maximised when the species in the sub-ring degrade more slowly than those outside, for a negative feedback subring. When the sub-ring forms a positive feedback loop, the optimal degradation rates are larger for the species in the sub-ring, provided the positive feedback is not too strong. By contrast, optimal degradation rates are smaller for the species in the sub-ring, when the positive feedback is very strong. Adding a positive feedback loop to a repressilator increases the probability of oscillations, provided the positive feedback is not too strong, whereas adding a negative feedback loop decreases the probability of oscillations. The work is illustrated with numerical simulations of example systems: an autoregulatory gene model in which transcription is downregulated by the protein dimer and three-species and four-species gene regulatory network examples.

Spatial Patterning from an Integrated Wnt/β-catenin and Notch/Delta Gene Circuit

Classically, the Wnt/β-catenin and Notch /Delta signaling pathways were thought to operate through separate mechanisms, performing distinct roles in tissue patterning. However, it has been shown that b-catenin activates transcription of Hesl, a signaling intermediate in the Notch /Delta pathway that controls its lateral inhibition mechanism. To investigate this non-canonical crosstalk mechanism, a new gene circuit, integrating the two pathways, is proposed and simulated in two-cell and multi-cell environments. This model also captures both Paneth cell- mediated and mesenchymal Wnt production. The simulations verify that the gene circuit is temporally bistable and capable of forming a pattern on a multi-cell grid. Last, the model exhibits a bifurcation based on the steady state concentration of Wnt and the relative amount of control b-catenin has over the Hesl promoter, providing a possible mechanism to explain why a homogeneous population of transit amplifying cells is observed directly above the more diverse stem niche.

Optimization-based synthesis of stochastic biocircuits with statistical specifications

Model-guided design has become a standard approach to engineering biomolecular circuits in synthetic biology. However, the stochastic nature of biomolecular reactions is often overlooked in the design process. As a result, cell-cell heterogeneity causes unexpected deviation of biocircuit behaviours from model predictions and requires additional iterations of design-build-test cycles. To enhance the design process of stochastic biocircuits, this paper presents a computational framework to systematically specify the level of intrinsic noise using well-defined metrics of statistics and design highly heterogeneous biocircuits based on the specifications. Specifically, we use descriptive statistics of population distributions as an intuitive specification language of stochastic biocircuits and develop an optimization-based computational tool that explores parameter configurations satisfying design requirements. Sensitivity analysis methods are also performed to ensure the robustness of a biocircuit design against extrinsic perturbations. These design tools are formulated with convex optimization programs to enable rigorous and efficient quantification of the statistics. We demonstrate these features by designing a stochastic negative feedback biocircuit that satisfies multiple statistical constraints and perform an in-depth study of noise propagation and regulation in negative feedback pathways.

Graph-facilitated resonant mode counting in stochastic interaction networks

Oscillations in dynamical systems are widely reported in multiple branches of applied mathematics. Critically, even a non-oscillatory deterministic system can produce cyclic trajectories when it is in a low copy number, stochastic regime. Common methods of finding parameter ranges for stochastically driven resonances, such as direct calculation, are cumbersome for any but the smallest networks. In this paper, we provide a systematic framework to efficiently determine the number of resonant modes and parameter ranges for stochastic oscillations relying on real root counting algorithms and graph theoretic methods. We argue that stochastic resonance is a network property by showing that resonant modes only depend on the squared Jacobian matrix J2, unlike deterministic oscillations which are determined by J By using graph theoretic tools, analysis of stochastic behaviour for larger interaction networks is simplified and stochastic dynamical systems with multiple resonant modes can be identified easily.

Multiple binding sites for transcriptional repressors can produce regular bursting and enhance noise suppression

Cells may control fluctuations in protein levels by means of negative autoregulation, where transcription factors bind DNA sites to repress their own production. Theoretical studies have assumed a single binding site for the repressor, while in most species it is found that multiple binding sites are arranged in clusters. We study a stochastic description of negative autoregulation with multiple binding sites for the repressor. We find that increasing the number of binding sites induces regular bursting of gene products. By tuning the threshold for repression, we show that multiple binding sites can also suppress fluctuations. Our results highlight possible roles for the presence of multiple binding sites of negative autoregulators.

Entraining Oscillations in the NF-κB Signaling System: With a Little Help from Noise

Two studies show that noise is a key ingredient of new mechanisms for entraining the NF-κB system.

Comparison of cellular oscillations driven by noise or deterministic mechanisms under cell-size scaling

Ultradian cycles are frequently observed in biological systems. They serve important roles in regulating, for example, cell fate and the development of the organism. Many mathematical models have been developed to analyze their behavior. Generally, these models can be classified into two classes: Deterministic models that generate oscillatory behavior by incorporating time delays or Hopf bifurcations, and stochastic models that generate oscillatory behavior by noise driven resonance. However, it is still unclear which of these two mechanisms applies to cellular oscillations. In this paper, we show through theoretical analysis and numerical simulation that we can distinguish which of these two mechanisms govern cellular oscillations, by measuring statistics of oscillation amplitudes for cells of different sizes. We found that, for oscillations driven deterministically, the normalized average amplitude is constant with respect to cell size, while the coefficient of variation of the amplitude scales with cell size with an exponent of -0.5. On the other hand, for oscillations driven stochastically, the coefficient of variation of the amplitude is constant with respect to cell size, while the normalized average amplitude scales with cell size with an exponent of -0.5. Our results provide a theoretical basis to discern whether a particular oscillatory behavior is governed by a deterministic or stochastic mechanism.

Heterodimer Autorepression Loop: A Robust and Flexible Pulse-Generating Genetic Module

We investigate the dynamics of the heterodimer autorepression loop (HAL), a small genetic module in which a protein A acts as an autorepressor and binds to a second protein B to form an AB dimer. For suitable values of the rate constants, the HAL produces pulses of A alternating with pulses of B. By means of analytical and numerical calculations, we show that the duration of A pulses is extremely robust against variation of the rate constants while the duration of the B pulses can be flexibly adjusted. The HAL is thus a minimal genetic module generating robust pulses with a tunable duration, an interesting property for cellular signaling.

Promoters Architecture-Based Mechanism for Noise-Induced Oscillations in a Single-Gene Circuit

It is well known that single-gene circuits with negative feedback loop can lead to oscillatory gene expression when they operate with time delay. In order to generate these oscillations many processes can contribute to properly timing such delay. Here we show that the time delay coming from the transitions between internal states of the cis-regulatory system (CRS) can drive sustained oscillations in an auto-repressive single-gene circuit operating in a small volume like a cell. We found that the cooperative binding of repressor molecules is not mandatory for a oscillatory behavior if there are enough binding sites in the CRS. These oscillations depend on an adequate balance between the CRS kinetic, and the synthesis/degradation rates of repressor molecules. This finding suggest that the multi-site CRS architecture can play a key role for oscillatory behavior of gene expression. Finally, our results can also help to synthetic biologists on the design of the promoters architecture for new genetic oscillatory circuits.

Mean field analysis of a spatial stochastic model of a gene regulatory network

A gene regulatory network may be defined as a collection of DNA segments which interact with each other indirectly through their RNA and protein products. Such a network is said to contain a negative feedback loop if its products inhibit gene transcription, and a positive feedback loop if a gene product promotes its own production. Negative feedback loops can create oscillations in mRNA and protein levels while positive feedback loops are primarily responsible for signal amplification. It is often the case in real biological systems that both negative and positive feedback loops operate in parameter regimes that result in low copy numbers of gene products. In this paper we investigate the spatio-temporal dynamics of a single feedback loop in a eukaryotic cell. We first develop a simplified spatial stochastic model of a canonical feedback system (either positive or negative). Using a Gillespie's algorithm, we compute sample trajectories and analyse their corresponding statistics. We then derive a system of equations that describe the spatio-temporal evolution of the stochastic means. Subsequently, we examine the spatially homogeneous case and compare the results of numerical simulations with the spatially explicit case. Finally, using a combination of steady-state analysis and data clustering techniques, we explore model behaviour across a subregion of the parameter space that is difficult to access experimentally and compare the parameter landscape of our spatio-temporal and spatially-homogeneous models.