A dual-feedback loop model of the mammalian circadian clock for multi-input control of circadian phase.

The molecular circadian clock is driven by interlocked transcriptional-translational feedback loops, producing oscillations in the expressions of genes and proteins to coordinate the timing of biological processes throughout the body. Modeling this system gives insight into the underlying processes driving oscillations in an activator-repressor architecture and allows us to make predictions about how to manipulate these oscillations. The knockdown or upregulation of different cellular components using small molecules can disrupt these rhythms, causing a phase shift, and we aim to determine the dosing of such molecules with a model-based control strategy. Mathematical models allow us to predict the phase response of the circadian clock to these interventions and time them appropriately but only if the model has enough physiological detail to describe these responses while maintaining enough simplicity for online optimization. We build a control-relevant, physiologically-based model of the two main feedback loops of the mammalian molecular clock, which provides sufficient detail to consider multi-input control. Our model captures experimentally observed peak to trough ratios, relative abundances, and phase differences in the model species, and we independently validate this model by showing that the in silico model reproduces much of the behavior that is observed in vitro under genetic knockout conditions. Because our model produces valid phase responses, it can be used in a model predictive control algorithm to determine inputs to shift phase. Our model allows us to consider multi-input control through small molecules that act on both feedback loops, and we find that changes to the parameters of the negative feedback loop are much stronger inputs for shifting phase. The strongest inputs predicted by this model provide targets for new experimental small molecules and suggest that the function of the positive feedback loop is to stabilize the oscillations while linking the circadian system to other clock-controlled processes.