Document Detail

Modeling and validating chronic pharmacological manipulation of circadian rhythms.
Jump to Full Text
MedLine Citation:
PMID:  23863866     Owner:  NLM     Status:  PubMed-not-MEDLINE    
Abstract/OtherAbstract:
Circadian rhythms can be entrained by a light-dark (LD) cycle and can also be reset pharmacologically, for example, by the CK1δ/ε inhibitor PF-670462. Here, we determine how these two independent signals affect circadian timekeeping from the molecular to the behavioral level. By developing a systems pharmacology model, we predict and experimentally validate that chronic CK1δ/ε inhibition during the earlier hours of a LD cycle can produce a constant stable delay of rhythm. However, chronic dosing later during the day, or in the presence of longer light intervals, is not predicted to yield an entrained rhythm. We also propose a simple method based on phase response curves (PRCs) that predicts the effects of a LD cycle and chronic dosing of a circadian drug. This work indicates that dosing timing and environmental signals must be carefully considered for accurate pharmacological manipulation of circadian phase.CPT: Pharmacometrics & Systems Pharmacology (2013) 2, e57; doi:10.1038/psp.2013.34; published online 17 July 2013.
Authors:
J K Kim; D B Forger; M Marconi; D Wood; A Doran; T Wager; C Chang; K M Walton
Related Documents :
4052496 - Dynamic and static hysteresis in crayfish stretch receptors.
4974746 - The structure and function of a slowly adapting touch corpuscle in hairy skin.
24167276 - A neuropeptide speeds circadian entrainment by reducing intercellular synchrony.
8903136 - Stochastic retinal mechanisms of light adaptation and gain control.
7863616 - Spatial scale shifts in amblyopia.
15610256 - Equivalence of information from single versus multiple frequency bioimpedance vector an...
Publication Detail:
Type:  Journal Article     Date:  2013-07-17
Journal Detail:
Title:  CPT: pharmacometrics & systems pharmacology     Volume:  2     ISSN:  2163-8306     ISO Abbreviation:  CPT Pharmacometrics Syst Pharmacol     Publication Date:  2013  
Date Detail:
Created Date:  2013-07-18     Completed Date:  2013-07-18     Revised Date:  2014-11-03    
Medline Journal Info:
Nlm Unique ID:  101580011     Medline TA:  CPT Pharmacometrics Syst Pharmacol     Country:  United States    
Other Details:
Languages:  eng     Pagination:  e57     Citation Subset:  -    
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): CPT Pharmacometrics Syst Pharmacol
Journal ID (iso-abbrev): CPT Pharmacometrics Syst Pharmacol
ISSN: 2163-8306
Publisher: Nature Publishing Group
Article Information
Download PDF
Copyright © 2013 American Society for Clinical Pharmacology and Therapeutics
open-access:
Received Day: 13 Month: 02 Year: 2013
Accepted Day: 07 Month: 05 Year: 2013
Print publication date: Month: 07 Year: 2013
Electronic publication date: Day: 17 Month: 07 Year: 2013
pmc-release publication date: Day: 1 Month: 7 Year: 2013
Volume: 2 Issue: 7
First Page: e57 Last Page:
PubMed Id: 23863866
ID: 3734602
Publisher Item Identifier: psp201334
DOI: 10.1038/psp.2013.34

Modeling and Validating Chronic Pharmacological Manipulation of Circadian Rhythms
J K Kim1
D B Forger12*
M Marconi3
D Wood3
A Doran4
T Wager5
C Chang6*
K M Walton37
1Department of Mathematics, University of Michigan, Ann Arbor, Michigan, USA
2Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA
3Neuroscience Research Unit, Worldwide Research & Development, Pfizer Inc, Groton, Connecticut, USA
4Pharmacokinetics, Dynamics and Metabolism, Worldwide Research & Development, Pfizer Inc, Groton, Connecticut, USA
5Neuroscience Medicinal Chemistry, Worldwide Research & Development, Pfizer Inc, Groton, Connecticut, USA
6Systems Modeling and Simulation, Pharmacokinetics, Dynamics and Metabolism, Worldwide Research & Development, Pfizer Inc, Groton, Connecticut, USA
7Current Address: Division of Pharmacotherapies and Medical Consequences of Drug Abuse, National Institute on Drug Abuse, NIH, Rockville, Maryland, USA
*(forger@umich.edu)
*(cheng.chang@pfizer.com)

The timing of activities such as waking, sleeping, body temperature, blood pressure, hormone expression, and feeding show circadian (daily) rhythms.1,2 These circadian rhythms are regulated by the master circadian clock in the suprachiasmatic nuclei where transcriptional activators CLOCK and BMAL1 drive the expression of repressors period (Per) and cryptochrome (Cry).3 This feedback system includes the PER1/2 proteins, which are phosphorylated by the CK1δ/ε, dimerize with the CRYs, then translocate to the nucleus to inhibit BMAL1/CLOCK and repress the transcription of Per and Cry.4 Further phosphorylation by CK1δ/ε signals PER degradation that releases BMAL1/CLOCK transcriptional inhibition and resumes transcription of Per and Cry.2

This endogenous timekeeping system can be synchronized to the earth's 24-h periodic environment through external cues, known as zeitgebers (e.g., light-dark (LD) cycle and temperature cycle).5,6 To maintain clock-environment synchrony, zeitgebers induce changes in the concentrations of the molecular components of the clock to levels consistent with the appropriate stage in the 24-h cycle. Misalignments of circadian timing with the external environment can cause significant physiological problems, such as jet lag, depression, insomnia, coronary heart disease, neurodegenerative disorders, and cancer.7 In particular, mood disorders and bipolar disorders appear to be tightly related to disrupted circadian rhythms.8,9,10,11 To treat the misalignment of circadian clocks with the external environment, pharmacological manipulation of circadian clocks has received much attention.12,13,14,15 We previously showed that acute dosing of PF-670462 (CK1δ/ε inhibitor) can delay circadian behavior, as well as re-establish a circadian rhythm in Vipr2-/- mice that are arrhythmic under dark-dark (DD) cycle or light-light (LL) cycle.16,17,18 To extend this work to real-life situations that proceed under LD cycles with seasonal variation, we need to study the effect of CK1δ/ε inhibition on circadian rhythms under different LD cycles.

Light and inhibition of CK1δ/ε simultaneously affect multiple components in the molecular feedback loops in circadian clocks.4,17,18,19,20,21 Light induces transcription of Per1 and Per2, whereas inhibition of CK1δ/ε decreases the degradation rate and the nuclear translocation rate of PER as well as the binding rate between PER and CRY.2 To understand these interactions systematically, mathematical modeling has played important roles. For instance, the correct function of the tau mutation in CK1ε was identified by the Forger–Peskin model, which was later confirmed experimentally.22,23

Here, we study how light stimuli and CK1δ/ε inhibition affect mammalian circadian timekeeping with a combination of experiments and simulations using a mathematical model of intracellular mammalian circadian clocks.24 We find that acute pharmacologic inhibition of CK1δ/ε via PF-670462 immediately delays all clock gene expression in suprachiasmatic nuclei and locomotor activity under LD cycles. The opposing actions of pharmacological delay and light can yield a constant stable delay of circadian behavior when CK1δ/ε is inhibited chronically under LD cycles. The occurrence and magnitude of a stable phase delay depend on dosing amount, dosing timing, and day lengths. We also find that complex behaviors induced by the LD chronic dosing can be predicted with phase response curves (PRCs) to light stimuli and CK1δ/ε inhibition. This work provides a way to determine a dosing strategy of chronic CK1δ/ε inhibition to treat the misalignment of circadian clocks by modulating the phase of circadian rhythms.


Results
The counteracting effects of light on multi-day dosing of PF-670462

We previously showed that the CK1δ/ε inhibitor, PF-670462 delays circadian phase under a DD cycle regardless of dosing timing.16,17,18 Here, we examine the effects of multi-day dosing of PF-670462 at CT11 (ZT11), the dosing timing that yielded the largest phase delay16 in both the DD and LD cycles (Figure 1a). CT and ZT refer to timing with respect to a circadian signal or external signal, respectively (e.g., for a nocturnal animal, time of the onset of activity and lights off define CT12 and ZT12, respectively).5 In a LD cycle, 10 and 32 mg/kg dosing induce phase delays of 1.3 and 2.9 h, respectively (Figure 1b and Supplementary Table S1 online). Treatment in DD, however, yields phase delays of 1.5 and 3.9 h at 10 and 32 mg/kg, respectively, which are greater in comparison to LD dosing (Figure 1b and Supplementary Table S1 online). These multi-day–dosing experiments indicated that a LD cycle partially counterbalances the effect of PF-670462 on circadian rhythms.

Chronic dosing of PF-670462 induces a constant stable delay under LD cycles

The 3-day dosing experiments indicates that the phase delays from multiple days of dosing are additive. Given the significant effect of light to partially counteract the delay in activity, we sought to determine whether daily dosing in LD would be continuously additive, as has been previously shown for extended DD dosing.18 Daily dosing of PF-670462 was performed in both the LD and DD cycles for up to 3 weeks (Figure 2a). The LD experiment was terminated after 2 weeks upon the observation of a constant daily phase delay, whereas the dosing in DD was continued for a third week. In DD, large accumulated delays of 7.6 and 13 h were detected at 10 and 32 mg/kg, respectively (Figure 2b,d and Supplementary Table S2 online). In a LD cycle, the reduced accumulated delays of 1.8 and 6 h were observed with 10 and 32 mg/kg, respectively (Figure 2c,d and Supplementary Table S2 online). This stable delay of LD chronic dosing, which was also observed in a previous study with rats,25 indicates that equilibrium has been reached by balancing light-induced phase advances and pharmacological-induced phase delays.

Effect of acute CK1δ/ε inhibition on gene expression and animal behavior under LD cycles

To extend our experimental results to different LD cycles or different dosing timings, we used an accurate and detailed mathematical model of the intracellular mammalian circadian clock,24 which accurately predicts many phenotypes of circadian mutations including the tau mutation in CK1ε. To use the molecular model, we first tested whether the effects of CK1δ/ε inhibition in animals entrained to a LD cycle were reflected both in gene expression and in whole animal locomotion since intracellular timekeeping is often not fully reflected in whole animal behavior.26

Mice housed in a 12:12 LD cycle were treated with PF-670462 at ZT11. Starting at 4 h following dosing, the hypothalamus, a region of the brain which includes the master circadian clock of mice, was harvested at 3-h intervals and the expression of eight clock genes was measured. Animals continued in the LD cycle until hypothalamus harvesting to test the effect of both the PF-670462 and LD cycle together. All the genes examined showed a dose-responsive delay in expression under LD cycles (Supplementary Figure S1 online and Supplementary Table S3 online). In particular, Rev-erbα and DBP showed the clearest shifts in the peak times of expression. While DBP showed no shift of peak time for 32 mg/kg dosing, other timepoints were shifted, indicating the overall phase shift might be <3 h, the time interval used to measure gene expression (Figure 3a).

We next examined how the delayed circadian gene expression is reflected in circadian behavior. Under treatment conditions identical to the gene expression studies, activity was measured by telemetry. Following treatment and one additional 12:12 LD cycle, reflecting the day of hypothalamus harvesting in the gene expression studies, the mice were transferred into DD to assess phase shifts (Figure 3b). Circadian behavior also showed a dose-response delay under a LD cycle (Supplementary Table S4 online). Delays of 1.3 and 6.7 h, for dosing of 32 and 100 mg/kg were found, respectively. Although the shift induced by 32 mg/kg is not significant (Supplementary Table S4 online), the behavioral shifts are remarkably consistent with the gene expression shifts of Rev-erbα and DBP (Figure 3c).

Taken together, these behavioral and gene expression results are consistent with a previous study showing that Rev-erbα and DBP expression best reflect altered locomotor activity in a natural LD cycle.27 The data also indicate that circadian behavior is strongly linked to the timing of clock genes when gene expression is altered by CK1δ/ε inhibition in an acute manner even under LD cycles.

A systems pharmacology model of the mammalian circadian clock including PF-670462

Since our experiments showed a tight relationship between molecular and behavioral rhythms in response to CK1δ/ε inhibition under a LD cycle (Figure 3), we hypothesized that our molecular model would give accurate predictions of both molecular timekeeping and whole animal behavior.24 Since this model does not include PF-670462, a multi-compartment pharmacokinetic/pharmacodynamic model was included into our model to describe disposition of PF-670462 and its interaction with CK1δ/ε28 (Supplementary Figure S2 online and Supplementary Methods online). This adds 11 parameters that describe the pharmacokinetics and pharmacodynamics of PF-670462 in the model (Supplementary Table S5 online). The values of these new parameters were estimated by fitting to experimental data via a simulated annealing (SA) method.29 See Methods and Supplementary Table S5 online for details of description and estimation of new parameters.

The estimated parameters reflect the pharmacokinetics and pharmacodynamics of PF-670462 accurately. Our parameters are fitted to the disposition profiles of PF-670462 in plasma and brain tissue following a single 32 mg/kg s.c. (Figure 4a). Furthermore, interestingly, the model correctly predicts CK1δ/ε occupancy, the fraction of bound CK1δ/ε by PF-670462 (Figure 4b), suggesting accurate prediction of the binding affinity between PF-670462 and CK1δ/ε.

Next, we compared simulations of our model with previous in vitro experimental studies that measured the effect of PF-670462 on clock gene expression.18 Our model successfully predicts dose-dependent period changes of clock gene expression in suprachiasmatic nuclei (Figure 4c). Matching previously published data,18 the model also predicts that period prolongation by CK1δ/ε inhibition is mainly due to the prolongation of the interval between the peaks of Per2 and Bmal1 (Figure 4d).

Finally, we compared the phase shifts of behavior produced by light and CK1δ/ε inhibition in the model and in the experiments. For this, we explored PRCs, which are measured by giving a stimulus (e.g., light) to circadian rhythms at different times and measuring the effect on the phase of rhythms.30,31 Previous studies have shown that PRCs of circadian rhythms to light pulses have both the advance and delay regions,32,33,34 but the PRC to dosing of PF-670462 shows only delays.16 That is, light stimuli can advance or delay circadian rhythms depending on the timing of stimuli, but a dose of PF-670462 always delays the phase, regardless of dosing timing. Our mathematical model successfully reproduced these two PRCs (Figure 4e) and was also able to reproduce both a constant stable delay induced by LD chronic dosing and a cumulative increasing delay induced by DD chronic dosing at ZT11 (Figures 2 and 4f). To test the reliability of estimated parameters, we also simulated PRCs to dosing of PF-670462 in the presence of parameter perturbations. Even with the significant perturbations of parameters, the model successfully produced PRCs that show only delays (Supplementary Figure S3 online). These simulations indicate that the model can accurately and robustly reproduce behavioral data in response to CK1δ/ε inhibition in a LD cycle.

Predictions about chronic dosing of a PF-670462

Our mathematical model was successfully able to reproduce experimental data on the effects of a light pulse or CK1δ/ε inhibition on circadian rhythms (Figure 4). With our mathematical model, we simulated chronic dosing of PF-670462 under various conditions. First, we investigated whether LD chronic dosing with 32 mg/kg of PF-670462 can induce a constant stable phase delay regardless of dosing timing. Based on our experiments showing that LD chronic dosing at ZT11 caused a constant stable phase delay (Figure 2c),16 one might assume that LD chronic dosing at other times might also lead to constant stable phase delays. However, the model surprisingly predicted that this was not the case. In fact, our model predicts that stable entrainment does not occur, i.e., no stable relationship between the LD cycle and circadian phase is achieved, when the inhibitor is applied during the early night (Figure 5a). This contrasts with late night dosing which is predicted to produce minimal phase shifts. The model also predicts that stable entrainment can occur during daytime dosing, however, the magnitude of the phase delay is predicted to vary greatly depending on the dosing timing (Figure 5a). Thus, we predict that dosing timing must be very carefully controlled to achieve a desired phase delay.

Since the timing of exact daily dosing is predicted to have a significant effect on the ability of the circadian clock to entrain to a LD cycle, we next wondered whether the accuracy of the dosing timing was also an important factor. To explore this, we simulated 32 mg/kg chronic dosing that varies somewhat but centers around either ZT2, ZT5, ZT8 or ZT11 under a 12:12 LD cycle (Figure 5b). While entrainment could still be seen in the presence of variations in dosing timing, the phase of the circadian clock was, as expected, less controlled particularly near ZT11 (Figure 5b). In fact, the variability in dosing timing often prevented entrainment from occurring at ZT11 (Figure 5b). These results indicate that LD chronic dosing in the morning or early afternoon is more likely to induce a stable and robust phase delay in the presence of a less-controlled dosing schedule.

Finally, the effect of the seasonal change of day length on chronic dosing was explored. Since light opposes the phase delays induced by chronic dosing and produces the stable entrainment (Figures 2c and 5a), we initially expected that increasing the amount of available light (e.g., to a 16:8 LD cycle) would allow for entrainment over a wider range of dosing times. However, our simulations showed the opposite results: as the light duration lengthens, entrainment is less likely to occur, and was lost for ZT11 dosing (Figure 5c). The ranges of dosing timings that induce the stable entrainment become narrower as day length increases (Figure 5d). These results indicate that dosing schedules should be adjusted according to the short and long day lengths, corresponding to winter and summer, to ensure stable entrainment.

Understanding complex behavior of LD chronic dosing via PRC analysis

We have shown that the interaction of two zeitgebers, light and dose of PF-670462 can lead to complex behavioral patterns, which depend on dosing amount, dosing timing, and day length (Figure 5). With a fixed relationship between the lighting schedule and the dosing timing, we can consider the light and dosing as part of one combined signal to shift the circadian clock. The effect of this combined signal can be understood through PRC analysis. As an example, we can consider dosing at ZT5 under a 12:12 LD cycle as a combined signal by a 12-h 300 lux light stimulus and subsequent CK1δ/ε inhibition occurring 5 h after light onset. To measure phase shifts induced by the combined signal at different phases, we simulated a PRC to the combined signal (Figure 6a). If the 12-h light pulse begins at circadian phase 0 (CT0), with subsequent dosing 5 h later, the PRC predicts a 1-h delay in phase due to the combined signal (Figure 6a). Thus, on the second day, the light onset begins at CT23. The PRC at CT23 predicts a phase delay of about 0.5 h after the second day of dosing. After several days, light onset begins near CT22, where the PRC is zero, no more phase shifting occurs, and entrainment is reached. This matches our ZT5 chronic-dosing simulation (Figure 5a).

Similar arguments can be used to predict the combined signal at different phase relationships between the dosing time (ZT2, ZT5, ZT8, ZT11, and ZT14) and LD cycle (Figure 6b). All of these PRCs correctly predicted whether the constant stable delays of the circadian rhythms occur and if a constant stable delay occurs, to what extent the phase is delayed. That is, the zero crossing of the PRCs correctly predicts the extent of phase delay with LD chronic dosing (Figures 5a and 6b). Moreover, the fact that the PRC with ZT14 dosing never takes the value of 0 explains why stable entrainment would not occur with LD chronic dosing at ZT14 (Figures 5a and 6b).

Next, to understand the effect of day length on chronic dosing, we simulated PRCs to a 300 lux light pulse with different durations and the subsequent 32 mg/kg PF-670462 dosing 11 h after light onset (ZT11). Interestingly, as light duration lengthens, the PRCs move downward and the PRCs to long light duration did not have zero crossings (Figure 6c). This explains why the ZT11 chronic dosing was not able to yield a stable entrainment during long days (Figure 5c).

The analysis of combined PRCs can also be used to predict the effect of free-running period from 24 h. If free-running period is longer than 24 h, the phase would be delayed every day, and this difference in period from 24 h must be made up by phase shifts from light and the inhibitor to achieve entrainment. We can incorporate this into our combined PRC analysis by subtracting the difference in free-running period from 24 h at all points. Due to the subtraction, the PRCs have less of a chance to have zeros, so stable entrainment is less likely to occur as the free-running period becomes longer, which matches simulations of the model (Supplementary Figure S4 online).

Estimation of the combined PRCs from two separate PRCs

We have shown that PRCs to a combined stimulus of light and PF-670462 dosing can be used to understand and predict phase changes caused by LD chronic dosing (Figure 6). However, experimentally measuring these combined PRCs for all combinations of different dosing timing and light duration require a tremendous amount of experimental work. Thus, we were curious whether the combined PRCs could be estimated from the already measured individual PRCs to light and PF-670462, which were successfully reproduced by our model (Figure 4e).16,34

To estimate the combined PRCs to a 12-h light pulse and ZT5 dosing, we added the PRC to 12-h light and a 5-h translated PRC to 32 mg/kg PF-670462 (Figure 7a). For example, the phase response to light at ZT0 and the phase response to the PF-670462 at CT5 are added. The summation of these two PRCs well matches the combined PRC with ZT5 dosing (Figure 7a). In a similar way, we added the light PRC and the PF-670462 PRC with time translations of 2, 8, 11, 14, 17, 20, and 23 h to estimate the combined PRCs to ZT2, ZT8, ZT11, ZT14, ZT17, ZT20, and ZT23 dosing, respectively (Figure 7b). The summed PRCs also nicely match with the combined PRCs and, in particular, correctly predict the zero crossings of combined PRCs. Thus, we can understand and predict LD chronic dosing-induced phase change by measuring the separate light PRC and PF-670462 PRC. For instance, as light duration lengthens, the advance region of light PRC of mice becomes smaller,34 which indicates that the summation of light PRCs and PF-670462 PRCs is more likely to be all delay or zeros of PRC is less likely to exist. This explains why stable entrainment was less likely to occur during chronic dosing as light duration lengthens (Figure 5c).


Discussion

We have employed gene expression quantitation, behavioral measures, and mathematical modeling together to study how the phases of circadian rhythms are modulated by two zeitgebers, the LD cycle, and pharmacological CK1δ/ε inhibition. By using acute inhibition of CK1δ/ε under a LD cycle, we have shown a direct and immediate relationship between the delay of circadian gene expression and circadian behavior (Figure 3). Our mathematical intracellular model also successfully reproduced the altered behaviors associated with CK1δ/ε inhibition in LD cycles (Figure 4).

The inhibition of CK1δ/ε over several days increased the phase shift in behavior. This behavioral shift is attenuated when the animal is maintained in LD (Figure 1). The counterbalancing of the delay by a LD environment is particularly dramatic when it occurs over multiple weeks. While chronic dosing of PF-670462 at ZT11 under DD cycle yielded a continuously accumulating delay, those under LD cycle induced a constant stable behavior shift in a dose-dependent manner, consistent with previous studies (Figure 2).18,25 The effect suggests that LD chronic dosing can be used to treat circadian rhythm sleep disorders, in particular advanced sleep phase disorder which is characterized by a several hour earlier sleep schedule than what is desired. Currently, treatment of advanced sleep phase disorder involves timed exposure to bright light in the evening (chronotherapy) or dosing of melatonin in the morning, each of which have their own limitations due to the lengthy process and sedating effect, respectively.13,35,36,37 These issues may be addressed by the LD chronic dosing of PF-670462 (Figure 2). Although there are side-effect concerns since CK1δ/ε are key regulators of diverse cellular growth and survival processes,38 a recent study showed that PF-670462 has only a modest effect on cell proliferation.39

To reliably use the LD chronic dosing of PF-670462 to modulate the phase of circadian rhythms, we explored the effect of the chronic dosing under various conditions. With an accurate and detailed mathematical model of mammalian circadian clock,24 we studied how dosing amount, dosing timing, and day length affect the chronic dosing-induced phase delay (Figure 5). Since chronic dosing at ZT11 under a 12:12 LD cycle showed a constant stable delay (Figure 2), dosing at other times was also expected to yield a constant stable delay since acute dosing at ZT11 induces the largest delay (Figure 4e).16 Moreover, considering the counterbalancing effects of PF-670462 and light on phase, a constant stable delay is expected as the day length becomes longer (Figure 1). However, simulations predict that the stable entrainment is less likely to occur when PF-670462 is dosed at early night and when day length becomes longer (Figure 5a,c). Furthermore, simulations predict that stable entrainment is less likely to occur as free-running period becomes longer (Supplementary Figure S4 online).

We found that these complex patterns can be understood by analyzing PRCs to a combination of light pulse and dosing of PF-670462. That is, the occurrence of a constant stable delay during LD chronic dosing depends on the existence of zeros of the combined PRC and to what extent phase delay can be estimated from the position of zeros (Figure 6). Moreover, we showed that these combined PRCs are successfully approximated by simple summation of a light PRC and a translated PF-670462 PRC (Figure 7). This result indicates that if we can measure a light PRC and a PF-670462 PRC for other organisms including humans, we can determine the optimized dosing timing and amount of LD chronic dosing to modulate the phase of circadian rhythms. Since the PRCs to light or PF-670462 may vary from species to species,32 the results of LD chronic dosing of PF-670462 may also vary from species to species. Testing whether the results of chronic dosing can be predicted from a PRC to light and a PRC to PF-670462 for other species will be a interesting future work. Similarly PRCs to other circadian drugs could be tested and incorporated into a similar model.

We demonstrate in our model that entrainment or the lack of entrainment is predicted to occur in response to LD chronic dosing. Periodically forced clocks can also show more complex behaviors including chaos.40 It would be an interesting future work to find whether these complex behaviors can be seen in the model or in experimental data.

Many studies have shown that bipolar disorders and mood disorders are associated with disrupted circadian rhythms.8,9,10,41 Specifically, phase misalignment is frequently present in mood disorders and the amplitude of circadian rhythms is often reduced.11,41,42 In particular, a recent study observed attenuated amplitudes in Rev-erbα and DBP gene expression in fibroblasts of bipolar disorder patients.43 Interestingly, our study showed that dosing of PF-670462 increases the amplitude of these genes' expression under a LD cycle (Figure 3a). This suggests that the chronic dosing of PF-670462 under a LD cycle could be a potential way to treat mood disorders on the basis of modulating the phase stably and increasing the amplitude of circadian rhythms. Indeed, recent studies showed that the chronic inhibition of CK1δ/ε rescues pathological behaviors in animal models of mania and alcoholism.44,45 Future work should test whether the chronic dosing of PF-670462 can treat these and other models of mood disorders.46,47,48,49,50


Methods
Modeling studies

Estimation of parameters. The values of newly added parameters are estimated by using SA method in two steps.29 In the first round, SA found 10's of parameter sets that simulate time course of PF-670462 in plasma and brain, matching experimental data (Figure 4a). Among these parameter sets, about 10 parameter sets also generated all delayed PRC to 50 mg/kg PF-670462 seen in the experimental data (Figure 4e).16 We used these parameter sets as initial parameter sets for the second round of SA. In the second round, SA found a final parameter set (Supplementary Table S5 online) that matches the chronic-dosing experimental data (Figures 2 and 4f) as well as generates all delayed PRC to 50 mg/kg PF-670462 (Figure 4e) and pharmacokinetic data of PF-670462 (Figure 4a). During SA, we did not fit simulated PRC to the experimentally measured PRC of rats since the original model is based on experimental data of mice. Nevertheless, the simulated PRC with the final parameter set is a good fit with the measured PRC of rats (Figure 4e). This indicates that mice and rats may have similar PRCs to PF-670462.

Simulation. All the simulation and parameter search were done with 150 × 8 Ghz CPU using MATHEMATICA 8.0 (Wolfram Research, Champaign, IL). The code of the model is available in MATHEMATICA format from the ModelDB (Access code: 148320).

Model description. Details can be found in Supplementary Methods online, Supplementary Tables S5–7 online, and Supplementary Figure S2 online.

Experimental studies. Details can be found in the Supplementary Data online and Supplementary Table s8 online.


Conflict of Interest

M.M., A.D., T.W., and C.C. are employees of Pfizer, Inc. The other authors declared no conflict of interest.


Author Contributions

J.K.K., D.B.F., T.W., C.C., and K.M.W. designed the research; J.K.K., D.B.F., M.M., D.W., A.D., and C.C. performed the research; J.K.K., D.B.F., M.M., D.W., A.D., T.W., C.C., and K.M.W. analyzed the data; J.K.K., D.B.F., C.C., and K.M.W. wrote the manuscript; J.K.K., D.B.F., M.M., D.W., A.D., T.W., C.C., and K.M.W. contributed new reagents/analytical tools.


Study Highlights


The authors thank Wendy Adamowicz, Emily Sylvain, Karin Wallace, and Katheleen Wood for excellent technical assistance; and James Offord, Robin Kleiman, Hugh Barton, Tristan Maurer, and Barbara Tate for discussions. This work was supported by the Human Frontiers of Science Program grant (RPG 24/2012), by the National Science Foundation grant DMS-1026317 for computing resources, and by Pfizer grant to University of Michigan, N011779.


References
Antle M.C.,Silver R.. Orchestrating time: arrangements of the brain circadian clock.Trends Neurosci28145151Year: 200515749168
Gallego M.,Virshup D.M.. Post-translational modifications regulate the ticking of the circadian clock.Nat. Rev. Mol. Cell Biol8139148Year: 200717245414
Ko C.H.,Takahashi J.S.. Molecular components of the mammalian circadian clock.Hum. Mol. Genet15 Spec No 2R271R277Year: 200616987893
Lee H.,Chen R.,Lee Y.,Yoo S.,Lee C.. Essential roles of CKIdelta and CKIepsilon in the mammalian circadian clock.Proc. Natl. Acad. Sci. U.S.A1062135921364Year: 200919948962
Dunlap J.C.,Loros J.J.,DeCoursey P.J.. Chronobiology: Biological TimekeepingSinauer Associates, Sunderland, MassachusettsYear: 2004
Oda G.A.,Friesen W.O.. Modeling two-oscillator circadian systems entrained by two environmental cycles.PLoS ONE6e23895Year: 201121886835
Sahar S.,Sassone-Corsi P.. Metabolism and cancer: the circadian clock connection.Nat. Rev. Cancer9886896Year: 200919935677
Kripke D.F.,Mullaney D.J.,Atkinson M.,Wolf S.. Circadian rhythm disorders in manic-depressives.Biol. Psychiatry13335351Year: 1978667233
Wehr T.A.,Turner E.H.,Shimada J.M.,Lowe C.H.,Barker C.,Leibenluft E.. Treatment of rapidly cycling bipolar patient by using extended bed rest and darkness to stabilize the timing and duration of sleep.Biol. Psychiatry43822828Year: 19989611672
Lenox R.H.,Gould T.D.,Manji H.K.. Endophenotypes in bipolar disorder.Am. J. Med. Genet114391406Year: 200211992561
McClung C.A.. Circadian rhythms and mood regulation: insights from pre-clinical models.Eur. Neuropsychopharmacol21suppl. 4S683S693Year: 201121835596
Hirota T.,Lewis W.G.,Liu A.C.,Lee J.W.,Schultz P.G.,Kay S.A.. A chemical biology approach reveals period shortening of the mammalian circadian clock by specific inhibition of GSK-3beta.Proc. Natl. Acad. Sci. U.S.A1052074620751Year: 200819104043
Pandi-Perumal S.R.,Trakht I.,Spence D.W.,Srinivasan V.,Dagan Y.,Cardinali D.P.. The roles of melatonin and light in the pathophysiology and treatment of circadian rhythm sleep disorders.Nat. Clin. Pract. Neurol4436447Year: 200818628753
Isojima Y.,et al. CKIepsilon/delta-dependent phosphorylation is a temperature-insensitive, period-determining process in the mammalian circadian clock.Proc. Natl. Acad. Sci. U.S.A1061574415749Year: 200919805222
Hirota T.,et al. Identification of small molecule activators of cryptochrome.Science33710941097Year: 201222798407
Badura L.,et al. An inhibitor of casein kinase I epsilon induces phase delays in circadian rhythms under free-running and entrained conditions.J. Pharmacol. Exp. Ther322730738Year: 200717502429
Walton K.M.,et al. Selective inhibition of casein kinase 1 epsilon minimally alters circadian clock period.J. Pharmacol. Exp. Ther330430439Year: 200919458106
Meng Q.J.,et al. Entrainment of disrupted circadian behavior through inhibition of casein kinase 1 (CK1) enzymes.Proc. Natl. Acad. Sci. U.S.A1071524015245Year: 201020696890
Reppert S.M.,Weaver D.R.. Coordination of circadian timing in mammals.Nature418935941Year: 200212198538
Wilsbacher L.D.,et al. Photic and circadian expression of luciferase in mPeriod1-luc transgenic mice in vivo.Proc. Natl. Acad. Sci. U.S.A99489494Year: 200211752392
Yan L.,Takekida S.,Shigeyoshi Y.,Okamura H.. Per1 and Per2 gene expression in the rat suprachiasmatic nucleus: circadian profile and the compartment-specific response to light.Neuroscience94141150Year: 199910613504
Forger D.B.,Peskin C.S.. A detailed predictive model of the mammalian circadian clock.Proc. Natl. Acad. Sci. U.S.A1001480614811Year: 200314657377
Gallego M.,Eide E.J.,Woolf M.F.,Virshup D.M.,Forger D.B.. An opposite role for tau in circadian rhythms revealed by mathematical modeling.Proc. Natl. Acad. Sci. U.S.A1031061810623Year: 200616818876
Kim J.K.,Forger D.B.. A mechanism for robust circadian timekeeping via stoichiometric balance.Mol. Syst. Biol8630Year: 201223212247
Sprouse J.,Reynolds L.,Kleiman R.,Tate B.,Swanson T.A.,Pickard G.E.. Chronic treatment with a selective inhibitor of casein kinase I delta/epsilon yields cumulative phase delays in circadian rhythms.Psychopharmacology (Berl.)210569576Year: 201020407760
Ko C.H.,et al. Emergence of noise-induced oscillations in the central circadian pacemaker.PLoS Biol8e1000513Year: 201020967239
Kiessling S.,Eichele G.,Oster H.. Adrenal glucocorticoids have a key role in circadian resynchronization in a mouse model of jet lag.J. Clin. Invest12026002609Year: 201020577050
Lévi F.,Altinok A.,Clairambault J.,Goldbeter A.. Implications of circadian clocks for the rhythmic delivery of cancer therapeutics.Philos. Trans. A. Math. Phys. Eng. Sci36635753598Year: 200818644767
Gonzalez O.R.,Küper C.,Jung K.,Naval P.C. Jr,Mendoza E.. Parameter estimation using Simulated Annealing for S-system models of biochemical networks.Bioinformatics23480486Year: 200717038344
Winfree A.T.. The Geometry of Biological TimeSpringer Verlag, New YorkYear: 1980
Glass L.,Winfree A.T.. Discontinuities in phase-resetting experiments.Am. J. Physiol246R251R258Year: 19846696148
Daan S.,Pittendrigh C.S.. Functional-analysis of circadian pacemakers in nocturnal rodents 2: variability of phase response curves.J. Comp. Physiol106253266Year: 1976
Kronauer R.E.,Forger D.B.,Jewett M.E.. Quantifying human circadian pacemaker response to brief, extended, and repeated light stimuli over the phototopic range.J. Biol. Rhythms14500515Year: 199910643747
Comas M.,Beersma D.G.,Spoelstra K.,Daan S.. Phase and period responses of the circadian system of mice (Mus musculus) to light stimuli of different duration.J. Biol. Rhythms21362372Year: 200616998156
Morgenthaler T.I.,Standards of Practice Committee of the American Academy of Sleep Medicineet al. Practice parameters for the clinical evaluation and treatment of circadian rhythm sleep disorders. An American Academy of Sleep Medicine report.Sleep3014451459Year: 200718041479
Gooley J.J.. Treatment of circadian rhythm sleep disorders with light.Ann. Acad. Med. Singap37669676Year: 200818797560
Barion A.. Circadian rhythm sleep disorders.Dis. Mon57423437Year: 201121929966
Price M.A.. CKI, there's more than one: casein kinase I family members in Wnt and Hedgehog signaling.Genes Dev20399410Year: 200616481469
Cheong J.K.,et al. IC261 induces cell cycle arrest and apoptosis of human cancer cells via CK1d/ε and Wnt/ß-catenin independent inhibition of mitotic spindle formation.Oncogene3025582569Year: 201121258417
Pikovsky A.,Rosenblum M.,Kurths J.. Synchronization: A Universal Concept in Nonlinear Sciences1st paperback edn. (Cambridge University Press, CambridgeYear: 2003
McCarthy M.J.,Welsh D.K.. Cellular circadian clocks in mood disorders.J. Biol. Rhythms27339352Year: 201223010657
Kennaway D.J.. Clock genes at the heart of depression.J. Psychopharmacol. (Oxford)24514Year: 201020663803
Yang S.,Van Dongen H.P.,Wang K.,Berrettini W.,Bucan M.. Assessment of circadian function in fibroblasts of patients with bipolar disorder.Mol. Psychiatry14143155Year: 200918301395
Arey R.,McClung C.A.. An inhibitor of casein kinase 1 epsilon/delta partially normalizes the manic-like behaviors of the ClockΔ19 mouse.Behav. Pharmacol23392396Year: 201222743604
Perreau-Lenz S.,et al. Inhibition of the casein-kinase-1-δ/ε prevents relapse-like alcohol drinking.Neuropsychopharmacology3721212131Year: 201222549116
Maier S.F.. Learned helplessness and animal models of depression.Prog. Neuropsychopharmacol. Biol. Psychiatry8435446Year: 19846385140
Spanagel R.,et al. The clock gene Per2 influences the glutamatergic system and modulates alcohol consumption.Nat. Med113542Year: 200515608650
Prickaerts J.,et al. Transgenic mice overexpressing glycogen synthase kinase 3beta: a putative model of hyperactivity and mania.J. Neurosci2690229029Year: 200616943560
Roybal K.,et al. Mania-like behavior induced by disruption of CLOCK.Proc. Natl. Acad. Sci. U.S.A10464066411Year: 200717379666
Hampp G.,et al. Regulation of monoamine oxidase A by circadian-clock components implies clock influence on mood.Curr. Biol18678683Year: 200818439826
Supplementary Material Supplementary Figure S1 Click here for additional data file (psp201334x1.png)

Supplementary Figure S2 Click here for additional data file (psp201334x2.png)

Supplementary Figure S3 Click here for additional data file (psp201334x3.png)

Supplementary Figure S4 Click here for additional data file (psp201334x4.png)

Supplementary Methods Supplementary Table S1 Supplementary Table S2 Supplementary Table S3 Supplementary Table S4 Supplementary Table S5 Supplementary Table S6 Supplementary Table S7 Supplementary Table S8

Article Categories:
  • Original Article


Previous Document:  Influence of a priori Information, Designs, and Undetectable Data on Individual Parameters Estimatio...
Next Document:  The impact of lobular and ductal breast cancer histology on the metastatic behavior and long term su...