Synchronization, clustering, and weak chimeras in a densely coupled transcription-based oscillator model for split circadian rhythms

Abstract

The synchronization dynamics for the circadian gene expression in the suprachiasmatic nucleus is investigated using a transcriptional circadian clock gene oscillator model. With global coupling in constant dark (DD) conditions, the model exhibits a one-cluster phase synchronized state, in dim light (dim LL), bistability between one- and two-cluster states and in bright LL, a two-cluster state. The two-cluster phase synchronized state, where some oscillator pairs synchronize in-phase, and some anti-phase, can explain the splitting of the circadian clock, i.e., generation of two bouts of daily activities with certain species, e.g., with hamsters. The one- and two-cluster states can be reached by transferring the animal from DD or bright LL to dim LL, i.e., the circadian synchrony has a memory effect. The stability of the one- and two-cluster states was interpreted analytically by extracting phase models from the ordinary differential equation models. In a modular network with two strongly coupled oscillator populations with weak intragroup coupling, with appropriate initial conditions, one group is synchronized to the one-cluster state and the other group to the two-cluster state, resulting in a weak-chimera state. Computational modeling suggests that the daily rhythms in sleep–wake depend on light intensity acting on bilateral networks of suprachiasmatic nucleus (SCN) oscillators. Addition of a network heterogeneity (coupling between the left and right SCN) allowed the system to exhibit chimera states. The simulations can guide experiments in the circadian rhythm research to explore the effect of light intensity on the complexities of circadian desynchronization.