Location: Nygren, Fiset, Firek, Clark, Lindblad, Clark, Giles, 1998 @ 51ef5ce71aa6 / nygren_fiset_firek_clark_lindblad_clark_giles_1998.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2010-02-11 14:05:06+13:00
Desc:
Updated model documentation, fixed e-notation and removed a redundant <ulink> tag from the citation.
Permanent Source URI:
https://staging.physiomeproject.org/workspace/nygren_fiset_firek_clark_lindblad_clark_giles_1998/rawfile/51ef5ce71aa6052f0685b6eb294c8828004b7a85/nygren_fiset_firek_clark_lindblad_clark_giles_1998.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 16/01/2008 at 16:52:25 using:

COR (0.9.31.857)
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:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="nygren_fiset_firek_clark_lindblad_clark_giles_1998" name="nygren_fiset_firek_clark_lindblad_clark_giles_1998">
   
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Mathematical model of an adult human atrial cell: the role of K+ currents in repolarization</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This CellML model runs in both PCEnv and COR to recreate the published results. The units have been checked and they are consistent.
          </para>

  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: We have developed a mathematical model of the human atria myocyte based on averaged voltage-clamp data recorded from isolated single myocytes. Our model consists of a Hodgkin-Huxley-type equivalent circuit for the sarcolemma, coupled with a fluid compartment model, which accounts for changes in ionic concentrations in the cytoplasm as well as in the sarcoplasmic reticulum. This formulation can reconstruct action potential data that are representative of recordings from a majority of human atrial cells in our laboratory and therefore provides a biophysically based account of the underlying ionic currents. This work is based in part on a previous model of the rabbit atrial myocyte published by our group and was motivated by differences in some of the repolarizing currents between human and rabbit atrium. We have therefore given particular attention to the sustained outward K+ current (I[sus]), which putatively has a prominent role in determining the duration of the human atrial action potential. Our results demonstrate that the action potential shape during the peak and plateau phases is determined primarily by transient outward K+ current, I(sus) and L-type Ca2+ current (I[Ca,L]) and that the role of I(sus) in the human atrial action potential can be modulated by the baseline sizes of I(Ca,L), I(sus) and the rapid delayed rectifier K+ current. As a result, our simulations suggest that the functional role of I(sus) can depend on the physiological/disease state of the cell.</para>

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

<para>
Mathematical model of an adult human atrial cell: the role of K+ currents in repolarization, Nygren A, Fiset C, Firek L, Clark JW, Lindblad DS, Clark RB and Giles WR, 1998,<emphasis>Circulation Research</emphasis>, 82, 63-81.  <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/9440706">PubMed ID: 9440706</ulink> 
</para>

<informalfigure float="0" id="fig_fluid_compartment_model">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>fluid compartment model of a human atrial cell</title>
    </objectinfo>
    <imagedata fileref="nygren_1998.png"/>
  </imageobject>
</mediaobject>
<caption>A fluid compartment model of the human atrial cell, including intracellular, cleft and extracellular spaces.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="per_millimolar_4">
      <unit units="millimolar" exponent="-4"/>
   </units>
   <units name="picoA_per_millimolar">
      <unit units="ampere" prefix="pico"/>
      <unit units="millimolar" exponent="-1"/>
   </units>
   <units name="picoA_per_millimolar_4">
      <unit units="ampere" prefix="pico"/>
      <unit units="millimolar" exponent="-4"/>
   </units>
   <units name="per_millimolar_second">
      <unit units="millimolar" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="picoA">
      <unit units="ampere" prefix="pico"/>
   </units>
   <units name="nanoS">
      <unit units="siemens" prefix="nano"/>
   </units>
   <units name="nanoF">
      <unit units="farad" prefix="nano"/>
   </units>
   <units name="nanolitre">
      <unit units="litre" prefix="nano"/>
   </units>
   <units name="nanolitre_per_second">
      <unit units="litre" prefix="nano"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millijoule_per_mole_kelvin">
      <unit units="joule" prefix="milli"/>
      <unit units="mole" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <component name="environment">
      <variable units="second" public_interface="out" name="time"/>
   </component>
   <component name="membrane">
      <variable units="millivolt" public_interface="out" name="V" initial_value="-74.2525"/>
      <variable units="millijoule_per_mole_kelvin" public_interface="out" name="R" initial_value="8314"/>
      <variable units="kelvin" public_interface="out" name="T" initial_value="306.15"/>
      <variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96487"/>
      <variable units="nanoF" name="Cm" initial_value="0.05"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="picoA" public_interface="in" name="i_Na"/>
      <variable units="picoA" public_interface="in" name="i_Ca_L"/>
      <variable units="picoA" public_interface="in" name="i_t"/>
      <variable units="picoA" public_interface="in" name="i_sus"/>
      <variable units="picoA" public_interface="in" name="i_K1"/>
      <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_B_Na"/>
      <variable units="picoA" public_interface="in" name="i_B_Ca"/>
      <variable units="picoA" public_interface="in" name="i_NaK"/>
      <variable units="picoA" public_interface="in" name="i_CaP"/>
      <variable units="picoA" public_interface="in" name="i_NaCa"/>
      <variable units="picoA" name="i_Stim"/>
      <variable units="second" name="stim_start" initial_value="0.1"/>
      <variable units="second" name="stim_end" initial_value="100000000"/>
      <variable units="second" name="stim_period" initial_value="1"/>
      <variable units="second" name="stim_duration" initial_value="0.006"/>
      <variable units="picoA" name="stim_amplitude" initial_value="-280"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Stim</ci>
            <piecewise>
               <piece>
                  <ci>stim_amplitude</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>time</ci>
                        <ci>stim_start</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>time</ci>
                        <ci>stim_end</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>time</ci>
                              <ci>stim_start</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>stim_start</ci>
                                    </apply>
                                    <ci>stim_period</ci>
                                 </apply>
                              </apply>
                              <ci>stim_period</ci>
                           </apply>
                        </apply>
                        <ci>stim_duration</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="picoA">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>i_Stim</ci>
                  <ci>i_Na</ci>
                  <ci>i_Ca_L</ci>
                  <ci>i_t</ci>
                  <ci>i_sus</ci>
                  <ci>i_K1</ci>
                  <ci>i_Kr</ci>
                  <ci>i_Ks</ci>
                  <ci>i_B_Na</ci>
                  <ci>i_B_Ca</ci>
                  <ci>i_NaK</ci>
                  <ci>i_CaP</ci>
                  <ci>i_NaCa</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current">
      <variable units="picoA" public_interface="out" name="i_Na"/>
      <variable units="millivolt" public_interface="out" name="E_Na"/>
      <variable units="nanolitre_per_second" name="P_Na" initial_value="0.0016"/>
      <variable units="second" 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="Na_c"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="dimensionless" private_interface="in" name="m"/>
      <variable units="dimensionless" private_interface="in" name="h1"/>
      <variable units="dimensionless" private_interface="in" name="h2"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>P_Na</ci>
                        <apply>
                           <power/>
                           <ci>m</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">0.9</cn>
                              <ci>h1</ci>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">0.1</cn>
                              <ci>h2</ci>
                           </apply>
                        </apply>
                        <ci>Na_c</ci>
                        <ci>V</ci>
                        <apply>
                           <power/>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <ci>E_Na</ci>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>V</ci>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Na_c</ci>
                     <ci>Na_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current_m_gate">
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.0032017"/>
      <variable units="dimensionless" name="m_infinity"/>
      <variable units="second" name="tau_m"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>m_infinity</ci>
                  <ci>m</ci>
               </apply>
               <ci>tau_m</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>m_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">27.12</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">8.21</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_m</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn type="e-notation" cellml:units="second">4.2<sep/>-5</cn>
				  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">25.57</cn>
                              </apply>
                              <cn cellml:units="millivolt">28.8</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn type="e-notation" cellml:units="second">2.4<sep/>-5</cn>
			</apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current_h1_gate">
      <variable units="dimensionless" public_interface="out" name="h1" initial_value="0.8814"/>
      <variable units="dimensionless" public_interface="out" name="h_infinity"/>
      <variable units="second" name="tau_h1"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h1</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h_infinity</ci>
                  <ci>h1</ci>
               </apply>
               <ci>tau_h1</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>h_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">63.6</cn>
                        </apply>
                        <cn cellml:units="millivolt">5.3</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_h1</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.03</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">35.1</cn>
                           </apply>
                           <cn cellml:units="millivolt">3.2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.0003</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current_h2_gate">
      <variable units="dimensionless" public_interface="out" name="h2" initial_value="0.8742"/>
      <variable units="second" name="tau_h2"/>
      <variable units="dimensionless" public_interface="in" name="h_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h2</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h_infinity</ci>
                  <ci>h2</ci>
               </apply>
               <ci>tau_h2</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_h2</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.12</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">35.1</cn>
                           </apply>
                           <cn cellml:units="millivolt">3.2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.003</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel">
      <variable units="picoA" public_interface="out" name="i_Ca_L"/>
      <variable units="nanoS" name="g_Ca_L" initial_value="6.75"/>
      <variable units="millivolt" name="E_Ca_app" initial_value="60"/>
      <variable units="dimensionless" name="f_Ca"/>
      <variable units="millimolar" name="k_Ca" initial_value="0.025"/>
      <variable units="second" 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="Ca_d"/>
      <variable units="dimensionless" private_interface="in" name="d_L"/>
      <variable units="dimensionless" private_interface="in" name="f_L_1"/>
      <variable units="dimensionless" private_interface="in" name="f_L_2"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca_L</ci>
            <apply>
               <times/>
               <ci>g_Ca_L</ci>
               <ci>d_L</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>f_Ca</ci>
                     <ci>f_L_1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>f_Ca</ci>
                     </apply>
                     <ci>f_L_2</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca_app</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_Ca</ci>
            <apply>
               <divide/>
               <ci>Ca_d</ci>
               <apply>
                  <plus/>
                  <ci>Ca_d</ci>
                  <ci>k_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_d_L_gate">
      <variable units="dimensionless" public_interface="out" name="d_L" initial_value="1.3005e-5"/>
      <variable units="dimensionless" name="d_L_infinity"/>
      <variable units="second" name="tau_d_L"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d_L</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>d_L_infinity</ci>
                  <ci>d_L</ci>
               </apply>
               <ci>tau_d_L</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>d_L_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">9</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">5.8</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_d_L</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.0027</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">35</cn>
                              </apply>
                              <cn cellml:units="millivolt">30</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.002</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f_L1_gate">
      <variable units="dimensionless" public_interface="out" name="f_L_1" initial_value="0.9986"/>
      <variable units="dimensionless" public_interface="out" name="f_L_infinity"/>
      <variable units="second" name="tau_f_L1"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_L_1</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_L_infinity</ci>
                  <ci>f_L_1</ci>
               </apply>
               <ci>tau_f_L1</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_L_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">27.4</cn>
                        </apply>
                        <cn cellml:units="millivolt">7.1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f_L1</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.161</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">40</cn>
                              </apply>
                              <cn cellml:units="millivolt">14.4</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.01</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f_L2_gate">
      <variable units="dimensionless" public_interface="out" name="f_L_2" initial_value="0.9986"/>
      <variable units="second" name="tau_f_L2"/>
      <variable units="dimensionless" public_interface="in" name="f_L_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f_L_2</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>f_L_infinity</ci>
                  <ci>f_L_2</ci>
               </apply>
               <ci>tau_f_L2</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_f_L2</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">1.3323</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">40</cn>
                              </apply>
                              <cn cellml:units="millivolt">14.2</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.0626</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_independent_transient_outward_K_current">
      <variable units="picoA" public_interface="out" name="i_t"/>
      <variable units="millivolt" public_interface="out" name="E_K"/>
      <variable units="nanoS" name="g_t" initial_value="7.5"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="millimolar" public_interface="in" name="K_c"/>
      <variable units="millimolar" public_interface="in" name="K_i"/>
      <variable units="dimensionless" private_interface="in" name="r"/>
      <variable units="dimensionless" private_interface="in" name="s"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_t</ci>
            <apply>
               <times/>
               <ci>g_t</ci>
               <ci>r</ci>
               <ci>s</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </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>K_c</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_independent_transient_outward_K_current_r_gate">
      <variable units="dimensionless" public_interface="out" name="r" initial_value="0.0010678"/>
      <variable units="second" name="tau_r"/>
      <variable units="dimensionless" name="r_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <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>
         <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/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">1</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">11</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_r</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.0035</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">30</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.0015</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_independent_transient_outward_K_current_s_gate">
      <variable units="dimensionless" public_interface="out" name="s" initial_value="0.949"/>
      <variable units="second" name="tau_s"/>
      <variable units="dimensionless" name="s_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>s_infinity</ci>
                  <ci>s</ci>
               </apply>
               <ci>tau_s</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>s_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">40.5</cn>
                        </apply>
                        <cn cellml:units="millivolt">11.5</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_s</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.4812</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">52.45</cn>
                              </apply>
                              <cn cellml:units="millivolt">14.97</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.01414</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sustained_outward_K_current">
      <variable units="picoA" public_interface="out" name="i_sus"/>
      <variable units="nanoS" name="g_sus" initial_value="2.75"/>
      <variable units="second" 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="dimensionless" private_interface="in" name="r_sus"/>
      <variable units="dimensionless" private_interface="in" name="s_sus"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_sus</ci>
            <apply>
               <times/>
               <ci>g_sus</ci>
               <ci>r_sus</ci>
               <ci>s_sus</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sustained_outward_K_current_r_sus_gate">
      <variable units="dimensionless" public_interface="out" name="r_sus" initial_value="0.00015949"/>
      <variable units="second" name="tau_r_sus"/>
      <variable units="dimensionless" name="r_sus_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r_sus</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>r_sus_infinity</ci>
                  <ci>r_sus</ci>
               </apply>
               <ci>tau_r_sus</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>r_sus_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">4.3</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">8</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_r_sus</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.009</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">5</cn>
                           </apply>
                           <cn cellml:units="millivolt">12</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.0005</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sustained_outward_K_current_s_sus_gate">
      <variable units="dimensionless" public_interface="out" name="s_sus" initial_value="0.9912"/>
      <variable units="second" name="tau_s_sus"/>
      <variable units="dimensionless" name="s_sus_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s_sus</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>s_sus_infinity</ci>
                  <ci>s_sus</ci>
               </apply>
               <ci>tau_s_sus</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>s_sus_infinity</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">0.4</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">20</cn>
                           </apply>
                           <cn cellml:units="millivolt">10</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="dimensionless">0.6</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_s_sus</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.047</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">60</cn>
                           </apply>
                           <cn cellml:units="millivolt">10</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.3</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="delayed_rectifier_K_currents">
      <variable units="picoA" public_interface="out" name="i_Ks"/>
      <variable units="picoA" public_interface="out" name="i_Kr"/>
      <variable units="nanoS" name="g_Ks" initial_value="1"/>
      <variable units="nanoS" name="g_Kr" initial_value="0.5"/>
      <variable units="second" 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="dimensionless" private_interface="in" name="n"/>
      <variable units="dimensionless" private_interface="in" name="p_a"/>
      <variable units="dimensionless" private_interface="in" name="p_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ks</ci>
            <apply>
               <times/>
               <ci>g_Ks</ci>
               <ci>n</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <ci>g_Kr</ci>
               <ci>p_a</ci>
               <ci>p_i</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="delayed_rectifier_K_currents_n_gate">
      <variable units="dimensionless" public_interface="out" name="n" initial_value="0.0048357"/>
      <variable units="second" name="tau_n"/>
      <variable units="dimensionless" name="n_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <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>
         <apply>
            <eq/>
            <ci>n_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">19.9</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">12.7</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_n</ci>
            <apply>
               <plus/>
               <cn cellml:units="second">0.7</cn>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.4</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">20</cn>
                              </apply>
                              <cn cellml:units="millivolt">20</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="delayed_rectifier_K_currents_pa_gate">
      <variable units="dimensionless" public_interface="out" name="p_a" initial_value="0.0001"/>
      <variable units="second" name="tau_p_a"/>
      <variable units="dimensionless" name="p_a_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>p_a</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>p_a_infinity</ci>
                  <ci>p_a</ci>
               </apply>
               <ci>tau_p_a</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>p_a_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">15</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">6</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_p_a</ci>
            <apply>
               <plus/>
               <cn cellml:units="second">0.03118</cn>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.21718</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">20.1376</cn>
                              </apply>
                              <cn cellml:units="millivolt">22.1996</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="delayed_rectifier_K_currents_pi_gate">
      <variable units="dimensionless" public_interface="out" name="p_i"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>p_i</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">55</cn>
                        </apply>
                        <cn cellml:units="millivolt">24</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="inward_rectifier">
      <variable units="picoA" public_interface="out" name="i_K1"/>
      <variable units="nanoS" name="g_K1" initial_value="3"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="millimolar" public_interface="in" name="K_c"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>g_K1</ci>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>K_c</ci>
                        <cn cellml:units="millimolar">1</cn>
                     </apply>
                     <cn cellml:units="dimensionless">0.4457</cn>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>V</ci>
                     <ci>E_K</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">1.5</cn>
                           <apply>
                              <plus/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <ci>E_K</ci>
                              </apply>
                              <cn cellml:units="millivolt">3.6</cn>
                           </apply>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="background_currents">
      <variable units="picoA" public_interface="out" name="i_B_Na"/>
      <variable units="picoA" public_interface="out" name="i_B_Ca"/>
      <variable units="nanoS" name="g_B_Na" initial_value="0.060599"/>
      <variable units="nanoS" name="g_B_Ca" initial_value="0.078681"/>
      <variable units="millivolt" name="E_Ca"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="millimolar" public_interface="in" name="Ca_c"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="second" 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>
         <apply>
            <eq/>
            <ci>i_B_Ca</ci>
            <apply>
               <times/>
               <ci>g_B_Ca</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_c</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable units="picoA" public_interface="out" name="i_NaK"/>
      <variable units="millimolar" name="k_NaK_K" initial_value="1"/>
      <variable units="millimolar" name="k_NaK_Na" initial_value="11"/>
      <variable units="picoA" name="i_NaK_max" initial_value="70.8253"/>
      <variable units="millimolar" public_interface="in" name="K_c"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>i_NaK_max</ci>
                              <ci>K_c</ci>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>K_c</ci>
                              <ci>k_NaK_K</ci>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <cn cellml:units="dimensionless">1.5</cn>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <cn cellml:units="dimensionless">1.5</cn>
                        </apply>
                        <apply>
                           <power/>
                           <ci>k_NaK_Na</ci>
                           <cn cellml:units="dimensionless">1.5</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>V</ci>
                     <cn cellml:units="millivolt">150</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">200</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sarcolemmal_calcium_pump_current">
      <variable units="picoA" public_interface="out" name="i_CaP"/>
      <variable units="picoA" name="i_CaP_max" initial_value="4"/>
      <variable units="millimolar" name="k_CaP" initial_value="0.0002"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_CaP</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_CaP_max</ci>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Ca_i</ci>
                  <ci>k_CaP</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Na_Ca_ion_exchanger_current">
      <variable units="picoA" public_interface="out" name="i_NaCa"/>
      <variable units="picoA_per_millimolar_4" name="k_NaCa" initial_value="0.0374842"/>
      <variable units="per_millimolar_4" name="d_NaCa" initial_value="0.0003"/>
      <variable units="dimensionless" name="gamma" initial_value="0.45"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millimolar" public_interface="in" name="Na_c"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_c"/>
      <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>k_NaCa</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <power/>
                           <ci>Na_i</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Ca_c</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>gamma</ci>
                                 <ci>F</ci>
                                 <ci>V</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <power/>
                           <ci>Na_c</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Ca_i</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>gamma</ci>
                                    <cn cellml:units="dimensionless">1</cn>
                                 </apply>
                                 <ci>V</ci>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <ci>d_NaCa</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <apply>
                              <power/>
                              <ci>Na_c</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                           <ci>Ca_i</ci>
                        </apply>
                        <apply>
                           <times/>
                           <apply>
                              <power/>
                              <ci>Na_i</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                           <ci>Ca_c</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_ion_concentrations">
      <variable units="millimolar" public_interface="out" name="Na_i" initial_value="8.5547"/>
      <variable units="millimolar" public_interface="out" name="Ca_i" initial_value="6.729e-5"/>
      <variable units="millimolar" public_interface="out" name="K_i" initial_value="129.435"/>
      <variable units="millimolar" public_interface="out" name="Ca_d" initial_value="7.2495e-5"/>
      <variable units="picoA" public_interface="out" name="phi_Na_en" initial_value="-1.68"/>
      <variable units="nanolitre" public_interface="out" name="Vol_i" initial_value="0.005884"/>
      <variable units="nanolitre" name="Vol_d"/>
      <variable units="picoA" name="i_di"/>
      <variable units="second" name="tau_di" initial_value="0.01"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="picoA" public_interface="in" name="i_Na"/>
      <variable units="picoA" public_interface="in" name="i_Ca_L"/>
      <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_Ks"/>
      <variable units="picoA" public_interface="in" name="i_Kr"/>
      <variable units="picoA" public_interface="in" name="i_K1"/>
      <variable units="picoA" public_interface="in" name="i_t"/>
      <variable units="picoA" public_interface="in" name="i_sus"/>
      <variable units="picoA" public_interface="in" name="i_CaP"/>
      <variable units="picoA" public_interface="in" name="i_B_Ca"/>
      <variable units="picoA" public_interface="in" name="i_up"/>
      <variable units="picoA" public_interface="in" name="i_rel"/>
      <variable units="per_second" public_interface="in" name="dOCdt"/>
      <variable units="per_second" public_interface="in" name="dOTCdt"/>
      <variable units="per_second" public_interface="in" name="dOTMgCdt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Vol_d</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.02</cn>
               <ci>Vol_i</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Na_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_Na</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>
                     <ci>phi_Na_en</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>Vol_i</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>K_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>i_t</ci>
                        <ci>i_sus</ci>
                        <ci>i_K1</ci>
                        <ci>i_Kr</ci>
                        <ci>i_Ks</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>i_NaK</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>Vol_i</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_i</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <apply>
                                    <minus/>
                                    <ci>i_di</ci>
                                 </apply>
                                 <ci>i_B_Ca</ci>
                                 <ci>i_CaP</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>i_NaCa</ci>
                              </apply>
                           </apply>
                           <ci>i_up</ci>
                        </apply>
                        <ci>i_rel</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>Vol_i</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="millimolar">0.08</cn>
                     <ci>dOTCdt</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millimolar">0.16</cn>
                     <ci>dOTMgCdt</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="millimolar">0.045</cn>
                     <ci>dOCdt</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_d</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_Ca_L</ci>
                     <ci>i_di</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>Vol_d</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_di</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>Ca_d</ci>
                     <ci>Ca_i</ci>
                  </apply>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>F</ci>
                  <ci>Vol_d</ci>
               </apply>
               <ci>tau_di</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_Ca_buffering">
      <variable units="dimensionless" name="O_C" initial_value="0.0275"/>
      <variable units="per_second" public_interface="out" name="dOCdt"/>
      <variable units="dimensionless" name="O_TC" initial_value="0.0133"/>
      <variable units="per_second" public_interface="out" name="dOTCdt"/>
      <variable units="dimensionless" name="O_TMgC" initial_value="0.1961"/>
      <variable units="per_second" public_interface="out" name="dOTMgCdt"/>
      <variable units="dimensionless" name="O_TMgMg" initial_value="0.7094"/>
      <variable units="millimolar" name="Mg_i" initial_value="2.5"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_C</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millimolar_second">200000</cn>
                  <ci>Ca_i</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>O_C</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">476</cn>
                  <ci>O_C</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dOCdt</ci>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_C</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_TC</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millimolar_second">78400</cn>
                  <ci>Ca_i</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>O_TC</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">392</cn>
                  <ci>O_TC</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dOTCdt</ci>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_TC</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_TMgC</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millimolar_second">200000</cn>
                  <ci>Ca_i</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>O_TMgC</ci>
                     </apply>
                     <ci>O_TMgMg</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">6.6</cn>
                  <ci>O_TMgC</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dOTMgCdt</ci>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_TMgC</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_TMgMg</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millimolar_second">2000</cn>
                  <ci>Mg_i</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>O_TMgC</ci>
                     </apply>
                     <ci>O_TMgMg</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">666</cn>
                  <ci>O_TMgMg</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="cleft_space_ion_concentrations">
      <variable units="millimolar" public_interface="out" name="Na_c" initial_value="130.011"/>
      <variable units="millimolar" public_interface="out" name="Ca_c" initial_value="1.8147"/>
      <variable units="millimolar" public_interface="out" name="K_c" initial_value="5.3581"/>
      <variable units="nanolitre" public_interface="in" name="Vol_i"/>
      <variable units="nanolitre" name="Vol_c"/>
      <variable units="second" name="tau_Na" initial_value="14.3"/>
      <variable units="second" name="tau_K" initial_value="10"/>
      <variable units="second" name="tau_Ca" initial_value="24.7"/>
      <variable units="millimolar" name="Na_b" initial_value="130"/>
      <variable units="millimolar" name="Ca_b" initial_value="1.8"/>
      <variable units="millimolar" name="K_b" initial_value="5.4"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="picoA" public_interface="in" name="i_Na"/>
      <variable units="picoA" public_interface="in" name="i_Ca_L"/>
      <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_Ks"/>
      <variable units="picoA" public_interface="in" name="i_Kr"/>
      <variable units="picoA" public_interface="in" name="i_K1"/>
      <variable units="picoA" public_interface="in" name="i_t"/>
      <variable units="picoA" public_interface="in" name="i_sus"/>
      <variable units="picoA" public_interface="in" name="i_CaP"/>
      <variable units="picoA" public_interface="in" name="i_B_Ca"/>
      <variable units="picoA" public_interface="in" name="phi_Na_en"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Vol_c</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.136</cn>
               <ci>Vol_i</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Na_c</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>Na_b</ci>
                     <ci>Na_c</ci>
                  </apply>
                  <ci>tau_Na</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <plus/>
                     <ci>i_Na</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>
                     <ci>phi_Na_en</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Vol_c</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>K_c</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>K_b</ci>
                     <ci>K_c</ci>
                  </apply>
                  <ci>tau_K</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>i_t</ci>
                        <ci>i_sus</ci>
                        <ci>i_K1</ci>
                        <ci>i_Kr</ci>
                        <ci>i_Ks</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>i_NaK</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Vol_c</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_c</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>Ca_b</ci>
                     <ci>Ca_c</ci>
                  </apply>
                  <ci>tau_Ca</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>i_Ca_L</ci>
                        <ci>i_B_Ca</ci>
                        <ci>i_CaP</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>i_NaCa</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>Vol_c</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_handling_by_the_SR">
      <variable units="picoA" public_interface="out" name="i_rel"/>
      <variable units="picoA" public_interface="out" name="i_up"/>
      <variable units="picoA" name="i_tr"/>
      <variable units="picoA" name="I_up_max" initial_value="2800"/>
      <variable units="millimolar" name="k_cyca" initial_value="0.0003"/>
      <variable units="millimolar" name="k_srca" initial_value="0.5"/>
      <variable units="dimensionless" name="k_xcs" initial_value="0.4"/>
      <variable units="picoA_per_millimolar" name="alpha_rel" initial_value="200000"/>
      <variable units="millimolar" name="Ca_rel" initial_value="0.6465"/>
      <variable units="millimolar" name="Ca_up" initial_value="0.6646"/>
      <variable units="nanolitre" name="Vol_up" initial_value="0.0003969"/>
      <variable units="nanolitre" name="Vol_rel" initial_value="4.41e-5"/>
      <variable units="per_second" name="r_act"/>
      <variable units="per_second" name="r_inact"/>
      <variable units="per_second" name="r_recov" initial_value="0.815"/>
      <variable units="dimensionless" name="O_Calse" initial_value="0.4369"/>
      <variable units="dimensionless" name="F1" initial_value="0.4284"/>
      <variable units="dimensionless" name="F2" initial_value="0.0028"/>
      <variable units="second" name="tau_tr" initial_value="0.01"/>
      <variable units="millimolar" name="k_rel_i" initial_value="0.0003"/>
      <variable units="millimolar" name="k_rel_d" initial_value="0.003"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_d"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_up</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>I_up_max</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <ci>k_cyca</ci>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <power/>
                              <ci>k_xcs</ci>
                              <cn cellml:units="dimensionless">2</cn>
                           </apply>
                           <ci>Ca_up</ci>
                        </apply>
                        <ci>k_srca</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>Ca_i</ci>
                        <ci>k_cyca</ci>
                     </apply>
                     <ci>k_cyca</ci>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>k_xcs</ci>
                        <apply>
                           <plus/>
                           <ci>Ca_up</ci>
                           <ci>k_srca</ci>
                        </apply>
                     </apply>
                     <ci>k_srca</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_tr</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>Ca_up</ci>
                     <ci>Ca_rel</ci>
                  </apply>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>F</ci>
                  <ci>Vol_rel</ci>
               </apply>
               <ci>tau_tr</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_rel</ci>
            <apply>
               <times/>
               <ci>alpha_rel</ci>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>F2</ci>
                     <apply>
                        <plus/>
                        <ci>F2</ci>
                        <cn cellml:units="dimensionless">0.25</cn>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>Ca_rel</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_Calse</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <cn cellml:units="per_millimolar_second">480</cn>
                  <ci>Ca_rel</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>O_Calse</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">400</cn>
                  <ci>O_Calse</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_rel</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>i_tr</ci>
                     <ci>i_rel</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>Vol_rel</ci>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="millimolar">31</cn>
                  <apply>
                     <diff/>
                     <bvar>
                        <ci>time</ci>
                     </bvar>
                     <ci>O_Calse</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_up</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>i_up</ci>
                  <ci>i_tr</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>Vol_up</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>F1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>r_recov</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>F1</ci>
                     </apply>
                     <ci>F2</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>r_act</ci>
                  <ci>F1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>F2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>r_act</ci>
                  <ci>F1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>r_inact</ci>
                  <ci>F2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>r_act</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">203.8</cn>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <apply>
                           <plus/>
                           <ci>Ca_i</ci>
                           <ci>k_rel_i</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">4</cn>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Ca_d</ci>
                        <apply>
                           <plus/>
                           <ci>Ca_d</ci>
                           <ci>k_rel_d</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">4</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>r_inact</ci>
            <apply>
               <plus/>
               <cn cellml:units="per_second">33.96</cn>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">339.6</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>Ca_i</ci>
                        <apply>
                           <plus/>
                           <ci>Ca_i</ci>
                           <ci>k_rel_i</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">4</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="sodium_current">
            <component_ref component="sodium_current_m_gate"/>
            <component_ref component="sodium_current_h1_gate"/>
            <component_ref component="sodium_current_h2_gate"/>
         </component_ref>
         <component_ref component="L_type_Ca_channel">
            <component_ref component="L_type_Ca_channel_d_L_gate"/>
            <component_ref component="L_type_Ca_channel_f_L1_gate"/>
            <component_ref component="L_type_Ca_channel_f_L2_gate"/>
         </component_ref>
         <component_ref component="Ca_independent_transient_outward_K_current">
            <component_ref component="Ca_independent_transient_outward_K_current_r_gate"/>
            <component_ref component="Ca_independent_transient_outward_K_current_s_gate"/>
         </component_ref>
         <component_ref component="sustained_outward_K_current">
            <component_ref component="sustained_outward_K_current_r_sus_gate"/>
            <component_ref component="sustained_outward_K_current_s_sus_gate"/>
         </component_ref>
         <component_ref component="delayed_rectifier_K_currents">
            <component_ref component="delayed_rectifier_K_currents_n_gate"/>
            <component_ref component="delayed_rectifier_K_currents_pa_gate"/>
            <component_ref component="delayed_rectifier_K_currents_pi_gate"/>
         </component_ref>
         <component_ref component="inward_rectifier"/>
         <component_ref component="background_currents"/>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="sarcolemmal_calcium_pump_current"/>
         <component_ref component="Na_Ca_ion_exchanger_current"/>
         <component_ref component="Ca_handling_by_the_SR"/>
         <component_ref component="intracellular_ion_concentrations"/>
         <component_ref component="intracellular_Ca_buffering"/>
         <component_ref component="cleft_space_ion_concentrations"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="sodium_current">
         <component_ref component="sodium_current_m_gate"/>
         <component_ref component="sodium_current_h1_gate"/>
         <component_ref component="sodium_current_h2_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
         <component_ref component="L_type_Ca_channel_d_L_gate"/>
         <component_ref component="L_type_Ca_channel_f_L1_gate"/>
         <component_ref component="L_type_Ca_channel_f_L2_gate"/>
      </component_ref>
      <component_ref component="Ca_independent_transient_outward_K_current">
         <component_ref component="Ca_independent_transient_outward_K_current_r_gate"/>
         <component_ref component="Ca_independent_transient_outward_K_current_s_gate"/>
      </component_ref>
      <component_ref component="sustained_outward_K_current">
         <component_ref component="sustained_outward_K_current_r_sus_gate"/>
         <component_ref component="sustained_outward_K_current_s_sus_gate"/>
      </component_ref>
      <component_ref component="delayed_rectifier_K_currents">
         <component_ref component="delayed_rectifier_K_currents_n_gate"/>
         <component_ref component="delayed_rectifier_K_currents_pa_gate"/>
         <component_ref component="delayed_rectifier_K_currents_pi_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="sodium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sustained_outward_K_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="delayed_rectifier_K_currents"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="inward_rectifier"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="background_currents"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sarcolemmal_calcium_pump_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="Na_Ca_ion_exchanger_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="Ca_handling_by_the_SR"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="intracellular_Ca_buffering"/>
      <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="cleft_space_ion_concentrations"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ca_L" variable_1="i_Ca_L"/>
   </connection>
   <connection>
      <map_components component_2="Ca_independent_transient_outward_K_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_t" variable_1="i_t"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="sustained_outward_K_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_sus" variable_1="i_sus"/>
   </connection>
   <connection>
      <map_components component_2="inward_rectifier" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K1" variable_1="i_K1"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="delayed_rectifier_K_currents" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
   </connection>
   <connection>
      <map_components component_2="background_currents" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_B_Na" variable_1="i_B_Na"/>
      <map_variables variable_2="i_B_Ca" variable_1="i_B_Ca"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="sodium_potassium_pump" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
   </connection>
   <connection>
      <map_components component_2="sarcolemmal_calcium_pump_current" component_1="membrane"/>
      <map_variables variable_2="i_CaP" variable_1="i_CaP"/>
   </connection>
   <connection>
      <map_components component_2="Na_Ca_ion_exchanger_current" component_1="membrane"/>
      <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="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="membrane"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="membrane"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="Ca_handling_by_the_SR" component_1="membrane"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sodium_current"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
   </connection>
   <connection>
      <map_components component_2="background_currents" component_1="sodium_current"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="sodium_current"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
      <map_variables variable_2="Na_c" variable_1="Na_c"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="i_Ca_L" variable_1="i_Ca_L"/>
      <map_variables variable_2="Ca_d" variable_1="Ca_d"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="i_Ca_L" variable_1="i_Ca_L"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="i_t" variable_1="i_t"/>
      <map_variables variable_2="K_i" variable_1="K_i"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="i_t" variable_1="i_t"/>
      <map_variables variable_2="K_c" variable_1="K_c"/>
   </connection>
   <connection>
      <map_components component_2="sustained_outward_K_current" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="inward_rectifier" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="delayed_rectifier_K_currents" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sustained_outward_K_current"/>
      <map_variables variable_2="i_sus" variable_1="i_sus"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="sustained_outward_K_current"/>
      <map_variables variable_2="i_sus" variable_1="i_sus"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="delayed_rectifier_K_currents"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="delayed_rectifier_K_currents"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="inward_rectifier"/>
      <map_variables variable_2="i_K1" variable_1="i_K1"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="inward_rectifier"/>
      <map_variables variable_2="i_K1" variable_1="i_K1"/>
      <map_variables variable_2="K_c" variable_1="K_c"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="background_currents"/>
      <map_variables variable_2="i_B_Na" variable_1="i_B_Na"/>
      <map_variables variable_2="i_B_Ca" variable_1="i_B_Ca"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="background_currents"/>
      <map_variables variable_2="i_B_Na" variable_1="i_B_Na"/>
      <map_variables variable_2="i_B_Ca" variable_1="i_B_Ca"/>
      <map_variables variable_2="Ca_c" variable_1="Ca_c"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
      <map_variables variable_2="K_c" variable_1="K_c"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sarcolemmal_calcium_pump_current"/>
      <map_variables variable_2="i_CaP" variable_1="i_CaP"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="sarcolemmal_calcium_pump_current"/>
      <map_variables variable_2="i_CaP" variable_1="i_CaP"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="Na_Ca_ion_exchanger_current"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_2="cleft_space_ion_concentrations" component_1="Na_Ca_ion_exchanger_current"/>
      <map_variables variable_2="Ca_c" variable_1="Ca_c"/>
      <map_variables variable_2="Na_c" variable_1="Na_c"/>
      <map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="Ca_handling_by_the_SR"/>
      <map_variables variable_2="i_rel" variable_1="i_rel"/>
      <map_variables variable_2="i_up" variable_1="i_up"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Ca_d" variable_1="Ca_d"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="cleft_space_ion_concentrations"/>
      <map_variables variable_2="phi_Na_en" variable_1="phi_Na_en"/>
      <map_variables variable_2="Vol_i" variable_1="Vol_i"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="intracellular_Ca_buffering"/>
      <map_variables variable_2="dOCdt" variable_1="dOCdt"/>
      <map_variables variable_2="dOTCdt" variable_1="dOTCdt"/>
      <map_variables variable_2="dOTMgCdt" variable_1="dOTMgCdt"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current_m_gate" component_1="sodium_current"/>
      <map_variables variable_2="m" variable_1="m"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current_h1_gate" component_1="sodium_current"/>
      <map_variables variable_2="h1" variable_1="h1"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current_h2_gate" component_1="sodium_current"/>
      <map_variables variable_2="h2" variable_1="h2"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current_h2_gate" component_1="sodium_current_h1_gate"/>
      <map_variables variable_2="h_infinity" variable_1="h_infinity"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_d_L_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="d_L" variable_1="d_L"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_f_L1_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="f_L_1" variable_1="f_L_1"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_f_L2_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="f_L_2" variable_1="f_L_2"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel_f_L2_gate" component_1="L_type_Ca_channel_f_L1_gate"/>
      <map_variables variable_2="f_L_infinity" variable_1="f_L_infinity"/>
   </connection>
   <connection>
      <map_components component_2="Ca_independent_transient_outward_K_current_r_gate" component_1="Ca_independent_transient_outward_K_current"/>
      <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="Ca_independent_transient_outward_K_current_s_gate" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="s" variable_1="s"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sustained_outward_K_current_r_sus_gate" component_1="sustained_outward_K_current"/>
      <map_variables variable_2="r_sus" variable_1="r_sus"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="sustained_outward_K_current_s_sus_gate" component_1="sustained_outward_K_current"/>
      <map_variables variable_2="s_sus" variable_1="s_sus"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="delayed_rectifier_K_currents_n_gate" component_1="delayed_rectifier_K_currents"/>
      <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="delayed_rectifier_K_currents_pa_gate" component_1="delayed_rectifier_K_currents"/>
      <map_variables variable_2="p_a" variable_1="p_a"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="delayed_rectifier_K_currents_pi_gate" component_1="delayed_rectifier_K_currents"/>
      <map_variables variable_2="p_i" variable_1="p_i"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>



<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#57a005cd-4723-485f-9953-3df032689207">
    <rdf:li rdf:resource="rdf:#626f0307-5bfa-4f7d-839f-0aaa6fc41c14"/>
    <rdf:li rdf:resource="rdf:#35eb2d78-e3dd-499e-9f1a-bf905737d68e"/>
    <rdf:li rdf:resource="rdf:#7c02ad2b-d7a2-4038-8ffe-16e34155cc7f"/>
    <rdf:li rdf:resource="rdf:#13eabd73-eacb-4de6-87b5-6bbbfd9641a7"/>
    <rdf:li rdf:resource="rdf:#07b76e21-f158-4fc0-ade9-159fb701b076"/>
    <rdf:li rdf:resource="rdf:#285aad83-b787-415e-9ec2-354a0670ec77"/>
    <rdf:li rdf:resource="rdf:#f3220b55-2e34-4a5e-9acf-95d82533df86"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#285aad83-b787-415e-9ec2-354a0670ec77">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#48869153-cc32-4475-b10f-9e3768192ba1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#79ea5fbc-292d-425f-87cc-e7a21c787865">
    <bqs:Pubmed_id>9440706</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#f3212125-87a8-4b25-a599-2c9da0dddf02"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Auckland Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#dcf648ec-6e68-46e9-9bec-d2c6f1cbe368"/>
    <dcterms:created rdf:resource="rdf:#16f2764e-66b0-462e-a724-8a9c972d8941"/>
    <dc:creator rdf:resource="rdf:#c561bcbd-7e00-41ad-a1d1-fb246d737870"/>
    <cmeta:modification rdf:resource="rdf:#ffb6683a-3b3a-4813-95e5-3f63e81370d8"/>
    <cmeta:modification rdf:resource="rdf:#77900f89-38c0-418d-a923-693b287b8194"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#48869153-cc32-4475-b10f-9e3768192ba1">
    <vCard:Given>R</vCard:Given>
    <vCard:Family>Clark</vCard:Family>
    <vCard:Other>B</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#07b76e21-f158-4fc0-ade9-159fb701b076">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#dbef145c-4ab5-4558-a8a1-d3a2279b8e58"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#caaf9885-8593-467a-a9b4-db00f1beca95">
    <vCard:Given>A</vCard:Given>
    <vCard:Family>Nygren</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#77900f89-38c0-418d-a923-693b287b8194">
    <dcterms:modified rdf:resource="rdf:#3f958671-f695-4a49-bb85-2b07f88d2121"/>
    <rdf:value>Dr. Daniel Ziemek noted that the model had two bugs in it. Those were, that the variable gCaL should be 6.75 instead of 3, and also that the stimulus amplitude should have been -280 pA instead of 400. These values were changed by Penny Noble, who sent James Lawson the file. James Lawson then made a few modifications to allow the file to be run in PCEnv as well as COR.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#bb57a8db-84c7-4611-bc07-65426853566c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#60c88f97-b0fd-4aec-82be-062e4940876f">
    <dcterms:W3CDTF>2006-03-31</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$8zT9U1">
    <ns7:boundIntervals rdf:resource="rdf:#$9zT9U1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4da186e3-cb6a-4dbb-8414-c46dcb500f51">
    <vCard:Given>W</vCard:Given>
    <vCard:Family>Giles</vCard:Family>
    <vCard:Other>R</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aede1c0c-6595-4ef9-9838-7c4ceb81a555">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>unknown</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#488e08a8-089a-4719-8320-93f001c3d366">
    <vCard:N rdf:resource="rdf:#0794bd8d-ce21-488f-9ac3-07bc84330de9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$azT9U1">
    <ns7:nonstandard-pointDensity>100000</ns7:nonstandard-pointDensity>
    <ns7:maximumStepSize>0.001</ns7:maximumStepSize>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#13006808-a969-4597-bf02-830b4414a57a">
    <vCard:N rdf:resource="rdf:#1eedce6f-2993-468e-8781-ec471834c270"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#19e3b15e-7b94-49ac-96c3-c0ee1974b7f8">
    <vCard:Given>J</vCard:Given>
    <vCard:Family>Clark</vCard:Family>
    <vCard:Other>W</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0d7ca1ab-10c1-41c3-a8c6-2c395cc312cf">
    <dc:creator rdf:resource="rdf:#3c10cda4-2b07-4f7f-a86a-8bc3a5327dd7"/>
    <rdf:value>This is the CellML description of Nygren et al's mathematical model of an adult human atrial cell.  They specifically study the role of K+ currents in repolarisation.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3c10cda4-2b07-4f7f-a86a-8bc3a5327dd7">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$9zT9U1">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$azT9U1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#16f2764e-66b0-462e-a724-8a9c972d8941">
    <dcterms:W3CDTF>2002-05-23T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b527a6b0-9b2d-429a-8f01-e2d23e5699d6">
    <dcterms:W3CDTF>2007-04-03</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3f958671-f695-4a49-bb85-2b07f88d2121">
    <dcterms:W3CDTF>2007-08-17T13:21:40+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bb76aeb3-65b2-4797-a8a8-1b1837c5962d">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#13eabd73-eacb-4de6-87b5-6bbbfd9641a7">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#19e3b15e-7b94-49ac-96c3-c0ee1974b7f8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#69bc2a9b-45f4-4721-b9a7-490af5afb098">
    <vCard:Given>James</vCard:Given>
    <vCard:Family>Lawson</vCard:Family>
    <vCard:Other>Richard</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b0f0d165-7554-462d-bd7f-0ef13c352365">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bb57a8db-84c7-4611-bc07-65426853566c">
    <vCard:N rdf:resource="rdf:#69bc2a9b-45f4-4721-b9a7-490af5afb098"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#85d66125-278c-416e-8668-fa6030c53ac1">
    <dc:creator rdf:resource="rdf:#7224713d-7360-405a-9f17-76261dfe470f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5e56c6a4-87fc-4516-a3b3-a72a4a11c66e">
    <vCard:Given>L</vCard:Given>
    <vCard:Family>Firek</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#05ddb03f-e086-413a-90f6-ab63087ab0b0">
    <vCard:FN>James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#626f0307-5bfa-4f7d-839f-0aaa6fc41c14">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#caaf9885-8593-467a-a9b4-db00f1beca95"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dfc6fc1f-b3ac-476c-b452-8e9c128f683c">
    <dcterms:modified rdf:resource="rdf:#60c88f97-b0fd-4aec-82be-062e4940876f"/>
    <rdf:value>added metadata</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#488e08a8-089a-4719-8320-93f001c3d366"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dcf648ec-6e68-46e9-9bec-d2c6f1cbe368">
    <dc:creator rdf:resource="rdf:#05ddb03f-e086-413a-90f6-ab63087ab0b0"/>
    <rdf:value>This model has been curated by Penny Noble and James Lawson and is known to run and produce the correct output in PCEnv and COR.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7c02ad2b-d7a2-4038-8ffe-16e34155cc7f">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#5e56c6a4-87fc-4516-a3b3-a72a4a11c66e"/>
  </rdf:Description>
  <rdf:Description rdf:about="#nygren_fiset_firek_clark_lindblad_clark_giles_1998">
    <ns7:simulation rdf:resource="rdf:#$8zT9U1"/>
    <bqs:reference rdf:resource="rdf:#79ea5fbc-292d-425f-87cc-e7a21c787865"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>k+ currents</rdf:li>
        <rdf:li>Atrial Myocyte</rdf:li>
        <rdf:li>atrial myocyte</rdf:li>
        <rdf:li>electrophysiology</rdf:li>
        <rdf:li>human</rdf:li>
        <rdf:li>cardiac</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <cmeta:comment rdf:resource="rdf:#0d7ca1ab-10c1-41c3-a8c6-2c395cc312cf"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f3212125-87a8-4b25-a599-2c9da0dddf02">
    <dc:creator rdf:resource="rdf:#57a005cd-4723-485f-9953-3df032689207"/>
    <dc:title>Mathematical Model of an Adult Human Atrial Cell: The Role of K+ Currents in Repolarization</dc:title>
    <bqs:volume>82</bqs:volume>
    <bqs:first_page>63</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#656c98f7-8ec6-427e-a9d3-9c1030c8f31b"/>
    <dcterms:issued rdf:resource="rdf:#1a65eb5a-90e4-4621-a6df-e49d16828515"/>
    <bqs:last_page>81</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#567cc0a5-77ad-4d73-95ed-8f8623b832a3">
    <dcterms:W3CDTF>2008-02-25T11:36:41+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#35eb2d78-e3dd-499e-9f1a-bf905737d68e">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#95b7576b-cf1d-4993-b720-a639139389cc"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0794bd8d-ce21-488f-9ac3-07bc84330de9">
    <vCard:Given>Peter</vCard:Given>
    <vCard:Family>Villiger</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#95b7576b-cf1d-4993-b720-a639139389cc">
    <vCard:Given>C</vCard:Given>
    <vCard:Family>Fiset</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c561bcbd-7e00-41ad-a1d1-fb246d737870">
    <vCard:ORG rdf:resource="rdf:#b0f0d165-7554-462d-bd7f-0ef13c352365"/>
    <vCard:EMAIL rdf:resource="rdf:#bb76aeb3-65b2-4797-a8a8-1b1837c5962d"/>
    <vCard:N rdf:resource="rdf:#6d762b96-d526-4072-9c09-ff187782fd6b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f3220b55-2e34-4a5e-9acf-95d82533df86">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#4da186e3-cb6a-4dbb-8414-c46dcb500f51"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <cmeta:modification rdf:resource="rdf:#dfc6fc1f-b3ac-476c-b452-8e9c128f683c"/>
    <dcterms:created rdf:resource="rdf:#b527a6b0-9b2d-429a-8f01-e2d23e5699d6"/>
    <dc:creator rdf:resource="rdf:#1ad72ab7-2286-44dd-adf9-4c6c215af53d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1a65eb5a-90e4-4621-a6df-e49d16828515">
    <dcterms:W3CDTF>1998-01-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1ad72ab7-2286-44dd-adf9-4c6c215af53d">
    <vCard:ORG rdf:resource="rdf:#548bcd7c-4ad0-44e1-9bd6-3859bc22748c"/>
    <vCard:EMAIL rdf:resource="rdf:#aede1c0c-6595-4ef9-9838-7c4ceb81a555"/>
    <vCard:N rdf:resource="rdf:#df29d364-fb6c-4335-abe1-fe926face5a4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#656c98f7-8ec6-427e-a9d3-9c1030c8f31b">
    <dc:title>Circulation Research</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dbef145c-4ab5-4558-a8a1-d3a2279b8e58">
    <vCard:Given>D</vCard:Given>
    <vCard:Family>Lindblad</vCard:Family>
    <vCard:Other>S</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ffb6683a-3b3a-4813-95e5-3f63e81370d8">
    <dcterms:modified rdf:resource="rdf:#567cc0a5-77ad-4d73-95ed-8f8623b832a3"/>
    <rdf:value>Units checked, curated</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#13006808-a969-4597-bf02-830b4414a57a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6d762b96-d526-4072-9c09-ff187782fd6b">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1eedce6f-2993-468e-8781-ec471834c270">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
</rdf:RDF>
</model>