Location: Leloup, Goldbeter, 2004 @ e12c13b3c70b / leloup_2004_1.1model.cellml

Author:
Dougal Cowan <devnull@localhost>
Date:
2011-11-04 09:46:20+13:00
Desc:
Updated images to fix the swapped K7/K8 arrows, noticed by Arne Gjuvsland.
Permanent Source URI:
http://models.cellml.org/workspace/leloup_goldbeter_2004/rawfile/e12c13b3c70b02aeb46795128b58f48326c20781/leloup_2004_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_2004_1.1model" name="leloup_2004_1.1model">
	<!-- documentation -->
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
		  <articleinfo>
			  <title>Modeling the mammalian circadian clock: sensitivity analysis and multiplicity of oscillatory mechanisms</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 2C from the published paper). The model needs to be run for 160 hours with a step size of 0.1 hours, and MP, MB and MC are plotted. 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. Ignore the first couple of oscillations as the model needs to stabilise. Also please note that in the version of the model presented here, parameter set 4 from the original paper has been used. The units have been checked and they are consistent.This file is known to run in OpenCell only.
					</para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
					<para>ABSTRACT: We extend the study of a computational model recently proposed for the mammalian circadian clock (Proc. Natl Acad. Sci. USA 100 (2003) 7051). The model, based on the intertwined positive and negative regulatory loops involving the Per, Cry, Bmal1, and Clock genes, can give rise to sustained circadian oscillations in conditions of continuous darkness. These limit cycle oscillations correspond to circadian rhythms autonomously generated by suprachiasmatic nuclei and by some peripheral tissues. By using different sets of parameter values producing circadian oscillations, we compare the effect of the various parameters and show that both the occurrence and the period of the oscillations are generally most sensitive to parameters related to synthesis or degradation of Bmal1 mRNA and BMAL1 protein. The mechanism of circadian oscillations relies on the formation of an inactive complex between PER and CRY and the activators CLOCK and BMAL1 that enhance Per and Cry expression. Bifurcation diagrams and computer simulations nevertheless indicate the possible existence of a second source of oscillatory behavior. Thus, sustained oscillations might arise from the sole negative autoregulation of Bmal1 expression. This second oscillatory mechanism may not be functional in physiological conditions, and its period need not necessarily be circadian. When incorporating the light-induced expression of the Per gene, the model accounts for entrainment of the oscillations by light-dark (LD) cycles. Long-term suppression of circadian oscillations by a single light pulse can occur in the model when a stable steady state coexists with a stable limit cycle. The phase of the oscillations upon entrainment in LD critically depends on the parameters that govern the level of CRY protein. Small changes in the parameters governing CRY levels can shift the peak in Per mRNA from the L to the D phase, or can prevent entrainment. The results are discussed in relation to physiological disorders of the sleep-wake cycle linked to perturbations of the human circadian clock, such as the familial advanced sleep phase syndrome or the non-24h sleep-wake syndrome.
					</para>
				<informalfigure float="0" id="fig_reaction_diagram">
					<mediaobject>
					  <imageobject>
						<objectinfo>
						  <title>model diagram</title>
						</objectinfo>
						<imagedata fileref="leloup_2004.png"/>
					  </imageobject>
					</mediaobject>
					<caption>Schematic diagram of the model for circadian oscillations in mammals including positive and negative feedback on the transcription of the Per, Cry and Bmal1 genes by their protein products.  Per, Cry and Bmal1 are transcribed in the nucleus and are then transferred to the cytosol where their proteins are translated and the mRNAs are degraded.  Light increases the rate of Per transcription.</caption>
				</informalfigure>
				<para>The original paper reference is cited below:</para>
				<para>Modeling the mammalian circadian clock: sensitivity analysis and multiplicity of oscillatory mechanisms, Jean-Christophe Leloup and Albert Goldbeter, 2004, 
				<emphasis>Journal of Theoretical Biology</emphasis>, 230, 541-562. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/15363675">PubMed ID: 15363675</ulink>
				</para>
			</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 name="vsP" private_interface="out" units="flux"/>
		<variable initial_value="1.1" name="vmP" private_interface="out" units="flux"/>
		<variable initial_value="0.01" name="kdmp" private_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.7" name="KAP" private_interface="out" units="nanomolar"/>
		<variable initial_value="0.3" name="KmP" private_interface="out" units="nanomolar"/>
		<variable name="vstot" public_interface="in" units="flux"/>
		<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"/>
	    <variable initial_value="0.5" name="MP_init" private_interface="out" units="nanomolar"/>
		<math xmlns="http://www.w3.org/1998/Math/MathML"> 
		  <apply>
			<eq/> 
			<ci> vsP </ci> 
			<ci> vstot </ci>
		  </apply>
		</math>
	</component>
	<component cmeta:id="MC" name="MC">
		<variable cmeta:id="MC_MC" name="MC" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable name="vsC" private_interface="out" units="flux"/>
		<variable initial_value="1.0" name="vmC" private_interface="out" units="flux"/>
		<variable initial_value="0.01" name="kdmc" private_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="1.0" name="KAC" private_interface="out" units="nanomolar"/>
		<variable initial_value="0.4" name="KmC" private_interface="out" units="nanomolar"/>
		<variable name="vstot" public_interface="in" units="flux"/>
		<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"/>
	    <variable initial_value="0.3" name="MC_init" private_interface="out" units="nanomolar"/>
		<math xmlns="http://www.w3.org/1998/Math/MathML"> 
		  <apply>
			<eq/> 
			<ci> vsC </ci> 
			<apply>
			  <times/>
			  <cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless"> 0.8 </cn>
			  <ci> vstot </ci>
			</apply>
		  </apply>
		</math>
	</component>
	<component cmeta:id="MB" name="MB">
		<variable cmeta:id="MB_MB" name="MB" private_interface="in" public_interface="out" units="nanomolar"/>
		<variable name="vsB" private_interface="out" units="flux"/>
		<variable initial_value="0.2" name="vmB" private_interface="out" units="flux"/>
		<variable initial_value="0.01" name="kdmb" private_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.8" name="KIB" private_interface="out" units="nanomolar"/>
		<variable name="vstot" public_interface="in" units="flux"/>
		<variable name="m" 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"/>
		<variable initial_value="3.1" name="MB_init" private_interface="out" units="nanomolar"/>
		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply>
				<eq/>
				<ci>vsB</ci>
				<apply>
					<times/>
					<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless">0.7</cn>
					<ci>vstot</ci>
				</apply>
			</apply>
		</math>
	<variable initial_value="0.4" name="KmB" private_interface="out" units="nanomolar"/></component>
	<component cmeta:id="PC" name="PC">
		<variable cmeta:id="PC_PC" name="PC" private_interface="in" public_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"/>
		<variable initial_value="0.1" name="PC_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="CC" name="CC">
		<variable cmeta:id="CC_CC" name="CC" private_interface="in" public_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"/>
		<variable initial_value="0.1" name="CC_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="PCP" name="PCP">
		<variable cmeta:id="PCP_PCP" name="PCP" private_interface="in" public_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"/>
		<variable initial_value="0.1" name="PCP_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="CCP" name="CCP">
		<variable cmeta:id="CCP_CCP" name="CCP" private_interface="in" public_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"/>
		<variable initial_value="0.1" name="CCP_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="PCC" name="PCC">
    <variable cmeta:id="PCC_PCC" initial_value="PCC_init" name="PCC" public_interface="out" 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"/>
	<variable initial_value="0.1" name="PCC_init" units="nanomolar"/>
	<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 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"/>
		<variable initial_value="0.1" name="PCN_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="PCCP" name="PCCP">
		<variable cmeta:id="PCCP_PCCP" name="PCCP" private_interface="in" public_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"/>
		<variable initial_value="0.1" name="PCCP_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="PCNP" name="PCNP">
		<variable cmeta:id="PCNP_PCNP" name="PCNP" private_interface="in" public_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"/>
		<variable initial_value="0.1" name="PCNP_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="BC" name="BC">
		<variable cmeta:id="BC_BC" initial_value="BC_init" name="BC" public_interface="out" 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"/>
		<variable initial_value="0.1" name="BC_init" units="nanomolar"/>
		<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 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"/>
		<variable initial_value="0.1" name="BCP_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="BN" name="BN">
		<variable cmeta:id="BN_BN" name="BN" private_interface="in" public_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"/>
		<variable initial_value="0.1" name="BN_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="BNP" name="BNP">
		<variable cmeta:id="BNP_BNP" name="BNP" private_interface="in" public_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"/>
		<variable initial_value="0.1" name="BNP_init" private_interface="out" units="nanomolar"/>
	</component>
	<component cmeta:id="IN" name="IN">
		<variable cmeta:id="IN_IN" initial_value="IN_init" name="IN" public_interface="out" 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"/>
		<variable initial_value="0.1" name="IN_init" units="nanomolar"/><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.2" 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.2" name="k4" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.4" name="k5" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.2" name="k6" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.5" name="k7" public_interface="out" units="second_order_rate_constant"/>
		<variable initial_value="0.1" name="k8" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.01" name="kdnc" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="0.01" name="kdn" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="1.0" name="kstot" units="first_order_rate_constant"/>
		<variable name="ksB" public_interface="out" units="first_order_rate_constant"/>
		<variable name="ksC" public_interface="out" units="first_order_rate_constant"/>
		<variable name="ksP" public_interface="out" units="first_order_rate_constant"/>
		<variable initial_value="4.0" name="m" public_interface="out" units="dimensionless"/>
		<variable initial_value="4.0" name="n" public_interface="out" units="dimensionless"/>
		<variable initial_value="1.0" name="vstot" public_interface="out" units="flux"/>
		<variable initial_value="0.3" name="Kd" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.3" name="Kdp" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.1" name="Kp" public_interface="out" units="nanomolar"/>
		<variable initial_value="0.6" name="Vphos" units="flux"/>
		<variable initial_value="1.0" name="V1B" public_interface="out" units="flux"/>
		<variable initial_value="0.6" name="V1C" public_interface="out" units="flux"/>
		<variable name="V1P" public_interface="out" units="flux"/>
		<variable name="V1PC" public_interface="out" units="flux"/>
		<variable initial_value="0.1" name="V2B" public_interface="out" units="flux"/>
		<variable initial_value="0.1" name="V2C" public_interface="out" units="flux"/>
		<variable initial_value="0.3" name="V2P" public_interface="out" units="flux"/>
		<variable initial_value="0.1" name="V2PC" public_interface="out" units="flux"/>
		<variable initial_value="1.0" name="V3B" public_interface="out" units="flux"/>
		<variable name="V3PC" public_interface="out" units="flux"/>
		<variable initial_value="0.2" name="V4B" public_interface="out" units="flux"/>
		<variable initial_value="0.1" name="V4PC" public_interface="out" units="flux"/>
		<variable initial_value="1.0" name="vdBC" public_interface="out" units="flux"/>
		<variable initial_value="0.5" name="vdBN" public_interface="out" units="flux"/>
		<variable initial_value="0.7" name="vdCC" public_interface="out" units="flux"/>
		<variable initial_value="0.8" name="vdIN" public_interface="out" units="flux"/>
		<variable initial_value="0.7" name="vdPC" public_interface="out" units="flux"/>
		<variable initial_value="1.0" name="vdPCC" public_interface="out" units="flux"/>
		<variable initial_value="1.0" name="vdPCN" public_interface="out" units="flux"/>
		<math xmlns="http://www.w3.org/1998/Math/MathML"> 
		  <apply>
			<eq/> 
			<ci> ksB </ci> 
			<ci> kstot </ci>
		  </apply>
		  <apply>
			<eq/> 
			<ci> ksC </ci> 
			<ci> kstot </ci>
		  </apply>
		  <apply>
			<eq/> 
			<ci> ksP </ci> 
			<apply>
			  <times/>
			  <cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless"> 0.5 </cn>
			  <ci> kstot </ci>
			</apply>
		  </apply>
		  <apply>
			<eq/> 
			<ci> V1P </ci> 
			<ci> Vphos </ci>
		  </apply>
		  <apply>
			<eq/> 
			<ci> V1PC </ci> 
			<ci> Vphos </ci>
		  </apply>
		  <apply>
			<eq/> 
			<ci> V3PC </ci> 
			<ci> Vphos </ci>
		  </apply>
		</math>
	</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="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="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="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="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"/>
		<map_variables variable_1="vstot" variable_2="vstot"/>
	</connection>
	<connection>
		<map_components component_1="MC" 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"/>
		<map_variables variable_1="vstot" variable_2="vstot"/>
	</connection>
	<connection>
		<map_components component_1="MB" component_2="BN"/>
		<map_variables variable_1="BN" variable_2="BN"/>
	</connection> 
	<connection>
		<map_components component_1="MB" component_2="model_parameters"/>
		<map_variables variable_1="m" variable_2="m"/>
		<map_variables variable_1="vstot" variable_2="vstot"/>
	</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="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="BN"/>
		<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="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_2004_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>29-12-2010</dcterms:W3CDTF>
			</dcterms:created>
		</rdf:Description>
		<rdf:Description rdf:about="#leloup_2004_1.1model">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Modeling the mammalian circadian clock: sensitivity analysis and multiplicity of oscillatory mechanisms</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>IN</rdf:li>
							<rdf:li>nuclear inactive complex between PER-CRY and CLOCK-BMAL1</rdf:li>
							<rdf:li>PCCP</rdf:li>
							<rdf:li>cytosolic phosphorylated PER-CRY protein complex</rdf:li>
							<rdf:li>BCP</rdf:li>
							<rdf:li>cytosolic phosphorylated BMAL1 protein</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>15363675</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>Goldbeter</vCard:Family>
										<vCard:Given>Albert</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
						</rdf:Seq>
					</dc:creator>
					<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Modeling the mammalian circadian clock: sensitivity analysis and multiplicity of oscillatory mechanisms</dc:title>
					<bqs:Journal rdf:parseType="Resource">
						<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Journal of Theoretical Biology</dc:title>
					</bqs:Journal>
					<bqs:volume>230</bqs:volume>
					<bqs:first_page>541</bqs:first_page>
					<bqs:last_page>562</bqs:last_page>
				</bqs:JournalArticle>
			</bqs:reference>
		</rdf:Description>
	</rdf:RDF>
</model>