Location: Dixit, Perelson, 2004 @ ca90068e1ced / dixit_perelson_2004.cellml

Author:
Geoff Nunns <gnunns1@jhu.edu>
Date:
2009-12-23 10:58:09+13:00
Desc:
Added value for parameter delta, noted presence of time delays in model and marked it as unsuitable for CellML.
Permanent Source URI:
https://models.cellml.org/workspace/dixit_perelson_2004/rawfile/ca90068e1ced53520ae00df07ef837155412803c/dixit_perelson_2004.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 17/12/2009 at 3:28:21 at p.m. using:

COR (0.9.31.1333)
Copyright 2002-2009 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="dixit_perelson_2004_version02" cmeta:id="dixit_perelson_2004_version02" 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>Modelling the Complex Patterns of Viral Load Decay under Antiretroviral Therapy</title>
				<author>
					<firstname>Catherine</firstname>
					<surname>Lloyd</surname>
					<affiliation>
						<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
            This model is valid CellML but can not be integrated as it is overconstrained. Some variables are not initialised but more importantly the model contains time delays, making it unsuitable for description using CellML. 
          </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
When treated with antiretroviral therapy, the plasma viral loads in HIV patients decline in three distinct phases:
</para>
				<itemizedlist>
					<listitem>
						<para>When the treatment is initially administered there is a 6 to 72 hour delay before the viral loads begin to decline;</para>
					</listitem>
					<listitem>
						<para>following this <emphasis>shoulder</emphasis> period viral load decline rapidly for about one week;</para>
					</listitem>
					<listitem>
						<para>finally, viral loads continue to decrease, but at a slower rate.</para>
					</listitem>
				</itemizedlist>
				<para>
In the Dixit and Perelson 2004 publication described here, the authors present a mathematical model that combines pharmacokinetics and viral dynamics and includes intracellular delay.  The model of viral dynamics is based on that published by <ulink url="${HTML_EXMPL_PERELSON_MODEL}">Perelson <emphasis>et al.</emphasis>
					</ulink> (see <xref linkend="fig_cell_diagram"/> below).  To accurately represent the pharmacokinetics they assume a two compartment model consisting of blood and cells.  The models of viral dynamics and pharmacokinetics are connected via a parameter which represents drug specific efficacy.  Because reverse transcriptase inhibitors (RTIs) and protease inhibitors (PIs) have very different effects on the duration of the intracellular delay, they have been modelled separately (see <xref linkend="fig_reaction_diagram1"/> and <xref linkend="fig_reaction_diagram2"/>).  </para>
				<para>
The complete original paper reference is cited below:
</para>
				<para>
Complex patterns of viral load decay under antiretroviral therapy: influence of pharmacokinetics and intracelllar delay, Narendra M. Dixit and Alan S. Perelson, 2004, 
            <emphasis>Journal of Theoretical Biology</emphasis>
          , 226, 95-109. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=14637059&amp;dopt=Abstract">PubMed ID: 14637059</ulink>
				</para>
				<informalfigure float="0" id="fig_cell_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>cell diagram</title>
							</objectinfo>
							<imagedata fileref="dixit_2004a.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic summary of the dynamics of HIV-1 infection <emphasis>in vivo</emphasis>.  This model is based on that published by Perelson <emphasis>et al.</emphasis> in 1996.</caption>
				</informalfigure>
				<informalfigure float="0" id="fig_reaction_diagram1">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>reaction diagram1</title>
							</objectinfo>
							<imagedata fileref="dixit_2004b.png"/>
						</imageobject>
					</mediaobject>
					<caption>Pharmacokinetics - the two compartment model for protease inhibitors.</caption>
				</informalfigure>
				<informalfigure float="0" id="fig_reaction_diagram2">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>reaction diagram2</title>
							</objectinfo>
							<imagedata fileref="dixit_2004c.png"/>
						</imageobject>
					</mediaobject>
					<caption>Pharmacokinetics - the two compartment model for reverse transcriptase inhibitors.  The main difference between the models is that RTIs must be phosphorylated within the cell in order to be in their active form.</caption>
				</informalfigure>
			</sect1>
		</article>
	</documentation>
	<rdf:RDF>
		<rdf:Description rdf:about="T">
			<dc:title>T</dc:title>
			<dcterms:alternative>uninfected target CD4 cells</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF>
		<rdf:Description rdf:about="T_">
			<dc:title>T_</dc:title>
			<dcterms:alternative>productively infected CD4 cells</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF>
		<rdf:Description rdf:about="VI">
			<dc:title>VI</dc:title>
			<dcterms:alternative>infectious viral load</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF>
		<rdf:Description rdf:about="VNI">
			<dc:title>VNI</dc:title>
			<dcterms:alternative>non-infectious viral load</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF>
		<rdf:Description rdf:about="epsilon_PI">
			<dc:title>instantaneous drug efficacy</dc:title>
			<dcterms:alternative>non-infectious viral load</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Seq rdf:about="rdf:#citationAuthorsSeq">
			<rdf:li rdf:resource="rdf:#author1Vcard"/>
			<rdf:li rdf:resource="rdf:#author2Vcard"/>
		</rdf:Seq>
		<rdf:Description rdf:about="rdf:#author2VcardN">
			<vCard:Given>Alan</vCard:Given>
			<vCard:Family>Perelson</vCard:Family>
			<vCard:Other>S</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#dce94775-74da-46a3-b668-f848194504e4">
			<vCard:Given>Catherine</vCard:Given>
			<vCard:Family>Lloyd</vCard:Family>
			<vCard:Other>May</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#cb2e90cf-8267-4746-a377-b69d2176da09">
			<dc:creator rdf:resource="rdf:#d8541a94-3dba-46fa-8de6-c74b880ef203"/>
			<rdf:value>This is the CellML description of Dixit and Perelson's 2nd 2004 mathematical model of viral dynamics in the presence of reverse transcriptase inhibitors.</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#8eed6c27-f49f-4f5f-bbab-b316a7c04b1c">
			<bqs:subject_type>keyword</bqs:subject_type>
			<rdf:value>
				<rdf:Bag>
					<rdf:li>immunology</rdf:li>
					<rdf:li>antiretroviral therapy</rdf:li>
					<rdf:li>hiv-1</rdf:li>
				</rdf:Bag>
			</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#5506b938-5664-4ec2-aceb-97e13c25d037">
			<dc:subject rdf:resource="rdf:#8eed6c27-f49f-4f5f-bbab-b316a7c04b1c"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#c97f69d1-92da-4b92-bd76-3feda7130f18">
			<dc:title>Journal of Theoretical Biology</dc:title>
		</rdf:Description>
		<rdf:Description rdf:about="#dixit_perelson_2004_version02">
			<dc:title>
        Dixit and Perelson's 2nd 2004 mathematical model of viral dynamics in the presence of reverse transcriptase inhibitors.
      </dc:title>
			<cmeta:comment xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#cb2e90cf-8267-4746-a377-b69d2176da09"/>
			<bqs:reference rdf:resource="rdf:#5506b938-5664-4ec2-aceb-97e13c25d037"/>
			<bqs:reference rdf:resource="rdf:#d35c0aad-add1-45a0-94f8-4ca7a3eba92c"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#author2Vcard">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#author2VcardN"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#8129c007-212b-4bbd-81a4-7f9c740266db">
			<vCard:ORG rdf:resource="rdf:#bb0e69ed-0f9d-439c-8017-0358f2424b9a"/>
			<vCard:EMAIL rdf:resource="rdf:#95098f57-c629-4452-9f43-3b433e7c4312"/>
			<vCard:N rdf:resource="rdf:#dce94775-74da-46a3-b668-f848194504e4"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#author1Vcard">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#author1VcardN"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#0952ec60-fecc-4011-ab80-0c0c2881bcc4">
			<dcterms:W3CDTF>2003-12-12</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#bb0e69ed-0f9d-439c-8017-0358f2424b9a">
			<vCard:Orgname>The University of Auckland</vCard:Orgname>
			<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#97693279-9a92-4dce-a9ae-278b507e4c97">
			<dcterms:W3CDTF>2004-01-07</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#author1VcardN">
			<vCard:Given>Narendra</vCard:Given>
			<vCard:Family>Dixit</vCard:Family>
			<vCard:Other>M</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#95098f57-c629-4452-9f43-3b433e7c4312">
			<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="">
			<dc:publisher>
        The University of Auckland, Bioengineering Institute
      </dc:publisher>
			<dcterms:created rdf:resource="rdf:#0952ec60-fecc-4011-ab80-0c0c2881bcc4"/>
			<dc:creator rdf:resource="rdf:#8129c007-212b-4bbd-81a4-7f9c740266db"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#d8541a94-3dba-46fa-8de6-c74b880ef203">
			<vCard:FN>Catherine Lloyd</vCard:FN>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#3ff37d96-bf85-4897-b91c-50d65c5fca8a">
			<dc:creator rdf:resource="rdf:#citationAuthorsSeq"/>
			<dc:title>Complex patterns of viral load decay under antiretroviral therapy: influence of pharmacokinetics and intracellular delay</dc:title>
			<bqs:volume>226</bqs:volume>
			<bqs:first_page>95</bqs:first_page>
			<bqs:Journal rdf:resource="rdf:#c97f69d1-92da-4b92-bd76-3feda7130f18"/>
			<dcterms:issued rdf:resource="rdf:#97693279-9a92-4dce-a9ae-278b507e4c97"/>
			<bqs:last_page>109</bqs:last_page>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#d35c0aad-add1-45a0-94f8-4ca7a3eba92c">
			<bqs:Pubmed_id>14637059</bqs:Pubmed_id>
			<bqs:JournalArticle rdf:resource="rdf:#3ff37d96-bf85-4897-b91c-50d65c5fca8a"/>
		</rdf:Description>
	</rdf:RDF>   <units xmlns="http://www.cellml.org/cellml/1.0#" name="ml">
      <unit units="litre" prefix="milli"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="mg">
      <unit units="gram" prefix="milli"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_ml">
      <unit units="litre" prefix="milli" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="mg_per_ml">
      <unit units="gram" prefix="milli"/>
      <unit units="litre" prefix="milli" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="day">
      <unit units="second" multiplier="86400"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="flux">
      <unit units="ml"/>
      <unit units="day" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="first_order_rate_constant">
      <unit units="day" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="second_order_rate_constant">
      <unit units="ml" exponent="-1"/>
      <unit units="day" exponent="-1"/>
   </units>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="environment">
      <variable name="time" units="day" public_interface="out"/>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="T">
      <variable name="T" units="per_ml" initial_value="1e6" public_interface="out"/>
      <variable name="lamda" units="second_order_rate_constant" initial_value="1e4"/>
      <variable name="d" units="first_order_rate_constant" initial_value="0.01"/>
      <variable name="k" units="flux" public_interface="in"/>
      <variable name="VI" units="per_ml" public_interface="in"/>
      <variable name="time" units="day" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>T</ci>
            </apply>
            <apply>
               <minus/>
               <ci>lamda</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>d</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k</ci>
                     <ci>T</ci>
                     <ci>VI</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="T_">
      <variable name="T_" units="per_ml" initial_value="1" public_interface="out"/>
      <variable name="tau" units="first_order_rate_constant" initial_value="1.5"/>
      <variable name="m" units="first_order_rate_constant" initial_value="0.01"/>
      <variable name="k" units="flux" public_interface="in"/>
      <variable name="VI" units="per_ml" public_interface="in"/>
      <variable name="T" units="per_ml" public_interface="in"/>
      <variable name="delta" units="first_order_rate_constant" public_interface="in"/>
      <variable name="time" units="day" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>T_</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k</ci>
                  <ci>T</ci>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <ci>tau</ci>
                  </apply>
                  <ci>VI</ci>
                  <apply>
                     <minus/>
                     <ci>time</ci>
                     <ci>tau</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>m</ci>
                        </apply>
                        <ci>tau</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>delta</ci>
                  <ci>T_</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="VI">
      <variable name="VI" units="per_ml" initial_value="1" public_interface="out"/>
      <variable name="N" units="dimensionless" public_interface="in"/>
      <variable name="c" units="first_order_rate_constant" public_interface="in"/>
      <variable name="epsilon_PI" units="dimensionless" public_interface="in"/>
      <variable name="T_" units="per_ml" public_interface="in"/>
      <variable name="delta" units="first_order_rate_constant" public_interface="in"/>
      <variable name="time" units="day" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>VI</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>N</ci>
                  <ci>delta</ci>
                  <ci>T_</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>epsilon_PI</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>c</ci>
                  <ci>VI</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="VNI">
      <variable name="VNI" units="per_ml" initial_value="2" public_interface="out"/>
      <variable name="N" units="dimensionless" public_interface="in"/>
      <variable name="c" units="first_order_rate_constant" public_interface="in"/>
      <variable name="epsilon_PI" units="dimensionless" public_interface="in"/>
      <variable name="T_" units="per_ml" public_interface="in"/>
      <variable name="delta" units="first_order_rate_constant" public_interface="in"/>
      <variable name="time" units="day" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>VNI</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>N</ci>
                  <ci>delta</ci>
                  <ci>T_</ci>
                  <ci>epsilon_PI</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>c</ci>
                  <ci>VNI</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="epsilon_PI">
      <variable name="epsilon_PI" units="dimensionless" public_interface="out"/>
      <variable name="IC50" units="mg_per_ml" initial_value="9e-7"/>
      <variable name="Cc" units="mg_per_ml" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>epsilon_PI</ci>
            <apply>
               <divide/>
               <ci>Cc</ci>
               <apply>
                  <plus/>
                  <ci>IC50</ci>
                  <ci>Cc</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Cb">
      <variable name="Cb" units="mg_per_ml" public_interface="out"/>
      <variable name="Vd" units="ml" initial_value="28000"/>
      <variable name="F" units="dimensionless" initial_value="1"/>
      <variable name="D" units="mg" initial_value="600"/>
      <variable name="ka" units="first_order_rate_constant" public_interface="in"/>
      <variable name="ke" units="first_order_rate_constant" public_interface="in"/>
      <variable name="time" units="day" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Cb</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>F</ci>
                     <ci>D</ci>
                  </apply>
                  <ci>Vd</ci>
               </apply>
               <apply>
                  <divide/>
                  <ci>ka</ci>
                  <apply>
                     <minus/>
                     <ci>ke</ci>
                     <ci>ka</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>ka</ci>
                        </apply>
                        <ci>time</ci>
                     </apply>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>ke</ci>
                        </apply>
                        <ci>time</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Cc">
      <variable name="Cc" units="mg_per_ml" initial_value="0" public_interface="out"/>
      <variable name="kacell" units="first_order_rate_constant" initial_value="24000"/>
      <variable name="kecell" units="first_order_rate_constant" initial_value="1.1"/>
      <variable name="time" units="day" public_interface="in"/>
      <variable name="Cx" units="mg_per_ml" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cc</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kacell</ci>
                  <ci>Cx</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>kecell</ci>
                  <ci>Cc</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Cx">
      <variable name="Cx" units="mg_per_ml" public_interface="out"/>
      <variable name="H" units="dimensionless" initial_value="0.052"/>
      <variable name="fb" units="dimensionless" initial_value="0.99"/>
      <variable name="Cb" units="mg_per_ml" public_interface="in"/>
      <variable name="Cc" units="mg_per_ml" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Cx</ci>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>fb</ci>
                        </apply>
                        <ci>H</ci>
                        <ci>Cb</ci>
                     </apply>
                     <ci>Cc</ci>
                  </apply>
                  <apply>
                     <gt/>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                              <ci>fb</ci>
                           </apply>
                           <ci>H</ci>
                           <ci>Cb</ci>
                        </apply>
                        <ci>Cc</ci>
                     </apply>
                     <cn cellml:units="mg_per_ml">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="mg_per_ml">0</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="kinetic_parameters">
      <variable name="N" units="dimensionless" initial_value="2500" public_interface="out"/>
      <variable name="c" units="first_order_rate_constant" initial_value="23" public_interface="out"/>
      <variable name="k" units="flux" initial_value="2.4e-8" public_interface="out"/>
      <variable name="delta" units="first_order_rate_constant" initial_value="0.01" public_interface="out"/>
      <variable name="ka" units="first_order_rate_constant" initial_value="14.64" public_interface="out"/>
      <variable name="ke" units="first_order_rate_constant" initial_value="6.86" public_interface="out"/>
   </component>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="T" 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="T_" 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="VI" 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="VNI" 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="Cb" 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="Cc" 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="T" component_2="T_"/>
      <map_variables variable_1="T" variable_2="T"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="T" component_2="VI"/>
      <map_variables variable_1="VI" variable_2="VI"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="T_" component_2="VI"/>
      <map_variables variable_1="VI" variable_2="VI"/>
      <map_variables variable_1="T_" variable_2="T_"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="T_" component_2="VNI"/>
      <map_variables variable_1="T_" variable_2="T_"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="T" component_2="kinetic_parameters"/>
      <map_variables variable_1="k" variable_2="k"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="T_" component_2="kinetic_parameters"/>
      <map_variables variable_1="delta" variable_2="delta"/>
      <map_variables variable_1="k" variable_2="k"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="VI" component_2="kinetic_parameters"/>
      <map_variables variable_1="delta" variable_2="delta"/>
      <map_variables variable_1="c" variable_2="c"/>
      <map_variables variable_1="N" variable_2="N"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="VI" component_2="epsilon_PI"/>
      <map_variables variable_1="epsilon_PI" variable_2="epsilon_PI"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="VNI" component_2="kinetic_parameters"/>
      <map_variables variable_1="delta" variable_2="delta"/>
      <map_variables variable_1="c" variable_2="c"/>
      <map_variables variable_1="N" variable_2="N"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="VNI" component_2="epsilon_PI"/>
      <map_variables variable_1="epsilon_PI" variable_2="epsilon_PI"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="epsilon_PI" component_2="Cc"/>
      <map_variables variable_1="Cc" variable_2="Cc"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Cb" component_2="kinetic_parameters"/>
      <map_variables variable_1="ka" variable_2="ka"/>
      <map_variables variable_1="ke" variable_2="ke"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Cc" component_2="Cx"/>
      <map_variables variable_1="Cx" variable_2="Cx"/>
      <map_variables variable_1="Cc" variable_2="Cc"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Cx" component_2="Cb"/>
      <map_variables variable_1="Cb" variable_2="Cb"/>
   </connection>
</model>