- Author:
- Catherine Lloyd <c.lloyd@auckland.ac.nz>
- Date:
- 2010-04-29 10:23:33+12:00
- Desc:
- Updated documention and curation status. Fixed RDF. Balanced units. Added initial conditions.
- Permanent Source URI:
- https://models.cellml.org/workspace/leloup_gonze_goldbeter_1999/rawfile/948987fe7cacf149d299752e5c7c7b096587f32a/leloup_gonze_goldbeter_1999_a.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: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#" cmeta:id="leloup_1999" name="leloup_1999">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Limit cycle models for circadian rhythms based on transcriptional regulation in Drosophila and Neurospora</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This CellML model runs in COR and OpenCell to reproduce the published results (figure 2a). The units are consistent. This particular version of the model has been translated from equations 1a-1j (Drosophila).
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: We examine theoretical models for circadian oscillations based on transcriptional regulation in Drosophila and Neurospora. For Drosophila, the molecular model is based on the negative feedback exerted on the expression of the per and tim genes by the complex formed between the PER and TIM proteins. For Neurospora, similarly, the model relies on the feedback exerted on the expression of the frq gene by its protein product FRQ. In both models, sustained rhythmic variations in protein and mRNA levels occur in continuous darkness, in the form of limit cycle oscillations. The effect of light on circadian rhythms is taken into account in the models by considering that it triggers degradation of the TIM protein in Drosophila, and frq transcription in Neurospora. When incorporating the control exerted by light at the molecular level, we show that the models can account for the entrainment of circadian rhythms by light-dark cycles and for the damping of the oscillations in constant light, though such damping occurs more readily in the Drosophila model. The models account for the phase shifts induced by light pulses and allow the construction of phase response curves. These compare well with experimental results obtained in Drosophila. The model for Drosophila shows that when applied at the appropriate phase, light pulses of appropriate duration and magnitude can permanently or transiently suppress circadian rhythmicity. We investigate the effects of the magnitude of light-induced changes on oscillatory behavior. Finally, we discuss the common and distinctive features of circadian oscillations in the two organisms.
</para>
<para>
The original paper reference is cited below:
</para>
<para>
Limit Cycle Models for Circadian Rhythms Based on Transcriptional Regulation in <emphasis>Drosophila</emphasis> and <emphasis>Neurospora</emphasis>, Jean-Christophe Leloup, Didier Gonze, and Albert Goldbeter, 1999, <emphasis>Journal of Biological Rhythms</emphasis>, 14, 433-448. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/10643740">PubMed ID: 10643740</ulink>
</para>
<informalfigure float="0" id="fig_reaction_diagram1">
<mediaobject>
<imageobject>
<objectinfo>
<title>reaction diagram</title>
</objectinfo>
<imagedata fileref="leloup_1999.png"/>
</imageobject>
</mediaobject>
<caption>Scheme for the model for circadian oscillations in <emphasis>Drosophila</emphasis> involving negative regulation of gene expression by the PER-TIM protein complex. And beneath this is the scheme for the model for circadian rhythms in <emphasis>Neurospora</emphasis>. This model is based on negative feedback exerted by the protein FRQ on the transcription of the <emphasis>frq</emphasis> gene.</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" cmeta:id="environment_time" name="time"/>
</component>
<!--
The following components describe all the reactants and products involved in
reactions.
-->
<component cmeta:id="MP" name="MP">
<rdf:RDF>
<rdf:Description rdf:about="#MP">
<dc:title>MP</dc:title>
<dcterms:alternative>Per mRNA</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="MP_MP" name="MP" initial_value="0.1"/>
<variable units="flux" name="vsP" initial_value="1"/>
<variable units="flux" name="vmP" initial_value="0.7"/>
<variable units="nanomolar" name="KmP" initial_value="0.2"/>
<variable units="nanomolar" name="KIP" initial_value="1.0"/>
<variable units="first_order_rate_constant" public_interface="in" name="kd"/>
<variable units="nanomolar" public_interface="in" name="CN"/>
<variable units="hour" public_interface="in" name="time"/>
<variable units="dimensionless" public_interface="in" name="n"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>MP</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> vsP </ci>
<apply>
<divide/>
<apply>
<power/>
<ci> KIP </ci>
<ci> n </ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci> KIP </ci>
<ci> n </ci>
</apply>
<apply>
<power/>
<ci> CN </ci>
<ci> n </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> vmP </ci>
<apply>
<divide/>
<ci> MP </ci>
<apply>
<plus/>
<ci> KmP </ci>
<ci> MP </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kd </ci>
<ci> MP </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="P0" name="P0">
<rdf:RDF>
<rdf:Description rdf:about="#P0">
<dc:title>P0</dc:title>
<dcterms:alternative>non-phopshorylated PER</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="P0_P0" name="P0" initial_value="0.1"/>
<variable units="first_order_rate_constant" name="ksP" initial_value="0.9"/>
<variable units="flux" public_interface="in" name="V1P"/>
<variable units="flux" public_interface="in" name="V2P"/>
<variable units="nanomolar" public_interface="in" name="K1P"/>
<variable units="nanomolar" public_interface="in" name="K2P"/>
<variable units="first_order_rate_constant" public_interface="in" name="kd"/>
<variable units="nanomolar" public_interface="in" name="P1"/>
<variable units="nanomolar" public_interface="in" name="MP"/>
<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>P0</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> ksP </ci>
<ci> MP </ci>
</apply>
<apply>
<times/>
<ci> V2P </ci>
<apply>
<divide/>
<ci> P1 </ci>
<apply>
<plus/>
<ci> K2P </ci>
<ci> P1 </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> V1P </ci>
<apply>
<divide/>
<ci> P0 </ci>
<apply>
<plus/>
<ci> K1P </ci>
<ci> P0 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kd </ci>
<ci> P0 </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="P1" name="P1">
<rdf:RDF>
<rdf:Description rdf:about="#P1">
<dc:title>P1</dc:title>
<dcterms:alternative>phopshorylated PER</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="P1_P1" name="P1" initial_value="0.1"/>
<variable units="flux" public_interface="in" name="V1P"/>
<variable units="flux" public_interface="in" name="V2P"/>
<variable units="flux" public_interface="in" name="V3P"/>
<variable units="flux" public_interface="in" name="V4P"/>
<variable units="nanomolar" public_interface="in" name="K1P"/>
<variable units="nanomolar" public_interface="in" name="K2P"/>
<variable units="nanomolar" public_interface="in" name="K3P"/>
<variable units="nanomolar" public_interface="in" name="K4P"/>
<variable units="first_order_rate_constant" public_interface="in" name="kd"/>
<variable units="nanomolar" public_interface="in" name="P0"/>
<variable units="nanomolar" public_interface="in" name="P2"/>
<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>P1</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> V1P </ci>
<apply>
<divide/>
<ci> P0 </ci>
<apply>
<plus/>
<ci> K1P </ci>
<ci> P0 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> V4P </ci>
<apply>
<divide/>
<ci> P2 </ci>
<apply>
<plus/>
<ci> K4P </ci>
<ci> P2 </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> V2P </ci>
<apply>
<divide/>
<ci> P1 </ci>
<apply>
<plus/>
<ci> K2P </ci>
<ci> P1 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> V3P </ci>
<apply>
<divide/>
<ci> P1 </ci>
<apply>
<plus/>
<ci> K3P </ci>
<ci> P1 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kd </ci>
<ci> P1 </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="P2" name="P2">
<rdf:RDF>
<rdf:Description rdf:about="#P2">
<dc:title>P2</dc:title>
<dcterms:alternative>phopshorylated PER</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="P2_P2" name="P2" initial_value="0.1"/>
<variable units="flux" name="vdP" initial_value="2"/>
<variable units="nanomolar" name="KdP" initial_value="0.2"/>
<variable units="nanomolar" public_interface="in" name="CN"/>
<variable units="nanomolar" cmeta:id="P2_Pt" name="Pt"/>
<variable units="nanomolar" public_interface="in" name="P0"/>
<variable units="flux" public_interface="in" name="V3P"/>
<variable units="flux" public_interface="in" name="V4P"/>
<variable units="second_order_rate_constant" public_interface="in" name="k3"/>
<variable units="first_order_rate_constant" public_interface="in" name="k4"/>
<variable units="nanomolar" public_interface="in" name="K3P"/>
<variable units="nanomolar" public_interface="in" name="K4P"/>
<variable units="first_order_rate_constant" public_interface="in" name="kd"/>
<variable units="nanomolar" public_interface="in" name="P1"/>
<variable units="nanomolar" public_interface="in" name="T2"/>
<variable units="nanomolar" public_interface="in" name="C"/>
<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>P2</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> V3P </ci>
<apply>
<divide/>
<ci> P1 </ci>
<apply>
<plus/>
<ci> K3P </ci>
<ci> P1 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> k4 </ci>
<ci> C </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> V4P </ci>
<apply>
<divide/>
<ci> P2 </ci>
<apply>
<plus/>
<ci> K4P </ci>
<ci> P2 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> k3 </ci>
<ci> P2 </ci>
<ci> T2 </ci>
</apply>
<apply>
<times/>
<ci> vdP </ci>
<apply>
<divide/>
<ci> P2 </ci>
<apply>
<plus/>
<ci> KdP </ci>
<ci> P2 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kd </ci>
<ci> P2 </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Pt</ci>
<apply>
<plus/>
<ci>P0</ci>
<ci>P1</ci>
<ci>P2</ci>
<ci>C</ci>
<ci>CN</ci>
</apply>
</apply>
</math>
</component>
<component cmeta:id="MT" name="MT">
<rdf:RDF>
<rdf:Description rdf:about="#MT">
<dc:title>MT</dc:title>
<dcterms:alternative>Tim mRNA</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="MT_MT" name="MT" initial_value="1.6"/>
<variable units="flux" name="vsT" initial_value="1.0"/>
<variable units="flux" name="vmT" initial_value="0.7"/>
<variable units="nanomolar" name="KmT" initial_value="0.2"/>
<variable units="nanomolar" name="KIT" initial_value="1.0"/>
<variable units="first_order_rate_constant" public_interface="in" name="kd"/>
<variable units="nanomolar" public_interface="in" name="CN"/>
<variable units="hour" public_interface="in" name="time"/>
<variable units="dimensionless" public_interface="in" name="n"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>MT</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> vsT </ci>
<apply>
<divide/>
<apply>
<power/>
<ci> KIT </ci>
<ci> n </ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci> KIT </ci>
<ci> n </ci>
</apply>
<apply>
<power/>
<ci> CN </ci>
<ci> n </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> vmT </ci>
<apply>
<divide/>
<ci> MT </ci>
<apply>
<plus/>
<ci> KmT </ci>
<ci> MT </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kd </ci>
<ci> MT </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="T0" name="T0">
<rdf:RDF>
<rdf:Description rdf:about="#T0">
<dc:title>T0</dc:title>
<dcterms:alternative>non-phopshorylated TIM</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="T0_T0" name="T0" initial_value="0.1"/>
<variable units="first_order_rate_constant" name="ksT" initial_value="0.9"/>
<variable units="flux" public_interface="in" name="V1T"/>
<variable units="flux" public_interface="in" name="V2T"/>
<variable units="nanomolar" public_interface="in" name="K1T"/>
<variable units="nanomolar" public_interface="in" name="K2T"/>
<variable units="first_order_rate_constant" public_interface="in" name="kd"/>
<variable units="nanomolar" public_interface="in" name="T1"/>
<variable units="nanomolar" public_interface="in" name="MT"/>
<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>T0</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> ksT </ci>
<ci> MT </ci>
</apply>
<apply>
<times/>
<ci> V2T </ci>
<apply>
<divide/>
<ci> T1 </ci>
<apply>
<plus/>
<ci> K2T </ci>
<ci> T1 </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> V1T </ci>
<apply>
<divide/>
<ci> T0 </ci>
<apply>
<plus/>
<ci> K1T </ci>
<ci> T0 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kd </ci>
<ci> T0 </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="T1" name="T1">
<rdf:RDF>
<rdf:Description rdf:about="#T1">
<dc:title>T1</dc:title>
<dcterms:alternative>phopshorylated TIM</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="T1_T1" name="T1" initial_value="0.1"/>
<variable units="flux" public_interface="in" name="V1T"/>
<variable units="flux" public_interface="in" name="V2T"/>
<variable units="flux" public_interface="in" name="V3T"/>
<variable units="flux" public_interface="in" name="V4T"/>
<variable units="nanomolar" public_interface="in" name="K1T"/>
<variable units="nanomolar" public_interface="in" name="K2T"/>
<variable units="nanomolar" public_interface="in" name="K3T"/>
<variable units="nanomolar" public_interface="in" name="K4T"/>
<variable units="first_order_rate_constant" public_interface="in" name="kd"/>
<variable units="nanomolar" public_interface="in" name="T0"/>
<variable units="nanomolar" public_interface="in" name="T2"/>
<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>T1</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> V1T </ci>
<apply>
<divide/>
<ci> T0 </ci>
<apply>
<plus/>
<ci> K1T </ci>
<ci> T0 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> V4T </ci>
<apply>
<divide/>
<ci> T2 </ci>
<apply>
<plus/>
<ci> K4T </ci>
<ci> T2 </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> V2T </ci>
<apply>
<divide/>
<ci> T1 </ci>
<apply>
<plus/>
<ci> K2T </ci>
<ci> T1 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> V3T </ci>
<apply>
<divide/>
<ci> T1 </ci>
<apply>
<plus/>
<ci> K3T </ci>
<ci> T1 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kd </ci>
<ci> T1 </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="T2" name="T2">
<rdf:RDF>
<rdf:Description rdf:about="#T2">
<dc:title>T2</dc:title>
<dcterms:alternative>phopshorylated TIM</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="T2_T2" name="T2" initial_value="0.1"/>
<variable units="flux" cmeta:id="T2_vdT" name="vdT" initial_value="2"/>
<variable units="nanomolar" name="KdT" initial_value="0.2"/>
<variable units="flux" public_interface="in" name="V3T"/>
<variable units="flux" public_interface="in" name="V4T"/>
<variable units="second_order_rate_constant" public_interface="in" name="k3"/>
<variable units="first_order_rate_constant" public_interface="in" name="k4"/>
<variable units="nanomolar" public_interface="in" name="K3T"/>
<variable units="nanomolar" public_interface="in" name="K4T"/>
<variable units="first_order_rate_constant" public_interface="in" name="kd"/>
<variable units="nanomolar" public_interface="in" name="T1"/>
<variable units="nanomolar" public_interface="in" name="P2"/>
<variable units="nanomolar" public_interface="in" name="C"/>
<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>T2</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> V3T </ci>
<apply>
<divide/>
<ci> T1 </ci>
<apply>
<plus/>
<ci> K3T </ci>
<ci> T1 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> k4 </ci>
<ci> C </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> V4T </ci>
<apply>
<divide/>
<ci> T2 </ci>
<apply>
<plus/>
<ci> K4T </ci>
<ci> T2 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> k3 </ci>
<ci> P2 </ci>
<ci> T2 </ci>
</apply>
<apply>
<times/>
<ci> vdT </ci>
<apply>
<divide/>
<ci> T2 </ci>
<apply>
<plus/>
<ci> KdT </ci>
<ci> T2 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kd </ci>
<ci> T2 </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="C" name="C">
<rdf:RDF>
<rdf:Description rdf:about="#C">
<dc:title>C</dc:title>
<dcterms:alternative>cytosolic PER-TIM complex</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="C_C" name="C" initial_value="0.1"/>
<variable units="first_order_rate_constant" name="kdC" initial_value="0.01"/>
<variable units="first_order_rate_constant" public_interface="in" name="k1"/>
<variable units="first_order_rate_constant" public_interface="in" name="k2"/>
<variable units="second_order_rate_constant" public_interface="in" name="k3"/>
<variable units="first_order_rate_constant" public_interface="in" name="k4"/>
<variable units="nanomolar" public_interface="in" name="T2"/>
<variable units="nanomolar" public_interface="in" name="P2"/>
<variable units="nanomolar" public_interface="in" name="CN"/>
<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>C</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> k3 </ci>
<ci> P2 </ci>
<ci> T2 </ci>
</apply>
<apply>
<times/>
<ci> k2 </ci>
<ci> CN </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> k4 </ci>
<ci> C </ci>
</apply>
<apply>
<times/>
<ci> k1 </ci>
<ci> C </ci>
</apply>
<apply>
<times/>
<ci> kdC </ci>
<ci> C </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="CN" name="CN">
<rdf:RDF>
<rdf:Description rdf:about="#CN">
<dc:title>CN</dc:title>
<dcterms:alternative>nuclear PER-TIM complex</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="nanomolar" public_interface="out" cmeta:id="CN_CN" name="CN" initial_value="1.25"/>
<variable units="first_order_rate_constant" name="kdN" initial_value="0.01"/>
<variable units="first_order_rate_constant" public_interface="in" name="k1"/>
<variable units="first_order_rate_constant" public_interface="in" name="k2"/>
<variable units="nanomolar" public_interface="in" name="C"/>
<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>CN</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k1 </ci>
<ci> C </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> k2 </ci>
<ci> CN </ci>
</apply>
<apply>
<times/>
<ci> kdN </ci>
<ci> CN </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!--
The following component defines the parameters of the model.
-->
<component name="parameters">
<variable units="first_order_rate_constant" public_interface="out" name="k1" initial_value="0.6"/>
<variable units="first_order_rate_constant" public_interface="out" name="k2" initial_value="0.2"/>
<variable units="second_order_rate_constant" public_interface="out" name="k3" initial_value="1.2"/>
<variable units="first_order_rate_constant" public_interface="out" name="k4" initial_value="0.6"/>
<variable units="first_order_rate_constant" public_interface="out" name="kd" initial_value="0.01"/>
<variable units="flux" public_interface="out" name="V1P" initial_value="8.0"/>
<variable units="flux" public_interface="out" name="V1T" initial_value="8.0"/>
<variable units="flux" public_interface="out" name="V2P" initial_value="1.0"/>
<variable units="flux" public_interface="out" name="V2T" initial_value="1.0"/>
<variable units="flux" public_interface="out" name="V3P" initial_value="8.0"/>
<variable units="flux" public_interface="out" name="V3T" initial_value="8.0"/>
<variable units="flux" public_interface="out" name="V4P" initial_value="1.0"/>
<variable units="flux" public_interface="out" name="V4T" initial_value="1.0"/>
<variable units="nanomolar" public_interface="out" name="K1P" initial_value="2.0"/>
<variable units="nanomolar" public_interface="out" name="K1T" initial_value="2.0"/>
<variable units="nanomolar" public_interface="out" name="K2P" initial_value="2.0"/>
<variable units="nanomolar" public_interface="out" name="K2T" initial_value="2.0"/>
<variable units="nanomolar" public_interface="out" name="K3P" initial_value="2.0"/>
<variable units="nanomolar" public_interface="out" name="K3T" initial_value="2.0"/>
<variable units="nanomolar" public_interface="out" name="K4P" initial_value="1"/>
<variable units="nanomolar" public_interface="out" name="K4T" initial_value="1"/>
<variable units="dimensionless" public_interface="out" name="n" initial_value="4.0"/>
</component>
<connection>
<map_components component_2="parameters" component_1="MP"/>
<map_variables variable_2="n" variable_1="n"/>
<map_variables variable_2="kd" variable_1="kd"/>
</connection>
<connection>
<map_components component_2="parameters" component_1="MT"/>
<map_variables variable_2="n" variable_1="n"/>
<map_variables variable_2="kd" variable_1="kd"/>
</connection>
<connection>
<map_components component_2="parameters" component_1="P0"/>
<map_variables variable_2="V1P" variable_1="V1P"/>
<map_variables variable_2="V2P" variable_1="V2P"/>
<map_variables variable_2="kd" variable_1="kd"/>
<map_variables variable_2="K1P" variable_1="K1P"/>
<map_variables variable_2="K2P" variable_1="K2P"/>
</connection>
<connection>
<map_components component_2="parameters" component_1="P1"/>
<map_variables variable_2="V1P" variable_1="V1P"/>
<map_variables variable_2="V2P" variable_1="V2P"/>
<map_variables variable_2="V3P" variable_1="V3P"/>
<map_variables variable_2="V4P" variable_1="V4P"/>
<map_variables variable_2="kd" variable_1="kd"/>
<map_variables variable_2="K1P" variable_1="K1P"/>
<map_variables variable_2="K2P" variable_1="K2P"/>
<map_variables variable_2="K3P" variable_1="K3P"/>
<map_variables variable_2="K4P" variable_1="K4P"/>
</connection>
<connection>
<map_components component_2="parameters" component_1="P2"/>
<map_variables variable_2="V3P" variable_1="V3P"/>
<map_variables variable_2="V4P" variable_1="V4P"/>
<map_variables variable_2="kd" variable_1="kd"/>
<map_variables variable_2="K3P" variable_1="K3P"/>
<map_variables variable_2="K4P" variable_1="K4P"/>
<map_variables variable_2="k3" variable_1="k3"/>
<map_variables variable_2="k4" variable_1="k4"/>
</connection>
<connection>
<map_components component_2="parameters" component_1="T0"/>
<map_variables variable_2="V1T" variable_1="V1T"/>
<map_variables variable_2="V2T" variable_1="V2T"/>
<map_variables variable_2="kd" variable_1="kd"/>
<map_variables variable_2="K1T" variable_1="K1T"/>
<map_variables variable_2="K2T" variable_1="K2T"/>
</connection>
<connection>
<map_components component_2="parameters" component_1="T1"/>
<map_variables variable_2="V1T" variable_1="V1T"/>
<map_variables variable_2="V2T" variable_1="V2T"/>
<map_variables variable_2="V3T" variable_1="V3T"/>
<map_variables variable_2="V4T" variable_1="V4T"/>
<map_variables variable_2="kd" variable_1="kd"/>
<map_variables variable_2="K1T" variable_1="K1T"/>
<map_variables variable_2="K2T" variable_1="K2T"/>
<map_variables variable_2="K3T" variable_1="K3T"/>
<map_variables variable_2="K4T" variable_1="K4T"/>
</connection>
<connection>
<map_components component_2="parameters" component_1="T2"/>
<map_variables variable_2="V3T" variable_1="V3T"/>
<map_variables variable_2="V4T" variable_1="V4T"/>
<map_variables variable_2="kd" variable_1="kd"/>
<map_variables variable_2="K3T" variable_1="K3T"/>
<map_variables variable_2="K4T" variable_1="K4T"/>
<map_variables variable_2="k3" variable_1="k3"/>
<map_variables variable_2="k4" variable_1="k4"/>
</connection>
<connection>
<map_components component_2="parameters" component_1="C"/>
<map_variables variable_2="k1" variable_1="k1"/>
<map_variables variable_2="k2" variable_1="k2"/>
<map_variables variable_2="k3" variable_1="k3"/>
<map_variables variable_2="k4" variable_1="k4"/>
</connection>
<connection>
<map_components component_2="parameters" component_1="CN"/>
<map_variables variable_2="k1" variable_1="k1"/>
<map_variables variable_2="k2" variable_1="k2"/>
</connection>
<connection>
<map_components component_2="CN" component_1="MP"/>
<map_variables variable_2="CN" variable_1="CN"/>
</connection>
<connection>
<map_components component_2="CN" component_1="MT"/>
<map_variables variable_2="CN" variable_1="CN"/>
</connection>
<connection>
<map_components component_2="MP" component_1="P0"/>
<map_variables variable_2="MP" variable_1="MP"/>
</connection>
<connection>
<map_components component_2="P1" component_1="P0"/>
<map_variables variable_2="P1" variable_1="P1"/>
<map_variables variable_2="P0" variable_1="P0"/>
</connection>
<connection>
<map_components component_2="P2" component_1="P1"/>
<map_variables variable_2="P2" variable_1="P2"/>
<map_variables variable_2="P1" variable_1="P1"/>
</connection>
<connection>
<map_components component_2="T2" component_1="P2"/>
<map_variables variable_2="T2" variable_1="T2"/>
<map_variables variable_2="P2" variable_1="P2"/>
</connection>
<connection>
<map_components component_2="C" component_1="P2"/>
<map_variables variable_2="C" variable_1="C"/>
<map_variables variable_2="P2" variable_1="P2"/>
</connection>
<connection>
<map_components component_2="MT" component_1="T0"/>
<map_variables variable_2="MT" variable_1="MT"/>
</connection>
<connection>
<map_components component_2="T1" component_1="T0"/>
<map_variables variable_2="T1" variable_1="T1"/>
<map_variables variable_2="T0" variable_1="T0"/>
</connection>
<connection>
<map_components component_2="T2" component_1="T1"/>
<map_variables variable_2="T2" variable_1="T2"/>
<map_variables variable_2="T1" variable_1="T1"/>
</connection>
<connection>
<map_components component_2="C" component_1="T2"/>
<map_variables variable_2="C" variable_1="C"/>
<map_variables variable_2="T2" variable_1="T2"/>
</connection>
<connection>
<map_components component_2="CN" component_1="C"/>
<map_variables variable_2="C" variable_1="C"/>
<map_variables variable_2="CN" variable_1="CN"/>
</connection>
<connection>
<map_components component_2="environment" component_1="MP"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="MT"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="P0"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="P1"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="P2"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="T0"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="T1"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="T2"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="C"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="CN"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="CN" component_1="P2"/>
<map_variables variable_2="CN" variable_1="CN"/>
</connection>
<connection>
<map_components component_2="P2" component_1="P0"/>
<map_variables variable_2="P0" variable_1="P0"/>
</connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Bag rdf:about="rdf:#2e2e8fb8-bb4b-496e-aa0f-d4061c62c2ba">
<rdf:li>circadian rhythms</rdf:li>
<rdf:li>signal transduction</rdf:li>
<rdf:li>drosophila</rdf:li>
</rdf:Bag>
<rdf:Seq rdf:about="rdf:#7da8af5e-7768-4b36-8330-8dcf3407b4c8">
<rdf:li rdf:resource="rdf:#b700623c-41da-4a37-ab0e-9a4af0f40668"/>
<rdf:li rdf:resource="rdf:#0a901c2a-c066-4e78-8cc1-9b940798b36b"/>
<rdf:li rdf:resource="rdf:#42cf7418-0743-4cba-8e83-95f091c5dede"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#1f9fa075-16b8-4a5a-915b-dae80cf14249">
<dcterms:W3CDTF>2003-09-24T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a3daf125-aa71-4b55-927e-1b34c8629dfb">
<vCard:Given>Albert</vCard:Given>
<vCard:Family>Goldbeter</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#923bfd0e-42c0-4d29-b031-bf5452611928">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:title>Limit Cycle Model for Circadian Rhythms in Drosophila and Neurospora (Drosophila model)</dc:title>
<dc:publisher>The University of Auckland, Auckland Bioengineering Institute</dc:publisher>
<cmeta:comment rdf:resource="rdf:#e1711b7c-891f-4a80-ab5a-21e0384cadd7"/>
<dcterms:created rdf:resource="rdf:#1f9fa075-16b8-4a5a-915b-dae80cf14249"/>
<dc:creator rdf:resource="rdf:#fcc7eb5e-1e73-406f-b11e-4c2ae5a6ebbe"/>
<cmeta:modification rdf:resource="rdf:#09e8c147-d534-4089-b129-3b42f92771b6"/>
<cmeta:modification rdf:resource="rdf:#87cc02ee-19c4-4b37-b155-aa9bbfcc9ca9"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#535a3d2b-57db-4085-a7c2-814d58e6fa1f">
<dcterms:W3CDTF>1999-12</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#72fd8854-409d-4535-80f5-f9922799f8be">
<bqs:Pubmed_id>10643740</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#bd0260e9-57a8-4df6-862e-612a3c11d625"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$RUs082">
<ns7:endingValue>100</ns7:endingValue>
<ns7:nonstandard-pointDensity>10000</ns7:nonstandard-pointDensity>
<ns7:maximumStepSize>0.1</ns7:maximumStepSize>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fcc7eb5e-1e73-406f-b11e-4c2ae5a6ebbe">
<vCard:ORG rdf:resource="rdf:#923bfd0e-42c0-4d29-b031-bf5452611928"/>
<vCard:EMAIL rdf:resource="rdf:#2a2663ff-79d2-431c-bb37-14c36ae5227e"/>
<vCard:N rdf:resource="rdf:#10c236f3-dca5-4336-926b-fd10da654f58"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#55ff19ef-a8ec-40f2-b8e0-c0ba6226c894">
<vCard:FN/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#62835309-4d36-4dbf-93b4-1396945103b0">
<vCard:Given>Jean-Christophe</vCard:Given>
<vCard:Family>Leloup</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2d6616c6-94da-4ace-95f9-9c3297743d0f">
<vCard:N rdf:resource="rdf:#4afe6db4-a1f0-40b0-8c0a-0c19b1d77d17"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#60918037-2dad-48eb-9023-ff585ed3c9df">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value rdf:resource="rdf:#2e2e8fb8-bb4b-496e-aa0f-d4061c62c2ba"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#bd0260e9-57a8-4df6-862e-612a3c11d625">
<dc:creator rdf:resource="rdf:#7da8af5e-7768-4b36-8330-8dcf3407b4c8"/>
<dc:title>Limit Cycle Models for Circadian Rhythms Based on Transcriptional Regulation in Drosophila and Neurospora</dc:title>
<bqs:volume>14</bqs:volume>
<bqs:first_page>433</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#67eb2311-6273-4465-9459-8760347e051b"/>
<dcterms:issued rdf:resource="rdf:#535a3d2b-57db-4085-a7c2-814d58e6fa1f"/>
<bqs:last_page>448</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$QUs082">
<rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
<rdf:first rdf:resource="rdf:#$RUs082"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#bd6761e1-7cb7-4bd8-b45a-9d1ded78f593">
<dc:subject rdf:resource="rdf:#60918037-2dad-48eb-9023-ff585ed3c9df"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0d2fc15c-ed23-4f60-95b1-b38283330348">
<dcterms:W3CDTF>2008-09-05T13:19:24+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2a2663ff-79d2-431c-bb37-14c36ae5227e">
<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:#f226f477-8730-46bd-9d2f-a97f3f899a25">
<vCard:Given>Didier</vCard:Given>
<vCard:Family>Gonze</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d321c9f4-3542-46a3-be9b-5d9729d110d2">
<dc:creator rdf:resource="rdf:#4e079499-d35f-4c5f-92b4-f600c0cb4a56"/>
<rdf:value>This is the CellML description of Leloup et al.'s 1999 mathematical model of circadian rhythms in Drosophila.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#09e8c147-d534-4089-b129-3b42f92771b6">
<dcterms:modified rdf:resource="rdf:#3db5e8b8-e91a-43cb-9fe6-d658657d171c"/>
<rdf:value>Corrected some initial conditions according to Dagmar Kohn's keen spotting!</rdf:value>
<cmeta:modifier rdf:resource="rdf:#2d6616c6-94da-4ace-95f9-9c3297743d0f"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b700623c-41da-4a37-ab0e-9a4af0f40668">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#62835309-4d36-4dbf-93b4-1396945103b0"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5fa52bda-ae3c-4e2f-a69d-fc8c025a5be5">
<vCard:N rdf:resource="rdf:#f28f9094-c62f-4d6c-918d-484258f9d922"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4afe6db4-a1f0-40b0-8c0a-0c19b1d77d17">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#87cc02ee-19c4-4b37-b155-aa9bbfcc9ca9">
<dcterms:modified rdf:resource="rdf:#0d2fc15c-ed23-4f60-95b1-b38283330348"/>
<rdf:value>Added equation for Pt to allow replication of publication figure 1.
Added simulation metadata.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#5fa52bda-ae3c-4e2f-a69d-fc8c025a5be5"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0a901c2a-c066-4e78-8cc1-9b940798b36b">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#f226f477-8730-46bd-9d2f-a97f3f899a25"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e1711b7c-891f-4a80-ab5a-21e0384cadd7">
<dc:creator rdf:resource="rdf:#55ff19ef-a8ec-40f2-b8e0-c0ba6226c894"/>
<rdf:value/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$PUs082">
<ns7:boundIntervals rdf:resource="rdf:#$QUs082"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f28f9094-c62f-4d6c-918d-484258f9d922">
<vCard:Given>James</vCard:Given>
<vCard:Family>Lawson</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4e079499-d35f-4c5f-92b4-f600c0cb4a56">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3db5e8b8-e91a-43cb-9fe6-d658657d171c">
<dcterms:W3CDTF>2009-04-14T13:51:23+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#42cf7418-0743-4cba-8e83-95f091c5dede">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#a3daf125-aa71-4b55-927e-1b34c8629dfb"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#10c236f3-dca5-4336-926b-fd10da654f58">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#67eb2311-6273-4465-9459-8760347e051b">
<dc:title>Journal of Biological Rhythms</dc:title>
</rdf:Description>
<rdf:Description rdf:about="#leloup_1999">
<ns7:simulation rdf:resource="rdf:#$PUs082"/>
<dc:title>Leloup et al.'s 1999 mathematical model of circadian rhythms in
Drosophila.</dc:title>
<cmeta:comment rdf:resource="rdf:#d321c9f4-3542-46a3-be9b-5d9729d110d2"/>
<bqs:reference rdf:resource="rdf:#72fd8854-409d-4535-80f5-f9922799f8be"/>
<bqs:reference rdf:resource="rdf:#bd6761e1-7cb7-4bd8-b45a-9d1ded78f593"/>
<cmeta:species>Drosophila</cmeta:species>
</rdf:Description>
</rdf:RDF>
</model>
