- Author:
- Catherine Lloyd <c.lloyd@auckland.ac.nz>
- Date:
- 2010-07-08 02:32:39+12:00
- Desc:
- Removed comments from metadata.
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/phillips_2007/rawfile/33e2a9fd5d8c7997657e869b1b04184636237049/phillips_2007.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_2007" name="phillips_2007">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>A quantitative model of sleep-wake dynamics based on the physiology of the brainstem ascending arousal system</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<sect1 id="sec_status">
<title>Model Status</title>
<para>
This CellML model runs in both OpenCell and COR to replicate figure 2 in the paper. The units have been checked and they are consistent.
</para>
</sect1>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
Abstract: A quantitative, physiology-based model of the ascending arousal system is developed, using continuum neuronal population modeling, which involves averaging properties such as firing rates across neurons in each population. The model includes the ventrolateral preoptic area (VLPO), where circadian and homeostatic drives enter the system, the monoaminergic and cholinergic nuclei of the ascending arousal system, and their interconnections. The human sleep-wake cycle is governed by the activities of these nuclei, which modulate the behavioral state of the brain via diffuse neuromodulatory projections. The model parameters are not free since they correspond to physiological observables. Approximate parameter bounds are obtained by requiring consistency with physiological and behavioral measures, and the model replicates the human sleep-wake cycle, with physiologically reasonable voltages and firing rates. Mutual inhibition between the wake-promoting monoaminergic group and sleep-promoting VLPO causes 'flip-flop' behavior, with most time spent in 2 stable steady states corresponding to wake and sleep, with transitions between them on a timescale of a few minutes. The model predicts hysteresis in the sleep-wake cycle, with a region of bistability of the wake and sleep states. Reducing the monoaminergic-VLPO mutual inhibition results in a smaller hysteresis loop. This makes the model more prone to wake-sleep transitions in both directions and makes the states less distinguishable, as in narcolepsy. The model behavior is robust across the constrained parameter ranges, but with sufficient flexibility to describe a wide range of observed phenomena.
</para>
<para>
The original paper reference is cited below:
</para>
<para>
A quantitative model of sleep-wake dynamics based on the physiology of the brainstem ascending arousal system, A.J.K. Phillips and P.A. Robinson, 2007,
<emphasis>Journal of Biological Rhythms</emphasis>, 22, 167-179. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/17440218">PubMed ID: 17440218</ulink>
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="phillips_2007.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram of the Phillips and Robinson 2007 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="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>
<units name="mV_per_nM">
<unit units="mV"/>
<unit units="nM" exponent="-1"/>
</units>
<units name="per_second">
<unit units="second" exponent="-1"/>
</units>
<component name="environment">
<variable units="hour" public_interface="out" name="time"/>
</component>
<component name="v" cmeta:id="v">
<rdf:RDF>
<rdf:Description rdf:about="#v">
<dc:title>v</dc:title>
<dcterms:alternative>VLPO</dcterms:alternative>
<dcterms:alternative>ventrolateral preoptic area</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="per_second" public_interface="out" name="Qv"/>
<variable units="mV" public_interface="out" name="Vv" initial_value="0.0"/>
<variable units="second" name="tau_v" initial_value="10.0"/>
<variable units="mV_second" name="v_vm" initial_value="-1.9"/>
<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>
<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>
</apply>
</math>
</component>
<component name="a" cmeta:id="a">
<rdf:RDF>
<rdf:Description rdf:about="#a">
<dc:title>a</dc:title>
<dcterms:alternative>ACh</dcterms:alternative>
<dcterms:alternative>acetylcholine group</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="per_second" name="Qa"/>
<variable units="mV" name="Va"/>
<variable units="mV" name="Vao" initial_value="1.0"/>
<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="3">
<apply>
<eq/>
<ci>Qa</ci>
<apply>
<divide/>
<ci>Qmax</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>Va</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/>
<ci>Va</ci>
<ci>Vao</ci>
</apply>
</math>
</component>
<component name="m" cmeta:id="m">
<rdf:RDF>
<rdf:Description rdf:about="#m">
<dc:title>m</dc:title>
<dcterms:alternative>MA</dcterms:alternative>
<dcterms:alternative>monoaminergic group</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<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.9"/>
<variable units="mV" name="v_maQao" initial_value="1.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="5">
<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="6">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Vm</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_maQao</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>
</apply>
</math>
</component>
<component name="H" cmeta:id="H">
<rdf:RDF>
<rdf:Description rdf:about="#H">
<dc:title>H</dc:title>
<dcterms:alternative>somnogen level</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nM" public_interface="out" name="H" initial_value="15.0"/>
<variable units="hour" name="chi" initial_value="10.8"/>
<variable units="nM_second" name="mu" initial_value="3.6"/>
<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="7">
<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" cmeta:id="D">
<rdf:RDF>
<rdf:Description rdf:about="#D">
<dc:title>D</dc:title>
<dcterms:alternative>total drive function</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="mV" public_interface="out" name="D"/>
<variable units="dimensionless" name="C"/>
<variable units="dimensionless" name="c0" initial_value="1.0"/>
<variable units="per_hour" name="omega"/>
<variable units="mV" name="v_vc" initial_value="-6.3"/>
<variable units="mV_per_nM" name="v_vh" initial_value="0.19"/>
<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="8">
<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="9">
<apply>
<eq/>
<ci>C</ci>
<apply>
<plus/>
<ci>c0</ci>
<apply>
<cos/>
<apply>
<times/>
<ci>omega</ci>
<ci>time</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="10">
<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="a"/>
<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="a" component_1="D"/>
<map_variables variable_2="D" variable_1="D"/>
</connection>
<connection>
<map_components component_2="a" 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#">
<rdf:Description rdf:about="">
<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>
<dcterms:created rdf:parseType="Resource">
<dcterms:W3CDTF>2009-11-09</dcterms:W3CDTF>
</dcterms:created>
</rdf:Description>
<rdf:Description rdf:about="#phillips_2007">
<dc:title>
The Phillips and Robinson 2007 quantitative model of sleep-wake dynamics
</dc:title>
<cmeta:comment rdf:parseType="Resource">
<rdf:value>
This is the CellML description of the Phillips and Robinson 2007 quantitative model of sleep-wake dynamics
</rdf:value>
<dc:creator rdf:parseType="Resource">
<vCard:FN>Catherine Lloyd</vCard:FN>
</dc:creator>
</cmeta:comment>
<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>17440218</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>
A quantitative model of sleep-wake dynamics based on the physiology of the brainstem ascending arousal system
</dc:title>
<dcterms:issued rdf:parseType="Resource">
<dcterms:W3CDTF>2007-04</dcterms:W3CDTF>
</dcterms:issued>
<bqs:Journal rdf:parseType="Resource">
<dc:title>Journal of Biological Rhythms</dc:title>
</bqs:Journal>
<bqs:volume>22</bqs:volume>
<bqs:first_page>167</bqs:first_page>
<bqs:last_page>179</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
</model>
