Rendering of the source text

<?xml version="1.0"?>
<!--
This CellML file was generated on 23/03/2010 at 3:47:27 at p.m. using:

COR (0.9.31.1207)
Copyright 2002-2010 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk

CellML 1.0 was used to generate this model
http://www.CellML.org/
-->
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="potter_model_2005" name="potter_model_2005">

	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Response to continuous and pulsatile PTH dosing: a mathematical model for parathyroid hormone receptor kinetics</title>
				<author>
					<firstname>Catherine</firstname>
					<surname>Lloyd</surname>
					<affiliation>
						<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
        This CellML version of the model has been checked in COR and OpenCell.  The units have been checked and are correct and are consistent.  The model runs in COR and Opencell to it recreate the published results. This version recreates pulsatile PTH dosing at high doses.   
				</para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
ABSTRACT: In this paper, we propose a mathematical model for parathyroid hormone receptor (PTH1R) kinetics, focusing on the receptor's response to PTH dosing to discern bone formation responses from bone resorption. The PTH1R is a major target for new osteoporosis treatments, as pulsatile PTH dosing has been shown to induce net bone formation in both animals and humans, and PTH(1-34) was recently FDA approved for the treatment of post-menopausal osteoporosis. PTH has also been shown to cause net bone loss when given continuously, so that the net action of PTH on bone is dependent on the dosing pattern. We have developed a simplified two-state receptor kinetics model for the PTH1R, based on the concepts of Segel et al., to distinguish the activity of active and inactive receptor and receptor-ligand complexes. The goal is to develop a plausible model of the minimal essential biological relationships necessary for understanding the responses to PTH dosing. A two-state model is able to effectively discriminate between continuous and pulsatile PTH dosing using the active species as surrogates for the downstream anabolic response. For continuous PTH dosing, the model predicts a desensitized system dominated by the inactive receptor and complex, consistent with downstream net bone loss that has been demonstrated experimentally. Using pulsatile PTH dosing, the model system predicts a highly sensitized state dominated by the active receptor and complex, corresponding to net bone formation. These results are consistent with the hypothesis that the kinetics of the receptor plays a critical role in the downstream effects of PTH dosing. Moreover, these results indicate that within a range of biologically relevant PTH doses, the two-state model is able to capture the differential behavior of the system for both continuous and pulsatile PTH dosing. The development of such a model provides a rational basis for developing more biologically extensive models that may support the design of optimal dosing strategies for PTH-based anti-osteoporosis treatments. Moreover, this model provides a unique starting point from which to design experiments investigating PTH receptor biology.
</para>
				<informalfigure float="0" id="fig_reaction_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>model diagram</title>
							</objectinfo>
							<imagedata fileref="potter_2005.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the two-state model for PTHR1 binding kinetics  PTH is secreted and/or dosed at a rate D and binds to the active (Ra) and inactive (Ri) forms of the receptor.  PTH-PTHR1 binding results in the formation of the active (Ca) and inactive (Ci) complexes.  In addition, there is a conversion between the active and inactive forms of the receptor and complexes, while unbound PTH is cleared from the system.</caption>
				</informalfigure>
				<para>
The original paper reference is cited below:
</para>
				<para>
Response to continuous and pulsatile PTH dosing: a mathematical model for parathyroid hormone receptor kinetics, Laura K. Potter, Larry D. Greller, Carolyn R. Cho, Mark E. Nuttall, George B. Stroup, Larry J. Suva, and Frank L. Tobin, 2005, <emphasis>Bone</emphasis>, 37, 159-169.  <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/15921971">PubMed ID: 15921971</ulink>
				</para>
				
			</sect1>
		</article>
	</documentation>
	
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#P">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">P</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">PTH</dcterms:alternative>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">parathyroid hormone</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#Ra">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Ra</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">active PTH1R</dcterms:alternative>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">active parathyroid hormone Type 1 receptor</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#Ri">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Ri</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">inactive PTH1R</dcterms:alternative>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">inactive parathyroid hormone Type 1 receptor</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#Ca">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Ca</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">active PTH-PTH1R complex</dcterms:alternative>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">active parathyroid hormone-parathyroid hormone Type 1 receptor</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#Ci">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Ci</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">inactive PTH-PTH1R complex</dcterms:alternative>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">inactive parathyroid hormone-parathyroid hormone Type 1 receptor</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#rho">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">rho</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">fraction of PTH1R in the active state</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Bag rdf:about="rdf:#4abfc8b0-f7d7-4525-9d28-b4cd292af49e">
			<rdf:li>parathyroid hormone</rdf:li>
			<rdf:li>parathyroid hormone receptor</rdf:li>
			<rdf:li>osteoporosis</rdf:li>
			<rdf:li>endocrine</rdf:li>
		</rdf:Bag>
		<rdf:Seq rdf:about="rdf:#aaa7b2e4-c9b6-42f8-8d83-4322ba940a07">
			<rdf:li rdf:resource="rdf:#2ea39e72-bd62-43e2-8f65-ee80e14c856b"/>
			<rdf:li rdf:resource="rdf:#f815bd64-c3ed-4287-8f50-c00e79f135af"/>
			<rdf:li rdf:resource="rdf:#59b34f3d-87d4-423b-86f7-d45eb5eeede5"/>
			<rdf:li rdf:resource="rdf:#1bb493a0-f279-4fbf-b9f0-5aac09d85626"/>
			<rdf:li rdf:resource="rdf:#956664a5-c801-4abf-bafb-67e4e78603ce"/>
			<rdf:li rdf:resource="rdf:#9bd771f2-8fbc-4305-ae14-53b1a0ffe316"/>
			<rdf:li rdf:resource="rdf:#774a71e8-3e1f-4d58-901a-234d5dffaee7"/>
		</rdf:Seq>
		<rdf:Description rdf:about="rdf:#1a66fc7a-ea7d-45e5-a24b-6452414ea60f">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Bone</dc:title>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#945fc318-18f9-4209-b808-8a7b2b6e25c9">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#875dfb28-dee8-4a20-82ae-c2c19740c37e"/>
			<rdf:value>This is a CellML description of Potter et al's 2005 mathematical model for parathyroid hormone receptor kinetics.</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="">
			<dc:title>This version recreates pulsatile PTH dosing at high doses</dc:title>
			<dc:publisher xmlns:dc="http://purl.org/dc/elements/1.1/">The University of Auckland, Auckland Bioengineering Institute</dc:publisher>
			<cmeta:comment rdf:resource="rdf:#6f5f9832-e954-4a71-82d8-c17188915067"/>
			<dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#807d61c7-7c86-47e0-89ef-c50674ed13d4"/>
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#27d22d11-068b-469f-8164-0e1252689b98"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#926284e1-da44-4000-b28b-a466e1dfbcb4">
			<dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#b565da34-bb98-4478-9c71-e16e5cfb58cc"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#774a71e8-3e1f-4d58-901a-234d5dffaee7">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#01cda416-ab54-4662-ae4e-fab70ffd3459"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#807d61c7-7c86-47e0-89ef-c50674ed13d4">
			<dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-07-23T00:00:00+00:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#d6c78163-a5e0-4bcc-8157-bcbff1f7e2df">
			<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
			<rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#35cfdd97-f2b5-4228-bdf8-be1a9eafbf5b">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Larry</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Suva</vCard:Family>
			<vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">J</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#875dfb28-dee8-4a20-82ae-c2c19740c37e">
			<vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine Lloyd</vCard:FN>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#14a133d2-18cf-4af9-a1dd-41e79c77b58e">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Larry</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Greller</vCard:Family>
			<vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">D</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="#potter_model_2005">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
        Potter et al's 2005 mathematical model for parathyroid hormone receptor kinetics.
      </dc:title>
			<cmeta:comment rdf:resource="rdf:#945fc318-18f9-4209-b808-8a7b2b6e25c9"/>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#22d69f71-1e94-4a10-a69a-a1df691d5fd5"/>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#926284e1-da44-4000-b28b-a466e1dfbcb4"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#947679f7-466d-43a8-94dd-65f29d4043e0">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#aaa7b2e4-c9b6-42f8-8d83-4322ba940a07"/>
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Response to continuous and pulsatile PTH dosing: a mathematical model for parathyroid hormone receptor kinetics</dc:title>
			<bqs:volume xmlns:bqs="http://www.cellml.org/bqs/1.0#">37</bqs:volume>
			<bqs:first_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">159</bqs:first_page>
			<bqs:Journal xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#1a66fc7a-ea7d-45e5-a24b-6452414ea60f"/>
			<dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#1f85f4be-dff7-4bf8-9858-993f8b0160f7"/>
			<bqs:last_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">169</bqs:last_page>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#b53fb134-71a9-4b1b-927c-64faaf958888">
			<vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1bb493a0-f279-4fbf-b9f0-5aac09d85626">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#2b495d83-8580-42a1-b0d9-4116a78b7ae8"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1f85f4be-dff7-4bf8-9858-993f8b0160f7">
			<dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2005-08-00 00:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#01cda416-ab54-4662-ae4e-fab70ffd3459">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Frank</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Tobin</vCard:Family>
			<vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">L</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#2ea39e72-bd62-43e2-8f65-ee80e14c856b">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#ade904da-2598-4ac1-acc4-1aa145efc829"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#ade904da-2598-4ac1-acc4-1aa145efc829">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Laura</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Potter</vCard:Family>
			<vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">K</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#4e6e4b06-9ef5-4d82-a044-56d196e0fd1a">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">George</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Stroup</vCard:Family>
			<vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">B</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#22d69f71-1e94-4a10-a69a-a1df691d5fd5">
			<bqs:Pubmed_id xmlns:bqs="http://www.cellml.org/bqs/1.0#">15921971</bqs:Pubmed_id>
			<bqs:JournalArticle xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#947679f7-466d-43a8-94dd-65f29d4043e0"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#fba1baec-7068-4233-b69f-1f5bdba6fa7f">
			<vCard:Orgname xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">The University of Auckland</vCard:Orgname>
			<vCard:Orgunit xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Auckland Bioengineering Institute</vCard:Orgunit>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#b565da34-bb98-4478-9c71-e16e5cfb58cc">
			<bqs:subject_type xmlns:bqs="http://www.cellml.org/bqs/1.0#">keyword</bqs:subject_type>
			<rdf:value rdf:resource="rdf:#4abfc8b0-f7d7-4525-9d28-b4cd292af49e"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#27d22d11-068b-469f-8164-0e1252689b98">
			<vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#fba1baec-7068-4233-b69f-1f5bdba6fa7f"/>
			<vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#d6c78163-a5e0-4bcc-8157-bcbff1f7e2df"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#cb392c88-80a9-4ac7-9e73-23a75b7e3cc3"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#f815bd64-c3ed-4287-8f50-c00e79f135af">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#14a133d2-18cf-4af9-a1dd-41e79c77b58e"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#59b34f3d-87d4-423b-86f7-d45eb5eeede5">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#bb96d1cc-7202-4f92-9e41-20fb1532f3bc"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#6f5f9832-e954-4a71-82d8-c17188915067">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#b53fb134-71a9-4b1b-927c-64faaf958888"/>
			<rdf:value/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#2b495d83-8580-42a1-b0d9-4116a78b7ae8">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Mark</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Nuttall</vCard:Family>
			<vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">E</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#cb392c88-80a9-4ac7-9e73-23a75b7e3cc3">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lloyd</vCard:Family>
			<vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">May</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#bb96d1cc-7202-4f92-9e41-20fb1532f3bc">
			<vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Carolyn</vCard:Given>
			<vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Cho</vCard:Family>
			<vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">R</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#956664a5-c801-4abf-bafb-67e4e78603ce">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#4e6e4b06-9ef5-4d82-a044-56d196e0fd1a"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#9bd771f2-8fbc-4305-ae14-53b1a0ffe316">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#35cfdd97-f2b5-4228-bdf8-be1a9eafbf5b"/>
		</rdf:Description>
	</rdf:RDF>   
	
   <units name="hour">
      <unit multiplier="3600" units="second"/>
   </units>
   <units name="picomolar">
      <unit prefix="pico" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="flux">
      <unit units="picomolar"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="first_order_rate_constant">
      <unit exponent="-1" units="second"/>
   </units>
   <units name="second_order_rate_constant">
      <unit exponent="-1" units="picomolar"/>
      <unit exponent="-1" units="second"/>
   </units>
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="second"/>
   </component>
   <component name="P">
      <variable initial_value="3" name="P" public_interface="out" units="picomolar"/>
      <variable name="k1" public_interface="in" units="second_order_rate_constant"/>
      <variable name="k1_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k2" public_interface="in" units="second_order_rate_constant"/>
      <variable name="k2_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="kcl" public_interface="in" units="first_order_rate_constant"/>
      <variable name="D" public_interface="in" units="flux"/>
      <variable name="Ca" public_interface="in" units="picomolar"/>
      <variable name="Ci" public_interface="in" units="picomolar"/>
      <variable name="Ra" public_interface="in" units="picomolar"/>
      <variable name="Ri" public_interface="in" units="picomolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k1_</ci>
                     <ci>Ca</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k2_</ci>
                     <ci>Ci</ci>
                  </apply>
                  <ci>D</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k1</ci>
                     <ci>Ra</ci>
                     <ci>P</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k2</ci>
                     <ci>Ri</ci>
                     <ci>P</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>kcl</ci>
                     <ci>P</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ra">
      <variable cmeta:id="Ra_Ra" initial_value="16.9" name="Ra" public_interface="out" units="picomolar"/>
      <variable name="k1" public_interface="in" units="second_order_rate_constant"/>
      <variable name="k1_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k3" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k3_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="Ca" public_interface="in" units="picomolar"/>
      <variable name="P" public_interface="in" units="picomolar"/>
      <variable name="Ri" public_interface="in" units="picomolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ra</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k1_</ci>
                     <ci>Ca</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k3</ci>
                     <ci>Ri</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k1</ci>
                     <ci>Ra</ci>
                     <ci>P</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k3_</ci>
                     <ci>Ra</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ri">
      <variable cmeta:id="Ri_Ri" initial_value="1.7" name="Ri" public_interface="out" units="picomolar"/>
      <variable name="k2" public_interface="in" units="second_order_rate_constant"/>
      <variable name="k2_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k3" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k3_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="Ci" public_interface="in" units="picomolar"/>
      <variable name="P" public_interface="in" units="picomolar"/>
      <variable name="Ra" public_interface="in" units="picomolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ri</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k2_</ci>
                     <ci>Ci</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k3_</ci>
                     <ci>Ra</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k2</ci>
                     <ci>Ri</ci>
                     <ci>P</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k3</ci>
                     <ci>Ri</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca">
      <variable cmeta:id="Ca_Ca" initial_value="4e-4" name="Ca" public_interface="out" units="picomolar"/>
      <variable name="k1" public_interface="in" units="second_order_rate_constant"/>
      <variable name="k1_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k4" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k4_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="Ci" public_interface="in" units="picomolar"/>
      <variable name="P" public_interface="in" units="picomolar"/>
      <variable name="Ra" public_interface="in" units="picomolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k1</ci>
                     <ci>Ra</ci>
                     <ci>P</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k4</ci>
                     <ci>Ci</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k1_</ci>
                     <ci>Ca</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k4_</ci>
                     <ci>Ca</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ci">
      <variable cmeta:id="Ci_Ci" initial_value="0.05" name="Ci" public_interface="out" units="picomolar"/>
      <variable name="k2" public_interface="in" units="second_order_rate_constant"/>
      <variable name="k2_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k4" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k4_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="Ca" public_interface="in" units="picomolar"/>
      <variable name="P" public_interface="in" units="picomolar"/>
      <variable name="Ri" public_interface="in" units="picomolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ci</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k2</ci>
                     <ci>Ri</ci>
                     <ci>P</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k4_</ci>
                     <ci>Ca</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k2_</ci>
                     <ci>Ci</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k4</ci>
                     <ci>Ci</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rho">
      <variable name="rho" units="dimensionless"/>
      <variable name="Ca" public_interface="in" units="picomolar"/>
      <variable name="Ci" public_interface="in" units="picomolar"/>
      <variable name="P" public_interface="in" units="picomolar"/>
      <variable name="Ra" public_interface="in" units="picomolar"/>
      <variable name="Ri" public_interface="in" units="picomolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>rho</ci>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <ci>Ra</ci>
                  <ci>Ca</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Ra</ci>
                  <ci>Ca</ci>
                  <ci>Ri</ci>
                  <ci>Ci</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="model_parameters">
      <variable cmeta:id="model_parameters_k1" initial_value="1e-6" name="k1" public_interface="out" units="second_order_rate_constant"/>
      <variable initial_value="1e-3" name="k1_" public_interface="out" units="first_order_rate_constant"/>
      <variable cmeta:id="model_parameters_k2" initial_value="1e-7" name="k2" public_interface="out" units="second_order_rate_constant"/>
      <variable initial_value="1e-3" name="k2_" public_interface="out" units="first_order_rate_constant"/>
      <variable cmeta:id="model_parameters_k3" initial_value="1e-3" name="k3" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="1e-4" name="k3_" public_interface="out" units="first_order_rate_constant"/>
      <variable cmeta:id="model_parameters_k4" initial_value="2e-3" name="k4" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.4" name="k4_" public_interface="out" units="first_order_rate_constant"/>
      <variable cmeta:id="model_parameters_kcl" initial_value="5e-3" name="kcl" public_interface="out" units="first_order_rate_constant"/>
      <variable name="D" public_interface="out" units="flux"/>
      <variable initial_value="0" name="De" units="flux"/>
      <variable name="Dd" units="flux"/>
      <variable initial_value="7.5" name="dmax" units="flux"/>
      <variable initial_value="0" name="dmin" units="flux"/>
      <variable initial_value="0.5" name="tau_on" units="hour"/>
      <variable initial_value="0.5" name="tau_off" units="hour"/>
      <variable name="cycle_length" units="hour"/>
      <variable initial_value="9" name="j" units="dimensionless"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="time_hour" units="hour"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>D</ci>
            <apply>
               <plus/>
               <ci>De</ci>
               <ci>Dd</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>cycle_length</ci>
            <apply>
               <plus/>
               <ci>tau_on</ci>
               <ci>tau_off</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>time_hour</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>time</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="hour">1</cn>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">3600</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Dd</ci>
            <piecewise>
               <piece>
                  <ci>dmin</ci>
                  <apply>
                     <geq/>
                     <apply>
                        <divide/>
                        <ci>time_hour</ci>
                        <ci>cycle_length</ci>
                     </apply>
                     <ci>j</ci>
                  </apply>
               </piece>
               <piece>
                  <ci>dmax</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <apply>
                           <minus/>
                           <ci>time_hour</ci>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <ci>time_hour</ci>
                                    <ci>cycle_length</ci>
                                 </apply>
                              </apply>
                              <ci>cycle_length</ci>
                           </apply>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="hour">0</cn>
                     </apply>
                     <apply>
                        <lt/>
                        <apply>
                           <minus/>
                           <ci>time_hour</ci>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <ci>time_hour</ci>
                                    <ci>cycle_length</ci>
                                 </apply>
                              </apply>
                              <ci>cycle_length</ci>
                           </apply>
                        </apply>
                        <ci>tau_on</ci>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <ci>dmin</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <apply>
                           <minus/>
                           <ci>time_hour</ci>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <ci>time_hour</ci>
                                    <ci>cycle_length</ci>
                                 </apply>
                              </apply>
                              <ci>cycle_length</ci>
                           </apply>
                        </apply>
                        <ci>tau_on</ci>
                     </apply>
                     <apply>
                        <lt/>
                        <apply>
                           <minus/>
                           <ci>time_hour</ci>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <ci>time_hour</ci>
                                    <ci>cycle_length</ci>
                                 </apply>
                              </apply>
                              <ci>cycle_length</ci>
                           </apply>
                        </apply>
                        <ci>cycle_length</ci>
                     </apply>
                  </apply>
               </piece>
            </piecewise>
         </apply>
      </math>
   </component>
   <connection>
      <map_components component_1="P" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Ra" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Ri" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Ca" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Ci" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="model_parameters" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Ra" component_2="P"/>
      <map_variables variable_1="P" variable_2="P"/>
      <map_variables variable_1="Ra" variable_2="Ra"/>
   </connection>
   <connection>
      <map_components component_1="Ri" component_2="P"/>
      <map_variables variable_1="P" variable_2="P"/>
      <map_variables variable_1="Ri" variable_2="Ri"/>
   </connection>
   <connection>
      <map_components component_1="Ca" component_2="P"/>
      <map_variables variable_1="P" variable_2="P"/>
      <map_variables variable_1="Ca" variable_2="Ca"/>
   </connection>
   <connection>
      <map_components component_1="Ci" component_2="P"/>
      <map_variables variable_1="P" variable_2="P"/>
      <map_variables variable_1="Ci" variable_2="Ci"/>
   </connection>
   <connection>
      <map_components component_1="P" component_2="model_parameters"/>
      <map_variables variable_1="k1" variable_2="k1"/>
      <map_variables variable_1="k1_" variable_2="k1_"/>
      <map_variables variable_1="k2" variable_2="k2"/>
      <map_variables variable_1="k2_" variable_2="k2_"/>
      <map_variables variable_1="kcl" variable_2="kcl"/>
      <map_variables variable_1="D" variable_2="D"/>
   </connection>
   <connection>
      <map_components component_1="Ca" component_2="Ra"/>
      <map_variables variable_1="Ra" variable_2="Ra"/>
      <map_variables variable_1="Ca" variable_2="Ca"/>
   </connection>
   <connection>
      <map_components component_1="Ri" component_2="Ra"/>
      <map_variables variable_1="Ra" variable_2="Ra"/>
      <map_variables variable_1="Ri" variable_2="Ri"/>
   </connection>
   <connection>
      <map_components component_1="Ra" component_2="model_parameters"/>
      <map_variables variable_1="k1" variable_2="k1"/>
      <map_variables variable_1="k1_" variable_2="k1_"/>
      <map_variables variable_1="k3" variable_2="k3"/>
      <map_variables variable_1="k3_" variable_2="k3_"/>
   </connection>
   <connection>
      <map_components component_1="Ci" component_2="Ri"/>
      <map_variables variable_1="Ri" variable_2="Ri"/>
      <map_variables variable_1="Ci" variable_2="Ci"/>
   </connection>
   <connection>
      <map_components component_1="Ri" component_2="model_parameters"/>
      <map_variables variable_1="k2" variable_2="k2"/>
      <map_variables variable_1="k2_" variable_2="k2_"/>
      <map_variables variable_1="k3" variable_2="k3"/>
      <map_variables variable_1="k3_" variable_2="k3_"/>
   </connection>
   <connection>
      <map_components component_1="Ca" component_2="Ci"/>
      <map_variables variable_1="Ci" variable_2="Ci"/>
      <map_variables variable_1="Ca" variable_2="Ca"/>
   </connection>
   <connection>
      <map_components component_1="Ca" component_2="model_parameters"/>
      <map_variables variable_1="k1" variable_2="k1"/>
      <map_variables variable_1="k1_" variable_2="k1_"/>
      <map_variables variable_1="k4" variable_2="k4"/>
      <map_variables variable_1="k4_" variable_2="k4_"/>
   </connection>
   <connection>
      <map_components component_1="Ci" component_2="model_parameters"/>
      <map_variables variable_1="k2" variable_2="k2"/>
      <map_variables variable_1="k2_" variable_2="k2_"/>
      <map_variables variable_1="k4" variable_2="k4"/>
      <map_variables variable_1="k4_" variable_2="k4_"/>
   </connection>
   <connection>
      <map_components component_1="rho" component_2="P"/>
      <map_variables variable_1="P" variable_2="P"/>
   </connection>
   <connection>
      <map_components component_1="rho" component_2="Ca"/>
      <map_variables variable_1="Ca" variable_2="Ca"/>
   </connection>
   <connection>
      <map_components component_1="rho" component_2="Ci"/>
      <map_variables variable_1="Ci" variable_2="Ci"/>
   </connection>
   <connection>
      <map_components component_1="rho" component_2="Ra"/>
      <map_variables variable_1="Ra" variable_2="Ra"/>
   </connection>
   <connection>
      <map_components component_1="rho" component_2="Ri"/>
      <map_variables variable_1="Ri" variable_2="Ri"/>
   </connection>
</model>