Location: Leloup, Goldbeter, 2003 @ 3f64359c54c2 / leloup_goldbeter_2003_1.1model.cellml

Author:
Mona Zhu <devnull@localhost>
Date:
2011-01-19 16:58:12+13:00
Desc:
Changed import locations
Permanent Source URI:
https://models.physiomeproject.org/workspace/leloup_goldbeter_2003/rawfile/3f64359c54c20162714e02f2361f792c1cf005dd/leloup_goldbeter_2003_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_goldbeter_2003_1.1model" name="leloup_goldbeter_2003_1.1model">
	<!-- documentation -->
	<documentation xmlns="http://cellml.org/tmp-documentation">
	<article>
		<articleinfo>
			<title>A Detailed Computational Model for the Mammalian Circadian Clock</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 CellML model runs in OpenCell to recreate the published results (figure 2A from the published paper, where KAC = 0.6nM). The 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 model needs to be run for 160 hours with a step size of 0.1 hours, and MP, MB and MC are plotted. Ignore the first couple of oscillations as the model needs to stabilise. Please note that the model presented here is the extended version of the model which includes equations to define the role of REV-ERB-alpha in inhibiting the transcription of Bmal1.  There are 19 equations and the parameters are listed as part of the Supplementary data set under Figure 8.</para>
		</section>
		<sect1 id="sec_structure">
		<title>Model Structure</title>

		<para>
		ABSTRACT: We present a computational model for the mammalian circadian clock based on the intertwined positive and negative regulatory loops involving the Per, Cry, Bmal1, Clock, and Rev-Erb alpha genes. In agreement with experimental observations, the model can give rise to sustained circadian oscillations in continuous darkness, characterized by an antiphase relationship between Per/Cry/Rev-Erbalpha and Bmal1 mRNAs. Sustained oscillations correspond to the rhythms autonomously generated by suprachiasmatic nuclei. For other parameter values, damped oscillations can also be obtained in the model. These oscillations, which transform into sustained oscillations when coupled to a periodic signal, correspond to rhythms produced by peripheral tissues. When incorporating the light-induced expression of the Per gene, the model accounts for entrainment of the oscillations by light-dark cycles. Simulations show that the phase of the oscillations can then vary by several hours with relatively minor changes in parameter values. Such a lability of the phase could account for physiological disorders related to circadian rhythms in humans, such as advanced or delayed sleep phase syndrome, whereas the lack of entrainment by light-dark cycles can be related to the non-24h sleep-wake syndrome. The model uncovers the possible existence of multiple sources of oscillatory behavior. Thus, in conditions where the indirect negative autoregulation of Per and Cry expression is inoperative, the model indicates the possibility that sustained oscillations might still arise from the negative autoregulation of Bmal1 expression.
		</para>

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

		<para>
		Toward a detailed computational model for the mammalian circadian clock, Jean-Christophe Leloup and Albert Goldbeter, 2003,<emphasis>PNAS</emphasis>, 100, 7051-7056.   <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=12775757&amp;dopt=Abstract">PubMed ID: 12775757</ulink>
		</para>

			<informalfigure float="0" id="fig_reaction_diagram">
				<mediaobject>
					<imageobject>
						<objectinfo>
						  <title>reaction diagram</title>
						</objectinfo>
					<imagedata fileref="leloup_2003.png"/>
					</imageobject>
				</mediaobject>
				<caption>Model for circadian oscillations in mammals involving interlocked negative and positive regulations of <emphasis>Per</emphasis>, <emphasis>Cry</emphasis>, <emphasis>Bmal1</emphasis>, and <emphasis>Rev-Erb</emphasis> genes by their protein products.</caption>
			</informalfigure>

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

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

	<units name="nanomolar">
		<unit prefix="nano" units="mole"/>
		<unit exponent="-1" units="litre"/>
	</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">
		<variable cmeta:id="MP_MP" name="MP" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="MP_init" private_interface="out" units="nanomolar"/>
		<variable initial_value="2.4" name="vsP" private_interface="out" units="flux"/>
		<variable initial_value="2.2" name="vmP" private_interface="out" units="flux"/>
		<variable initial_value="0.02" name="kdmp" private_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.6" name="KAP" private_interface="out" units="nanomolar"/>
		<variable initial_value="0.3" name="KmP" private_interface="out" units="nanomolar"/>
		<variable name="n" private_interface="out" public_interface="in" units="dimensionless"/>
		<variable name="BN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="MC" name="MC">  
		<variable cmeta:id="MC_MC" name="MC" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="1.2" name="MC_init" private_interface="out" units="nanomolar"/>
		<variable initial_value="2.2" name="vsC" private_interface="out" units="flux"/>
		<variable initial_value="2.0" name="vmC" private_interface="out" units="flux"/>
		<variable initial_value="0.02" name="kdmc" private_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.6" name="KAC" private_interface="out" units="nanomolar"/>
		<variable initial_value="0.4" name="KmC" private_interface="out" units="nanomolar"/>
		<variable name="n" private_interface="out" public_interface="in" units="dimensionless"/>
		<variable name="BN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="MB" name="MB">  
		<variable cmeta:id="MB_MB" name="MB" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="9" name="MB_init" private_interface="out" units="nanomolar"/>
		<variable initial_value="1.8" name="vsB" private_interface="out" units="flux"/>
		<variable initial_value="1.3" name="vmB" private_interface="out" units="flux"/>
		<variable initial_value="0.02" name="kdmb" private_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="2.2" name="KIB" private_interface="out" units="nanomolar"/>
		<variable initial_value="0.4" name="KmB" private_interface="out" units="nanomolar"/>
		<variable name="m" private_interface="out" public_interface="in" units="dimensionless"/>
		<variable name="RN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="MR" name="MR">
		<variable cmeta:id="MR_MR" name="MR" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="1.5" name="MR_init" private_interface="out" units="nanomolar"/>
		<variable initial_value="1.6" name="vsR" private_interface="out" units="flux"/>
		<variable initial_value="1.6" name="vmR" private_interface="out" units="flux"/>
		<variable initial_value="0.02" name="kdmr" private_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.6" name="KAR" private_interface="out" units="nanomolar"/>
		<variable initial_value="0.4" name="KmR" private_interface="out" units="nanomolar"/>
		<variable name="h" private_interface="out" public_interface="in" units="dimensionless"/>
		<variable name="BN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="PC" name="PC">
		<variable cmeta:id="PC_PC" name="PC" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="PC_init" private_interface="out" units="nanomolar"/>
		<variable name="ksP" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/>
		<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="kdn" private_interface="out" public_interface="in" 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="PCP" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="PCC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="CC" 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"/>
	</component>
	<component cmeta:id="CC" name="CC">
		<variable cmeta:id="CC_CC" name="CC" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="CC_init" private_interface="out" units="nanomolar"/>
		<variable name="ksC" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/>
		<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="kdnc" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V1C" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V2C" private_interface="out" public_interface="in" units="flux"/>
		<variable name="CCP" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="PC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="PCC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="MC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="RC" name="RC">
		<variable cmeta:id="RC_RC" initial_value="0.1" name="RC" public_interface="out" units="nanomolar"/>
		<variable name="ksR" public_interface="in" units="first_order_rate_constant"/>
		<variable name="Kd" public_interface="in" units="nanomolar"/>
		<variable name="k9" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k10" public_interface="in" units="first_order_rate_constant"/> 
		<variable name="kdn" public_interface="in" units="first_order_rate_constant"/>
		<variable name="vdRC" public_interface="in" units="flux"/>
		<variable name="RN" public_interface="in" units="nanomolar"/>
		<variable name="MR" 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> RC </ci> 
			</apply>
			<apply>
			  <minus/>
			  <apply>
				<plus/>
				<apply>
				  <times/>
				  <ci> ksR </ci>
				  <ci> MR </ci>
				</apply>
				<apply>
				  <times/>
				  <ci> k10 </ci>
				  <ci> RN </ci>
				</apply>
			  </apply>
			  <apply>
				<plus/>
				<apply>
				  <times/>
				  <ci> k9 </ci>
				  <ci> RC </ci>
				</apply>
				<apply>
				  <times/>
				  <ci> vdRC </ci> 
				  <apply>
					<divide/>
					<ci> RC </ci>
					<apply>
					  <plus/>
					  <ci> Kd </ci>
					  <ci> RC </ci>
					</apply>
				  </apply>
				</apply>
				<apply>
				  <times/>
				  <ci> kdn </ci>
				  <ci> RC </ci>
				</apply>
			  </apply>
			</apply>
		  </apply>
		</math>
	</component>
	<component cmeta:id="PCP" name="PCP">
		<variable cmeta:id="PCP_PCP" name="PCP" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="PCP_init" private_interface="out" units="nanomolar"/>
		<variable name="vdPC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/> 
		<variable name="Kd" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kdn" private_interface="out" public_interface="in" 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="PC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="CCP" name="CCP">
		<variable cmeta:id="CCP_CCP" name="CCP" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="CCP_init" private_interface="out" units="nanomolar"/>
		<variable name="vdCC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/> 
		<variable name="Kd" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kdn" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V1C" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V2C" private_interface="out" public_interface="in" units="flux"/>
		<variable name="CC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="PCC" name="PCC">
		<variable cmeta:id="PCC_PCC" initial_value="PCC_init" name="PCC" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="PCC_init" units="nanomolar"/>
		<variable name="Kp" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" public_interface="in" units="nanomolar"/>
		<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="kdn" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V1PC" public_interface="in" units="flux"/>
		<variable name="V2PC" public_interface="in" units="flux"/>
		<variable name="PCCP" public_interface="in" units="nanomolar"/>
		<variable name="PCN" public_interface="in" units="nanomolar"/>
		<variable name="PC" public_interface="in" units="nanomolar"/>
		<variable name="CC" 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> PCC </ci> 
			</apply>
			<apply>
			  <minus/>
			  <apply>
				<plus/>
				<apply>
				  <times/>
				  <ci> V2PC </ci> 
				  <apply>
					<divide/>
					<ci> PCCP </ci>
					<apply>
					  <plus/>
					  <ci> Kdp </ci>
					  <ci> PCCP </ci>
					</apply>
				  </apply>
				</apply>
				<apply>
				  <times/>
				  <ci> k3 </ci>
				  <ci> PC </ci>
				  <ci> CC </ci>
				</apply>
				<apply>
				  <times/>
				  <ci> k2 </ci>
				  <ci> PCN </ci>
				</apply>
			  </apply>
			  <apply>
				<plus/>
				<apply>
				  <times/>
				  <ci> V1PC </ci> 
				  <apply>
					<divide/>
					<ci> PCC </ci>
					<apply>
					  <plus/>
					  <ci> Kp </ci>
					  <ci> PCC </ci>
					</apply>
				  </apply>
				</apply>
				<apply>
				  <times/>
				  <ci> k4 </ci>
				  <ci> PCC </ci>
				</apply>
				<apply>
				  <times/>
				  <ci> k1 </ci>
				  <ci> PCC </ci>
				</apply>
				<apply>
				  <times/>
				  <ci> kdn </ci>
				  <ci> PCC </ci>
				</apply> 
			  </apply>
			</apply>
		  </apply>
		</math>
	</component>
	<component cmeta:id="PCN" name="PCN">
		<variable cmeta:id="PCN_PCN" name="PCN" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="PCN_init" private_interface="out" units="nanomolar"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="k1" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k2" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k7" private_interface="out" public_interface="in" units="second_order_rate_constant"/>
		<variable name="k8" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="kdn" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V3PC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V4PC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="PCNP" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="PCC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="BN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="IN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="RN" name="RN">
		<variable cmeta:id="RN_RN" initial_value="RN_init" name="RN" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="RN_init" units="nanomolar"/>
		<variable name="Kd" public_interface="in" units="nanomolar"/>
		<variable name="k9" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k10" public_interface="in" units="first_order_rate_constant"/>
		<variable name="kdn" public_interface="in" units="first_order_rate_constant"/>
		<variable name="vdRN" public_interface="in" units="flux"/>
		<variable name="RC" public_interface="in" units="nanomolar"/>
		<variable name="time" public_interface="in" units="hour"/>
		<variable initial_value="1" name="one" units="nanomolar"/>
		<math xmlns="http://www.w3.org/1998/Math/MathML"> 
		  <apply>
			<eq/> 
			<apply>
			  <diff/> 
			  <bvar>
				<ci>time</ci>
			  </bvar> 
			  <ci> RN </ci> 
			</apply>
			<apply>
			  <minus/>
			  <apply>
				<times/>
				<ci> k9 </ci>
				<ci> RC </ci>
			  </apply>
			  <apply>
				<plus/>  
				<apply>
				  <times/>
				  <ci> k10 </ci>
				  <ci> RN </ci>
				</apply>
				<apply>
				  <times/>
				  <ci> vdRN</ci> 
				  <apply>
					<divide/>
					<ci> RN </ci>
					<apply>
					  <plus/>
					  <ci> Kd </ci>
					  <ci> RN </ci>
					</apply>
				  </apply>
				</apply>
				<apply>
				  <times/>
				  <ci> kdn </ci>
				  <ci> RN </ci>
				</apply>
			  </apply>
			</apply>
		  </apply>
		</math>
	</component>
	<component cmeta:id="PCCP" name="PCCP">
		<variable cmeta:id="PCCP_PCCP" name="PCCP" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="PCCP_init" private_interface="out" units="nanomolar"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kd" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kdn" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V1PC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V2PC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="vdPCC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="PCC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="PCNP" name="PCNP">
		<variable cmeta:id="PCNP_PCNP" name="PCNP" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="PCNP_init" private_interface="out" units="nanomolar"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kd" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kdn" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V3PC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V4PC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="vdPCN" private_interface="out" public_interface="in" units="flux"/>
		<variable name="PCN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="BC" name="BC">
		<variable cmeta:id="BC_BC" initial_value="BC_init" name="BC" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="BC_init" units="nanomolar"/>
		<variable name="Kp" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" public_interface="in" units="nanomolar"/>
		<variable name="ksB" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k5" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k6" public_interface="in" units="first_order_rate_constant"/>
		<variable name="kdn" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V1B" public_interface="in" units="flux"/>
		<variable name="V2B" public_interface="in" units="flux"/>
		<variable name="MB" public_interface="in" units="nanomolar"/>
		<variable name="BCP" public_interface="in" units="nanomolar"/>
		<variable name="BN" 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>BC</ci>
				</apply>
				<apply>
					<minus/>
					<apply>
						<plus/>
						<apply>
							<times/>
							<ci>V2B</ci>
							<apply>
								<divide/>
								<ci>BCP</ci>
								<apply>
									<plus/>
									<ci>Kdp</ci>
									<ci>BCP</ci>
								</apply>
							</apply>
						</apply>
						<apply>
							<times/>
							<ci>k6</ci>
							<ci>BN</ci>
						</apply>
						<apply>
							<times/>
							<ci>ksB</ci>
							<ci>MB</ci>
						</apply>
					</apply>
					<apply>
						<plus/>
						<apply>
							<times/>
							<ci>V1B</ci>
							<apply>
								<divide/>
								<ci>BC</ci>
								<apply>
									<plus/>
									<ci>Kp</ci>
									<ci>BC</ci>
								</apply>
							</apply>
						</apply>
						<apply>
							<times/>
							<ci>k5</ci>
							<ci>BC</ci>
						</apply>
						<apply>
							<times/>
							<ci>kdn</ci>
							<ci>BC</ci>
						</apply>
					</apply>
				</apply>
			</apply>
		</math>
	</component>
	<component cmeta:id="BCP" name="BCP">
		<variable cmeta:id="BCP_BCP" name="BCP" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="BCP_init" private_interface="out" units="nanomolar"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kd" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kdn" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V1B" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V2B" private_interface="out" public_interface="in" units="flux"/>
		<variable name="vdBC" private_interface="out" public_interface="in" units="flux"/>
		<variable name="BC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="BN" name="BN">
		<variable cmeta:id="BN_BN" name="BN" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="BN_init" private_interface="out" units="nanomolar"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="k5" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k6" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k7" private_interface="out" public_interface="in" units="second_order_rate_constant"/>
		<variable name="k8" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="kdn" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V3B" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V4B" private_interface="out" public_interface="in" units="flux"/>
		<variable name="BNP" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="BC" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="PCN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="IN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="BNP" name="BNP"> 
		<variable cmeta:id="BNP_BNP" name="BNP" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="BNP_init" private_interface="out" units="nanomolar"/>
		<variable name="Kp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kdp" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="Kd" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="kdn" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
		<variable name="V3B" private_interface="out" public_interface="in" units="flux"/>
		<variable name="V4B" private_interface="out" public_interface="in" units="flux"/>
		<variable name="vdBN" private_interface="out" public_interface="in" units="flux"/>
		<variable name="BN" private_interface="out" public_interface="in" units="nanomolar"/>
		<variable name="time" private_interface="out" public_interface="in" units="hour"/>
	</component>
	<component cmeta:id="IN" name="IN"> 
		<variable cmeta:id="IN_IN" initial_value="IN_init" name="IN" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="IN_init" units="nanomolar"/>
		<variable name="Kd" public_interface="in" units="nanomolar"/>
		<variable name="kdn" public_interface="in" units="first_order_rate_constant"/>
		<variable name="k7" public_interface="in" units="second_order_rate_constant"/>
		<variable name="k8" public_interface="in" units="first_order_rate_constant"/>
		<variable name="vdIN" public_interface="in" units="flux"/>
		<variable name="PCN" public_interface="in" units="nanomolar"/>
		<variable name="BN" 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> IN </ci> 
			</apply>
			<apply>
			  <minus/>
			  <apply>
				<times/>
				<ci> k7 </ci>
				<ci> BN </ci>
				<ci> PCN </ci>
			  </apply>
			  <apply>
				<plus/>
				<apply>
				  <times/>
				  <ci> k8 </ci>
				  <ci> IN </ci>
				</apply> 
				<apply>
				  <times/>
				  <ci> vdIN </ci> 
				  <apply>
					<divide/>
					<ci> IN </ci>
					<apply>
					  <plus/>
					  <ci> Kd </ci>
					  <ci> IN </ci>
					</apply>
				  </apply>
				</apply>
				<apply>
				  <times/>
				  <ci> kdn </ci>
				  <ci> IN </ci>
				</apply> 
			  </apply>
			</apply>
		  </apply>
		</math>
	</component>
	<component name="model_parameters"> 
		<variable initial_value="0.8" name="k1" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.4" name="k2" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.8" name="k3" public_interface="out" units="second_order_rate_constant"/>
		<variable initial_value="0.4" name="k4" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.8" name="k5" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.4" name="k6" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="1.0" name="k7" public_interface="out" units="second_order_rate_constant"/>
		<variable initial_value="0.2" name="k8" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.8" name="k9" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.4" name="k10" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.02" name="kdnc" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.02" name="kdn" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.32" name="ksB" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="3.2" name="ksC" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="1.2" name="ksP" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="1.7" name="ksR" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="2.0" name="m" public_interface="out" units="dimensionless"/>
		<variable initial_value="2.0" name="n" public_interface="out" units="dimensionless"/>
		<variable initial_value="2.0" name="h" public_interface="out" units="dimensionless"/>
		<variable initial_value="0.3" name="Kd" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="Kdp" public_interface="out" units="nanomolar"/>
		<variable initial_value="1.006" name="Kp" public_interface="out" units="nanomolar"/>
		<variable initial_value="1.4" name="V1B" public_interface="out" units="flux"/>
		<variable initial_value="1.2" name="V1C" public_interface="out" units="flux"/>
		<variable initial_value="9.6" name="V1P" public_interface="out" units="flux"/>
		<variable initial_value="2.4" name="V1PC" public_interface="out" units="flux"/>
		<variable initial_value="0.2" name="V2B" public_interface="out" units="flux"/>
		<variable initial_value="0.2" name="V2C" public_interface="out" units="flux"/>
		<variable initial_value="0.6" name="V2P" public_interface="out" units="flux"/>
		<variable initial_value="0.2" name="V2PC" public_interface="out" units="flux"/>
		<variable initial_value="1.4" name="V3B" public_interface="out" units="flux"/>
		<variable initial_value="2.4" name="V3PC" public_interface="out" units="flux"/>
		<variable initial_value="0.4" name="V4B" public_interface="out" units="flux"/>
		<variable initial_value="0.2" name="V4PC" public_interface="out" units="flux"/>
		<variable initial_value="3.0" name="vdBC" public_interface="out" units="flux"/>
		<variable initial_value="3.0" name="vdBN" public_interface="out" units="flux"/>
		<variable initial_value="1.4" name="vdCC" public_interface="out" units="flux"/>
		<variable initial_value="1.6" name="vdIN" public_interface="out" units="flux"/>
		<variable initial_value="3.4" name="vdPC" public_interface="out" units="flux"/>
		<variable initial_value="1.4" name="vdPCC" public_interface="out" units="flux"/>
		<variable initial_value="1.4" name="vdPCN" public_interface="out" units="flux"/>
		<variable initial_value="4.4" name="vdRC" public_interface="out" units="flux"/>
		<variable initial_value="0.8" name="vdRN" public_interface="out" units="flux"/>
	</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="MC">
				<component_ref component="M_C"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="MB">
				<component_ref component="M_B"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="MR">
				<component_ref component="M_R"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="PC">
				<component_ref component="P_C"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="CC">
				<component_ref component="C_C"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="PCP">
				<component_ref component="P_CP"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="CCP">
				<component_ref component="C_CP"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="PCN">
				<component_ref component="PC_N"/>
			</component_ref>
	</group>
	
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="PCCP">
				<component_ref component="PC_CP"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="PCNP">
				<component_ref component="PC_NP"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="BCP">
				<component_ref component="B_CP"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="BN">
				<component_ref component="B_N"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="BNP">
				<component_ref component="B_NP"/>
			</component_ref>
	</group>
	
	
	<!-- imports -->
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="mRNAs.cellml">
		<component component_ref="interface" name="M_P"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="mRNAs.cellml">
		<component component_ref="interface" name="M_C"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="mRNAs.cellml">
		<component component_ref="interface" name="M_B"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="mRNAs.cellml">
		<component component_ref="interface" name="M_R"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law1.cellml">
		<component component_ref="interface" name="P_C"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law1.cellml">
		<component component_ref="interface" name="C_C"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law2.cellml">
		<component component_ref="interface" name="P_CP"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law2.cellml">
		<component component_ref="interface" name="C_CP"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law2.cellml">
		<component component_ref="interface" name="PC_CP"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law2.cellml">
		<component component_ref="interface" name="PC_NP"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law2.cellml">
		<component component_ref="interface" name="B_CP"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law2.cellml">
		<component component_ref="interface" name="B_NP"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law3.cellml">
		<component component_ref="interface" name="PC_N"/>
	</import>
	<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="law3.cellml">
		<component component_ref="interface" name="B_N"/>
	</import>
	
	
  
	<!-- connections -->
	<connection>
		<map_components component_1="MP" component_2="environment"/>
		<map_variables variable_1="time" variable_2="time"/>
	</connection>

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

	<connection>
		<map_components component_1="MP" component_2="BN"/>
		<map_variables variable_1="BN" variable_2="BN"/>
	</connection> 

	<connection>
		<map_components component_1="MP" component_2="model_parameters"/>
		<map_variables variable_1="n" variable_2="n"/>
	</connection>

	<connection>
		<map_components component_1="MC" component_2="BN"/>
		<map_variables variable_1="BN" variable_2="BN"/>
	</connection>

	<connection>
		<map_components component_1="MR" component_2="BN"/>
		<map_variables variable_1="BN" variable_2="BN"/>
	</connection>  

	<connection>
		<map_components component_1="MC" component_2="model_parameters"/>
		<map_variables variable_1="n" variable_2="n"/>
	</connection>

	<connection>
		<map_components component_1="MR" component_2="model_parameters"/>
		<map_variables variable_1="h" variable_2="h"/>
	</connection>

	<connection>
		<map_components component_1="MB" component_2="RN"/>
		<map_variables variable_1="RN" variable_2="RN"/>
	</connection> 

	<connection>
		<map_components component_1="MB" component_2="model_parameters"/>
		<map_variables variable_1="m" variable_2="m"/>
	</connection>

	<connection>
		<map_components component_1="PC" component_2="MP"/>
		<map_variables variable_1="MP" variable_2="MP"/>
	</connection> 

	<connection>
		<map_components component_1="PC" component_2="PCP"/>
		<map_variables variable_1="PC" variable_2="PC"/>
		<map_variables variable_1="PCP" variable_2="PCP"/>
	</connection> 

	<connection>
		<map_components component_1="PC" component_2="PCC"/>
		<map_variables variable_1="PC" variable_2="PC"/>
		<map_variables variable_1="PCC" variable_2="PCC"/>
	</connection>

	<connection>
		<map_components component_1="PC" component_2="CC"/>
		<map_variables variable_1="PC" variable_2="PC"/>
		<map_variables variable_1="CC" variable_2="CC"/>
	</connection>  

	<connection>
		<map_components component_1="PC" component_2="model_parameters"/>
		<map_variables variable_1="ksP" variable_2="ksP"/>
		<map_variables variable_1="V1P" variable_2="V1P"/>
		<map_variables variable_1="V2P" variable_2="V2P"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="k3" variable_2="k3"/>
		<map_variables variable_1="k4" variable_2="k4"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="CC" component_2="MC"/>
		<map_variables variable_1="MC" variable_2="MC"/>
	</connection> 

	<connection>
		<map_components component_1="CC" component_2="CCP"/>
		<map_variables variable_1="CC" variable_2="CC"/>
		<map_variables variable_1="CCP" variable_2="CCP"/>
	</connection> 

	<connection>
		<map_components component_1="CC" component_2="PCC"/>
		<map_variables variable_1="CC" variable_2="CC"/>
		<map_variables variable_1="PCC" variable_2="PCC"/>
	</connection>

	<connection>
		<map_components component_1="CC" component_2="model_parameters"/>
		<map_variables variable_1="ksC" variable_2="ksC"/>
		<map_variables variable_1="V1C" variable_2="V1C"/>
		<map_variables variable_1="V2C" variable_2="V2C"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="k3" variable_2="k3"/>
		<map_variables variable_1="k4" variable_2="k4"/>
		<map_variables variable_1="kdnc" variable_2="kdnc"/>
	</connection>

	<connection>
		<map_components component_1="PCP" component_2="model_parameters"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="V1P" variable_2="V1P"/>
		<map_variables variable_1="V2P" variable_2="V2P"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="vdPC" variable_2="vdPC"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="CCP" component_2="model_parameters"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="V1C" variable_2="V1C"/>
		<map_variables variable_1="V2C" variable_2="V2C"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="vdCC" variable_2="vdCC"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="PCC" component_2="PCCP"/>
		<map_variables variable_1="PCC" variable_2="PCC"/>
		<map_variables variable_1="PCCP" variable_2="PCCP"/>
	</connection> 

	<connection>
		<map_components component_1="PCC" component_2="PCN"/>
		<map_variables variable_1="PCC" variable_2="PCC"/>
		<map_variables variable_1="PCN" variable_2="PCN"/>
	</connection>

	<connection>
		<map_components component_1="PCC" component_2="model_parameters"/>
		<map_variables variable_1="k1" variable_2="k1"/>
		<map_variables variable_1="k2" variable_2="k2"/>
		<map_variables variable_1="V1PC" variable_2="V1PC"/>
		<map_variables variable_1="V2PC" variable_2="V2PC"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="k3" variable_2="k3"/>
		<map_variables variable_1="k4" variable_2="k4"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="PCN" component_2="BN"/>
		<map_variables variable_1="PCN" variable_2="PCN"/>
		<map_variables variable_1="BN" variable_2="BN"/>
	</connection>

	<connection>
		<map_components component_1="PCN" component_2="IN"/>
		<map_variables variable_1="PCN" variable_2="PCN"/>
		<map_variables variable_1="IN" variable_2="IN"/>
	</connection>  

	<connection>
		<map_components component_1="PCNP" component_2="PCN"/>
		<map_variables variable_1="PCNP" variable_2="PCNP"/>
		<map_variables variable_1="PCN" variable_2="PCN"/>
	</connection>

	<connection>
		<map_components component_1="PCN" component_2="model_parameters"/>
		<map_variables variable_1="k1" variable_2="k1"/>
		<map_variables variable_1="k2" variable_2="k2"/>
		<map_variables variable_1="V3PC" variable_2="V3PC"/>
		<map_variables variable_1="V4PC" variable_2="V4PC"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="k7" variable_2="k7"/>
		<map_variables variable_1="k8" variable_2="k8"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="PCCP" component_2="model_parameters"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="V1PC" variable_2="V1PC"/>
		<map_variables variable_1="V2PC" variable_2="V2PC"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="vdPCC" variable_2="vdPCC"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="PCNP" component_2="model_parameters"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="V3PC" variable_2="V3PC"/>
		<map_variables variable_1="V4PC" variable_2="V4PC"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="vdPCN" variable_2="vdPCN"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="BC" component_2="MB"/>
		<map_variables variable_1="MB" variable_2="MB"/>
	</connection> 

	<connection>
		<map_components component_1="BC" component_2="BCP"/>
		<map_variables variable_1="BC" variable_2="BC"/>
		<map_variables variable_1="BCP" variable_2="BCP"/>
	</connection> 

	<connection>
		<map_components component_1="BC" component_2="BN"/>
		<map_variables variable_1="BC" variable_2="BC"/>
		<map_variables variable_1="BN" variable_2="BN"/>
	</connection>

	<connection>
		<map_components component_1="BC" component_2="model_parameters"/>
		<map_variables variable_1="ksB" variable_2="ksB"/>
		<map_variables variable_1="V1B" variable_2="V1B"/>
		<map_variables variable_1="V2B" variable_2="V2B"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="k5" variable_2="k5"/>
		<map_variables variable_1="k6" variable_2="k6"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="BCP" component_2="model_parameters"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="V1B" variable_2="V1B"/>
		<map_variables variable_1="V2B" variable_2="V2B"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="vdBC" variable_2="vdBC"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="BN" component_2="BNP"/>
		<map_variables variable_1="BN" variable_2="BN"/>
		<map_variables variable_1="BNP" variable_2="BNP"/>
	</connection> 

	<connection>
		<map_components component_1="IN" component_2="BN"/>
		<map_variables variable_1="IN" variable_2="IN"/>
		<map_variables variable_1="BN" variable_2="BN"/>
	</connection>

	<connection>
		<map_components component_1="BN" component_2="model_parameters"/>
		<map_variables variable_1="k7" variable_2="k7"/>
		<map_variables variable_1="k8" variable_2="k8"/>
		<map_variables variable_1="V3B" variable_2="V3B"/>
		<map_variables variable_1="V4B" variable_2="V4B"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="k5" variable_2="k5"/>
		<map_variables variable_1="k6" variable_2="k6"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="BNP" component_2="model_parameters"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="V3B" variable_2="V3B"/>
		<map_variables variable_1="V4B" variable_2="V4B"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="vdBN" variable_2="vdBN"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="IN" component_2="model_parameters"/>
		<map_variables variable_1="k7" variable_2="k7"/>
		<map_variables variable_1="k8" variable_2="k8"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="vdIN" variable_2="vdIN"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</connection>

	<connection>
		<map_components component_1="RC" component_2="RN"/>
		<map_variables variable_1="RN" variable_2="RN"/>
		<map_variables variable_1="RC" variable_2="RC"/>
	</connection>

	<connection>
		<map_components component_1="RC" component_2="MR"/>
		<map_variables variable_1="MR" variable_2="MR"/>
	</connection> 

	<connection>
		<map_components component_1="RC" component_2="model_parameters"/>
		<map_variables variable_1="k9" variable_2="k9"/>
		<map_variables variable_1="k10" variable_2="k10"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="vdRC" variable_2="vdRC"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
		<map_variables variable_1="ksR" variable_2="ksR"/>
	</connection>

	<connection>
		<map_components component_1="RN" component_2="model_parameters"/>
		<map_variables variable_1="k9" variable_2="k9"/>
		<map_variables variable_1="k10" variable_2="k10"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="vdRN" variable_2="vdRN"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
	</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="kdm" variable_2="kdmp"/>
		<map_variables variable_1="B" variable_2="KAP"/>
		<map_variables variable_1="Km" variable_2="KmP"/>
		<map_variables variable_1="c" variable_2="n"/>
		<map_variables variable_1="A" variable_2="BN"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="M_init" variable_2="MP_init"/>
	</connection>
	<connection>
		<map_components component_1="M_C" component_2="MC"/>
		<map_variables variable_1="M" variable_2="MC"/>
		<map_variables variable_1="vs" variable_2="vsC"/>
		<map_variables variable_1="vm" variable_2="vmC"/>
		<map_variables variable_1="kdm" variable_2="kdmc"/>
		<map_variables variable_1="B" variable_2="KAC"/>
		<map_variables variable_1="Km" variable_2="KmC"/>
		<map_variables variable_1="c" variable_2="n"/>
		<map_variables variable_1="A" variable_2="BN"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="M_init" variable_2="MC_init"/>
	</connection>
	<connection>
		<map_components component_1="M_B" component_2="MB"/>
		<map_variables variable_1="M" variable_2="MB"/>
		<map_variables variable_1="vs" variable_2="vsB"/>
		<map_variables variable_1="vm" variable_2="vmB"/>
		<map_variables variable_1="kdm" variable_2="kdmb"/>
		<map_variables variable_1="B" variable_2="RN"/>
		<map_variables variable_1="Km" variable_2="KmB"/>
		<map_variables variable_1="c" variable_2="m"/>
		<map_variables variable_1="A" variable_2="KIB"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="M_init" variable_2="MB_init"/>
	</connection>
	<connection>
		<map_components component_1="M_R" component_2="MR"/>
		<map_variables variable_1="M" variable_2="MR"/>
		<map_variables variable_1="vs" variable_2="vsR"/>
		<map_variables variable_1="vm" variable_2="vmR"/>
		<map_variables variable_1="kdm" variable_2="kdmr"/>
		<map_variables variable_1="B" variable_2="KAR"/>
		<map_variables variable_1="Km" variable_2="KmR"/>
		<map_variables variable_1="c" variable_2="h"/>
		<map_variables variable_1="A" variable_2="BN"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="M_init" variable_2="MR_init"/>
	</connection>
	<connection>
		<map_components component_1="P_C" component_2="PC"/>
		<map_variables variable_1="YC" variable_2="PC"/>
		<map_variables variable_1="ks" variable_2="ksP"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="ka" variable_2="k4"/>
		<map_variables variable_1="kb" variable_2="k3"/>
		<map_variables variable_1="kd" variable_2="kdn"/>
		<map_variables variable_1="V1" variable_2="V1P"/>
		<map_variables variable_1="V2" variable_2="V2P"/>
		<map_variables variable_1="YCP" variable_2="PCP"/>
		<map_variables variable_1="A" variable_2="PCC"/>
		<map_variables variable_1="M" variable_2="MP"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="B" variable_2="CC"/>
		<map_variables variable_1="YC_init" variable_2="PC_init"/>
	</connection>
	<connection>
		<map_components component_1="C_C" component_2="CC"/>
		<map_variables variable_1="YC" variable_2="CC"/>
		<map_variables variable_1="ks" variable_2="ksC"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="ka" variable_2="k4"/>
		<map_variables variable_1="kb" variable_2="k3"/>
		<map_variables variable_1="kd" variable_2="kdnc"/>
		<map_variables variable_1="V1" variable_2="V1C"/>
		<map_variables variable_1="V2" variable_2="V2C"/>
		<map_variables variable_1="YCP" variable_2="CCP"/>
		<map_variables variable_1="A" variable_2="PCC"/>
		<map_variables variable_1="M" variable_2="MC"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="B" variable_2="PC"/>
		<map_variables variable_1="YC_init" variable_2="CC_init"/>
	</connection>
	<connection>
		<map_components component_1="P_CP" component_2="PCP"/>
		<map_variables variable_1="YCP" variable_2="PCP"/>
		<map_variables variable_1="vd" variable_2="vdPC"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
		<map_variables variable_1="V1" variable_2="V1P"/>
		<map_variables variable_1="V2" variable_2="V2P"/>
		<map_variables variable_1="YC" variable_2="PC"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="YCP_init" variable_2="PCP_init"/>
	</connection>
	<connection>
		<map_components component_1="C_CP" component_2="CCP"/>
		<map_variables variable_1="YCP" variable_2="CCP"/>
		<map_variables variable_1="vd" variable_2="vdCC"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
		<map_variables variable_1="V1" variable_2="V1C"/>
		<map_variables variable_1="V2" variable_2="V2C"/>
		<map_variables variable_1="YC" variable_2="CC"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="YCP_init" variable_2="CCP_init"/>
	</connection>
	<connection>
		<map_components component_1="PC_CP" component_2="PCCP"/>
		<map_variables variable_1="YCP" variable_2="PCCP"/>
		<map_variables variable_1="vd" variable_2="vdPCC"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
		<map_variables variable_1="V1" variable_2="V1PC"/>
		<map_variables variable_1="V2" variable_2="V2PC"/>
		<map_variables variable_1="YC" variable_2="PCC"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="YCP_init" variable_2="PCCP_init"/>
	</connection>
	<connection>
		<map_components component_1="PC_NP" component_2="PCNP"/>
		<map_variables variable_1="YCP" variable_2="PCNP"/>
		<map_variables variable_1="vd" variable_2="vdPCN"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
		<map_variables variable_1="V1" variable_2="V3PC"/>
		<map_variables variable_1="V2" variable_2="V4PC"/>
		<map_variables variable_1="YC" variable_2="PCN"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="YCP_init" variable_2="PCNP_init"/>
	</connection>
	<connection>
		<map_components component_1="B_CP" component_2="BCP"/>
		<map_variables variable_1="YCP" variable_2="BCP"/>
		<map_variables variable_1="vd" variable_2="vdBC"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
		<map_variables variable_1="V1" variable_2="V1B"/>
		<map_variables variable_1="V2" variable_2="V2B"/>
		<map_variables variable_1="YC" variable_2="BC"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="YCP_init" variable_2="BCP_init"/>
	</connection>
	<connection>
		<map_components component_1="B_NP" component_2="BNP"/>
		<map_variables variable_1="YCP" variable_2="BNP"/>
		<map_variables variable_1="vd" variable_2="vdBN"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="Kd" variable_2="Kd"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
		<map_variables variable_1="V1" variable_2="V3B"/>
		<map_variables variable_1="V2" variable_2="V4B"/>
		<map_variables variable_1="YC" variable_2="BN"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="YCP_init" variable_2="BNP_init"/>
	</connection>
	<connection>
		<map_components component_1="PC_N" component_2="PCN"/>
		<map_variables variable_1="Y" variable_2="PCN"/>
		<map_variables variable_1="Kp" variable_2="Kp"/>
		<map_variables variable_1="Kdp" variable_2="Kdp"/>
		<map_variables variable_1="kc" variable_2="k7"/>
		<map_variables variable_1="ke" variable_2="k8"/>
		<map_variables variable_1="kb" variable_2="k1"/>
		<map_variables variable_1="ka" variable_2="k2"/>
		<map_variables variable_1="kdn" variable_2="kdn"/>
		<map_variables variable_1="V1" variable_2="V3PC"/>
		<map_variables variable_1="V2" variable_2="V4PC"/>
		<map_variables variable_1="YP" variable_2="PCNP"/>
		<map_variables variable_1="E" variable_2="IN"/>
		<map_variables variable_1="A" variable_2="PCC"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="C" variable_2="BN"/>
		<map_variables variable_1="Y_init" variable_2="PCN_init"/>
	</connection>
	
	
	
	<connection>
			<map_components component_1="B_N" component_2="BN"/>
			<map_variables variable_1="Y" variable_2="BN"/>
			<map_variables variable_1="Kp" variable_2="Kp"/>
			<map_variables variable_1="Kdp" variable_2="Kdp"/>
			<map_variables variable_1="kc" variable_2="k7"/>
			<map_variables variable_1="ke" variable_2="k8"/>
			<map_variables variable_1="kb" variable_2="k5"/>
			<map_variables variable_1="ka" variable_2="k6"/>
			<map_variables variable_1="kdn" variable_2="kdn"/>
			<map_variables variable_1="V1" variable_2="V3B"/>
			<map_variables variable_1="V2" variable_2="V4B"/>
			<map_variables variable_1="YP" variable_2="BNP"/>
			<map_variables variable_1="E" variable_2="IN"/>
			<map_variables variable_1="A" variable_2="BC"/>
			<map_variables variable_1="time" variable_2="time"/>
			<map_variables variable_1="C" variable_2="PCN"/>
			<map_variables variable_1="Y_init" variable_2="BN_init"/>
	</connection>
   
   	<!-- metadata -->
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="leloup_goldbeter_2003_1.1model">
			<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>18-01-2011</dcterms:W3CDTF>
			</dcterms:created>
		</rdf:Description>
		<rdf:Description rdf:about="#leloup_goldbeter_2003_1.1model">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Toward a detailed computational model for the mammalian circadian clock</dc:title>
			<cmeta:comment rdf:parseType="Resource">
				<rdf:value>This is a CellML description of Leloup and Goldbeter's 2004 mathematical model of the mammalian circadian clock.
				</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>
			<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>Bmal1 mRNA</rdf:li>
							<rdf:li>signal transduction</rdf:li>
							<rdf:li>circadian rhythms</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>12775757</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>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/">Toward a detailed computational model for the mammalian circadian clock</dc:title>
					<bqs:volume>100</bqs:volume>
					<bqs:first_page>7051</bqs:first_page>
					<bqs:last_page>7056</bqs:last_page>
				</bqs:JournalArticle>
			</bqs:reference>
		</rdf:Description>
	</rdf:RDF>
</model>