<?xml version="1.0"?>
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.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 prefix="milli" units="volt"/>
</units>
<units name="mV_second">
<unit units="mV"/>
<unit units="second"/>
</units>
<units name="hour">
<unit multiplier="3600" units="second"/>
</units>
<units name="per_hour">
<unit exponent="-1" units="hour"/>
</units>
<units name="second_per_hour">
<unit units="second"/>
<unit exponent="-1" units="hour"/>
</units>
<units name="nM">
<unit prefix="nano" units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="nM_second">
<unit units="nM"/>
<unit units="second"/>
</units>
<units name="mV_per_nM">
<unit units="mV"/>
<unit exponent="-1" units="nM"/>
</units>
<units name="per_second">
<unit exponent="-1" units="second"/>
</units>
<component name="environment">
<variable cmeta:id="environment_time" name="time" public_interface="out" units="hour"/>
</component>
<component cmeta:id="v" name="v">
<variable name="Qv" public_interface="out" units="per_second"/>
<variable cmeta:id="v_Vv" initial_value="0.0" name="Vv" public_interface="out" units="mV"/>
<variable initial_value="10.0" name="tau_v" units="second"/>
<variable initial_value="-1.9" name="v_vm" units="mV_second"/>
<variable name="Qm" public_interface="in" units="per_second"/>
<variable name="Qmax" public_interface="in" units="per_second"/>
<variable name="D" public_interface="in" units="mV"/>
<variable name="theta" public_interface="in" units="mV"/>
<variable name="sigma" public_interface="in" units="mV"/>
<variable name="time" public_interface="in" units="hour"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="1">
<apply>
<eq/>
<ci>Qv</ci>
<apply>
<divide/>
<ci>Qmax</ci>
<apply>
<plus/>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second_per_hour">3600</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="a" name="a">
<variable name="Qa" units="per_second"/>
<variable cmeta:id="a_Va" name="Va" units="mV"/>
<variable initial_value="1.0" name="Vao" units="mV"/>
<variable name="Qmax" public_interface="in" units="per_second"/>
<variable name="D" public_interface="in" units="mV"/>
<variable name="theta" public_interface="in" units="mV"/>
<variable name="sigma" public_interface="in" units="mV"/>
<variable name="time" public_interface="in" units="hour"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="3">
<apply>
<eq/>
<ci>Qa</ci>
<apply>
<divide/>
<ci>Qmax</ci>
<apply>
<plus/>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 cmeta:id="m" name="m">
<variable name="Qm" public_interface="out" units="per_second"/>
<variable cmeta:id="m_Vm" initial_value="0.0" name="Vm" public_interface="out" units="mV"/>
<variable initial_value="10.0" name="tau_m" units="second"/>
<variable initial_value="-1.9" name="v_mv" units="mV_second"/>
<variable initial_value="1.0" name="v_maQao" units="mV"/>
<variable name="Qmax" public_interface="in" units="per_second"/>
<variable name="theta" public_interface="in" units="mV"/>
<variable name="sigma" public_interface="in" units="mV"/>
<variable name="Qv" public_interface="in" units="per_second"/>
<variable name="time" public_interface="in" units="hour"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="5">
<apply>
<eq/>
<ci>Qm</ci>
<apply>
<divide/>
<ci>Qmax</ci>
<apply>
<plus/>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second_per_hour">3600</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="H" name="H">
<variable initial_value="15.0" name="H" public_interface="out" units="nM"/>
<variable initial_value="10.8" name="chi" units="hour"/>
<variable initial_value="3.6" name="mu" units="nM_second"/>
<variable name="Qm" public_interface="in" units="per_second"/>
<variable name="time" public_interface="in" units="hour"/>
<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 cmeta:id="D" name="D">
<variable cmeta:id="D_D" name="D" public_interface="out" units="mV"/>
<variable name="C" units="dimensionless"/>
<variable initial_value="1.0" name="c0" units="dimensionless"/>
<variable name="omega" units="per_hour"/>
<variable initial_value="-6.3" name="v_vc" units="mV"/>
<variable initial_value="0.19" name="v_vh" units="mV_per_nM"/>
<variable name="H" public_interface="in" units="nM"/>
<variable name="time" public_interface="in" units="hour"/>
<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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
<pi/>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="hour">24</cn>
</apply>
</apply>
</math>
</component>
<component name="model_parameters">
<variable initial_value="100.0" name="Qmax" public_interface="out" units="per_second"/>
<variable initial_value="10.0" name="theta" public_interface="out" units="mV"/>
<variable initial_value="3.0" name="sigma" public_interface="out" units="mV"/>
</component>
<connection>
<map_components component_1="v" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="a" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="m" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="H" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="D" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="D" component_2="v"/>
<map_variables variable_1="D" variable_2="D"/>
</connection>
<connection>
<map_components component_1="m" component_2="v"/>
<map_variables variable_1="Qm" variable_2="Qm"/>
<map_variables variable_1="Qv" variable_2="Qv"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="v"/>
<map_variables variable_1="Qmax" variable_2="Qmax"/>
<map_variables variable_1="theta" variable_2="theta"/>
<map_variables variable_1="sigma" variable_2="sigma"/>
</connection>
<connection>
<map_components component_1="D" component_2="a"/>
<map_variables variable_1="D" variable_2="D"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="a"/>
<map_variables variable_1="Qmax" variable_2="Qmax"/>
<map_variables variable_1="theta" variable_2="theta"/>
<map_variables variable_1="sigma" variable_2="sigma"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="m"/>
<map_variables variable_1="Qmax" variable_2="Qmax"/>
<map_variables variable_1="theta" variable_2="theta"/>
<map_variables variable_1="sigma" variable_2="sigma"/>
</connection>
<connection>
<map_components component_1="m" component_2="H"/>
<map_variables variable_1="Qm" variable_2="Qm"/>
</connection>
<connection>
<map_components component_1="H" component_2="D"/>
<map_variables variable_1="H" variable_2="H"/>
</connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"><rdf:Description rdf:nodeID="n1"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n2"/></N></rdf:Description><rdf:Description rdf:about="#D"><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">total drive function</alternative><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">D</title></rdf:Description><rdf:Description rdf:nodeID="n3"><Orgname xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Auckland Bioengineering Institute</Orgname></rdf:Description><rdf:Description rdf:about="#a"><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">acetylcholine group</alternative><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">ACh</alternative><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">a</title></rdf:Description><rdf:Description rdf:about="#v"><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">VLPO</alternative><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">ventrolateral preoptic area</alternative><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">v</title></rdf:Description><rdf:Description rdf:nodeID="n4"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n5"/></N></rdf:Description><rdf:Description rdf:nodeID="n6"><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Bag"/></rdf:type><rdf:_3 xml:lang="en">diurnal cycle</rdf:_3><rdf:_1 xml:lang="en">sleep</rdf:_1><rdf:_2 xml:lang="en">neurobiology</rdf:_2><rdf:_4 xml:lang="en">circadian rhythms</rdf:_4></rdf:Description><rdf:Description rdf:nodeID="n7"><endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">48</endingValue></rdf:Description><rdf:Description rdf:nodeID="n8"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2009-11-09</W3CDTF></rdf:Description><rdf:Description rdf:nodeID="n2"><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">K</Other><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">J</Other><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">A</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Phillips</Family></rdf:Description><rdf:Description rdf:nodeID="n9"><boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n10"/></boundIntervals></rdf:Description><rdf:Description rdf:nodeID="n11"><subject_type xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">keyword</subject_type><rdf:value><rdf:Description rdf:nodeID="n6"/></rdf:value></rdf:Description><rdf:Description rdf:nodeID="n12"><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n13"/></creator><rdf:value xml:lang="en">
This is the CellML description of the Phillips and Robinson 2007 quantitative model of sleep-wake dynamics
</rdf:value></rdf:Description><rdf:Description rdf:nodeID="n14"><ORG xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n3"/></ORG><EMAIL xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n15"/></EMAIL><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n16"/></N></rdf:Description><rdf:Description rdf:nodeID="n15"><rdf:type><rdf:Description rdf:about="http://imc.org/vCard/3.0#internet"/></rdf:type><rdf:value xml:lang="en">c.lloyd@auckland.ac.nz</rdf:value></rdf:Description><rdf:Description rdf:nodeID="n5"><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">A</Other><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">P</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Robinson</Family></rdf:Description><rdf:Description rdf:nodeID="n17"><JournalArticle xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n18"/></JournalArticle><Pubmed_id xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">17440218</Pubmed_id></rdf:Description><rdf:Description rdf:nodeID="n19"><subject xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n11"/></subject></rdf:Description><rdf:Description rdf:nodeID="n13"><FN xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Catherine Lloyd</FN></rdf:Description><rdf:Description rdf:nodeID="n16"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Catherine</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Lloyd</Family></rdf:Description><rdf:Description rdf:about="#phillips_2007"><cmeta:comment><rdf:Description rdf:nodeID="n12"/></cmeta:comment><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
The Phillips and Robinson 2007 quantitative model of sleep-wake dynamics
</title><simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n9"/></simulation><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n19"/></reference><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n17"/></reference></rdf:Description><rdf:Description rdf:nodeID="n20"><Person xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n4"/></Person></rdf:Description><rdf:Description rdf:nodeID="n21"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Journal of Biological Rhythms</title></rdf:Description><rdf:Description rdf:about="#H"><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">somnogen level</alternative><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">H</title></rdf:Description><rdf:Description rdf:nodeID="n10"><rdf:first><rdf:Description rdf:nodeID="n7"/></rdf:first><rdf:rest><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/></rdf:rest></rdf:Description><rdf:Description rdf:nodeID="n18"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
A quantitative model of sleep-wake dynamics based on the physiology of the brainstem ascending arousal system
</title><issued xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:nodeID="n22"/></issued><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n23"/></creator><Journal xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n21"/></Journal><first_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">167</first_page><last_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">179</last_page><volume xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">22</volume></rdf:Description><rdf:Description rdf:nodeID="n24"><Person xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n1"/></Person></rdf:Description><rdf:Description rdf:about=""><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n14"/></creator><created xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:nodeID="n8"/></created></rdf:Description><rdf:Description rdf:about="#m"><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">monoaminergic group</alternative><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">MA</alternative><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">m</title></rdf:Description><rdf:Description rdf:nodeID="n23"><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Seq"/></rdf:type><rdf:_1><rdf:Description rdf:nodeID="n24"/></rdf:_1><rdf:_2><rdf:Description rdf:nodeID="n20"/></rdf:_2></rdf:Description><rdf:Description rdf:nodeID="n22"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2007-04</W3CDTF></rdf:Description></rdf:RDF></model>
