Introduction

All organisms, including microbes, plants, and animals, use an endogenous timekeeping system, namely the circadian clock (CC), which helps organisms to adapt to the 24-h cycle of the earth to control their daily physiology and behavior rhythms. Molecular pacemakers inside organisms drive the CC. In mammals, the coordination of essential behavioral, hormonal, and other physiological rhythms throughout the body relies on the CC1. It is also known that the circadian clock modulates cognitive activities2,3,4,5 and is linked to mood disorders6. In Drosophila, the CC controls the timing of eclosion and courtship, the period of rest and activity, and the timing of feeding; it also influences temperature preference7,8. Despite the differences in molecular components of the CCs, their features, organization, and the molecular mechanism that generate rhythmicity are very alike across organisms9.

Environmental cues such as light, food, and temperature can modulate the rhythmicity of the CC10. It is also known that the CC is susceptible to external magnetic fields (MFs). In the 1960s, Brown et al.11 found that small changes in the intensity of Earth’s MF synchronize the CCs of fiddler crabs and other organisms. Since then, the effects of external MF on the CC have been observed in multiple studies12,13,14,15,16,17,18,19,20,21,22,23,24. Similarly, Yoshii et al.25 have shown the effects of static MFs on the CC of Drosophila and found that exposure to these fields exhibited enhanced slowing of clock rhythms in the presence of blue light, with a maximal alteration at \(300 \,\mu \hbox {T}\), and reduced effects at both lower and slightly higher field strengths. However, the exact mechanism behind this phenomenon is still mostly unknown.

Additionally, a growing body of evidence points to the circadian cycles as a target for bipolar disorder treatments26,27. Bipolar disorder is correlated with disruptions in circadian rhythms26,27 and abnormalities in oxidative stress28,29,30,31,32,33,34,35. Lithium is the first-line treatment for bipolar disorders36,37, yet the exact mechanisms and pathways underlying this treatment are under debate. It has been shown that lithium treatment for hyperactivity in rats is isotope dependent38. Lithium has two stable isotopes, \({}^{6}{{\mathrm{Li}}}\) and \({}^{7}{{\mathrm{Li}}}\), which have different nuclear spin angular momentum, \(I_6=1\) and \(I_7=3/2\), respectively. In a recent study, it has also been proposed that lithium affects hyperactivity via the clock center channel in the brain39. This study further predicted that the magnetic field would influence the potency of lithium treatment. Furthermore, Dokucu et al.40 showed that in Drosophila lithium lengthens the period of the CC. Here, based on these findings, we propose a mechanism that can explain both the MF effects and lithium effects on Drosophila’s CC.

In the CC of Drosophila, the CLOCK (CLK) and CYCLE (CYC) transcription factors form a heterodimeric complex and promote the period (per) and timeless (tim) transcription mRNAs, which result in the assembly of the PERIOD (PER) and TIMELESS (TIM) proteins in the cytoplasm41, shown in Fig. 1a. During the night, PER and TIM accumulate and form a heterodimer. The TIM/PER complex enters the nucleus and then promotes the phosphorylation of CLK/CYC, which inhibits the promotion of the per and tim mRNAs. During the day, TIM and PER are gradually degraded, and consequently, CLK/CYC are released from repression to start a new cycle.

In this process, light activation of cryptochrome (CRY) protein is critical for the rhythmicity of the CC. CRYs regulate growth and development in plants; they also act as photo-receptors in some animal’s CC, where they are necessary components of the circadian clock42,43. In Drosophila’s CC, upon light absorption, CRY undergoes a conformational change that allows it to bind TIM7,10,44 which results in the degradation of TIM and hence resetting the clock, see Fig. 1a. In Drosophila, CRYs contain the flavin adenine dinucleotide (FAD) cofactor, which is the photoactive cofactor. Upon blue-light absorption, FAD can go through various redox states. In insects, this process produces the anionic semiquinone FAD, \(\hbox {FAD}^{{\varvec{\cdot }}-}\), and reactive oxygen species (ROS), which are thought to be the key signaling states for initiating TIM degradation45,46. In mammals, CRY’s do not bind FAD. However, CRY’s are essential for the development of intercellular networks in the suprachiasmatic nucleus (SCN), a circadian pacemaker in the brain, that subserve coherent rhythm expression; the network synchronizes cellular oscillators and corrects errors47.

It has been known for many years that migratory birds use Earth’s magnetic field for finding their way during migrations. Later, it was proposed that the radical pair mechanism (RPM) could be the key for birds’ magnetoreception48. Ritz et al.49 proposed that the candidate protein for such a mechanism could be the CRY in the retina of birds. Ever since, there have been extensive studies on that hypothesis, and to date, it is the most promising model for avian magnetoreception50,51,52 in birds, sharks, sea turtles, monarch butterflies, fruit flies, etc. These models are based on the study of the dynamics of the created pair of radicals which can be in a superposition of singlet (S) and triplet (T) states53, depending on the parent molecule’s spin configuration54. The key elements in such reactions are radical molecules—transient molecules with an odd number of electrons in their outer molecular shell. Protons, neutrons, and electrons possess spin angular momentum, an inherently quantum characteristic. In a simple picture, quantum spins are like tiny magnets; any other spins or magnetic field in the vicinity could alter their states. In the framework of RPM for avian magnetoreception, it is thought that in CRY RPs can be in the form of anionic semiquinone FAD radical (\(\hbox {FAD}^{{\varvec{\cdot }}-}\)) and terminal tryptophan radical (\(\hbox {TrpH}^{{\varvec{\cdot }}+}\))48,55,56,57. It is also well-known that the superoxide radical, \(\hbox {O}_2^{{\varvec{\cdot }}-}\), can be an alternative partner for the flavin radical58,59,60,61. It has also been proposed that RPs can play important roles in the magnetosensitivity of Drosophila’s CC21,25,62.

Applied magnetic fields can also influence oxidative stress in the presence of CRY63. Moreover, the CC rhythmicity is associated with an endogenous rhythm in the generation of ROS64. Thus it seems pertinent to explore the connection between magnetic field effects and ROS role in the CC. It has recently been proposed that radical pairs (RPs) could play roles in other brain functions. Dufor et al.65 propose that weak MFs activate cellular signaling cascade in neural circuits via acting through CRY, most likely by modulating the state of RPs. The authors concluded that the presence of CRY is critical in axon outgrowth under low-intensity repetitive transcranial magnetic stimulation (rTMS). It has also been suggested that RP may help explain xenon-induced anesthesia66 and the lithium effects on mania39. It, therefore, seems RPs could play critical roles in the functionalities of the brain in general and the CC in particular.

Figure 1
figure 1

(a) Simplified models of the circadian clock feedback loop in Drosophila. CLOCK (CLK) and CYCLE (CYC) promote the tim and per genes. PER and TIM first accumulate in the cytoplasm and then enter into the nucleus to block their gene transcription. Upon light, absorption CRY binds to TIM, and this results in the degradation of TIM67,68. (b) Flavinsemiquinone, \(\hbox {FADH}^{{\varvec{\cdot }}}\), and superoxide \(\hbox {O}_2^{{\varvec{\cdot }}-}\) radical pair in CRY, considered in the RPM model in the present work. The radical pair undergoes interconversion between singlet and triplet states.

The circadian oscillations in Drosophila can be modeled by incorporating the formation of a complex between the PER and TIM proteins and introducing negative feedback loops69, which are the key to the rhythmicity of PER and TIM and their mRNA transcription. The models can be described by a set of a few kinetic equations70. However, modeling Drosophila CC70 can be further simplified into two nonlinear equations67. Furthermore, Player et al.71 show that quantum effects such as magnetic field effects and hyperfine interaction of radical pairs can be introduced to the chemical oscillator by considering the quantum effects on the corresponding reaction rates.

Here, we propose that the RPM could be the underlying mechanism behind the lithium treatment effects and MF effects on Drosophila’s CC. MF via the Zeeman interaction and lithium nucleus via HFIs modulate the recombination dynamics of singlet-triplet interconversion in the naturally occurring RPs in the [\(\hbox {FADH}^{{\varvec{\cdot }}}\) ... \(\hbox {O}_2^{{\varvec{\cdot }}-}\)] complex, shown in Fig. 1b, and hence influence the period of the CC.

In the following, we review the quantitative experimental results for the effects of applied magnetic field25 and lithium40 on the period of Drosophila’s CC. Next, we briefly describe the quantum spin dynamics for the radical pair model where the magnetic field effects and the HFIs are relevant. Moving on, we present our singlet yield calculation for the RP system, inspired by the CRY-based model of birds’ avian magnetoreception72. Later we use a simple model for the mathematical presentation of Drosophila’s CC, following the work on Tyson et al.67. Then we introduce the quantum effect to the period of the CC model, and we show the consistency of our model’s predictions and the experimental findings on the magnetic field and lithium treatment effects. Finally, we discuss new predictions for experiments.

Results

Magnetic field and lithium treatment effects on circadian clock and RPM

Results from prior experiments

Here, we focus on the effects of static MF on Drosophila’s CC observed by Yoshii et al.25. The authors conducted experiments to observe the effects of static magnetic fields with different intensities, [0, 150, 300, 500] \(\mu \hbox {T}\), on changes in the period of Drosophila’s CC under blue light illumination, shown in Table 1. These magnetic fields are, excepting the control of \(0\,\mu \hbox {T}\), approximately 3, 6, and 10 times stronger than natural magnetic fields, respectively. That observation revealed that the period alterations significantly depended on the strength of the magnetic field such that the period change reached a maximum of \(0.522\pm 0.072 \,\hbox {h}\) at \(300 \,\mu \hbox {T}\). In this experiment, the geomagnetic field was shielded, and the arrhythmic flies were excluded from the analysis. We also consider the results of the experiment conducted by Dokucu et al.40 observing the effects of chronic lithium administration on Drosophila’s CC for a range of doses [0, 300] mM. It was shown that lithium treatment lengthens the CC with a maximum prolongation of \(0.7 \pm 0.217\) h at 30 mM of lithium compared to zero lithium intake, see Table 2. In that work, the lethality of lithium up to 30 mM was relatively low until the end of the experiments. Here, we consider that 30 mM is the optimal concentration of lithium where all RPs interact with lithium atoms. We assume that the lithium administered in that work was in its natural abundance, 92.5% and 7.5% of \({}^7{{\mathrm{Li}}}\) and \({}^6{{\mathrm{Li}}}\), respectively. Here we will refer to the natural lithium as Li. In our model here, \(0\,\mu \hbox {T}\) of MF and 0 mM of lithium are our control sets for MF and lithium effects on the CC.

Table 1 Period Changes in the free-running rhythm of Drosophila after application of magnetic fields (MFs) under blue light illumination and lithium administration, taken from the work of Yoshii et al.25.
Table 2 Period in the free-running rhythm of Drosophila for zero and 30 mM intake of lithium, taken from the work of Dokucu et al.40.

RPM model

We develop an RP model to reproduce static MFs and lithium administration effects on the rhythmicity of Drosophila’s CC observed in Ref.25 and Ref.40, respectively. Taking into account the facts that the CC is associated with oxidative stress levels under light exposure45,46,64,73,74 and applied MF75,76, and the CC is affected by lithium intake, we propose that the applied magnetic field interacts with the spins of RPs on FADH and superoxide, and the nuclear spin of lithium modulates the spin state of the radical on superoxide. The correlated spins of RP are assumed to be in the [\(\hbox {FADH}^{{\varvec{\cdot }}}\) ... \(\hbox {O}_2^{{\varvec{\cdot }}-}\)] form, where the unpaired electron on each molecule couples to the nuclear spins in the corresponding molecule, see Fig. 1b. In [\(\hbox {FADH}^{{\varvec{\cdot }}}\) ... \(\hbox {O}_2^{{\varvec{\cdot }}-}\)].

We consider a simplified system in which the unpaired electron is coupled to the flavin’s nitrogen nucleus with an isotropic HF coupling constant (HFCCs) of \(431.3 \,\mu \hbox {T}\)77. In this model, for simplicity, we consider only Zeeman and HF interactions48,78. Following the work of Hore72, the anisotropic components of the hyperfine interactions are excluded, which are only relevant when the radicals are aligned and immobilized79. The RPs are assumed to have the g-values of a free electron. The Hamiltonian for the RP system reads as follows:

$$\begin{aligned} {\hat{H}}=\omega {\hat{S}}_{A_{z}}+a_{A} \hat{{\mathbf {S}}}_A.\hat{{\mathbf {I}}}_{A}+\omega {\hat{S}}_{B_{z}}+a_{B} \hat{{\mathbf {S}}}_B.\hat{{\mathbf {I}}}_{B}, \end{aligned}$$
(1)

where \(\hat{{\mathbf {S}}}_A\) and \(\hat{{\mathbf {S}}}_B\) are the spin operators of radical electron A and B, respectively, \(\hat{{\mathbf {I}}}_A\) is the nuclear spin operator of the isoalloxazine nitrogen of \(\hbox {FADH}^{{\varvec{\cdot }}}\), similar to Refs.39,72, \(\hat{{\mathbf {I}}}_B\) is the nuclear spin operator of the Li nucleus, \(a_{A}\) and \(a_{B}\) are HFCCs, taken from39,77, and \(\omega \) is the Larmor precession frequency of the electrons due to the Zeeman effect. Ref72 uses two nitrogen atoms, while for the purpose of the computational cost we use only the largest one without loss of generality. Of note, oxygen has a zero nuclear spin and thus its HFCC equals zero, (\(a_{B}=0\)), however in the model for lithium effects \(a_{B}\) corresponds to the nuclear spin of lithium. We assumed that the RPs start off from singlet states (see the Discussion section).

Singlet yield calculation

The singlet yield resulting from the radical pair mechanism can be obtained by solving the Liouville-von Neumann equation for the spin state of the radical pair throughout the reaction. Using the eigenvalues and eigenvectors of the Hamiltonian, the ultimate singlet yield, \(\Phi _S\), for periods much greater than the RP lifetime72 has the following form:

$$\begin{aligned} \Phi _S=\frac{1}{4}-\frac{k}{4(k+r)}+\frac{1}{M}\sum _{m=1}^{4M} \sum _{n=1}^{4M}|\langle m|{\hat{P}}^S |n\rangle |^2 \frac{ k(k+r)}{(k+r)^2+(\omega _m-\omega _n)^2}, \end{aligned}$$
(2)

where \(M = M_A M_B\), \(M_{X} =\prod \limits _{i}^{N_X} I_{iX}(I_{iX}+1)\), is the nuclear spin multiplicity, \({\hat{P}}^S\) is the singlet projection operator, \(|m\rangle \) and \(|n\rangle \) are eigenstates of \({\hat{H}}\) with corresponding eigenenergies of \(\omega _m\) and \(\omega _n\), respectively, k is the RP reaction rate, and r is the RP spin-coherence lifetime rate (relaxation rate).

Here we look at the sensitivity of the singlet yield to changes in the strength of the external magnetic field for the [\(\hbox {FADH}^{{\varvec{\cdot }}}\) ... \(\hbox {O}_2^{{\varvec{\cdot }}-}\)] radical complex. Figure 2 illustrates the dependence of the singlet yield of the [\(\hbox {FADH}^{{\varvec{\cdot }}}\) ... \(\hbox {O}_2^{{\varvec{\cdot }}-}\)] complex on external magnetic field B with a maximum yield in [280–360] \(\mu \)T for \(k=4\times 10^{7}\,\hbox {s}^{-1}\) and \(r=3\times 10^{7}\, \hbox {s}^{-1}\) with \(a_{1A}=431.3\,\mu \hbox {T}\). In our model, the magnetic dependence of singlet yield is the foundation of the magnetic sensitivity of the circadian clock. Using the singlet yield, we can reproduce the experimental finding on the effects of applied MF25 and lithium administration40 on the period of the circadian clock of Drosophila, as we discuss below. It is worth mentioning that the singlet-product of the RP system in [\(\hbox {FADH}^{{\varvec{\cdot }}}\) ... \(\hbox {O}_2^{{\varvec{\cdot }}-}\)] is \(\hbox {H}_2\hbox {O}_2\)80, which is the major ROS in redox regulation of biological activities and signaling81.

Figure 2
figure 2

The dependence of the singlet yield of the [\(\hbox {FADH}^{{\varvec{\cdot }}}\) ... \(\hbox {O}_2^{{\varvec{\cdot }}-}\)] complex on external magnetic field B for \(a_{1A}=431.3\,\mu \hbox {T}\), reaction rate k, and relaxation rate r. The singlet yield reaches a minimum value of 48.45% in [280–360] \(\mu \hbox {T}\) (see the inset).

Circadian clock model

We use a simple mathematical model for the circadian clock of Drosophila, following the work of Tyson et al.67. Despite its simplicity, the model is very well known69,82,83,84 and captures the most important part of the clock’s function. In this model, PER monomers are rapidly phosphorylated and degraded, whereas PER/TIM dimers are less susceptible to proteolysis, shown in Fig. 1a. In this context, it is also assumed that the cytoplasmic and nuclear pools of dimeric protein are in rapid equilibrium. With these considerations, it is possible to write the mathematical model in two coupled equations as follows:

$$\begin{aligned} \frac{dM(t)}{dt}= & {} \frac{v_m}{1+(P_t(t)(1-q(t))/2P_{crit})^2}-k_m M(t), \end{aligned}$$
(3)
$$\begin{aligned} \frac{dP_t(t)}{dt}= & {} v_d M(t)-\frac{k_{p1} P_t(t) q(t)+k_{p2} P_t(t)}{J_p+P_t(t)}-k_{p3} P_t(t), \end{aligned}$$
(4)

where \(q(t)=\frac{2}{1+\sqrt{1+8K_{eq}P_t(t)}}\), \(P_t(t)\) and M(t) are the total protein and the mRNA concentrations, respectively. For the descriptions and values of the parameters, see Table 3. In this simple model, \(k_{p3}\) represents the role of CRY’s light activation and hence proteolysis of protein. By solving Eqs. 3 and 4, we obtain the oscillation of protein and mRNA concentrations. Figure 3 shows the explicit time-dependence of protein and mRNA concentrations and the parametric representation of the chemical oscillator limit cycle for Drosophila’s CC. To obtain the period of the clock, we take the average differences between successive peaks and likewise troughs of either \(P_t(t)\) or M(t) by keeping track of when the derivative is zero.

Table 3 Parameter values for the circadian clock of Drosophila, taken from the work of Tyson et al.67. \(\hbox {C}_m\) and \(\hbox {C}_p\) are characteristic concentrations for mRNA and protein, respectively.
Figure 3
figure 3

(a) Explicit time-dependence of the concentrations of protein [red] and mRNA [blue] and (b) Parametric representations of oscillations in the concentrations of protein and mRNA, shown as a limit cycle in Drosophila’s circadian clock model using Eqs. (3) and (4), and the parameters from Table 3.

Effects of singlet yield change on circadian clock

The effects of applied magnetic fields and hyperfine interactions can be introduced to the chemical oscillator of the circadian clock by modifying the rate \(k_f\)71, following the work of Player et al., see Methods. In the CC Eqs. (3) and (4) the corresponding rate is \(k_{p3}\), which represents the role of CRY’s light activation and hence proteolysis of protein and is 0.1 \(\hbox {h}^{-1}\) for the natural cycle of the clock. Hence for the occasions with no singlet yields effects, this value must be retained. The singlet yield effects on \(k_{p3}\) can be written as follows:

$$\begin{aligned} k'_{p3} \propto k_{p3} \frac{\Phi _S'}{\Phi _S}, \end{aligned}$$
(5)

where \(k'_{p3}\), \(\Phi _S\), and \(\Phi '_S\) are the modified rate constant \(k_{p3}\), the singlet yield with no quantum effects, and the singlet yield resulted from quantum effects due to the Zeeman and/or hyperfine interactions, respectively.

Based on the above considerations, here, we calculate the explicit effects of an applied magnetic field and the hyperfine interactions on the period of the CC. Using Eqs. (3), (4), and (5), we explored the parameter space of relaxation rate r and recombination rate k in order to find allowed regions for which our model can reproduce both experimental findings of static MF of \(300 \,\mu \hbox {T}\)25 and 30 mM of lithium40 effects on Drosophila’s CC, which respectively lengthen the clock’s period by \(0.224\pm 0.068 \,\hbox {h}\) and \(0.567\pm 0.11 \,\hbox {h}\). The results are shown in Fig. 4. We find an allowed region where the model reproduces both experiments, see Fig. 4. The parameters for calculating the period of the circadian clock are taken from Table 3. As discussed above, \(k_{p3}\) corresponds to the degradation of TIM due to blue light exposure. For the MF effects under blue light illumination, we set \(k_{p3}=0.085\) \(\hbox {h}^{-1}\) to obtain the control period of the circadian clock \(25.8\pm 0.14 \,\hbox {h}\) under blue light illumination observed in Ref.25. Figure 5 shows the effects of lithium on the rhythmicity of CC, such that \({}^7{{\mathrm{Li}}}\) lengthens the period of the clock longer than \({}^6{{\mathrm{Li}}}\). For the effects of lithium on the circadian clock, the geomagnetic field of \(50\,\mu \hbox {T}\) is taken into account. Figure 6 shows the effects of \(300 \,\mu \hbox {T}\) MF on the CC.

Figure 4
figure 4

The RPM model can reproduce both magnetic field and lithium effects. The comparison between period changes due to applied magnetic fields measured in the experiment25, \(\Delta _{expt-MF}\), and obtained by the RPM model, \(\Delta _{RPM-MF}\), where \(\Delta _{MF}\) is the difference between period changes at 300 \(\mu \)T and 0 \(\mu \)T, \(\Delta _{MF}=period_{300 \mu T}-period_{ 0 \mu T}\). The solid blue line indicates \(\Delta _{expt-MF}-\Delta _{RPM-MF}=0\,\hbox {h}\) and the dashed blue line indicates the region where \(|\Delta _{expt-MF}-\Delta _{RPM-MF}|\le std_{expt-MF}=0.089\,\hbox {h}\). The difference between period changes due to the lithium administration measured in the experiment40, \(\Delta _{expt-Li}\), and obtained by the RPM model, \(\Delta _{RPM-Li}\) is presented by red lines. The solid red line indicates \(\Delta _{expt-Li}-\Delta _{RPM-Li}=0\) h and the dashed red line indicates the region where \(|\Delta _{expt-Li}-\Delta _{RPM-Li}|\le std_{expt-Li}=0.214\,\hbox {h}\), \(\Delta _{Li}=period_{30 mM}-period_{ 0 mM}\). The green shaded color indicates the regions where the RPM model can reproduce both magnetic field25 and lithium40 effects on Drosophila’s CC. The parameters for calculating the period of the circadian clock are taken from Table 3, except that for the MF effects under blue light illumination \(k_{p3}=0.085\,\hbox {h}^{-1}\).

Figure 5
figure 5

Lithium effects on the circadian clock are reproduced by the RPM model. (a) Explicit time-dependence of the concentrations of protein [the solid lines] and mRNA [the dashed lines] and (b) Parametric representations of oscillations in the concentrations of protein and mRNA, in Drosophila’s circadian clock model using the parameters from Table 3. The black, red, blue and purple colors indicate zero-lithium, \({}^{6}{\mathrm{Li}}\), \({}^{7}{{\mathrm{Li}}}\), and Li, respectively. Lithium administration prolongs the period of the clock, such that \({}^{7}{{\mathrm{Li}}}\) has more potency than \({}^{6}{{\mathrm{Li}}}\).

Figure 6
figure 6

Magnetic field effects on the circadian clock are reproduced by the RPM model. (a) Explicit time-dependence of the concentrations of protein [the solid lines] and mRNA [the dashed lines] and (b) Parametric representations of oscillations in the concentrations of protein and mRNA, in Drosophila’s circadian clock model using the parameters from Table 3, except \(k_{p3}=0.085\,\hbox {h}^{-1}\). The black and green colors indicate zero-MF and \(300\,\mu \hbox {T}\) MF effects, respectively.

The model here reproduces the dependence of the CC’s period on the applied MF’s strength and Li administration, shown in Fig. 7. The model predicts that further increases in the intensity of the MF would shorten the period of the clock significantly. For the cases considering MF effects solely, for both the experimental data and the RPM model, the period reaches a maximum between \(0\,\mu \hbox {T}\) and \(500 \,\mu \hbox {T}\) and exhibits reduced effects at both lower and slightly higher field strengths, shown in Fig. 7a. For the cases of \({}^6{{\mathrm{Li}}}\) or without lithium intake, the largest prolongation of the period occurs in the same range of magnetic field as well, shown in Fig. 7b. Another prediction of the model is that \({}^7{{\mathrm{Li}}}\) prolongs the clock’s period stronger than \({}^6{{\mathrm{Li}}}\), which has a smaller spin compared to \({}^7{{\mathrm{Li}}}\), see Fig. 7b. In this model, in the cases where lithium effects are considered, the geomagnetic effects of \(50 \,\mu \hbox {T}\) are also considered. For the comparison between our model and the experimental data on the lithium effects, we assume that natural lithium was administered in the experiment40. Following Ref.39, we considered \(a_{{{}^7{{\mathrm{Li}}}}}=-224.4\,\mu \hbox {T}\). Note that moderate changes to the HFCC of lithium do not impact our results. Figure 7 shows that the dependence of the period on applied MFs and lithium effects calculated by the RPM model used in the present work is consistent with the experimental observations. We compare the maximum lengthening of the period in both the RPM model and experimental data25, \(\Delta _{RPM}=0.154\,\hbox {h}\) and \(\Delta _{expt}=0.22\pm 0.089\,\hbox {h}\), respectively, where \(\Delta =period_{max}-period_{min}\). The results from the RPM fall into the uncertainty of the experimental data, \(|\Delta _{expt}-\Delta _{RPM}|\le std_{expt}\).

Figure 7
figure 7

The dependence of the period of Drosophila’s circadian clock calculated by the RPM model on the static magnetic field strength B without (a) and with (b) lithium effects for \(a_A = 431.3\) \(\mu \)T, \(a_B = a_{{{}^7{{\mathrm{Li}}}}} =-224.4\,\mu \hbox {T}\), relaxation rate \(r = 3\times 10^7\, \hbox {s}^{-1}\), and reaction rate \(k= 4\times 10^7\,\hbox {s}^{-1}\). Higher magnetic field intensities shorten the period of the circadian clock. For the case without lithium effects (a), the applied magnetic field lengthens the period of the clock to a maximum in [280–360] \(\mu \)T and reduces effects at both lower and higher field strengths. The comparison between the dependence of the period on applied magnetic field calculated by the RPM model [black line in the inset of plot (a)] and the experimental findings [green dots with error-bars] of Ref.25. (b) The effects of Li [purple], \({}^6{{\mathrm{Li}}}\) [red], \({}^7{{\mathrm{Li}}}\) [blue], and zero Li [black]. The inset indicates the comparison between the effects of Li on the period of the clock calculated by the RPM model [purple line] and the experimental findings [orange dots with error-bars] of Ref.40. The results from the RPM fits into the uncertainty of the experimental data, such that \(|\Delta _{expt}-\Delta _{RPM}|\le std_{expt}\).

Discussion

In this project, we aimed to probe whether a RP model can explain the experimental findings for both the effects of static magnetic field25 and lithium40 on the circadian clock in Drosophila. We showed how the quantum effects affect the rates, which then yields a change in the period of the clock. This is a significant step forward compared to the previous studies on xenon anesthesia66 and the lithium effects on hyperactivity39, where the quantum effects were correlated to experimental findings without explicitly modeling the related chemical reaction networks. With a set of reasonable parameters, our model reproduces the experimental findings, as shown in Figs. 4, and 7. In addition, this strengthens the previously proposed explanation for the effects of lithium on hyperactivity39 via the circadian clock.

We proposed that applied magnetic fields and nuclear spins of lithium influence the spin state of naturally occurring [\(\hbox {FADH}^{{\varvec{\cdot }}}\) ... \(\hbox {O}_2^{{\varvec{\cdot }}-}\)] radical pairs in the circadian clock. This is inspired by the observations that the Drosophila circadian clock is altered by external magnetic fields24,25, which is accompanied by modulations in the ROS level64,75, and by lithium administration40. Let us note that it has also been suggested that lithium exerts its effects by inhibiting Glycogen synthase kinase-3 (GSK-3)85,86 . However, while the presence of RPs is a natural explanation for magnetic field effects, their existence in GSK-3 requires experimental support.

The suggested [\(\hbox {FADH}^{{\varvec{\cdot }}}\)...\(\hbox {O}_2^{{\varvec{\cdot }}-}\)] radical pairs depend on a form of FAD in the Drosophila CRY photocycle that is thought to be unusual. However, Baik et al. showed the possibility that light-activated CRY in Drosophila neurons express a \(\hbox {FADH}^{{\varvec{\cdot }}}\) neutral semiquinone state87. They conclude that further investigation would shed more light on whether the red light responses are due to Drosophila CRY signaling or opsin signaling or a combination of both. On the other hand, Chandrasekaran et al. reported that Drosophila CRY mutants that can form \(\hbox {FADH}^{{\varvec{\cdot }}}\) are incapable of light-activated conformational changes88. However, more recent studies show that alternative signaling pathways may be involved in the magnetosensitivity of the Drosophila circadian clock89. Thus, other forms of radical pairs, e.g. [\(\hbox {FAD}^{{\varvec{\cdot }}-}\) ... \(\hbox {TrpH}^{+.}\)], could also be considered instead of [\(\hbox {FADH}^{{\varvec{\cdot }}}\)...\(\hbox {O}_2^{{\varvec{\cdot }}-}\)] to explain the magnetosensitivity of the Drosophila circadian clock.

Of note, there is a large body of evidence that ROS are involved in the context of magnetosensing and the circadian clock modulations63,64,73,74,76,90. it has been shown that oscillating magnetic fields at Zeeman resonance can influence the biological production of ROS in vivo, indicating coherent S-T mixing in the ROS formation80. Additionally, it has been observed that extremely low frequency pulsed electromagnetic fields cause defense mechanisms in human osteoblasts via induction of \(\hbox {O}_2^{{\varvec{\cdot }}-}\) and \(\hbox {H}_2\hbox {O}_2\)75. Sherrard et al.63 observed that weak pulsed electromagnetic fields (EMFs) stimulate the rapid accumulation of ROS, where the presence of CRY was required63. The authors of that work concluded that modulation of intracellular ROS via CR represents a general response to weak EMFs. Further, Sheppard et al.90 demonstrated that MFs of a few millitesla can indeed influence transfer reactions in Drosophila CRY. It has also been shown that illumination of Drosophila CRY results in the enzymatic conversion of molecular oxygen to transient formation of superoxide \(\hbox {O}_2^{{\varvec{\cdot }}-}\) and accumulation of hydrogen peroxide \(\hbox {H}_2\hbox {O}_2\) in the nucleus of insect cell cultures73. These findings indicate the light-driven electron transfer to the flavin in CRY signaling74.

The feasibility for the \(\hbox {O}_2^{{\varvec{\cdot }}-}\) radical to be involved in the RPM is a matter of debate in this scenario due to its likely fast spin relaxation rate r. Because of fast molecular rotation, the spin relaxation lifetime of \(\hbox {O}_2^{{\varvec{\cdot }}-}\) is thought to be on the orders of 1 ns91,92. Nonetheless, it has also been pointed out that this fast spin relaxation can be decreased on account of its biological environment. Additionally, Kattnig et al.93,94 proposed that scavenger species around \(\hbox {O}_2^{{\varvec{\cdot }}-}\) can also reduce its fast spin relaxation. Moreover, in such a model, the effects of exchange and dipolar interactions can also be minimized.

It is often assumed that in the RP complexes involving superoxide are formed in triplet states, as opposed to the case considered here. This is because the ground state of the oxygen molecule is a triplet state. The initial state for RP formation could also be its excited singlet state, which is a biologically relevant ROS95,96,97. Further, the transition of the initial RP state from triplet to singlet could also take place due to spin-orbit coupling98,99.

Our model predicts that increasing the intensity of the applied magnetic field will shorten the period of the clock. This is a significant new prediction of our model that would be very interesting to check. The isotopic-dependence of the period is another prediction of our present model, such that \({}^7{{\mathrm{Li}}}\) lengthens the period of the clock longer than \({}^6{{\mathrm{Li}}}\).

The circadian clock not only controls the rhythms of the biological processes, but it also has intimate connections to other vital processes in the body100 and particularly in the brain101. It has been suggested that environmental perturbations in the circadian period could increase the risk of selected cancers and hence the circadian clock could be a therapeutic target for cancer risks102. It also appears that the way drugs function depends on the circadian clock103,104. Notably, it has been shown that the circadian clock is vital for maintaining the anti-oxidative defense105. Moreover, it has been suggested that the circadian clock could be a new potential target for anti-aging106,107 and neurodegenerative disorders therapeutics108. Thus this project also paves a potential path to study other functionalities of the body and the brain connected to the circadian clock in the light of the RPM.

To sum up, our results suggest that quantum effects may underlie the magnetic field and lithium effects on the circadian clock. A similar mechanism is likely to be at the heart of magnetoreception in animals109, xenon-induced anesthesia66, and lithium treatment for mania39. Our work is thus another piece of evidence that quantum effects may play essential roles in the brain’s functionalities110,111,112,113,114,115,116,117,118,119,120,121.

Methods

Quantum effects and chemical oscillator

The effects of applied magnetic fields and hyperfine interactions can be introduced, following the work of Player et al.71, by assuming that the FAD signaling of CRY to TIM and hence the TIM degradation in the CC process proceed by RPs, shown in Eq. 6:

(6)

where rate constants \(k_1\) and \(k_{-1}\) are assumed to conserve electron spin. External magnetic fields and the hyperfine interactions can influence the overall rate of production of \(\hbox {FADH}^{{\varvec{\cdot }}}\) by altering the extent and timing of coherent singlet/triplet interconversion in RP and so changing the probability that it reacts to give \(\hbox {FADH}^{{\varvec{\cdot }}}\) rather than returning to \(\hbox {FAD}^*\). Based on a spin dynamics calculation, one can describe the effect of applied magnetic fields and the HFIs on the kinetics of the CC simply by modifying the rate constant \(k_{p3}\)122,123,124,125 which corresponds to the degradation of TIM in Eqs. (3) and (4). The spin dynamics of the RP in Eq. (6) can be written as follows:

$$\begin{aligned} \frac{d{\hat{\rho }}(t)}{dt}=-{\hat{L}} {\hat{\rho }}(t)+k_1 \frac{{\hat{P}}^S}{M}=-i[{\hat{H}},{\hat{\rho }}(t)]-\frac{k_{-1}}{2} \{{\hat{P}}^S,{\hat{\rho }}(t)\}-k_f {\hat{\rho }}(t)+k_1 \frac{{\hat{P}}^S}{M}, \end{aligned}$$
(7)

where \({\hat{L}}\) is the Liouvillian, [..., ...] and \(\{...,...\}\) are the commutator and anti-commutator operators, \({\hat{\rho }}(t)\) is the spin density operator of the RP system; its trace, \(Tr[{\hat{\rho }}(t)]\), equals the concentration of RPs divided by the fixed concentration of \(\hbox {FAD}^*\) in Eq. 6. As RPs are short-lived intermediates, their concentrations are very low, and hence one can obtain the steady-state solutions as follows (see Ref.71):

$$\begin{aligned} \Phi =\frac{k_f}{M}Tr[{\hat{L}}^{-1}{\hat{P}}^S], \end{aligned}$$
(8)

where \(\Phi \) is the singlet yield of the RPM. It is, therefore, possible to introduce the singlet yield of the RPM to the chemical reaction by modifying the rate \(k_f\). In the CC Eqs. (3) and (4) the corresponding rate is \(k_{p3}\), which is \(0.1\, \hbox {h}^{-1}\) for the natural cycle of the clock without blue light illumination and \(0.085\,\hbox {h}^{-1}\) for blue light illumination. The singlet yield effects on \(k_{p3}\) can be written as follows:

$$\begin{aligned} k'_{p3} \propto k_{p3} \frac{\Phi '_S}{\Phi _S}, \end{aligned}$$
(9)

where \(k'_{p3}\), \(\Phi _S\), and \(\Phi '_S\) are the modified rate constant \(k_{p3}\), the singlet yield with no quantum effects, and the singlet yield resulted from quantum effects due to the Zeeman and/or hyperfine interactions, respectively.