- Author:
- Catherine Lloyd <c.lloyd@auckland.ac.nz>
- Date:
- 2009-11-05 17:39:00+13:00
- Desc:
- Added model files to the new workspace.
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/2ed/rawfile/49397de30eac1520542c630afc344e430d0f6278/phillips_sleep_deprivation_2008.cellml
<?xml version='1.0' encoding='utf-8'?>
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="phillips_2008" name="phillips_2008">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Sleep deprivation in a quantitative physiologically based model of the ascending arousal system</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<sect1 id="sec_status">
<title>Model Status</title>
<para>
This CellML model replicates the sleep deprivation version of the original Phillips and Robinson 2008 model. It runs in both PCEnv and COR to replicate figure 5 in the paper (using the new 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.
</para>
</sect1>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
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.
</para>
<para>
The original paper reference is cited below:
</para>
<para>
Sleep deprivation in a quantitative physiologically based model of the ascending arousal system, A.J.K. Phillips and P.A. Robinson, 2008,
<emphasis>Journal of Theoretical Biology</emphasis>, 255, 413-423. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&cmd=Retrieve&dopt=AbstractPlus&list_uids=18805427&query_hl=1&itool=pubmed_docsum">PubMed ID: 18805427</ulink>
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="phillips_2008.png"/>
</imageobject>
</mediaobject>
<caption>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.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<units name="mV">
<unit units="volt" prefix="milli"/>
</units>
<units name="mV_second">
<unit units="mV"/>
<unit units="second"/>
</units>
<units name="mV_per_second">
<unit units="mV"/>
<unit units="second" exponent="-1"/>
</units>
<units name="mV_per_hour">
<unit units="mV"/>
<unit units="hour" exponent="-1"/>
</units>
<units name="mV_per_nM">
<unit units="mV"/>
<unit units="nM" exponent="-1"/>
</units>
<units name="per_second">
<unit units="second" exponent="-1"/>
</units>
<units name="hour">
<unit units="second" multiplier="3600"/>
</units>
<units name="per_hour">
<unit units="hour" exponent="-1"/>
</units>
<units name="second_per_hour">
<unit units="second"/>
<unit units="hour" exponent="-1"/>
</units>
<units name="nM">
<unit units="mole" prefix="nano"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="nM_second">
<unit units="nM"/>
<unit units="second"/>
</units>
<component name="environment">
<variable units="hour" public_interface="out" name="time"/>
</component>
<component name="v">
<variable units="per_second" public_interface="out" name="Qv"/>
<variable units="mV" public_interface="out" name="Vv" initial_value="-12.5"/>
<variable units="second" name="tau_v" initial_value="10.0"/>
<variable units="mV_second" name="v_vm" initial_value="-2.1"/>
<variable units="hour" name="awake_start" initial_value="167.04"/>
<variable units="hour" name="awake_end" initial_value="240.0"/>
<variable units="per_second" public_interface="in" name="Qm"/>
<variable units="per_second" public_interface="in" name="Qmax"/>
<variable units="mV" public_interface="in" name="D"/>
<variable units="mV" public_interface="in" name="theta"/>
<variable units="mV" public_interface="in" name="sigma"/>
<variable units="hour" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="1">
<apply>
<eq/>
<ci>Qv</ci>
<apply>
<divide/>
<ci>Qmax</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vv</ci>
<ci>theta</ci>
</apply>
</apply>
<ci>sigma</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="2">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Vv</ci>
</apply>
<piecewise>
<piece>
<cn cellml:units="mV_per_hour">0.0</cn>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<ci>awake_start</ci>
</apply>
<apply>
<lt/>
<ci>time</ci>
<ci>awake_end</ci>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>v_vm</ci>
<ci>Qm</ci>
</apply>
<ci>D</ci>
</apply>
<ci>Vv</ci>
</apply>
<apply>
<divide/>
<ci>tau_v</ci>
<cn cellml:units="second_per_hour">3600</cn>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="m">
<variable units="per_second" public_interface="out" name="Qm"/>
<variable units="mV" public_interface="out" name="Vm" initial_value="0.0"/>
<variable units="second" name="tau_m" initial_value="10.0"/>
<variable units="mV_second" name="v_mv" initial_value="-1.8"/>
<variable units="mV" name="A" initial_value="1.3"/>
<variable units="hour" name="awake_start" initial_value="160.0"/>
<variable units="hour" name="awake_end" initial_value="240.0"/>
<variable units="per_second" public_interface="in" name="Qmax"/>
<variable units="mV" public_interface="in" name="theta"/>
<variable units="mV" public_interface="in" name="sigma"/>
<variable units="per_second" public_interface="in" name="Qv"/>
<variable units="hour" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="3">
<apply>
<eq/>
<ci>Qm</ci>
<apply>
<divide/>
<ci>Qmax</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>Vm</ci>
<ci>theta</ci>
</apply>
</apply>
<ci>sigma</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="4">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Vm</ci>
</apply>
<piecewise>
<piece>
<cn cellml:units="mV_per_hour">0.0</cn>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<ci>awake_start</ci>
</apply>
<apply>
<lt/>
<ci>time</ci>
<ci>awake_end</ci>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>A</ci>
<apply>
<times/>
<ci>v_mv</ci>
<ci>Qv</ci>
</apply>
</apply>
<ci>Vm</ci>
</apply>
<apply>
<divide/>
<ci>tau_m</ci>
<cn cellml:units="second_per_hour">3600</cn>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="H">
<variable units="nM" public_interface="out" name="H" initial_value="15.0"/>
<variable units="hour" name="chi" initial_value="45.0"/>
<variable units="nM_second" name="mu" initial_value="4.4"/>
<variable units="per_second" public_interface="in" name="Qm"/>
<variable units="hour" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="5">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>H</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>mu</ci>
<ci>Qm</ci>
</apply>
<ci>H</ci>
</apply>
<ci>chi</ci>
</apply>
</apply>
</math>
</component>
<component name="D">
<variable units="mV" public_interface="out" name="D"/>
<variable units="dimensionless" name="C"/>
<variable units="dimensionless" name="c0" initial_value="4.5"/>
<variable units="per_hour" name="omega"/>
<variable units="hour" name="alpha" initial_value="0"/>
<variable units="mV" name="v_vc" initial_value="-2.9"/>
<variable units="mV_per_nM" name="v_vh" initial_value="1.0"/>
<variable units="nM" public_interface="in" name="H"/>
<variable units="hour" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="6">
<apply>
<eq/>
<ci>D</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>v_vc</ci>
<ci>C</ci>
</apply>
<apply>
<times/>
<ci>v_vh</ci>
<ci>H</ci>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="7">
<apply>
<eq/>
<ci>C</ci>
<apply>
<plus/>
<ci>c0</ci>
<apply>
<cos/>
<apply>
<times/>
<ci>omega</ci>
<apply>
<minus/>
<ci>time</ci>
<ci>alpha</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="8">
<apply>
<eq/>
<ci>omega</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<pi/>
</apply>
<cn cellml:units="hour">24</cn>
</apply>
</apply>
</math>
</component>
<component name="model_parameters">
<variable units="per_second" public_interface="out" name="Qmax" initial_value="100.0"/>
<variable units="mV" public_interface="out" name="theta" initial_value="10.0"/>
<variable units="mV" public_interface="out" name="sigma" initial_value="3.0"/>
</component>
<connection>
<map_components component_2="environment" component_1="v"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="m"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="H"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="D"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="v" component_1="D"/>
<map_variables variable_2="D" variable_1="D"/>
</connection>
<connection>
<map_components component_2="v" component_1="m"/>
<map_variables variable_2="Qm" variable_1="Qm"/>
<map_variables variable_2="Qv" variable_1="Qv"/>
</connection>
<connection>
<map_components component_2="v" component_1="model_parameters"/>
<map_variables variable_2="Qmax" variable_1="Qmax"/>
<map_variables variable_2="theta" variable_1="theta"/>
<map_variables variable_2="sigma" variable_1="sigma"/>
</connection>
<connection>
<map_components component_2="m" component_1="model_parameters"/>
<map_variables variable_2="Qmax" variable_1="Qmax"/>
<map_variables variable_2="theta" variable_1="theta"/>
<map_variables variable_2="sigma" variable_1="sigma"/>
</connection>
<connection>
<map_components component_2="H" component_1="m"/>
<map_variables variable_2="Qm" variable_1="Qm"/>
</connection>
<connection>
<map_components component_2="D" component_1="H"/>
<map_variables variable_2="H" variable_1="H"/>
</connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
<foo xmlns="http://www.example.org/foo/bar/baz" />
<!--
The following RDF block contains metadata that applies to this document
as a whole, as indicated by the empty about attribute on the
<rdf:Description> element.
-->
<rdf:Description rdf:about="">
<!--
The Model Builder Metadata. The Dublin Core "creator" element is used
to indicate the person who translated the model into CellML.
-->
<dc:creator rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Lloyd</vCard:Family>
<vCard:Given>Catherine</vCard:Given>
</vCard:N>
<vCard:EMAIL rdf:parseType="Resource">
<rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
</vCard:EMAIL>
<vCard:ORG rdf:parseType="Resource">
<vCard:Orgname>Auckland Bioengineering Institute</vCard:Orgname>
</vCard:ORG>
</dc:creator>
<!--
The Creation Date metadata. This is the date on which the model
was translated into CellML.
-->
<dcterms:created rdf:parseType="Resource">
<dcterms:W3CDTF>2009-10-29</dcterms:W3CDTF>
</dcterms:created>
<!--
The Modification History metadata. This lists the changes that have been
made to the document, who made the changes, and when they were made.
-->
</rdf:Description>
<!--
The following metadata refers to the model itself, as indicated by the
reference to the ID "phillips_2008", which is
declared on the <model> element.
-->
<rdf:Description rdf:about="#phillips_2008">
<!-- A human readable name for the model. -->
<dc:title>
The Phillips and Robinson 2008 quantitative model of sleep deprivation
</dc:title>
<!-- A comment regarding the model. -->
<cmeta:comment rdf:parseType="Resource">
<rdf:value>
This is the CellML description of the Phillips and Robinson 2008 quantitative model of sleep deprivation
</rdf:value>
<!-- The creator of the comment. -->
<dc:creator rdf:parseType="Resource">
<vCard:FN>Catherine Lloyd</vCard:FN>
</dc:creator>
</cmeta:comment>
<!-- Keyword(s) -->
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>sleep</rdf:li>
<rdf:li>neurobiology</rdf:li>
<rdf:li>diurnal cycle</rdf:li>
<rdf:li>circadian rhythms</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<bqs:reference rdf:parseType="Resource">
<bqs:Pubmed_id>18805427</bqs:Pubmed_id>
<bqs:JournalArticle rdf:parseType="Resource">
<dc:creator>
<rdf:Seq>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Phillips</vCard:Family>
<vCard:Given>A</vCard:Given>
<vCard:Other>J</vCard:Other>
<vCard:Other>K</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Robinson</vCard:Family>
<vCard:Given>P</vCard:Given>
<vCard:Other>A</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
</rdf:Seq>
</dc:creator>
<dc:title>
Sleep deprivation in a quantitative physiologically based model of the ascending arousal system
</dc:title>
<dcterms:issued rdf:parseType="Resource">
<dcterms:W3CDTF>2008-12-21</dcterms:W3CDTF>
</dcterms:issued>
<bqs:Journal rdf:parseType="Resource">
<dc:title>Journal of Theoretical Biology</dc:title>
</bqs:Journal>
<bqs:volume>255</bqs:volume>
<bqs:first_page>413</bqs:first_page>
<bqs:last_page>423</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
</model>
