Location: Michailova 2005 @ 6c54bb01d941 / michailova_2005.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2010-01-20 14:37:27+13:00
Desc:
dded katp images in ai svg and png format
Permanent Source URI:
https://staging.physiomeproject.org/workspace/michailova_saucerman_belik_mcculloch_2005/rawfile/6c54bb01d941abfefb85930761320d246d74d99d/michailova_2005.cellml

<?xml version='1.0' encoding='utf-8'?>

<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="michailova_2005" name="michailova_2005">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modeling regulation of cardiac KATP and L-type Ca2+ currents by ATP, ADP, and Mg2+</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
          This CellML model runs in both COR and PCEnv.  The units have been checked and they are consistent.  Unfortunately the CellML model is still not quite reproducing the published results.  We will continue to curate the model and hope to get it fully working.  
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: Changes in cytosolic free Mg(2+) and adenosine nucleotide phosphates affect cardiac excitability and contractility. To investigate how modulation by Mg(2+), ATP, and ADP of K(ATP) and L-type Ca(2+) channels influences excitation-contraction coupling, we incorporated equations for intracellular ATP and MgADP regulation of the K(ATP) current and MgATP regulation of the L-type Ca(2+) current in an ionic-metabolic model of the canine ventricular myocyte. The new model: 1), quantitatively reproduces a dose-response relationship for the effects of changes in ATP on K(ATP) current, 2), simulates effects of ADP in modulating ATP sensitivity of K(ATP) channel, 3), predicts activation of Ca(2+) current during rapid increase in MgATP, and 4), demonstrates that decreased ATP/ADP ratio with normal total Mg(2+) or increased free Mg(2+) with normal ATP and ADP activate K(ATP) current, shorten action potential, and alter ionic currents and intracellular Ca(2+) signals. The model predictions are in agreement with experimental data measured under normal and a variety of pathological conditions.
</para>

<informalfigure float="0" id="katp_channel_diagram">
<mediaobject>
  <imageobject>
    <imagedata fileref="michailova_2005_katp.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic representation of Kir6.2 and SUR2A subunit stoichiometry in cardiac octameric KATP channel. It is assumed that the single ATP molecule is sufficient to close the channel and that the binding of 2MgADP molecules to one SUR2A subunit increases the channel open probability.</caption>
</informalfigure>

<informalfigure float="0" id="cell_diagram">
<mediaobject>
  <imageobject>
    <imagedata fileref="michailova_2005_fullcell.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the full cell model illustrating Ca2+ and Mg2+ buffering and transport by ATP and ADP, adenine nucleotides regulation of ionic channels and pumps, and electrophysiology in a ventricular myocyte.</caption>
</informalfigure>

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

<para>
Modeling Regulation of Cardiac KATP and L-type Ca2+ Currents by ATP, ADP, and Mg2+, Anushka Michailova, Jeffrey Saucerman, Mary Ellen Belik, and Andrew D. McCulloch, 2005, <emphasis>Biophysical Journal</emphasis>. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/15738467">PubMed ID: 18339738</ulink>
</para>

</sect1>
</article>
</documentation>

   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="second_per_millisecond">
      <unit units="second"/>
      <unit units="second" prefix="milli" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="per_millivolt">
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="per_millivolt_second">
      <unit units="millivolt" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="milliS_per_microF">
      <unit units="siemens" prefix="milli"/>
      <unit units="farad" prefix="micro" exponent="-1"/>
   </units>
   <units name="microF_per_cm2">
      <unit units="farad" prefix="micro"/>
      <unit units="metre" prefix="centi" exponent="-2"/>
   </units>
   <units name="microA_per_microF">
      <unit units="ampere" prefix="micro"/>
      <unit units="farad" prefix="micro" exponent="-1"/>
   </units>
   <units name="nA_microA">
      <unit units="ampere" prefix="nano"/>
      <unit units="ampere" prefix="micro" exponent="-1"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="millimolar_per_second">
      <unit units="millimolar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_millimolar_second">
      <unit units="millimolar" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millimolar3_per_second">
      <unit units="millimolar" exponent="-3"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millimolar4_per_second">
      <unit units="millimolar" exponent="-4"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="joule_per_mole_kelvin">
      <unit units="joule"/>
      <unit units="mole" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   <units name="coulomb_per_millimole">
      <unit units="coulomb"/>
      <unit units="mole" prefix="milli" exponent="-1"/>
   </units>
   <units name="cm_per_second">
      <unit units="metre" prefix="centi"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="cm_per_farad_second">
      <unit units="metre" prefix="centi"/>
      <unit units="farad" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="cm2">
      <unit units="metre" prefix="centi" exponent="2"/>
   </units>
   <units name="microlitre">
      <unit units="litre" prefix="micro"/>
   </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="-96.1638"/>
      <variable units="joule_per_mole_kelvin" public_interface="out" name="R" initial_value="8.314472"/>
      <variable units="kelvin" public_interface="out" name="T" initial_value="310"/>
      <variable units="coulomb_per_millimole" public_interface="out" name="F" initial_value="96.4853415"/>
      <variable units="microF_per_cm2" public_interface="out" name="C_sc" initial_value="0.001"/>
	  <variable units="microA_per_microF" public_interface="out" name="i_Stim"/>
      
      <variable units="second" public_interface="in" name="time"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Na"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca_K"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Kr"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ks"/>
      <variable units="microA_per_microF" public_interface="in" name="i_to1"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K1"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Kp"/>
      <variable units="microA_per_microF" public_interface="in" name="i_NaCa"/>
      <variable units="microA_per_microF" public_interface="in" name="i_NaK"/>
      <variable units="microA_per_microF" public_interface="in" name="i_p_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca_b"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Na_b"/>
	  <variable units="microA_per_microF" public_interface="in" name="i_K_ATP"/>
      
	  <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.002"/>
      <variable units="microA_per_microF" name="stim_amplitude" initial_value="-21.1268"/>
	  
      <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="microA_per_microF">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="microF_per_cm2">1</cn>
                  </apply>
                  <cn cellml:units="nA_microA">1</cn>
                  <apply>
                     <plus/>
                     <ci>i_Na</ci>
                     <ci>i_Ca</ci>
                     <ci>i_Ca_K</ci>
                     <ci>i_Kr</ci>
                     <ci>i_Ks</ci>
                     <ci>i_to1</ci>
                     <ci>i_K1</ci>
                     <ci>i_Kp</ci>
                     <ci>i_NaCa</ci>
                     <ci>i_NaK</ci>
                     <ci>i_p_Ca</ci>
                     <ci>i_Na_b</ci>
                     <ci>i_Ca_b</ci>
					 <ci>i_K_ATP</ci>
                     <ci>i_Stim</ci>
                  </apply>
               </apply>
               <ci>C_sc</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="fast_sodium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Na"/>
      <variable units="millivolt" public_interface="out" name="E_Na"/>
      <variable units="milliS_per_microF" name="g_Na" initial_value="12.8"/>
      <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_o"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_millimole" 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="h"/>
      <variable units="dimensionless" private_interface="in" name="j"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Na_o</ci>
                     <ci>Na_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <ci>j</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.0328302"/>
      <variable units="per_second" name="alpha_m"/>
      <variable units="per_second" name="beta_m"/>
      <variable units="millivolt" name="E0_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/>
            <ci>E0_m</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn cellml:units="millivolt">47.13</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <cn cellml:units="per_second">1000</cn>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">0.1</cn>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millivolt">0.005</cn>
                           <ci>E0_m</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_m</ci>
                     </apply>
                     <cn cellml:units="millivolt">0.00001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt_second">320</cn>
                        <ci>E0_m</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1</cn>
                              </apply>
                              <ci>E0_m</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">80</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">11</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>alpha_m</ci>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>m</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>beta_m</ci>
                        <ci>m</ci>
                     </apply>
                  </apply>
                  <apply>
                     <geq/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">90</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="per_second">0</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h_gate">
      <variable units="dimensionless" public_interface="out" name="h" initial_value="0.988354"/>
      <variable units="per_second" name="alpha_h"/>
      <variable units="per_second" name="beta_h"/>
      <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>alpha_h</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <cn cellml:units="per_second">135</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <cn cellml:units="millivolt">80</cn>
                              <ci>V</ci>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">6.8</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="per_second">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <piecewise>
               <piece>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_second">3560</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.079</cn>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="per_second">310000</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.35</cn>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1000</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="second">0.13</cn>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">10.66</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="millivolt">11.1</cn>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_h</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>h</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_h</ci>
                  <ci>h</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_j_gate">
      <variable units="dimensionless" public_interface="out" name="j" initial_value="0.99254"/>
      <variable units="per_second" name="alpha_j"/>
      <variable units="per_second" name="beta_j"/>
      <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>alpha_j</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">1000</cn>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="per_millivolt_second">127140</cn>
                                 <apply>
                                    <exp/>
                                    <apply>
                                       <times/>
                                       <cn cellml:units="per_millivolt">0.2444</cn>
                                       <ci>V</ci>
                                    </apply>
                                 </apply>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="per_millivolt_second">0.00003474</cn>
                                 <apply>
                                    <exp/>
                                    <apply>
                                       <times/>
                                       <apply>
                                          <minus/>
                                          <cn cellml:units="per_millivolt">0.04391</cn>
                                       </apply>
                                       <ci>V</ci>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">37.78</cn>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.311</cn>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">79.23</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="per_second">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_j</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_second">121.2</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.01052</cn>
                              </apply>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1378</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">40.14</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_second">300</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.0000002535</cn>
                              </apply>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">32</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>j</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_j</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>j</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_j</ci>
                  <ci>j</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_activating_delayed_rectifiyer_K_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Kr"/>
      <variable units="millivolt" public_interface="out" name="E_K"/>
      <variable units="milliS_per_microF" name="g_Kr" initial_value="0.0034"/>
      <variable units="dimensionless" name="f_K_o"/>
      <variable units="dimensionless" name="R_V"/>
      <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="K_o"/>
      <variable units="millimolar" public_interface="in" name="K_i"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="dimensionless" private_interface="in" name="X_kr"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <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_o</ci>
                     <ci>K_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_K_o</ci>
            <apply>
               <root/>
               <apply>
                  <divide/>
                  <ci>K_o</ci>
                  <cn cellml:units="millimolar">4</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>R_V</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">1.4945</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millivolt">0.0446</cn>
                           <ci>V</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Kr</ci>
            <apply>
               <times/>
               <ci>g_Kr</ci>
               <ci>f_K_o</ci>
               <ci>R_V</ci>
               <ci>X_kr</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rapid_activating_delayed_rectifiyer_K_current_X_kr_gate">
      <variable units="dimensionless" public_interface="out" name="X_kr" initial_value="0.51"/>
      <variable units="dimensionless" name="K12"/>
      <variable units="dimensionless" name="K21"/>
      <variable units="dimensionless" name="X_kr_inf"/>
      <variable units="second" name="tau_X_kr"/>
      <variable units="dimensionless" name="tau_factor" initial_value="1"/>
      <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>K12</ci>
            <apply>
               <exp/>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">5.495</cn>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt">0.1691</cn>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K21</ci>
            <apply>
               <exp/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">7.677</cn>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt">0.0128</cn>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>X_kr_inf</ci>
            <apply>
               <divide/>
               <ci>K12</ci>
               <apply>
                  <plus/>
                  <ci>K12</ci>
                  <ci>K21</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_X_kr</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="second">0.001</cn>
                  <apply>
                     <plus/>
                     <ci>K12</ci>
                     <ci>K21</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>tau_factor</ci>
                  <cn cellml:units="second">0.027</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>X_kr</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>X_kr_inf</ci>
                  <ci>X_kr</ci>
               </apply>
               <ci>tau_X_kr</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_activating_delayed_rectifiyer_K_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Ks"/>
      <variable units="milliS_per_microF" name="g_Ks" initial_value="0.0027134"/>
      <variable units="millivolt" name="E_Ks"/>
      <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="K_o"/>
      <variable units="millimolar" public_interface="in" name="K_i"/>
      <variable units="millimolar" public_interface="in" name="Na_o"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="dimensionless" private_interface="in" name="X_ks"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Ks</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>K_o</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.01833</cn>
                           <ci>Na_o</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_i</ci>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.01833</cn>
                           <ci>Na_i</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ks</ci>
            <apply>
               <times/>
               <ci>g_Ks</ci>
               <apply>
                  <power/>
                  <ci>X_ks</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ks</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_activating_delayed_rectifiyer_K_current_X_ks_gate">
      <variable units="dimensionless" public_interface="out" name="X_ks" initial_value="0.264"/>
      <variable units="second" name="tau_X_ks"/>
      <variable units="dimensionless" name="X_ks_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/>
            <ci>X_ks_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">24.7</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">13.6</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_X_ks</ci>
            <apply>
               <divide/>
               <cn cellml:units="second">0.001</cn>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.0000719</cn>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">10</cn>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.148</cn>
                              </apply>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">10</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.000131</cn>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">10</cn>
                        </apply>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.0687</cn>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">10</cn>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>X_ks</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>X_ks_infinity</ci>
                  <ci>X_ks</ci>
               </apply>
               <ci>tau_X_ks</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_to1"/>
      <variable units="milliS_per_microF" name="g_to1" initial_value="0.23815"/>
      <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="X_to1"/>
      <variable units="dimensionless" private_interface="in" name="Y_to1"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_to1</ci>
            <apply>
               <times/>
               <ci>g_to1</ci>
               <ci>X_to1</ci>
               <ci>Y_to1</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_potassium_current_X_to1_gate">
      <variable units="dimensionless" public_interface="out" name="X_to1" initial_value="2.63"/>
      <variable units="per_second" name="alpha_X_to1"/>
      <variable units="per_second" name="beta_X_to1"/>
      <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>alpha_X_to1</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">45.16</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt">0.03577</cn>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_X_to1</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">98.9</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="per_millivolt">0.06237</cn>
                     </apply>
                     <ci>V</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>X_to1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_X_to1</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>X_to1</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_X_to1</ci>
                  <ci>X_to1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_potassium_current_Y_to1_gate">
      <variable units="dimensionless" public_interface="out" name="Y_to1" initial_value="0.99"/>
      <variable units="per_second" name="alpha_Y_to1"/>
      <variable units="per_second" name="beta_Y_to1"/>
      <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>alpha_Y_to1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">5.415</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">33.5</cn>
                           </apply>
                        </apply>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.051335</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">33.5</cn>
                              </apply>
                           </apply>
                           <cn cellml:units="millivolt">5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_Y_to1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">5.415</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">33.5</cn>
                        </apply>
                        <cn cellml:units="millivolt">5</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.051335</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">33.5</cn>
                           </apply>
                           <cn cellml:units="millivolt">5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Y_to1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_Y_to1</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>Y_to1</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_Y_to1</ci>
                  <ci>Y_to1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_K1"/>
      <variable units="milliS_per_microF" name="g_K1" initial_value="2.8"/>
      <variable units="millimolar" name="K_mK1" initial_value="13"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" private_interface="out" name="R"/>
      <variable units="kelvin" public_interface="in" private_interface="out" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" private_interface="out" name="F"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="E_K"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
      <variable units="dimensionless" private_interface="in" name="K1_infinity_V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>g_K1</ci>
                     <ci>K1_infinity_V</ci>
                     <ci>K_o</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>K_o</ci>
                     <ci>K_mK1</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current_K1_gate">
      <variable units="dimensionless" public_interface="out" name="K1_infinity_V"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>K1_infinity_V</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">2</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">1.5</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <ci>E_K</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="plateau_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Kp"/>
      <variable units="milliS_per_microF" name="g_Kp" initial_value="0.002216"/>
      <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="Kp_V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Kp</ci>
            <apply>
               <times/>
               <ci>g_Kp</ci>
               <ci>Kp_V</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="plateau_potassium_current_Kp_gate">
      <variable units="dimensionless" public_interface="out" name="Kp_V"/>
      <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>Kp_V</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">7.488</cn>
                           <ci>V</ci>
                        </apply>
                        <cn cellml:units="millivolt">5.98</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Na_Ca_exchanger">
      <variable units="microA_per_microF" public_interface="out" name="i_NaCa"/>
      <variable units="millimolar" name="K_mCa" initial_value="1.38"/>
      <variable units="millimolar" name="K_mNa" initial_value="87.5"/>
      <variable units="microA_per_microF" name="K_NaCa" initial_value="0.3"/>
      <variable units="dimensionless" name="K_sat" initial_value="0.2"/>
      <variable units="dimensionless" name="eta" initial_value="0.35"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millimolar" public_interface="in" name="Na_o"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_o"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>K_NaCa</ci>
                     <cn cellml:units="dimensionless">5000</cn>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <apply>
                           <power/>
                           <ci>K_mNa</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Na_o</ci>
                           <cn cellml:units="dimensionless">3</cn>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_mCa</ci>
                        <ci>Ca_o</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>K_sat</ci>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <times/>
                                    <apply>
                                       <minus/>
                                       <ci>eta</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>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>eta</ci>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <power/>
                        <ci>Na_i</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                     <ci>Ca_o</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>eta</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>
                        <power/>
                        <ci>Na_o</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="sodium_potassium_pump">
      <variable units="microA_per_microF" public_interface="out" name="i_NaK"/>
	  
	  <variable units="microA_per_microF" name="I_NaK" initial_value="0.693"/>
      <variable units="dimensionless" name="f_NaK"/>
      <variable units="millimolar" name="K_mNa_i" initial_value="10"/>
      <variable units="millimolar" name="K_mK_o" initial_value="1.5"/>
      <variable units="dimensionless" name="sigma"/>
	  
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Na_i"/>
      <variable units="millimolar" public_interface="in" name="Na_o"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_NaK</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.1245</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">0.1</cn>
                              </apply>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0365</cn>
                     <ci>sigma</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>sigma</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <cn cellml:units="dimensionless">7</cn>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <ci>Na_o</ci>
                        <cn cellml:units="millimolar">67.3</cn>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>i_NaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>I_NaK</ci>
                        <ci>f_NaK</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <ci>K_mNa_i</ci>
                              <ci>Na_i</ci>
                           </apply>
                           <cn cellml:units="dimensionless">1.5</cn>
                        </apply>
                     </apply>
                  </apply>
                  <ci>K_o</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_o</ci>
                  <ci>K_mK_o</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="sarcolemmal_calcium_pump">
      <variable units="microA_per_microF" public_interface="out" name="i_p_Ca"/>
      
	  <variable units="microA_per_microF" name="i_p_Ca_winslow"/>
	  <variable units="millimolar" name="K_mpCa" initial_value="0.00005"/>
      <variable units="microA_per_microF" name="I_pCa" initial_value="0.05"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
	  <variable units="millimolar" public_interface="in" name="MgATP_i"/>
	  <variable units="millimolar" public_interface="in" name="MgATP_i0"/>
	  
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_p_Ca</ci>
            <apply>
               <times/>
			   <apply>
                  <divide/>
				  <ci>MgATP_i</ci>
				  <ci>MgATP_i0</ci>
			   </apply>
			   <ci>i_p_Ca_winslow</ci>
			</apply>
		 </apply>
		 
		 <apply>
            <eq/>
            <ci>i_p_Ca_winslow</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>I_pCa</ci>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_mpCa</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="calcium_background_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Ca_b"/>
      <variable units="milliS_per_microF" name="g_Cab" initial_value="0.0003842"/>
      <variable units="millivolt" name="E_Ca"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_o"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_o</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_b</ci>
            <apply>
               <times/>
               <ci>g_Cab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_background_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Na_b"/>
      <variable units="milliS_per_microF" name="g_Nab" initial_value="0.0031"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na_b</ci>
            <apply>
               <times/>
               <ci>g_Nab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="L_type_Ca_current">
      <variable units="microA_per_microF" public_interface="out" name="i_Ca"/>
      <variable units="microA_per_microF" public_interface="out" name="i_Ca_K"/>
      
	  <variable units="cm_per_second" name="P_Ca" initial_value="3.594e-4"/>
      <variable units="cm_per_second" name="p_prime_k" initial_value="6.658e-7"/>
      <variable units="dimensionless" name="f_Ca"/>
	  <variable units="dimensionless" name="P_rom"/>
      <variable units="millimolar" name="k_MgATP_ss" initial_value="1.4"/>
      <variable units="per_second" name="alpha"/>
      <variable units="per_second" name="beta"/>
	  <variable units="per_second" name="gamma"/>
      <variable units="per_second" name="v_gamma"/>
      <variable units="per_second" name="alpha_prime"/>
      <variable units="per_second" name="beta_prime"/>
      <variable units="dimensionless" name="a" initial_value="2"/>
      <variable units="dimensionless" name="b" initial_value="2"/>
      <variable units="per_second" name="g" initial_value="2000"/>
      <variable units="per_second" name="f" initial_value="300"/>
      <variable units="per_second" name="v_omega"/>
	  <variable units="per_second" name="omega" initial_value="10"/>
	  <variable units="dimensionless" name="x" />
      <variable units="dimensionless" name="z" initial_value="0.1"/>
	  <variable units="dimensionless" name="v" initial_value="0.1"/>
      <variable units="dimensionless" name="w" initial_value="0.1"/>
      
	  <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="microF_per_cm2" public_interface="in" name="C_sc"/>
      <variable units="millimolar" public_interface="in" name="Ca_ss"/>
	  <variable units="millimolar" public_interface="in" name="MgATP_ss"/>
      <variable units="millimolar" public_interface="in" name="Ca_o"/>
      <variable units="millimolar" public_interface="in" name="K_o"/>
      <variable units="millimolar" public_interface="in" name="K_i"/>
      <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      
	  <variable units="dimensionless" private_interface="in" name="y"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca</ci>
            <apply>
               <times/>
               <ci>f_Ca</ci>
			   <apply>
                  <divide/>
                  <ci>P_Ca</ci>
				  <ci>C_sc</ci>
			   </apply>
			   <apply>
                  <divide/>
                  <apply>
                     <times/> 
					 <cn cellml:units="dimensionless">2</cn>
					 <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>
                  <divide/>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
						<cn cellml:units="millimolar">0.001</cn>
						<apply>
                           <exp/>
						   <apply>
                              <divide/>
							  <apply>
                                 <times/>
								 <cn cellml:units="dimensionless">2</cn>
					             <ci>V</ci> 	
							  </apply>
                              <apply>
                                 <times/> 
					             <ci>R</ci>
					             <ci>T</ci>
                              </apply>
			               </apply>	
						</apply>
					 </apply>
                     <apply>
                        <times/>
						<cn cellml:units="dimensionless">0.34</cn>
                        <ci>Ca_o</ci>
                     </apply>
				  </apply>
                  <apply>
                     <minus/>
					 <apply>
                        <exp/>
						<apply>
                           <divide/>
						   <apply>
                              <times/>
						      <cn cellml:units="dimensionless">2</cn>
					          <ci>V</ci> 	
						   </apply>
                           <apply>
                              <times/> 
					          <ci>R</ci>
					          <ci>T</ci>
                           </apply>
			            </apply>
					 </apply>
					 <cn cellml:units="dimensionless">1.0</cn>	
				  </apply>
			   </apply>	
			   <ci>P_rom</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>i_Ca_K</ci>
            <apply>
               <times/>
			   <ci>f_Ca</ci>
			   <apply>
                  <divide/>
                  <ci>p_prime_k</ci>
                  <ci>C_sc</ci>
			   </apply>	
			   <ci>P_rom</ci>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>v</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha</ci>
				  <apply>
                     <minus/>
					 <cn cellml:units="dimensionless">1</cn>
					 <ci>v</ci> 	
			      </apply>
               </apply>
			   <apply>
                  <times/>
                  <ci>beta</ci>
                  <ci>v</ci>
			   </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>w</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_prime</ci>
				  <apply>
                     <minus/>
					 <cn cellml:units="dimensionless">1</cn>
					 <ci>w</ci> 	
			      </apply>
               </apply>
			   <apply>
                  <times/>
                  <ci>beta_prime</ci>
                  <ci>w</ci>
			   </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>x</ci>
            <apply>
               <divide/>
			   <ci>f</ci>
               <apply>
                  <plus/>
                  <ci>f</ci>
                  <ci>g</ci>
			   </apply>	
			</apply>
         </apply>
		 
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>z</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>v_omega</ci>
				  <apply>
                     <minus/>
					 <cn cellml:units="dimensionless">1</cn>
					 <ci>z</ci> 	
			      </apply>
               </apply>
			   <apply>
                  <times/>
                  <ci>v_gamma</ci>
                  <ci>z</ci>
			   </apply>
            </apply>
         </apply>
		 
         <apply>
            <eq/>
            <ci>alpha</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">400</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">2</cn>
                     </apply>
                     <cn cellml:units="millivolt">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">50</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">2</cn>
                        </apply>
                     </apply>
                     <cn cellml:units="millivolt">13</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>alpha_prime</ci>
            <apply>
               <times/>
               <ci>alpha</ci>
               <ci>a</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>beta_prime</ci>
            <apply>
               <divide/>
               <ci>beta</ci>
               <ci>b</ci>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>gamma</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="per_second">103.75</cn>
                  <ci>Ca_ss</ci>
               </apply>
               <cn cellml:units="millimolar">1</cn>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>P_rom</ci>
            <apply>
               <times/>
			   <apply>
                  <power/>
                  <ci>v</ci>
                  <cn cellml:units="dimensionless">4</cn>
               </apply>
			   <ci>y</ci>
			   <ci>z</ci>
               <apply>
                  <divide/>
                  <ci>f</ci>
				  <apply>
                     <plus/>
                     <ci>f</ci>
                     <ci>g</ci>
			      </apply>	
			   </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>f_Ca</ci>
            <apply>
               <divide/>
			   <cn cellml:units="dimensionless">1</cn>
			   <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
				  <apply>
                     <power/>
					 <apply>
                        <divide/>
                        <ci>k_MgATP_ss</ci>
						<ci>MgATP_ss</ci>
			         </apply>
					 <cn cellml:units="dimensionless">2.6</cn>
                  </apply>	
			   </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>v_omega</ci>
            <apply>
               <times/>
			   <ci>omega</ci>
               <apply>
                  <plus/>
				  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>b</ci>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">0</cn>
			            </apply>
					 </apply>
                     <apply>
                        <power/>
                        <ci>w</ci>
                        <cn cellml:units="dimensionless">0</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>w</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">0</cn>
			            </apply>
					 </apply>
				  </apply>
                  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>b</ci>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
			            </apply>
					 </apply>
                     <apply>
                        <power/>
                        <ci>w</ci>
                        <cn cellml:units="dimensionless">1</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>w</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">1</cn>
			            </apply>
					 </apply>
				  </apply>
                  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>b</ci>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">2</cn>
			            </apply>
					 </apply>
                     <apply>
                        <power/>
                        <ci>w</ci>
                        <cn cellml:units="dimensionless">2</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>w</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">2</cn>
			            </apply>
					 </apply>
				  </apply>
                  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>b</ci>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">3</cn>
			            </apply>
					 </apply>
                     <apply>
                        <power/>
                        <ci>w</ci>
                        <cn cellml:units="dimensionless">3</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>w</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">3</cn>
			            </apply>
					 </apply>
				  </apply>	
                  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>b</ci>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
			            </apply>
					 </apply>
                     <apply>
                        <power/>
                        <ci>w</ci>
                        <cn cellml:units="dimensionless">4</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>w</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">4</cn>
			            </apply>
					 </apply>
				  </apply>						  
			   </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>v_gamma</ci>
            <apply>
               <times/>
			   <ci>gamma</ci>
               <apply>
                  <plus/>
				  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>a</ci>
                        <cn cellml:units="dimensionless">0</cn>
			         </apply>
                     <apply>
                        <power/>
                        <ci>v</ci>
                        <cn cellml:units="dimensionless">0</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>v</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">0</cn>
			            </apply>
					 </apply>
				  </apply>
                  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>a</ci>
                        <cn cellml:units="dimensionless">1</cn>
			         </apply>
                     <apply>
                        <power/>
                        <ci>v</ci>
                        <cn cellml:units="dimensionless">1</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>v</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">1</cn>
			            </apply>
					 </apply>
				  </apply>
                  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>a</ci>
                        <cn cellml:units="dimensionless">2</cn>
			         </apply>
                     <apply>
                        <power/>
                        <ci>v</ci>
                        <cn cellml:units="dimensionless">2</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>v</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">2</cn>
			            </apply>
					 </apply>
				  </apply>
                  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>a</ci>
                        <cn cellml:units="dimensionless">3</cn>
			         </apply>
                     <apply>
                        <power/>
                        <ci>v</ci>
                        <cn cellml:units="dimensionless">3</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>v</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">3</cn>
			            </apply>
					 </apply>
				  </apply>	
                  <apply>
                     <times/>
					 <apply>
                        <power/>
                        <ci>a</ci>
                        <cn cellml:units="dimensionless">4</cn>
			         </apply>
                     <apply>
                        <power/>
                        <ci>v</ci>
                        <cn cellml:units="dimensionless">4</cn>
			         </apply>
                     <apply>
                        <power/>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">1</cn>
						   <ci>v</ci>
						</apply>
                        <apply>
                           <minus/>
						   <cn cellml:units="dimensionless">4</cn>
						   <cn cellml:units="dimensionless">4</cn>
			            </apply>
					 </apply>
				  </apply>	
				  <apply>
                     <minus/>
					 <apply>
                        <times/>
						<apply>
                           <power/>
						   <ci>a</ci>
                           <cn cellml:units="dimensionless">4</cn>
			            </apply>
						<apply>
                           <power/>
						   <ci>v</ci>
                           <cn cellml:units="dimensionless">4</cn>
			            </apply>
						<ci>x</ci>
				     </apply>
				  </apply>
			   </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="L_type_Ca_current_y_gate">
      <variable units="dimensionless" public_interface="out" name="y" initial_value="0.798"/>
      <variable units="dimensionless" name="y_infinity"/>
      <variable units="second" name="tau_y"/>
      <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>y_infinity</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">0.8</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">12.5</cn>
                           </apply>
                           <cn cellml:units="millivolt">5</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="dimensionless">0.2</cn>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>tau_y</ci>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <cn cellml:units="second">20</cn>
                  <apply>
                     <divide/>
                     <cn cellml:units="second">600</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">9.5</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="dimensionless">1000</cn>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>y</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>y_infinity</ci>
                  <ci>y</ci>
               </apply>
               <ci>tau_y</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="ATP_sensitive_potassium_current">
	  <variable name="i_K_ATP" public_interface="out" units="microA_per_microF"/>
      
	  <variable name="g_K_ATP" units="milliS_per_microF"/>
	  <variable name="G_K_ATP" initial_value="0.05" units="milliS_per_microF"/>
	  <variable name="f_K_ATP" units="dimensionless"/>
	  <variable name="f_ATP" units="dimensionless"/>
	  <variable name="k_ATP" initial_value="0.6" units="millimolar"/>
	  <variable name="k_MgADP" initial_value="0.4" units="millimolar"/>
	  <variable name="K_o_normal" initial_value="4.5" units="millimolar"/>
	  <variable name="f_MgADP_" units="dimensionless"/>
	  <variable name="f_MgADP" units="dimensionless"/>
	  <variable name="c_MgADP" units="dimensionless"/>
	  <variable name="go" initial_value="0.08" units="dimensionless"/>
	  <variable name="gd" initial_value="0.89" units="dimensionless"/>
      
	  <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
	  <variable name="K_o" public_interface="in" units="millimolar"/>
	  <variable name="ATP_i" public_interface="in" units="millimolar"/>
	  <variable name="MgADP_i" public_interface="in" units="millimolar"/>
      <variable name="time" public_interface="in" units="second"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML" id="1">
         <apply>
            <eq/>
            <ci>i_K_ATP</ci>
            <apply>
               <times/>
               <ci>g_K_ATP</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
	  </math>
       
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="2">	   
		 <apply>
            <eq/>
            <ci>g_K_ATP</ci>
            <apply>
               <times/>
               <ci>G_K_ATP</ci>
			   <ci>f_K_ATP</ci>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>K_o</ci>
					 <ci>K_o_normal</ci>
                  </apply>
				  <cn cellml:units="dimensionless">0.24</cn>
               </apply>
            </apply>
         </apply>
      </math>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML" id="3">
	     <apply>
            <eq/>
            <ci>f_ATP</ci>
            <apply>
               <power/>
			   <apply>
                  <minus/>
				  <cn cellml:units="dimensionless">1.0</cn>
				  <apply>
                     <divide/>
                     <ci>ATP_i</ci>
					 <apply>
                        <plus/>
                        <ci>ATP_i</ci>
			            <ci>k_ATP</ci>
                     </apply>
				  </apply>
			   </apply>
               <cn cellml:units="dimensionless">4</cn>
            </apply>
         </apply>
      </math>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML" id="4">
	     <apply>
            <eq/>
            <ci>f_MgADP_</ci>
            <apply>
               <power/>
			   <apply>
                  <divide/>
                  <ci>MgADP_i</ci>
			      <apply>
                     <plus/>
                     <ci>MgADP_i</ci>
			         <ci>k_MgADP</ci>
                  </apply>
			   </apply>
               <cn cellml:units="dimensionless">2</cn>
            </apply>
         </apply>
      </math>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML" id="5">
	     <apply>
            <eq/>
            <ci>c_MgADP</ci>
            <apply>
               <minus/>
			   <cn cellml:units="dimensionless">1</cn>
			   <ci>f_MgADP_</ci>
			</apply>
         </apply>
      </math>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML" id="6">
	     <apply>
            <eq/>
            <ci>f_MgADP</ci>
            <apply>
               <minus/>
			   <cn cellml:units="dimensionless">1</cn>
			   <apply>
                  <power/>
			      <ci>c_MgADP</ci>
				  <cn cellml:units="dimensionless">4</cn>
               </apply>
			</apply>
         </apply>
      </math>
	  
	  <math xmlns="http://www.w3.org/1998/Math/MathML" id="7">
	     <apply>
            <eq/>
            <ci>f_K_ATP</ci>
            <apply>
               <plus/>
			   <apply>
                  <times/>
				  <ci>go</ci>
				  <ci>f_ATP</ci>
				  <apply>
                     <minus/>
			         <cn cellml:units="dimensionless">1</cn>
					 <ci>f_MgADP</ci>
			      </apply>
			   </apply>
			   <apply>
                  <times/>
				  <ci>gd</ci>
				  <ci>f_ATP</ci>
                  <ci>f_MgADP</ci>
			   </apply>
			</apply>
         </apply>
      </math>
   </component>
   
   <component name="RyR_channel">
      <variable units="millimolar_per_second" public_interface="out" name="J_rel"/>
      <variable units="per_second" name="v1" initial_value="1800"/>
      <variable units="millimolar4_per_second" name="k_a_plus" initial_value="1.215e13"/>
      <variable units="per_second" name="k_a_minus" initial_value="576"/>
      <variable units="millimolar3_per_second" name="k_b_plus" initial_value="4.05e9"/>
      <variable units="per_second" name="k_b_minus" initial_value="1930"/>
      <variable units="per_second" name="k_c_plus" initial_value="100"/>
      <variable units="per_second" name="k_c_minus" initial_value="0.8"/>
      <variable units="dimensionless" name="P_O1" initial_value="0"/>
      <variable units="dimensionless" name="P_O2" initial_value="0"/>
      <variable units="dimensionless" name="P_C1" initial_value="0.47"/>
      <variable units="dimensionless" name="P_C2" initial_value="0.53"/>
      <variable units="dimensionless" name="n" initial_value="4"/>
      <variable units="dimensionless" name="m" initial_value="3"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millimolar" public_interface="in" name="Ca_ss"/>
      <variable units="millimolar" public_interface="in" name="Ca_JSR"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_C1</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>k_a_plus</ci>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Ca_ss</ci>
                     <ci>n</ci>
                  </apply>
                  <ci>P_C1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_a_minus</ci>
                  <ci>P_O1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_O1</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>k_a_plus</ci>
                     <apply>
                        <power/>
                        <ci>Ca_ss</ci>
                        <ci>n</ci>
                     </apply>
                     <ci>P_C1</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>k_a_minus</ci>
                        <ci>P_O1</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>k_b_plus</ci>
                        <apply>
                           <power/>
                           <ci>Ca_ss</ci>
                           <ci>m</ci>
                        </apply>
                        <ci>P_O1</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>k_c_plus</ci>
                        <ci>P_O1</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_b_minus</ci>
                  <ci>P_O2</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_c_minus</ci>
                  <ci>P_C2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_O2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_b_plus</ci>
                  <apply>
                     <power/>
                     <ci>Ca_ss</ci>
                     <ci>m</ci>
                  </apply>
                  <ci>P_O1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_b_minus</ci>
                  <ci>P_O2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_C2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_c_plus</ci>
                  <ci>P_O1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_c_minus</ci>
                  <ci>P_C2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_rel</ci>
            <apply>
               <times/>
               <ci>v1</ci>
               <apply>
                  <plus/>
                  <ci>P_O1</ci>
                  <ci>P_O2</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>Ca_JSR</ci>
                  <ci>Ca_ss</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="SERCA2a_pump">
      <variable units="millimolar_per_second" public_interface="out" name="J_up"/>
	  
	  <variable units="millimolar_per_second" name="J_up_winslow"/>
      <variable units="millimolar" name="K_fb" initial_value="0.000168"/>
      <variable units="millimolar" name="K_rb" initial_value="3.29"/>
      <variable units="dimensionless" name="fb"/>
      <variable units="dimensionless" name="rb"/>
      <variable units="millimolar_per_second" name="Vmaxf" initial_value="0.0813"/>
      <variable units="millimolar_per_second" name="Vmaxr" initial_value="0.318"/>
      <variable units="dimensionless" name="K_SR" initial_value="1"/>
      <variable units="dimensionless" name="N_fb" initial_value="1.2"/>
      <variable units="dimensionless" name="N_rb" initial_value="1"/>
      <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_NSR"/>
	  <variable units="millimolar" public_interface="in" name="MgATP_i"/>
	  <variable units="millimolar" public_interface="in" name="MgATP_i0"/>
	  
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>fb</ci>
            <apply>
               <power/>
               <apply>
                  <divide/>
                  <ci>Ca_i</ci>
                  <ci>K_fb</ci>
               </apply>
               <ci>N_fb</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>rb</ci>
            <apply>
               <power/>
               <apply>
                  <divide/>
                  <ci>Ca_NSR</ci>
                  <ci>K_rb</ci>
               </apply>
               <ci>N_rb</ci>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>J_up</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <ci>MgATP_i</ci>
				  <ci>MgATP_i0</ci>
			   </apply>
               <ci>J_up_winslow</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>J_up_winslow</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>K_SR</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>Vmaxf</ci>
                        <ci>fb</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>Vmaxr</ci>
                        <ci>rb</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>fb</ci>
                  <ci>rb</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="intracellular_Ca_fluxes">
      <variable units="millimolar_per_second" public_interface="out" name="J_tr"/>
      <variable units="millimolar_per_second" public_interface="out" name="J_xfer"/>
      <variable units="millimolar_per_second" public_interface="out" name="J_trpn"/>
	  
      <variable units="second" name="tau_tr" initial_value="0.0005747"/>
      <variable units="second" name="tau_xfer" initial_value="0.0267"/>
      <variable units="millimolar" name="HTRPNCa" initial_value="0.98"/>
      <variable units="millimolar" name="LTRPNCa" initial_value="0.078"/>
      <variable units="millimolar_per_second" name="J_HTRPNCa"/>
      <variable units="millimolar_per_second" name="J_LTRPNCa"/>
      <variable units="dimensionless" name="HTRPN_tot" initial_value="0.14"/>
      <variable units="dimensionless" name="LTRPN_tot" initial_value="0.07"/>
      <variable units="per_millimolar_second" name="k_htrpn_plus" initial_value="20000"/>
      <variable units="per_second" name="k_htrpn_minus" initial_value="0.066"/>
      <variable units="per_millimolar_second" name="k_ltrpn_plus" initial_value="40000"/>
      <variable units="per_second" name="k_ltrpn_minus" initial_value="40"/>
      
	  <variable units="second" public_interface="in" name="time"/>
      <variable units="millimolar" public_interface="in" name="Ca_ss"/>
      <variable units="millimolar" public_interface="in" name="Ca_i"/>
      <variable units="millimolar" public_interface="in" name="Ca_NSR"/>
      <variable units="millimolar" public_interface="in" name="Ca_JSR"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_tr</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Ca_NSR</ci>
                  <ci>Ca_JSR</ci>
               </apply>
               <ci>tau_tr</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_xfer</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Ca_ss</ci>
                  <ci>Ca_i</ci>
               </apply>
               <ci>tau_xfer</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_trpn</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>HTRPN_tot</ci>
                  <ci>J_HTRPNCa</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>LTRPN_tot</ci>
                  <ci>J_LTRPNCa</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_HTRPNCa</ci>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>HTRPNCa</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>HTRPNCa</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_htrpn_plus</ci>
                  <ci>Ca_i</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="millimolar">1</cn>
                     <ci>HTRPNCa</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_htrpn_minus</ci>
                  <ci>HTRPNCa</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_LTRPNCa</ci>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>LTRPNCa</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>LTRPNCa</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_ltrpn_plus</ci>
                  <ci>Ca_i</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="millimolar">1</cn>
                     <ci>LTRPNCa</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_ltrpn_minus</ci>
                  <ci>LTRPNCa</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="intracellular_ion_concentrations">
      <variable units="millimolar" public_interface="out" name="Na_i" initial_value="10"/>
      <variable units="millimolar" public_interface="out" name="Ca_i" initial_value="8.464E-5"/>
      <variable units="millimolar" public_interface="out" name="K_i" initial_value="159.48"/>
      <variable units="millimolar" public_interface="out" name="Ca_ss" initial_value="1.315E-4"/>
      <variable units="millimolar" public_interface="out" name="Ca_JSR" initial_value="0.2616"/>
      <variable units="millimolar" public_interface="out" name="Ca_NSR" initial_value="0.2620"/>
      
	  <variable units="cm2" name="A_cap" initial_value="0.0001534"/>
      <variable units="microlitre" name="V_myo" initial_value="0.00002584"/>
      <variable units="microlitre" name="V_JSR" initial_value="0.00000016"/>
      <variable units="microlitre" name="V_NSR" initial_value="0.0000021"/>
      <variable units="microlitre" name="V_ss" initial_value="0.0000000012"/>
      <variable units="millimolar" name="K_mCMDN" initial_value="0.00238"/>
      <variable units="millimolar" name="K_mEGTA" initial_value="0.00015"/>
      <variable units="millimolar" name="K_mCSQN" initial_value="0.8"/>
      <variable units="millimolar" name="CMDN_tot" initial_value="0.05"/>
      <variable units="millimolar" name="EGTA_tot" initial_value="0"/>
      <variable units="millimolar" name="CSQN_tot" initial_value="15"/>
      <variable units="dimensionless" name="beta_i"/>
      <variable units="dimensionless" name="beta_SS"/>
      <variable units="dimensionless" name="beta_JSR"/>
	  
      <variable units="second" public_interface="in" name="time"/>
      <variable units="coulomb_per_millimole" public_interface="in" name="F"/>
      <variable units="microF_per_cm2" public_interface="in" name="C_sc"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Na"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Na_b"/>
      <variable units="microA_per_microF" public_interface="in" name="i_NaCa"/>
      <variable units="microA_per_microF" public_interface="in" name="i_NaK"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca_K"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Kr"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ks"/>
      <variable units="microA_per_microF" public_interface="in" name="i_K1"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Kp"/>
      <variable units="microA_per_microF" public_interface="in" name="i_to1"/>
      <variable units="microA_per_microF" public_interface="in" name="i_p_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="i_Ca_b"/>
	  <variable units="microA_per_microF" public_interface="in" name="i_K_ATP"/>
	  <variable units="millimolar_per_second" public_interface="in" name="J_up"/>
      <variable units="millimolar_per_second" public_interface="in" name="J_rel"/>
      <variable units="millimolar_per_second" public_interface="in" name="J_xfer"/>
      <variable units="millimolar_per_second" public_interface="in" name="J_trpn"/>
      <variable units="millimolar_per_second" public_interface="in" name="J_tr"/>
	  <variable units="per_millimolar_second" public_interface="in" name="k_plus_CaATP"/>
      <variable units="per_second" public_interface="in" name="k_minus_CaATP"/>
      <variable units="per_millimolar_second" public_interface="in" name="k_plus_CaADP"/>
      <variable units="per_second" public_interface="in" name="k_minus_CaADP"/>
      <variable units="millimolar" public_interface="in" name="CaADP_i"/>
      <variable units="millimolar" public_interface="in" name="CaADP_ss"/>
      <variable units="millimolar" public_interface="in" name="CaATP_ss"/>
      <variable units="millimolar" public_interface="in" name="ATP_i"/>
      <variable units="millimolar" public_interface="in" name="ADP_i"/>
      <variable units="millimolar" public_interface="in" name="ADP_ss"/>
      <variable units="millimolar" public_interface="in" name="ATP_ss"/>
      <variable units="millimolar" public_interface="in" name="CaATP_i"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_i</ci>
            </apply>
            <apply>
               <times/>
               <ci>beta_i</ci>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <ci>J_xfer</ci>
                     <apply>
                        <plus/>
                        <ci>J_up</ci>
                        <ci>J_trpn</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
					 <apply>
                        <times/>
                        <apply>
                           <minus/>
						   <apply>
                              <plus/>
                              <ci>i_p_Ca</ci>
                              <ci>i_Ca_b</ci>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">2</cn>
                              <ci>i_NaCa</ci>
                           </apply>
						</apply>
						<apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>A_cap</ci>
                              <ci>C_sc</ci>
                           </apply>
                           <apply>
                               <times/>
                               <cn cellml:units="dimensionless">2</cn>
                               <ci>V_myo</ci>
                               <ci>F</ci>
						   </apply>
						</apply>
                     </apply>
                  </apply>
				  <apply>
                     <minus/>
					 <apply>
                        <plus/>
						<apply>
                           <times/>
                           <ci> k_minus_CaATP </ci>
                           <ci> CaATP_i </ci>
                        </apply>
						<apply>
                           <times/>
                           <ci> k_minus_CaADP </ci>
                           <ci> CaADP_i </ci>
                        </apply>
					 </apply>
					 <apply>
                        <plus/>
						<apply>
                           <times/>
                           <ci> k_plus_CaATP </ci>
                           <ci> Ca_i </ci>
						   <ci> ATP_i </ci>
                        </apply>
						<apply>
                           <times/>
                           <ci> k_plus_CaADP </ci>
                           <ci> Ca_i </ci>
						   <ci> ADP_i </ci>
                        </apply>
					 </apply>
				  </apply>
               </apply>
            </apply>
         </apply>
	  
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Na_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>i_Na</ci>
                     <ci>i_Na_b</ci>
                     <apply>
                        <times/>
                        <ci>i_NaCa</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                     <apply>
                        <times/>
                        <ci>i_NaK</ci>
                        <cn cellml:units="dimensionless">3</cn>
                     </apply>
                  </apply>
                  <ci>A_cap</ci>
                  <cn cellml:units="microF_per_cm2">1</cn>
               </apply>
               <apply>
                  <times/>
                  <ci>V_myo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>K_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>i_Ca_K</ci>
                     <ci>i_Kr</ci>
                     <ci>i_Ks</ci>
                     <ci>i_K1</ci>
                     <ci>i_Kp</ci>
                     <ci>i_to1</ci>
					 <ci>i_K_ATP</ci>
                     <apply>
                        <times/>
                        <ci>i_NaK</ci>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <ci>A_cap</ci>
                  <cn cellml:units="microF_per_cm2">1</cn>
               </apply>
               <apply>
                  <times/>
                  <ci>V_myo</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_i</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>CMDN_tot</ci>
                        <ci>K_mCMDN</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>K_mCMDN</ci>
                           <ci>Ca_i</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>EGTA_tot</ci>
                        <ci>K_mEGTA</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>K_mEGTA</ci>
                           <ci>Ca_i</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_SS</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>CMDN_tot</ci>
                        <ci>K_mCMDN</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>K_mCMDN</ci>
                           <ci>Ca_ss</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>EGTA_tot</ci>
                        <ci>K_mEGTA</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>K_mEGTA</ci>
                           <ci>Ca_ss</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_JSR</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>CSQN_tot</ci>
                        <ci>K_mCSQN</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>K_mCSQN</ci>
                           <ci>Ca_JSR</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      
	     <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_ss</ci>
            </apply>
            <apply>
               <times/>
               <ci>beta_SS</ci>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>J_rel</ci>
                           <ci>V_JSR</ci>
                        </apply>
                        <ci>V_ss</ci>
                     </apply>
					 <apply>
                        <times/>
                        <ci>k_minus_CaATP</ci>
						<ci>CaATP_ss</ci>
				     </apply>
					 <apply>
                        <times/>
                        <ci>k_minus_CaADP</ci>
						<ci>CaADP_ss</ci>
				     </apply>
                  </apply>
				  <apply>
                     <plus/> 
					 <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>J_xfer</ci>
                           <ci>V_myo</ci>
                        </apply>
                        <ci>V_ss</ci>
                     </apply>
				     <apply>
                        <times/>
                        <ci>i_Ca</ci>
						<apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>A_cap</ci>
                              <cn cellml:units="microF_per_cm2">1</cn>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">2</cn>
                              <ci>V_ss</ci>
                              <ci>F</ci>
                           </apply>
                        </apply>
                     </apply>
				     <apply>
                        <times/>
                        <ci>k_plus_CaATP</ci>
						<ci>Ca_ss</ci>
					    <ci>ATP_ss</ci>
                     </apply>
				     <apply>
                        <times/>
                        <ci>k_plus_CaADP</ci>
						<ci>Ca_ss</ci>
						<ci>ADP_ss</ci>
                     </apply>
				  </apply>
               </apply>
            </apply>
         </apply>
	  
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_JSR</ci>
            </apply>
            <apply>
               <times/>
               <ci>beta_JSR</ci>
               <apply>
                  <minus/>
                  <ci>J_tr</ci>
                  <ci>J_rel</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_NSR</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>J_up</ci>
                     <ci>V_myo</ci>
                  </apply>
                  <ci>V_NSR</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>J_tr</ci>
                     <ci>V_JSR</ci>
                  </apply>
                  <ci>V_NSR</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="Ca_and_Mg_buffering_by_ATP">
    <variable units="millimolar" public_interface="out" name="CaATP_ss" initial_value="0.25E-3"/>
    <variable units="millimolar" public_interface="out" name="ATP_ss"/>
    <variable units="millimolar" public_interface="out" name="MgATP_i" initial_value="6.4395"/>
	<variable units="millimolar" public_interface="out" name="MgADP_i" initial_value="0.298E-2"/>
	<variable units="millimolar" public_interface="out" name="MgATP_i0" initial_value="2.888"/>
    <variable units="per_millimolar_second" public_interface="out" name="k_plus_CaATP" initial_value="225000.0"/>
    <variable units="per_second" public_interface="out" name="k_minus_CaATP" initial_value="45000.0"/>
    <variable units="per_millimolar_second" public_interface="out" name="k_plus_CaADP" initial_value="125000.0"/>
    <variable units="per_second" public_interface="out" name="k_minus_CaADP" initial_value="193500"/>
    <variable units="millimolar" public_interface="out" name="CaADP_i" initial_value="0.11E-6"/>
    <variable units="millimolar" public_interface="out" name="CaADP_ss" initial_value="0.13E-6"/>
    <variable units="millimolar" public_interface="out" name="ATP_i"/>
    <variable units="millimolar" public_interface="out" name="ADP_i"/>
    <variable units="millimolar" public_interface="out" name="ADP_ss"/>
    <variable units="millimolar" public_interface="out" name="CaATP_i" initial_value="0.237E-3"/>
    <variable units="millimolar" public_interface="out" name="MgATP_ss" initial_value="6.4395"/>
	
    <variable units="millimolar" name="Mg_ss" initial_value="1.0"/>
    <variable units="millimolar" name="Mg_i" initial_value="1.0"/>
    <variable units="millimolar" name="MgADP_ss" initial_value="0.298E-2"/>
    <variable units="millimolar" name="ATP_tot" initial_value="7.0"/>
    <variable units="per_millimolar_second" name="k_plus_MgATP" initial_value="125000.0"/>
    <variable units="per_second" name="k_minus_MgATP" initial_value="10875.0"/>
    <variable units="millimolar_per_second" name="Jxfer_CaATP"/>
    <variable units="millimolar_per_second" name="Jxfer_MgATP"/>
    <variable units="millimolar_per_second" name="Jxfer_Mg"/>
    <variable units="second" name="tau_xfer_CaATP" initial_value="0.0534"/>
    <variable units="second" name="tau_xfer_MgATP" initial_value="0.0534"/>
    <variable units="second" name="tau_xfer_Mg" initial_value="0.0267"/>
    <variable units="millimolar" name="ADP_tot" initial_value="0.005"/>
    <variable units="per_millimolar_second" name="k_plus_MgADP" initial_value="125000.0"/>
    <variable units="per_second" name="k_minus_MgADP" initial_value="84500.0"/>
    <variable units="millimolar_per_second" name="Jxfer_CaADP"/>
    <variable units="millimolar_per_second" name="Jxfer_MgADP"/>
    <variable units="second" name="tau_xfer_CaADP" initial_value="0.0534"/>
    <variable units="second" name="tau_xfer_MgADP" initial_value="0.0534"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="millimolar" public_interface="in" name="Ca_ss"/>
    <variable units="millimolar" public_interface="in" name="Ca_i"/>
    <variable units="microlitre" public_interface="in" name="V_myo"/>
    <variable units="microlitre" public_interface="in" name="V_ss"/>
  
    <math xmlns="http://www.w3.org/1998/Math/MathML">    
      <apply id="ATP_ss_calculation">
        <eq/>
        <ci> ATP_ss </ci>
        <apply>
          <minus/>
          <ci> ATP_tot </ci>
          <apply>
            <plus/>
            <ci> CaATP_ss </ci>
            <ci> MgATP_ss </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="CaATP_ss_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> CaATP_ss </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_plus_CaATP </ci>
            <ci> Ca_ss </ci>
            <ci> ATP_ss </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> Jxfer_CaATP </ci>
              <apply>
                <divide/>
                <ci> V_myo </ci>
                <ci> V_ss </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> k_minus_CaATP </ci>
              <ci> CaATP_ss </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="MgATP_ss_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> MgATP_ss </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_plus_MgATP </ci>
            <ci> Mg_ss </ci>
            <ci> ATP_ss </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> Jxfer_MgATP </ci>
              <apply>
                <divide/>
                <ci> V_myo </ci>
                <ci> V_ss </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> k_minus_MgATP </ci>
              <ci> MgATP_ss </ci>
            </apply>
          </apply>
        </apply>
      </apply>
     
      <apply id="ATP_i_calculation">
        <eq/>
        <ci> ATP_i </ci>
        <apply>
          <minus/>
          <ci> ATP_tot </ci>
          <apply>
            <plus/>
            <ci> CaATP_i </ci>
            <ci> MgATP_i </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="CaATP_i_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> CaATP_i </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> Jxfer_CaATP </ci>
            <apply>
              <times/>
              <ci> k_plus_CaATP </ci>
              <ci> Ca_i </ci>
              <ci> ATP_i </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_minus_CaATP </ci>
            <ci> CaATP_i </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="MgATP_i_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> MgATP_i </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> Jxfer_MgATP </ci>
            <apply>
              <times/>
              <ci> k_plus_MgATP </ci>
              <ci> Mg_i </ci>
              <ci> ATP_i </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_minus_MgATP </ci>
            <ci> MgATP_i </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Jxfer_CaATP_calculation">
        <eq/>
        <ci> Jxfer_CaATP </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> CaATP_ss </ci>
            <ci> CaATP_i </ci>
          </apply>
          <ci> tau_xfer_CaATP </ci>
        </apply>
      </apply>
      
      <apply id="Jxfer_MgATP_calculation">
        <eq/>
        <ci> Jxfer_MgATP </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> MgATP_ss </ci>
            <ci> MgATP_i </ci>
          </apply>
          <ci> tau_xfer_MgATP </ci>
        </apply>
      </apply>
      
      <apply id="ADP_ss_calculation">
        <eq/>
        <ci> ADP_ss </ci>
        <apply>
          <minus/>
          <ci> ADP_tot </ci>
          <apply>
            <plus/>
            <ci> CaADP_ss </ci>
            <ci> MgADP_ss </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="CaADP_ss_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> CaADP_ss </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_plus_CaADP </ci>
            <ci> Ca_ss </ci>
            <ci> ADP_ss </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> Jxfer_CaADP </ci>
              <apply>
                <divide/>
                <ci> V_myo </ci>
                <ci> V_ss </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> k_minus_CaADP </ci>
              <ci> CaADP_ss </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="MgADP_ss_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> MgADP_ss </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_plus_MgADP </ci>
            <ci> Mg_ss </ci>
            <ci> ADP_ss </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> Jxfer_MgADP </ci>
              <apply>
                <divide/>
                <ci> V_myo </ci>
                <ci> V_ss </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> k_minus_MgADP </ci>
              <ci> MgADP_ss </ci>
            </apply>
          </apply>
        </apply>
      </apply>
     
      <apply id="ADP_i_calculation">
        <eq/>
        <ci> ADP_i </ci>
        <apply>
          <minus/>
          <ci> ADP_tot </ci>
          <apply>
            <plus/>
            <ci> CaADP_i </ci>
            <ci> MgADP_i </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="CaADP_i_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> CaADP_i </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> Jxfer_CaADP </ci>
            <apply>
              <times/>
              <ci> k_plus_CaADP </ci>
              <ci> Ca_i </ci>
              <ci> ADP_i </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_minus_CaADP </ci>
            <ci> CaADP_i </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="MgADP_i_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> MgADP_i </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> Jxfer_MgADP </ci>
            <apply>
              <times/>
              <ci> k_plus_MgADP </ci>
              <ci> Mg_i </ci>
              <ci> ADP_i </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_minus_MgADP </ci>
            <ci> MgADP_i </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Jxfer_CaADP_calculation">
        <eq/>
        <ci> Jxfer_CaADP </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> CaADP_ss </ci>
            <ci> CaADP_i </ci>
          </apply>
          <ci> tau_xfer_CaADP </ci>
        </apply>
      </apply>
      
      <apply id="Jxfer_MgADP_calculation">
        <eq/>
        <ci> Jxfer_MgADP </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> MgADP_ss </ci>
            <ci> MgADP_i </ci>
          </apply>
          <ci> tau_xfer_MgADP </ci>
        </apply>
      </apply>
      
      <apply id="Mg_ss_calculation">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Mg_ss </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> k_minus_MgATP </ci>
              <ci> MgATP_ss </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_minus_MgADP </ci>
              <ci> MgADP_ss </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> k_plus_MgATP </ci>
              <ci> Mg_ss </ci>
              <ci> ATP_ss </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_plus_MgADP </ci>
              <ci> Mg_ss </ci>
              <ci> ADP_ss </ci>
            </apply>
            <apply>
              <times/>
              <ci> Jxfer_Mg </ci>
              <apply>
                <divide/>
                <ci> V_myo </ci>
                <ci> V_ss </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Mg_i_calculation">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Mg_i </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> Jxfer_Mg </ci>
            <apply>
              <times/>
              <ci> k_minus_MgATP </ci>
              <ci> MgATP_i </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_minus_MgADP </ci>
              <ci> MgADP_i </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> k_plus_MgATP </ci>
              <ci> Mg_i </ci>
              <ci> ATP_i </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_plus_MgADP </ci>
              <ci> Mg_i </ci>
              <ci> ADP_i </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Jxfer_Mg_calculation">
        <eq/>
        <ci> Jxfer_Mg </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> Mg_ss </ci>
            <ci> Mg_i </ci>
          </apply>
          <ci> tau_xfer_Mg </ci>
        </apply>
      </apply>
    </math>
  </component>
   
   <component name="extracellular_ion_concentrations">
      <variable units="millimolar" public_interface="out" name="Na_o" initial_value="138"/>
      <variable units="millimolar" public_interface="out" name="Ca_o" initial_value="2"/>
      <variable units="millimolar" public_interface="out" name="K_o" initial_value="4"/>
   </component>
   
   <component name="model_parameters">
    <variable units="microlitre" public_interface="out" name="V_myo" initial_value="0.00002584"/>
    <variable units="microlitre" public_interface="out" name="V_ss" initial_value="0.0000000012"/>
  </component>
   
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_h_gate"/>
         <component_ref component="fast_sodium_current_j_gate"/>
      </component_ref>
      <component_ref component="rapid_activating_delayed_rectifiyer_K_current">
         <component_ref component="rapid_activating_delayed_rectifiyer_K_current_X_kr_gate"/>
      </component_ref>
      <component_ref component="slow_activating_delayed_rectifiyer_K_current">
         <component_ref component="slow_activating_delayed_rectifiyer_K_current_X_ks_gate"/>
      </component_ref>
      <component_ref component="transient_outward_potassium_current">
         <component_ref component="transient_outward_potassium_current_X_to1_gate"/>
         <component_ref component="transient_outward_potassium_current_Y_to1_gate"/>
      </component_ref>
      <component_ref component="time_independent_potassium_current">
         <component_ref component="time_independent_potassium_current_K1_gate"/>
      </component_ref>
      <component_ref component="plateau_potassium_current">
         <component_ref component="plateau_potassium_current_Kp_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_current">
         <component_ref component="L_type_Ca_current_y_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="fast_sodium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="rapid_activating_delayed_rectifiyer_K_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="slow_activating_delayed_rectifiyer_K_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="transient_outward_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="plateau_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="Na_Ca_exchanger"/>
      <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"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_background_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="sodium_background_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="L_type_Ca_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="ATP_sensitive_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="RyR_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="SERCA2a_pump"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="intracellular_Ca_fluxes"/>
      <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="Ca_and_Mg_buffering_by_ATP"/>
    <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="ATP_sensitive_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K_ATP" variable_1="i_K_ATP"/>
   </connection>
   <connection>
      <map_components component_2="fast_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="rapid_activating_delayed_rectifiyer_K_current" 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="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="slow_activating_delayed_rectifiyer_K_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="transient_outward_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_to1" variable_1="i_to1"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" 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="plateau_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Kp" variable_1="i_Kp"/>
   </connection>
   <connection>
      <map_components component_2="Na_Ca_exchanger" 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="sodium_potassium_pump" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
      <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="sarcolemmal_calcium_pump" component_1="membrane"/>
      <map_variables variable_2="i_p_Ca" variable_1="i_p_Ca"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ca_b" variable_1="i_Ca_b"/>
      <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_background_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Na_b" variable_1="i_Na_b"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Ca" variable_1="i_Ca"/>
      <map_variables variable_2="i_Ca_K" variable_1="i_Ca_K"/>
      <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="C_sc" variable_1="C_sc"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="membrane"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="C_sc" variable_1="C_sc"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="fast_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="sodium_background_current" component_1="fast_sodium_current"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_ion_concentrations" component_1="fast_sodium_current"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
   </connection>
   <connection>
      <map_components component_2="rapid_activating_delayed_rectifiyer_K_current" component_1="intracellular_ion_concentrations"/>
      <map_variables variable_2="i_Kr" variable_1="i_Kr"/>
      <map_variables variable_2="K_i" variable_1="K_i"/>
   </connection>
   <connection>
      <map_components component_2="rapid_activating_delayed_rectifiyer_K_current" component_1="extracellular_ion_concentrations"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="rapid_activating_delayed_rectifiyer_K_current" component_1="transient_outward_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="rapid_activating_delayed_rectifiyer_K_current" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="rapid_activating_delayed_rectifiyer_K_current" component_1="plateau_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="rapid_activating_delayed_rectifiyer_K_current" component_1="ATP_sensitive_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="ATP_sensitive_potassium_current" component_1="extracellular_ion_concentrations"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="slow_activating_delayed_rectifiyer_K_current" component_1="intracellular_ion_concentrations"/>
      <map_variables variable_2="i_Ks" variable_1="i_Ks"/>
      <map_variables variable_2="K_i" variable_1="K_i"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
   </connection>
   <connection>
      <map_components component_2="slow_activating_delayed_rectifiyer_K_current" component_1="extracellular_ion_concentrations"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="transient_outward_potassium_current"/>
      <map_variables variable_2="i_to1" variable_1="i_to1"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="i_K1" variable_1="i_K1"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_ion_concentrations" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="plateau_potassium_current"/>
      <map_variables variable_2="i_Kp" variable_1="i_Kp"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="Na_Ca_exchanger"/>
      <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="extracellular_ion_concentrations" component_1="Na_Ca_exchanger"/>
      <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="Na_i" variable_1="Na_i"/>
      <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_ion_concentrations" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="Na_o" variable_1="Na_o"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sarcolemmal_calcium_pump"/>
      <map_variables variable_2="i_p_Ca" variable_1="i_p_Ca"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="calcium_background_current"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="i_Ca_b" variable_1="i_Ca_b"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_ion_concentrations" component_1="calcium_background_current"/>
      <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="sodium_background_current"/>
      <map_variables variable_2="i_Na_b" variable_1="i_Na_b"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="L_type_Ca_current"/>
      <map_variables variable_2="i_Ca_K" variable_1="i_Ca_K"/>
      <map_variables variable_2="i_Ca" variable_1="i_Ca"/>
      <map_variables variable_2="K_i" variable_1="K_i"/>
      <map_variables variable_2="Ca_ss" variable_1="Ca_ss"/>
   </connection>
   <connection>
      <map_components component_2="extracellular_ion_concentrations" component_1="L_type_Ca_current"/>
      <map_variables variable_2="K_o" variable_1="K_o"/>
      <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="RyR_channel"/>
      <map_variables variable_2="J_rel" variable_1="J_rel"/>
      <map_variables variable_2="Ca_ss" variable_1="Ca_ss"/>
      <map_variables variable_2="Ca_JSR" variable_1="Ca_JSR"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="SERCA2a_pump"/>
      <map_variables variable_2="J_up" variable_1="J_up"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Ca_NSR" variable_1="Ca_NSR"/>
   </connection>
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="intracellular_Ca_fluxes"/>
      <map_variables variable_2="J_tr" variable_1="J_tr"/>
      <map_variables variable_2="J_xfer" variable_1="J_xfer"/>
      <map_variables variable_2="J_trpn" variable_1="J_trpn"/>
      <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
      <map_variables variable_2="Ca_NSR" variable_1="Ca_NSR"/>
      <map_variables variable_2="Ca_ss" variable_1="Ca_ss"/>
      <map_variables variable_2="Ca_JSR" variable_1="Ca_JSR"/>
   </connection>
   
   <connection>
    <map_components component_2="ATP_sensitive_potassium_current" component_1="Ca_and_Mg_buffering_by_ATP"/>
    <map_variables variable_2="ATP_i" variable_1="ATP_i"/>
	<map_variables variable_2="MgADP_i" variable_1="MgADP_i"/>
   </connection>
   
   <connection>
    <map_components component_2="L_type_Ca_current" component_1="Ca_and_Mg_buffering_by_ATP"/>
    <map_variables variable_2="MgATP_ss" variable_1="MgATP_ss"/>
   </connection>
  
   <connection>
      <map_components component_2="intracellular_ion_concentrations" component_1="ATP_sensitive_potassium_current"/>
      <map_variables variable_2="i_K_ATP" variable_1="i_K_ATP"/>
   </connection>
   
  <connection>
    <map_components component_2="sarcolemmal_calcium_pump" component_1="Ca_and_Mg_buffering_by_ATP"/>
    <map_variables variable_2="MgATP_i" variable_1="MgATP_i"/>
	<map_variables variable_2="MgATP_i0" variable_1="MgATP_i0"/>
  </connection>
  
  <connection>
    <map_components component_2="SERCA2a_pump" component_1="Ca_and_Mg_buffering_by_ATP"/>
    <map_variables variable_2="MgATP_i" variable_1="MgATP_i"/>
	<map_variables variable_2="MgATP_i0" variable_1="MgATP_i0"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_and_Mg_buffering_by_ATP" component_1="model_parameters"/>
    <map_variables variable_2="V_myo" variable_1="V_myo"/>
    <map_variables variable_2="V_ss" variable_1="V_ss"/>
  </connection>
   
  <connection>
    <map_components component_2="intracellular_ion_concentrations" component_1="Ca_and_Mg_buffering_by_ATP"/>
    <map_variables variable_2="CaATP_ss" variable_1="CaATP_ss"/>
    <map_variables variable_2="CaATP_i" variable_1="CaATP_i"/>
    <map_variables variable_2="CaADP_ss" variable_1="CaADP_ss"/>
    <map_variables variable_2="CaADP_i" variable_1="CaADP_i"/>
    <map_variables variable_2="Ca_ss" variable_1="Ca_ss"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
    <map_variables variable_2="ATP_ss" variable_1="ATP_ss"/>
    <map_variables variable_2="ATP_i" variable_1="ATP_i"/>
    <map_variables variable_2="ADP_ss" variable_1="ADP_ss"/>
    <map_variables variable_2="ADP_i" variable_1="ADP_i"/>
    <map_variables variable_2="k_plus_CaATP" variable_1="k_plus_CaATP"/>
    <map_variables variable_2="k_minus_CaATP" variable_1="k_minus_CaATP"/>
    <map_variables variable_2="k_plus_CaADP" variable_1="k_plus_CaADP"/>
    <map_variables variable_2="k_minus_CaADP" variable_1="k_minus_CaADP"/>
  </connection>
  
   <connection>
      <map_components component_2="fast_sodium_current_m_gate" component_1="fast_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="fast_sodium_current_h_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="h" variable_1="h"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_j_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="j" variable_1="j"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="rapid_activating_delayed_rectifiyer_K_current_X_kr_gate" component_1="rapid_activating_delayed_rectifiyer_K_current"/>
      <map_variables variable_2="X_kr" variable_1="X_kr"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="slow_activating_delayed_rectifiyer_K_current_X_ks_gate" component_1="slow_activating_delayed_rectifiyer_K_current"/>
      <map_variables variable_2="X_ks" variable_1="X_ks"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="transient_outward_potassium_current_X_to1_gate" component_1="transient_outward_potassium_current"/>
      <map_variables variable_2="X_to1" variable_1="X_to1"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="transient_outward_potassium_current_Y_to1_gate" component_1="transient_outward_potassium_current"/>
      <map_variables variable_2="Y_to1" variable_1="Y_to1"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current_K1_gate" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
      <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="E_K" variable_1="E_K"/>
      <map_variables variable_2="K1_infinity_V" variable_1="K1_infinity_V"/>
   </connection>
   <connection>
      <map_components component_2="plateau_potassium_current_Kp_gate" component_1="plateau_potassium_current"/>
      <map_variables variable_2="Kp_V" variable_1="Kp_V"/>
      <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_current_y_gate" component_1="L_type_Ca_current"/>
      <map_variables variable_2="y" variable_1="y"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   
   <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
    <rdf:Description rdf:about="">
      
      <dc:creator rdf:parseType="Resource">
        <vCard:N rdf:parseType="Resource">
          <vCard:Family>Lloyd</vCard:Family>
          <vCard:Given>Catherine</vCard:Given>
		  <vCard:Other>May</vCard:Other>
        </vCard:N>
        <vCard:EMAIL rdf:parseType="Resource">
          <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
          <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
        </vCard:EMAIL>
        <vCard:ORG rdf:parseType="Resource">
          <vCard:Orgname>The University of Auckland</vCard:Orgname>
          <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
        </vCard:ORG>
      </dc:creator>
    
      <dcterms:created rdf:parseType="Resource">
        <dcterms:W3CDTF>2010-01-20</dcterms:W3CDTF>
      </dcterms:created>
     
    </rdf:Description>
   
    <rdf:Description rdf:about="#michailova_2005">
      
	  <cmeta:bio_entity>ventricular myocyte</cmeta:bio_entity>
      
      <!--  Keyword(s) -->
      <bqs:reference rdf:parseType="Resource">
        <dc:subject rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>ventricular myocyte</rdf:li>
              <rdf:li>calcium dynamics</rdf:li>
              <rdf:li>electrophysiology</rdf:li>
              <rdf:li>cardiac</rdf:li>
              <rdf:li>buffering</rdf:li>
              <rdf:li>metabolism</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>

      <bqs:reference rdf:parseType="Resource">
        <bqs:Pubmed_id>15738467</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator>
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Michailova</vCard:Family>
                    <vCard:Given>A</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Saucerman</vCard:Family>
                    <vCard:Given>J</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Belik</vCard:Family>
                    <vCard:Given>M</vCard:Given>
					<vCard:Other>E</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>McCulloch</vCard:Family>
                    <vCard:Given>A</vCard:Given>
					<vCard:Other>D</vCard:Other>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title>
            Modeling regulation of cardiac KATP and L-type Ca2+ currents by ATP, ADP, and Mg2+
          </dc:title>
          <dcterms:issued rdf:parseType="Resource">
            <dcterms:W3CDTF>2005-03</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title>Biophysical Journal</dc:title>
          </bqs:Journal>
		  <bqs:volume>88</bqs:volume>
          <bqs:first_page>2234</bqs:first_page>
          <bqs:last_page>2249</bqs:last_page>
        </bqs:JournalArticle>
      </bqs:reference>
    </rdf:Description>
  </rdf:RDF>
</model>