Location: Phillips 2007 @ 72ec675bef4b / phillips_2007.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2009-11-09 13:19 +1300
Desc:
Added files to the new workspace.
Permanent Source URI:
http://models.cellml.org/workspace/phillips_2007/rawfile/72ec675bef4b533604f518b0ce5fd7ce9ebdfab8/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>Bioengineering Institute, University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<sect1 id="sec_status">
				<title>Model Status</title>
				<para>
     This CellML model runs in both PCEnv 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/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=17440218&amp;query_hl=1&amp;itool=pubmed_docsum">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">
    <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">
    <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">
    <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">
    <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">
    <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#">
    <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-11-09</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_2007", which is    
      declared on the <model> element.
    -->
    <rdf:Description rdf:about="#phillips_2007">
      <!-- A human readable name for the model. -->
      <dc:title>
        The Phillips and Robinson 2007 quantitative model of sleep-wake dynamics
      </dc:title>
      
      <!-- A comment regarding the model. -->
      <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>
        <!-- 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>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>