Rendering of the source text

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  st_hilaire_model_2007.xml

CREATED :  24th September 2007

LAST MODIFIED : 24th September 2007

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.1 Specification.

DESCRIPTION :  This file contains a CellML description of St Hilaire et al.'s 2007 physiologically based mathematical model of melatonin including ocular light suppression and interactions with the circadian pacemaker.

CHANGES:  
  
--><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="st_hilaire_model_2007" name="st_hilaire_model_2007_version01">

<documentation xmlns="http://cellml.org/tmp-documentation">

<article>
  <articleinfo>
  <title>A physiologically based mathematical model of melatonin including ocular light suppression and interactions with the circadian pacemaker</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
        This CellML version of the model has been checked in COR and PCEnv.  Currently it runs in PCEnv but it does not recreate the published results.  Once the model author loads the original model onto their website (<ulink url="http://sleep.med.harvard.edu/people/faculty/225/klerman">found here</ulink>) we'll check the CellML code against this original version and will hopefully be able to fix the CellML model.
          </para>
  </section>
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
Since it is impossible to directly measure the activity of the circadian pacemaker in humans, alternative markers of the circadian rhythms must be used, such as plasma melatonin, core body temperature and cortisol levels.  Of these markers, it has been demonstrated that plasma melatonin is the most reliable, and its presence in both blood and saliva makes it more accessible for assay.  Furthermore, unlike cortisol melatonin is relatively unaffected by sleep-wake state.
</para>

<para>
Several mathematical models of circadian rhythms have been developed based on plasma melatonin data.  These employ a number of different methods to estimate circadian phase and amplitude from the measured melatonin data.  However, the accuracy of these methods comes into question since melatonin levels can be subject to inter-individual variation and can also be influenced by environmental factors - especially light.  Furthermore, the underlying physiology of melatonin secretion and clearance is often ignored.
</para>

<para>
In the study described here, Melissa St Hilaire <emphasis>et al.</emphasis> address these issues by developing a physiologically based mathematical model of melatonin, including ocular light suppression and interactions with the circadian pacemaker.  Their model is an extension of the previously published model of melatonin rhythms by Brown <emphasis>et al.</emphasis> (also described in CellML).  A saliva component has been added, and the effect of light on melatonin levels has been incorporated into the model by including equations based on those published by Kronauer <emphasis>et al.</emphasis>. 
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="sthilaire_2007.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the components and variables in a mathematical model of human melatonin.  Melatonin is synthesised in the pineal gland where its precursor tryptophan is converted to serotonin, and then to melatonin.  The rate determining step in this synthesis pathway is the governed by the activity of N-acetylserotonin (NAT) (here represented by A), and the activity of this enzyme is controlled by the suprachiasmatic nuclei of the hypothalamus, the main mammalian circadian pacemaker.  Melatonin is not stored in the pineal gland.  Once synthesised melatonin diffuses into the systemic circulatory system.  From the plasma it is absorbed by various tissues and partially bind with albumin to be subsequently cleared by the liver.  In addition, some melatonin infuses from the plasma into the saliva, from which it can also be cleared.</caption>
</informalfigure>

<para>
The complete original paper reference is cited below:
</para>

<para>
<ulink url="http://www.blackwell-synergy.com/doi/abs/10.1111/j.1600-079X.2007.00477.x">A physiologically based mathematical model of melatonin including ocular light suppression and interactions with the circadian pacemaker</ulink>, Melissa A. St Hilaire, Claude Gronfier, Jamie M. Zeitzer, and Elizabeth B. Klerman, 2007, <ulink url="http://www.blackwell-synergy.com/toc/jpi/43/3">
            <emphasis>Journal of Pineal Research</emphasis>
          </ulink>, 43, 294-304.  (<ulink url="http://www.blackwell-synergy.com/doi/full/10.1111/j.1600-079X.2007.00477.x">Full text</ulink> and PDF versions of the article are available to journal subscribers on the <emphasis>Journal of Pineal Research</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=17803528&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 17803528</ulink>
</para>


</sect1>
</article>
</documentation> 



  
  <units name="minute">
    <unit units="second" multiplier="60.0"/>
  </units>  

  <units name="hour">
    <unit units="second" multiplier="3600.0"/>
  </units>

  <units name="picomolar">
    <unit units="mole" prefix="pico"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="flux">
    <unit units="picomolar"/>
    <unit units="minute" exponent="-1"/>
  </units>

  <units name="first_order_rate_constant">
    <unit units="minute" exponent="-1"/>
  </units>
  
  
  <component name="environment">
    <variable units="hour" public_interface="out" name="time"/>
  </component>
 
  
  <component cmeta:id="H1" name="H1">
     
    <variable units="picomolar" public_interface="out" name="H1" initial_value="0.5"/>
    
    <variable units="flux" public_interface="in" name="A"/>
    <variable units="first_order_rate_constant" public_interface="in" name="betaIP"/>
    <variable units="hour" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> H1 </ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> betaIP </ci>
              <ci> H1 </ci>
            </apply>
          </apply>
          <ci> A </ci>
        </apply>
      </apply>  
    </math> 
  </component>

  <component cmeta:id="H2" name="H2">
     
    <variable units="picomolar" public_interface="out" name="H2" initial_value="0.5"/>
    
    <variable units="picomolar" public_interface="in" name="H1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="betaIP"/>
    <variable units="first_order_rate_constant" public_interface="in" name="betaIS"/>
    <variable units="first_order_rate_constant" public_interface="in" name="betaCP"/>
    <variable units="hour" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> H2 </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> betaIP </ci>
            <ci> H1 </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> betaCP </ci>
              <ci> H2 </ci>
            </apply>
            <apply>
              <times/>
              <ci> betaIS </ci>
              <ci> H2 </ci>
            </apply>
          </apply>
        </apply>
      </apply>  
    </math> 
  </component>

  <component cmeta:id="H3" name="H3">
     
    <variable units="picomolar" public_interface="out" name="H3" initial_value="0.5"/>
    
    <variable units="picomolar" public_interface="in" name="H2"/>
    <variable units="first_order_rate_constant" public_interface="in" name="betaIS"/>
    <variable units="first_order_rate_constant" public_interface="in" name="betaCS"/>
    <variable units="hour" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> H3 </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> betaIS </ci>
            <ci> H2 </ci>
          </apply>
          <apply>
            <times/>
            <ci> betaCS </ci>
            <ci> H3 </ci>
          </apply>
        </apply>
      </apply>  
    </math> 
  </component>

  <component cmeta:id="A" name="A">
     
    <variable units="flux" public_interface="out" name="A"/>
    <variable units="flux" public_interface="out" name="A_"/>
    
    <variable units="flux" name="Amax" initial_value="4.834"/>
    <variable units="dimensionless" name="m" initial_value="7.0"/>
    <variable units="hour" name="t_on" initial_value="21.93"/>
    <variable units="hour" name="t_off" initial_value="29.867"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="alpha"/>
    <variable units="first_order_rate_constant" public_interface="in" name="lamda"/>
    <variable units="dimensionless" public_interface="in" name="B_"/>
    <variable units="hour" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> A </ci> 
        </apply>
        <piecewise>
          <piece>
            <apply>
              <times/>
              <ci> Amax </ci> 
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <exp/>
                    <apply>
                      <times/>
                      <apply>
                        <minus/>
                        <ci> lamda </ci>
                      </apply>
                      <apply>
                        <minus/>
                        <ci> time </ci>
                        <ci> t_on </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <exp/>
                    <apply>
                      <times/>
                      <apply>
                        <minus/>
                        <ci> lamda </ci>
                      </apply>
                      <apply>
                        <minus/>
                        <ci> t_off </ci>
                        <ci> t_on </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <or/>
              <apply>
                <lt/>
                <ci> time </ci>
                <ci> t_off </ci>
              </apply>
              <apply>
                <geq/>
                <ci> time </ci>
                <ci> t_on </ci>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <times/>
              <ci> Amax </ci> 
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <ci> alpha </ci>
                  </apply>
                  <apply>
                    <minus/>
                    <ci> time </ci>
                    <ci> t_off </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <geq/>
              <ci> time </ci>
              <ci> t_off </ci>
            </apply>
          </piece>
        </piecewise>
      </apply>

      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> A_ </ci> 
        </apply>
        <apply>
          <times/>
          <ci> A </ci> 
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <ci> m </ci>
              <ci> B_ </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component>

  <component cmeta:id="Process_L" name="Process_L">
     
    <variable units="dimensionless" public_interface="out" name="B_"/>

    <variable units="dimensionless" name="n" initial_value="7.0"/>
    <variable units="dimensionless" name="p" initial_value="0.6"/>
    <variable units="dimensionless" name="I"/>
    <variable units="dimensionless" name="I0" initial_value="9500.0"/>
    <variable units="dimensionless" name="G" initial_value="37.0"/>
    <variable units="first_order_rate_constant" name="alph_0" initial_value="0.1"/>
    <variable units="first_order_rate_constant" name="beta" initial_value="0.007"/>
    <variable units="first_order_rate_constant" name="alph"/>

    <variable units="hour" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> alph </ci> 
        <apply>
          <times/>
          <ci> alph_0 </ci>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> I </ci>
              <ci> I0 </ci>
            </apply>
            <ci> p </ci>
          </apply>
          <apply>
            <divide/>
            <ci> I </ci>
            <apply>
              <plus/>
              <ci> I </ci>
              <cn cellml:units="dimensionless"> 100.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply> 

      <apply>
        <eq/> 
        <ci> B_ </ci> 
        <apply>
          <times/>
          <ci> G </ci>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> n </ci>
          </apply>
          <ci> alph </ci>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> n </ci> 
        </apply>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 60.0 </cn>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> alph </ci>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> n </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> beta </ci>
              <ci> n </ci>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> I </ci> 
        </apply>
        <piecewise>
          <piece>
            <cn cellml:units="dimensionless"> 0.0 </cn>
            <apply>
              <and/>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="hour"> 0.0 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="hour"> 4.75 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="dimensionless"> 9500.0 </cn>
            <apply>
              <and/>
              <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="hour"> 4.75 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="hour"> 11.2 </cn>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="dimensionless"> 0.0 </cn>
          </otherwise>
        </piecewise>
      </apply>    
    </math> 
  </component>

  <component cmeta:id="Process_P" name="Process_P">
     
    
    <variable units="dimensionless" name="x"/>
    <variable units="dimensionless" name="xc"/>
    <variable units="dimensionless" name="B"/>
    <variable units="dimensionless" name="k" initial_value="0.55"/>
    <variable units="dimensionless" name="mu" initial_value="0.13"/>
    <variable units="dimensionless" name="q" initial_value="0.333333"/>
    <variable units="dimensionless" name="tau_x" initial_value="24.2"/>
    
    <variable units="dimensionless" public_interface="in" name="B_"/>
    <variable units="hour" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> x </ci> 
        </apply>
        <apply>
          <times/>
          <apply>
            <divide/>
            <pi/>
            <cn cellml:units="dimensionless"> 12.0 </cn>
          </apply>
          <apply>
            <plus/>
            <ci> xc </ci>
            <apply>
              <times/>
              <ci> mu </ci>
              <apply>
                <minus/>
                <apply>
                  <plus/>
                  <apply>
                    <times/>
                    <apply>
                      <divide/>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                      <cn cellml:units="dimensionless"> 3.0 </cn>  
                    </apply>
                    <ci> x </ci> 
                  </apply>
                  <apply>
                    <times/>
                    <apply>
                      <divide/>
                      <cn cellml:units="dimensionless"> 4.0 </cn>
                      <cn cellml:units="dimensionless"> 3.0 </cn>  
                    </apply>
                    <apply>
                      <power/>
                      <ci> x </ci>
                      <cn cellml:units="dimensionless"> 3.0 </cn>  
                    </apply> 
                  </apply>
                </apply>
                <apply>
                  <times/>
                  <apply>
                    <divide/>
                    <cn cellml:units="dimensionless"> 256.0 </cn>
                    <cn cellml:units="dimensionless"> 105.0 </cn>  
                  </apply>
                  <apply>
                    <power/>
                    <ci> x </ci>
                    <cn cellml:units="dimensionless"> 7.0 </cn>  
                  </apply> 
                </apply>
              </apply>
            </apply>
            <ci> B </ci> 
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> xc </ci> 
        </apply>
        <apply>
          <times/>
          <apply>
            <divide/>
            <pi/>
            <cn cellml:units="dimensionless"> 12.0 </cn>
          </apply>
          <apply>
            <plus/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> q </ci>
                <ci> B </ci>
                <ci> xc </ci>
              </apply>
              <apply>
                <times/>
                <ci> x </ci>
                <apply>
                  <power/>
                  <apply>
                    <divide/>
                    <cn cellml:units="dimensionless"> 24.0 </cn>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 0.99729 </cn>  
                      <ci> tau_x </ci> 
                    </apply>
                  </apply>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> k </ci>
              <ci> B </ci>
            </apply>
          </apply>
        </apply>
      </apply>    

      <apply>
        <eq/> 
        <ci> B </ci> 
        <apply>
          <times/>
          <ci> B_ </ci>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.4 </cn>
              <ci> x </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.4 </cn>
              <ci> xc </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math> 
  </component>
  
  <component name="model_parameters">
    <variable units="first_order_rate_constant" public_interface="out" name="betaIP" initial_value="0.047"/>
    <variable units="first_order_rate_constant" public_interface="out" name="betaCP" initial_value="0.025"/>
    <variable units="first_order_rate_constant" public_interface="out" name="betaIS" initial_value="0.065"/>
    <variable units="first_order_rate_constant" public_interface="out" name="betaCS" initial_value="0.028"/>
    <variable units="first_order_rate_constant" public_interface="out" name="alpha" initial_value="0.067"/>
    <variable units="first_order_rate_constant" public_interface="out" name="lamda" initial_value="0.588"/>
  </component>


  <connection>
    <map_components component_2="environment" component_1="H1"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="H2"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="H3"/>
    <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="Process_L"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Process_P"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="A" component_1="H1"/>
    <map_variables variable_2="A" variable_1="A"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="H1"/>
    <map_variables variable_2="betaIP" variable_1="betaIP"/>
  </connection>

  <connection>
    <map_components component_2="H1" component_1="H2"/>
    <map_variables variable_2="H1" variable_1="H1"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="H2"/>
    <map_variables variable_2="betaIP" variable_1="betaIP"/>
    <map_variables variable_2="betaIS" variable_1="betaIS"/>
    <map_variables variable_2="betaCP" variable_1="betaCP"/>
  </connection>

  <connection>
    <map_components component_2="H2" component_1="H3"/>
    <map_variables variable_2="H2" variable_1="H2"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="H3"/>
    <map_variables variable_2="betaIS" variable_1="betaIS"/>
    <map_variables variable_2="betaCS" variable_1="betaCS"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="A"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
    <map_variables variable_2="lamda" variable_1="lamda"/>
  </connection>

  <connection>
    <map_components component_2="Process_L" component_1="A"/>
    <map_variables variable_2="B_" variable_1="B_"/>
  </connection>

  <connection>
    <map_components component_2="Process_L" component_1="Process_P"/>
    <map_variables variable_2="B_" variable_1="B_"/>
  </connection>

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#7e9c2aff-aa14-4227-8c2e-d4365841c3da">
    <rdf:li>melatonin</rdf:li>
    <rdf:li>signal transduction</rdf:li>
    <rdf:li>pineal gland</rdf:li>
    <rdf:li>circadian</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#5e944931-8600-4e00-94b6-642fc2624d82">
    <rdf:li rdf:resource="rdf:#b7f1354d-469d-4e94-8ec8-5bcabcad2668"/>
    <rdf:li rdf:resource="rdf:#8af59fd6-0024-47c6-bc27-50949c4fac28"/>
    <rdf:li rdf:resource="rdf:#1641c416-3ee1-49da-a7d5-d7d81c74b7b2"/>
    <rdf:li rdf:resource="rdf:#a462fede-8368-4da9-9c8a-17418ec61aee"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#4f1994f1-fc8e-454d-959f-09cb11ab846d">
    <vCard:Given>Jamie</vCard:Given>
    <vCard:Family>Zeitzer</vCard:Family>
    <vCard:Other>M</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#ca1e7ed8-358e-4d6b-ae2e-4d67cb736954"/>
    <dcterms:created rdf:resource="rdf:#9336916b-9277-41a4-a64e-ab02485f634d"/>
    <dc:creator rdf:resource="rdf:#08c348f7-428a-4b55-922e-cc9ed3e5f4e0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b7f1354d-469d-4e94-8ec8-5bcabcad2668">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#17ab9aa7-3e59-4bb1-8cbf-cd521438370c"/>
  </rdf:Description>
  <rdf:Description rdf:about="/A">
    <dcterms:alternative>NAT activity</dcterms:alternative>
    <dc:title>A</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b6093df1-3c7a-4790-b0ed-90bf19fbba81">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#17ab9aa7-3e59-4bb1-8cbf-cd521438370c">
    <vCard:Given>Melissa</vCard:Given>
    <vCard:Family>St Hilaire</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#38dc65b6-c2b2-407c-9346-08ae99c4188f">
    <dc:creator rdf:resource="rdf:#5e944931-8600-4e00-94b6-642fc2624d82"/>
    <dc:title>A physiologically based mathematical model of melatonin including ocular light suppression and interactions with the circadian pacemaker</dc:title>
    <bqs:volume>43</bqs:volume>
    <bqs:first_page>294</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#ff0d6155-dd86-47b3-adda-381d4c9b1c02"/>
    <dcterms:issued rdf:resource="rdf:#2820d4e0-a024-4092-a56a-5aa93c4cc574"/>
    <bqs:last_page>304</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9336916b-9277-41a4-a64e-ab02485f634d">
    <dcterms:W3CDTF>2007-09-25T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#32c57989-fc37-4009-be68-c25cb4b1e2de">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8cb2686c-7390-44b0-b1e6-517a8b5c42b2">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2820d4e0-a024-4092-a56a-5aa93c4cc574">
    <dcterms:W3CDTF>2007-10</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d35c4dd3-a020-4059-b73b-6d42121b8464">
    <vCard:Given>Elizabeth</vCard:Given>
    <vCard:Family>Klerman</vCard:Family>
    <vCard:Other>B</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="/H3">
    <dcterms:alternative>melatonin concentration in the saliva</dcterms:alternative>
    <dc:title>H3</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="/H2">
    <dcterms:alternative>melatonin concentration in the plasma</dcterms:alternative>
    <dc:title>H2</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="/H1">
    <dcterms:alternative>melatonin concentration in the pineal gland</dcterms:alternative>
    <dc:title>H1</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1d4921bd-39f3-48ae-9e18-7a0c47d79598">
    <bqs:Pubmed_id>17803528</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#38dc65b6-c2b2-407c-9346-08ae99c4188f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1760e46d-08e2-47d8-af45-2e2d5f450d68">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#7e9c2aff-aa14-4227-8c2e-d4365841c3da"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ca1e7ed8-358e-4d6b-ae2e-4d67cb736954">
    <dc:creator rdf:resource="rdf:#f4192fc2-1e2e-4c28-8f64-d7ee3e3856bc"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f4192fc2-1e2e-4c28-8f64-d7ee3e3856bc">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#456de80b-a59c-4d17-9b8c-d3189ecee684">
    <dc:subject rdf:resource="rdf:#1760e46d-08e2-47d8-af45-2e2d5f450d68"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d3d54c5e-c57d-4b0b-aac5-e3f57bd4e37f">
    <vCard:Given>Claude</vCard:Given>
    <vCard:Family>Gronfier</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="#st_hilaire_model_2007">
    <dc:title>
        St Hilaire et al.'s 2007 physiologically based mathematical model of melatonin including ocular light supression and interactions with the circadian pacemaker.
      </dc:title>
    <cmeta:bio_entity>pineal gland</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#60e7da91-7d16-4989-8e16-1e70316ff4e6"/>
    <bqs:reference rdf:resource="rdf:#1d4921bd-39f3-48ae-9e18-7a0c47d79598"/>
    <bqs:reference rdf:resource="rdf:#456de80b-a59c-4d17-9b8c-d3189ecee684"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#60e7da91-7d16-4989-8e16-1e70316ff4e6">
    <dc:creator rdf:resource="rdf:#b6093df1-3c7a-4790-b0ed-90bf19fbba81"/>
    <rdf:value>This is the CellML description of St Hilaire et al.'s 2007 physiologically based mathematical model of melatonin including ocular light supression and interactions with the circadian pacemaker.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="/Process_P">
    <dcterms:alternative>dynamic circadian pacemaker</dcterms:alternative>
    <dc:title>Process_P</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="/Process_L">
    <dcterms:alternative>the physiological process by which light initiates a chemical reaction within the photo-pigments of the retinal photoreceptors</dcterms:alternative>
    <dc:title>Process_L</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1641c416-3ee1-49da-a7d5-d7d81c74b7b2">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#4f1994f1-fc8e-454d-959f-09cb11ab846d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8af59fd6-0024-47c6-bc27-50949c4fac28">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#d3d54c5e-c57d-4b0b-aac5-e3f57bd4e37f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a462fede-8368-4da9-9c8a-17418ec61aee">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#d35c4dd3-a020-4059-b73b-6d42121b8464"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ff0d6155-dd86-47b3-adda-381d4c9b1c02">
    <dc:title>Journal of Pineal Research</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#08c348f7-428a-4b55-922e-cc9ed3e5f4e0">
    <vCard:ORG rdf:resource="rdf:#387a4217-18cd-4ebc-a505-d31f1b418f2a"/>
    <vCard:EMAIL rdf:resource="rdf:#32c57989-fc37-4009-be68-c25cb4b1e2de"/>
    <vCard:N rdf:resource="rdf:#8cb2686c-7390-44b0-b1e6-517a8b5c42b2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#387a4217-18cd-4ebc-a505-d31f1b418f2a">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
</rdf:RDF>
</model>