Stochastic simulation of a model for circadian rhythms in plants

Circadian clocks allow living organisms to anticipate and adapt to the daily variations of the environment. The interlocked feedback loops of the transcription factors network in the plant clock generate oscillations with expression peaks at specific times of the day. In this work, we explore the effect of molecular noise on the behavior of the plant circadian clock through numerical simulations. The influence of system size, photoperiod, and mutations of clock genes on the robustness of the oscillations are discussed. Our simulations show that the oscillations remain robust when the mRNA and protein levels are in the range of a few hundreds molecules. Entrainment by light-dark cycles enhances the robustness compared to constant conditions. Multiple light inputs and inter-cellular coupling also contribute to the robustness of the oscillations. The comparison between deterministic and stochastic simulations of single and double mutants shows that stochasticity does not qualitatively affect the behaviour of mutants but that they do not have the same robustness to noise. Finally, the model shows that noise can induce transitions between two limit cycles in a birhythmic clock mutant.