Rendering of the source text

<?xml version="1.0"?>

<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="warren_2010" name="warren_2010">

	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Mathematical modelling of calcium wave propagation in mammalian airway epithelium: evidence for regenerative ATP release</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 both COR and OpenCell to replicate the results from the original published model.  The units have been checked and they balance.  We'd like to thank the original model author Nic Warren for taking the time to develop the CellML version of his model.
				</para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
ABSTRACT: Airway epithelium has been shown to exhibit intracellular calcium waves after mechanical stimulation. Two classes of mechanism have been proposed to explain calcium wave propagation: diffusion through gap junctions of the intracellular messenger inositol 1,4,5-trisphosphate (IP3), and diffusion of paracrine extracellular messengers such as ATP. We have used single cell recordings of airway epithelium to parametreize a model of an airway epithelial cell. This was then incorporated into a spatial model of a cell culture where both mechanisms for calcium wave propagation are possible. It is shown that a decreasing return on the radius of Ca2+ wave propagation is achieved as the amount of ATP released from the stimulated cell increases. It is therefore shown that for a Ca2+ wave to propagate large distances, a significant fraction of the intracellular ATP pool would be required to be released. Further to this, the radial distribution of maximal calcium response from the stimulated cell does not produce the same flat profile of maximal calcium response seen in experiential studies. This suggests that an additional mechanism is important in Ca2+ wave propagation, such as regenerative release of ATP from cells downstream of the stimulated cell.
</para>
				<para>
The original paper reference is cited below:
</para>
				<para>Mathematical modelling of calcium wave propagation in mammalian airway epithelium: evidence for regenerative ATP release, N.J. Warren, M.H. Tawhai and E.J. Crampin, 2010,
						<emphasis>Experimental Physiology</emphasis> 95, 232-249. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/19700517">PubMed ID: 19700517</ulink>
				</para>
				
				<informalfigure float="0" id="diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>diagram</title>
							</objectinfo>
							<imagedata fileref="warren_2010.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic representation of model components and subsystems: Mechanical stimulation triggers the release of ATP from the stimulated cell as well as direct opening of plasma membrane Ca2+ channels (PMCCs). The release of ATP results in the production of IP3 after binding to cell-surface P2Y2 receptors. The IP3 can diffuse through the cytosol to initiate the release of Ca2+ from the ER or diffuse through gap junctions to initiate Ca2+ release from the ER of an adjacent cell. 
</caption>
				</informalfigure>
			</sect1>
		</article>
	</documentation>

	<units name="micron">
      <unit prefix="micro" units="meter"/>
   </units>
   <units name="volume">
      <unit exponent="3" units="micron"/>
   </units>
   <units name="milliV">
      <unit prefix="milli" units="volt"/>
   </units>
   <units name="millijoule_per_mole_kelvin">
      <unit prefix="milli" units="joule"/>
      <unit exponent="-1" units="mole"/>
      <unit exponent="-1" units="kelvin"/>
   </units>
   <units name="area">
      <unit exponent="2" units="micron"/>
   </units>
   <units name="capacitance">
      <unit prefix="pico" units="farad"/>
      <unit exponent="-1" units="area"/>
   </units>
   <units name="current">
      <unit prefix="-21" units="ampere"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="ampere"/>
      <unit units="second"/>
      <unit exponent="-1" units="mole"/>
   </units>
   <units name="conductance">
      <unit prefix="-18" units="siemens"/>
   </units>
   <units name="conductancepervolume">
      <unit units="conductance"/>
      <unit exponent="-1" units="volume"/>
   </units>
   <units name="conductancesec2permicroM2volume">
      <unit units="conductance"/>
      <unit exponent="-1" units="volume"/>
      <unit exponent="2" units="second"/>
      <unit exponent="-2" units="microM"/>
   </units>
   <units name="mM">
      <unit prefix="milli" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="mMvolume">
      <unit units="mM"/>
      <unit units="volume"/>
   </units>
   <units name="mMvolumepersec">
      <unit units="mM"/>
      <unit units="volume"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="microMvolumepersec">
      <unit units="microM"/>
      <unit units="volume"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="microMpervolumesec">
      <unit units="microM"/>
      <unit exponent="-1" units="volume"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="microM2persec">
      <unit exponent="2" units="microM"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="microM">
      <unit prefix="micro" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="per_microM">
      <unit exponent="-1" units="microM"/>
   </units>
   <units name="microMpersec">
      <unit units="microM"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="microMVolume">
      <unit units="microM"/>
      <unit units="volume"/>
   </units>
   <units name="Volumepersec">
      <unit exponent="-1" units="second"/>
      <unit units="volume"/>
   </units>
   <units name="mmol">
      <unit prefix="milli" units="mole"/>
   </units>
   <units name="per_mmol">
      <unit exponent="-1" units="mmol"/>
   </units>
   <units name="per_sec">
      <unit exponent="-1" units="second"/>
   </units>
   <units name="per_secmicroM">
      <unit exponent="-1" units="second"/>
      <unit exponent="-1" units="microM"/>
   </units>
   <units name="micromolar">
      <unit prefix="micro" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="flux">
      <unit units="microM"/>
      <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="micromolar"/>
      <unit exponent="-1" units="second"/>
   </units>
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="second"/>
   </component>
   <component name="volume_i">
      <variable initial_value="1450" name="volume_i0" public_interface="out" units="volume"/>
      <variable name="volume_i" public_interface="out" units="volume"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>volume_i</ci>
            <ci>volume_i0</ci>
         </apply>
      </math>
   </component>
   <component name="volume_er">
      <variable initial_value="4.9" name="gamma" units="dimensionless"/>
      <variable name="volume_i" public_interface="in" units="volume"/>
      <variable name="volume_er" public_interface="out" units="volume"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>volume_er</ci>
            <apply>
               <divide/>
               <ci>volume_i</ci>
               <ci>gamma</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="volumeCa_i">
      <variable initial_value="75.8335890358239" name="volumeCa_i" public_interface="out" units="microMVolume"/>
      <variable name="volume_i" public_interface="in" units="volume"/>
      <variable name="Jipt" public_interface="in" units="microMvolumepersec"/>
      <variable name="Jryr" public_interface="in" units="microMvolumepersec"/>
      <variable name="Jer" public_interface="in" units="microMvolumepersec"/>
      <variable name="Jserca" public_interface="in" units="microMvolumepersec"/>
      <variable name="Jin" public_interface="in" units="microMvolumepersec"/>
      <variable name="Jpm" public_interface="in" units="microMvolumepersec"/>
      <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>volumeCa_i</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>Jipt</ci>
                     <ci>Jryr</ci>
                     <ci>Jer</ci>
                  </apply>
                  <ci>Jserca</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>Jin</ci>
                  <ci>Jpm</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_i">
      <variable name="Ca_i" public_interface="out" units="microM"/>
      <variable name="volume_i" public_interface="in" units="volume"/>
      <variable name="volumeCa_i" public_interface="in" units="microMVolume"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Ca_i</ci>
            <apply>
               <divide/>
               <ci>volumeCa_i</ci>
               <ci>volume_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="volumeCa_er">
      <variable initial_value="1199.73747802144" name="volumeCa_er" public_interface="out" units="microMVolume"/>
      <variable name="volume_er" public_interface="in" units="volume"/>
      <variable name="Jipt" public_interface="in" units="microMvolumepersec"/>
      <variable name="Jryr" public_interface="in" units="microMvolumepersec"/>
      <variable name="Jer" public_interface="in" units="microMvolumepersec"/>
      <variable name="Jserca" public_interface="in" units="microMvolumepersec"/>
      <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>volumeCa_er</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>Jipt</ci>
                     <ci>Jryr</ci>
                     <ci>Jer</ci>
                  </apply>
                  <ci>Jserca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_er">
      <variable cmeta:id="Ca_er_Ca_er" name="Ca_er" public_interface="out" units="microM"/>
      <variable name="volume_er" public_interface="in" units="volume"/>
      <variable name="volumeCa_er" public_interface="in" units="microMVolume"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Ca_er</ci>
            <apply>
               <divide/>
               <ci>volumeCa_er</ci>
               <ci>volume_er</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="volume_iIP3">
      <variable initial_value="-3.28287434380887e-93" name="volume_iIP3" public_interface="out" units="microMVolume"/>
      <variable name="volume_i" public_interface="in" units="volume"/>
      <variable name="J_ip3P" public_interface="in" units="microMpersec"/>
      <variable name="J_ip3D" public_interface="in" units="microMpersec"/>
      <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>volume_iIP3</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>J_ip3P</ci>
                  <ci>J_ip3D</ci>
               </apply>
               <ci>volume_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="IP3">
      <variable cmeta:id="IP3_IP3" name="IP3" public_interface="out" units="microM"/>
      <variable name="volume_i" public_interface="in" units="volume"/>
      <variable name="volume_iIP3" public_interface="in" units="microMVolume"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>IP3</ci>
            <apply>
               <divide/>
               <ci>volume_iIP3</ci>
               <ci>volume_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="PIP2">
      <variable initial_value="1" name="PIP2" public_interface="out" units="dimensionless"/>
      <variable initial_value="5" name="PIP2_Total" units="microM"/>
      <variable name="J_ip3P" public_interface="in" units="microMpersec"/>
      <variable name="J_ip3R" public_interface="in" units="per_sec"/>
      <variable name="IPX" public_interface="in" units="dimensionless"/>
      <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>PIP2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>J_ip3R</ci>
                  <ci>IPX</ci>
               </apply>
               <apply>
                  <divide/>
                  <ci>J_ip3P</ci>
                  <ci>PIP2_Total</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="IPX">
      <variable name="IPX" public_interface="out" units="dimensionless"/>
      <variable name="PIP2" public_interface="in" units="dimensionless"/>
      <variable name="IP3" public_interface="in" units="microM"/>
      <variable initial_value="25" name="PIP2_Total" units="microM"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>IPX</ci>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>IP3</ci>
                     <ci>PIP2_Total</ci>
                  </apply>
               </apply>
               <ci>PIP2</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="J_ip3R">
      <variable cmeta:id="J_ip3R_J_ip3R" name="J_ip3R" public_interface="out" units="per_sec"/>
      <variable initial_value="10" name="J_ip3R0" units="per_sec"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_ip3R</ci>
            <ci>J_ip3R0</ci>
         </apply>
      </math>
   </component>
   <component name="receptor">
      <variable cmeta:id="receptor_ATP_e" initial_value="0" name="ATP_e" public_interface="out" units="microM"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="V_IP3" public_interface="out" units="microMpersec"/>
      <variable initial_value="0" name="MechanicalStimulation" units="dimensionless"/>
      <variable name="V_mech" public_interface="out" units="dimensionless"/>
      <variable initial_value="1" name="V_mech0" units="dimensionless"/>
      <variable initial_value="2.6" name="C_1" public_interface="out" units="microMpersec"/>
      <variable initial_value="0.02" name="C_2" public_interface="out" units="microM"/>
      <variable initial_value="2.3" name="C_3" public_interface="out" units="microMpersec"/>
      <variable initial_value="10" name="C_4" public_interface="out" units="microM"/>
      <variable initial_value="0.001" name="V_ATP" units="microMpersec"/>
      <variable initial_value="0.0001" name="K_ATP" units="microM"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_mech</ci>
            <piecewise>
               <piece>
                  <ci>V_mech0</ci>
                  <apply>
                     <and/>
                     <apply>
                        <gt/>
                        <ci>time</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">10</cn>
                     </apply>
                     <apply>
                        <lt/>
                        <ci>time</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">25</cn>
                     </apply>
                     <apply>
                        <gt/>
                        <ci>MechanicalStimulation</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>V_IP3</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>C_1</ci>
                     <ci>ATP_e</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>C_2</ci>
                     <ci>ATP_e</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>C_3</ci>
                     <apply>
                        <power/>
                        <ci>ATP_e</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <power/>
                        <ci>C_4</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     </apply>
                     <apply>
                        <power/>
                        <ci>ATP_e</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ATP_e</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>V_ATP</ci>
                  </apply>
                  <ci>ATP_e</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_ATP</ci>
                  <ci>ATP_e</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="J_ip3P">
      <variable name="J_ip3P" public_interface="out" units="microMpersec"/>
      <variable name="V_IP3" public_interface="in" units="microMpersec"/>
      <variable name="PIP2" public_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_ip3P</ci>
            <apply>
               <times/>
               <ci>V_IP3</ci>
               <ci>PIP2</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="J_ip3D">
      <variable name="J_ip3D" public_interface="out" units="microMpersec"/>
      <variable name="IP3" public_interface="in" units="microM"/>
      <variable name="volume_i0" public_interface="in" units="volume"/>
      <variable initial_value="0" name="Beta_1" units="microMpersec"/>
      <variable initial_value="36" name="Beta_2" units="microMpersec"/>
      <variable name="Ca_i" public_interface="in" units="microM"/>
      <variable name="F_Ca" units="dimensionless"/>
      <variable initial_value="0.3" name="K_rc" units="microM"/>
      <variable initial_value="0.01" name="K_IP3" units="microM"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>F_Ca</ci>
            <apply>
               <divide/>
               <ci>Ca_i</ci>
               <apply>
                  <plus/>
                  <ci>K_rc</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_ip3D</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <ci>Beta_1</ci>
                     <apply>
                        <times/>
                        <ci>Beta_2</ci>
                        <apply>
                           <power/>
                           <ci>F_Ca</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <ci>IP3</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_IP3</ci>
                  <ci>IP3</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Jipt">
      <variable name="Jipt" public_interface="out" units="microMvolumepersec"/>
      <variable name="Ca_i" public_interface="in" units="microM"/>
      <variable name="Ca_er" public_interface="in" units="microM"/>
      <variable name="Pipt" units="dimensionless"/>
      <variable name="O" public_interface="in" units="dimensionless"/>
      <variable name="volume_i0" public_interface="in" units="volume"/>
      <variable initial_value="28" name="Jipt0" units="per_sec"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Pipt</ci>
            <apply>
               <power/>
               <ci>O</ci>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">4</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jipt</ci>
            <apply>
               <times/>
               <ci>Jipt0</ci>
               <ci>volume_i0</ci>
               <ci>Pipt</ci>
               <apply>
                  <minus/>
                  <ci>Ca_er</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Jryr">
      <variable initial_value="0.3722419" name="Ka" units="microM"/>
      <variable initial_value="0.636005" name="Kb" units="microM"/>
      <variable initial_value="0.0571428" name="Kc" units="dimensionless"/>
      <variable name="Ca_i" public_interface="in" units="microM"/>
      <variable name="Ca_er" public_interface="in" units="microM"/>
      <variable name="volume_i0" public_interface="in" units="volume"/>
      <variable cmeta:id="Jryr_Jryr" name="Jryr" public_interface="out" units="microMvolumepersec"/>
      <variable initial_value="0" name="Jryr0" units="per_sec"/>
      <variable name="Pryr" units="dimensionless"/>
      <variable name="W_inf" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>W_inf</ci>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Ka</ci>
                        <ci>Ca_i</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">4</cn>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <ci>Kb</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>Kc</ci>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <ci>Kb</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Ka</ci>
                        <ci>Ca_i</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">4</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Pryr</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>W_inf</ci>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>Ka</ci>
                           <ci>Ca_i</ci>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">4</cn>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>Ca_i</ci>
                           <ci>Kb</ci>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>Kc</ci>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <ci>Kb</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Ka</ci>
                        <ci>Ca_i</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">4</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jryr</ci>
            <apply>
               <times/>
               <ci>Jryr0</ci>
               <ci>volume_i0</ci>
               <apply>
                  <minus/>
                  <ci>Ca_er</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Jer">
      <variable cmeta:id="Jer_Jer" name="Jer" public_interface="out" units="microMvolumepersec"/>
      <variable name="Ca_i" public_interface="in" units="microM"/>
      <variable name="Ca_er" public_interface="in" units="microM"/>
      <variable name="volume_i0" public_interface="in" units="volume"/>
      <variable initial_value="0.0035" name="Jer0" units="per_sec"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Jer</ci>
            <apply>
               <times/>
               <ci>Jer0</ci>
               <ci>volume_i0</ci>
               <apply>
                  <minus/>
                  <ci>Ca_er</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Jserca">
      <variable cmeta:id="Jserca_Jserca" name="Jserca" public_interface="out" units="microMvolumepersec"/>
      <variable name="Ca_i" public_interface="in" units="microM"/>
      <variable name="Ca_er" public_interface="in" units="microM"/>
      <variable name="volume_i0" public_interface="in" units="volume"/>
      <variable initial_value="0.09" name="Vserca" units="microM2persec"/>
      <variable initial_value="0.04" name="Kserca" units="microM"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Jserca</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>volume_i0</ci>
                        <ci>Vserca</ci>
                        <apply>
                           <power/>
                           <ci>Ca_i</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <power/>
                           <ci>Kserca</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Ca_i</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
               <ci>Ca_er</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Jin">
      <variable name="Jin" public_interface="out" units="microMvolumepersec"/>
      <variable name="Ca_i" public_interface="in" units="microM"/>
      <variable name="Ca_er" public_interface="in" units="microM"/>
      <variable initial_value="-39" name="Vm" units="milliV"/>
      <variable name="V_Ca" units="milliV"/>
      <variable initial_value="8314" name="R" units="millijoule_per_mole_kelvin"/>
      <variable initial_value="310" name="T" units="kelvin"/>
      <variable initial_value="96485" name="F" units="coulomb_per_mole"/>
      <variable initial_value="2" name="z_Ca" units="dimensionless"/>
      <variable initial_value="1300" name="Ca_ext" units="microM"/>
      <variable initial_value="0.6" name="P_pm_1" units="conductancepervolume"/>
      <variable initial_value="28" name="P_pm_2" units="conductancepervolume"/>
      <variable name="Iin_1" units="current"/>
      <variable name="Iin_2" units="current"/>
      <variable name="Jin_1" units="microMvolumepersec"/>
      <variable name="Jin_2" units="microMvolumepersec"/>
      <variable name="volume_i0" public_interface="in" units="volume"/>
      <variable name="V_mech" public_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>z_Ca</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_ext</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Iin_1</ci>
            <apply>
               <times/>
               <ci>volume_i0</ci>
               <ci>P_pm_1</ci>
               <apply>
                  <minus/>
                  <ci>V_Ca</ci>
                  <ci>Vm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Iin_2</ci>
            <apply>
               <times/>
               <ci>volume_i0</ci>
               <ci>P_pm_2</ci>
               <ci>V_mech</ci>
               <apply>
                  <minus/>
                  <ci>V_Ca</ci>
                  <ci>Vm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jin_1</ci>
            <apply>
               <divide/>
               <ci>Iin_1</ci>
               <apply>
                  <times/>
                  <ci>z_Ca</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jin_2</ci>
            <apply>
               <divide/>
               <ci>Iin_2</ci>
               <apply>
                  <times/>
                  <ci>z_Ca</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jin</ci>
            <apply>
               <plus/>
               <ci>Jin_1</ci>
               <ci>Jin_2</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Jpm">
      <variable cmeta:id="Jpm_Jpm" name="Jpm" public_interface="out" units="microMvolumepersec"/>
      <variable name="Ca_i" public_interface="in" units="microM"/>
      <variable initial_value="0.072" name="Vpm" units="microMpersec"/>
      <variable initial_value="0.6" name="Kpm" units="microM"/>
      <variable name="volume_i0" public_interface="in" units="volume"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Jpm</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Vpm</ci>
                  <ci>volume_i0</ci>
                  <apply>
                     <power/>
                     <ci>Ca_i</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Kpm</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Ca_i</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="O">
      <variable initial_value="-2.42013500345263e-89" name="O" public_interface="out" units="dimensionless"/>
      <variable name="S" units="dimensionless"/>
      <variable name="I1" public_interface="in" units="dimensionless"/>
      <variable name="I2" public_interface="in" units="dimensionless"/>
      <variable name="IP3" public_interface="in" units="micromolar"/>
      <variable name="k1" public_interface="in" units="per_secmicroM"/>
      <variable name="k1_" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k2" public_interface="in" units="first_order_rate_constant"/>
      <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>O</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k1</ci>
                  <ci>IP3</ci>
                  <ci>S</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k1_</ci>
                     <ci>O</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k2</ci>
                     <ci>O</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>S</ci>
            <apply>
               <minus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>O</ci>
                  <ci>I1</ci>
                  <ci>I2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I1">
      <variable initial_value="2.24343309680416e-81" name="I1" public_interface="out" units="dimensionless"/>
      <variable name="O" public_interface="in" units="dimensionless"/>
      <variable name="k2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k3" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k4" public_interface="in" units="first_order_rate_constant"/>
      <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>I1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k2</ci>
                  <ci>O</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <ci>k3</ci>
                     <ci>k4</ci>
                  </apply>
                  <ci>I1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I2">
      <variable initial_value="1.01276768045278e-40" name="I2" public_interface="out" units="dimensionless"/>
      <variable name="I1" public_interface="in" units="dimensionless"/>
      <variable name="k4" public_interface="in" units="first_order_rate_constant"/>
      <variable name="k5" public_interface="in" units="first_order_rate_constant"/>
      <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>I2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k4</ci>
                  <ci>I1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k5</ci>
                  <ci>I2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="constants">
      <variable name="k1" public_interface="out" units="per_secmicroM"/>
      <variable initial_value="0.88" name="k1_" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.5" name="k2" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.5" name="k3" public_interface="out" units="first_order_rate_constant"/>
      <variable name="k4" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.02" name="k5" public_interface="out" units="first_order_rate_constant"/>
      <variable name="IP3" public_interface="in" units="microM"/>
      <variable initial_value="40" name="alpha_1" units="per_secmicroM"/>
      <variable initial_value="0.8" name="beta_1" units="microM"/>
      <variable initial_value="0.06" name="alpha_4" units="first_order_rate_constant"/>
      <variable initial_value="0.01" name="beta_4" units="microM"/>
      <variable name="Ca_i" public_interface="in" units="microM"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>k1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>alpha_1</ci>
                  <apply>
                     <power/>
                     <ci>Ca_i</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>beta_1</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Ca_i</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k4</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>alpha_4</ci>
                  <ci>IP3</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>beta_4</ci>
                  <ci>IP3</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <connection>
      <map_components component_1="receptor" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="receptor" component_2="J_ip3P"/>
      <map_variables variable_1="V_IP3" variable_2="V_IP3"/>
   </connection>
   <connection>
      <map_components component_1="receptor" component_2="Jin"/>
      <map_variables variable_1="V_mech" variable_2="V_mech"/>
   </connection>
   <connection>
      <map_components component_1="Jer" component_2="volume_i"/>
      <map_variables variable_1="volume_i0" variable_2="volume_i0"/>
   </connection>
   <connection>
      <map_components component_1="Jserca" component_2="volume_i"/>
      <map_variables variable_1="volume_i0" variable_2="volume_i0"/>
   </connection>
   <connection>
      <map_components component_1="Jipt" component_2="volume_i"/>
      <map_variables variable_1="volume_i0" variable_2="volume_i0"/>
   </connection>
   <connection>
      <map_components component_1="Jipt" component_2="O"/>
      <map_variables variable_1="O" variable_2="O"/>
   </connection>
   <connection>
      <map_components component_1="Jipt" component_2="Ca_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="Jipt" component_2="Ca_er"/>
      <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
   </connection>
   <connection>
      <map_components component_1="Jryr" component_2="volume_i"/>
      <map_variables variable_1="volume_i0" variable_2="volume_i0"/>
   </connection>
   <connection>
      <map_components component_1="J_ip3D" component_2="volume_i"/>
      <map_variables variable_1="volume_i0" variable_2="volume_i0"/>
   </connection>
   <connection>
      <map_components component_1="Jin" component_2="volume_i"/>
      <map_variables variable_1="volume_i0" variable_2="volume_i0"/>
   </connection>
   <connection>
      <map_components component_1="Jpm" component_2="volume_i"/>
      <map_variables variable_1="volume_i0" variable_2="volume_i0"/>
   </connection>
   <connection>
      <map_components component_1="IP3" component_2="O"/>
      <map_variables variable_1="IP3" variable_2="IP3"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_i" component_2="Jipt"/>
      <map_variables variable_1="Jipt" variable_2="Jipt"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_i" component_2="Jryr"/>
      <map_variables variable_1="Jryr" variable_2="Jryr"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_i" component_2="Jer"/>
      <map_variables variable_1="Jer" variable_2="Jer"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_i" component_2="Jserca"/>
      <map_variables variable_1="Jserca" variable_2="Jserca"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_i" component_2="Jpm"/>
      <map_variables variable_1="Jpm" variable_2="Jpm"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_i" component_2="Jin"/>
      <map_variables variable_1="Jin" variable_2="Jin"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_er" component_2="Jipt"/>
      <map_variables variable_1="Jipt" variable_2="Jipt"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_er" component_2="Jryr"/>
      <map_variables variable_1="Jryr" variable_2="Jryr"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_er" component_2="Jer"/>
      <map_variables variable_1="Jer" variable_2="Jer"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_er" component_2="Jserca"/>
      <map_variables variable_1="Jserca" variable_2="Jserca"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_i" component_2="volume_i"/>
      <map_variables variable_1="volume_i" variable_2="volume_i"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_i" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_i" component_2="Ca_i"/>
      <map_variables variable_1="volumeCa_i" variable_2="volumeCa_i"/>
   </connection>
   <connection>
      <map_components component_1="volume_i" component_2="Ca_i"/>
      <map_variables variable_1="volume_i" variable_2="volume_i"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_er" component_2="volume_er"/>
      <map_variables variable_1="volume_er" variable_2="volume_er"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_er" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="O" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="volumeCa_er" component_2="Ca_er"/>
      <map_variables variable_1="volumeCa_er" variable_2="volumeCa_er"/>
   </connection>
   <connection>
      <map_components component_1="volume_er" component_2="Ca_er"/>
      <map_variables variable_1="volume_er" variable_2="volume_er"/>
   </connection>
   <connection>
      <map_components component_1="volume_er" component_2="volume_i"/>
      <map_variables variable_1="volume_i" variable_2="volume_i"/>
   </connection>
   <connection>
      <map_components component_1="volume_iIP3" component_2="J_ip3P"/>
      <map_variables variable_1="J_ip3P" variable_2="J_ip3P"/>
   </connection>
   <connection>
      <map_components component_1="volume_iIP3" component_2="J_ip3D"/>
      <map_variables variable_1="J_ip3D" variable_2="J_ip3D"/>
   </connection>
   <connection>
      <map_components component_1="volume_iIP3" component_2="IP3"/>
      <map_variables variable_1="volume_iIP3" variable_2="volume_iIP3"/>
   </connection>
   <connection>
      <map_components component_1="volume_i" component_2="IP3"/>
      <map_variables variable_1="volume_i" variable_2="volume_i"/>
   </connection>
   <connection>
      <map_components component_1="volume_iIP3" component_2="volume_i"/>
      <map_variables variable_1="volume_i" variable_2="volume_i"/>
   </connection>
   <connection>
      <map_components component_1="volume_iIP3" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="PIP2" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="PIP2" component_2="J_ip3P"/>
      <map_variables variable_1="J_ip3P" variable_2="J_ip3P"/>
      <map_variables variable_1="PIP2" variable_2="PIP2"/>
   </connection>
   <connection>
      <map_components component_1="PIP2" component_2="IPX"/>
      <map_variables variable_1="IPX" variable_2="IPX"/>
      <map_variables variable_1="PIP2" variable_2="PIP2"/>
   </connection>
   <connection>
      <map_components component_1="PIP2" component_2="J_ip3R"/>
      <map_variables variable_1="J_ip3R" variable_2="J_ip3R"/>
   </connection>
   <connection>
      <map_components component_1="IPX" component_2="IP3"/>
      <map_variables variable_1="IP3" variable_2="IP3"/>
   </connection>
   <connection>
      <map_components component_1="J_ip3D" component_2="Ca_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="J_ip3D" component_2="IP3"/>
      <map_variables variable_1="IP3" variable_2="IP3"/>
   </connection>
   <connection>
      <map_components component_1="Jryr" component_2="Ca_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="Jryr" component_2="Ca_er"/>
      <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
   </connection>
   <connection>
      <map_components component_1="Jer" component_2="Ca_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="Jer" component_2="Ca_er"/>
      <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
   </connection>
   <connection>
      <map_components component_1="Jserca" component_2="Ca_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="Jserca" component_2="Ca_er"/>
      <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
   </connection>
   <connection>
      <map_components component_1="Jin" component_2="Ca_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="Jin" component_2="Ca_er"/>
      <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
   </connection>
   <connection>
      <map_components component_1="Jpm" component_2="Ca_i"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="IP3" component_2="constants"/>
      <map_variables variable_1="IP3" variable_2="IP3"/>
   </connection>
   <connection>
      <map_components component_1="I1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="I2" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="O" component_2="I1"/>
      <map_variables variable_1="O" variable_2="O"/>
      <map_variables variable_1="I1" variable_2="I1"/>
   </connection>
   <connection>
      <map_components component_1="O" component_2="I2"/>
      <map_variables variable_1="I2" variable_2="I2"/>
   </connection>
   <connection>
      <map_components component_1="O" component_2="constants"/>
      <map_variables variable_1="k1" variable_2="k1"/>
      <map_variables variable_1="k1_" variable_2="k1_"/>
      <map_variables variable_1="k2" variable_2="k2"/>
   </connection>
   <connection>
      <map_components component_1="I1" component_2="constants"/>
      <map_variables variable_1="k2" variable_2="k2"/>
      <map_variables variable_1="k3" variable_2="k3"/>
      <map_variables variable_1="k4" variable_2="k4"/>
   </connection>
   <connection>
      <map_components component_1="I2" component_2="I1"/>
      <map_variables variable_1="I1" variable_2="I1"/>
   </connection>
   <connection>
      <map_components component_1="I2" component_2="constants"/>
      <map_variables variable_1="k5" variable_2="k5"/>
      <map_variables variable_1="k4" variable_2="k4"/>
   </connection>
   <connection>
      <map_components component_1="Ca_i" component_2="constants"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   
   <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
    
	<rdf:Description rdf:about="">
      <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
        <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
          <vCard:Family>Warren</vCard:Family>
          <vCard:Given>Nic</vCard:Given>
		  <vCard:Other>J</vCard:Other>
        </vCard:N>
        <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
          <rdf:value>n.warren@auckland@auckland.ac.nz</rdf:value>
          <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
        </vCard:EMAIL>
        <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
          <vCard:Orgname>The University of Auckland</vCard:Orgname>
          <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
        </vCard:ORG>
      </dc:creator>
      
       <dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
        <dcterms:W3CDTF>2009-12-10</dcterms:W3CDTF>
      </dcterms:created>
     
    </rdf:Description>
    
    <rdf:Description rdf:about="#warren_2010">
      
	  <cmeta:species>mammalia</cmeta:species>
	  
      <cmeta:bio_entity>lung</cmeta:bio_entity>
	  <cmeta:bio_entity>epithelial cell</cmeta:bio_entity>
      
      <!--  Keyword(s) -->
      <bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
        <dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>electrophysiology</rdf:li>
			  <rdf:li>lung</rdf:li>
			  <rdf:li>epithelia</rdf:li>
			  <rdf:li>calcium dynamics</rdf:li>
			  <rdf:li>signal transduction</rdf:li>
			  <rdf:li>mechanotransduction</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>

      <bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
        <bqs:Pubmed_id>19700517</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/">
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Warren</vCard:Family>
                    <vCard:Given>N</vCard:Given>
					<vCard:Other>J</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Tawhai</vCard:Family>
                    <vCard:Given>M</vCard:Given>
					<vCard:Other>H</vCard:Other>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Crampin</vCard:Family>
                    <vCard:Given>E</vCard:Given>
					<vCard:Other>J</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
            Mathematical modelling of calcium wave propagation in mammalian airway epithelium: evidence for regenerative ATP release
          </dc:title>
          <dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
            <dcterms:W3CDTF>2010-01</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Experimental Physiology</dc:title>
          </bqs:Journal>
		  <bqs:volume>95</bqs:volume>
          <bqs:first_page>232</bqs:first_page>
          <bqs:last_page>249</bqs:last_page>
        </bqs:JournalArticle>
      </bqs:reference>
    </rdf:Description>
  </rdf:RDF>
</model>