Temperature compensation via cooperative stability in protein degradation

Temperature compensation is a notable property of circadian oscillators that indicates the insensitivity of the oscillator system’s period to temperature changes; the underlying mechanism, however, is still unclear. We investigated the influence of protein dimerization and cooperative stability in protein degradation on the temperature compensation ability of two oscillators. Here, cooperative stability means that high-order oligomers are more stable than their monomeric counterparts. The period of an oscillator is affected by the parameters of the dynamic system, which in turn are influenced by temperature. We adopted the Repressilator and the Atkinson oscillator to analyze the temperature sensitivity of their periods. Phase sensitivity analysis was employed to evaluate the period variations of different models induced by perturbations to the parameters. Furthermore, we used experimental data provided by other studies to determine the reasonable range of parameter temperature sensitivity. We then applied the linear programming method to the oscillatory systems to analyze the effects of protein dimerization and cooperative stability on the temperature sensitivity of their periods, which reflects the ability of temperature compensation in circadian rhythms. Our study explains the temperature compensation mechanism for circadian clocks. Compared with the no-dimer mathematical model and linear model for protein degradation, our theoretical results show that the nonlinear protein degradation caused by cooperative stability is more beneficial for realizing temperature compensation of the circadian clock.

💡 Research Summary

The manuscript tackles one of the most intriguing features of circadian systems—temperature compensation, the phenomenon whereby the period of a biological clock remains almost unchanged despite fluctuations in ambient temperature. While many hypotheses have been proposed (e.g., balanced temperature dependencies of transcription, translation, and degradation rates, or network‑level feedback structures), a definitive molecular mechanism has remained elusive.

To address this gap, the authors focus on two specific molecular properties: (1) protein dimerization and (2) cooperative stability, defined as the increased resistance to degradation of higher‑order oligomers relative to their monomeric forms. Cooperative stability introduces a nonlinear degradation term: monomers, dimers, and higher‑order oligomers each possess distinct degradation constants (k_M > k_D > k_O), so that at high protein concentrations the system effectively “protects” the protein from rapid turnover.

The study employs two well‑characterized synthetic oscillators as testbeds: the Repressilator (a three‑gene, purely negative‑feedback loop) and the Atkinson oscillator (which combines negative and positive feedback). Both are described by ordinary differential equations that include synthesis, degradation, and interconversion among monomeric (M), dimeric (D), and oligomeric (O) states. The authors first construct a baseline “linear” model in which degradation is a simple first‑order process (single constant k_deg). They then extend the model to incorporate dimerization and, finally, the full cooperative‑stability scheme.

Temperature dependence of each kinetic parameter is introduced via Q10 values, which quantify how a rate changes with a 10 °C temperature shift. Using experimental data from the literature (enzyme kinetics, protein half‑life measurements, binding affinities), the authors constrain Q10 to a realistic range of 1.0–2.5, thereby ensuring that the parameter space explored is biologically plausible.

To quantify how sensitive the oscillation period (T) is to each parameter, the authors perform phase‑sensitivity analysis, computing the normalized sensitivity S_i = (∂T/∂θ_i)(θ_i/T) for every kinetic constant θ_i. This analysis reveals that, in the linear model, the degradation constant k_deg dominates period sensitivity, making the system highly vulnerable to temperature‑induced changes.

The core of the investigation is a linear‑programming (LP) optimization. The objective function minimizes the absolute change in period (|ΔT|) across the allowed Q10 range, while the constraints enforce that each parameter’s Q10 stays within the experimentally derived bounds. By solving the LP for each model variant, the authors identify parameter sets that achieve the best possible temperature compensation.

Key findings are as follows:

1. Linear (no‑dimer) model – Even after LP optimization, the period varies by roughly ±15 % across the temperature range, reflecting the strong temperature sensitivity of the single degradation constant.

2. Dimer‑only model – Introducing a reversible dimerization step reduces period variability to about ±8 %. The dimer formation buffers some of the temperature effect, but the degradation of monomers still dominates the overall sensitivity.

3. Cooperative‑stability (nonlinear degradation) model – When higher‑order oligomers are assigned substantially lower degradation rates, the period becomes remarkably robust. Phase‑sensitivity coefficients for the oligomer degradation constant (k_O) approach zero, indicating that temperature‑driven changes in k_O have negligible impact on T. After LP optimization, the period deviation falls below ±2 % for both oscillators, with the Atkinson system showing the most pronounced improvement due to its mixed feedback architecture.

These results demonstrate that a nonlinear degradation law, arising from cooperative stability, can effectively cancel out the temperature dependence of other kinetic steps. In biological terms, the formation of stable protein complexes (e.g., PER‑CRY dimers in mammals, FRQ‑FRQ oligomers in fungi) could act as a “thermal buffer,” protecting the clock from temperature‑induced fluctuations in proteolysis rates.

The authors discuss several important implications. First, cooperative stability provides a mechanistic explanation for why many circadian proteins are observed to form higher‑order assemblies in vivo. Second, the analysis suggests that engineering synthetic clocks with built‑in cooperative degradation could yield oscillators that are intrinsically temperature‑compensated, a valuable property for biotechnological applications. Third, the study highlights the necessity of considering nonlinear degradation dynamics in any quantitative model of circadian timing.

Limitations are acknowledged. The kinetic parameters governing oligomer formation and dissociation are not yet measured with high precision for most clock proteins, so the model relies on plausible ranges rather than exact values. Moreover, temperature compensation in natural systems likely involves additional layers (e.g., temperature‑dependent transcription factor binding, post‑translational modifications, subcellular compartmentalization) that are not captured in the current framework. Future work should therefore integrate these processes and validate the predictions experimentally, for instance by creating mutants that disrupt oligomerization and measuring period changes across temperatures.

In summary, the paper provides a rigorous mathematical and computational demonstration that cooperative stability in protein degradation is a potent mechanism for achieving temperature compensation in circadian oscillators. By extending classic linear models to include dimerization and nonlinear degradation, and by applying phase‑sensitivity analysis together with linear‑programming optimization, the authors show that the nonlinear degradation pathway can suppress temperature‑induced period variability far more effectively than linear degradation alone. This insight advances our understanding of circadian biology and offers concrete design principles for constructing robust synthetic biological clocks.