Location: Dokos, Celler, Lovell, 1996 @ 58a19cab7f0b / dokos_celler_lovell_1996.cellml

Author:
Dougal Cowan <devnull@localhost>
Date:
2010-05-07 11:35:44+12:00
Desc:
fixed a highlight path
Permanent Source URI:
https://models.cellml.org/workspace/dokos_celler_lovell_1996/rawfile/58a19cab7f0b07a8fee3e0d8d8cca1059b8909b9/dokos_celler_lovell_1996.cellml

<?xml version="1.0"?>
<!--
This CellML file was generated on 29/10/2007 at 18:25:00 using:

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

CellML 1.0 was used to generate this cellular 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="dokos_1996" cmeta:id="dokos_1996" 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>A Single Cell Mathematical Model of the Ion Currents Underlying Sinoatrial Node Pacemaker Activity</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 COR and OpenCell to recreate the published results. The units have been checked and they are consistent.
    </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: The ionic currents underlying autorhythmicity of the mammalian sinoatrial node and their wider contribution to each phase of the action potential have been investigated in this study using a new single cell mathematical model. The new model provides a review and update of existing formulations of sinoatrial node membrane currents, derived from a wide range of electrophysiological data available in the literature. Simulations of spontaneous activity suggest that the dominant mechanism underlying pacemaker depolarisation is the inward background Na+ current, ib,Na. In contrast to previous models, the decay of the delayed rectifying K+ current, iK, was insignificant during this phase. Despite the presence of a pseudo-outward background current throughout the pacemaker range of potentials (Na-K pump+leak currents), the hyperpolarisation-activated current i(f) was not essential to pacemaker activity. A closer inspection of the current-voltage characteristics of the model revealed that the "instantaneous" time-independent current was inward for holding potentials in the pacemaker range, which rapidly became outward within 2 ms due to the inactivation of the L-type Ca2+ current, iCa,L. This suggests that reports in the literature in which the net background current is outward throughout the pacemaker range of potentials may be exaggerated. The magnitudes of the action potential overshoot and the maximum diastolic potential were determined largely by the reversal potentials of iCa,L and iK respectively. The action potential was sustained by the incomplete deactivation of iCa,L and the Na-Ca exchanger, iNaCa. Despite the incorporation of "square-root" activation by [K]o of all K+ currents, the model was unable to correctly simulate the response to elevated [K]o.</para>

<para>
The original paper reference is cited below:
</para>

<para>
Ion Currents Underlying Sinoatrial Node Pacemaker Activity: A New Single Cell Mathematical Model, Socrates Dokos, Branko Celler, and Nigel Lovell, 1996, 
            <emphasis>Theoretical Journal of Biology</emphasis>
          , 181, 245-272.  <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/8869126">PubMed ID: 8869126</ulink> 

</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram</title>
    </objectinfo>
    <imagedata fileref="dokos_1996.png"/>
  </imageobject>

</mediaobject>
<caption>A schematic diagram of the Dokos <emphasis>et al.</emphasis> 1996 mathematical model of the SA node cell.  Sodium, calcium and potassium ions are exchanged between the intracellular and extracellular environments through channels, the sodium-potassium pump and the sodium-calcium exchanger.  Calcium is transferred between the cytosol and the sarcoplasmic reticulum (SR), and between the local regions of the SR.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
   <units name="per_second">
      <unit exponent="-1" units="second"/>
   </units>
   <units name="per_millimolar_second">
      <unit exponent="-1" units="second"/>
      <unit exponent="-1" units="millimolar"/>
   </units>
   <units name="millivolt">
      <unit prefix="milli" units="volt"/>
   </units>
   <units name="per_millivolt">
      <unit exponent="-1" prefix="milli" units="volt"/>
   </units>
   <units name="per_millivolt_second">
      <unit exponent="-1" units="millivolt"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="nanoS">
      <unit prefix="nano" units="siemens"/>
   </units>
   <units name="nanoF">
      <unit prefix="nano" units="farad"/>
   </units>
   <units name="picoA_per_millimolar">
      <unit prefix="pico" units="ampere"/>
      <unit exponent="-1" units="millimolar"/>
   </units>
   <units name="picoA">
      <unit prefix="pico" units="ampere"/>
   </units>
   <units name="mA_pA">
      <unit prefix="milli" units="ampere"/>
      <unit exponent="-1" prefix="pico" units="ampere"/>
   </units>
   <units name="millimolar">
      <unit prefix="milli" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="joule_per_kilomole_kelvin">
      <unit units="joule"/>
      <unit exponent="-1" prefix="kilo" units="mole"/>
      <unit exponent="-1" units="kelvin"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit exponent="-1" units="mole"/>
   </units>
   <units name="microlitre">
      <unit prefix="micro" units="litre"/>
   </units>
   <units name="litre_micrometre3">
      <unit units="litre"/>
      <unit exponent="-1" units="microlitre"/>
   </units>
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="second"/>
   </component>
   <component name="constants">
      <variable initial_value="8314.472" name="R" public_interface="out" units="joule_per_kilomole_kelvin"/>
      <variable initial_value="310" name="T" public_interface="out" units="kelvin"/>
      <variable initial_value="96485.3415" name="F" public_interface="out" units="coulomb_per_mole"/>
   </component>
   <component name="membrane">
      <variable cmeta:id="membrane_E" initial_value="-64.9" name="E" public_interface="out" units="millivolt"/>
      <variable initial_value="3.2e-5" name="C" units="nanoF"/>
      <variable name="i_tot" units="picoA"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="i_CaL" public_interface="in" units="picoA"/>
      <variable name="i_CaT" public_interface="in" units="picoA"/>
      <variable name="i_Na" public_interface="in" units="picoA"/>
      <variable name="i_K" public_interface="in" units="picoA"/>
      <variable name="i_f" public_interface="in" units="picoA"/>
      <variable name="i_p" public_interface="in" units="picoA"/>
      <variable name="i_NaCa" public_interface="in" units="picoA"/>
      <variable name="i_bNa" public_interface="in" units="picoA"/>
      <variable name="i_bK" public_interface="in" units="picoA"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_tot</ci>
            <apply>
               <plus/>
               <ci>i_CaL</ci>
               <ci>i_CaT</ci>
               <ci>i_Na</ci>
               <ci>i_K</ci>
               <ci>i_f</ci>
               <ci>i_p</ci>
               <ci>i_NaCa</ci>
               <ci>i_bNa</ci>
               <ci>i_bK</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>E</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>i_tot</ci>
               </apply>
               <ci>C</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="reversal_potentials">
      <variable name="E_Ca" public_interface="out" units="millivolt"/>
      <variable name="E_Na" public_interface="out" units="millivolt"/>
      <variable name="E_K" public_interface="out" units="millivolt"/>
      <variable name="R" public_interface="in" units="joule_per_kilomole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="Cai" public_interface="in" units="millimolar"/>
      <variable name="Cao" public_interface="in" units="millimolar"/>
      <variable name="Nai" public_interface="in" units="millimolar"/>
      <variable name="Nao" public_interface="in" units="millimolar"/>
      <variable name="Ki" public_interface="in" units="millimolar"/>
      <variable name="Ko" public_interface="in" units="millimolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Nao</ci>
                     <ci>Nai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <ci>Ki</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_current">
      <variable cmeta:id="L_type_calcium_current_i_CaL" name="i_CaL" public_interface="out" units="picoA"/>
      <variable initial_value="0.4" name="g_CaL" units="nanoS"/>
      <variable name="E_Ca" public_interface="in" units="millivolt"/>
      <variable name="time" private_interface="out" public_interface="in" units="second"/>
      <variable name="E" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="Cai" private_interface="out" public_interface="in" units="millimolar"/>
      <variable name="dL" private_interface="in" units="dimensionless"/>
      <variable name="fL" private_interface="in" units="dimensionless"/>
      <variable name="fL2" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_CaL</ci>
            <apply>
               <times/>
               <ci>g_CaL</ci>
               <ci>dL</ci>
               <ci>fL</ci>
               <ci>fL2</ci>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <ci>E</ci>
                     <ci>E_Ca</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">75</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_current_d_gate">
      <variable initial_value="0.0001" name="dL" public_interface="out" units="dimensionless"/>
      <variable name="dL_infinity" units="dimensionless"/>
      <variable initial_value="0.002" name="tau_dL" units="second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>dL_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>E</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">6.6</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">6.6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>dL</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>dL_infinity</ci>
                  <ci>dL</ci>
               </apply>
               <ci>tau_dL</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_current_f_gate">
      <variable initial_value="0.1505" name="fL" public_interface="out" units="dimensionless"/>
      <variable name="fL_infinity" units="dimensionless"/>
      <variable name="tau_fL" units="second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>fL_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>E</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">25</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_fL</ci>
            <apply>
               <plus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">0.031</cn>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>E</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">37.6</cn>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">8.1</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>fL</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>fL_infinity</ci>
                  <ci>fL</ci>
               </apply>
               <ci>tau_fL</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_current_f2_gate">
      <variable initial_value="0.219" name="fL2" public_interface="out" units="dimensionless"/>
      <variable initial_value="3" name="alpha_fL2" units="per_second"/>
      <variable initial_value="40000" name="beta_fL2" units="per_millimolar_second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="Cai" public_interface="in" units="millimolar"/>
      <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>fL2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_fL2</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>fL2</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_fL2</ci>
                  <ci>Cai</ci>
                  <ci>fL2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_calcium_current">
      <variable cmeta:id="T_type_calcium_current_i_CaT" name="i_CaT" public_interface="out" units="picoA"/>
      <variable initial_value="0.085" name="g_CaT" units="nanoS"/>
      <variable name="E_Ca" public_interface="in" units="millivolt"/>
      <variable name="time" private_interface="out" public_interface="in" units="second"/>
      <variable name="E" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="dT" private_interface="in" units="dimensionless"/>
      <variable name="fT" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_CaT</ci>
            <apply>
               <times/>
               <ci>g_CaT</ci>
               <ci>dT</ci>
               <ci>fT</ci>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <ci>E</ci>
                     <ci>E_Ca</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">75</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_calcium_current_d_gate">
      <variable initial_value="0.001" name="dT" public_interface="out" units="dimensionless"/>
      <variable name="dT_infinity" units="dimensionless"/>
      <variable name="tau_dT" units="second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>dT_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>E</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">23</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">6.1</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_dT</ci>
            <apply>
               <plus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">0.0006</cn>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">0.0054</cn>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.03</cn>
                           <apply>
                              <plus/>
                              <ci>E</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">100</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>dT</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>dT_infinity</ci>
                  <ci>dT</ci>
               </apply>
               <ci>tau_dT</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_calcium_current_f_gate">
      <variable initial_value="0.1328" name="fT" public_interface="out" units="dimensionless"/>
      <variable name="fT_infinity" units="dimensionless"/>
      <variable name="tau_fT" units="second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>fT_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>E</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">75</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">6.6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_fT</ci>
            <apply>
               <plus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">0.001</cn>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">0.04</cn>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.08</cn>
                           <apply>
                              <plus/>
                              <ci>E</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">65</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>fT</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>fT_infinity</ci>
                  <ci>fT</ci>
               </apply>
               <ci>tau_fT</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable cmeta:id="fast_sodium_current_i_Na" name="i_Na" public_interface="out" units="picoA"/>
      <variable initial_value="0.25" name="g_Na" units="nanoS"/>
      <variable name="E_Na" public_interface="in" units="millivolt"/>
      <variable name="time" private_interface="out" public_interface="in" units="second"/>
      <variable name="E" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="m" private_interface="in" units="dimensionless"/>
      <variable name="h" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>E</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable initial_value="0.0139" name="m" public_interface="out" units="dimensionless"/>
      <variable name="alpha_m" units="per_second"/>
      <variable name="beta_m" units="per_second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt_second">200</cn>
                  <apply>
                     <plus/>
                     <ci>E</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">34.3</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.09</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>E</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">34.3</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">8000</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.15</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>E</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">56.2</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_m</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>m</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_m</ci>
                  <ci>m</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h_gate">
      <variable initial_value="0.0087" name="h" public_interface="out" units="dimensionless"/>
      <variable name="alpha_h" units="per_second"/>
      <variable name="beta_h" units="per_second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">32.4</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.14</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>E</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">93.4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">709</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">4.2</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.06</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>E</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">45.4</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_h</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>h</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_h</ci>
                  <ci>h</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="delayed_rectifying_potassium_current">
      <variable cmeta:id="delayed_rectifying_potassium_current_i_K" name="i_K" public_interface="out" units="picoA"/>
      <variable name="i_KK" public_interface="out" units="picoA"/>
      <variable name="i_KNa" public_interface="out" units="picoA"/>
      <variable initial_value="0.00026" name="Kk" units="picoA_per_millimolar"/>
      <variable initial_value="0.035" name="P_KNa" units="dimensionless"/>
      <variable name="R" public_interface="in" units="joule_per_kilomole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="Ko" public_interface="in" units="millimolar"/>
      <variable name="Ki" public_interface="in" units="millimolar"/>
      <variable name="Nao" public_interface="in" units="millimolar"/>
      <variable name="Nai" public_interface="in" units="millimolar"/>
      <variable name="time" private_interface="out" public_interface="in" units="second"/>
      <variable name="E" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="x" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_KK</ci>
            <apply>
               <times/>
               <ci>x</ci>
               <ci>Kk</ci>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">1</cn>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.59</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>Ki</ci>
                  <apply>
                     <times/>
                     <ci>Ko</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>E</ci>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_KNa</ci>
            <apply>
               <times/>
               <ci>x</ci>
               <ci>Kk</ci>
               <ci>P_KNa</ci>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">1</cn>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.59</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>Nai</ci>
                  <apply>
                     <times/>
                     <ci>Nao</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>E</ci>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K</ci>
            <apply>
               <plus/>
               <ci>i_KK</ci>
               <ci>i_KNa</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="delayed_rectifying_potassium_current_x_gate">
      <variable initial_value="0.5682" name="x" public_interface="out" units="dimensionless"/>
      <variable name="x_infinity" units="dimensionless"/>
      <variable name="tau_x" units="second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>x_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>E</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">25.1</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">7.4</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_x</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">17</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.0398</cn>
                           <ci>E</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.211</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.051</cn>
                           </apply>
                           <ci>E</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>x_infinity</ci>
                  <ci>x</ci>
               </apply>
               <ci>tau_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarising_activated_current">
      <variable cmeta:id="hyperpolarising_activated_current_i_f" name="i_f" public_interface="out" units="picoA"/>
      <variable name="i_fNa" public_interface="out" units="picoA"/>
      <variable name="i_fK" public_interface="out" units="picoA"/>
      <variable initial_value="10.3" name="Kmf" units="millimolar"/>
      <variable initial_value="0.0081" name="g_fNa" units="nanoS"/>
      <variable initial_value="0.0135" name="g_fK" units="nanoS"/>
      <variable name="E_Na" public_interface="in" units="millivolt"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
      <variable name="Ko" public_interface="in" units="millimolar"/>
      <variable name="time" private_interface="out" public_interface="in" units="second"/>
      <variable name="E" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="y" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_fNa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>y</ci>
                     <apply>
                        <power/>
                        <ci>Ko</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.83</cn>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <power/>
                        <ci>Ko</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.83</cn>
                     </apply>
                     <apply>
                        <power/>
                        <ci>Kmf</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.83</cn>
                     </apply>
                  </apply>
               </apply>
               <ci>g_fNa</ci>
               <apply>
                  <minus/>
                  <ci>E</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_fK</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>y</ci>
                     <apply>
                        <power/>
                        <ci>Ko</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.83</cn>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <power/>
                        <ci>Ko</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.83</cn>
                     </apply>
                     <apply>
                        <power/>
                        <ci>Kmf</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.83</cn>
                     </apply>
                  </apply>
               </apply>
               <ci>g_fK</ci>
               <apply>
                  <minus/>
                  <ci>E</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_f</ci>
            <apply>
               <plus/>
               <ci>i_fK</ci>
               <ci>i_fNa</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarising_activated_current_y_gate">
      <variable initial_value="0.0287" name="y" public_interface="out" units="dimensionless"/>
      <variable name="alpha_y" units="per_second"/>
      <variable name="beta_y" units="per_second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_y</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt_second">0.36</cn>
                  <apply>
                     <plus/>
                     <ci>E</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">137.8</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.066</cn>
                        <apply>
                           <plus/>
                           <ci>E</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">137.8</cn>
                        </apply>
                     </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_y</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt_second">0.1</cn>
                  <apply>
                     <plus/>
                     <ci>E</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">76.3</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.21</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>E</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">76.3</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>y</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_y</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>y</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_y</ci>
                  <ci>y</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable cmeta:id="sodium_potassium_pump_i_p" name="i_p" public_interface="out" units="picoA"/>
      <variable initial_value="40" name="KmNa" units="millimolar"/>
      <variable initial_value="1" name="KmK" units="millimolar"/>
      <variable initial_value="0.226" name="i_pmax" units="picoA"/>
      <variable name="Ko" public_interface="in" units="millimolar"/>
      <variable name="Nai" public_interface="in" units="millimolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_p</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>i_pmax</ci>
                           <ci>Nai</ci>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>Nai</ci>
                           <ci>KmNa</ci>
                        </apply>
                     </apply>
                     <ci>Ko</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>Ko</ci>
                     <ci>KmK</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>E</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">211</cn>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_calcium_exchange_current">
      <variable cmeta:id="sodium_calcium_exchange_current_i_NaCa" name="i_NaCa" public_interface="out" units="picoA"/>
      <variable initial_value="4" name="kNaCa" units="picoA"/>
      <variable name="x1" units="dimensionless"/>
      <variable name="x2" units="dimensionless"/>
      <variable name="x3" units="dimensionless"/>
      <variable name="x4" units="dimensionless"/>
      <variable name="k41" units="dimensionless"/>
      <variable name="k34" units="dimensionless"/>
      <variable name="k23" units="dimensionless"/>
      <variable name="k21" units="dimensionless"/>
      <variable name="k32" units="dimensionless"/>
      <variable name="k43" units="dimensionless"/>
      <variable name="k12" units="dimensionless"/>
      <variable name="k14" units="dimensionless"/>
      <variable initial_value="0.1369" name="Qci" units="dimensionless"/>
      <variable initial_value="0.4315" name="Qn" units="dimensionless"/>
      <variable initial_value="0" name="Qco" units="dimensionless"/>
      <variable initial_value="26.44" name="K3ni" units="millimolar"/>
      <variable initial_value="0.0207" name="Kci" units="millimolar"/>
      <variable initial_value="395.3" name="K1ni" units="millimolar"/>
      <variable initial_value="2.289" name="K2ni" units="millimolar"/>
      <variable initial_value="26.44" name="Kcni" units="millimolar"/>
      <variable initial_value="4.663" name="K3no" units="millimolar"/>
      <variable initial_value="1628" name="K1no" units="millimolar"/>
      <variable initial_value="561.4" name="K2no" units="millimolar"/>
      <variable initial_value="3.663" name="Kco" units="millimolar"/>
      <variable name="do" units="dimensionless"/>
      <variable name="di" units="dimensionless"/>
      <variable name="Cao" public_interface="in" units="millimolar"/>
      <variable name="Nao" public_interface="in" units="millimolar"/>
      <variable name="Cai" public_interface="in" units="millimolar"/>
      <variable name="Nai" public_interface="in" units="millimolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="R" public_interface="in" units="joule_per_kilomole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>do</ci>
            <apply>
               <plus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <divide/>
                  <ci>Cao</ci>
                  <ci>Kco</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Kco</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>Qco</ci>
                           <ci>E</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <ci>Nao</ci>
                  <ci>K1no</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <power/>
                     <ci>Nao</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K1no</ci>
                     <ci>K2no</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <power/>
                     <ci>Nao</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K1no</ci>
                     <ci>K2no</ci>
                     <ci>K3no</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k32</ci>
            <apply>
               <exp/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>Qn</ci>
                     <ci>E</ci>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k23</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <apply>
                           <power/>
                           <ci>Nao</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>K1no</ci>
                           <ci>K2no</ci>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <power/>
                           <ci>Nao</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>K1no</ci>
                           <ci>K2no</ci>
                           <ci>K3no</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Qn</ci>
                           </apply>
                           <ci>E</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>do</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k21</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Kco</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Qco</ci>
                           </apply>
                           <ci>E</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>do</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k34</ci>
            <apply>
               <divide/>
               <ci>Nao</ci>
               <apply>
                  <plus/>
                  <ci>K3no</ci>
                  <ci>Nao</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>di</ci>
            <apply>
               <plus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <divide/>
                  <ci>Cai</ci>
                  <ci>Kci</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Cai</ci>
                     <ci>Kci</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Qci</ci>
                           </apply>
                           <ci>E</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>Cai</ci>
                     <ci>Nai</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Kci</ci>
                     <ci>Kcni</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <ci>Nai</ci>
                  <ci>K1ni</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <power/>
                     <ci>Nai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K1ni</ci>
                     <ci>K2ni</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <power/>
                     <ci>Nai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K1ni</ci>
                     <ci>K2ni</ci>
                     <ci>K3ni</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k41</ci>
            <apply>
               <exp/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <ci>Qn</ci>
                     </apply>
                     <ci>E</ci>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k14</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <apply>
                           <power/>
                           <ci>Nai</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>K1ni</ci>
                           <ci>K2ni</ci>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <power/>
                           <ci>Nai</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        </apply>
                        <apply>
                           <times/>
                           <ci>K1ni</ci>
                           <ci>K2ni</ci>
                           <ci>K3ni</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>Qn</ci>
                           <ci>E</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>di</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k12</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Cai</ci>
                     <ci>Kci</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Qci</ci>
                           </apply>
                           <ci>E</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>di</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k43</ci>
            <apply>
               <divide/>
               <ci>Nai</ci>
               <apply>
                  <plus/>
                  <ci>K3ni</ci>
                  <ci>Nai</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>x1</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>k41</ci>
                  <ci>k34</ci>
                  <apply>
                     <plus/>
                     <ci>k23</ci>
                     <ci>k21</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k21</ci>
                  <ci>k32</ci>
                  <apply>
                     <plus/>
                     <ci>k43</ci>
                     <ci>k41</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>x2</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>k32</ci>
                  <ci>k43</ci>
                  <apply>
                     <plus/>
                     <ci>k14</ci>
                     <ci>k12</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k41</ci>
                  <ci>k12</ci>
                  <apply>
                     <plus/>
                     <ci>k34</ci>
                     <ci>k32</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>x3</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>k14</ci>
                  <ci>k43</ci>
                  <apply>
                     <plus/>
                     <ci>k23</ci>
                     <ci>k21</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k12</ci>
                  <ci>k23</ci>
                  <apply>
                     <plus/>
                     <ci>k43</ci>
                     <ci>k41</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>x4</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>k23</ci>
                  <ci>k34</ci>
                  <apply>
                     <plus/>
                     <ci>k14</ci>
                     <ci>k12</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k14</ci>
                  <ci>k21</ci>
                  <apply>
                     <plus/>
                     <ci>k34</ci>
                     <ci>k32</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>kNaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>x2</ci>
                        <ci>k21</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>x1</ci>
                        <ci>k12</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <ci>x1</ci>
                  <ci>x2</ci>
                  <ci>x3</ci>
                  <ci>x4</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="background_sodium_current">
      <variable cmeta:id="background_sodium_current_i_bNa" name="i_bNa" public_interface="out" units="picoA"/>
      <variable initial_value="0.00024" name="g_Nab" units="nanoS"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="E_Na" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_bNa</ci>
            <apply>
               <times/>
               <ci>g_Nab</ci>
               <apply>
                  <minus/>
                  <ci>E</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="background_potassium_current">
      <variable cmeta:id="background_potassium_current_i_bK" name="i_bK" public_interface="out" units="picoA"/>
      <variable initial_value="0.00007" name="KbK" units="picoA_per_millimolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="E" public_interface="in" units="millivolt"/>
      <variable name="R" public_interface="in" units="joule_per_kilomole_kelvin"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="Ki" public_interface="in" units="millimolar"/>
      <variable name="Ko" public_interface="in" units="millimolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_bK</ci>
            <apply>
               <times/>
               <ci>KbK</ci>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">1</cn>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.41</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>Ki</ci>
                  <apply>
                     <times/>
                     <ci>Ko</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>E</ci>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sarcoplasmic_reticulum_kinetics">
      <variable cmeta:id="sarcoplasmic_reticulum_kinetics_i_up" name="i_up" public_interface="out" units="picoA"/>
      <variable cmeta:id="sarcoplasmic_reticulum_kinetics_i_tr" name="i_tr" public_interface="out" units="picoA"/>
      <variable cmeta:id="sarcoplasmic_reticulum_kinetics_i_rel" name="i_rel" public_interface="out" units="picoA"/>
      <variable name="V_i" public_interface="in" units="microlitre"/>
      <variable name="V_rel" public_interface="out" units="microlitre"/>
      <variable name="V_up" public_interface="out" units="microlitre"/>
      <variable initial_value="0.0212" name="i_up_max" units="picoA"/>
      <variable initial_value="0.0005" name="KmCaup" units="millimolar"/>
      <variable initial_value="0.001" name="KmCarel" units="millimolar"/>
      <variable initial_value="0.005" name="tau_rel" units="second"/>
      <variable initial_value="0.4" name="tau_tr" units="second"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="Cai" public_interface="in" units="millimolar"/>
      <variable name="Caup" public_interface="in" units="millimolar"/>
      <variable name="Carel" public_interface="in" units="millimolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_rel</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.006</cn>
               <ci>V_i</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>V_up</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.014</cn>
               <ci>V_i</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_up</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_up_max</ci>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>KmCaup</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_tr</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                     <ci>V_rel</ci>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                     <ci>tau_tr</ci>
                  </apply>
               </apply>
               <ci>Caup</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_rel</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                        <ci>V_rel</ci>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                        <ci>tau_rel</ci>
                     </apply>
                  </apply>
                  <ci>Carel</ci>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Cai</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>KmCarel</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ion_concentrations">
      <variable initial_value="0.000034" name="Cai" public_interface="out" units="millimolar"/>
      <variable initial_value="2.0004" name="Cao" public_interface="out" units="millimolar"/>
      <variable initial_value="0.5832" name="Caup" public_interface="out" units="millimolar"/>
      <variable initial_value="0.1101" name="Carel" public_interface="out" units="millimolar"/>
      <variable initial_value="7.4994" name="Nai" public_interface="out" units="millimolar"/>
      <variable initial_value="139.9929" name="Nao" public_interface="out" units="millimolar"/>
      <variable initial_value="140.0073" name="Ki" public_interface="out" units="millimolar"/>
      <variable initial_value="5.4243" name="Ko" public_interface="out" units="millimolar"/>
      <variable initial_value="2.5e-6" name="V_i" public_interface="out" units="microlitre"/>
      <variable name="V_e" units="microlitre"/>
      <variable initial_value="0.1" name="tau_b" units="second"/>
      <variable initial_value="140" name="Nab" units="millimolar"/>
      <variable initial_value="2" name="Cab" units="millimolar"/>
      <variable initial_value="5.4" name="Kb" units="millimolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="i_bNa" public_interface="in" units="picoA"/>
      <variable name="i_fNa" public_interface="in" units="picoA"/>
      <variable name="i_Na" public_interface="in" units="picoA"/>
      <variable name="i_p" public_interface="in" units="picoA"/>
      <variable name="i_NaCa" public_interface="in" units="picoA"/>
      <variable name="i_KNa" public_interface="in" units="picoA"/>
      <variable name="i_KK" public_interface="in" units="picoA"/>
      <variable name="i_fK" public_interface="in" units="picoA"/>
      <variable name="i_bK" public_interface="in" units="picoA"/>
      <variable name="i_CaL" public_interface="in" units="picoA"/>
      <variable name="i_CaT" public_interface="in" units="picoA"/>
      <variable name="i_up" public_interface="in" units="picoA"/>
      <variable name="i_rel" public_interface="in" units="picoA"/>
      <variable name="i_tr" public_interface="in" units="picoA"/>
      <variable name="V_rel" public_interface="in" units="microlitre"/>
      <variable name="V_up" public_interface="in" units="microlitre"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_e</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.2</cn>
               <ci>V_i</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Nai</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>i_bNa</ci>
                     <ci>i_fNa</ci>
                     <ci>i_Na</ci>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        <ci>i_p</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        <ci>i_NaCa</ci>
                     </apply>
                     <ci>i_KNa</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>F</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                  <ci>V_i</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Nao</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                     <apply>
                        <plus/>
                        <ci>i_bNa</ci>
                        <ci>i_fNa</ci>
                        <ci>i_Na</ci>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                           <ci>i_p</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                           <ci>i_NaCa</ci>
                        </apply>
                        <ci>i_KNa</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>F</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                     <ci>V_e</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>Nab</ci>
                     <ci>Nao</ci>
                  </apply>
                  <ci>tau_b</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ki</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>i_KK</ci>
                           <ci>i_fK</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           <ci>i_p</ci>
                        </apply>
                     </apply>
                     <ci>i_bK</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>F</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                  <ci>V_i</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ko</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                     <apply>
                        <plus/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>i_KK</ci>
                              <ci>i_fK</ci>
                           </apply>
                           <apply>
                              <times/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                              <ci>i_p</ci>
                           </apply>
                        </apply>
                        <ci>i_bK</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>F</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                     <ci>V_e</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>Kb</ci>
                     <ci>Ko</ci>
                  </apply>
                  <ci>tau_b</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cai</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>i_CaL</ci>
                           <ci>i_CaT</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           <ci>i_NaCa</ci>
                        </apply>
                     </apply>
                     <ci>i_up</ci>
                     <apply>
                        <minus/>
                        <ci>i_rel</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>F</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                  <ci>V_i</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cao</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>i_CaL</ci>
                           <ci>i_CaT</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           <ci>i_NaCa</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     <ci>F</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                     <ci>V_e</ci>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>Cab</ci>
                     <ci>Cao</ci>
                  </apply>
                  <ci>tau_b</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Caup</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                  <apply>
                     <minus/>
                     <ci>i_up</ci>
                     <ci>i_tr</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                  <ci>V_up</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Carel</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mA_pA">1</cn>
                  <apply>
                     <minus/>
                     <ci>i_tr</ci>
                     <ci>i_rel</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="litre_micrometre3">1</cn>
                  <ci>V_rel</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="L_type_calcium_current">
            <component_ref component="L_type_calcium_current_d_gate"/>
            <component_ref component="L_type_calcium_current_f_gate"/>
            <component_ref component="L_type_calcium_current_f2_gate"/>
         </component_ref>
         <component_ref component="T_type_calcium_current">
            <component_ref component="T_type_calcium_current_d_gate"/>
            <component_ref component="T_type_calcium_current_f_gate"/>
         </component_ref>
         <component_ref component="fast_sodium_current">
            <component_ref component="fast_sodium_current_m_gate"/>
            <component_ref component="fast_sodium_current_h_gate"/>
         </component_ref>
         <component_ref component="delayed_rectifying_potassium_current">
            <component_ref component="delayed_rectifying_potassium_current_x_gate"/>
         </component_ref>
         <component_ref component="hyperpolarising_activated_current">
            <component_ref component="hyperpolarising_activated_current_y_gate"/>
         </component_ref>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="sodium_calcium_exchange_current"/>
         <component_ref component="background_sodium_current"/>
         <component_ref component="background_potassium_current"/>
         <component_ref component="sarcoplasmic_reticulum_kinetics"/>
         <component_ref component="ion_concentrations"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="L_type_calcium_current">
         <component_ref component="L_type_calcium_current_d_gate"/>
         <component_ref component="L_type_calcium_current_f_gate"/>
         <component_ref component="L_type_calcium_current_f2_gate"/>
      </component_ref>
      <component_ref component="T_type_calcium_current">
         <component_ref component="T_type_calcium_current_d_gate"/>
         <component_ref component="T_type_calcium_current_f_gate"/>
      </component_ref>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_h_gate"/>
      </component_ref>
      <component_ref component="delayed_rectifying_potassium_current">
         <component_ref component="delayed_rectifying_potassium_current_x_gate"/>
      </component_ref>
      <component_ref component="hyperpolarising_activated_current">
         <component_ref component="hyperpolarising_activated_current_y_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="membrane" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="T_type_calcium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="delayed_rectifying_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchange_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="background_sodium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="background_potassium_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="sarcoplasmic_reticulum_kinetics" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="ion_concentrations" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="constants"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="T" variable_2="T"/>
   </connection>
   <connection>
      <map_components component_1="delayed_rectifying_potassium_current" component_2="constants"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="T" variable_2="T"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchange_current" component_2="constants"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="T" variable_2="T"/>
   </connection>
   <connection>
      <map_components component_1="background_potassium_current" component_2="constants"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="T" variable_2="T"/>
   </connection>
   <connection>
      <map_components component_1="sarcoplasmic_reticulum_kinetics" component_2="constants"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="ion_concentrations" component_2="constants"/>
      <map_variables variable_1="F" variable_2="F"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="membrane"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="i_CaL" variable_2="i_CaL"/>
   </connection>
   <connection>
      <map_components component_1="T_type_calcium_current" component_2="membrane"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="i_CaT" variable_2="i_CaT"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="membrane"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
   </connection>
   <connection>
      <map_components component_1="delayed_rectifying_potassium_current" component_2="membrane"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="i_K" variable_2="i_K"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="membrane"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="i_f" variable_2="i_f"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="membrane"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="i_p" variable_2="i_p"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchange_current" component_2="membrane"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_1="background_sodium_current" component_2="membrane"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="i_bNa" variable_2="i_bNa"/>
   </connection>
   <connection>
      <map_components component_1="background_potassium_current" component_2="membrane"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="i_bK" variable_2="i_bK"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="ion_concentrations"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
      <map_variables variable_1="Cao" variable_2="Cao"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="reversal_potentials"/>
      <map_variables variable_1="E_Ca" variable_2="E_Ca"/>
   </connection>
   <connection>
      <map_components component_1="T_type_calcium_current" component_2="reversal_potentials"/>
      <map_variables variable_1="E_Ca" variable_2="E_Ca"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="reversal_potentials"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="reversal_potentials"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="background_sodium_current" component_2="reversal_potentials"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="ion_concentrations"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
      <map_variables variable_1="i_CaL" variable_2="i_CaL"/>
   </connection>
   <connection>
      <map_components component_1="delayed_rectifying_potassium_current" component_2="ion_concentrations"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
      <map_variables variable_1="i_KK" variable_2="i_KK"/>
      <map_variables variable_1="i_KNa" variable_2="i_KNa"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="ion_concentrations"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
      <map_variables variable_1="i_p" variable_2="i_p"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchange_current" component_2="ion_concentrations"/>
      <map_variables variable_1="Nao" variable_2="Nao"/>
      <map_variables variable_1="Nai" variable_2="Nai"/>
      <map_variables variable_1="Cao" variable_2="Cao"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_1="background_potassium_current" component_2="ion_concentrations"/>
      <map_variables variable_1="Ki" variable_2="Ki"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
      <map_variables variable_1="i_bK" variable_2="i_bK"/>
   </connection>
   <connection>
      <map_components component_1="sarcoplasmic_reticulum_kinetics" component_2="ion_concentrations"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
      <map_variables variable_1="Caup" variable_2="Caup"/>
      <map_variables variable_1="Carel" variable_2="Carel"/>
      <map_variables variable_1="i_tr" variable_2="i_tr"/>
      <map_variables variable_1="i_up" variable_2="i_up"/>
      <map_variables variable_1="i_rel" variable_2="i_rel"/>
      <map_variables variable_1="V_up" variable_2="V_up"/>
      <map_variables variable_1="V_rel" variable_2="V_rel"/>
      <map_variables variable_1="V_i" variable_2="V_i"/>
   </connection>
   <connection>
      <map_components component_1="T_type_calcium_current" component_2="ion_concentrations"/>
      <map_variables variable_1="i_CaT" variable_2="i_CaT"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="ion_concentrations"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="ion_concentrations"/>
      <map_variables variable_1="i_fNa" variable_2="i_fNa"/>
      <map_variables variable_1="i_fK" variable_2="i_fK"/>
      <map_variables variable_1="Ko" variable_2="Ko"/>
   </connection>
   <connection>
      <map_components component_1="background_sodium_current" component_2="ion_concentrations"/>
      <map_variables variable_1="i_bNa" variable_2="i_bNa"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="L_type_calcium_current_d_gate"/>
      <map_variables variable_1="dL" variable_2="dL"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="E" variable_2="E"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="L_type_calcium_current_f_gate"/>
      <map_variables variable_1="fL" variable_2="fL"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="E" variable_2="E"/>
   </connection>
   <connection>
      <map_components component_1="L_type_calcium_current" component_2="L_type_calcium_current_f2_gate"/>
      <map_variables variable_1="fL2" variable_2="fL2"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="E" variable_2="E"/>
      <map_variables variable_1="Cai" variable_2="Cai"/>
   </connection>
   <connection>
      <map_components component_1="T_type_calcium_current" component_2="T_type_calcium_current_d_gate"/>
      <map_variables variable_1="dT" variable_2="dT"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="E" variable_2="E"/>
   </connection>
   <connection>
      <map_components component_1="T_type_calcium_current" component_2="T_type_calcium_current_f_gate"/>
      <map_variables variable_1="fT" variable_2="fT"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="E" variable_2="E"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_m_gate"/>
      <map_variables variable_1="m" variable_2="m"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="E" variable_2="E"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_h_gate"/>
      <map_variables variable_1="h" variable_2="h"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="E" variable_2="E"/>
   </connection>
   <connection>
      <map_components component_1="delayed_rectifying_potassium_current" component_2="delayed_rectifying_potassium_current_x_gate"/>
      <map_variables variable_1="x" variable_2="x"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="E" variable_2="E"/>
   </connection>
   <connection>
      <map_components component_1="hyperpolarising_activated_current" component_2="hyperpolarising_activated_current_y_gate"/>
      <map_variables variable_1="y" variable_2="y"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="E" variable_2="E"/>
   </connection>




<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <rdf:Seq rdf:about="rdf:#02c5e67d-7325-49df-8c31-eefc2f5d7393">
    <rdf:li rdf:resource="rdf:#eca5cc77-0860-4b43-b2f1-32bcfa90e38b"/>
    <rdf:li rdf:resource="rdf:#06e3604c-9f6b-4d50-904c-30a9599863e6"/>
    <rdf:li rdf:resource="rdf:#89b5388f-9e91-4149-86c5-a69c743cf225"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#86c035d1-4adc-4e18-bc4c-eb46192720f1">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">N</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lovell</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ace57437-aeea-46fd-910e-93aaa6711e62">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#edd1c5c4-0ad8-4b0f-a37c-92d081a800b5"/>
    <rdf:value>cmeta:id's have been added to several variables to allow creation of a PCEnv session file.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#76e798f0-ee41-408d-abc7-d16d76454ffa"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#acfbbe34-f9dd-4218-90f6-01bdd8000e08">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>unknown</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#89b5388f-9e91-4149-86c5-a69c743cf225">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#86c035d1-4adc-4e18-bc4c-eb46192720f1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#29b2e279-f828-4d28-869a-4bb483ec2507">
    <bqs:Pubmed_id xmlns:bqs="http://www.cellml.org/bqs/1.0#">8869126</bqs:Pubmed_id>
    <bqs:JournalArticle xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#8df3cc09-f780-42e6-be09-2781df6cdff1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#34d6402c-debe-4537-9e31-731d53191d6d">
    <vCard:Orgname xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">The University of Auckland</vCard:Orgname>
    <vCard:Orgunit xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5ccea1c1-0664-4858-b820-2785fcd42d9a">
    <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#98e79f29-dc05-486b-8f82-a4523678164a"/>
    <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#acfbbe34-f9dd-4218-90f6-01bdd8000e08"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#8745b654-170e-45dc-8715-8bf5e193b97e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#75aa486a-f48b-4405-b5e2-dc3a2e627175">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#e80f779c-f090-43b2-b7da-a34a6f2fe1a3"/>
    <rdf:value>This model has been validated by Penny Noble of Oxford University and is known to run in COR and PCEnv. A PCEnv session file is also associated with this version.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7d1b853b-fee1-434d-b811-84fcdfd20b1a">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">1996-08-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f02d9c27-566e-48ad-a14b-4d7146644dc5">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#90a8dba4-4417-4f57-a81f-f68fe07f3bb9"/>
    <rdf:value>added cmeta:id's for some variables to allow referencing by PCEnv</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#20a6f99f-d94a-4f9e-89e5-5fc770b3d9b6"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8df3cc09-f780-42e6-be09-2781df6cdff1">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#02c5e67d-7325-49df-8c31-eefc2f5d7393"/>
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Ion currents underlying sinoatrial node pacemaker activity: a new single cell mathematical model</dc:title>
    <bqs:volume xmlns:bqs="http://www.cellml.org/bqs/1.0#">181(3)</bqs:volume>
    <bqs:first_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">245</bqs:first_page>
    <bqs:Journal xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#1724b578-45c5-4380-8246-fc2dd34bc2b2"/>
    <dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#7d1b853b-fee1-434d-b811-84fcdfd20b1a"/>
    <bqs:last_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">272</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#eca5cc77-0860-4b43-b2f1-32bcfa90e38b">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#f2c2757e-e7ed-4cbb-90a6-8ffaf67f96d9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1724b578-45c5-4380-8246-fc2dd34bc2b2">
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Journal of Theroretical Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a574a2ba-22e8-4c75-a2ec-c37859024a4b">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#261db724-b5a5-42bc-abed-2388a95a6e0a"/>
    <rdf:value>This is the CellML description of Dokos et al's 1996 mathematical 
model of the ion currents underlying sinoatrial node pacemaker 
activity. Please note that equations for reversal potentials were not           included in the original published model and have had to be added.  
Also, a current called i_KACh was included in the equations for 
external and internal potassium ion concentration calculation.  This 
current had not been previously mentioned in the paper, and as i_KACh 
is usually referred to in ischemic tissue, and also the background 
potassium current i_bK had not been included in the equation, I am 
assuming that this was a typo.  The model authors have been contacted 
and we are awaiting their reply.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#edd1c5c4-0ad8-4b0f-a37c-92d081a800b5">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-06-19T16:01:50+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8745b654-170e-45dc-8715-8bf5e193b97e">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher xmlns:dc="http://purl.org/dc/elements/1.1/">The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#75aa486a-f48b-4405-b5e2-dc3a2e627175"/>
    <dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#d54dfe01-5253-4e09-ba29-2c899b8f43c3"/>
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#34e891fc-e1e0-4625-8934-c6f5a8c82ad0"/>
    <cmeta:modification rdf:resource="rdf:#f02d9c27-566e-48ad-a14b-4d7146644dc5"/>
    <cmeta:modification rdf:resource="rdf:#5c2070b9-7348-4845-990d-be21a6a4dd27"/>
    <cmeta:modification rdf:resource="rdf:#6ed5da74-dec8-406f-ba93-911f93fd6f69"/>
    <cmeta:modification rdf:resource="rdf:#ace57437-aeea-46fd-910e-93aaa6711e62"/>
  </rdf:Description>
  <rdf:Description rdf:about="#dokos_1996">
    <ns7:simulation xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" rdf:resource="rdf:#$pk2nE"/>
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">unknown</dc:title>
    <cmeta:comment rdf:resource="rdf:#a574a2ba-22e8-4c75-a2ec-c37859024a4b"/>
    <bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#29b2e279-f828-4d28-869a-4bb483ec2507"/>
<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>cardiac</rdf:li>
        <rdf:li>electrophysiology</rdf:li>
        <rdf:li>Sinoatrial Node</rdf:li>
        <rdf:li>pacemaker</rdf:li>
        <rdf:li>sinoatrial node</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7f774b4b-e9f8-40f9-b84a-9e59c5d426b1">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Penny</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Noble</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#24727d09-7250-4a4d-bbf0-addb6c67ec23">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Penny</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Noble</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e80f779c-f090-43b2-b7da-a34a6f2fe1a3">
    <vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#994bbbe8-76fe-4d5c-9bf0-e3afec3f89b4">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">B</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Celler</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8c05249a-2613-443c-8234-fa601710f804">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lawson</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Richard</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c2070b9-7348-4845-990d-be21a6a4dd27">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#16a51c5e-1c96-4dfc-8c81-48e872f171c6"/>
    <rdf:value>This model has been validated by Penny Noble</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#c0f222ac-f1c1-4aaa-821b-43f5dd20ae0d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fd71e851-3895-44f3-8e03-176c535ca656">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lloyd</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6ed5da74-dec8-406f-ba93-911f93fd6f69">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#50bd5b2a-64b6-4aaf-beeb-4ec42c7b5865"/>
    <rdf:value>Units checked, curated.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#5abcd52a-3b14-407a-9b90-d2d5a0f21f1a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$pk2nE">
    <ns7:boundIntervals xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" rdf:resource="rdf:#$qk2nE"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d54dfe01-5253-4e09-ba29-2c899b8f43c3">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2003-04-10T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#90a8dba4-4417-4f57-a81f-f68fe07f3bb9">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-12-19T14:19:43+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#34e891fc-e1e0-4625-8934-c6f5a8c82ad0">
    <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#34d6402c-debe-4537-9e31-731d53191d6d"/>
    <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#5f4f4c45-bf1b-405d-a10b-abdc66064fe6"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#fd71e851-3895-44f3-8e03-176c535ca656"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$rk2nE">
    <ns7:endingValue xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#">2</ns7:endingValue>
    <ns7:nonstandard-pointDensity xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#">100000</ns7:nonstandard-pointDensity>
    <ns7:maximumStepSize xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#">0.1</ns7:maximumStepSize>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f2c2757e-e7ed-4cbb-90a6-8ffaf67f96d9">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">S</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Dokos</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fc49bcaf-1457-4f79-a984-5a0b7afc1205">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lawson</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Richard</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#16a51c5e-1c96-4dfc-8c81-48e872f171c6">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-06-19T15:55:00+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5f4f4c45-bf1b-405d-a10b-abdc66064fe6">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$qk2nE">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$rk2nE"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#06e3604c-9f6b-4d50-904c-30a9599863e6">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#994bbbe8-76fe-4d5c-9bf0-e3afec3f89b4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#98e79f29-dc05-486b-8f82-a4523678164a">
    <vCard:Orgname xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Orgname>
    <vCard:Orgunit xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5abcd52a-3b14-407a-9b90-d2d5a0f21f1a">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#7f774b4b-e9f8-40f9-b84a-9e59c5d426b1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#76e798f0-ee41-408d-abc7-d16d76454ffa">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#8c05249a-2613-443c-8234-fa601710f804"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#4c91e375-d489-498d-91a9-6ad412eb2cb4"/>
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#5ccea1c1-0664-4858-b820-2785fcd42d9a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#50bd5b2a-64b6-4aaf-beeb-4ec42c7b5865">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-12-19T14:19:12+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#20a6f99f-d94a-4f9e-89e5-5fc770b3d9b6">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#fc49bcaf-1457-4f79-a984-5a0b7afc1205"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4c91e375-d489-498d-91a9-6ad412eb2cb4">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2006-01-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#261db724-b5a5-42bc-abed-2388a95a6e0a">
    <vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c0f222ac-f1c1-4aaa-821b-43f5dd20ae0d">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#24727d09-7250-4a4d-bbf0-addb6c67ec23"/>
  </rdf:Description>
</rdf:RDF>
</model>