Location: Liu, Hu, Peng, Liu, 1999 @ 82357f460156 / liu_hu_peng_liu_1999.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2010-04-08 15:27:12+12:00
Desc:
Made minor changes to the documentation.
Permanent Source URI:
https://staging.physiomeproject.org/workspace/liu_hu_peng_liu_1999/rawfile/82357f4601563d1f97485d83e06a8b507147fd03/liu_hu_peng_liu_1999.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 30/03/2010 at 12:44:59 at p.m. using:

COR (0.9.31.1371)
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:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="liu_model_1999" cmeta:id="liu_model_1999" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Mathematical modeling and qualitative analysis of insulin therapies</title>
				<author>
					<firstname>Catherine</firstname>
					<surname>Lloyd</surname>
					<affiliation>
						<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
				This CellML model runs in PCenv, COR and OpenCell to recreate the published results. The initial conditions were not explicitly defined in the paper, but were extrapolated from Table 1 and Figure 2. The model was then run until it reached a steady state, and new intial values were taken from the steady state. 
				</para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
ABSTRACT:  We propose a dynamical model for the pulsatile secretion of the hypothalamo-pituitary-adrenal axis. This model takes into account both the binding of hormone with proteins in the plasma and tissues and mutual interactions of the hormones in the system. The deductions from 
this model are in good agreement with experiment results.
</para>
				<informalfigure float="0" id="fig_reaction_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>model diagram</title>
							</objectinfo>
							<imagedata fileref="liu_1999.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the hormone flow and interactions between the different components of the hypothalamic-pituitary-adrenal (HPA) axis which are described by the mathematical model.  CRH represents corticotropin releasing hormone, ACTH represents corticotropin, and CBG is corticosteroid binding globulin.  Red arrows represent hormone disposal, while black arrows represent flow between the different glands and plasma. + represents a stimulatory effect while - represents an inhibitory effect.</caption>
				</informalfigure>
				
				<para>
The original paper reference is cited below:
</para>
				<para>
A dynamical model for the pulsatile secretion of the hypothalamo-pituitary-adrenal axis, Yi-Wei Liu, Zhi-Hong Hu, Jian-Hua Peng, and Bing-Zheng Liu, 1999, Mathematical and Computer Modelling, 29 (Issue 4), 103-110.  (No Pubmed ID)  
</para>
			</sect1>
		</article>
	</documentation>
	
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
		<rdf:Description rdf:about="">
			<dc:creator rdf:parseType="Resource">
				<vCard:N rdf:parseType="Resource">
					<vCard:Family>Lloyd</vCard:Family>
					<vCard:Given>Catherine</vCard:Given>
					<vCard:Other>May</vCard:Other>
				</vCard:N>
				<vCard:EMAIL rdf:parseType="Resource">
					<rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
					<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
				</vCard:EMAIL>
				<vCard:ORG rdf:parseType="Resource">
					<vCard:Orgname>The University of Auckland</vCard:Orgname>
					<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
				</vCard:ORG>
			</dc:creator>
			<dcterms:created rdf:parseType="Resource">
				<dcterms:W3CDTF>2007-07-13</dcterms:W3CDTF>
			</dcterms:created>
			<dc:publisher>
        The University of Auckland, Auckland Bioengineering Institute
      </dc:publisher>
		</rdf:Description>
		<rdf:Description rdf:about="#liu_model_1999">
			<dc:title>
        Liu et al's 1999 mathematical model for the pulsatile secretion of the hypothalamic-pituitary-adrenal axis.
      </dc:title>
			<cmeta:comment xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:parseType="Resource">
				<rdf:value>
          This is a CellML description of Liu et al's 1999 mathematical model for the pulsatile secretion of the hypothalamic-pituitary-adrenal axis.
        </rdf:value>
				<dc:creator rdf:parseType="Resource">
					<vCard:FN>Catherine Lloyd</vCard:FN>
				</dc:creator>
			</cmeta:comment>
			<bqs:reference rdf:parseType="Resource">
				<dc:subject rdf:parseType="Resource">
					<bqs:subject_type>keyword</bqs:subject_type>
					<rdf:value>
						<rdf:Bag>
							<rdf:li>cortisol</rdf:li>
							<rdf:li>endocrine</rdf:li>
							<rdf:li>adrenal cortex</rdf:li>
							<rdf:li>HPA axis</rdf:li>
							<rdf:li>hypothalamus</rdf:li>
							<rdf:li>hpa axis</rdf:li>
							<rdf:li>anterior pituitary</rdf:li>
						</rdf:Bag>
					</rdf:value>
				</dc:subject>
			</bqs:reference>
			<bqs:reference rdf:parseType="Resource">
				<bqs:JournalArticle rdf:parseType="Resource">
					<dc:creator>
						<rdf:Seq>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N rdf:parseType="Resource">
										<vCard:Family>Liu</vCard:Family>
										<vCard:Given>Yi-Wei</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N rdf:parseType="Resource">
										<vCard:Family>Hu</vCard:Family>
										<vCard:Given>Zhi-Hong</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N rdf:parseType="Resource">
										<vCard:Family>Peng</vCard:Family>
										<vCard:Given>Jian-Hua</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N rdf:parseType="Resource">
										<vCard:Family>Liu</vCard:Family>
										<vCard:Given>Bing-Zheng</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
						</rdf:Seq>
					</dc:creator>
					<dc:title>
            A dynamical model for the pulsatile secretion of the hypothalamo-pituitary-adrenal axis
          </dc:title>
					<dcterms:issued rdf:parseType="Resource">
						<dcterms:W3CDTF>1999-02</dcterms:W3CDTF>
					</dcterms:issued>
					<bqs:Journal rdf:parseType="Resource">
						<dc:title>
              Mathematical and Computer Modelling
            </dc:title>
					</bqs:Journal>
					<bqs:volume>29</bqs:volume>
					<bqs:issue>4</bqs:issue>
					<bqs:first_page>103</bqs:first_page>
					<bqs:last_page>110</bqs:last_page>
				</bqs:JournalArticle>
			</bqs:reference>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/">
		<rdf:Description rdf:about="#x1">
			<dc:title>x1</dc:title>
			<dcterms:alternative>CRH</dcterms:alternative>
			<dcterms:alternative>corticotropin-releasing hormone</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/">
		<rdf:Description rdf:about="#x2">
			<dc:title>x2</dc:title>
			<dcterms:alternative>ACTH</dcterms:alternative>
			<dcterms:alternative>corticotropin</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/">
		<rdf:Description rdf:about="#x3">
			<dc:title>x3</dc:title>
			<dcterms:alternative>free cortisol</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/">
		<rdf:Description rdf:about="#x4">
			<dc:title>x4</dc:title>
			<dcterms:alternative>CBG-bound cortisol</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/">
		<rdf:Description rdf:about="#x5">
			<dc:title>x5</dc:title>
			<dcterms:alternative>albumin-bound cortisol</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	
	<units xmlns="http://www.cellml.org/cellml/1.0#" name="minute">
      <unit units="second" multiplier="60"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="microg_l">
      <unit units="gram" prefix="micro"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_microg_l">
      <unit units="microg_l" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_microg_l2">
      <unit units="microg_l" exponent="-2"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="flux">
      <unit units="microg_l"/>
      <unit units="minute" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="first_order_rate_constant">
      <unit units="minute" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="second_order_rate_constant">
      <unit units="microg_l" exponent="-1"/>
      <unit units="minute" exponent="-1"/>
   </units>
   
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="environment">
      <variable name="time" units="minute" public_interface="out"/>
   </component>
   
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="x1">
      <variable name="x1" units="microg_l" initial_value="0.01067" public_interface="out"/>
      <variable name="lambda_1" units="first_order_rate_constant" public_interface="in"/>
      <variable name="x3" units="microg_l" public_interface="in"/>
      <variable name="a1" units="flux" public_interface="in"/>
      <variable name="a2" units="flux" public_interface="in"/>
      <variable name="a3" units="first_order_rate_constant" public_interface="in"/>
      <variable name="a4" units="second_order_rate_constant" public_interface="in"/>
      <variable name="a5" units="per_microg_l" public_interface="in"/>
      <variable name="a6" units="per_microg_l2" public_interface="in"/>
      <variable name="a7" units="per_microg_l" public_interface="in"/>
      <variable name="a8" units="per_microg_l2" public_interface="in"/>
      <variable name="time" units="minute" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>a1</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>a2</ci>
                        <apply>
                           <times/>
                           <ci>a3</ci>
                           <ci>x1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a4</ci>
                           <apply>
                              <power/>
                              <ci>x1</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>a5</ci>
                           <ci>x1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a6</ci>
                           <apply>
                              <power/>
                              <ci>x1</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a7</ci>
                           <ci>x3</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a8</ci>
                           <apply>
                              <power/>
                              <ci>x3</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_1</ci>
                  <ci>x1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="x2">
      <variable name="x2" units="microg_l" initial_value="0.04665" public_interface="out"/>
      <variable name="lambda_2" units="first_order_rate_constant" public_interface="in"/>
      <variable name="x1" units="microg_l" public_interface="in"/>
      <variable name="x3" units="microg_l" public_interface="in"/>
      <variable name="a9" units="flux" public_interface="in"/>
      <variable name="a10" units="first_order_rate_constant" public_interface="in"/>
      <variable name="a11" units="second_order_rate_constant" public_interface="in"/>
      <variable name="a12" units="per_microg_l" public_interface="in"/>
      <variable name="a13" units="per_microg_l2" public_interface="in"/>
      <variable name="a14" units="per_microg_l" public_interface="in"/>
      <variable name="a15" units="per_microg_l2" public_interface="in"/>
      <variable name="time" units="minute" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <plus/>
                     <ci>a9</ci>
                     <apply>
                        <times/>
                        <ci>a10</ci>
                        <ci>x1</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>a11</ci>
                        <apply>
                           <power/>
                           <ci>x1</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <ci>a12</ci>
                        <ci>x1</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>a13</ci>
                        <apply>
                           <power/>
                           <ci>x1</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>a14</ci>
                        <ci>x3</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>a15</ci>
                        <apply>
                           <power/>
                           <ci>x3</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_2</ci>
                  <ci>x2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="x3">
      <variable name="x3" units="microg_l" initial_value="6.51" public_interface="out"/>
      <variable name="lambda_3_" units="first_order_rate_constant" public_interface="in"/>
      <variable name="x1" units="microg_l" public_interface="in"/>
      <variable name="x2" units="microg_l" public_interface="in"/>
      <variable name="x4" units="microg_l" public_interface="in"/>
      <variable name="x5" units="microg_l" public_interface="in"/>
      <variable name="a16" units="flux" public_interface="in"/>
      <variable name="a17" units="first_order_rate_constant" public_interface="in"/>
      <variable name="a18" units="second_order_rate_constant" public_interface="in"/>
      <variable name="a19" units="first_order_rate_constant" public_interface="in"/>
      <variable name="a20" units="second_order_rate_constant" public_interface="in"/>
      <variable name="a21" units="per_microg_l" public_interface="in"/>
      <variable name="a22" units="per_microg_l2" public_interface="in"/>
      <variable name="a23" units="per_microg_l" public_interface="in"/>
      <variable name="a24" units="per_microg_l2" public_interface="in"/>
      <variable name="a25" units="first_order_rate_constant" public_interface="in"/>
      <variable name="a26" units="first_order_rate_constant" public_interface="in"/>
      <variable name="time" units="minute" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x3</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>a16</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>a17</ci>
                           <ci>x1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a18</ci>
                           <apply>
                              <power/>
                              <ci>x1</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a19</ci>
                           <ci>x2</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a20</ci>
                           <apply>
                              <power/>
                              <ci>x2</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>a21</ci>
                           <ci>x1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a22</ci>
                           <apply>
                              <power/>
                              <ci>x1</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a23</ci>
                           <ci>x2</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>a24</ci>
                           <apply>
                              <power/>
                              <ci>x2</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>a25</ci>
                     <ci>x4</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>a26</ci>
                     <ci>x5</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_3_</ci>
                  <ci>x3</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="x4">
      <variable name="x4" units="microg_l" initial_value="60.61" public_interface="out"/>
      <variable name="lambda_4_" units="first_order_rate_constant" public_interface="in"/>
      <variable name="x3" units="microg_l" public_interface="in"/>
      <variable name="a27" units="first_order_rate_constant" public_interface="in"/>
      <variable name="time" units="minute" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x4</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>a27</ci>
                  <ci>x3</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_4_</ci>
                  <ci>x4</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="x5">
      <variable name="x5" units="microg_l" initial_value="12.61" public_interface="out"/>
      <variable name="lambda_5_" units="first_order_rate_constant" public_interface="in"/>
      <variable name="x3" units="microg_l" public_interface="in"/>
      <variable name="a28" units="first_order_rate_constant" public_interface="in"/>
      <variable name="time" units="minute" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x5</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>a28</ci>
                  <ci>x3</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_5_</ci>
                  <ci>x5</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="model_parameters">
      <variable name="a1" units="flux" initial_value="0.000017" public_interface="out"/>
      <variable name="a2" units="flux" initial_value="0.0023" public_interface="out"/>
      <variable name="a3" units="first_order_rate_constant" initial_value="0.6" public_interface="out"/>
      <variable name="a4" units="second_order_rate_constant" initial_value="45" public_interface="out"/>
      <variable name="a5" units="per_microg_l" initial_value="36" public_interface="out"/>
      <variable name="a6" units="per_microg_l2" initial_value="216" public_interface="out"/>
      <variable name="a7" units="per_microg_l" initial_value="0.28" public_interface="out"/>
      <variable name="a8" units="per_microg_l2" initial_value="0.36" public_interface="out"/>
      <variable name="a9" units="flux" initial_value="0.0003" public_interface="out"/>
      <variable name="a10" units="first_order_rate_constant" initial_value="0.18" public_interface="out"/>
      <variable name="a11" units="second_order_rate_constant" initial_value="150" public_interface="out"/>
      <variable name="a12" units="per_microg_l" initial_value="18" public_interface="out"/>
      <variable name="a13" units="per_microg_l2" initial_value="460" public_interface="out"/>
      <variable name="a14" units="per_microg_l" initial_value="0.46" public_interface="out"/>
      <variable name="a15" units="per_microg_l2" initial_value="0.1" public_interface="out"/>
      <variable name="a16" units="flux" initial_value="0.04" public_interface="out"/>
      <variable name="a17" units="first_order_rate_constant" initial_value="150" public_interface="out"/>
      <variable name="a18" units="second_order_rate_constant" initial_value="3800" public_interface="out"/>
      <variable name="a19" units="first_order_rate_constant" initial_value="57" public_interface="out"/>
      <variable name="a20" units="second_order_rate_constant" initial_value="2600" public_interface="out"/>
      <variable name="a21" units="per_microg_l" initial_value="200" public_interface="out"/>
      <variable name="a22" units="per_microg_l2" initial_value="9400" public_interface="out"/>
      <variable name="a23" units="per_microg_l" initial_value="10" public_interface="out"/>
      <variable name="a24" units="per_microg_l2" initial_value="320" public_interface="out"/>
      <variable name="a25" units="first_order_rate_constant" initial_value="0.04" public_interface="out"/>
      <variable name="a26" units="first_order_rate_constant" initial_value="0.00097" public_interface="out"/>
      <variable name="a27" units="first_order_rate_constant" initial_value="0.57" public_interface="out"/>
      <variable name="a28" units="first_order_rate_constant" initial_value="0.0017" public_interface="out"/>
      <variable name="lambda_1" units="first_order_rate_constant" initial_value="0.059" public_interface="out"/>
      <variable name="lambda_2" units="first_order_rate_constant" initial_value="0.028" public_interface="out"/>
      <variable name="lambda_3_" units="first_order_rate_constant" public_interface="out"/>
      <variable name="lambda_4_" units="first_order_rate_constant" public_interface="out"/>
      <variable name="lambda_5_" units="first_order_rate_constant" public_interface="out"/>
	  
	  <variable name="lambda_3" units="first_order_rate_constant" initial_value="0.0986"/>
      <variable name="lambda_4" units="first_order_rate_constant" initial_value="0.024"/>
      <variable name="lambda_5" units="first_order_rate_constant" initial_value="3e-5"/>
   
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>lambda_3_</ci>
            <apply>
               <plus/>
               <ci>lambda_3</ci>
               <ci>a27</ci>
               <ci>a28</ci>
            </apply>
         </apply>
      </math>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>lambda_4_</ci>
            <apply>
               <plus/>
               <ci>lambda_4</ci>
               <ci>a25</ci>
            </apply>
         </apply>
      </math>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>lambda_5_</ci>
            <apply>
               <plus/>
               <ci>lambda_5</ci>
               <ci>a26</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x2" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x3" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x4" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x5" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x1" component_2="x3"/>
      <map_variables variable_1="x1" variable_2="x1"/>
      <map_variables variable_1="x3" variable_2="x3"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x1" component_2="x2"/>
      <map_variables variable_1="x1" variable_2="x1"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x1" component_2="model_parameters"/>
      <map_variables variable_1="a1" variable_2="a1"/>
      <map_variables variable_1="a2" variable_2="a2"/>
      <map_variables variable_1="a3" variable_2="a3"/>
      <map_variables variable_1="a4" variable_2="a4"/>
      <map_variables variable_1="a5" variable_2="a5"/>
      <map_variables variable_1="a6" variable_2="a6"/>
      <map_variables variable_1="a7" variable_2="a7"/>
      <map_variables variable_1="a8" variable_2="a8"/>
      <map_variables variable_1="lambda_1" variable_2="lambda_1"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x2" component_2="x3"/>
      <map_variables variable_1="x2" variable_2="x2"/>
      <map_variables variable_1="x3" variable_2="x3"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x2" component_2="model_parameters"/>
      <map_variables variable_1="a9" variable_2="a9"/>
      <map_variables variable_1="a10" variable_2="a10"/>
      <map_variables variable_1="a11" variable_2="a11"/>
      <map_variables variable_1="a12" variable_2="a12"/>
      <map_variables variable_1="a13" variable_2="a13"/>
      <map_variables variable_1="a14" variable_2="a14"/>
      <map_variables variable_1="a15" variable_2="a15"/>
      <map_variables variable_1="lambda_2" variable_2="lambda_2"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x4" component_2="x3"/>
      <map_variables variable_1="x4" variable_2="x4"/>
      <map_variables variable_1="x3" variable_2="x3"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x5" component_2="x3"/>
      <map_variables variable_1="x5" variable_2="x5"/>
      <map_variables variable_1="x3" variable_2="x3"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x3" component_2="model_parameters"/>
      <map_variables variable_1="a16" variable_2="a16"/>
      <map_variables variable_1="a17" variable_2="a17"/>
      <map_variables variable_1="a18" variable_2="a18"/>
      <map_variables variable_1="a19" variable_2="a19"/>
      <map_variables variable_1="a20" variable_2="a20"/>
      <map_variables variable_1="a21" variable_2="a21"/>
      <map_variables variable_1="a22" variable_2="a22"/>
      <map_variables variable_1="a23" variable_2="a23"/>
      <map_variables variable_1="a24" variable_2="a24"/>
      <map_variables variable_1="a25" variable_2="a25"/>
      <map_variables variable_1="a26" variable_2="a26"/>
      <map_variables variable_1="lambda_3_" variable_2="lambda_3_"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x4" component_2="model_parameters"/>
      <map_variables variable_1="a27" variable_2="a27"/>
      <map_variables variable_1="lambda_4_" variable_2="lambda_4_"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="x5" component_2="model_parameters"/>
      <map_variables variable_1="a28" variable_2="a28"/>
      <map_variables variable_1="lambda_5_" variable_2="lambda_5_"/>
   </connection>
</model>