Location: Kurata, Hisatome, Imanishi, Shibamoto, 2002 @ 5cc3253726ef / kurata_hisatome_imanishi_shibamoto_2002.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2009-06-17 14:49:49+12:00
Desc:
committing version05 of kurata_hisatome_imanishi_shibamoto_2002
Permanent Source URI:
https://models.physiomeproject.org/workspace/kurata_hisatome_imanishi_shibamoto_2002/rawfile/5cc3253726eff00620b85907cfd5d4e7200f63d9/kurata_hisatome_imanishi_shibamoto_2002.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 07/05/2008 at 17:28:48 using:

COR (0.9.31.927)
Copyright 2002-2008 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="http://www.cellml.org/cellml/1.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:RDF="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:NS1="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:NS2="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:NS3="http://purl.org/dc/elements/1.1/" xmlns:NS4="http://www.cellml.org/metadata/1.0#" xmlns:NS5="http://www.cellml.org/metadata/simulation/1.0#" xmlns:NS6="http://purl.org/dc/terms/" xmlns:NS7="http://imc.org/vCard/3.0#" cmeta:id="kurata_hisatome_imanishi_shibamoto_2006_version02" name="kurata_hisatome_imanishi_shibamoto_2002_version02">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Improved Mathematical Model for the Primary Pacemaker Cell</title>
  <author>
    <firstname>Penny</firstname>
          <surname>Noble</surname>
    <affiliation>
      <shortaffil>Cardiac Electrophysiology Group, Oxford University</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This model was created from Version 01 by Penny Noble of Oxford University. Version 02 is known to run in COR and PCEnv. A PCEnv session is also associated with this model.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
The pacemaker activity of the sinoatrial (SA) node initiates the spontaneous beating of the heart.  Over the past decade a large body of data on the ionic currents underlying this pacemaker activity has been elucidated in patch-clamp experiments.  On the basis of this data, several mathematical models describing the pace maker activity of a single rabbit SA node cell have been developed.  These include:
<itemizedlist>
  <listitem>
              <para>
                <ulink url="${HTML_EXMPL_D_SAN_MODEL}">Demir <emphasis>et al.</emphasis> Sinoatrial Node Model, 1994;</ulink>
              </para>
            </listitem>
  <listitem>
              <para>
                <ulink url="${HTML_EXMPL_DOKOS_MODEL}">Modelling the Ion Currents Underlying Sinoatrial Node Pacemaker Activity, Dokos <emphasis>et al.</emphasis>, 1996;</ulink>
              </para>
            </listitem>
  <listitem>
              <para>
                <ulink url="${HTML_EXMPL_DOKOS_MODEL_II}">A Model of Sinoatrial Node Vagal Control, Dokos <emphasis>et al.</emphasis>, 1996;</ulink>
              </para>
            </listitem>
  <listitem>
              <para>
                <ulink url="${HTML_EXMPL_D99_SAN_MODEL}">Demir <emphasis>et al.</emphasis> Sinoatrial Node Model, 1999;</ulink>
              </para>
            </listitem>
  <listitem>
              <para>
                <ulink url="${HTML_EXMPL_Z_SAN_MODEL}">Zhang <emphasis>et al.</emphasis> Sinoatrial Node Model, 2000; and</ulink>
              </para>
            </listitem>
  <listitem>
              <para>
                <ulink url="${HTML_EXMPL_B_SAN_MODEL}">Boyett <emphasis>et al.</emphasis> Sinoatrial Node Model, 2001.</ulink>
              </para>
            </listitem>
</itemizedlist>
</para>

<para>
With the elucidation of more experimental data, these models have become increasingly complex, with equations that describe calcium buffering, intracellular compartmentation, and additional ionic currents through channels, pumps and exchangers.  In the study by Kurata <emphasis>et al.</emphasis> described here, the aim is to develop an improved mathematical model of a single primary pacemaker cell of the rabbit SA node, in order to investigate the mechanisms underlying pacemaker generation.  On the basis of the most recent experimental data, Kurata <emphasis>et al.</emphasis> updated the previously published models (see list above) in several ways:
</para>
<itemizedlist>
  <listitem>
            <para>A sustained inward current has been added (<emphasis>I<subscript>st</subscript>
              </emphasis>).</para>
          </listitem>
  <listitem>
            <para>The voltage- and Ca<superscript>2+</superscript>-dependent activation kinetics of the L-type calcium channel current (<emphasis>I<subscript>Ca,L</subscript>
              </emphasis>) have been reformulated.</para>
          </listitem>
  <listitem>
            <para>Expressions for the activation kinetics of the rapidly activating delayed rectifier potassium current (<emphasis>I<subscript>Kr</subscript>
              </emphasis>) have been updated.</para>
          </listitem>
  <listitem>
            <para>Revised kinetic equations for two 4-AP sensitive currents (<emphasis>I<subscript>to</subscript>
              </emphasis> and <emphasis>I<subscript>sus</subscript>
              </emphasis>) have been incorporated.</para>
          </listitem>
  <listitem>
            <para>Voltage- and concentration-dependent kinetics of the Na<superscript>+</superscript>-K<superscript>+</superscript> pump current (<emphasis>I<subscript>NaK</subscript>
              </emphasis>) have been reformulated.</para>
          </listitem>
  <listitem>
            <para>The subsarcolemmal space (see <xref linkend="fig_cell_diagram"/> below) as a diffusion barrier for intracellular Ca<superscript>2+</superscript> has been added.</para>
          </listitem>
</itemizedlist>

<para>
The model was validated by comparing simulation results with experimental data, and also by comparing them with those of previous mathematical models.  The authors concluded that their model represented a significant improvement over previous models because it can:
<itemizedlist>
  <listitem>
              <para>Simulate whole cell voltage-clamp data for <emphasis>I<subscript>Ca,L</subscript>
                </emphasis>, <emphasis>I<subscript>Kr</subscript>
                </emphasis>, and <emphasis>I<subscript>st</subscript>
                </emphasis>;</para>
            </listitem>
  <listitem>
              <para>Reproduce the wave shapes of spontaneous action potentials and ionic currents during action potential clamp recordings; and</para>
            </listitem>
  <listitem>
              <para>Mimic the effects of channel blockers or Ca <superscript>2+</superscript> buffers on pacemaker activity more accurately than the previous models.</para>
            </listitem>
</itemizedlist>
</para>

<para>
The model has been described here in CellML (the raw CellML description of the Kurata <emphasis>et al.</emphasis> 2002 model can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>).     
</para>

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

<para>
<ulink url="http://ajpheart.physiology.org/cgi/content/abstract/283/5/H2074">Dynamical description of sinoatrial node pacemaking: improved mathematical model for primary pacemaker cell</ulink>, Yasutaka Kurata, Ichiro Hisatome, Sunao Imanishi, and Toshishige Shibamoto, 2002, <ulink url="http://ajpheart.physiology.org/">
            <emphasis>American Journal of Physiology</emphasis>
          </ulink>, 283, H2074-H2101.  (<ulink url="http://ajpheart.physiology.org/cgi/content/full/283/5/H2074">Full text (HTML)</ulink> and <ulink url="http://ajpheart.physiology.org/cgi/reprint/283/5/H2074.pdf">PDF</ulink> versions of the article are available on the <emphasis>American Journal of Physiology</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=12384487&amp;dopt=Abstract">PubMed ID: 12384487</ulink>
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>reaction diagram</title>
    </objectinfo>
    <imagedata fileref="kurata_2002.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram depicting the intracellular compartments for Ca<superscript>2+</superscript>.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
   <units name="per_millisecond">
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="millimolar_per_millisecond">
      <unit units="millimolar"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="per_millimolar_millisecond">
      <unit units="millimolar" exponent="-1"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="millisecond">
      <unit units="second" prefix="milli"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="per_millivolt">
      <unit units="millivolt" exponent="-1"/>
   </units>
   <units name="per_millivolt_millisecond">
      <unit units="millivolt" exponent="-1"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="nanoS">
      <unit units="siemens" prefix="nano"/>
   </units>
   <units name="nanoS_per_picoF">
      <unit units="siemens" prefix="nano"/>
      <unit units="farad" prefix="pico" exponent="-1"/>
   </units>
   <units name="picoF">
      <unit units="farad" prefix="pico"/>
   </units>
   <units name="per_picoF">
      <unit units="farad" prefix="pico" exponent="-1"/>
   </units>
   <units name="picoA_per_millimolar_picoF">
      <unit units="ampere" prefix="pico"/>
      <unit units="millimolar" exponent="-1"/>
      <unit units="farad" prefix="pico" exponent="-1"/>
   </units>
   <units name="picoA">
      <unit units="ampere" prefix="pico"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="attojoule_per_millimole_kelvin">
      <unit units="joule" prefix="atto"/>
      <unit units="mole" prefix="milli" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="femtocoulomb_per_millimole">
      <unit units="coulomb" prefix="femto"/>
      <unit units="mole" prefix="milli" exponent="-1"/>
   </units>
   <component name="environment">
      <variable units="millisecond" public_interface="out" name="time"/>
   </component>
   <component name="membrane">
      <variable units="millivolt" public_interface="out" name="V" initial_value="-58.600291137693"/>
      <variable units="attojoule_per_millimole_kelvin" public_interface="out" name="R" initial_value="8314400000000000"/>
      <variable units="kelvin" public_interface="out" name="T" initial_value="310.15"/>
      <variable units="femtocoulomb_per_millimole" public_interface="out" name="F" initial_value="96485000000000000"/>
      <variable units="picoF" public_interface="out" name="Cm" initial_value="32"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="picoA" public_interface="in" name="i_CaT"/>
      <variable units="picoA" public_interface="in" name="i_CaL"/>
      <variable units="picoA" public_interface="in" name="i_Kr"/>
      <variable units="picoA" public_interface="in" name="i_Ks"/>
      <variable units="picoA" public_interface="in" name="i_to"/>
      <variable units="picoA" public_interface="in" name="i_sus"/>
      <variable units="picoA" public_interface="in" name="i_h"/>
      <variable units="picoA" public_interface="in" name="i_st"/>
      <variable units="picoA" public_interface="in" name="i_b_Na"/>
      <variable units="picoA" public_interface="in" name="i_NaCa"/>
      <variable units="picoA" public_interface="in" name="i_NaK"/>
      <variable units="picoA" public_interface="in" name="i_K_ACh"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_CaL</ci>
                     <ci>i_CaT</ci>
                     <ci>i_Kr</ci>
                     <ci>i_Ks</ci>
                     <ci>i_to</ci>
                     <ci>i_sus</ci>
                     <ci>i_h</ci>
                     <ci>i_st</ci>
                     <ci>i_b_Na</ci>
                     <ci>i_K_ACh</ci>
                     <ci>i_NaK</ci>
                     <ci>i_NaCa</ci>
                  </apply>
               </apply>
               <cn cellml:units="picoF">1</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="reversal_potentials">
      <variable units="millivolt" public_interface="out" name="E_Na"/>
      <variable units="millivolt" public_interface="out" name="E_K"/>
      <variable units="attojoule_per_millimole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="femtocoulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="millimolar" public_interface="in" name="Ki"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <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_channel_current">
      <variable units="picoA" public_interface="out" name="i_CaL"/>
      <variable units="nanoS" name="g_CaL" initial_value="0.58"/>
      <variable units="millivolt" name="E_CaL" initial_value="45"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millimolar" public_interface="in" private_interface="out" name="Ca_sub"/>
      <variable units="dimensionless" private_interface="in" name="d"/>
      <variable units="dimensionless" private_interface="in" name="f"/>
      <variable units="dimensionless" private_interface="in" name="fCa"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_CaL</ci>
            <apply>
               <times/>
               <ci>g_CaL</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_CaL</ci>
               </apply>
               <ci>d</ci>
               <ci>f</ci>
               <ci>fCa</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_channel_current_d_gate">
      <variable units="dimensionless" public_interface="out" name="d" initial_value="0.000602055134"/>
      <variable units="dimensionless" name="d_infinity"/>
      <variable units="millisecond" name="tau_d"/>
      <variable units="per_millisecond" name="alpha_d"/>
      <variable units="per_millisecond" name="beta_d"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">14.1</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_d</ci>
                  <ci>beta_d</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_d</ci>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt_millisecond">0.02839</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">35</cn>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">35</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">2.5</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt_millisecond">0.0849</cn>
                     <ci>V</ci>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">4.808</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_d</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millivolt_millisecond">0.01143</cn>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">5</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">5</cn>
                        </apply>
                        <cn cellml:units="millivolt">2.5</cn>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_infinity</ci>
                  <ci>d</ci>
               </apply>
               <ci>tau_d</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_channel_current_f_gate">
      <variable units="dimensionless" public_interface="out" name="f" initial_value="0.626999773853"/>
      <variable units="dimensionless" name="f_infinity"/>
      <variable units="millisecond" name="tau_f"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f</ci>
            <apply>
               <plus/>
               <cn cellml:units="millisecond">44.3</cn>
               <apply>
                  <times/>
                  <cn cellml:units="millisecond">257.1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">32.5</cn>
                              </apply>
                              <cn cellml:units="millivolt">13.9</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_infinity</ci>
                  <ci>f</ci>
               </apply>
               <ci>tau_f</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_calcium_channel_current_fCa_gate">
      <variable units="dimensionless" public_interface="out" name="fCa" initial_value="0.589580408056"/>
      <variable units="per_millisecond" name="alpha_fCa"/>
      <variable units="per_millimolar_millisecond" name="beta_fCa" initial_value="60"/>
      <variable units="dimensionless" name="fCa_infinity"/>
      <variable units="millisecond" name="tau_fCa"/>
      <variable units="millimolar" name="Km_fCa" initial_value="0.00035"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millimolar" public_interface="in" name="Ca_sub"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_fCa</ci>
            <apply>
               <times/>
               <ci>Km_fCa</ci>
               <ci>beta_fCa</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>fCa_infinity</ci>
            <apply>
               <divide/>
               <ci>Km_fCa</ci>
               <apply>
                  <plus/>
                  <ci>Km_fCa</ci>
                  <ci>Ca_sub</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_fCa</ci>
            <apply>
               <divide/>
               <ci>fCa_infinity</ci>
               <ci>alpha_fCa</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>fCa</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>fCa_infinity</ci>
                  <ci>fCa</ci>
               </apply>
               <ci>tau_fCa</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_calcium_channel_current">
      <variable units="picoA" public_interface="out" name="i_CaT"/>
      <variable units="nanoS" name="g_CaT" initial_value="0.458"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millivolt" name="E_CaT" initial_value="45"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="d"/>
      <variable units="dimensionless" private_interface="in" name="f"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_CaT</ci>
            <apply>
               <times/>
               <ci>g_CaT</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_CaT</ci>
               </apply>
               <ci>d</ci>
               <ci>f</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_calcium_channel_current_d_gate">
      <variable units="dimensionless" public_interface="out" name="d" initial_value="0.004571884917"/>
      <variable units="dimensionless" name="d_infinity"/>
      <variable units="millisecond" name="tau_d"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>d_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">26.3</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d</ci>
            <apply>
               <divide/>
               <cn cellml:units="millisecond">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">1.068</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">26.3</cn>
                           </apply>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">1.068</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">26.3</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_infinity</ci>
                  <ci>d</ci>
               </apply>
               <ci>tau_d</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="T_type_calcium_channel_current_f_gate">
      <variable units="dimensionless" public_interface="out" name="f" initial_value="0.249637570396"/>
      <variable units="dimensionless" name="f_infinity"/>
      <variable units="millisecond" name="tau_f"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">61.7</cn>
                        </apply>
                        <cn cellml:units="millivolt">5.6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f</ci>
            <apply>
               <divide/>
               <cn cellml:units="millisecond">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0153</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">61.7</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">83.3</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.015</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">61.7</cn>
                           </apply>
                           <cn cellml:units="millivolt">15.38</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_infinity</ci>
                  <ci>f</ci>
               </apply>
               <ci>tau_f</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapidly_activating_delayed_rectifier_potassium_current">
      <variable units="picoA" public_interface="out" name="i_Kr"/>
      <variable units="nanoS" name="g_Kr"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="dimensionless" private_interface="in" name="paS"/>
      <variable units="dimensionless" private_interface="in" name="paF"/>
      <variable units="dimensionless" private_interface="in" name="piy"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_Kr</ci>
            <apply>
               <times/>
               <cn cellml:units="nanoS">0.025</cn>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <cn cellml:units="millimolar">1</cn>
                  </apply>
                  <cn cellml:units="dimensionless">0.59</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <ci>g_Kr</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.6</cn>
                     <ci>paF</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.4</cn>
                     <ci>paS</ci>
                  </apply>
               </apply>
               <ci>piy</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapidly_activating_delayed_rectifier_potassium_current_pa_gate">
      <variable units="dimensionless" public_interface="out" name="paS" initial_value="0.629323128348"/>
      <variable units="dimensionless" public_interface="out" name="paF" initial_value="0.3493633709533"/>
      <variable units="dimensionless" name="pa_infinity"/>
      <variable units="millisecond" name="tau_paS"/>
      <variable units="millisecond" name="tau_paF"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>pa_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">23.2</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">10.6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_paS</ci>
            <apply>
               <divide/>
               <cn cellml:units="millisecond">0.84655354</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0042</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">17</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.00015</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">21.6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_paF</ci>
            <apply>
               <divide/>
               <cn cellml:units="millisecond">0.84655354</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0372</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">15.9</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.00096</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">22.5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>paS</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>pa_infinity</ci>
                  <ci>paS</ci>
               </apply>
               <ci>tau_paS</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>paF</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>pa_infinity</ci>
                  <ci>paF</ci>
               </apply>
               <ci>tau_paF</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapidly_activating_delayed_rectifier_potassium_current_pi_gate">
      <variable units="dimensionless" public_interface="out" name="piy" initial_value="0.852396631172"/>
      <variable units="dimensionless" name="pi_infinity"/>
      <variable units="millisecond" name="tau_pi"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>pi_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">28.6</cn>
                        </apply>
                        <cn cellml:units="millivolt">17.1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_pi</ci>
            <apply>
               <divide/>
               <cn cellml:units="millisecond">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">54.645</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.656</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">106.157</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>piy</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>pi_infinity</ci>
                  <ci>piy</ci>
               </apply>
               <ci>tau_pi</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slowly_activating_delayed_rectifier_potassium_current">
      <variable units="picoA" public_interface="out" name="i_Ks"/>
      <variable units="nanoS" name="g_Ks" initial_value="0.0259"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millivolt" name="E_Ks"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millimolar" public_interface="in" name="Ki"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="attojoule_per_millimole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="femtocoulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="dimensionless" private_interface="in" name="n"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Ks</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>Ko</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>Nao</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Ki</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.12</cn>
                           <ci>Nai</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ks</ci>
            <apply>
               <times/>
               <ci>g_Ks</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ks</ci>
               </apply>
               <apply>
                  <power/>
                  <ci>n</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slowly_activating_delayed_rectifier_potassium_current_n_gate">
      <variable units="dimensionless" public_interface="out" name="n" initial_value="0.054409723782"/>
      <variable units="dimensionless" name="n_infinity"/>
      <variable units="millisecond" name="tau_n"/>
      <variable units="per_millisecond" name="alpha_n"/>
      <variable units="per_millisecond" name="beta_n"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>n_infinity</ci>
            <apply>
               <divide/>
               <ci>alpha_n</ci>
               <apply>
                  <plus/>
                  <ci>alpha_n</ci>
                  <ci>beta_n</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_n</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_n</ci>
                  <ci>beta_n</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_n</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_millisecond">0.014</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">40</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">9</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_n</ci>
            <apply>
               <times/>
               <cn cellml:units="per_millisecond">0.001</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">45</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>n</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>n_infinity</ci>
                  <ci>n</ci>
               </apply>
               <ci>tau_n</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="AP_sensitive_currents">
      <variable units="picoA" public_interface="out" name="i_to"/>
      <variable units="picoA" public_interface="out" name="i_sus"/>
      <variable units="nanoS" name="g_to" initial_value="0.18"/>
      <variable units="nanoS" name="g_sus" initial_value="0.02"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="q"/>
      <variable units="dimensionless" private_interface="in" name="r"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_to</ci>
            <apply>
               <times/>
               <ci>g_to</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
               <ci>q</ci>
               <ci>r</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_sus</ci>
            <apply>
               <times/>
               <ci>g_sus</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
               <ci>r</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="AP_sensitive_currents_q_gate">
      <variable units="dimensionless" public_interface="out" name="q" initial_value="0.531446952485"/>
      <variable units="dimensionless" name="q_infinity"/>
      <variable units="millisecond" name="tau_q"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>q_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">49</cn>
                        </apply>
                        <cn cellml:units="millivolt">13</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_q</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.6</cn>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">65.17</cn>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millisecond">0.57</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="per_millivolt">0.08</cn>
                                 </apply>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">44</cn>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millisecond">0.065</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="per_millivolt">0.1</cn>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">45.93</cn>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="millisecond">10.1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>q</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>q_infinity</ci>
                  <ci>q</ci>
               </apply>
               <ci>tau_q</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="AP_sensitive_currents_r_gate">
      <variable units="dimensionless" public_interface="out" name="r" initial_value="0.005550489445"/>
      <variable units="dimensionless" name="r_infinity"/>
      <variable units="millisecond" name="tau_r"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>r_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">19.3</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">15</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_r</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.66</cn>
               <cn cellml:units="dimensionless">1.4</cn>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <cn cellml:units="millisecond">15.59</cn>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">1.037</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="per_millivolt">0.09</cn>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">30.61</cn>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.369</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="per_millivolt">0.12</cn>
                                 </apply>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">23.84</cn>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="millisecond">2.98</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>r_infinity</ci>
                  <ci>r</ci>
               </apply>
               <ci>tau_r</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarisation_activated_current">
      <variable units="picoA" public_interface="out" name="i_h"/>
      <variable units="picoA" public_interface="out" name="i_h_Na"/>
      <variable units="picoA" public_interface="out" name="i_h_K"/>
      <variable units="nanoS" name="g_h_Na" initial_value="0.1437375"/>
      <variable units="nanoS" name="g_h_K" initial_value="0.2312625"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="y"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_h_Na</ci>
            <apply>
               <times/>
               <ci>g_h_Na</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
               <apply>
                  <power/>
                  <ci>y</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_h_K</ci>
            <apply>
               <times/>
               <ci>g_h_K</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
               <apply>
                  <power/>
                  <ci>y</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_h</ci>
            <apply>
               <plus/>
               <ci>i_h_Na</ci>
               <ci>i_h_K</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarisation_activated_current_y_gate">
      <variable units="dimensionless" public_interface="out" name="y" initial_value="0.067156687129"/>
      <variable units="dimensionless" name="y_infinity"/>
      <variable units="millisecond" name="tau_y"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>y_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">64</cn>
                        </apply>
                        <cn cellml:units="millivolt">13.5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_y</ci>
            <apply>
               <divide/>
               <cn cellml:units="millisecond">0.7166529</cn>
               <apply>
                  <plus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">386.9</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">45.302</cn>
                     </apply>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">73.08</cn>
                        </apply>
                        <cn cellml:units="millivolt">19.231</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>y</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>y_infinity</ci>
                  <ci>y</ci>
               </apply>
               <ci>tau_y</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sustained_inward_current">
      <variable units="picoA" public_interface="out" name="i_st"/>
      <variable units="nanoS" name="g_st" initial_value="0.015"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millivolt" name="E_st" initial_value="37.4"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="qa"/>
      <variable units="dimensionless" private_interface="in" name="qi"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_st</ci>
            <apply>
               <times/>
               <ci>g_st</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_st</ci>
               </apply>
               <ci>qa</ci>
               <ci>qi</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sustained_inward_current_qa_gate">
      <variable units="dimensionless" public_interface="out" name="qa" initial_value="0.426018100136"/>
      <variable units="dimensionless" name="qa_infinity"/>
      <variable units="millisecond" name="tau_qa"/>
      <variable units="per_millisecond" name="alpha_qa"/>
      <variable units="per_millisecond" name="beta_qa"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>qa_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">57</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_qa</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>alpha_qa</ci>
                  <ci>beta_qa</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_qa</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">0.15</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">11</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">0.2</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">700</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_qa</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">16</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">8</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">15</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">50</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>qa</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>qa_infinity</ci>
                  <ci>qa</ci>
               </apply>
               <ci>tau_qa</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sustained_inward_current_qi_gate">
      <variable units="dimensionless" public_interface="out" name="qi" initial_value="0.333330378068"/>
      <variable units="dimensionless" name="qi_infinity"/>
      <variable units="millisecond" name="tau_qi"/>
      <variable units="per_millisecond" name="alpha_qi"/>
      <variable units="per_millisecond" name="beta_qi"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>qi_infinity</ci>
            <apply>
               <divide/>
               <ci>alpha_qi</ci>
               <apply>
                  <plus/>
                  <ci>alpha_qi</ci>
                  <ci>beta_qi</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_qi</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">6.65</cn>
               <apply>
                  <plus/>
                  <ci>alpha_qi</ci>
                  <ci>beta_qi</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_qi</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">3100</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">13</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millisecond">700</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">70</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_qi</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="millisecond">95</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <cn cellml:units="millivolt">10</cn>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="millisecond">50</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <cn cellml:units="millivolt">700</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <cn cellml:units="per_millisecond">0.000229</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                           </apply>
                           <cn cellml:units="millivolt">5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>qi</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>qi_infinity</ci>
                  <ci>qi</ci>
               </apply>
               <ci>tau_qi</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_dependent_background_current">
      <variable units="picoA" public_interface="out" name="i_b_Na"/>
      <variable units="nanoS" name="g_b_Na" initial_value="0.0054"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_Na</ci>
            <apply>
               <times/>
               <ci>g_b_Na</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="background_muscarinic_potassium_channel_current">
      <variable units="picoA" public_interface="out" name="i_K_ACh"/>
      <variable units="picoA" name="g_K_ACh"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millimolar" public_interface="in" name="Ki"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="attojoule_per_millimole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="femtocoulomb_per_millimole" public_interface="in" name="F"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>g_K_ACh</ci>
            <apply>
               <times/>
               <cn cellml:units="picoA">0.0011</cn>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <cn cellml:units="millimolar">1</cn>
                  </apply>
                  <cn cellml:units="dimensionless">0.41</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_K_ACh</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>g_K_ACh</ci>
                  <apply>
                     <minus/>
                     <ci>Ki</ci>
                     <apply>
                        <times/>
                        <ci>Ko</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                 </apply>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millimolar">1</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump_current">
      <variable units="picoA" public_interface="out" name="i_NaK"/>
      <variable units="millimolar" name="Km_Kp" initial_value="1.4"/>
      <variable units="millimolar" name="Km_Nap" initial_value="14"/>
      <variable units="picoA" name="i_NaK_max" initial_value="3.6"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaK</ci>
            <apply>
               <times/>
               <ci>i_NaK_max</ci>
               <apply>
                  <power/>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>Km_Kp</ci>
                           <ci>Ko</ci>
                        </apply>
                        <cn cellml:units="dimensionless">1.2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
               <apply>
                  <power/>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <power/>
                        <apply>
                           <divide/>
                           <ci>Km_Nap</ci>
                           <ci>Nai</ci>
                        </apply>
                        <cn cellml:units="dimensionless">1.3</cn>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
               <apply>
                  <power/>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <ci>E_Na</ci>
                                 </apply>
                                 <cn cellml:units="millivolt">120</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">30</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_calcium_exchange_current">
      <variable units="picoA" public_interface="out" name="i_NaCa"/>
      <variable units="picoA" name="kNaCa" initial_value="125"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="dimensionless" name="x1"/>
      <variable units="dimensionless" name="x2"/>
      <variable units="dimensionless" name="x3"/>
      <variable units="dimensionless" name="x4"/>
      <variable units="dimensionless" name="k41"/>
      <variable units="dimensionless" name="k34"/>
      <variable units="dimensionless" name="k23"/>
      <variable units="dimensionless" name="k21"/>
      <variable units="dimensionless" name="k32"/>
      <variable units="dimensionless" name="k43"/>
      <variable units="dimensionless" name="k12"/>
      <variable units="dimensionless" name="k14"/>
      <variable units="dimensionless" name="Qci" initial_value="0.1369"/>
      <variable units="dimensionless" name="Qn" initial_value="0.4315"/>
      <variable units="dimensionless" name="Qco" initial_value="0"/>
      <variable units="millimolar" name="K3ni" initial_value="26.44"/>
      <variable units="millimolar" name="Kci" initial_value="0.0207"/>
      <variable units="millimolar" name="K1ni" initial_value="395.3"/>
      <variable units="millimolar" name="K2ni" initial_value="2.289"/>
      <variable units="millimolar" name="Kcni" initial_value="26.44"/>
      <variable units="millimolar" name="K3no" initial_value="4.663"/>
      <variable units="millimolar" name="K1no" initial_value="1628"/>
      <variable units="millimolar" name="K2no" initial_value="561.4"/>
      <variable units="millimolar" name="Kco" initial_value="3.663"/>
      <variable units="dimensionless" name="do"/>
      <variable units="dimensionless" name="di"/>
      <variable units="millimolar" public_interface="in" name="Cao"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="millimolar" public_interface="in" name="Ca_sub"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="attojoule_per_millimole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="femtocoulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <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>
         <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>k43</ci>
            <apply>
               <divide/>
               <ci>Nai</ci>
               <apply>
                  <plus/>
                  <ci>K3ni</ci>
                  <ci>Nai</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k12</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Ca_sub</ci>
                     <ci>Kci</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Qci</ci>
                           </apply>
                           <ci>V</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>k14</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <ci>Nai</ci>
                           <ci>K1ni</ci>
                        </apply>
                        <ci>Nai</ci>
                     </apply>
                     <ci>K2ni</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Nai</ci>
                        <ci>K3ni</ci>
                     </apply>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>Qn</ci>
                           <ci>V</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>di</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k41</ci>
            <apply>
               <exp/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <ci>Qn</ci>
                     </apply>
                     <ci>V</ci>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>di</ci>
            <apply>
               <plus/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Ca_sub</ci>
                     <ci>Kci</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>Qci</ci>
                              </apply>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <ci>Nai</ci>
                        <ci>Kcni</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Nai</ci>
                     <ci>K1ni</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <ci>Nai</ci>
                           <ci>K2ni</ci>
                        </apply>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <ci>Nai</ci>
                              <ci>K3ni</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </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>k21</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Kco</ci>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>Qco</ci>
                           <ci>V</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>k23</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <ci>Nao</ci>
                           <ci>K1no</ci>
                        </apply>
                        <ci>Nao</ci>
                     </apply>
                     <ci>K2no</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Nao</ci>
                        <ci>K3no</ci>
                     </apply>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Qn</ci>
                           </apply>
                           <ci>V</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <ci>do</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k32</ci>
            <apply>
               <exp/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>Qn</ci>
                     <ci>V</ci>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>do</ci>
            <apply>
               <plus/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Kco</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>Qco</ci>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>Nao</ci>
                     <ci>K1no</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Nao</ci>
                        <ci>K2no</ci>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Nao</ci>
                        <ci>K3no</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_calcium_dynamics">
      <variable units="millimolar_per_millisecond" public_interface="out" name="j_Ca_dif"/>
      <variable units="millimolar_per_millisecond" public_interface="out" name="j_rel"/>
      <variable units="millimolar_per_millisecond" public_interface="out" name="j_up"/>
      <variable units="millimolar_per_millisecond" public_interface="out" name="j_tr"/>
      <variable units="millisecond" name="tau_dif_Ca" initial_value="0.04"/>
      <variable units="millisecond" name="tau_tr" initial_value="60"/>
      <variable units="millimolar" name="K_rel" initial_value="0.0012"/>
      <variable units="millimolar_per_millisecond" name="P_up" initial_value="0.005"/>
      <variable units="per_millisecond" name="P_rel" initial_value="0.5"/>
      <variable units="millimolar" name="K_up" initial_value="0.0006"/>
      <variable units="millimolar" public_interface="in" name="Ca_up"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <variable units="millimolar" public_interface="in" name="Ca_sub"/>
      <variable units="millimolar" public_interface="in" name="Ca_rel"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>j_Ca_dif</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Ca_sub</ci>
                  <ci>Cai</ci>
               </apply>
               <ci>tau_dif_Ca</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>j_rel</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>P_rel</ci>
                  <apply>
                     <minus/>
                     <ci>Ca_rel</ci>
                     <ci>Ca_sub</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>K_rel</ci>
                        <ci>Ca_sub</ci>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>j_up</ci>
            <apply>
               <divide/>
               <ci>P_up</ci>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>K_up</ci>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>j_tr</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Ca_up</ci>
                  <ci>Ca_rel</ci>
               </apply>
               <ci>tau_tr</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_ion_concentrations">
      <variable units="millimolar" public_interface="out" name="Cai" initial_value="0.000312494921"/>
      <variable units="millimolar" public_interface="out" name="Cao" initial_value="2"/>
      <variable units="millimolar" public_interface="out" name="Ca_up" initial_value="1.462338380106"/>
      <variable units="millimolar" public_interface="out" name="Ca_rel" initial_value="0.296742023718"/>
      <variable units="millimolar" public_interface="out" name="Ca_sub" initial_value="0.00019074741"/>
      <variable units="millimolar" public_interface="out" name="Nai" initial_value="9.438646305915"/>
      <variable units="millimolar" public_interface="out" name="Nao" initial_value="140"/>
      <variable units="millimolar" public_interface="out" name="Ki" initial_value="139.984146485614"/>
      <variable units="millimolar" public_interface="out" name="Ko" initial_value="5.4"/>
      <variable units="millimolar" public_interface="out" name="Mgi" initial_value="2.5"/>
      <variable units="litre" name="V_i" initial_value="0.0000000000015835"/>
      <variable units="litre" name="V_rel" initial_value="0.0000000000000042223"/>
      <variable units="litre" name="V_up" initial_value="0.000000000000040816"/>
      <variable units="litre" name="V_sub" initial_value="0.000000000000035098"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="femtocoulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="picoF" public_interface="in" name="Cm"/>
      <variable units="picoA" public_interface="in" name="i_b_Na"/>
      <variable units="picoA" public_interface="in" name="i_CaL"/>
      <variable units="picoA" public_interface="in" name="i_CaT"/>
      <variable units="picoA" public_interface="in" name="i_NaCa"/>
      <variable units="millimolar_per_millisecond" public_interface="in" name="j_rel"/>
      <variable units="millimolar_per_millisecond" public_interface="in" name="j_Ca_dif"/>
      <variable units="millimolar_per_millisecond" public_interface="in" name="j_tr"/>
      <variable units="millimolar_per_millisecond" public_interface="in" name="j_up"/>
      <variable units="picoA" public_interface="in" name="i_h_Na"/>
      <variable units="picoA" public_interface="in" name="i_st"/>
      <variable units="picoA" public_interface="in" name="i_NaK"/>
      <variable units="picoA" public_interface="in" name="i_Kr"/>
      <variable units="picoA" public_interface="in" name="i_Ks"/>
      <variable units="picoA" public_interface="in" name="i_to"/>
      <variable units="picoA" public_interface="in" name="i_sus"/>
      <variable units="picoA" public_interface="in" name="i_h_K"/>
      <variable units="picoA" public_interface="in" name="i_K_ACh"/>
      <variable units="dimensionless" public_interface="in" name="TMC_tot"/>
      <variable units="dimensionless" public_interface="in" name="CM_tot"/>
      <variable units="dimensionless" public_interface="in" name="TC_tot"/>
      <variable units="dimensionless" public_interface="in" name="CQ_tot"/>
      <variable units="millimolar_per_millisecond" public_interface="in" name="delta_fTMC"/>
      <variable units="millimolar_per_millisecond" public_interface="in" name="delta_fCMi"/>
      <variable units="millimolar_per_millisecond" public_interface="in" name="delta_fCMs"/>
      <variable units="millimolar_per_millisecond" public_interface="in" name="delta_fTC"/>
      <variable units="millimolar_per_millisecond" public_interface="in" name="delta_fCQ"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Nai</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>i_h_Na</ci>
                        <ci>i_st</ci>
                        <ci>i_b_Na</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">3</cn>
                           <ci>i_NaK</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">3</cn>
                           <ci>i_NaCa</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="picoF">1</cn>
                  <ci>F</ci>
                  <apply>
                     <plus/>
                     <ci>V_i</ci>
                     <ci>V_sub</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ki</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>i_Kr</ci>
                        <ci>i_Ks</ci>
                        <ci>i_to</ci>
                        <ci>i_sus</ci>
                        <ci>i_h_K</ci>
                        <ci>i_K_ACh</ci>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <ci>i_NaK</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="picoF">1</cn>
                  <ci>F</ci>
                  <apply>
                     <plus/>
                     <ci>V_i</ci>
                     <ci>V_sub</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cai</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>j_Ca_dif</ci>
                        <ci>V_sub</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>j_up</ci>
                        <ci>V_up</ci>
                     </apply>
                  </apply>
                  <ci>V_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>CM_tot</ci>
                     <ci>delta_fCMi</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>TC_tot</ci>
                     <ci>delta_fTC</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>TMC_tot</ci>
                     <ci>delta_fTMC</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_sub</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>i_CaL</ci>
                                    <ci>i_CaT</ci>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">2</cn>
                                    <ci>i_NaCa</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <ci>Cm</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="picoF">1</cn>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>F</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>j_rel</ci>
                        <ci>V_rel</ci>
                     </apply>
                  </apply>
                  <ci>V_sub</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>j_Ca_dif</ci>
                  <apply>
                     <times/>
                     <ci>CM_tot</ci>
                     <ci>delta_fCMs</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_up</ci>
            </apply>
            <apply>
               <minus/>
               <ci>j_up</ci>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>j_tr</ci>
                     <ci>V_rel</ci>
                  </apply>
                  <ci>V_up</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_rel</ci>
            </apply>
            <apply>
               <minus/>
               <ci>j_tr</ci>
               <apply>
                  <plus/>
                  <ci>j_rel</ci>
                  <apply>
                     <times/>
                     <ci>CQ_tot</ci>
                     <ci>delta_fCQ</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_buffering">
      <variable units="dimensionless" public_interface="out" name="TC_tot" initial_value="0.031"/>
      <variable units="dimensionless" public_interface="out" name="TMC_tot" initial_value="0.062"/>
      <variable units="dimensionless" public_interface="out" name="CM_tot" initial_value="0.045"/>
      <variable units="dimensionless" public_interface="out" name="CQ_tot" initial_value="10"/>
      <variable units="millimolar_per_millisecond" public_interface="out" name="delta_fTC"/>
      <variable units="millimolar_per_millisecond" public_interface="out" name="delta_fTMC"/>
      <variable units="millimolar_per_millisecond" public_interface="out" name="delta_fCMs"/>
      <variable units="millimolar_per_millisecond" public_interface="out" name="delta_fCMi"/>
      <variable units="millimolar_per_millisecond" public_interface="out" name="delta_fCQ"/>
      <variable units="millimolar_per_millisecond" name="delta_fTMM"/>
      <variable units="millimolar" name="fTMM" initial_value="0.350600895635"/>
      <variable units="millimolar" name="fCMi" initial_value="0.116947220413"/>
      <variable units="millimolar" name="fCMs" initial_value="0.074631965653"/>
      <variable units="millimolar" name="fTC" initial_value="0.059206293446"/>
      <variable units="millimolar" name="fTMC" initial_value="0.602955114871"/>
      <variable units="millimolar" name="fCQ" initial_value="0.260317260703"/>
      <variable units="per_millimolar_millisecond" name="kf_TC" initial_value="88.8"/>
      <variable units="per_millimolar_millisecond" name="kf_TMM" initial_value="2.277"/>
      <variable units="per_millimolar_millisecond" name="kf_TMC" initial_value="227.7"/>
      <variable units="per_millimolar_millisecond" name="kf_CM" initial_value="227.7"/>
      <variable units="per_millimolar_millisecond" name="kf_CQ" initial_value="0.534"/>
      <variable units="per_millisecond" name="kb_TC" initial_value="0.446"/>
      <variable units="per_millisecond" name="kb_TMC" initial_value="0.00751"/>
      <variable units="per_millisecond" name="kb_TMM" initial_value="0.751"/>
      <variable units="per_millisecond" name="kb_CM" initial_value="0.542"/>
      <variable units="per_millisecond" name="kb_CQ" initial_value="0.445"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="millimolar" public_interface="in" name="Mgi"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <variable units="millimolar" public_interface="in" name="Ca_sub"/>
      <variable units="millimolar" public_interface="in" name="Ca_rel"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>fTC</ci>
            </apply>
            <ci>delta_fTC</ci>
         </apply>
         <apply>
            <eq/>
            <ci>delta_fTC</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kf_TC</ci>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="millimolar">1</cn>
                     <ci>fTC</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>kb_TC</ci>
                  <ci>fTC</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>fTMC</ci>
            </apply>
            <ci>delta_fTMC</ci>
         </apply>
         <apply>
            <eq/>
            <ci>delta_fTMC</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kf_TMC</ci>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="millimolar">1</cn>
                     <apply>
                        <plus/>
                        <ci>fTMC</ci>
                        <ci>fTMM</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>kb_TMC</ci>
                  <ci>fTMC</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>fTMM</ci>
            </apply>
            <ci>delta_fTMM</ci>
         </apply>
         <apply>
            <eq/>
            <ci>delta_fTMM</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kf_TMM</ci>
                  <ci>Mgi</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="millimolar">1</cn>
                     <apply>
                        <plus/>
                        <ci>fTMC</ci>
                        <ci>fTMM</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>kb_TMM</ci>
                  <ci>fTMM</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>fCMi</ci>
            </apply>
            <ci>delta_fCMi</ci>
         </apply>
         <apply>
            <eq/>
            <ci>delta_fCMi</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kf_CM</ci>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="millimolar">1</cn>
                     <ci>fCMi</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>kb_CM</ci>
                  <ci>fCMi</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>fCMs</ci>
            </apply>
            <ci>delta_fCMs</ci>
         </apply>
         <apply>
            <eq/>
            <ci>delta_fCMs</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kf_CM</ci>
                  <ci>Ca_sub</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="millimolar">1</cn>
                     <ci>fCMs</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>kb_CM</ci>
                  <ci>fCMs</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>fCQ</ci>
            </apply>
            <ci>delta_fCQ</ci>
         </apply>
         <apply>
            <eq/>
            <ci>delta_fCQ</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kf_CQ</ci>
                  <ci>Ca_rel</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="millimolar">1</cn>
                     <ci>fCQ</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>kb_CQ</ci>
                  <ci>fCQ</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="L_type_calcium_channel_current">
            <component_ref component="L_type_calcium_channel_current_d_gate"/>
            <component_ref component="L_type_calcium_channel_current_f_gate"/>
            <component_ref component="L_type_calcium_channel_current_fCa_gate"/>
         </component_ref>
         <component_ref component="T_type_calcium_channel_current">
            <component_ref component="T_type_calcium_channel_current_d_gate"/>
            <component_ref component="T_type_calcium_channel_current_f_gate"/>
         </component_ref>
         <component_ref component="rapidly_activating_delayed_rectifier_potassium_current">
            <component_ref component="rapidly_activating_delayed_rectifier_potassium_current_pa_gate"/>
            <component_ref component="rapidly_activating_delayed_rectifier_potassium_current_pi_gate"/>
         </component_ref>
         <component_ref component="slowly_activating_delayed_rectifier_potassium_current">
            <component_ref component="slowly_activating_delayed_rectifier_potassium_current_n_gate"/>
         </component_ref>
         <component_ref component="AP_sensitive_currents">
            <component_ref component="AP_sensitive_currents_q_gate"/>
            <component_ref component="AP_sensitive_currents_r_gate"/>
         </component_ref>
         <component_ref component="hyperpolarisation_activated_current">
            <component_ref component="hyperpolarisation_activated_current_y_gate"/>
         </component_ref>
         <component_ref component="sustained_inward_current">
            <component_ref component="sustained_inward_current_qa_gate"/>
            <component_ref component="sustained_inward_current_qi_gate"/>
         </component_ref>
         <component_ref component="sodium_dependent_background_current"/>
         <component_ref component="background_muscarinic_potassium_channel_current"/>
         <component_ref component="sodium_potassium_pump_current"/>
         <component_ref component="sodium_calcium_exchange_current"/>
         <component_ref component="intracellular_calcium_dynamics"/>
         <component_ref component="intracellular_ion_concentrations"/>
         <component_ref component="calcium_buffering"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="L_type_calcium_channel_current">
         <component_ref component="L_type_calcium_channel_current_d_gate"/>
         <component_ref component="L_type_calcium_channel_current_f_gate"/>
         <component_ref component="L_type_calcium_channel_current_fCa_gate"/>
      </component_ref>
      <component_ref component="T_type_calcium_channel_current">
         <component_ref component="T_type_calcium_channel_current_d_gate"/>
         <component_ref component="T_type_calcium_channel_current_f_gate"/>
      </component_ref>
      <component_ref component="rapidly_activating_delayed_rectifier_potassium_current">
         <component_ref component="rapidly_activating_delayed_rectifier_potassium_current_pa_gate"/>
         <component_ref component="rapidly_activating_delayed_rectifier_potassium_current_pi_gate"/>
      </component_ref>
      <component_ref component="slowly_activating_delayed_rectifier_potassium_current">
         <component_ref component="slowly_activating_delayed_rectifier_potassium_current_n_gate"/>
      </component_ref>
      <component_ref component="AP_sensitive_currents">
         <component_ref component="AP_sensitive_currents_q_gate"/>
         <component_ref component="AP_sensitive_currents_r_gate"/>
      </component_ref>
      <component_ref component="hyperpolarisation_activated_current">
         <component_ref component="hyperpolarisation_activated_current_y_gate"/>
      </component_ref>
      <component_ref component="sustained_inward_current">
         <component_ref component="sustained_inward_current_qa_gate"/>
         <component_ref component="sustained_inward_current_qi_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_2="environment" component_1="membrane"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="L_type_calcium_channel_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="T_type_calcium_channel_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="slowly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="AP_sensitive_currents"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="hyperpolarisation_activated_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sustained_inward_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_dependent_background_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="background_muscarinic_potassium_channel_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_calcium_exchange_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_potassium_pump_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="intracellular_ion_concentrations"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_buffering"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="reversal_potentials"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="T" variable_1="T"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="L_type_calcium_channel_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_CaL" variable_1="i_CaL"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="T_type_calcium_channel_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_CaT" variable_1="i_CaT"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="slowly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="AP_sensitive_currents"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_to" variable_1="i_to"/>
      <map_variables variable_2="i_sus" variable_1="i_sus"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="hyperpolarisation_activated_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_h" variable_1="i_h"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="sustained_inward_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_st" variable_1="i_st"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="sodium_dependent_background_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="background_muscarinic_potassium_channel_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K_ACh" variable_1="i_K_ACh"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="sodium_potassium_pump_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="sodium_calcium_exchange_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="intracellular_ion_concentrations"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="Cm" variable_1="Cm"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="reversal_potentials"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
   </connection>
   <connection>
      <map_components component_2="rapidly_activating_delayed_rectifier_potassium_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="AP_sensitive_currents" component_1="reversal_potentials"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="reversal_potentials" component_1="hyperpolarisation_activated_current"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="sodium_dependent_background_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="sodium_potassium_pump_current" component_1="reversal_potentials"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="L_type_calcium_channel_current"/>
      <map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
      <map_variables variable_2="i_CaL" variable_1="i_CaL"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="T_type_calcium_channel_current"/>
      <map_variables variable_2="i_CaT" variable_1="i_CaT"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="slowly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="AP_sensitive_currents"/>
      <map_variables variable_2="i_to" variable_1="i_to"/>
      <map_variables variable_2="i_sus" variable_1="i_sus"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="hyperpolarisation_activated_current"/>
      <map_variables variable_2="i_h_Na" variable_1="i_h_Na"/>
      <map_variables variable_2="i_h_K" variable_1="i_h_K"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sustained_inward_current"/>
      <map_variables variable_2="i_st" variable_1="i_st"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sodium_dependent_background_current"/>
      <map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="background_muscarinic_potassium_channel_current"/>
      <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="i_K_ACh" variable_1="i_K_ACh"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sodium_potassium_pump_current"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sodium_calcium_exchange_current"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
      <map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="intracellular_calcium_dynamics"/>
      <map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
      <map_variables variable_2="Ca_up" variable_1="Ca_up"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
      <map_variables variable_2="j_Ca_dif" variable_1="j_Ca_dif"/>
      <map_variables variable_2="j_up" variable_1="j_up"/>
      <map_variables variable_2="j_rel" variable_1="j_rel"/>
      <map_variables variable_2="j_tr" variable_1="j_tr"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="calcium_buffering"/>
      <map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
      <map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="Mgi" variable_1="Mgi"/>
      <map_variables variable_2="TC_tot" variable_1="TC_tot"/>
      <map_variables variable_2="TMC_tot" variable_1="TMC_tot"/>
      <map_variables variable_2="CM_tot" variable_1="CM_tot"/>
      <map_variables variable_2="CQ_tot" variable_1="CQ_tot"/>
      <map_variables variable_2="delta_fCMi" variable_1="delta_fCMi"/>
      <map_variables variable_2="delta_fCMs" variable_1="delta_fCMs"/>
      <map_variables variable_2="delta_fTMC" variable_1="delta_fTMC"/>
      <map_variables variable_2="delta_fTC" variable_1="delta_fTC"/>
      <map_variables variable_2="delta_fCQ" variable_1="delta_fCQ"/>
   </connection>
   <connection>
      <map_components component_2="L_type_calcium_channel_current_d_gate" component_1="L_type_calcium_channel_current"/>
      <map_variables variable_2="d" variable_1="d"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_calcium_channel_current_f_gate" component_1="L_type_calcium_channel_current"/>
      <map_variables variable_2="f" variable_1="f"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_calcium_channel_current_fCa_gate" component_1="L_type_calcium_channel_current"/>
      <map_variables variable_2="fCa" variable_1="fCa"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
   </connection>
   <connection>
      <map_components component_2="T_type_calcium_channel_current_d_gate" component_1="T_type_calcium_channel_current"/>
      <map_variables variable_2="d" variable_1="d"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="T_type_calcium_channel_current_f_gate" component_1="T_type_calcium_channel_current"/>
      <map_variables variable_2="f" variable_1="f"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="rapidly_activating_delayed_rectifier_potassium_current_pa_gate" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="paF" variable_1="paF"/>
      <map_variables variable_2="paS" variable_1="paS"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="rapidly_activating_delayed_rectifier_potassium_current_pi_gate" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="piy" variable_1="piy"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="slowly_activating_delayed_rectifier_potassium_current_n_gate" component_1="slowly_activating_delayed_rectifier_potassium_current"/>
      <map_variables variable_2="n" variable_1="n"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="AP_sensitive_currents_q_gate" component_1="AP_sensitive_currents"/>
      <map_variables variable_2="q" variable_1="q"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="AP_sensitive_currents_r_gate" component_1="AP_sensitive_currents"/>
      <map_variables variable_2="r" variable_1="r"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="hyperpolarisation_activated_current_y_gate" component_1="hyperpolarisation_activated_current"/>
      <map_variables variable_2="y" variable_1="y"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sustained_inward_current_qa_gate" component_1="sustained_inward_current"/>
      <map_variables variable_2="qa" variable_1="qa"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sustained_inward_current_qi_gate" component_1="sustained_inward_current"/>
      <map_variables variable_2="qi" variable_1="qi"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#4e7eb78e-506d-42db-ae89-95af5cf22e4a">
    <rdf:li rdf:resource="rdf:#a4f239de-cae4-4771-93e4-36cae1f80104"/>
    <rdf:li rdf:resource="rdf:#bb38323f-44d0-497f-bad3-9a0f638f0cd1"/>
    <rdf:li rdf:resource="rdf:#0632c177-07ba-4185-a6ed-98d9c540e9c0"/>
    <rdf:li rdf:resource="rdf:#a16b9599-a60b-4d6d-91a3-6c7f03f887d2"/>
  </rdf:Seq>
  <rdf:Description rdf:about="">
    <dc:publisher/>
    <cmeta:comment rdf:resource="rdf:#ce8a4b31-1d53-4c25-af24-3d1c47e0737a"/>
    <dcterms:created rdf:resource="rdf:#52f55f93-32da-4cfc-9417-6812a46b8719"/>
    <dc:creator rdf:resource="rdf:#ef4eebea-cfdb-452e-b8b1-409705d845bf"/>
    <cmeta:modification rdf:resource="rdf:#b27c7efb-dfa1-4ed4-95af-e10caf46a28c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#52f55f93-32da-4cfc-9417-6812a46b8719">
    <dcterms:W3CDTF>2008-05-08T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b27c7efb-dfa1-4ed4-95af-e10caf46a28c">
    <dcterms:modified rdf:resource="rdf:#58edd3a2-b5a6-480a-90e3-f63f8f4922da"/>
    <rdf:value>Curation of units.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#8fd6885b-b079-4a8d-9ed3-ff0599ac7956"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a2e8aec2-8234-4020-99a9-e71204fd9b79">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ce8a4b31-1d53-4c25-af24-3d1c47e0737a">
    <dc:creator rdf:resource="rdf:#599a52b6-bb36-4be6-9964-81b44645ac00"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a16b9599-a60b-4d6d-91a3-6c7f03f887d2">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#32e4216e-118c-4ee5-b74a-80267e5c6bdb"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8fd6885b-b079-4a8d-9ed3-ff0599ac7956">
    <vCard:N rdf:resource="rdf:#6ce2ba4b-85dc-46af-8b8b-9e14b236aae3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#659ea48f-5eb5-4872-a128-4e5a7445e1ee">
    <dc:creator rdf:resource="rdf:#4e7eb78e-506d-42db-ae89-95af5cf22e4a"/>
    <dc:title>Dynamical description of sinoatrial node pacemaking: improved mathematical model for primary pacemaker cell</dc:title>
    <bqs:volume>283</bqs:volume>
    <bqs:first_page>2074</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#2030f9b7-8e59-4df2-9df2-15f151f85b3c"/>
    <dcterms:issued rdf:resource="rdf:#d5de7ad3-8979-49b2-93a5-b4ad967a1a51"/>
    <bqs:last_page>2101</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="#kurata_hisatome_imanishi_shibamoto_2006_version02">
    <bqs:reference rdf:resource="rdf:#996c404c-7126-4a04-9634-cd5a09e99615"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject 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>
    <cmeta:comment rdf:resource="rdf:#e2a13f7b-ac1d-4019-a68b-3850fcb56eae"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d5de7ad3-8979-49b2-93a5-b4ad967a1a51">
    <dcterms:W3CDTF>2002-11-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bb38323f-44d0-497f-bad3-9a0f638f0cd1">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#d20a1024-60ed-4afd-8526-aad9e6686ce0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e2a13f7b-ac1d-4019-a68b-3850fcb56eae">
    <dc:creator rdf:resource="rdf:#a2e8aec2-8234-4020-99a9-e71204fd9b79"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d20a1024-60ed-4afd-8526-aad9e6686ce0">
    <vCard:Given>I</vCard:Given>
    <vCard:Family>Hisatome</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0632c177-07ba-4185-a6ed-98d9c540e9c0">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#79ca0021-882f-43ec-a764-453779d76d08"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ef4eebea-cfdb-452e-b8b1-409705d845bf">
    <vCard:ORG rdf:resource="rdf:#2b7056bd-c110-421b-953b-b83e49e71f02"/>
    <vCard:EMAIL rdf:resource="rdf:#e060641a-97bd-400e-b5f7-58dc50b3a430"/>
    <vCard:N rdf:resource="rdf:#5c63855c-81a3-48ad-aed6-93977249c9e3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2030f9b7-8e59-4df2-9df2-15f151f85b3c">
    <dc:title>Am J Physiol Heart Circ Physiol</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#58edd3a2-b5a6-480a-90e3-f63f8f4922da">
    <dcterms:W3CDTF>2008-05-08T04:35:11+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#32e4216e-118c-4ee5-b74a-80267e5c6bdb">
    <vCard:Given>T</vCard:Given>
    <vCard:Family>Shibamoto</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2b7056bd-c110-421b-953b-b83e49e71f02">
    <vCard:Orgname>University of Oxford</vCard:Orgname>
    <vCard:Orgunit/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5e038ad6-49e7-42db-ad02-5109c177fe3e">
    <vCard:Given>Y</vCard:Given>
    <vCard:Family>Kurata</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#996c404c-7126-4a04-9634-cd5a09e99615">
    <bqs:Pubmed_id>12384487</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#659ea48f-5eb5-4872-a128-4e5a7445e1ee"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#599a52b6-bb36-4be6-9964-81b44645ac00">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c63855c-81a3-48ad-aed6-93977249c9e3">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#79ca0021-882f-43ec-a764-453779d76d08">
    <vCard:Given>S</vCard:Given>
    <vCard:Family>Imanishi</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a4f239de-cae4-4771-93e4-36cae1f80104">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#5e038ad6-49e7-42db-ad02-5109c177fe3e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6ce2ba4b-85dc-46af-8b8b-9e14b236aae3">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e060641a-97bd-400e-b5f7-58dc50b3a430">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>penny.noble@dpag.ox.ac.uk</rdf:value>
  </rdf:Description>
</rdf:RDF>
</model>