Disruption in the regulation of casein kinase 2 in circadian rhythm leads to pathological states: cancer, diabetes and neurodegenerative disorders

Abstract

Introduction:

Circadian rhythm maintains the sleep–wake cycle in biological systems. Various biological activities are regulated and modulated by the circadian rhythm, disruption of which can result in onset of diseases. Robust rhythms of phosphorylation profiles and abundances of PERIOD (PER) proteins are thought to be the master keys that drive circadian clock functions. The role of casein kinase 2 (CK2) in circadian rhythm via its direct interactions with the PER protein has been extensively studied; however, the exact mechanism by which it affects circadian rhythms at the molecular level is not known.

Methods:

Here, we propose an extended circadian rhythm model in Drosophila that incorporates the crosstalk between the PER protein and CK2. We studied the regulatory role of CK2 in the dynamics of PER proteins involved in circadian rhythm using the stochastic simulation algorithm.

Results:

We observed that variations in the concentration of CK2 in the circadian rhythm model modulates the PER protein dynamics at different cellular states, namely, active, weakly active, and rhythmic death. These oscillatory states may correspond to distinct pathological cellular states of the living system. We find molecular noise at the expression level of CK2 to switch normal circadian rhythm to any of the three above-mentioned circadian oscillatory states. Our results suggest that the concentration levels of CK2 in the system has a strong impact on its dynamics, which is reflected in the time evolution of PER protein.

Discussion:

We believe that our findings can contribute towards understanding the molecular mechanisms of circadian dysregulation in pathways driven by the PER mutant genes and their pathological states, including cancer, obesity, diabetes, neurodegenerative disorders, and socio-psychological disease.

Introduction

Life is a manifestation of rhythms. One such rhythm is circadian, an endogenous biological process with an oscillation of an approximately 24-h cycle (Dalchau et al., 2018; Lane et al., 2023). The circadian rhythm is present in a wide range of living systems, such as plants, animals, fungi, and cyanobacteria (Goldbeter, 2002; Gonze et al., 2002; Albrecht, 2010; Lane et al., 2023). Although these rhythms are self-sustained and endogenous, they can be entrained to the local environment by external factors (zeitgebers), such as daylight, temperature, and molecular fluctuations (Zeng et al., 1996; Yu and Hardin, 2006). These environmental and metabolic stimuli (e.g., dietary intake) regulate biological processes, such as the sleep–wake cycle, energy metabolism, hormone and immunological responses, and cell proliferation (Koronowski, 2021; Lee, 2021).

A panoply of researchers has reported that changes in the rhythmic process may lead to different diseases (Zeng et al., 1996; Foster et al., 2013; Tan et al., 2018; Mokhlesi et al., 2019). Onset of diseases, such as mood and sleep disorders, cancer, obesity, and diabetes, are all strongly linked with disturbances in these rhythms caused by factors, such as chronic jet lag, eating late at night, sleep deprivation, variations in sunlight, and modifications in hormone regulation (including progesterone and testosterone) (Amaral et al., 2014; Lee, 2021). The expression and activity of various oncogenes and tumor suppressors, in both tumor tissues and the host, are extensively altered by environmental and genetic disruptions in circadian rhythms. Such alterations lead to the incidence and progression of cancer (Papagiannakopoulos et al., 2016; Lee et al., 2019). Circadian disturbances can influence the immunological and metabolic functions of the host, favoring permissive tumor microenvironments in different types of cancer (Aiello et al., 2020; Hadadi et al., 2020). It is further known that disruption of the molecular clock in skeletal muscle promotes insulin resistance and obesity (Dyar et al., 2014). Although it is known that disrupted circadian rhythms affect metabolism, its impact on patients with type 2 diabetes (T2D) is not well-understood.

Circadian rhythm is caused by a genetic regulatory negative feedback loop that involves several clock genes and proteins in the biochemical reaction model (Gerstner and Yin, 2010; Duong et al., 2011). An essential enzyme known to be at the heart of self-sustaining circadian clocks in fungi, plants, and animals is casein kinase 2 (CK2) (Akten et al., 2003; Filhol and Cochet, 2009). CK2 is a ubiquitous eukaryotic protein kinase present in both the nucleus and cytoplasm (Litchfield, 2003; Tsuchiya et al., 2009). CK2 contributes to a wide variety of physiological functions, complex cellular processes (including DNA repair and cell cycle control), and regulation of cell viability (Lin et al., 2002; Ueda et al., 2005). It destabilizes and phosphorylates the TIMELESS (TIM) and PERIOD (PER) proteins in Drosophila, which subsequently suppress the transcriptional activity of the CLOCK (CLK) gene (Mizoguchi et al., 2006; Szabó et al., 2013). The CLK activator is directly targeted by CK2. It is commonly found in tetrameric complexes formed from catalytic subunits (α and/or α’) and two regulatory β subunits and is traditionally categorized as a messenger-independent protein threonine/serine kinase (Lu et al., 2011; Lee and Kim, 2014). Furthermore, CK2 phosphorylates TIM and PER in vitro, indicating that it has an impact on these proteins (Edery et al., 1994; Kloss et al., 2001; Miyazaki et al., 2004). It has also been observed that CK2 phosphorylates TIM protein to a lesser extent than PER. These findings support the notion that CK2 directly regulates TIM and PER. Despite the several reports that have been published regarding the role of CK2 in circadian rhythm of Drosophila, the exact mechanism by which it regulates circadian rhythms at the molecular level is still unclear. Moreover, the dynamics of clock proteins modulated by CK2 is still an ongoing research question and needs to be systematically studied.

In the present study, we used a stochastic approach to examine the dynamic behavior of circadian rhythm driven by CK2 and thereby to discern the regulatory role of CK2 at a molecular level. We extended the biochemical pathway model for circadian rhythm in Drosophila by incorporating various possible interactions of CK2 with clock proteins. We present the method to simulate the developed biochemical network model and we discuss the numerical simulation results and interpretations.

Materials and methods

The extended model of the Drosophila circadian pathway, which incorporates the impact of CK2 and describes the stochastic simulation algorithm used for the simulation of biochemical reaction network, is proposed in the following sub-sections.

Description of the circadian CK2 model

Prompted by published experimental reports, which indicate that CK2 protein interacts with clock proteins in circadian rhythm (Yu et al., 2006; Szabó et al., 2013), we incorporated the molecular signaling pathways of CK2 in our model. The proposed schematic circadian–CK2 integrative model (Figure 1) is an extension of the circadian rhythm model in Drosophila, reported by Gonze et al. (2002) and Albrecht (2010). This model is based on the inhibition of a nuclear clock protein () at the level of transcription of its gene into mRNA () (Albrecht, 2010). mRNA is synthesized in the nucleus and is transferred to the cytosol, where it accumulates at a maximum rate of and gets ubiquitinated by enzymes (E1, E2, E3, and E4) with a rate of . In our model, the spatial temporal fluctuations in the concentrations of various forms of the nuclear clock protein () or cytosolic regulatory protein ( and ) is governed by the biochemical pathways. The rate of synthesis of protein is proportional to the formation of and is characterized by an apparent first order rate constant . Parameters and denote the maximum rate(s) and Michaelis constant(s) of the phosphatase and kinase involved in the reversible phosphorylation of into and into , respectively. The fully phosphorylated state is degraded by CK2 and is transported into the nucleus at a rate characterized by the apparent first-order rate constant Transport of the into the cytosol is characterized by the apparent first-order rate constant . It is reported that the negative feedback is exerted by the nuclear clock protein on gene transcription. CK2 phosphorylates the clock protein, which may form three complexes, namely , and due to three different forms of the available proteins. Synthesis of is assumed to occur with a rate constant of and subsequently, the dissociation of this complex is assumed to occur with a rate constant of . In a similar manner, it is hypothesized that complex formation of occurs with a rate constant of and its dissociation follows with a rate constant of . Further, complex formation of CK2 − P2 is considered to occur with a rate constant of and subsequently, the dissociation of this complex is assumed to occur with rate constants of . CK2 is ubiquitinated at a rate constant of . Synthesis of CK2 in the network is assumed to occur at the rate constant of . CK2 plays a significant role in various processes of the organism such as biological clocks. Table 1 provides the list of clock proteins linked to the integrated model and Table 2 presents the list of the biochemical reactions and propensity functions (probability of reaction). The rate constants associated with the proposed model are presented in Table 3.

Figure 1

Technique for simulation of the biochemical reaction pathway

Complex dynamical processes, governed by a set of well-defined reaction channels, are generally noise-induced stochastic processes due to random molecular interactions in the system (origin of intrinsic noise) and continuous interaction of the system with random environmental fluctuations (origin of extrinsic noise) (Doob, 1942; Gillespie, 1977; Miyazaki et al., 2004). The system with clock protein variables, whose population vector is defined by, , where is the transpose of the vector, which undergo reaction channels, is given by,

where is the set of classical rate constants. The stochastic rate constant of any of its reactions can be expressed in terms of classical rate constant as , where V is the system size and is the state change parameter of its reaction (Gillespie, 1977; Gillespie, 2000). The trajectories of the variables provided by birth and death processes due to molecular interaction given by the Equation 1 can be traced by solving the master equation, for the rate of change of configurational probability as a function of time. The master equation is denoted as

where and are the transition probabilities of the two configurational states and { ′}. Solving the master Equation 2 analytically for complex biological process is extremely difficult. However, numerical solutions for the master Equation 2 can be derived using stochastic simulation algorithm (SSA) (Gillespie, 1977; Gillespie, 2000) that is based on the theoretical foundations developed by Doob (1942) and initially proposed by Kendall (1950). The stochastic simulation implements a Monte Carlo algorithm that provides the exact numerical solution by considering every possible interaction in the system (Gillespie, 1977; Singh et al., 2018). This algorithm is in deed a non-spatial, individual-based analog of the master Equation 2, which is constructed on the physical basis of molecular collision in each reaction channel at a certain constant temperature.

This SSA is based on two crucial independent random processes, namely reaction fire and reaction time. These two processes are implemented in this algorithm by generating two statistically independent random numbers, namely and such that the reaction time is computed using , where , is the propensity function given by , where is the number of possible molecular combinations of reaction, and the reaction will fire when it satisfies . Intrinsic noise () associated with the clock protein dynamics in the system is inversely proportional to the square root of the systems size (i.e., (Gillespie, 2000; Singh et al., 2018; Sharma et al., 2019; Singh et al., 2021).

Algorithm to calculate permutation entropy: the Bandt and Pompe approach

Permutation entropy can be used to measure the complexity of a system associated with the dynamics of the system’s variables (Bandt and Pompe, 2002). The basic algorithm for calculating permutation entropy ) of a time series is as follows:

Consider a dynamical variable of a system given by the discrete time series , where is the finite total number of discrete time elements in the time series data. We define an embedding dimension ), preferably to represent the data in consecutive patterns of the size dimension . For a particular value of , there are possible permuted sequences of inequalities of sequence elements. If we take , we will have arrangements of permutations given by,

, , , , , , where, Then, we calculate the probability for each single permutation sequence, which is the ratio of number of values for a particular sequence of permutations to the total number of all possible permutations for the embedding dimension, in the time series data. Shannon entropy for a sequence of perturbations can then be calculated by Equation 3,

and permutation entropy of embedded dimension (d) by the sum of these entropies as mentioned in Equation 4,

where . The mapped permutation entropy spectrum of time series is indicated by , and exhibits a behavior similar to that of Lyapunov spectrum of the same time series (Bodine et al., 2001).

Results

We performed large-scale numerical simulations of the proposed circadian–CK2 integrated model developed using stochastic simulation algorithm (SSA) (Gillespie, 1977). We demonstrate the SSA with a primary focus on the regulatory role of CK2 in the dynamics of clock proteins involved in circadian rhythm. We also demonstrate the impact of molecular noise owing to molecular integration of environment within the system (Goldbeter et al., 2017).

Modulation of clock proteins (, ) by CK2

We first present the results of how regulates the dynamics of clock protein (, ). Since the population of the protein in the system of size (at fixed is proportional to the rate of synthesis of this protein in the system (), variation in might cause changes in the interaction rate of other clock proteins in the reaction network model. Hence, we looked for variations in the dynamics of the clock proteins, driven by via changes in the values for . First, we allowed all three phosphorylation events in the cytosolic PER proteins (well-known clock proteins), with to occur (Figure 1), and found that for small values of , prominent oscillation (active state) in the nuclear PER protein is exhibited with time period of oscillation, h and amplitude of oscillation (Figure 2A upper panel). Further increase in suppressed the oscillation ( allowing increased and decreased significantly, which could correspond to weak circadian activity (Crino, 2011). This increase in the time period was due to increased phosphorylation owing to an increase in interaction of CK2 with the PER protein (), which was evident from the experimental reports of Cao et al. (2009). Further, weak circadian activity may cause various diseases, such as aging of the brain, metabolic dysfunction, dementia, and cancer-related disorders (Rao et al., 2002). When is sufficiently large (), the dynamics showed both oscillation and amplitude death scenarios (Rao et al., 2002). This state may correspond to circadian rhythm death (Rao et al., 2002). These three important states in the dynamics of of the circadian rhythm induced by CK2 can be shown in two-dimensional space of parameter as ( and ) (Figure 2B). Here, the parameters and are mean values of amplitude and time period for a range of hours. In the phase diagram presented in Figure 2B, all three circadian states are seen to be clearly demarcated.

Figure 2

We delineated the measure of complexity in the three derived states (active, weak activity, and rhythmic death) driven by by calculating permutation entropy ( of the dynamics of at the three respective states (Figure 2C). Upon considering the average permutation entropies of active, weak activity, and rhythmic death as and the simulation results demonstrated that . The three circadian rhythm states (active, weak activity, and rhythmic death) could easily be classified using (Figure 2D). In this phase diagram, each point is the average of permutation spectrum for each value of with error bars. Hence, the results indicated that could be used as a parameter to classify various states of the circadian rhythm.

The dynamics of exhibited similarity to that of at fixed system size and variable concentration of CK2. The dynamics of for various values of for constant values of system size is presented in Figure 3A. We measured the permutation entropies of MP at the three states under the situation of fixed system size and of variations of the CK2 concentration. In the case of permutation entropies of MP at the three states (Figure 3B), we made observations similar to that of (as depicted in Figure 2B).

Figure 3

We propose that and can be useful for clinicians and medical practitioners to classify the various states of circadian rhythm. These results suggest that CK2 plays a very dynamic role in the regulation of circadian rhythm in living systems. The change in rhythmic properties may alter the physiological processes of the organism leading to various diseases.

Impact of CK2 configurational interactions with the per proteins on circadian rhythm

Figure 4 presents the various possible interactions of CK2 with each of the PER proteins , and to form complexes , , and . All simulations were performed for the same values of and . used in the previous simulation depicted in Figure 3. The results demonstrated that if interacts individually with any one of the PER proteins via formation of a complex, dynamics exhibits both active (for smaller values of ) and weak activity (for larger values of ) states (Figure 4A). However, the system needs significantly large values of to obtain the rhythmic death state (data not shown). Further, configurational interaction of with is much more sensitive in driving the circadian rhythm states, as compared to those with or . We then calculated and for each time series when CK2 was allowed to interact with each one of , and in the circadian–CK2 model (Figure 4B). Each curve in and (Figure 4B) illustrates that the three circadian rhythm states could be observed. We also found that interaction of CK2 with is more sensitive than that with either or because the three circadian states are found at smaller values of in the case of .

Figure 4

The three models, namely model1, model2, and model3, as presented in Figure 4, are patterns of weak activity when CK2 systematically interacts with any one of the PER proteins (); these three patterns can also be considered as pathological states of the corresponding configuration. This rationale is acceptable because the patterns of CK2 activities (namely amplitude and time period) are significantly changed due to drastic increase in phosphorylation of PER protein(s) with CK2 (Cao et al., 2009). This drastic change in circadian rhythm may cause various diseases ranging from socio-psychological diseases (Wulff et al., 2010) and metabolic syndrome (Maury et al., 2010), to various types of cancer (Savvidis and Koutsilieris, 2012). Excess phosphorylation of CK2 with PER protein switches the pathological state to rhythmic death pattern which could be apoptosis signature.

We then calculated permutation entropy ( of the time series of as a function of for individual interactions of CK2 with (Figure 4C). We found that the three circadian states can be detected distinctly for each time series, as in the previous simulation, and can be used as a parameter to detect these states which can be useful to clinicians and medical practitioners (Figure 4D).

The dynamics of MP (Figure 5A) exhibited similarity to that of (as seen in Figure 4A). We also measured the permutation entropies of MP at fixed CK2 concentration and variation in system size. The observations from permutation entropies of MP at the three dynamical states of MP (Figure 5B) were similar to those of .

Figure 5

Noise can regulate the circadian states

Noise is an inherent property associated with the dynamics of any natural system (Glass, 2001; Bhadana et al., 2019, 2021). We studied the dynamics of for four possible configurational interactions of CK2 with the cytosolic PER proteins as we did in the previous simulation by keeping fixed at (which is associated with regular pattern of active circadian rhythm), and changing the strength of the noise with four values for as and (Figure 6). In all the configurational interactions of with the PER protein, we observed that an increase in noise (i.e., decrease in the value of ) can alter the dynamics in each of the three circadian states. This indicates that noise can regulate the dynamics of circadian rhythm.

Figure 6

Generally, the size of an animal cell (of both the proliferation and non-proliferation types) can change due to various reasons (Lloyd, 2013), such as cell cycle progression and cell growth (Conlon and Raff, 2003), pathological states in muscle cells and starvation (Bodine et al., 2001), tumor or cancer progression (Lum et al., 2005), defects in synaptic wiring/rewiring in neurons (Crino, 2011), and manipulation of extracellular signals to prevent apoptosis (Lum et al., 2005). Animal cells can have size variability of up to of the normal cell size (Lum et al., 2005), and this change in the cell size can drastically affect molecular crowding within the cell (Feig et al., 2017). This variation in the system’s size is reflected in the dynamics of the system’s variables as internal noise fluctuation () (Gillespie, 2000). Since the circadian rhythm system is governed by well-defined reaction channels (Table 2), this change in molecular crowding could cause two significant impacts on the system. Firstly, it allows changes in the rate of interaction of the clock proteins in the system. This change in molecular interaction leads to changes in the internal noise associated with the system reflected in the dynamics of the constituting variables and of other clock proteins in the system. Hence, excess noise may destroy the signal associated with the system variables and may lead to external collapse of the system itself. Secondly, this molecular crowding may trigger changes in molecular traffic in the system; traffic jam due to excess molecular crowding may make the system unable to perform normal functions. In such a scenario, the amplitude is minimized, and the system faces death with infinitely large time period (Figure 6). Similar behavior can be observed in the other three configurational interaction models of CK2 with each of and particularly, in the interaction of CK2 with wherein the noise was slightly more sensitive (Figure 6A).

Results of calculations to derive the amplitudes and time periods of dynamics as a function of for four different configurational interaction of with and are presented in Figure 6B. The results indicated that noise, measured by could drive the system to the three distinct states, namely active, weak activity, and rhythmic death. Hence, it can be inferred that noise is an important parameter that can trigger the system at various rhythmic states and can regulate the dynamics of the system. Permutation entropy spectra of corresponding to all the four models as a function of values are presented in Figure 6C. We delineated the measure of complexity in three different derived states driven by by calculating permutation entropy of the dynamics of for all possible configurational interaction of with cytosolic PER proteins (, and . Further, it can be observed that the measure of as a function of can distinctly classify the three different circadian rhythm states (Figure 6C). Upon denoting and as permutation entropies corresponding to the circadian states of active, weak activity, and rhythmic death, it can be inferred from the results that (Figure 6D).

A similar behavior can be seen in the case of (Figure 7A) and (Figure 7B) dynamics. The transition from sustained oscillation fluctuations (active state) to amplitude death (rhythmic death) via the weak activity state can be clearly seen in the two-dimensional plot (, ) as presented in Figure 7C.

Figure 7

Circadian per-induced cellular pathways that can trigger pathological states

We then identified cellular pathways that could be triggered by PER mutants (P₀, P_1, and P₂) in the circadian rhythm model (Figure 8). These cellular pathways can be grouped into three categories, cancer pathways (Savvidis and Koutsilieris, 2012; Papagiannakopoulos et al., 2016; Lee et al., 2019), socio–psychological pathways (Foster et al., 2013), and metabolic pathways (Maury et al., 2010). PER genes with their own rhythms generally regulate expressions of other genes participating in various important cellular pathways, significant changes of which may lead to various diseases. Since disruptions in the dynamics of and in circadian rhythm may lead to drastic changes in the dynamics and mechanisms of these closely interacting pathways and to various pathological states, they can cause various diseases. Dysfunctional PER proteins disturb key biological functions, such as cell proliferation, DNA damage, cell cycle, and apoptosis, resulting in various type of cancers (Deibel et al., 2015).

Figure 8

According to various studies, there may exist a tissue-dependent relationship between insulin resistance (or type 2 diabetes, T2D) and dysregulated molecular clock activity (Hart, 2013; Gabriel et al., 2021). In a time-course experiment (Jakubowicz et al., 2017) with white adipose tissue biopsies from people with healthy weight or obesity or T2D, no variation in the rhythm and amplitude of core-clock (PER1, PER2, PER3, DBP, BMAL1, and CRY2), metabolic (PGC1), and clock-related (REVERB) genes could be observed in biopsy tissues (Jakubowicz et al., 2017). On the other hand, when the sleep–wake cycle and dietary factors were controlled, it was seen that the amplitude oscillations of core-clock genes and a number of rhythmic genes are decreased in adipose tissues from patients with T2D, as compared to lean and healthy individuals (Stenvers et al., 2019). The mRNA expression of BMAL1, PER1, PER2, and PER3 in leukocytes was observed to be lower in non-diabetic individuals, as compared to those with diabetes (Hart, 2013; Stenvers et al., 2019). Additionally, the expression levels of PER1, PER3, and BMAL1 molecular-clock genes in leukocytes obtained from patients with T2D was seen inversely correlated with hemoglobin A1C (HbA1c) levels, indicating an association between insulin resistance and T2D (Gabriel et al., 2021). Expression of PER2, PER3, and CRY2 mRNA were also significantly associated with plasma HbA1c levels and islet insulin content in pancreatic islets from healthy people and from those with T2D (Gabriel et al., 2021). Uncertainties still exist regarding the fundamental mechanisms that regulate metabolic rhythmicity, especially regarding whether rhythmicity is lost in T2D. Mutations in clock genes were initially associated with glucose homeostasis (Rudic et al., 2004) and later with hyperinsulinemia, hyperglycaemia, and obesity in murine models (Turek et al., 2005). Numerous single nucleotide polymorphisms (SNPs) in human clock genes, including rs18012602, rs4580704, rs4864584, rs3749474 and rs1464490, have also been associated with obesity (Sookoian et al., 2008), hyperglycaemia, and a higher prevalence of T2D (Sookoian et al., 2008).

Aberrations in the PER protein disrupts multiple biological functions, such as age-related hydroxy methylation, which causes increased risk for socio–psychological diseases (e.g., learning and memory impairment, and Alzheimer’s disease) (Liu and Chang, 2017). Significant changes in the dynamics of and in the circadian rhythm and noise in the system may lead to drastic changes in the dynamics and mechanisms of pathways and to various pathological states, causing diseases corresponding to these pathways. Therefore, it is crucial to identify pathways that are significantly affected by changes in the and dynamics.

On the other hand, patients having any one or more of the diseases mentioned in Figure 8 can become healthy by maintaining a proper circadian rhythm (Foster et al., 2013), and the process could possibly be termed as a reversal of pathological state. However, even though this may be possible for some diseases, one needs to categorically study these possibilities from both the mathematical modeling and experimental points of view. In accordance with these postulations, previous reports have also recommended studies on genes/proteins in the circadian rhythm pathways as potential drug targets (Wood et al., 2009).

Discussion

Circadian rhythm is one of the most important biological rhythms that can regulate and interfere with various biological processes, such as those involved in cell regeneration, hormone production, and controlling brain activity. In the present work, we studied the impact of and stochastic noise on a circadian rhythm model. It has been reported in literature that CK2 promotes the progressive phosphorylation of clock protein that leads to the rapid degradation of hyperphosphorylated isoforms by the ubiquitin–proteasome pathway. Our results suggest that the presence of CK2 in the system has a strong impact on its dynamics, as reflected in the time evolution of the nuclear PER protein and mRNA. We found that CK2 drives three distinct circadian rhythm states, namely active, weak activity, and rhythmic death. The active state corresponds to the time period of oscillation () in the rhythm variable ( etc) as approximately 22–24 h, with optimal amplitude (A). When is larger than that of the active state, it is termed as weak activity, whereas the circadian state corresponding to , where is known as rhythmic death (Konopka and Benzer, 1971; van Soest et al., 2020).

Noise is an inherent property of any natural system. An interesting and important role of noise in a system is its capability to control the behavior of the system. The size of an animal cell can be changed due to various intracellular and extracellular factors (Lloyd, 2013); this variation in size is reflected in internal noise fluctuation in the system’s dynamics. In our study, we observed that noise can trigger the three circadian states and can control the behavior of the system. Our observations on the impact of noise (due to defects in biochemical reactions) on the circadian states support the hypothesis that the circadian clock is highly sensitive to CK2 activity.

Radical changes in circadian rhythm can occur due to various reasons, such as genetic defects, irregular work shifts, and aging (Jakubowicz et al., 2017). Change in the circadian state leading to weak activity may result in various diseases, especially cancer. This is because variations in the PER gene (denoted as in drosophila and as in human) may affect disruption in the cell cycle via the cMyc pathway with PER, ATK pathway with PER1 and estrogen signal, and HER with PER mutants; these pathways are known to be signature for progression of cancer in organs such as breast (Konopka and Benzer, 1971). Further, since binding of PER to androgen receptor (AR) causes the inhibition of AR transcriptional activity, the disruption of circadian rhythm may cause prostate cancer (Cao et al., 2009). Apart from prostate cancer, defects in circadian rhythm may cause various other cancer types (Savvidis and Koutsilieris, 2012), such as colorectal cancer due to the PER2–ATM–Chk1/Chk2 pathway (Mazzoccoli et al., 2011).

Disruptions in circadian rhythm may also cause several other disorders, such as psychiatric and neurodegenerative diseases (Wulff et al., 2010), jet lag disorder, and mental illness (Foster et al., 2013), and illnesses due to aging (Deibel et al., 2015). This could be because circadian rhythm is associated with various important cellular pathways. Some of the circadian genes/proteins are found to preserve cellular stability; for instance, PER1 is experimentally found to be anti-apoptotic in nature. Jet lag refers to misalignment of body’s internal clock with the local time at the destination. The jet lag phenomenon often occurs when flying across two or more time zones. The symptoms include sleeping problems, impaired thinking, hampered physical function and stomach problems. In rare instances, jet lag leads to sleep paralysis and seizures. Jet lag, in the case of people taking frequent long-haul flights, can be a long-term problem. Chronic circadian rhythm disruption can raise the risk of chronic disorders such as diabetes, depression and cancer (Foster et al., 2013). Hence, one must maintain proper circadian rhythm in their day-to-day life.

Conclusion

We studied the model for circadian rhythm using stochastic simulation algorithm, and examined the behavior of the amplitude, time period and permutation entropy of PER proteins to identify three distinct circadian states, namely active, weak activity, and rhythmic death all driven by CK2 protein. The interaction between changes in the PER gene expression under certain conditions illustrate the need for mathematical models to understand the underlying processes. Although a full representation of biological systems is hard to achieve owing to modeling limitations, the present study might help to understand how complex oscillatory dynamics occur at the molecular level. However, experimental data are needed to validate such phenomena at molecular level of the circadian clock.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1

  AielloI.FedeleM. M.RománF.MarpeganL.CaldartC.ChiesaJ. J.et al. (2020). Circadian disruption promotes tumor-immune microenvironment remodeling favoring tumor cell proliferation. Sci. Adv.6:eaaz4530. doi: 10.1126/sciadv.aaz4530

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AktenB.JauchE.GenovaG. K.KimE. Y.EderyI.RaabeT.et al. (2003). A role for CK2 in the Drosophila circadian oscillator. Nat. Neurosci.6, 251–257. doi: 10.1038/nn1007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  AlbrechtU: The circadian clock, Vol. 12. New York, NY: Springer Science & Business Media. (2010).

  • Google Scholar
  • View reference in article

• 4

  AmaralF. G.CastrucciA.Cipolla-NetoJ.PoletiniM. O.MendezN.RichterH. G.et al. (2014). Environmental control of biological rhythms: effects on development, fertility and metabolism. J. Neuroendocrinol.26, 603–612. doi: 10.1111/jne.12144

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BandtC.PompeB. (2002). Permutation entropy: a natural complexity measure for time series. Phys. Rev. Lett.88:174102. doi: 10.1103/PhysRevLett.88.174102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BhadanaJ.ChanuA. L.MalikM. Z.SinghR. B. (2021). Noise and delay can shape distribution functions in stochastic reaction dynamics. Nonlinear Dyn.105, 797–811. doi: 10.1007/s11071-021-06643-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BhadanaJ.MalikM. Z.SinghR. B. (2019). Universality in stochastic enzymatic futile cycle. Appl. Math. Model.74, 658–667. doi: 10.1016/j.apm.2019.05.008

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  BodineS. C.StittT. N.GonzalezM.KlineW. O.StoverG. L.BauerleinR.et al. (2001). Akt/mTOR pathway is a crucial regulator of skeletal muscle hypertrophy and can prevent muscle atrophy in vivo. Nat. Cell Biol.3, 1014–1019. doi: 10.1038/ncb1101-1014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  CaoQ.GeryS.DashtiA.YinD.ZhouY.GuJ.et al. (2009). A role for the clock gene per1 in prostate cancer. Cancer Res.69, 7619–7625. doi: 10.1158/0008-5472.CAN-08-4199

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  ConlonI.RaffM. (2003). Differences in the way a mammalian cell and yeast cells coordinate cell growth and cell-cycle progression. J. Biol.2:7. doi: 10.1186/1475-4924-2-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  CrinoP. B. (2011). mTOR: a pathogenic signaling pathway in developmental brain malformations. Trends Mol. Med.17, 734–742. doi: 10.1016/j.molmed.2011.07.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  DalchauN.SzépG.Hernansaiz-BallesterosR.BarnesC. P.CardelliL.PhillipsA.et al. (2018). Computing with biological switches and clocks. Nat. Comput.17, 761–779. doi: 10.1007/s11047-018-9686-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  DeibelS. H.ZelinskiE. L.KeeleyR. J.KovalchukO.McDonaldR. J. (2015). Epigenetic alterations in the suprachiasmatic nucleus and hippocampus contribute to age-related cognitive decline. Oncotarget6, 23181–23203. doi: 10.18632/oncotarget.4036

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  DoobJ. L. (1942). Topics in the theory of Markoff chains. Trans. Am. Math. Soc.52, 37–64. doi: 10.1090/S0002-9947-1942-0006633-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  DuongH. A.RoblesM. S.KnuttiD.WeitzC. J. (2011). A molecular mechanism for circadian clock negative feedback. Science332, 1436–1439. doi: 10.1126/science.1196766

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  DyarK. A.CiciliotS.WrightL. E.BiensøR. S.TagliazucchiG. M.PatelV. R.et al. (2014). Muscle insulin sensitivity and glucose metabolism are controlled by the intrinsic muscle clock. Mol. Metabolism3, 29–41. doi: 10.1016/j.molmet.2013.10.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  EderyI.ZwiebelL. J.DembinskaM. E.RosbashM. (1994). Temporal phosphorylation of the Drosophila period protein. Proc. Natl. Acad. Sci. U. S. A.91, 2260–2264. doi: 10.1073/pnas.91.6.2260

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  FeigM.YuI.WangP.-H.NawrockiG.SugitaY. (2017). Crowding in cellular environments at an atomistic level from computer simulations. J. Phys. Chem. B121, 8009–8025. doi: 10.1021/acs.jpcb.7b03570

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  FilholO.CochetC. (2009). Cellular functions of protein kinase CK2: a dynamic affair. Cell. Mol. Life Sci.66, 1830–1839. doi: 10.1007/s00018-009-9151-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  FosterR. G.PeirsonS. N.WulffK.WinnebeckE.VetterC.RoennebergT. (2013). Sleep and circadian rhythm disruption in social jetlag and mental illness. Prog. Mol. Biol. Transl. Sci.119, 325–346. doi: 10.1016/B978-0-12-396971-2.00011-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  GabrielB. M.AltıntaşA.SmithJ. A.Sardon-PuigL.ZhangX.BasseA. L.et al. (2021). Disrupted circadian oscillations in type 2 diabetes are linked to altered rhythmic mitochondrial metabolism in skeletal muscle. Sci. Adv.7:eabi9654. doi: 10.1126/sciadv.abi9654

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  GerstnerJ. R.YinJ. C. (2010). Circadian rhythms and memory formation. Nat. Rev. Neurosci.11, 577–588. doi: 10.1038/nrn2881

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  GillespieD. T. (1977). Exact stochastic simulation of coupled chemical reactions. J. Phys. Chem.81, 2340–2361. doi: 10.1021/j100540a008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  GillespieD. T. (2000). The chemical Langevin equation. J. Chem. Phys.113, 297–306. doi: 10.1063/1.481811

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  GlassL. (2001). Synchronization and rhythmic processes in physiology. Nature410, 277–284. doi: 10.1038/35065745

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  GoldbeterA. (2002). Computational approaches to cellular rhythms. Nature420, 238–245. doi: 10.1038/nature01259

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  GoldbeterA. A.TaylorL.WakafZ.VasudevanS. R.FosterR. G. (2017). The genetics of circadian rhythms, sleep and health. Hum. Mol. Genet.26, R128–R138. doi: 10.1093/hmg/ddx240

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  GonzeD.HalloyJ.GoldbeterA. (2002). Deterministic versus stochastic models for circadian rhythms. J. Biol. Phys.28, 637–653. doi: 10.1023/A:1021286607354

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  HadadiE.TaylorW.LiX.-M.AslanY.VilloteM.RivièreJ.et al. (2020). Chronic circadian disruption modulates breast cancer stemness and immune microenvironment to drive metastasis in mice. Nat. Commun.11:3193. doi: 10.1038/s41467-020-16890-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  HartG. W. (2013). How sugar tunes your clock. Cell Metab.17, 155–156. doi: 10.1016/j.cmet.2013.01.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  JakubowiczD.WainsteinJ.LandauZ.RazI.AhrenB.ChapnikN.et al. (2017). Influences of breakfast on clock gene expression and postprandial glycemia in healthy individuals and individuals with diabetes: a randomized clinical trial. Diabetes Care40, 1573–1579. doi: 10.2337/dc16-2753

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  KendallD. (1950). An artificial realization of a simple “birth-and-death” process. J. R. Stat. Soc. Ser. B Methodol.12, 116–119. doi: 10.1111/j.2517-6161.1950.tb00048.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  KlossB.RothenfluhA.YoungM. W.SaezL. (2001). Phosphorylation of period is influenced by cycling physical associations of double-time, period, and timeless in the Drosophila clock. Neuron30, 699–706. doi: 10.1016/S0896-6273(01)00320-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  KonopkaR. J.BenzerS. (1971). Clock mutants of Drosophila melanogaster. Proc. Natl. Acad. Sci. U. S. A.68, 2112–2116. doi: 10.1073/pnas.68.9.2112

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  KoronowskiK. B. (2021). Sassone-Corsi P: communicating clocks shape circadian homeostasis. Science371:eabd0951. doi: 10.1126/science.abd0951

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  LaneJ. M.QianJ.MignotE.RedlineS.ScheerF. A.SaxenaR. (2023). Genetics of circadian rhythms and sleep in human health and disease. Nat. Rev. Genet.24, 4–20. doi: 10.1038/s41576-022-00519-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  LeeY. (2021). Roles of circadian clocks in cancer pathogenesis and treatment. Exp. Mol. Med.53, 1529–1538. doi: 10.1038/s12276-021-00681-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  LeeE.KimE. Y. (2014). A role for timely nuclear translocation of clock repressor proteins in setting circadian clock speed. Exp. Neurobiol.23, 191–199. doi: 10.5607/en.2014.23.3.191

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  LeeY.LahensN. F.ZhangS.BedontJ.FieldJ. M.SehgalA. (2019). G1/S cell cycle regulators mediate effects of circadian dysregulation on tumor growth and provide targets for timed anticancer treatment. PLoS Biol.17:e3000228. doi: 10.1371/journal.pbio.3000228

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  LinJ. M.KilmanV. L.KeeganK.PaddockB.Emery-LeM.RosbashM.et al. (2002). A role for casein kinase 2alpha in the Drosophila circadian clock. Nature420, 816–820. doi: 10.1038/nature01235

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  LitchfieldD. W. (2003). Protein kinase CK2: structure, regulation and role in cellular decisions of life and death. Biochem. J.369, 1–15. doi: 10.1042/bj20021469

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  LiuF.ChangH. C. (2017). Physiological links of circadian clock and biological clock of aging. Protein Cell8, 477–488. doi: 10.1007/s13238-016-0366-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  LloydA. C. (2013). The regulation of cell size. Cells154, 1194–1205. doi: 10.1016/j.cell.2013.08.053

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  LuS. X.LiuH.KnowlesS. M.LiJ.MaL.TobinE. M.et al. (2011). A role for protein kinase casein kinase2 α-subunits in the Arabidopsis circadian clock. Plant Physiol.157, 1537–1545. doi: 10.1104/pp.111.179846

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  LumJ. J.BauerD. E.KongM.HarrisM. H.LiC.LindstenT.et al. (2005). Growth factor regulation of autophagy and cell survival in the absence of apoptosis. Cells120, 237–248. doi: 10.1016/j.cell.2004.11.046

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  MauryE.RamseyK. M.BassJ. (2010). Circadian rhythms and metabolic syndrome: from experimental genetics to human disease. Circ. Res.106, 447–462. doi: 10.1161/CIRCRESAHA.109.208355

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  MazzoccoliG.PanzaA.ValvanoM. R.PalumboO.CarellaM.PazienzaV.et al. (2011). Clock gene expression levels and relationship with clinical and pathological features in colorectal cancer patients. Chronobiol. Int.28, 841–851. doi: 10.3109/07420528.2011.615182

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  MiyazakiK.NagaseT.MesakiM.NarukawaJ.OharaO.IshidaN. (2004). Phosphorylation of clock protein PER1 regulates its circadian degradation in normal human fibroblasts. Biochem. J.380, 95–103. doi: 10.1042/BJ20031308

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  MizoguchiT.PutterillJ.OhkoshiY. (2006). Kinase and phosphatase: the cog and spring of the circadian clock. Int. Rev. Cytol.250, 47–72. doi: 10.1016/S0074-7696(06)50002-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  MokhlesiB.TempleK. A.TjadenA. H.EdelsteinS. L.UtzschneiderK. M.NadeauK. J.et al. (2019). Association of self-reported sleep and circadian measures with glycemia in adults with prediabetes or recently diagnosed untreated type 2 diabetes. Diabetes Care42, 1326–1332. doi: 10.2337/dc19-0298

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  PapagiannakopoulosT.BauerM. R.DavidsonS. M.HeimannM.SubbarajL.BhutkarA.et al. (2016). Circadian rhythm disruption promotes lung tumorigenesis. Cell Metab.24, 324–331. doi: 10.1016/j.cmet.2016.07.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  RaoC. V.WolfD. M.ArkinA. P. (2002). Control, exploitation and tolerance of intracellular noise. Nature420, 231–237. doi: 10.1038/nature01258

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  RudicR. D.McNamaraP.CurtisA. M.BostonR. C.PandaS.HogeneschJ. B.et al. (2004). BMAL1 and CLOCK, two essential components of the circadian clock, are involved in glucose homeostasis. PLoS Biol.2:e377. doi: 10.1371/journal.pbio.0020377

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  SavvidisC.KoutsilierisM. (2012). Circadian rhythm disruption in cancer biology. Mol. Med.18, 1249–1260. doi: 10.2119/molmed.2012.00077

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  SharmaS. K.HaobijamD.SinghS. S.MalikM. Z.SinghR. B. (2019). Neuronal communication: stochastic neuron dynamics and multi-synchrony states. AEU Int. J. Electr. Commun.100, 75–85. doi: 10.1016/j.aeue.2019.01.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  SinghS. N.ChanuA. L.MalikM. Z.SinghR. B. (2021). Interplay of cellular states: role of delay as control mechanism. Phys A Stat. Mech. Appl.572:125869. doi: 10.1016/j.physa.2021.125869

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  SinghS. S.HaobijamD.MalikM. Z.IshratR.SinghR. B. (2018). Fractal rules in brain networks: signatures of self-organization. J. Theor. Biol.437, 58–66. doi: 10.1016/j.jtbi.2017.09.014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 58

  SookoianS.GemmaC.GianottiT. F.BurgueñoA.CastañoG.PirolaC. J. (2008). Genetic variants of clock transcription factor are associated with individual susceptibility to obesity. Am. J. Clin. Nutr.87, 1606–1615. doi: 10.1093/ajcn/87.6.1606

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 59

  StenversD. J.JongejanA.AtiqiS.VreijlingJ. P.LimonardE. J.EndertE.et al. (2019). Diurnal rhythms in the white adipose tissue transcriptome are disturbed in obese individuals with type 2 diabetes compared with lean control individuals. Diabetologia62, 704–716. doi: 10.1007/s00125-019-4813-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  SzabóÁ.PapinC.ZornD.PonienP.WeberF.RaabeT.et al. (2013). The CK2 kinase stabilizes CLOCK and represses its activity in the Drosophila circadian oscillator. PLoS Biol.11:e1001645. doi: 10.1371/journal.pbio.1001645

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 61

  TanX.ChapmanC. D.CedernaesJ.BenedictC. (2018). Association between long sleep duration and increased risk of obesity and type 2 diabetes: a review of possible mechanisms. Sleep Med. Rev.40, 127–134. doi: 10.1016/j.smrv.2017.11.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 62

  TsuchiyaY.AkashiM.MatsudaM.GotoK.MiyataY.NodeK.et al. (2009). Involvement of the protein kinase CK2 in the regulation of mammalian circadian rhythms. Sci. Signal.2:ra26. doi: 10.1126/scisignal.2000305

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  TurekF. W.JoshuC.KohsakaA.LinE.IvanovaG.McDearmonE.et al. (2005). Obesity and metabolic syndrome in circadian clock mutant mice. Science308, 1043–1045. doi: 10.1126/science.1108750

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 64

  UedaH. R.HayashiS.ChenW.SanoM.MachidaM.ShigeyoshiY.et al. (2005). System-level identification of transcriptional circuits underlying mammalian circadian clocks. Nat. Genet.37, 187–192. doi: 10.1038/ng1504

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 65

  van SoestI.Del OlmoM.SchmalC.HerzelH. (2020). Nonlinear phenomena in models of the circadian clock. J. R. Soc. Interface17:20200556. doi: 10.1098/rsif.2020.0556

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 66

  WoodP. A.YangX.HrusheskyW. J. (2009). Clock genes and cancer. Integr. Cancer Ther.8, 303–308. doi: 10.1177/1534735409355292

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 67

  WulffK.GattiS.WettsteinJ. G.FosterR. G. (2010). Sleep and circadian rhythm disruption in psychiatric and neurodegenerative disease. Nat. Rev. Neurosci.11, 589–599. doi: 10.1038/nrn2868

  • CrossRef
  • Google Scholar
  • View reference in article

• 68

  YuW.HardinP. E. (2006). Circadian oscillators of Drosophila and mammals. J. Cell Sci.119, 4793–4795. doi: 10.1242/jcs.03174

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 69

  YuW.ZhengH.HoulJ. H.DauwalderB.HardinP. E. (2006). PER-dependent rhythms in CLK phosphorylation and E-box binding regulate circadian transcription. Genes Dev.20, 723–733. doi: 10.1101/gad.1404406

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 70

  ZengH.QianZ.MyersM. P.RosbashM. (1996). A light-entrainment mechanism for the Drosophila circadian clock. Nature380, 129–135. doi: 10.1038/380129a0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article