Quantification of protein abundance and interaction defines a mechanism for operation of the circadian clock

  1. Alex A Koch
  2. James S Bagnall
  3. Nicola J Smyllie
  4. Nicola Begley
  5. Antony D Adamson
  6. Jennifer L Fribourgh
  7. David G Spiller
  8. Qing-Jun Meng
  9. Carrie L Partch
  10. Korbinian Strimmer
  11. Thomas A House
  12. Michael H Hastings
  13. Andrew SI Loudon  Is a corresponding author
  1. Faculty of Biology, Medicine and Health, University of Manchester, United Kingdom
  2. MRC Laboratory of Molecular Biology, United Kingdom
  3. Department of Chemistry and Biochemistry, University of California, Santa Cruz, United States
  4. Department of Mathematics, University of Manchester, United Kingdom
7 figures, 4 tables and 1 additional file

Figures

Figure 1 with 2 supplements
Short-lived DNA binding of BMAL1 and CLOCK.

(A) Schematic representation of parameters regulating CLOCK:BMAL1 dimers binding to target DNA sites. (B) NIH/3T3 cells are either singularly or sequentially transduced to express fluorescent …

Figure 1—figure supplement 1
Ectopically expressed mRNA is the major form in a lentivirus transduced system.

(A) NIH/3T3 cells were stained for CLOCK mRNA by single-molecule fluorescent in situ hybridisation (smFISH). (B) Mature mRNA was counted from images of many single cells to determine the mRNA …

Figure 1—figure supplement 2
Binding plays a significant role in BMAL1 mobility.

Fluorescence recovery after photobleaching for NIH/3T3 cells transduced with tagRFP::CLOCK and BMAL1:EGFP. (A) Imaging protocol was performed on BMAL1::EGFP signal. Regions of photobleaching are …

Figure 2 with 2 supplements
Live-cell interaction measurements demonstrate BMAL1 and CLOCK mobility is regulated by dimerisation and DNA binding.

(A) Schematic of confocal volume used in FCCS with corresponding photon count traces. Interaction may be seen by correlation between both channels. Representative auto- and cross- correlation data …

Figure 2—figure supplement 1
Anomalous diffusion best fits protein movement.

(A) Representative normal and anomalous model fits for BMAL1::EGFP FCS data sets for cells transduced with lentivirus to express BMAL1, CLOCK or control fluorescent proteins. (B) Summary data sets …

Figure 2—figure supplement 2
Fluorescent BMAL1 and CLOCK proteins behave similarly when colours are swapped.

(A) Confocal images of NIH/3T3-LV2 EGFP::CLOCK-BMAL1::tagRFP cells. (B) Average auto- and cross- correlation curves shown as mean (line) and standard deviation (error envelope) for NIH/3T3 cells …

Figure 3 with 3 supplements
A rhythmic and strong interaction observed between slow-diffusing BMAL1 and CRY1 facilitates repression.

(A) Schematic of triple-labelled mice from which primary lung fibroblasts were isolated (B) Confocal images of two cells shown for Venus::BMAL1 and CRY1::mRuby3 over time. FCS determined measurement …

Figure 3—figure supplement 1
BMAL1 concentration and DNA binding parameters minimally vary across cell types.

(A) Confocal microscopy images of primary cultures isolated from Venus::BMAL1 mice. (B) Characteristic bound time calculated from FRAP measurements of primary cultures from A. Chondrocyte data was …

Figure 3—figure supplement 2
Generation of CRY1::mRuby3 mouse line.

(A) Genotyping of pups by PCR. Images are PCR reactions run on Qiaxcel with red arrows indicating correct HDR product size, blue asterisks indicate mice in which InDels resulting from NHEJ are …

Figure 3—figure supplement 3
Triple endogenous labelled mice used to assay rhythms in SCN and peripheral lung fibroblasts.

(A) Schematic representation of newly made transgenic mouse engineered to express CRY1::mRuby3, Venus::BMAL1, and PER2::LUC. (B) Confocal microscopy image of SCN organotypic slice expressing …

Figure 4 with 1 supplement
PER2 modulates CRY1 mobility via a high-affinity association.

(A) Confocal images of transduced NIH/3T3 cells that either solo- or co- express PER2 and CRY1. (B) FCS data showing diffusion for PER2 and CRY1 in solo- and co-expressed conditions (n = 165, 174, …

Figure 4—figure supplement 1
CRY1 mobility is affected by co-expression with PER2.

(A) Average auto- and cross- correlation curves shown as mean (line) and standard deviation (error envelope) for NIH/3T3 cells transduced to express EGFP::PER2 and CRY1::tagRFP. Measurements were …

Figure 5 with 1 supplement
PER2 acts via CRY1 to mediate rhythmic displacement of CLOCK:BMAL1 from DNA.

(A) Schematic representation of model topology used for the deterministic model of CLOCK:BMAL1 DNA binding. (B) Primary lung fibroblasts from BMAL1 x CRY1 x PER2 mice were synchronised with …

Figure 5—figure supplement 1
ODE model of CLOCK:BMAL1 DNA binding using measured inputs and modelled perturbations.

(A) Normalised model simulated FRAP using ODE model by counting recovery of site bound BMAL1 species after removal. (B) Model determined residence time of DNA bound CLOCK:BMAL1 across different …

Circadian proteins operate within an optimal range to modulate E-Box binding.

Sensitivity analysis of the deterministic binding model showing relationship of measured parameters (bottom) against model for occupancy of active BMAL1:CLOCK on target sites (top). (A) Changing …

Figure 7 with 1 supplement
Mathematical modelling demonstrates dual function of PER:CRY mediated repression.

Stochastic binding model outputs using parameters corresponding to T28, T32 or T40 post dexamethasone BMAL1 x CRY1 data sets. (A) Shows a promoter corresponding to the average binding rate of …

Figure 7—figure supplement 1
Stochastic binding model using experimentally measured parameters (A) Stochastic model showing the average binding (with SD) of CLOCK:BMAL1 bound target sites using input measurements from all time points for both WT and without PER2 simulations.

(B) The time to visit every E-Box site once for T28 showing fit. (C) Model simulation plots showing CLOCK:BMAL1 visits to a single promoter over time. Red and blue lines show simulations using …

Tables

Table 1
Summary of ordinary differential equation model parameters.

Model fit χ2=7.46.

Input parameters
ParameterUnitDescriptionValue±SD
KD(C:B - E-Box)nMCLOCK:BMAL1 - E-Box dissociation constant10 (Huang et al., 2012)
KD(C-B)nMCLOCK - BMAL1 dissociation constant147.6 ± 9.8
KD(B-C1)nMBMAL1 - CRY1 dissociation constantTime-point dependent, Figure 3F
KD(C1-P2)nMCRY1 - PER2 dissociation constant81.8 ± 4.9
Fitted parameters
ParameterUnitDescriptionValue±SD (Inverse Hessian eigenvalue of fit)
kONnm-1s-1CLOCK:BMAL1 - DNA binding on rate0.027±1.034(1.96)
dONnm-1s-1BMAL1 - CRY1 binding rate0.237±1.003(3.42×10-2)
aONnm-1s-1PER2 - CRY1 binding rate6.34±1.00(6.79×10-8)
ROFFs-1CLOCK:BMAL1:CRY1:PER2 - DNA unbinding rate(1.23±0.33)×101(1.00)
bONnm-1s-1CLOCK - BMAL1 binding rate9.17±1.29(1.00)
Derived parameters
ParameterUnitDescriptionValue±SD
kOFFs-1CLOCK:BMAL1 - DNA binding unbinding ratekON×KD(C:B - E-Box)=0.27±10.34
bOFFs-1CLOCK - BMAL1 unbinding ratebON×KD(C-B)=(1.35±0.21)×103
dOFFnm-1s-1BMAL1 - CRY1 unbinding rateTime-point dependent,dON×KD(B:C1)
aOFFs-1PER2 - CRY1 unbinding rateaON×KD(C1-P2)=(5.19±0.88)×102
Table 2
BMAL1 ChIP reports.
No.TissueBMAL1 peaksReference
1Liver2049Rey et al., 2011
2Liver5952Koike et al., 2012
3U2OS2001Wu et al., 2017
4PECS2026Oishi et al., 2017
5Liver4813Beytebiere et al., 2019
6Kidney4034Beytebiere et al., 2019
7Heart2520Beytebiere et al., 2019
8NIH3T34740Chiou et al., 2016
9Skeletal muscle2787Dyar et al., 2018
Mean average3436
Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Genetic reagent (M. musclus)C57BL/6 Venus::BMAL1Yang et al., 2020Venus sequence inserted before BMAL1 start codon.
Genetic reagent (M. musclus)C57BL/6 Cry1::mRuby3This paperCRY1 stop codon replaced with mRuby3
Genetic reagent (M. musclus)C57BL/6 Venus::BMAL1 x CRY1::mRuby3This paperCrossed from Venus::BMAL1 and CRY1::mRuby3 mice
Cell line (M. musculus)NIH/3T3ATCCCRL-1658
Transfected construct (M. musculus)pLNT-NLS::EGFPVector Builder VB900119-0501njqLentiviral construct to express nuclear EGFP.
Transfected construct (M. musculus)pLNT-BMAL1::EGFP or pLNT-BMAL1::RFPThis paperNCBI reference: NM_007489.4Lentiviral construct to express fluorescent BMAL1.
Transfected construct (M. musculus)pLNT-BMAL1-L95E::EGFPThis paperNCBI reference: NM_007489.4Lentiviral construct to express fluorescent BMAL1 L95E mutant.
Transfected construct (M. musculus)pLNT-BMAL1-V435R::EGFPThis paperNCBI reference: NM_007489.4Lentiviral construct to express fluorescent BMAL1 V435R mutant.
Transfected construct (M. musculus)pLNT-EGFP::CLOCK or pLNT-RFP::CLOCKThis paperNCBI reference: NM_007715.6Lentiviral construct to express fluorescent CLOCK.
Transfected construct (M. musculus)pLNT-EGFP::PER2This paperNCBI reference: NM_011066Lentiviral construct to express fluorescent PER2.
Transfected construct (M. musculus)pLNT-CRY1::RFPThis paperNCBI reference: NM_007771.3Lentiviral construct to express fluorescent CRY1.
Chemical compound, drugDexamethasoneSigma AldrichD4902
Software, algorithmGraphPad PrismGraphPad PrismVersion 9
Software, algorithmFCCS analysis pipelineThis paperhttps://github.com/LoudonLab/FcsAnalysisPipeline,(copy archived at swh:1:rev:b12e9007ed7f8a033485e57c8605e27c67df74f1; Koch, 2021)
Table 3
Stochastic model reactions and propensities.

Counter for arrivals by CLOCK:BMAL1 (CB) without CRY1 (C1) to previously unbound sites S converting them to S0 given by ACB as well as counters for marked site M binding represented by BX, and …

No.ReactionPropensity
1CB+SCBS+ACB(kON/Ω)CBS
2CB+S0CBS(kON/Ω)CBS0
3CBSCB+S0kOFFCBS
4CBC1+SCBC1S(kON/Ω)CBC1S
5CBC1SCBC1+SkOFFCBC1S
6CBC1+S0CBC1S0(kON/Ω)CBC1S0
7CBC1S0CBC1+S0kOFFCBC1S0
8CB+MCBM+BCB(kON/Ω)CBM
9CBMCB+M+UCBkOFFCBM
10CBC1+MCBC1M+BCBC1(kON/Ω)CBC1M
11CBC1MCBC1+M+UCBC1kOFFCBC1M

Additional files

Download links