A stochastic oscillator model simulates the entrainment of vertebrate cellular clocks by light

The circadian clock is a cellular mechanism that synchronizes various biological processes with respect to the time of the day. While much progress has been made characterizing the molecular mechanisms underlying this clock, it is less clear how external light cues influence the dynamics of the core clock mechanism and thereby entrain it with the light-dark cycle. Zebrafish-derived cell cultures possess clocks that are directly light-entrainable, thus providing an attractive laboratory model for circadian entrainment. Here, we have developed a stochastic oscillator model of the zebrafish circadian clock, which accounts for the core clock negative feedback loop, light input, and the proliferation of single-cell oscillator noise into population-level luminescence recordings. The model accurately predicts the entrainment dynamics observed in bioluminescent clock reporter assays upon exposure to a wide range of lighting conditions. Furthermore, we have applied the model to obtain refitted parameter sets for cell cultures exposed to a variety of pharmacological treatments and predict changes in single-cell oscillator parameters. Our work paves the way for model-based, large-scale screens for genetic or pharmacologically-induced modifications to the entrainment of circadian clock function. Author summaryThe circadian clock is a key, cell-autonomous timing mechanism that is encountered in most organisms. It is entrained by environmental lighting conditions and in turn temporally coordinates most aspects of physiology according to the time of day. Cell lines derived from zebrafish are attractive experimental models for studying how clocks are entrained by light since they possess clocks that respond directly to light stimuli. Here we describe a mathematical model for the behavior of the circadian clock in zebrafish cell lines during exposure to a range of lighting conditions. Using this model, we can determine how different pharmacological treatments may affect the entrainment dynamics of the clock and the degree of synchronization of individual cells circadian clocks in bioluminescent clock reporter assays. Our current model is mathematically simple and thus easy to apply and extend in future studies.