Location: Ueda, Hagiwara, Kitano, 2001 @ 27118d5f2f43 / ueda_hagiwara_kitano_2001.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2009-06-17 16:16:04+12:00
Desc:
committing version03 of ueda_hagiwara_kitano_2001
Permanent Source URI:
https://staging.physiomeproject.org/workspace/ueda_hagiwara_kitano_2001/rawfile/27118d5f2f432dbda8a9c1d59aa4e8d92f184542/ueda_hagiwara_kitano_2001.cellml

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

CREATED :  23rd September 2003

LAST MODIFIED : 23rd September 2003

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/1/02 CellML Metadata 1.0 Specification.

DESCRIPTION :  This file contains a CellML description of Ueda et al.'s mathematical model of circadian rhythms in Drosophila.

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#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" name="ueda_hagiwara_kitano_2001_version01" cmeta:id="ueda_hagiwara_kitano_2001_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Interlocked Feedback Model of Drosophila Circadian Rhythm</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 model is known to run in both PCEnv and COR to replicate the published results.  The units have been checked and they are consistent.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
Many living organisms, from bacteria to plants, insects to mammals, display circadian rhythms.  These are spontaneously sustained oscillations with a period close to 24 hours.  Even in the absence of environmental cues, such as the light changes associated with day and night, organisms have been shown to retain their circadian behaviour, therefore suggesting that the rhythms are endogenous.  Experiments with <emphasis>Drosophila</emphasis> (fruit fly), <emphasis>Neurospora</emphasis> (fungus), cyanobacteria, plants and mammals have improved our understanding of the molecular mechanisms underlying circadian rhythms.  It seems that they rely on a negative feedback on gene expression.  A number of genes involved in circadian rhythms have been identified.  These include: <emphasis>Per</emphasis> (period), <emphasis>Tim</emphasis> (timeless), <emphasis>Cyc</emphasis> (cycle), <emphasis>Clk</emphasis> (clock), <emphasis>dbt</emphasis> (doubletime), and <emphasis>Cry</emphasis> (cryptochrome).  Among these genes, <emphasis>per</emphasis>, <emphasis>tim</emphasis> and <emphasis>Clk</emphasis> are rhythmically expressed.  <emphasis>per</emphasis> and <emphasis>tim</emphasis> mRNA levels peak in the evening, where as <emphasis>Clk</emphasis> mRNA levels peak late at night or early in the morning. 
</para>

<para>
Analysis of <emphasis>per</emphasis> and <emphasis>tim</emphasis> oscillatory expression has revealed a negative feedback loop (<emphasis>
            <emphasis>per-tim</emphasis> feedback loop</emphasis> in the figure below).  Several mathematical models of this negative feedback loop have been developed, and model simulations have shown that a delayed negative feedback loop is sufficient to produce spontaneous oscillations.  However, these models fail to take into account the <emphasis>Clk</emphasis> negative feedback loop.  Based on experimental analyses of <emphasis>Clk</emphasis> expression, Ueda <emphasis>et al.</emphasis> have developed a different type of mathematical model which they call an <emphasis>interlocked feedback model</emphasis> (as opposed to a <emphasis>simple model</emphasis>).  The interlocked feedback model is composed of two coupled negative feedback loops: a per-tim feedback loop which is activated by CLK-CYC and inhibited by PER-TIM, and a Clk feedback loop, which is inhibited by CLK-CYC and reactivated by PER-TIM. 
</para>

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

<para>
<ulink url="http://linkinghub.elsevier.com/retrieve/pii/S0022519300922264">Robust Oscillations within the Interlocked Feedback Model of <emphasis>Drosophila</emphasis> Circadian Rhythm</ulink>, Hiroki R. Ueda, Masatoshi Hagiwara, and Hiroaki Kitano, 2001, <ulink url="http://www.sciencedirect.com/science/journal/00225193">
            <emphasis>Journal of Theoretical Biology</emphasis>
          </ulink>, 210, 401-406.  (<ulink url="http://www.sciencedirect.com/science?_ob=MImg&amp;_imagekey=B6WMD-457VD36-24-1&amp;_cdi=6932&amp;_user=140507&amp;_orig=search&amp;_coverDate=06%2F21%2F2001&amp;_sk=997899995&amp;view=c&amp;wchp=dGLzVlz-zSkWz&amp;md5=6924ecbd484d35a98e1e70ee47751f1f&amp;ie=/sdarticle.pdf">A PDF</ulink> version of the article is available on the <emphasis>Journal of Theoretical Biology</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/11403560">PubMed ID: 11403460</ulink>
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>reaction diagram</title>
    </objectinfo>
    <imagedata fileref="ueda_2001.png"/>
  </imageobject>
</mediaobject>
<caption>The interlocked feedback model of <emphasis>Drosophila</emphasis> circadian rhythm.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
  
  
  <units name="nanomolar">
    <unit units="mole" prefix="nano"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="hour">
    <unit units="second" multiplier="3600.0"/>
  </units>
  
  <units name="flux">
    <unit units="nanomolar"/>
    <unit units="hour" exponent="-1"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="hour" exponent="-1"/>
  </units>
  
  <units name="second_order_rate_constant">
    <unit units="nanomolar" exponent="-1"/>
    <unit units="hour" exponent="-1"/>
  </units>
  
  
  <component name="environment">
    <variable units="hour" public_interface="out" name="time" cmeta:id="environment_time"/>
  </component>
  
  
  <component name="Per_m" cmeta:id="Per_m">
    <rdf:RDF>
      <rdf:Description rdf:about="Per_m">
        <dc:title>Per_m</dc:title>
        <dcterms:alternative>Per mRNA</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="Per_m_Per_m" name="Per_m" initial_value="0.5"/>
    
    <variable units="dimensionless" name="B1" initial_value="0.0"/>
    <variable units="flux" name="C1" initial_value="0.0"/> 
    <variable units="flux" name="S1" initial_value="1.45"/>
    <variable units="flux" name="D1" initial_value="0.94"/> 
    <variable units="nanomolar" name="L1" initial_value="0.3"/>
    <variable units="nanomolar" name="R1" initial_value="1.02"/>
    <variable units="nanomolar" name="A1" initial_value="0.45"/> 
    
    <variable units="nanomolar" public_interface="in" name="PT_n"/>
    <variable units="nanomolar" public_interface="in" name="CC_n"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    <variable units="dimensionless" public_interface="in" name="a"/>
    <variable units="dimensionless" public_interface="in" name="r"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Per_m</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> C1 </ci>
            <apply>
              <times/>
              <ci> S1 </ci>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> CC_n </ci>
                      <ci> A1 </ci>
                    </apply>
                    <ci> a </ci>
                  </apply>
                  <ci> B1 </ci>
                </apply>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> PT_n </ci>
                      <ci> R1 </ci>
                    </apply>
                    <ci> r </ci>
                  </apply>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> CC_n </ci>
                      <ci> A1 </ci>
                    </apply>
                    <ci> a </ci>
                  </apply>
                  <ci> B1 </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> D1 </ci>
              <apply>
                <divide/>
                <ci> Per_m </ci>
                <apply>
                  <plus/>
                  <ci> L1 </ci>
                  <ci> Per_m </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> Per_m </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Per_c" cmeta:id="Per_c">
    <rdf:RDF>
      <rdf:Description rdf:about="Per_c">
        <dc:title>Per_c</dc:title>
        <dcterms:alternative>cytosolic Per protein</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" name="Per_c" initial_value="0.6"/>
    
    <variable units="first_order_rate_constant" name="S2" initial_value="0.48"/>
    <variable units="first_order_rate_constant" name="D2" initial_value="0.44"/> 
    <variable units="nanomolar" name="L2" initial_value="0.2"/>
    <variable units="nanomolar" name="Dbt_c" initial_value="1.0"/>
    
    <variable units="second_order_rate_constant" public_interface="in" name="V1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="V2"/>
    <variable units="nanomolar" public_interface="in" name="Per_m"/>
    <variable units="nanomolar" public_interface="in" name="Tim_c"/>
    <variable units="nanomolar" public_interface="in" name="PT_c"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Per_c</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> S2 </ci>
              <ci> Per_m </ci>
            </apply>
            <apply>
              <times/>
              <ci> V2 </ci>
              <ci> PT_c </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> V1 </ci>
              <ci> Per_c </ci> 
              <ci> Tim_c </ci>
            </apply>
            <apply>
              <times/>
              <ci> D2 </ci>
              <ci> Dbt_c </ci>
              <apply>
                <divide/>
                <ci> Per_c </ci>
                <apply>
                  <plus/>
                  <ci> L2 </ci>
                  <ci> Per_c </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> Per_c </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Tim_m" cmeta:id="Tim_m">
    <rdf:RDF>
      <rdf:Description rdf:about="Tim_m">
        <dc:title>Tim_m</dc:title>
        <dcterms:alternative>Tim mRNA</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="Tim_m_Tim_m" name="Tim_m" initial_value="0.7"/>
    
    <variable units="dimensionless" name="B2" initial_value="0.0"/>
    <variable units="flux" name="C2" initial_value="0.0"/> 
    <variable units="flux" name="S3" initial_value="1.45"/>
    <variable units="flux" name="D3" initial_value="0.94"/> 
    <variable units="nanomolar" name="L3" initial_value="0.3"/>
    <variable units="nanomolar" name="R2" initial_value="1.02"/>
    <variable units="nanomolar" name="A2" initial_value="0.45"/> 
    
    <variable units="nanomolar" public_interface="in" name="PT_n"/>
    <variable units="nanomolar" public_interface="in" name="CC_n"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    <variable units="dimensionless" public_interface="in" name="a"/>
    <variable units="dimensionless" public_interface="in" name="r"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Tim_m</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> C2 </ci>
            <apply>
              <times/>
              <ci> S3 </ci>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> CC_n </ci>
                      <ci> A2 </ci>
                    </apply>
                    <ci> a </ci>
                  </apply>
                  <ci> B2 </ci>
                </apply>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> PT_n </ci>
                      <ci> R2 </ci>
                    </apply>
                    <ci> r </ci>
                  </apply>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> CC_n </ci>
                      <ci> A2 </ci>
                    </apply>
                    <ci> a </ci>
                  </apply>
                  <ci> B2 </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> D3 </ci>
              <apply>
                <divide/>
                <ci> Tim_m </ci>
                <apply>
                  <plus/>
                  <ci> L3 </ci>
                  <ci> Tim_m </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> Tim_m </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Tim_c" cmeta:id="Tim_c">
    <rdf:RDF>
      <rdf:Description rdf:about="Tim_c">
        <dc:title>Tim_c</dc:title>
        <dcterms:alternative>cytosolic Tim protein</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="Tim_c_Tim_c" name="Tim_c" initial_value="0.8"/>
    
    <variable units="first_order_rate_constant" name="S4" initial_value="0.48"/>
    <variable units="flux" name="D4" initial_value="0.44"/> 
    <variable units="nanomolar" name="L4" initial_value="0.2"/>
    
    <variable units="second_order_rate_constant" public_interface="in" name="V1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="V2"/>
    <variable units="nanomolar" public_interface="in" name="Tim_m"/>
    <variable units="nanomolar" public_interface="in" name="Per_c"/>
    <variable units="nanomolar" public_interface="in" name="PT_c"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Tim_c</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> S4 </ci>
              <ci> Tim_m </ci>
            </apply>
            <apply>
              <times/>
              <ci> V2 </ci>
              <ci> PT_c </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> V1 </ci>
              <ci> Per_c </ci> 
              <ci> Tim_c </ci>
            </apply>
            <apply>
              <times/>
              <ci> D4 </ci>
              <apply>
                <divide/>
                <ci> Tim_c </ci>
                <apply>
                  <plus/>
                  <ci> L4 </ci>
                  <ci> Tim_c </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> Tim_c </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="PT_c" cmeta:id="PT_c">
    <rdf:RDF>
      <rdf:Description rdf:about="PT_c">
        <dc:title>PT_c</dc:title>
        <dcterms:alternative>cytosolic Per-Tim protein complex</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="PT_c_PT_c" name="PT_c" initial_value="0.9"/>
    
    <variable units="flux" name="D5" initial_value="0.44"/> 
    <variable units="nanomolar" name="L5" initial_value="0.2"/>
    
    <variable units="nanomolar" public_interface="in" name="K1"/>
    <variable units="nanomolar" public_interface="in" name="K2"/>
    <variable units="flux" public_interface="in" name="T1"/>  
    <variable units="flux" public_interface="in" name="T2"/>
    <variable units="second_order_rate_constant" public_interface="in" name="V1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="V2"/>
    <variable units="nanomolar" public_interface="in" name="Tim_c"/>
    <variable units="nanomolar" public_interface="in" name="Per_c"/>
    <variable units="nanomolar" public_interface="in" name="PT_n"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>PT_c</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> V1 </ci>
              <ci> Per_c </ci> 
              <ci> Tim_c </ci>
            </apply>
            <apply>
              <times/>
              <ci> T2 </ci>
              <apply>
                <divide/>
                <ci> PT_n </ci>
                <apply>
                  <plus/>
                  <ci> K2 </ci>
                  <ci> PT_n </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> V2 </ci>
              <ci> PT_c </ci>
            </apply>
            <apply>
              <times/>
              <ci> T1 </ci>
              <apply>
                <divide/>
                <ci> PT_c </ci>
                <apply>
                  <plus/>
                  <ci> K1 </ci>
                  <ci> PT_c </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D5 </ci>
              <apply>
                <divide/>
                <ci> PT_c </ci>
                <apply>
                  <plus/>
                  <ci> L5 </ci>
                  <ci> PT_c </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> PT_c </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="PT_n" cmeta:id="PT_n">
    <rdf:RDF>
      <rdf:Description rdf:about="PT_n">
        <dc:title>PT_n</dc:title>
        <dcterms:alternative>nuclear Per-Tim protein complex</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="PT_n_PT_n" name="PT_n" initial_value="1.0"/>
    
    <variable units="flux" name="D6" initial_value="0.29"/> 
    <variable units="nanomolar" name="L6" initial_value="0.2"/>
    
    <variable units="nanomolar" public_interface="in" name="K1"/>
    <variable units="nanomolar" public_interface="in" name="K2"/>
    <variable units="flux" public_interface="in" name="T1"/>  
    <variable units="flux" public_interface="in" name="T2"/>
    <variable units="nanomolar" public_interface="in" name="PT_c"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>PT_n</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> T1 </ci>
            <apply>
              <divide/>
              <ci> PT_c </ci>
              <apply>
                <plus/>
                <ci> K1 </ci>
                <ci> PT_c </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> T2 </ci>
              <apply>
                <divide/>
                <ci> PT_n </ci>
                <apply>
                  <plus/>
                  <ci> K2 </ci>
                  <ci> PT_n </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D6 </ci>
              <apply>
                <divide/>
                <ci> PT_n </ci>
                <apply>
                  <plus/>
                  <ci> L6 </ci>
                  <ci> PT_n </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> PT_n </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
    
  <component name="Clk_m" cmeta:id="Clk_m">
    <rdf:RDF>
      <rdf:Description rdf:about="Clk_m">
        <dc:title>Clk_m</dc:title>
        <dcterms:alternative>Clk mRNA</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="Clk_m_Clk_m" name="Clk_m" initial_value="0.1"/>
    
    <variable units="dimensionless" name="B3" initial_value="0.6"/>
    <variable units="flux" name="C3" initial_value="0.0"/> 
    <variable units="flux" name="S5" initial_value="1.63"/>
    <variable units="flux" name="D7" initial_value="0.54"/> 
    <variable units="nanomolar" name="L7" initial_value="0.13"/>
    <variable units="nanomolar" name="R3" initial_value="0.89"/>
    <variable units="nanomolar" name="A3" initial_value="0.8"/> 
    
    <variable units="nanomolar" public_interface="in" name="CC_n"/>
    <variable units="nanomolar" public_interface="in" name="PT_n"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    <variable units="dimensionless" public_interface="in" name="a"/>
    <variable units="dimensionless" public_interface="in" name="r"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Clk_m</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> C3 </ci>
            <apply>
              <times/>
              <ci> S5 </ci>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> PT_n </ci>
                      <ci> A3 </ci>
                    </apply>
                    <ci> a </ci>
                  </apply>
                  <ci> B3 </ci>
                </apply>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> CC_n </ci>
                      <ci> R3 </ci>
                    </apply>
                    <ci> r </ci>
                  </apply>
                  <apply>
                    <power/>
                    <apply>
                      <divide/>
                      <ci> PT_n </ci>
                      <ci> A3 </ci>
                    </apply>
                    <ci> a </ci>
                  </apply>
                  <ci> B3 </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> D7 </ci>
              <apply>
                <divide/>
                <ci> Clk_m </ci>
                <apply>
                  <plus/>
                  <ci> L7 </ci>
                  <ci> Clk_m </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> Clk_m </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Clk_c" cmeta:id="Clk_c">
    <rdf:RDF>
      <rdf:Description rdf:about="Clk_c">
        <dc:title>Clk_c</dc:title>
        <dcterms:alternative>cytosolic Clk protein</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="Clk_c_Clk_c" name="Clk_c" initial_value="0.2"/>
    
    <variable units="first_order_rate_constant" name="S6" initial_value="0.47"/>
    <variable units="flux" name="D8" initial_value="0.6"/> 
    <variable units="nanomolar" name="L8" initial_value="0.2"/>
    
    <variable units="second_order_rate_constant" public_interface="in" name="V3"/>
    <variable units="first_order_rate_constant" public_interface="in" name="V4"/>
    <variable units="nanomolar" public_interface="in" name="Clk_m"/>
    <variable units="nanomolar" public_interface="in" cmeta:id="Cyc_c_Cyc_c" name="Cyc_c"/>
    <variable units="nanomolar" public_interface="in" name="CC_c"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Clk_c</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> S6 </ci>
              <ci> Clk_m </ci>
            </apply>
            <apply>
              <times/>
              <ci> V4 </ci>
              <ci> CC_c </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> V3 </ci>
              <ci> Clk_c </ci> 
              <ci> Cyc_c </ci>
            </apply>
            <apply>
              <times/>
              <ci> D8 </ci>
              <apply>
                <divide/>
                <ci> Clk_c </ci>
                <apply>
                  <plus/>
                  <ci> L8 </ci>
                  <ci> Clk_c </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> Clk_c </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Cyc_c" cmeta:id="Cyc_c">
    <rdf:RDF>
      <rdf:Description rdf:about="Cyc_c">
        <dc:title>Cyc_c</dc:title>
        <dcterms:alternative>cytosolic Cyc protein</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="Cyc_c_Cyc_c" name="Cyc_c" initial_value="1.0"/>
  </component>
  
  <component name="CC_c" cmeta:id="CC_c">
    <rdf:RDF>
      <rdf:Description rdf:about="CC_c">
        <dc:title>CC_c</dc:title>
        <dcterms:alternative>cytosolic Clk-Cyc protein complex</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="CC_c_CC_c" name="CC_c" initial_value="0.3"/>
    
    <variable units="flux" name="D9" initial_value="0.6"/> 
    <variable units="nanomolar" name="L9" initial_value="0.2"/>
    
    <variable units="nanomolar" public_interface="in" name="K3"/>
    <variable units="nanomolar" public_interface="in" name="K4"/>
    <variable units="flux" public_interface="in" name="T3"/>  
    <variable units="flux" public_interface="in" name="T4"/>
    <variable units="second_order_rate_constant" public_interface="in" name="V3"/>
    <variable units="first_order_rate_constant" public_interface="in" name="V4"/>
    <variable units="nanomolar" public_interface="in" name="Cyc_c"/>
    <variable units="nanomolar" public_interface="in" name="Clk_c"/>
    <variable units="nanomolar" public_interface="in" name="CC_n"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>CC_c</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> V3 </ci>
              <ci> Clk_c </ci> 
              <ci> Cyc_c </ci>
            </apply>
            <apply>
              <times/>
              <ci> T4 </ci>
              <apply>
                <divide/>
                <ci> CC_n </ci>
                <apply>
                  <plus/>
                  <ci> K4 </ci>
                  <ci> CC_n </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> V4 </ci>
              <ci> CC_c </ci>
            </apply>
            <apply>
              <times/>
              <ci> T3 </ci>
              <apply>
                <divide/>
                <ci> CC_c </ci>
                <apply>
                  <plus/>
                  <ci> K3 </ci>
                  <ci> CC_c </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D9 </ci>
              <apply>
                <divide/>
                <ci> CC_c </ci>
                <apply>
                  <plus/>
                  <ci> L9 </ci>
                  <ci> CC_c </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> CC_c </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="CC_n" cmeta:id="CC_n">
    <rdf:RDF>
      <rdf:Description rdf:about="CC_n">
        <dc:title>CC_n</dc:title>
        <dcterms:alternative>nuclear Clk-Cyc protein complex</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="nanomolar" public_interface="out" cmeta:id="CC_n_CC_n" name="CC_n" initial_value="0.4"/>
    
    <variable units="flux" name="D10" initial_value="0.3"/> 
    <variable units="nanomolar" name="L10" initial_value="0.2"/>
    
    <variable units="nanomolar" public_interface="in" name="K3"/>
    <variable units="nanomolar" public_interface="in" name="K4"/>
    <variable units="flux" public_interface="in" name="T3"/>  
    <variable units="flux" public_interface="in" name="T4"/>
    <variable units="nanomolar" public_interface="in" name="CC_c"/>
    <variable units="hour" public_interface="in" name="time"/>
    <variable units="first_order_rate_constant" public_interface="in" name="D0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>CC_n</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> T3 </ci>
            <apply>
              <divide/>
              <ci> CC_c </ci>
              <apply>
                <plus/>
                <ci> K3 </ci>
                <ci> CC_c </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> T4 </ci>
              <apply>
                <divide/>
                <ci> CC_n </ci>
                <apply>
                  <plus/>
                  <ci> K4 </ci>
                  <ci> CC_n </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D10 </ci>
              <apply>
                <divide/>
                <ci> CC_n </ci>
                <apply>
                  <plus/>
                  <ci> L10 </ci>
                  <ci> CC_n </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> D0 </ci>
              <ci> CC_n </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="parameters">
    <variable units="nanomolar" public_interface="out" name="K1" initial_value="2.0"/>
    <variable units="nanomolar" public_interface="out" name="K2" initial_value="2.0"/>
    <variable units="nanomolar" public_interface="out" name="K3" initial_value="2.0"/>
    <variable units="nanomolar" public_interface="out" name="K4" initial_value="2.0"/>
    <variable units="second_order_rate_constant" public_interface="out" name="V1" initial_value="1.45"/>
    <variable units="first_order_rate_constant" public_interface="out" name="V2" initial_value="1.45"/>
    <variable units="second_order_rate_constant" public_interface="out" name="V3" initial_value="1.63"/>
    <variable units="first_order_rate_constant" public_interface="out" name="V4" initial_value="1.63"/>
    <variable units="flux" public_interface="out" name="T1" initial_value="1.73"/>
    <variable units="flux" public_interface="out" name="T2" initial_value="0.72"/>
    <variable units="flux" public_interface="out" name="T3" initial_value="1.63"/>
    <variable units="flux" public_interface="out" name="T4" initial_value="0.52"/>
    <variable units="dimensionless" public_interface="out" name="a" initial_value="1.0"/>
    <variable units="dimensionless" public_interface="out" name="r" initial_value="4.0"/>
    <variable units="first_order_rate_constant" public_interface="out" name="D0" initial_value="0.012"/>
  </component>
  
  <connection>
    <map_components component_2="parameters" component_1="Per_m"/>
    <map_variables variable_2="a" variable_1="a"/>
    <map_variables variable_2="r" variable_1="r"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="parameters" component_1="Per_c"/>
    <map_variables variable_2="V1" variable_1="V1"/>
    <map_variables variable_2="V2" variable_1="V2"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="parameters" component_1="Tim_m"/>
    <map_variables variable_2="a" variable_1="a"/>
    <map_variables variable_2="r" variable_1="r"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="parameters" component_1="Tim_c"/>
    <map_variables variable_2="V1" variable_1="V1"/>
    <map_variables variable_2="V2" variable_1="V2"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="parameters" component_1="PT_c"/>
    <map_variables variable_2="V1" variable_1="V1"/>
    <map_variables variable_2="V2" variable_1="V2"/>
    <map_variables variable_2="T1" variable_1="T1"/>
    <map_variables variable_2="T2" variable_1="T2"/>
    <map_variables variable_2="K1" variable_1="K1"/>
    <map_variables variable_2="K2" variable_1="K2"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="parameters" component_1="PT_n"/>
    <map_variables variable_2="T1" variable_1="T1"/>
    <map_variables variable_2="T2" variable_1="T2"/>
    <map_variables variable_2="K1" variable_1="K1"/>
    <map_variables variable_2="K2" variable_1="K2"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="parameters" component_1="Clk_m"/>
    <map_variables variable_2="a" variable_1="a"/>
    <map_variables variable_2="r" variable_1="r"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="parameters" component_1="Clk_c"/>
    <map_variables variable_2="V3" variable_1="V3"/>
    <map_variables variable_2="V4" variable_1="V4"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="parameters" component_1="CC_c"/>
    <map_variables variable_2="V3" variable_1="V3"/>
    <map_variables variable_2="V4" variable_1="V4"/>
    <map_variables variable_2="T3" variable_1="T3"/>
    <map_variables variable_2="T4" variable_1="T4"/>
    <map_variables variable_2="K3" variable_1="K3"/>
    <map_variables variable_2="K4" variable_1="K4"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="parameters" component_1="CC_n"/>
    <map_variables variable_2="T3" variable_1="T3"/>
    <map_variables variable_2="T4" variable_1="T4"/>
    <map_variables variable_2="K3" variable_1="K3"/>
    <map_variables variable_2="K4" variable_1="K4"/>
    <map_variables variable_2="D0" variable_1="D0"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Per_m"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Per_c"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Tim_m"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Tim_c"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="PT_c"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="PT_n"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Clk_m"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Clk_c"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="CC_c"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="CC_n"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="CC_n" component_1="Per_m"/>
    <map_variables variable_2="CC_n" variable_1="CC_n"/>
  </connection>
  
  <connection>
    <map_components component_2="PT_n" component_1="Per_m"/>
    <map_variables variable_2="PT_n" variable_1="PT_n"/>
  </connection>
  
  <connection>
    <map_components component_2="Per_m" component_1="Per_c"/>
    <map_variables variable_2="Per_m" variable_1="Per_m"/>
  </connection>
  
  <connection>
    <map_components component_2="Tim_c" component_1="Per_c"/>
    <map_variables variable_2="Tim_c" variable_1="Tim_c"/>
    <map_variables variable_2="Per_c" variable_1="Per_c"/>
  </connection>
  
  <connection>
    <map_components component_2="PT_c" component_1="Per_c"/>
    <map_variables variable_2="PT_c" variable_1="PT_c"/>
    <map_variables variable_2="Per_c" variable_1="Per_c"/>
  </connection>
  
  <connection>
    <map_components component_2="CC_n" component_1="Tim_m"/>
    <map_variables variable_2="CC_n" variable_1="CC_n"/>
  </connection>
  
  <connection>
    <map_components component_2="PT_n" component_1="Tim_m"/>
    <map_variables variable_2="PT_n" variable_1="PT_n"/>
  </connection>
  
  <connection>
    <map_components component_2="PT_c" component_1="Tim_c"/>
    <map_variables variable_2="PT_c" variable_1="PT_c"/>
    <map_variables variable_2="Tim_c" variable_1="Tim_c"/>
  </connection>
  
  <connection>
    <map_components component_2="Tim_m" component_1="Tim_c"/>
    <map_variables variable_2="Tim_m" variable_1="Tim_m"/>
  </connection>
  
  <connection>
    <map_components component_2="PT_n" component_1="PT_c"/>
    <map_variables variable_2="PT_n" variable_1="PT_n"/>
    <map_variables variable_2="PT_c" variable_1="PT_c"/>
  </connection>
  
  <connection>
    <map_components component_2="PT_n" component_1="Clk_m"/>
    <map_variables variable_2="PT_n" variable_1="PT_n"/>
  </connection>
  
  <connection>
    <map_components component_2="CC_n" component_1="Clk_m"/>
    <map_variables variable_2="CC_n" variable_1="CC_n"/>
  </connection>
  
  <connection>
    <map_components component_2="Clk_m" component_1="Clk_c"/>
    <map_variables variable_2="Clk_m" variable_1="Clk_m"/>
  </connection>
  
  <connection>
    <map_components component_2="CC_c" component_1="Clk_c"/>
    <map_variables variable_2="CC_c" variable_1="CC_c"/>
    <map_variables variable_2="Clk_c" variable_1="Clk_c"/>
  </connection>
  
  <connection>
    <map_components component_2="Cyc_c" component_1="Clk_c"/>
    <map_variables variable_2="Cyc_c" variable_1="Cyc_c"/>
  </connection>
  
  <connection>
    <map_components component_2="Cyc_c" component_1="CC_c"/>
    <map_variables variable_2="Cyc_c" variable_1="Cyc_c"/>
  </connection>
  
  <connection>
    <map_components component_2="CC_n" component_1="CC_c"/>
    <map_variables variable_2="CC_n" variable_1="CC_n"/>
    <map_variables variable_2="CC_c" variable_1="CC_c"/>
  </connection>
  




<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <rdf:Bag rdf:about="rdf:#e31d78d3-89ba-4fd4-98c7-cc84cb7e442a">
    <rdf:li>circadian rhythms</rdf:li>
    <rdf:li>signal transduction</rdf:li>
    <rdf:li>oscillator</rdf:li>
    <rdf:li>circadian rhythm</rdf:li>
    <rdf:li>drosophila</rdf:li>
  </rdf:Bag>
  <rdf:Description rdf:about="rdf:#e31d78d3-89ba-4fd4-98c7-cc84cb7e442a"/>
  <rdf:Seq rdf:about="rdf:#adf93871-147a-49b7-8283-38f74fdad713">
    <rdf:li rdf:resource="rdf:#8e791b2b-0a31-46c6-ab25-9a82455bc05d"/>
    <rdf:li rdf:resource="rdf:#a5b81f9f-5948-4f07-8ca9-ee4d47cbbc7f"/>
    <rdf:li rdf:resource="rdf:#22ec0d5f-57ec-4ca8-9d69-7ce26cff4e56"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#60b96a2c-86e8-4db7-8183-46af24007ee1">
    <dcterms:W3CDTF>2001-06-21</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a51328f1-60f0-4536-810f-5f1130912397">
    <bqs:Pubmed_id>11403560</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#83d21e6c-2e49-4935-b790-4a3091922805"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$KY9Jj2">
    <ns7:boundIntervals rdf:resource="rdf:#$LY9Jj2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b3ca88cc-c02a-457a-8365-2e4a53ea287f">
    <vCard:Given>Hiroki</vCard:Given>
    <vCard:Family>Ueda</vCard:Family>
    <vCard:Other>R</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8f850fd5-4c1d-41e5-a7dd-e8d602534978">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#be517f4c-43cd-4510-a553-538de75df537">
    <vCard:Given>Hiroaki</vCard:Given>
    <vCard:Family>Kitano</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1eabd2fe-564a-4580-9127-20555a9242f6">
    <vCard:N rdf:resource="rdf:#937b3879-b2a0-47a7-a2ba-6513fe1ec065"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f9d6e962-db27-4f58-b85f-9cbe239a0e77">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2f22e793-d182-4565-954a-7792fe487f64">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a5b81f9f-5948-4f07-8ca9-ee4d47cbbc7f">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#5ae91878-fc85-479c-98c0-42ed45d94128"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0275ac6a-e378-4352-91e4-adc44722562c">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f6430b16-99b0-4f5e-a5aa-0966f356e5e1">
    <dc:creator rdf:resource="rdf:#2f22e793-d182-4565-954a-7792fe487f64"/>
    <rdf:value>This CellML model is known to run in both PCEnv and COR to replicate the published results.  The units have been checked and they are consistent.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$MY9Jj2">
    <ns7:endingValue>1000</ns7:endingValue>
    <ns7:nonstandard-pointDensity>10000</ns7:nonstandard-pointDensity>
    <ns7:maximumStepSize>0.1</ns7:maximumStepSize>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#80e5a22b-def2-4541-b546-9f3d95fa5b53">
    <dc:subject rdf:resource="rdf:#4d506f4a-f8ed-45f1-a371-c5f477f04372"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a27dff20-77a9-437f-9e36-46278e79d7bd">
    <dcterms:W3CDTF>2008-10-14T15:10:39+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8e791b2b-0a31-46c6-ab25-9a82455bc05d">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#b3ca88cc-c02a-457a-8365-2e4a53ea287f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5ae91878-fc85-479c-98c0-42ed45d94128">
    <vCard:Given>Masatoshi</vCard:Given>
    <vCard:Family>Hagiwara</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$LY9Jj2">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$MY9Jj2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4d506f4a-f8ed-45f1-a371-c5f477f04372">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#e31d78d3-89ba-4fd4-98c7-cc84cb7e442a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#22ec0d5f-57ec-4ca8-9d69-7ce26cff4e56">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#be517f4c-43cd-4510-a553-538de75df537"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#f6430b16-99b0-4f5e-a5aa-0966f356e5e1"/>
    <dcterms:created rdf:resource="rdf:#58bf27b5-7384-41e6-8f1d-0fc082cea50a"/>
    <dc:creator rdf:resource="rdf:#fe7855a9-24ac-4a34-b5aa-f49136820039"/>
    <cmeta:modification rdf:resource="rdf:#da4a6185-e2bd-4ef7-b383-6597522953e7"/>
    <cmeta:modification rdf:resource="rdf:#9b2b0aec-554b-4345-b89d-95902ef55786"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#da4a6185-e2bd-4ef7-b383-6597522953e7">
    <dcterms:modified rdf:resource="rdf:#eccb8b25-eb7e-4bc7-93a2-dd789480091e"/>
    <rdf:value>Added simulation metadata and cmeta:id's to some variables</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#1eabd2fe-564a-4580-9127-20555a9242f6"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#83d21e6c-2e49-4935-b790-4a3091922805">
    <dc:creator rdf:resource="rdf:#adf93871-147a-49b7-8283-38f74fdad713"/>
    <dc:title>Robust Oscillations within the Interlocked Feedback Model of Drosophila Circadian Rhythm</dc:title>
    <bqs:volume>210</bqs:volume>
    <bqs:first_page>401</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#f3f6bda1-6d0b-48d6-8c24-7e5b786c3a6e"/>
    <dcterms:issued rdf:resource="rdf:#60b96a2c-86e8-4db7-8183-46af24007ee1"/>
    <bqs:last_page>406</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#eccb8b25-eb7e-4bc7-93a2-dd789480091e">
    <dcterms:W3CDTF>2008-10-20T14:24:23+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#58337fc9-fba8-41f4-b1ca-f4cd92949d37">
    <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:#9b2b0aec-554b-4345-b89d-95902ef55786">
    <dcterms:modified rdf:resource="rdf:#a27dff20-77a9-437f-9e36-46278e79d7bd"/>
    <rdf:value>Added initial conditions to all the substrates (these were based on the values used in the corresponding SBML model in the BioModels database), and fixed a typo in the equation for dPTc/dt such that the first term in this equation now reads V1*Perc*Timc rather than V1*Perc*PTc.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#57ea31f7-5c8a-495c-a2e9-6c8ce2fd1858"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#937b3879-b2a0-47a7-a2ba-6513fe1ec065">
    <vCard:Given>James</vCard:Given>
    <vCard:Family>Lawson</vCard:Family>
    <vCard:Other>Richard</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f3f6bda1-6d0b-48d6-8c24-7e5b786c3a6e">
    <dc:title>Journal of Theoretical Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4adad4e7-8fce-44b0-af9c-a04168908601">
    <dc:creator rdf:resource="rdf:#992dde9a-388e-42b2-b701-85f1984698c3"/>
    <rdf:value>This is the CellML description of Ueda et al.'s mathematical model of circadian rhythms in Drosophila.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fe7855a9-24ac-4a34-b5aa-f49136820039">
    <vCard:ORG rdf:resource="rdf:#8f850fd5-4c1d-41e5-a7dd-e8d602534978"/>
    <vCard:EMAIL rdf:resource="rdf:#58337fc9-fba8-41f4-b1ca-f4cd92949d37"/>
    <vCard:N rdf:resource="rdf:#0275ac6a-e378-4352-91e4-adc44722562c"/>
  </rdf:Description>
  <rdf:Description rdf:about="#ueda_hagiwara_kitano_2001_version01">
    <ns7:simulation rdf:resource="rdf:#$KY9Jj2"/>
    <dc:title>Ueda et al.'s mathematical model of circadian rhythms in Drosophila.</dc:title>
    <cmeta:comment rdf:resource="rdf:#4adad4e7-8fce-44b0-af9c-a04168908601"/>
    <bqs:reference rdf:resource="rdf:#80e5a22b-def2-4541-b546-9f3d95fa5b53"/>
    <bqs:reference rdf:resource="rdf:#a51328f1-60f0-4536-810f-5f1130912397"/>
    <cmeta:species>Drosophila</cmeta:species>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#58bf27b5-7384-41e6-8f1d-0fc082cea50a">
    <dcterms:W3CDTF>2003-09-23T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#57ea31f7-5c8a-495c-a2e9-6c8ce2fd1858">
    <vCard:N rdf:resource="rdf:#f9d6e962-db27-4f58-b85f-9cbe239a0e77"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#992dde9a-388e-42b2-b701-85f1984698c3">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
</rdf:RDF>
</model>