Location: Leloup, Gonze, Goldbeter, 1999 @ bfaac0e80b23 / leloup_1999_1.1model.cellml

Author:
Mona Zhu <devnull@localhost>
Date:
2011-01-19 17:09:20+13:00
Desc:
Changed import locations
Permanent Source URI:
https://models.cellml.org/workspace/leloup_gonze_goldbeter_1999/rawfile/bfaac0e80b23726ffe05b02f98b3d1d01a2ee3b7/leloup_1999_1.1model.cellml

<?xml version="1.0"?>
<model xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.1#" cmeta:id="leloup_1999_1.1model" name="leloup_1999_1.1model">
	<!-- documentation -->
	<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>Mona</firstname>
			  <surname>Zhu</surname>
		<affiliation>
		  <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
		</affiliation>
	  </author>
	</articleinfo>
	  <section id="sec_status">
		<title>Model Status</title>
		<para>
				This model is a 1.1 modification of the orignal 1.0 model; it uses imports that are reusable in defining equations in the model. The units are consistent. This particular version of the model has been translated from equations 1a-1j. This file is known to run in OpenCell only. (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_1.1model.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 -->
	<units name="nanomolar">
		<unit prefix="nano" units="mole"/>
		<unit exponent="-1" units="litre"/>
	</units>

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

	<units name="flux">
		<unit units="nanomolar"/>
		<unit exponent="-1" units="hour"/>
	</units>

	<units name="first_order_rate_constant">
		<unit exponent="-1" units="hour"/>
	</units>

	<units name="second_order_rate_constant">
		<unit exponent="-1" units="nanomolar"/>
		<unit exponent="-1" units="hour"/>
	</units>

	<!-- components -->
	<component name="environment">
		<variable cmeta:id="environment_time" name="time" public_interface="out" units="hour"/>
	</component>
	<component cmeta:id="MP" name="MP">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		  <rdf:Description rdf:about="#MP">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">MP</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">Per mRNA</dcterms:alternative>  
		  </rdf:Description>
		</rdf:RDF> 
		<variable cmeta:id="MP_MP" name="MP" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="1" name="vsP" private_interface="out" units="flux"/>
		<variable initial_value="0.7" name="vmP" private_interface="out" units="flux"/> 
		<variable initial_value="0.2" name="KmP" private_interface="out" units="nanomolar"/>
		<variable initial_value="1.0" name="KIP" private_interface="out" units="nanomolar"/> 
		<variable name="kd" private_interface="out" public_interface="in" units="first_order_rate_constant"/> 
		<variable name="CN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
		<variable name="n" private_interface="out" public_interface="in" units="dimensionless"/>
		<variable initial_value="0.1" name="MP_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="P0" name="P0">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		  <rdf:Description rdf:about="#P0">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">P0</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">non-phopshorylated PER</dcterms:alternative>  
		  </rdf:Description>
		</rdf:RDF>
		<variable cmeta:id="P0_P0" name="P0" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.9" name="ksP" private_interface="out" units="first_order_rate_constant"/>
		<variable name="V1P" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V2P" private_interface="out" public_interface="in" units="flux"/>
		<variable name="K1P" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K2P" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kd" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="P1" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="MP" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
		<variable initial_value="0.1" name="P0_init" private_interface="out" units="nanomolar"/>
	</component> 
	<component cmeta:id="P1" name="P1">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
			<rdf:Description rdf:about="#P1">
				<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">P1</dc:title>
				<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">phopshorylated PER</dcterms:alternative>  
			</rdf:Description>
		</rdf:RDF>
		<variable cmeta:id="P1_P1" name="P1" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable name="V1P" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V2P" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V3P" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V4P" private_interface="out" public_interface="in" units="flux"/>
		<variable name="K1P" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K2P" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K3P" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K4P" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kd" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="P0" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="P2" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
		<variable initial_value="0.1" name="P1_init" private_interface="out" units="nanomolar"/>
	</component> 
	<component cmeta:id="P2" name="P2">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#P2">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">P2</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">phopshorylated PER</dcterms:alternative>  
		</rdf:Description>
		</rdf:RDF>
		<variable cmeta:id="P2_P2" name="P2" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="2" name="vdP" private_interface="out" units="flux"/>
		<variable initial_value="0.2" name="KdP" private_interface="out" units="nanomolar"/>
		<variable name="CN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable cmeta:id="P2_Pt" name="Pt" private_interface="in" units="nanomolar"/>
		<variable name="P0" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="V3P" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V4P" private_interface="out" public_interface="in" units="flux"/>
		<variable name="k3" private_interface="out" public_interface="in" units="second_order_rate_constant"/>
		<variable name="k4" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="K3P" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K4P" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kd" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="P1" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="T2" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="C" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
		<variable initial_value="0.1" name="P2_init" private_interface="out" units="nanomolar"/>
	</component> 
	<component cmeta:id="MT" name="MT">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		  <rdf:Description rdf:about="#MT">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">MT</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">Tim mRNA</dcterms:alternative>  
		  </rdf:Description>
		</rdf:RDF> 
		<variable cmeta:id="MT_MT" name="MT" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="1.0" name="vsT" private_interface="out" units="flux"/>
		<variable initial_value="0.7" name="vmT" private_interface="out" units="flux"/> 
		<variable initial_value="0.2" name="KmT" private_interface="out" units="nanomolar"/>
		<variable initial_value="1.0" name="KIT" private_interface="out" units="nanomolar"/> 
		<variable name="kd" private_interface="out" public_interface="in" units="first_order_rate_constant"/> 
		<variable name="CN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
		<variable name="n" private_interface="out" public_interface="in" units="dimensionless"/>
		<variable initial_value="1.6" name="MT_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="T0" name="T0">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		  <rdf:Description rdf:about="#T0">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">T0</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">non-phopshorylated TIM</dcterms:alternative>  
		  </rdf:Description>
		</rdf:RDF>
		<variable cmeta:id="T0_T0" name="T0" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.9" name="ksT" private_interface="out" units="first_order_rate_constant"/>
		<variable name="V1T" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V2T" private_interface="out" public_interface="in" units="flux"/>
		<variable name="K1T" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K2T" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kd" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="T1" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="MT" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
		<variable initial_value="0.1" name="T0_init" private_interface="out" units="nanomolar"/></component> 
  
	<component cmeta:id="T1" name="T1">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		  <rdf:Description rdf:about="#T1">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">T1</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">phopshorylated TIM</dcterms:alternative>  
		  </rdf:Description>
		</rdf:RDF>
		<variable cmeta:id="T1_T1" name="T1" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable name="V1T" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V2T" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V3T" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V4T" private_interface="out" public_interface="in" units="flux"/>
		<variable name="K1T" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K2T" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K3T" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K4T" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kd" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="T0" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="T2" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
		<variable initial_value="0.1" name="T1_init" private_interface="out" units="nanomolar"/>
	</component> 
	<component cmeta:id="T2" name="T2">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		  <rdf:Description rdf:about="#T2">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">T2</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">phopshorylated TIM</dcterms:alternative>  
		  </rdf:Description>
		</rdf:RDF>
		<variable cmeta:id="T2_T2" name="T2" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable cmeta:id="T2_vdT" initial_value="2" name="vdT" private_interface="out" units="flux"/>
		<variable initial_value="0.2" name="KdT" private_interface="out" units="nanomolar"/>
		<variable name="V3T" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V4T" private_interface="out" public_interface="in" units="flux"/>
		<variable name="k3" private_interface="out" public_interface="in" units="second_order_rate_constant"/>
		<variable name="k4" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="K3T" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="K4T" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kd" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="T1" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="P2" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="C" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
		<variable initial_value="0.1" name="T2_init" private_interface="out" units="nanomolar"/>
		<variable name="CN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Tt" private_interface="in" units="nanomolar"/>
		<variable name="T0" private_interface="out" public_interface="in" units="nanomolar"/>
	</component> 
	<component cmeta:id="C" name="C">
		<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		  <rdf:Description rdf:about="#C">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">C</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">cytosolic PER-TIM complex</dcterms:alternative>  
		  </rdf:Description>
		</rdf:RDF>
		<variable cmeta:id="C_C" initial_value="0.1" name="C" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.01" name="kdC" units="first_order_rate_constant"/>
		<variable name="k1" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k2" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k3" public_interface="in" units="second_order_rate_constant"/>
		<variable name="k4" public_interface="in" units="first_order_rate_constant"/>
		<variable name="T2" public_interface="in" units="nanomolar"/>
		<variable name="P2" public_interface="in" units="nanomolar"/>
		<variable name="CN" public_interface="in" units="nanomolar"/>
		<variable name="time" public_interface="in" units="hour"/>
		<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 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		  <rdf:Description rdf:about="#CN">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">CN</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">nuclear PER-TIM complex</dcterms:alternative>  
		  </rdf:Description>
		</rdf:RDF>
		<variable cmeta:id="CN_CN" initial_value="1.25" name="CN" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.01" name="kdN" units="first_order_rate_constant"/>
		<variable name="k1" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k2" public_interface="in" units="first_order_rate_constant"/>
		<variable name="C" public_interface="in" units="nanomolar"/>
		<variable name="time" public_interface="in" units="hour"/>
		<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 initial_value="0.6" name="k1" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.2" name="k2" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="1.2" name="k3" public_interface="out" units="second_order_rate_constant"/>
		<variable initial_value="0.6" name="k4" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.01" name="kd" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="8.0" name="V1P" public_interface="out" units="flux"/>
		<variable initial_value="8.0" name="V1T" public_interface="out" units="flux"/>
		<variable initial_value="1.0" name="V2P" public_interface="out" units="flux"/>
		<variable initial_value="1.0" name="V2T" public_interface="out" units="flux"/>
		<variable initial_value="8.0" name="V3P" public_interface="out" units="flux"/>
		<variable initial_value="8.0" name="V3T" public_interface="out" units="flux"/>
		<variable initial_value="1.0" name="V4P" public_interface="out" units="flux"/>
		<variable initial_value="1.0" name="V4T" public_interface="out" units="flux"/>
		<variable initial_value="2.0" name="K1P" public_interface="out" units="nanomolar"/>
		<variable initial_value="2.0" name="K1T" public_interface="out" units="nanomolar"/>
		<variable initial_value="2.0" name="K2P" public_interface="out" units="nanomolar"/>
		<variable initial_value="2.0" name="K2T" public_interface="out" units="nanomolar"/>
		<variable initial_value="2.0" name="K3P" public_interface="out" units="nanomolar"/>
		<variable initial_value="2.0" name="K3T" public_interface="out" units="nanomolar"/>
		<variable initial_value="1" name="K4P" public_interface="out" units="nanomolar"/>
		<variable initial_value="1" name="K4T" public_interface="out" units="nanomolar"/>
		<variable initial_value="4.0" name="n" public_interface="out" units="dimensionless"/>
	</component>
	
	<!-- groups -->
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="MP">
				<component_ref component="M_P"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="P0">
				<component_ref component="P_0"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="P1">
				<component_ref component="P_1"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="P2">
				<component_ref component="P_2"/>
				<component_ref component="P_t"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="MT">
				<component_ref component="M_T"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="T0">
				<component_ref component="T_0"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="T1">
				<component_ref component="T_1"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="T2">
				<component_ref component="T_2"/>
				<component_ref component="T_t"/>
			</component_ref>
	</group>
	
	<!-- imports -->
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law1.cellml">
		<component component_ref="interface" name="M_P"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law2.cellml">
		<component component_ref="interface" name="P_0"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law3.cellml">
		<component component_ref="interface" name="P_1"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law4.cellml">
		<component component_ref="interface" name="P_2"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law1.cellml">
		<component component_ref="interface" name="M_T"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law2.cellml">
		<component component_ref="interface" name="T_0"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law3.cellml">
		<component component_ref="interface" name="T_1"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law4.cellml">
		<component component_ref="interface" name="T_2"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="total.cellml">
		<component component_ref="interface" name="P_t"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="total.cellml">
		<component component_ref="interface" name="T_t"/>
	</import>
	
	<!-- connections -->
	<connection>
		<map_components component_1="MP" component_2="parameters"/>
		<map_variables variable_1="n" variable_2="n"/>
		<map_variables variable_1="kd" variable_2="kd"/>
	</connection>
	<connection>
		<map_components component_1="MT" component_2="parameters"/>
		<map_variables variable_1="n" variable_2="n"/>
		<map_variables variable_1="kd" variable_2="kd"/>
	</connection>
	<connection>
		<map_components component_1="P0" component_2="parameters"/>
		<map_variables variable_1="V1P" variable_2="V1P"/>
		<map_variables variable_1="V2P" variable_2="V2P"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="K1P" variable_2="K1P"/>
		<map_variables variable_1="K2P" variable_2="K2P"/>
	</connection>
	<connection>
		<map_components component_1="P1" component_2="parameters"/>
		<map_variables variable_1="V1P" variable_2="V1P"/>
		<map_variables variable_1="V2P" variable_2="V2P"/>
		<map_variables variable_1="V3P" variable_2="V3P"/>
		<map_variables variable_1="V4P" variable_2="V4P"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="K1P" variable_2="K1P"/>
		<map_variables variable_1="K2P" variable_2="K2P"/>
		<map_variables variable_1="K3P" variable_2="K3P"/>
		<map_variables variable_1="K4P" variable_2="K4P"/>
	</connection>
	<connection>
		<map_components component_1="P2" component_2="parameters"/>
		<map_variables variable_1="V3P" variable_2="V3P"/>
		<map_variables variable_1="V4P" variable_2="V4P"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="K3P" variable_2="K3P"/>
		<map_variables variable_1="K4P" variable_2="K4P"/>
		<map_variables variable_1="k3" variable_2="k3"/>
		<map_variables variable_1="k4" variable_2="k4"/>
	</connection>
	<connection>
		<map_components component_1="T0" component_2="parameters"/>
		<map_variables variable_1="V1T" variable_2="V1T"/>
		<map_variables variable_1="V2T" variable_2="V2T"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="K1T" variable_2="K1T"/>
		<map_variables variable_1="K2T" variable_2="K2T"/>
	</connection>
	<connection>
		<map_components component_1="T1" component_2="parameters"/>
		<map_variables variable_1="V1T" variable_2="V1T"/>
		<map_variables variable_1="V2T" variable_2="V2T"/>
		<map_variables variable_1="V3T" variable_2="V3T"/>
		<map_variables variable_1="V4T" variable_2="V4T"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="K1T" variable_2="K1T"/>
		<map_variables variable_1="K2T" variable_2="K2T"/>
		<map_variables variable_1="K3T" variable_2="K3T"/>
		<map_variables variable_1="K4T" variable_2="K4T"/>
	</connection>
	<connection>
		<map_components component_1="T2" component_2="parameters"/>
		<map_variables variable_1="V3T" variable_2="V3T"/>
		<map_variables variable_1="V4T" variable_2="V4T"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="K3T" variable_2="K3T"/>
		<map_variables variable_1="K4T" variable_2="K4T"/>
		<map_variables variable_1="k3" variable_2="k3"/>
		<map_variables variable_1="k4" variable_2="k4"/>
	</connection>
	<connection>
		<map_components component_1="C" component_2="parameters"/>
		<map_variables variable_1="k1" variable_2="k1"/>
		<map_variables variable_1="k2" variable_2="k2"/>
		<map_variables variable_1="k3" variable_2="k3"/>
		<map_variables variable_1="k4" variable_2="k4"/>
	</connection>
	<connection>
		<map_components component_1="CN" component_2="parameters"/>
		<map_variables variable_1="k1" variable_2="k1"/>
		<map_variables variable_1="k2" variable_2="k2"/>
	</connection>
	<connection>
		<map_components component_1="MP" component_2="CN"/>
		<map_variables variable_1="CN" variable_2="CN"/>
	</connection>
	<connection>
		<map_components component_1="MT" component_2="CN"/>
		<map_variables variable_1="CN" variable_2="CN"/>
	</connection>
	<connection>
		<map_components component_1="P0" component_2="MP"/>
		<map_variables variable_1="MP" variable_2="MP"/>
	</connection>
	<connection>
		<map_components component_1="P0" component_2="P1"/>
		<map_variables variable_1="P1" variable_2="P1"/>
		<map_variables variable_1="P0" variable_2="P0"/>
	</connection>
	<connection>
		<map_components component_1="P1" component_2="P2"/>
		<map_variables variable_1="P2" variable_2="P2"/>
		<map_variables variable_1="P1" variable_2="P1"/>
	</connection>
	<connection>
		<map_components component_1="P2" component_2="T2"/>
		<map_variables variable_1="T2" variable_2="T2"/>
		<map_variables variable_1="P2" variable_2="P2"/>
	</connection>
	<connection>
		<map_components component_1="P2" component_2="C"/>
		<map_variables variable_1="C" variable_2="C"/>
		<map_variables variable_1="P2" variable_2="P2"/>
	</connection>
	<connection>
		<map_components component_1="T0" component_2="MT"/>
		<map_variables variable_1="MT" variable_2="MT"/>
	</connection>
	<connection>
		<map_components component_1="T0" component_2="T1"/>
		<map_variables variable_1="T1" variable_2="T1"/>
		<map_variables variable_1="T0" variable_2="T0"/>
	</connection>
	<connection>
		<map_components component_1="T1" component_2="T2"/>
		<map_variables variable_1="T2" variable_2="T2"/>
		<map_variables variable_1="T1" variable_2="T1"/>
	</connection>
	<connection>
		<map_components component_1="T2" component_2="C"/>
		<map_variables variable_1="C" variable_2="C"/>
		<map_variables variable_1="T2" variable_2="T2"/>
	</connection>
	<connection>
		<map_components component_1="C" component_2="CN"/>
		<map_variables variable_1="C" variable_2="C"/>
		<map_variables variable_1="CN" variable_2="CN"/>
	</connection>
	<connection>
		<map_components component_1="MP" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="MT" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="P0" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="P1" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="P2" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="T0" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="T1" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="T2" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="C" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="CN" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>
	<connection>
		<map_components component_1="P2" component_2="CN"/>
		<map_variables variable_1="CN" variable_2="CN"/>
	</connection>
	<connection>
		<map_components component_1="P0" component_2="P2"/>
		<map_variables variable_1="P0" variable_2="P0"/>
	</connection>
	<connection>
		<map_components component_1="M_P" component_2="MP"/>
		<map_variables variable_1="M" variable_2="MP"/>
		<map_variables variable_1="vs" variable_2="vsP"/>
		<map_variables variable_1="vm" variable_2="vmP"/>
		<map_variables variable_1="Km" variable_2="KmP"/>
		<map_variables variable_1="KI" variable_2="KIP"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="CN" variable_2="CN"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="n" variable_2="n"/>
		<map_variables variable_1="M_init" variable_2="MP_init"/>
	</connection>
	<connection>
		<map_components component_1="P_0" component_2="P0"/>
		<map_variables variable_1="Y0" variable_2="P0"/>
		<map_variables variable_1="ks" variable_2="ksP"/>
		<map_variables variable_1="V1" variable_2="V1P"/>
		<map_variables variable_1="V2" variable_2="V2P"/>
		<map_variables variable_1="K1" variable_2="K1P"/>
		<map_variables variable_1="K2" variable_2="K2P"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="Y1" variable_2="P1"/>
		<map_variables variable_1="M" variable_2="MP"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="Y0_init" variable_2="P0_init"/>
	</connection>
	<connection>
		<map_components component_1="P_1" component_2="P1"/>
		<map_variables variable_1="Y1" variable_2="P1"/>
		<map_variables variable_1="V1" variable_2="V1P"/>
		<map_variables variable_1="V2" variable_2="V2P"/>
		<map_variables variable_1="V3" variable_2="V3P"/>
		<map_variables variable_1="V4" variable_2="V4P"/>
		<map_variables variable_1="K1" variable_2="K1P"/>
		<map_variables variable_1="K2" variable_2="K2P"/>
		<map_variables variable_1="K3" variable_2="K3P"/>
		<map_variables variable_1="K4" variable_2="K4P"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="Y0" variable_2="P0"/>
		<map_variables variable_1="Y2" variable_2="P2"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="Y1_init" variable_2="P1_init"/>
	</connection>
	<connection>
		<map_components component_1="P_2" component_2="P2"/>
		<map_variables variable_1="Y2" variable_2="P2"/>
		<map_variables variable_1="vd" variable_2="vdP"/>
		<map_variables variable_1="KdP" variable_2="KdP"/>
		<map_variables variable_1="V3" variable_2="V3P"/>
		<map_variables variable_1="V4" variable_2="V4P"/>
		<map_variables variable_1="k3" variable_2="k3"/>
		<map_variables variable_1="k4" variable_2="k4"/>
		<map_variables variable_1="K3" variable_2="K3P"/>
		<map_variables variable_1="K4" variable_2="K4P"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="Y1" variable_2="P1"/>
		<map_variables variable_1="Z2" variable_2="T2"/>
		<map_variables variable_1="C" variable_2="C"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="Y2_init" variable_2="P2_init"/>
	</connection>
	<connection>
		<map_components component_1="M_T" component_2="MT"/>
		<map_variables variable_1="M" variable_2="MT"/>
		<map_variables variable_1="vs" variable_2="vsT"/>
		<map_variables variable_1="vm" variable_2="vmT"/>
		<map_variables variable_1="Km" variable_2="KmT"/>
		<map_variables variable_1="KI" variable_2="KIT"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="CN" variable_2="CN"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="n" variable_2="n"/>
		<map_variables variable_1="M_init" variable_2="MT_init"/>
	</connection>
	<connection>
		<map_components component_1="T_0" component_2="T0"/>
		<map_variables variable_1="Y0" variable_2="T0"/>
		<map_variables variable_1="ks" variable_2="ksT"/>
		<map_variables variable_1="V1" variable_2="V1T"/>
		<map_variables variable_1="V2" variable_2="V2T"/>
		<map_variables variable_1="K1" variable_2="K1T"/>
		<map_variables variable_1="K2" variable_2="K2T"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="Y1" variable_2="T1"/>
		<map_variables variable_1="M" variable_2="MT"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="Y0_init" variable_2="T0_init"/>
	</connection>
	<connection>
		<map_components component_1="T_1" component_2="T1"/>
		<map_variables variable_1="Y1" variable_2="T1"/>
		<map_variables variable_1="V1" variable_2="V1T"/>
		<map_variables variable_1="V2" variable_2="V2T"/>
		<map_variables variable_1="V3" variable_2="V3T"/>
		<map_variables variable_1="V4" variable_2="V4T"/>
		<map_variables variable_1="K1" variable_2="K1T"/>
		<map_variables variable_1="K2" variable_2="K2T"/>
		<map_variables variable_1="K3" variable_2="K3T"/>
		<map_variables variable_1="K4" variable_2="K4T"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="Y0" variable_2="T0"/>
		<map_variables variable_1="Y2" variable_2="T2"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="Y1_init" variable_2="T1_init"/>
	</connection>
	<connection>
		<map_components component_1="T_2" component_2="T2"/>
		<map_variables variable_1="Y2" variable_2="T2"/>
		<map_variables variable_1="vd" variable_2="vdT"/>
		<map_variables variable_1="KdP" variable_2="KdT"/>
		<map_variables variable_1="V3" variable_2="V3T"/>
		<map_variables variable_1="V4" variable_2="V4T"/>
		<map_variables variable_1="k3" variable_2="k3"/>
		<map_variables variable_1="k4" variable_2="k4"/>
		<map_variables variable_1="K3" variable_2="K3T"/>
		<map_variables variable_1="K4" variable_2="K4T"/>
		<map_variables variable_1="kd" variable_2="kd"/>
		<map_variables variable_1="Y1" variable_2="T1"/>
		<map_variables variable_1="Z2" variable_2="P2"/>
		<map_variables variable_1="C" variable_2="C"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="Y2_init" variable_2="T2_init"/>
	</connection>
	<connection>
		<map_components component_1="P_t" component_2="P2"/>
		<map_variables variable_1="Y0" variable_2="P0"/>
		<map_variables variable_1="Y1" variable_2="P1"/>
		<map_variables variable_1="C" variable_2="C"/>
		<map_variables variable_1="CN" variable_2="CN"/>
		<map_variables variable_1="Yt" variable_2="Pt"/>
	</connection>
	<connection>
		<map_components component_1="T_t" component_2="T2"/>
		<map_variables variable_1="Y0" variable_2="T0"/>
		<map_variables variable_1="Y1" variable_2="T1"/>
		<map_variables variable_1="C" variable_2="C"/>
		<map_variables variable_1="CN" variable_2="CN"/>
		<map_variables variable_1="Yt" variable_2="Tt"/>
	</connection>
	<connection>
		<map_components component_1="T2" component_2="CN"/>
		<map_variables variable_1="CN" variable_2="CN"/>
	</connection>
	<connection>
		<map_components component_1="T2" component_2="T0"/>
		<map_variables variable_1="T0" variable_2="T0"/>
	</connection>
	<connection>
		<map_components component_1="P_t" component_2="P_2"/>
		<map_variables variable_1="Y2" variable_2="Y2"/>
	</connection>
	<connection>
		<map_components component_1="T_t" component_2="T_2"/>
		<map_variables variable_1="Y2" variable_2="Y2"/>
	</connection>
   
   	<!-- metadata -->
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
				<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<vCard:Family>Zhu</vCard:Family>
					<vCard:Given>Xin Yue</vCard:Given>
					<vCard:Other>Mona</vCard:Other>
				</vCard:N>
				<vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<rdf:value>xzhu055@aucklanduni.ac.nz</rdf:value>
					<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
				</vCard:EMAIL>
				<vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<vCard:Orgname>The University of Auckland</vCard:Orgname>
					<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
				</vCard:ORG>
			</dc:creator>
			<dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
				<dcterms:W3CDTF>06-01-2011</dcterms:W3CDTF>
			</dcterms:created>
		</rdf:Description>
		<rdf:Description rdf:about="#leloup_1999_1.1">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Limit Cycle Models for Circadian Rhythms Based on Transcriptional Regulation in Drosophila and Neurospora</dc:title>
			<cmeta:comment rdf:parseType="Resource">
				<rdf:value>This is the CellML description of Leloup et al.'s 1999 mathematical model of circadian rhythms in Drosophila.
				</rdf:value>
				<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
					<vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Mona Zhu</vCard:FN>
				</dc:creator>
			</cmeta:comment>
			<cmeta:species>Drosophila</cmeta:species>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.1#" rdf:parseType="Resource">
				<dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
					<bqs:subject_type>keyword</bqs:subject_type>
					<rdf:value>
						<rdf:Bag>
							<rdf:li>Circadian Rhythms</rdf:li>
							<rdf:li>Signal Transduction</rdf:li>
							<rdf:li>Drosophilia</rdf:li>
						</rdf:Bag>
					</rdf:value>
				</dc:subject>
			</bqs:reference>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.1#" rdf:parseType="Resource">
				<bqs:Pubmed_id>10643740</bqs:Pubmed_id>
				<bqs:JournalArticle rdf:parseType="Resource">
					<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/">
						<rdf:Seq>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Leloup</vCard:Family>
										<vCard:Given>Jean-Christophe</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Gonze</vCard:Family>
										<vCard:Given>Didier</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Albert</vCard:Family>
										<vCard:Given>Goldbeter</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
						</rdf:Seq>
					</dc:creator>
					<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Limit Cycle Models for Circadian Rhythms Based on Transcriptional Regulation in Drosophila and Neurospora</dc:title>
					<bqs:Journal rdf:parseType="Resource">
						<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Journal of Biological Rhythms</dc:title>
					</bqs:Journal>
					<bqs:volume>14</bqs:volume>
					<bqs:first_page>433</bqs:first_page>
					<bqs:last_page>448</bqs:last_page>
				</bqs:JournalArticle>
			</bqs:reference>
		</rdf:Description>
	</rdf:RDF>
</model>