Sleep deprivation in a quantitative physiologically based model of the ascending arousal system
Catherine
Lloyd
Bioengineering Institute, University of Auckland
Model Status
This CellML model replicates the "baseline" version of the original Phillips and Robinson 2008 model. It runs in both PCEnv and COR to replicate figure 2 in the paper (using the old nominal values from table one). The units have been checked and they are consistent. Note that the model has to be run for at least 3 cycles for it to stabilse into steady oscillations.
Model Structure
Abstract: A physiologically based quantitative model of the human ascending arousal system is used to study sleep deprivation after being calibrated on a small set of experimentally based criteria. The model includes the sleep-wake switch of mutual inhibition between nuclei which use monoaminergic neuromodulators, and the ventrolateral preoptic area. The system is driven by the circadian rhythm and sleep homeostasis. We use a small number of experimentally derived criteria to calibrate the model for sleep deprivation, then investigate model predictions for other experiments, demonstrating the scope of application. Calibration gives an improved parameter set, in which the form of the homeostatic drive is better constrained, and its weighting relative to the circadian drive is increased. Within the newly constrained parameter ranges, the model predicts repayment of sleep debt consistent with experiment in both quantity and distribution, asymptoting to a maximum repayment for very long deprivations. Recovery is found to depend on circadian phase, and the model predicts that it is most efficient to recover during normal sleeping phases of the circadian cycle, in terms of the amount of recovery sleep required. The form of the homeostatic drive suggests that periods of wake during recovery from sleep deprivation are phases of relative recovery, in the sense that the homeostatic drive continues to converge toward baseline levels. This undermines the concept of sleep debt, and is in agreement with experimentally restricted recovery protocols. Finally, we compare our model to the two-process model, and demonstrate the power of physiologically based modeling by correctly predicting sleep latency times following deprivation from experimental data.
The original paper reference is cited below:
Sleep deprivation in a quantitative physiologically based model of the ascending arousal system, A.J.K. Phillips and P.A. Robinson, 2008,
Journal of Theoretical Biology, 255, 413-423. PubMed ID: 18805427
model diagram
Schematic diagram of the Phillips and Robinson 2008 sleep model. MA is the monoaminergic group, ACh is acetylcholine, VLPO represents the ventrolateral preoptic area, and D is the drive input which consists of circadian and homeostatic components.
$\mathrm{Qv}=\frac{\mathrm{Qmax}}{1+e^{\frac{-(\mathrm{Vv}-\mathrm{theta})}{\mathrm{sigma}}}}$
$\frac{d \mathrm{Vv}}{d \mathrm{time}}=\frac{\mathrm{v\_vm}\mathrm{Qm}+D-\mathrm{Vv}}{\frac{\mathrm{tau\_v}}{3600}}$
$\mathrm{Qm}=\frac{\mathrm{Qmax}}{1+e^{\frac{-(\mathrm{Vm}-\mathrm{theta})}{\mathrm{sigma}}}}$
$\frac{d \mathrm{Vm}}{d \mathrm{time}}=\frac{A+\mathrm{v\_mv}\mathrm{Qv}-\mathrm{Vm}}{\frac{\mathrm{tau\_m}}{3600}}$
$\frac{d H}{d \mathrm{time}}=\frac{\mathrm{mu}\mathrm{Qm}-H}{\mathrm{chi}}$
$D=\mathrm{v\_vc}C+\mathrm{v\_vh}H$
$C=\mathrm{c0}+\cos (\mathrm{omega}(\mathrm{time}-\mathrm{alpha}))$
$\mathrm{omega}=\frac{2\pi}{24}$
CatherineLloydkeyword400neurobiologysleepdiurnal cyclecircadian rhythms2009-10-29Auckland Bioengineering Institute
This is the CellML description of the Phillips and Robinson 2008 quantitative model of sleep deprivation
18805427c.lloyd@auckland.ac.nzSleep deprivation in a quantitative physiologically based model of the ascending arousal system (Baseline Model)PARobinson2008-12-21Catherine LloydAKJPhillips
The Phillips and Robinson 2008 quantitative model of sleep deprivation
Sleep deprivation in a quantitative physiologically based model of the ascending arousal system
423413255Journal of Theoretical Biology