Stochastic models for circadian rhythms: effect of molecular noise on periodic and chaotic behaviour

Circadian rhythms are endogenous oscillations that occur with a period close to 24 h in nearly all living organisms. These rhythms originate from the negative autoregulation of gene expression. Deterministic models based on such genetic regulatory processes account for the occurrence of circadian rhythms in constant environmental conditions (e.g., constant darkness), for entrainment of these rhythms by light-dark cycles, and for their phase-shifting by light pulses. When the numbers of protein and mRNA molecules involved in the oscillations are small, as may occur in cellular conditions, it becomes necessary to resort to stochastic simulations to assess the influence of molecular noise on circadian oscillations. We address the effect of molecular noise by considering the stochastic version of a deterministic model previously proposed for circadian oscillations of the PER and TIM proteins and their mRNAs in Drosophila. The model is based on repression of the per and tim genes by a complex between the PER and TIM proteins. Numerical simulations of the stochastic version of the model are performed by means of the Gillespie method. The predictions of the stochastic approach compare well with those of the deterministic model with respect both to sustained oscillations of the limit cycle type and to the influence of the proximity from a bifurcation point beyond which the system evolves to stable steady state. Stochastic simulations indicate that robust circadian oscillations can emerge at the cellular level even when the maximum numbers of mRNA and protein molecules involved in the oscillations are of the order of only a few tens or hundreds. The stochastic model also reproduces the evolution to a strange attractor in conditions where the deterministic PER-TIM model admits chaotic behaviour. The difference between periodic oscillations of the limit cycle type and aperiodic oscillations (i.e. chaos) persists in the presence of molecular noise, as shown by means of Poincaré sections. The progressive obliteration of periodicity observed as the number of molecules decreases can thus be distinguished from the aperiodicity originating from chaotic dynamics. As long as the numbers of molecules involved in the oscillations remain sufficiently large (of the order of a few tens or hundreds, or more), stochastic models therefore provide good agreement with the predictions of the deterministic model for circadian rhythms.     

Stable stochastic dynamics in yeast cell cycle

Chemical reactions in cells are subject to intense stochastic fluctuations. An important question is how the fundamental physiological behavior of the cell is kept stable against those noisy perturbations. In this study, a stochastic model of the cell cycle of budding yeast was constructed to analyze the effects of noise on the cell-cycle oscillation. The model predicts intense noise in levels of mRNAs and proteins, and the simulated protein levels explain the observed statistical tendency of noise in populations of synchronous and asynchronous cells. Despite intense noise in levels of proteins and mRNAs, the cell cycle is stable enough to bring the largely perturbed cells back to the physiological cyclic oscillation. The model shows that consecutively appearing fixed points are the origin of this stability of the cell cycle.     

Deterministic Versus Stochastic Models for Circadian Rhythms

Circadian rhythms which occur with a period close to 24 h in nearly all living organisms originate from the negative autoregulation of gene expression.Deterministic models based on genetic regulatory processes account for theoccurrence of circadian rhythms in constant environmental conditions (e.g.constant darkness), for entrainment of these rhythms by light-dark cycles, and for their phase-shifting by light pulses. At low numbers of protein and mRNA molecules, it becomes necessary to resort to stochastic simulations to assess the influence of molecular noise on circadian oscillations. We address the effect of molecular noise by considering two stochastic versions of a core model for circadian rhythms. The deterministic version of this core modelwas previously proposed for circadian oscillations of the PER protein in Drosophila and of the FRQ protein in Neurospora. In the first, non-developed version of the stochastic model, we introduce molecular noise without decomposing the deterministic mechanism into detailed reaction steps while in the second, developed version we carry out such a detailed decomposition. Numerical simulations of the two stochastic versions of the model are performed by means of the Gillespie method. We compare the predictions of the deterministic approach with those of the two stochastic models, with respect both to sustained oscillations of the limit cycle type and to the influence of the proximity of a bifurcation point beyond which the system evolves to a stable steady state. The results indicate that robust circadian oscillations can occur even when the numbers of mRNA and nuclear protein involved in the oscillatory mechanism are reduced to a few tens orhundreds, respectively. The non-developed and developed versions of the stochastic model yield largely similar results and provide good agreement with the predictions of the deterministic model for circadian rhythms.     

Effect of periodic stimulus on a neuronal diffusion model with signal-dependent noise

To relate the noise intensity with a periodically modulated input signal in a single neuron stochastic model we introduce a diffusion model with both time modulated drift and diffusion coefficient. Such a model is the continuous version of a Stein model with time oscillating frequencies for the Poisson processes describing the inputs impinging on the neuron. We focus here on some aspects of the resonance phenomenon for such a model. We compare the corresponding interspike interval distribution with the analogous distribution for a model sharing the same parameter values, but with constant noise intensity. Examples with two different levels for this noise intensity are discussed. The enhancement of the height of the peaks in the interspike interval distribution appearing at the modulation period, the improvement of the phase locking behavior and an enlargement of the noise ranges where a resonance like behavior arises are the main features observed in the considered cases.     

A stochastic limit cycle oscillator model of the EEG

We present an empirical model of the electroencephalogram (EEG) signal based on the construction of a stochastic limit cycle oscillator using Itô calculus. This formulation, where the noise influences actually interact with the dynamics, is substantially different from the usual definition of measurement noise. Analysis of model data is compared with actual EEG data using both traditional methods and modern techniques from nonlinear time series analysis. The model demonstrates visually displayed patterns and statistics that are similar to actual EEG data. In addition, the nonlinear mechanisms underlying the dynamics of the model do not manifest themselves in nonlinear time series analysis, paralleling the situation with real, non-pathological EEG data. This modeling exercise suggests that the EEG is optimally described by stochastic limit cycle behavior.     

Persistence and extinction of a stochastic single-specie model under regime switching in a polluted environment II

This is a continuation of our paper [Liu, M., Wang, K., 2010. Persistence and extinction of a stochastic single-specie model under regime switching in a polluted environment, J. Theor. Biol. 264, 934-944]. Taking both white noise and colored noise into account, a stochastic single-species model under regime switching in a polluted environment is studied. Sufficient conditions for extinction, stochastic nonpersistence in the mean, stochastic weak persistence and stochastic permanence are established. The threshold between stochastic weak persistence and extinction is obtained. The results show that a different type of noise has a different effect on the survival results.     

Stochastic and regulatory role of chromatin silencing in genomic response to environmental changes

Phenotypic diversity and fidelity can be balanced by controlling stochastic molecular mechanisms. Epigenetic silencing is one that has a critical role in stress response. Here we show that in yeast, incomplete silencing increases stochastic noise in gene expression, probably owing to unstable chromatin structure. Telomere position effect is suggested as one mechanism. Expression diversity in a population achieved in this way may render a subset of cells to readily respond to various acute stresses. By contrast, strong silencing tends to suppress noisy expression of genes, in particular those involved in life cycle control. In this regime, chromatin may act as a noise filter for precisely regulated responses to environmental signals that induce huge phenotypic changes such as a cell fate transition. These results propose modulation of chromatin stability as an important determinant of environmental adaptation and cellular differentiation.     

Engineered internal noise stochastic resonator in gene network: a model study

Based on a genetic bistable switch model coupled with a gene oscillator model, we have constructed a mesoscopic stochastic model for the coupled synthetic gene network, and studied how internal noise would influence the oscillation of such a system. We found that the state-to-state transitions can occur if the internal noise is taken into account, and the performance of resulting oscillation can reach a maximum in a certain internal noise level, which indicates the occurrence of internal noise stochastic resonance (SR) and makes the coupled gene network work as a stochastic resonator. The potential role of such an effect on gene expression systems is also discussed.     

Overcompensatory population dynamic responses to environmental stochasticity

1. To quantify the interactions between density-dependent, population regulation and density-independent limitation, we studied the time-series dynamics of an experimental laboratory insect microcosm system in which both environmental noise and resource limitation were manipulated. 2. A hierarchical Bayesian state-space approach is presented through which it is feasible to capture all sources of uncertainty, including observation error to accurately quantify the density dependence operating on the dynamics. 3. The regulatory processes underpinning the dynamics of two different bruchid beetles (Callosobruchus maculatus and Callosobruchus chinensis) are principally determined by environmental conditions, with fluctuations in abundance explained in terms of changes in overcompensatory dynamics and stochastic processes. 4. A general, stochastic population model is developed to explore the link between abundance fluctuations and the interaction between density dependence and noise. Taking account of time-lags in population regulation can substantially increase predicted population fluctuations resulting from underlying noise processes.     

Slow Noise in the Period of a Biological Oscillator Underlies Gradual Trends and Abrupt Transitions in Phasic Relationships in Hybrid Neural Networks

In order to study the ability of coupled neural oscillators to synchronize in the presence of intrinsic as opposed to synaptic noise, we constructed hybrid circuits consisting of one biological and one computational model neuron with reciprocal synaptic inhibition using the dynamic clamp. Uncoupled, both neurons fired periodic trains of action potentials. Most coupled circuits exhibited qualitative changes between one-to-one phase-locking with fairly constant phasic relationships and phase slipping with a constant progression in the phasic relationships across cycles. The phase resetting curve (PRC) and intrinsic periods were measured for both neurons, and used to construct a map of the firing intervals for both the coupled and externally forced (PRC measurement) conditions. For the coupled network, a stable fixed point of the map predicted phase locking, and its absence produced phase slipping. Repetitive application of the map was used to calibrate different noise models to simultaneously fit the noise level in the measurement of the PRC and the dynamics of the hybrid circuit experiments. Only a noise model that added history-dependent variability to the intrinsic period could fit both data sets with the same parameter values, as well as capture bifurcations in the fixed points of the map that cause switching between slipping and locking. We conclude that the biological neurons in our study have slowly-fluctuating stochastic dynamics that confer history dependence on the period. Theoretical results to date on the behavior of ensembles of noisy biological oscillators may require re-evaluation to account for transitions induced by slow noise dynamics.     

Diversity and noise effects in a model of homeostatic regulation of the sleep-wake cycle

Recent advances in sleep neurobiology have allowed development of physiologically based mathematical models of sleep regulation that account for the neuronal dynamics responsible for the regulation of sleep-wake cycles and allow detailed examination of the underlying mechanisms. Neuronal systems in general, and those involved in sleep regulation in particular, are noisy and heterogeneous by their nature. It has been shown in various systems that certain levels of noise and diversity can significantly improve signal encoding. However, these phenomena, especially the effects of diversity, are rarely considered in the models of sleep regulation. The present paper is focused on a neuron-based physiologically motivated model of sleep-wake cycles that proposes a novel mechanism of the homeostatic regulation of sleep based on the dynamics of a wake-promoting neuropeptide orexin. Here this model is generalized by the introduction of intrinsic diversity and noise in the orexin-producing neurons, in order to study the effect of their presence on the sleep-wake cycle. A simple quantitative measure of the quality of a sleep-wake cycle is introduced and used to systematically study the generalized model for different levels of noise and diversity. The model is shown to exhibit a clear diversity-induced resonance: that is, the best wake-sleep cycle turns out to correspond to an intermediate level of diversity at the synapses of the orexin-producing neurons. On the other hand, only a mild evidence of stochastic resonance is found, when the level of noise is varied. These results show that disorder, especially in the form of quenched diversity, can be a key-element for an efficient or optimal functioning of the homeostatic regulation of the sleep-wake cycle. Furthermore, this study provides an example of a constructive role of diversity in a neuronal system that can be extended beyond the system studied here.     

Investigating the two-moment characterisation of subcellular biochemical networks

While ordinary differential equations (ODEs) form the conceptual framework for modelling many cellular processes, specific situations demand stochastic models to capture the influence of noise. The most common formulation of stochastic models for biochemical networks is the chemical master equation (CME). While stochastic simulations are a practical way to realise the CME, analytical approximations offer more insight into the influence of noise. Towards that end, the two-moment approximation (2MA) is a promising addition to the established analytical approaches including the chemical Langevin equation (CLE) and the related linear noise approximation (LNA). The 2MA approach directly tracks the mean and (co)variance which are coupled in general. This coupling is not obvious in CME and CLE and ignored by LNA and conventional ODE models. We extend previous derivations of 2MA by allowing (a) non-elementary reactions and (b) relative concentrations. Often, several elementary reactions are approximated by a single step. Furthermore, practical situations often require the use of relative concentrations. We investigate the applicability of the 2MA approach to the well-established fission yeast cell cycle model. Our analytical model reproduces the clustering of cycle times observed in experiments. This is explained through multiple resettings of M-phase promoting factor (MPF), caused by the coupling between mean and (co)variance, near the G2/M transition.     

The Interplay of Intrinsic and Extrinsic Bounded Noises in Biomolecular Networks

After being considered as a nuisance to be filtered out, it became recently clear that biochemical noise plays a complex role, often fully functional, for a biomolecular network. The influence of intrinsic and extrinsic noises on biomolecular networks has intensively been investigated in last ten years, though contributions on the co-presence of both are sparse. Extrinsic noise is usually modeled as an unbounded white or colored gaussian stochastic process, even though realistic stochastic perturbations are clearly bounded. In this paper we consider Gillespie-like stochastic models of nonlinear networks, i.e. the intrinsic noise, where the model jump rates are affected by colored bounded extrinsic noises synthesized by a suitable biochemical state-dependent Langevin system. These systems are described by a master equation, and a simulation algorithm to analyze them is derived. This new modeling paradigm should enlarge the class of systems amenable at modeling. We investigated the influence of both amplitude and autocorrelation time of a extrinsic Sine-Wiener noise on: the Michaelis-Menten approximation of noisy enzymatic reactions, which we show to be applicable also in co-presence of both intrinsic and extrinsic noise, a model of enzymatic futile cycle and a genetic toggle switch. In and we show that the presence of a bounded extrinsic noise induces qualitative modifications in the probability densities of the involved chemicals, where new modes emerge, thus suggesting the possible functional role of bounded noises.     

Summation of Spatiotemporal Input Patterns in Leaky Integrate-and-Fire Neurons: Application to Neurons in the Cochlear Nucleus Receiving Converging Auditory Nerve Fiber Input

The response of leaky integrate-and-fire neurons is analyzed for periodic inputs whose phases vary with their spatial location. The model gives the relationship between the spatial summation distance and the degree of phase locking of the output spikes (i.e., locking to the periodic stochastic inputs, measured by the synchronization index). The synaptic inputs are modeled as an inhomogeneous Poisson process, and the analysis is carried out in the Gaussian approximation. The model has been applied to globular bushy cells of the cochlear nucleus, which receive converging inputs from auditory nerve fibers that originate at neighboring sites in the cochlea. The model elucidates the roles played by spatial summation and coincidence detection, showing how synchronization decreases with an increase in both frequency and spatial spread of inputs. It also shows under what conditions an enhancement of synchronization of the output relative to the input takes place.     

Monte Carlo simulation of biospecific interactions

Numerical simulations of the stochastic time evolution of biospecific interactions are described and show that when molecular populations are large, time course predictions match those obtained using a deterministic expression. When population size is decreased the effects of stochastic noise become apparent. The significance of stochastic noise in sensitive binding-based assay systems suggests an immediate need for models of this type.