Location: Iribe, Kohl, Noble, 2006 @ fc297fb5d36b / iribe_kohl_noble_2006.cellml

Author:
James Lawson <j.lawson@auckland.ac.nz>
Date:
2009-08-14 16:28:05+12:00
Desc:
added cmeta id's to cellml file to allow session to point to variables
Permanent Source URI:
https://models.cellml.org/workspace/iribe_kohl_noble_2006/rawfile/fc297fb5d36b2be7622a2040d7571f71d2e05d02/iribe_kohl_noble_2006.cellml

<?xml version="1.0"?>
<!--
This CellML file was generated on 27/11/2007 at 13:26:38 using:

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

CellML 1.0 was used to generate this cellular model
http://www.CellML.org/
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xml:base="file:///home/jlaw060/cellml_repository/pmr2/iribe_kohl_noble_2006/iribe_kohl_noble_2006.cellml" cmeta:id="iribe_kohl_noble_2006" name="iribe_kohl_noble_2006">
  <documentation xmlns="http://cellml.org/tmp-documentation">
    <article>
      <section id="sec_status">
        <title>Model Status</title>
        <para>
This model was created by Penny Noble of Oxford University.  This file is known to run COR and OpenCell. A PCEnv session file is also associated with this model.
          </para>
<para>
This model is valid CellML and has consistent units.
</para>
      </section>
      <sect1 id="sec_structure">
        <title>Model Structure</title>
        <para>
This model was based on the paper:
</para>
        <para>
Modulatory effect of calmodulin-dependent kinase II (CaMKII) on sarcoplasmic reticulum Ca2+ handling and interval-force relations: a modelling study.
</para>
        <para>
by Iribe G, Kohl P, Noble D
</para>
        <para>
Philos Transact A Math Phys Eng Sci. 2006 May 15;364(1842):1107-33
</para>
        <para>
ABSTRACT:
</para>
        <para>
We hypothesize that slow inactivation of Ca2+/calmodulin-dependent kinase II (CaMKII) and its modulatory effect on sarcoplasmic reticulum (SR) Ca2+ handling are important for various interval-force (I-F) relations, in particular for the beat interval dependency in transient alternans during the decay of post-extrasystolic potentiation. We have developed a mathematical model of a single cardiomyocyte to integrate various I-F relations, including alternans, by incorporating a conceptual CaMKII kinetics model into the SR Ca2+ handling model. Our model integrates I-F relations, such as the beat interval-dependent twitch force duration, restitution and potentiation, positive staircase phenomenon and alternans. We found that CaMKII affects more or less all I-F relations, and it is a key factor for integration of the various I-F relations in our model. Alternans arises, in the model, out of a steep relation between SR Ca2+ load and release, owing to SR load-dependent changes in the releasability of Ca2+ via the ryanodine receptor. Beat interval-dependent CaMKII activity, owing to its kinetic properties and amplifying effect on SR Ca2+ load dependency of Ca2+ release, replicated the beat interval dependency of alternans, as observed experimentally. Additionally, our model enabled reproduction of the effects of various interventions on alternans, such as the slowing or accelerating of Ca2+ release and/or uptake. We conclude that a slow time-dependent factor, represented in the model by CaMKII, is important for the integration of I-F relations, including alternans, and that our model offers a useful tool for further analysis of the roles of integrative Ca2+ handling in myocardial I-F relations.
</para>
      </sect1>
    </article>
  </documentation>
   <units name="micrometre">
      <unit prefix="micro" units="metre"/>
   </units>
   <units name="per_second">
      <unit exponent="-1" units="second"/>
   </units>
   <units name="millivolt">
      <unit prefix="milli" units="volt"/>
   </units>
   <units name="per_millivolt">
      <unit exponent="-1" prefix="milli" units="volt"/>
   </units>
   <units name="per_millivolt_second">
      <unit exponent="-1" units="millivolt"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units base_units="yes" name="arbitary"/>
   <units name="microS">
      <unit prefix="micro" units="siemens"/>
   </units>
   <units name="nanoS_per_cm2">
      <unit prefix="nano" units="siemens"/>
      <unit exponent="-2" prefix="centi" units="metre"/>
   </units>
   <units name="microF">
      <unit prefix="micro" units="farad"/>
   </units>
   <units name="nanoA">
      <unit prefix="nano" units="ampere"/>
   </units>
   <units name="nanoA_per_millimolar">
      <unit prefix="nano" units="ampere"/>
      <unit exponent="-1" units="millimolar"/>
   </units>
   <units name="nanoA_per_millimolar4">
      <unit prefix="nano" units="ampere"/>
      <unit exponent="-4" units="millimolar"/>
   </units>
   <units name="millimolar">
      <unit prefix="milli" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="micromolar">
      <unit prefix="micro" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="millimolar_per_second">
      <unit units="millimolar"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="per_millimolar_per_second">
      <unit exponent="-1" units="millimolar"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="per_millimolar2_per_second">
      <unit exponent="-2" units="millimolar"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="millijoule_per_mole_kelvin">
      <unit prefix="milli" units="joule"/>
      <unit exponent="-1" units="mole"/>
      <unit exponent="-1" units="kelvin"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit exponent="-1" units="mole"/>
   </units>
   <units name="per_nanocoulomb">
      <unit exponent="-1" prefix="nano" units="coulomb"/>
   </units>
   <units name="cm_per_second">
      <unit prefix="centi" units="metre"/>
      <unit exponent="-1" units="second"/>
   </units>
   <units name="cm2">
      <unit exponent="2" prefix="centi" units="metre"/>
   </units>
   <units name="microlitre">
      <unit prefix="micro" units="litre"/>
   </units>
   <units name="N_per_mm2">
      <unit units="newton"/>
      <unit exponent="-2" prefix="milli" units="metre"/>
   </units>
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="second"/>
   </component>
   <component name="cell_parameters">
      <variable initial_value="8314.472" name="R" public_interface="out" units="millijoule_per_mole_kelvin"/>
      <variable initial_value="310" name="T" public_interface="out" units="kelvin"/>
      <variable initial_value="96485.3415" name="F" public_interface="out" units="coulomb_per_mole"/>
      <variable initial_value="9.5e-5" name="Cm" public_interface="out" units="microF"/>
      <variable initial_value="1.6404e-5" name="v_i" public_interface="out" units="microlitre"/>
      <variable initial_value="3.3477e-6" name="v_SR" public_interface="out" units="microlitre"/>
      <variable initial_value="140" name="Na_o" public_interface="out" units="millimolar"/>
      <variable initial_value="4" name="K_o" public_interface="out" units="millimolar"/>
      <variable initial_value="2" name="Ca_o" public_interface="out" units="millimolar"/>
   </component>
   <component name="membrane_potential">
      <variable cmeta:id="membrane_potential_V" initial_value="-92.849333" name="V" public_interface="out" units="millivolt"/>
      <variable name="Cm" public_interface="in" units="microF"/>
      <variable name="i_Na" public_interface="in" units="nanoA"/>
      <variable name="i_K1" public_interface="in" units="nanoA"/>
      <variable name="i_to" public_interface="in" units="nanoA"/>
      <variable name="i_K" public_interface="in" units="nanoA"/>
      <variable name="i_Ca_L" public_interface="in" units="nanoA"/>
      <variable name="i_NaK" public_interface="in" units="nanoA"/>
      <variable name="i_NaCa" public_interface="in" units="nanoA"/>
      <variable name="i_b_Na" public_interface="in" units="nanoA"/>
      <variable name="i_b_K" public_interface="in" units="nanoA"/>
      <variable name="i_b_Ca" public_interface="in" units="nanoA"/>
      <variable name="i_Stim" units="nanoA"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable initial_value="0" name="stim_start" units="second"/>
      <variable initial_value="1000" name="stim_end" units="second"/>
      <variable initial_value="0.5" name="stim_period" units="second"/>
      <variable initial_value="0.002" name="stim_duration" units="second"/>
      <variable initial_value="-4" name="stim_amplitude" units="nanoA"/>
      <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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="nanoA">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  </apply>
                  <ci>Cm</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>i_Na</ci>
                  <ci>i_b_Na</ci>
                  <ci>i_K1</ci>
                  <ci>i_K</ci>
                  <ci>i_to</ci>
                  <ci>i_b_K</ci>
                  <ci>i_Ca_L</ci>
                  <ci>i_b_Ca</ci>
                  <ci>i_NaCa</ci>
                  <ci>i_NaK</ci>
                  <ci>i_Stim</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="reversal_potentials">
      <variable name="E_Na" public_interface="out" units="millivolt"/>
      <variable name="E_K" public_interface="out" units="millivolt"/>
      <variable name="E_Ca" public_interface="out" units="millivolt"/>
      <variable name="E_mh" public_interface="out" units="millivolt"/>
      <variable name="K_o" public_interface="in" units="millimolar"/>
      <variable name="Na_o" public_interface="in" units="millimolar"/>
      <variable name="K_i" public_interface="in" units="millimolar"/>
      <variable name="Na_i" public_interface="in" units="millimolar"/>
      <variable name="Ca_o" public_interface="in" units="millimolar"/>
      <variable name="Ca_i" public_interface="in" units="millimolar"/>
      <variable name="R" public_interface="in" units="millijoule_per_mole_kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <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>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>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.5</cn>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ca_o</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>E_mh</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>Na_o</ci>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.12</cn>
                           <ci>K_o</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>Na_i</ci>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.12</cn>
                           <ci>K_i</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable name="i_Na" public_interface="out" units="nanoA"/>
      <variable initial_value="2.5" name="g_Na" units="microS"/>
      <variable name="E_mh" public_interface="in" units="millivolt"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="m" private_interface="in" units="dimensionless"/>
      <variable name="h" private_interface="in" units="dimensionless"/>
      <variable name="time" private_interface="out" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_mh</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable initial_value="0.0013809" name="m" public_interface="out" units="dimensionless"/>
      <variable name="alpha_m" units="per_second"/>
      <variable name="beta_m" units="per_second"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable initial_value="1e-5" name="delta_m" units="millivolt"/>
      <variable name="E0_m" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_m</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">41</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">2000</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_m</ci>
                     </apply>
                     <ci>delta_m</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt_second">200</cn>
                        <ci>E0_m</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.1</cn>
                              </apply>
                              <ci>E0_m</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">8000</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.056</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">66</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_m</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>m</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_m</ci>
                  <ci>m</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h_gate">
      <variable initial_value="0.99569" name="h" public_interface="out" units="dimensionless"/>
      <variable name="alpha_h" units="per_second"/>
      <variable name="beta_h" units="per_second"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_h</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">20</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.125</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">75</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">2000</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">320</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">75</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_h</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>h</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_h</ci>
                  <ci>h</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_dependent_rectifier_potassium_current">
      <variable name="i_K" public_interface="out" units="nanoA"/>
      <variable initial_value="1" name="i_Kmax" units="nanoA"/>
      <variable name="K_i" public_interface="in" units="millimolar"/>
      <variable name="K_o" public_interface="in" units="millimolar"/>
      <variable name="R" public_interface="in" units="millijoule_per_mole_kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="x" private_interface="in" units="dimensionless"/>
      <variable name="time" private_interface="out" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_Kmax</ci>
                  <ci>x</ci>
                  <apply>
                     <minus/>
                     <ci>K_i</ci>
                     <apply>
                        <times/>
                        <ci>K_o</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                 </apply>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">140</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_dependent_rectifier_potassium_current_x_gate">
      <variable initial_value="5.1127e-2" name="x" public_interface="out" units="dimensionless"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="alpha_x" units="per_second"/>
      <variable name="beta_x" units="per_second"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_x</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">0.5</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.0826</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.057</cn>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_x</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">1.3</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.06</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">20</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.04</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">20</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>x</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_x</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>x</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_x</ci>
                  <ci>x</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="time_independent_potassium_current">
      <variable name="i_K1" public_interface="out" units="nanoA"/>
      <variable initial_value="10" name="K_mk1" units="millimolar"/>
      <variable initial_value="1" name="g_K1" units="microS"/>
      <variable name="K_o" public_interface="in" units="millimolar"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
      <variable name="R" public_interface="in" units="millijoule_per_mole_kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>g_K1</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>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                           <ci>F</ci>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <ci>E_K</ci>
                              </apply>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current">
      <variable name="i_to" public_interface="out" units="nanoA"/>
      <variable initial_value="0.005" name="g_to" units="microS"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="s" private_interface="in" units="dimensionless"/>
      <variable name="r" private_interface="in" units="dimensionless"/>
      <variable name="time" private_interface="out" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_to</ci>
            <apply>
               <times/>
               <ci>g_to</ci>
               <ci>s</ci>
               <ci>r</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_s_gate">
      <variable initial_value="0.95854" name="s" public_interface="out" units="dimensionless"/>
      <variable name="alpha_s" units="per_second"/>
      <variable name="beta_s" units="per_second"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha_s</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">0.033</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">17</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_s</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">33</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt">0.125</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_s</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>s</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_s</ci>
                  <ci>s</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="transient_outward_current_r_gate">
      <variable initial_value="1.5185e-8" name="r" public_interface="out" units="dimensionless"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r</ci>
            </apply>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">333</cn>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <plus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
                                 </apply>
                              </apply>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <ci>r</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel">
      <variable cmeta:id="L_type_Ca_channel_i_Ca_L" name="i_Ca_L" public_interface="out" units="nanoA"/>
      <variable name="i_Ca_L_Na" public_interface="out" units="nanoA"/>
      <variable name="i_Ca_L_K" public_interface="out" units="nanoA"/>
      <variable name="i_Ca_L_Ca" public_interface="out" units="nanoA"/>
      <variable initial_value="0.25" name="P_Ca_L_Ca" units="nanoA_per_millimolar"/>
      <variable name="R" public_interface="in" units="millijoule_per_mole_kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="Ca_o" public_interface="in" units="millimolar"/>
      <variable name="Ca_i" private_interface="out" public_interface="in" units="millimolar"/>
      <variable name="K_o" public_interface="in" units="millimolar"/>
      <variable name="K_i" public_interface="in" units="millimolar"/>
      <variable name="Na_o" public_interface="in" units="millimolar"/>
      <variable name="Na_i" public_interface="in" units="millimolar"/>
      <variable name="d" private_interface="in" units="dimensionless"/>
      <variable name="f" private_interface="in" units="dimensionless"/>
      <variable name="time" private_interface="out" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca_L_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">4</cn>
                     <ci>d</ci>
                     <ci>f</ci>
                     <ci>P_Ca_L_Ca</ci>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                     </apply>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                                 </apply>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Ca_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">100</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Ca_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                              </apply>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_K</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.002</cn>
                     <ci>d</ci>
                     <ci>f</ci>
                     <ci>P_Ca_L_Ca</ci>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                     </apply>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <apply>
                                       <minus/>
                                       <ci>V</ci>
                                       <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                                    </apply>
                                 </apply>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>K_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.01</cn>
                     <ci>d</ci>
                     <ci>f</ci>
                     <ci>P_Ca_L_Ca</ci>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                     </apply>
                     <ci>F</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <apply>
                                       <minus/>
                                       <ci>V</ci>
                                       <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                                    </apply>
                                 </apply>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>Na_i</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Na_o</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>V</ci>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">50</cn>
                                 </apply>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>i_Ca_L</ci>
            <apply>
               <plus/>
               <ci>i_Ca_L_Ca</ci>
               <ci>i_Ca_L_K</ci>
               <ci>i_Ca_L_Na</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_d_gate">
      <variable initial_value="1.7908e-8" name="d" public_interface="out" units="dimensionless"/>
      <variable name="alpha_d" units="per_second"/>
      <variable name="beta_d" units="per_second"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="E0_d" units="millivolt"/>
      <variable initial_value="3" name="speed_d" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_d</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">24</cn>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_d</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <ci>speed_d</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">120</cn>
                  </apply>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">0.00001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>speed_d</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt_second">30</cn>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>E0_d</ci>
                              </apply>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_d</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <ci>speed_d</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">120</cn>
                  </apply>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_d</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">0.00001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>speed_d</ci>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt_second">12</cn>
                        </apply>
                        <ci>E0_d</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_d</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>d</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_d</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>d</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_d</ci>
                  <ci>d</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="L_type_Ca_channel_f_gate">
      <variable initial_value="1" name="f" public_interface="out" units="dimensionless"/>
      <variable name="alpha_f" units="per_second"/>
      <variable name="beta_f" units="per_second"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable initial_value="0.5" name="speed_f" units="dimensionless"/>
      <variable name="E0_f" units="millivolt"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>E0_f</ci>
            <apply>
               <plus/>
               <ci>V</ci>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">34</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_f</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <ci>speed_f</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">25</cn>
                  </apply>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>E0_f</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">0.00001</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>speed_f</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millivolt_second">6.25</cn>
                        <ci>E0_f</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <minus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        </apply>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <ci>E0_f</ci>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_f</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>speed_f</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">50</cn>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>E0_f</ci>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>f</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_f</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>f</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_f</ci>
                  <ci>f</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_calcium_exchanger">
      <variable name="i_NaCa" public_interface="out" units="nanoA"/>
      <variable initial_value="0.0005" name="i_NaCa_max" units="nanoA_per_millimolar4"/>
      <variable initial_value="0.5" name="gamma" units="dimensionless"/>
      <variable name="Na_i" public_interface="in" units="millimolar"/>
      <variable name="Na_o" public_interface="in" units="millimolar"/>
      <variable name="Ca_i" public_interface="in" units="millimolar"/>
      <variable name="Ca_o" public_interface="in" units="millimolar"/>
      <variable name="R" public_interface="in" units="millijoule_per_mole_kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaCa</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>i_NaCa_max</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>gamma</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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Ca_o</ci>
                     </apply>
                     <apply>
                        <times/>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>gamma</ci>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Ca_i</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>Ca_i</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">0.0069</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_potassium_pump">
      <variable name="i_NaK" public_interface="out" units="nanoA"/>
      <variable initial_value="1.36" name="i_NaK_max" units="nanoA"/>
      <variable initial_value="1" name="K_mK" units="millimolar"/>
      <variable initial_value="21.7" name="K_mNa" units="millimolar"/>
      <variable name="K_o" public_interface="in" units="millimolar"/>
      <variable name="Na_i" public_interface="in" units="millimolar"/>
      <variable name="R" public_interface="in" units="millijoule_per_mole_kelvin"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="T" public_interface="in" units="kelvin"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>i_NaK_max</ci>
                              <ci>K_o</ci>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>K_mK</ci>
                              <ci>K_o</ci>
                           </apply>
                        </apply>
                        <ci>Na_i</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_mNa</ci>
                        <ci>Na_i</ci>
                     </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1245</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" 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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.0353</cn>
                     <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>
      </math>
   </component>
   <component name="calcium_background_current">
      <variable name="i_b_Ca" public_interface="out" units="nanoA"/>
      <variable initial_value="0.00025" name="g_b_Ca" units="microS"/>
      <variable name="E_Ca" public_interface="in" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_Ca</ci>
            <apply>
               <times/>
               <ci>g_b_Ca</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="potassium_background_current">
      <variable name="i_b_K" public_interface="out" units="nanoA"/>
      <variable initial_value="0.0006" name="g_b_K" units="microS"/>
      <variable name="E_K" public_interface="in" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_K</ci>
            <apply>
               <times/>
               <ci>g_b_K</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_background_current">
      <variable name="i_b_Na" public_interface="out" units="nanoA"/>
      <variable initial_value="0.0006" name="g_b_Na" units="microS"/>
      <variable name="E_Na" public_interface="in" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_b_Na</ci>
            <apply>
               <times/>
               <ci>g_b_Na</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="CaMKII_factor">
      <variable initial_value="1.028" name="F_CaMK" public_interface="out" units="dimensionless"/>
      <variable name="Inf_CaMK" units="dimensionless"/>
      <variable initial_value="0.8" name="Tau_CaMK" units="second"/>
      <variable name="Cmdn_Ca" public_interface="in" units="millimolar"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Inf_CaMK</ci>
            <apply>
               <divide/>
               <ci>Cmdn_Ca</ci>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">0.00005</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>F_CaMK</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Inf_CaMK</ci>
                  <ci>F_CaMK</ci>
               </apply>
               <ci>Tau_CaMK</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="RyR">
      <variable name="j_rel" public_interface="out" units="millimolar_per_second"/>
      <variable name="K_rel" units="per_second"/>
      <variable name="F_rel" units="dimensionless"/>
      <variable name="Ca_i" public_interface="in" units="millimolar"/>
      <variable name="Ca_SR" public_interface="in" units="millimolar"/>
      <variable name="i_Ca_L" public_interface="in" units="nanoA"/>
      <variable name="F_CaMK" public_interface="in" units="dimensionless"/>
      <variable initial_value="500" name="K_rel_max" units="per_second"/>
      <variable initial_value="0.25089" name="F_SRCa_RyR" units="millimolar"/>
      <variable initial_value="0.05" name="Tau_SRCa_RyR" units="second"/>
      <variable name="N_CaMK" units="dimensionless"/>
      <variable initial_value="1" name="gain_k1" units="dimensionless"/>
      <variable initial_value="1" name="gain_k2" units="dimensionless"/>
      <variable initial_value="1" name="gain_k3" units="dimensionless"/>
      <variable initial_value="1" name="gain_k4" units="dimensionless"/>
      <variable name="k_1" units="per_second"/>
      <variable name="k_2" units="per_second"/>
      <variable name="k_3" units="per_second"/>
      <variable name="k_4" units="per_second"/>
      <variable initial_value="0.5268" name="F_1" units="dimensionless"/>
      <variable initial_value="8.7508e-6" name="F_2" units="dimensionless"/>
      <variable name="F_3" units="dimensionless"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable initial_value="0" name="K_leak_rate" units="per_second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>N_CaMK</ci>
            <apply>
               <power/>
               <apply>
                  <divide/>
                  <ci>F_CaMK</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.7</cn>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k_1</ci>
            <apply>
               <times/>
               <ci>gain_k1</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millimolar2_per_second">30625000</cn>
                     <apply>
                        <power/>
                        <ci>Ca_i</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_nanocoulomb">245</cn>
                     <ci>i_Ca_L</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k_2</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>gain_k2</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">450</cn>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">0.36</cn>
                     <ci>Ca_SR</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k_3</ci>
            <apply>
               <times/>
               <ci>gain_k3</ci>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">1.885</cn>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>F_SRCa_RyR</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">0.22</cn>
                  </apply>
                  <ci>N_CaMK</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k_4</ci>
            <apply>
               <times/>
               <ci>gain_k4</ci>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">1.8</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>F_1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>k_3</ci>
                     <ci>F_3</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_4</ci>
                     <ci>F_1</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_1</ci>
                  <ci>F_1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>F_2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_1</ci>
                  <ci>F_1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_2</ci>
                  <ci>F_2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>F_3</ci>
            <apply>
               <minus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>F_1</ci>
                  <ci>F_2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>F_rel</ci>
            <apply>
               <power/>
               <apply>
                  <divide/>
                  <ci>F_2</ci>
                  <apply>
                     <plus/>
                     <ci>F_2</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.25</cn>
                  </apply>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>F_SRCa_RyR</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Ca_SR</ci>
                  <ci>F_SRCa_RyR</ci>
               </apply>
               <ci>Tau_SRCa_RyR</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K_rel</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>K_rel_max</ci>
                  <ci>F_SRCa_RyR</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>F_SRCa_RyR</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">0.2</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>j_rel</ci>
            <apply>
               <times/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>K_rel</ci>
                     <ci>F_rel</ci>
                  </apply>
                  <ci>K_leak_rate</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>Ca_SR</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="SERCA">
      <variable name="j_up" public_interface="out" units="millimolar_per_second"/>
      <variable name="F_CaMK" public_interface="in" units="dimensionless"/>
      <variable initial_value="0.292" name="V_max_f" units="millimolar_per_second"/>
      <variable initial_value="0.391" name="V_max_r" units="millimolar_per_second"/>
      <variable name="f_b" units="dimensionless"/>
      <variable name="r_b" units="dimensionless"/>
      <variable name="Ca_i" public_interface="in" units="millimolar"/>
      <variable name="Ca_SR" public_interface="in" units="millimolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_b</ci>
            <apply>
               <power/>
               <apply>
                  <divide/>
                  <ci>Ca_i</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">0.00024</cn>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>r_b</ci>
            <apply>
               <power/>
               <apply>
                  <divide/>
                  <ci>Ca_SR</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">1.64</cn>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>j_up</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>F_CaMK</ci>
                     <ci>V_max_f</ci>
                     <ci>f_b</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>V_max_r</ci>
                     <ci>r_b</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <ci>f_b</ci>
                  <ci>r_b</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calmodulin">
      <variable initial_value="3.9636e-6" name="Cmdn_Ca" public_interface="out" units="millimolar"/>
      <variable name="Ca_i" public_interface="in" units="millimolar"/>
      <variable initial_value="0.02" name="Cmdn_tot" units="millimolar"/>
      <variable initial_value="10000" name="alpha_cmdn" units="per_millimolar_per_second"/>
      <variable initial_value="500" name="beta_cmdn" units="per_second"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="dCmdn_Ca_dtime" public_interface="out" units="millimolar_per_second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>dCmdn_Ca_dtime</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_cmdn</ci>
                  <apply>
                     <minus/>
                     <ci>Cmdn_tot</ci>
                     <ci>Cmdn_Ca</ci>
                  </apply>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_cmdn</ci>
                  <ci>Cmdn_Ca</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cmdn_Ca</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_cmdn</ci>
                  <apply>
                     <minus/>
                     <ci>Cmdn_tot</ci>
                     <ci>Cmdn_Ca</ci>
                  </apply>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_cmdn</ci>
                  <ci>Cmdn_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="troponin">
      <variable initial_value="2.7661e-4" name="Trpn_Ca" public_interface="out" units="millimolar"/>
      <variable name="Ca_i" public_interface="in" units="millimolar"/>
      <variable initial_value="0.07" name="Trpn_tot" public_interface="out" units="millimolar"/>
      <variable initial_value="80000" name="alpha_trpn" public_interface="out" units="per_millimolar_per_second"/>
      <variable initial_value="200" name="beta_trpn" public_interface="out" units="per_second"/>
      <variable name="Force_norm" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="dTrpn_Ca_dtime" public_interface="out" units="millimolar_per_second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>dTrpn_Ca_dtime</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_trpn</ci>
                  <apply>
                     <minus/>
                     <ci>Trpn_tot</ci>
                     <ci>Trpn_Ca</ci>
                  </apply>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>beta_trpn</ci>
                        <apply>
                           <plus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                           <apply>
                              <times/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                              <apply>
                                 <minus/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                                 <ci>Force_norm</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
                  <ci>Trpn_Ca</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Trpn_Ca</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_trpn</ci>
                  <apply>
                     <minus/>
                     <ci>Trpn_tot</ci>
                     <ci>Trpn_Ca</ci>
                  </apply>
                  <ci>Ca_i</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>beta_trpn</ci>
                        <apply>
                           <plus/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                           <apply>
                              <times/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                              <apply>
                                 <minus/>
                                 <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                                 <ci>Force_norm</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  </apply>
                  <ci>Trpn_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_calcium_concentration">
      <variable initial_value="9.91e-6" name="Ca_i" public_interface="out" units="millimolar"/>
      <variable name="i_Ca_L_Ca" public_interface="in" units="nanoA"/>
      <variable name="i_b_Ca" public_interface="in" units="nanoA"/>
      <variable name="i_NaCa" public_interface="in" units="nanoA"/>
      <variable name="Cmdn_Ca" public_interface="in" units="millimolar"/>
      <variable name="Trpn_Ca" public_interface="in" units="millimolar"/>
      <variable name="j_rel" public_interface="in" units="millimolar_per_second"/>
      <variable name="j_up" public_interface="in" units="millimolar_per_second"/>
      <variable name="v_SR" public_interface="in" units="microlitre"/>
      <variable name="v_i" public_interface="in" units="microlitre"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="time" public_interface="in" units="second"/>
      <variable name="dCmdn_Ca_dtime" public_interface="in" units="millimolar_per_second"/>
      <variable name="dTrpn_Ca_dtime" public_interface="in" units="millimolar_per_second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_i</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <apply>
                        <minus/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>i_Ca_L_Ca</ci>
                                    <ci>i_b_Ca</ci>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                                    <ci>i_NaCa</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                              <ci>v_i</ci>
                              <ci>F</ci>
                           </apply>
                        </apply>
                        <ci>j_up</ci>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>j_rel</ci>
                           <ci>v_SR</ci>
                        </apply>
                        <ci>v_i</ci>
                     </apply>
                  </apply>
                  <ci>dCmdn_Ca_dtime</ci>
               </apply>
               <ci>dTrpn_Ca_dtime</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="SR_calcium_concentration">
      <variable initial_value="0.24886" name="Ca_SR" public_interface="out" units="millimolar"/>
      <variable name="j_rel" public_interface="in" units="millimolar_per_second"/>
      <variable name="j_up" public_interface="in" units="millimolar_per_second"/>
      <variable name="v_SR" public_interface="in" units="microlitre"/>
      <variable name="v_i" public_interface="in" units="microlitre"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_SR</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>j_up</ci>
                     <ci>v_i</ci>
                  </apply>
                  <ci>v_SR</ci>
               </apply>
               <ci>j_rel</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_sodium_concentration">
      <variable initial_value="5.8041" name="Na_i" public_interface="out" units="millimolar"/>
      <variable name="i_Na" public_interface="in" units="nanoA"/>
      <variable name="i_b_Na" public_interface="in" units="nanoA"/>
      <variable name="i_Ca_L_Na" public_interface="in" units="nanoA"/>
      <variable name="i_NaCa" public_interface="in" units="nanoA"/>
      <variable name="i_NaK" public_interface="in" units="nanoA"/>
      <variable name="v_i" public_interface="in" units="microlitre"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Na_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_Na</ci>
                     <ci>i_b_Na</ci>
                     <ci>i_Ca_L_Na</ci>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        <ci>i_NaCa</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        <ci>i_NaK</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>v_i</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="intracellular_potassium_concentration">
      <variable initial_value="138.22" name="K_i" public_interface="out" units="millimolar"/>
      <variable name="i_K1" public_interface="in" units="nanoA"/>
      <variable name="i_K" public_interface="in" units="nanoA"/>
      <variable name="i_to" public_interface="in" units="nanoA"/>
      <variable name="i_b_K" public_interface="in" units="nanoA"/>
      <variable name="i_Ca_L_K" public_interface="in" units="nanoA"/>
      <variable name="i_NaK" public_interface="in" units="nanoA"/>
      <variable name="v_i" public_interface="in" units="microlitre"/>
      <variable name="F" public_interface="in" units="coulomb_per_mole"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>K_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>i_K1</ci>
                        <ci>i_K</ci>
                        <ci>i_to</ci>
                        <ci>i_b_K</ci>
                        <ci>i_Ca_L_K</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        <ci>i_NaK</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>v_i</ci>
                  <ci>F</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Force">
      <variable cmeta:id="Force_Force" name="Force" units="N_per_mm2"/>
      <variable initial_value="0.1" name="zeta" units="N_per_mm2"/>
      <variable name="Force_max" units="dimensionless"/>
      <variable name="phi_SL" units="dimensionless"/>
      <variable name="P_1_max" units="dimensionless"/>
      <variable name="P_2_max" units="dimensionless"/>
      <variable name="P_3_max" units="dimensionless"/>
      <variable name="sigma_paths" units="dimensionless"/>
      <variable initial_value="0.99917" name="N_0" units="dimensionless"/>
      <variable initial_value="9.8593e-5" name="P_0" units="dimensionless"/>
      <variable initial_value="1.3331e-4" name="P_1" units="dimensionless"/>
      <variable initial_value="2.3505e-4" name="P_2" units="dimensionless"/>
      <variable initial_value="1.5349e-4" name="P_3" units="dimensionless"/>
      <variable name="N_1" units="dimensionless"/>
      <variable name="alpha_tm" units="per_second"/>
      <variable initial_value="40" name="beta_tm" units="per_second"/>
      <variable name="K_tm" units="dimensionless"/>
      <variable name="N_tm" units="dimensionless"/>
      <variable initial_value="2.15" name="SL" units="micrometre"/>
      <variable name="SL_norm" units="dimensionless"/>
      <variable name="f_01" units="per_second"/>
      <variable name="f_12" units="per_second"/>
      <variable name="f_23" units="per_second"/>
      <variable name="g_01" units="per_second"/>
      <variable name="g_12" units="per_second"/>
      <variable name="g_23" units="per_second"/>
      <variable initial_value="10" name="f_XB" units="per_second"/>
      <variable initial_value="30" name="g_XB" units="per_second"/>
      <variable name="Trpn_Ca" public_interface="in" units="millimolar"/>
      <variable name="Trpn_tot" public_interface="in" units="millimolar"/>
      <variable name="alpha_trpn" public_interface="in" units="per_millimolar_per_second"/>
      <variable name="beta_trpn" public_interface="in" units="per_second"/>
      <variable name="Force_norm" public_interface="out" units="dimensionless"/>
      <variable name="time" public_interface="in" units="second"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>f_01</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
               <ci>f_XB</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_12</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">10</cn>
               <ci>f_XB</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_23</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">7</cn>
               <ci>f_XB</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>g_01</ci>
            <apply>
               <times/>
               <ci>g_XB</ci>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>SL_norm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>g_12</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
               <ci>g_XB</ci>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>SL_norm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>g_23</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
               <ci>g_XB</ci>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>SL_norm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>SL_norm</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>SL</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">1.7</cn>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">0.7</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>N_tm</ci>
            <apply>
               <plus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3.5</cn>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2.5</cn>
                  <ci>SL_norm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K_tm</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <divide/>
                        <ci>beta_trpn</ci>
                        <ci>alpha_trpn</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">0.0017</cn>
                        <apply>
                           <times/>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millimolar">0.0009</cn>
                           <ci>SL_norm</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>alpha_tm</ci>
            <apply>
               <times/>
               <ci>beta_tm</ci>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>Trpn_Ca</ci>
                     <apply>
                        <times/>
                        <ci>Trpn_tot</ci>
                        <ci>K_tm</ci>
                     </apply>
                  </apply>
                  <ci>N_tm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>N_0</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>beta_tm</ci>
                     <ci>P_0</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>alpha_tm</ci>
                     <ci>N_0</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>g_01</ci>
                  <ci>N_1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_0</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>beta_tm</ci>
                        <ci>f_01</ci>
                     </apply>
                  </apply>
                  <ci>P_0</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>alpha_tm</ci>
                  <ci>N_0</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>g_01</ci>
                  <ci>P_1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_1</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>beta_tm</ci>
                        <ci>f_12</ci>
                        <ci>g_01</ci>
                     </apply>
                  </apply>
                  <ci>P_1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>alpha_tm</ci>
                  <ci>N_1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>f_01</ci>
                  <ci>P_0</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>g_12</ci>
                  <ci>P_2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_2</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>f_23</ci>
                        <ci>g_12</ci>
                     </apply>
                  </apply>
                  <ci>P_2</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>f_12</ci>
                  <ci>P_1</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>g_23</ci>
                  <ci>P_3</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_3</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>g_23</ci>
                  </apply>
                  <ci>P_3</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>f_23</ci>
                  <ci>P_2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>N_1</ci>
            <apply>
               <minus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>N_0</ci>
                  <ci>P_0</ci>
                  <ci>P_1</ci>
                  <ci>P_2</ci>
                  <ci>P_3</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>sigma_paths</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>g_XB</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">2</cn>
                  <ci>g_XB</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">3</cn>
                  <ci>g_XB</ci>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_01</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">2</cn>
                  <ci>g_XB</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">3</cn>
                  <ci>g_XB</ci>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_01</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_12</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">3</cn>
                  <ci>g_XB</ci>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_01</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_12</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_23</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_1_max</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_01</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">2</cn>
                  <ci>g_XB</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">3</cn>
                  <ci>g_XB</ci>
               </apply>
               <ci>sigma_paths</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_2_max</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_01</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_12</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">3</cn>
                  <ci>g_XB</ci>
               </apply>
               <ci>sigma_paths</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P_3_max</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_01</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_12</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
                  <ci>f_23</ci>
               </apply>
               <ci>sigma_paths</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Force_max</ci>
            <apply>
               <plus/>
               <ci>P_1_max</ci>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>P_2_max</ci>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                  <ci>P_3_max</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>phi_SL</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>SL</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">0.6</cn>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">1.4</cn>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>SL</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">1.7</cn>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>SL</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">2</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <and/>
                     <apply>
                        <gt/>
                        <ci>SL</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">2</cn>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>SL</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">2.2</cn>
                     </apply>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">3.6</cn>
                        <ci>SL</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">1.4</cn>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <gt/>
                        <ci>SL</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">2.2</cn>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>SL</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="micrometre">2.3</cn>
                     </apply>
                  </apply>
               </piece>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>Force_norm</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>phi_SL</ci>
                  <apply>
                     <plus/>
                     <ci>P_1</ci>
                     <ci>N_1</ci>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                        <ci>P_2</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        <ci>P_3</ci>
                     </apply>
                  </apply>
               </apply>
               <ci>Force_max</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Force</ci>
            <apply>
               <times/>
               <ci>zeta</ci>
               <ci>Force_norm</ci>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="Force">
         <component_ref component="troponin"/>
      </component_ref>
      <component_ref component="CaMKII_factor">
         <component_ref component="membrane_potential"/>
         <component_ref component="RyR"/>
         <component_ref component="SERCA"/>
      </component_ref>
      <component_ref component="intracellular_calcium_concentration">
         <component_ref component="SR_calcium_concentration"/>
         <component_ref component="intracellular_sodium_concentration"/>
         <component_ref component="intracellular_potassium_concentration"/>
         <component_ref component="calmodulin"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
         <component_ref component="L_type_Ca_channel_d_gate"/>
         <component_ref component="L_type_Ca_channel_f_gate"/>
      </component_ref>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_h_gate"/>
      </component_ref>
      <component_ref component="time_dependent_rectifier_potassium_current">
         <component_ref component="time_dependent_rectifier_potassium_current_x_gate"/>
      </component_ref>
      <component_ref component="transient_outward_current">
         <component_ref component="transient_outward_current_s_gate"/>
         <component_ref component="transient_outward_current_r_gate"/>
      </component_ref>
      <component_ref component="sodium_background_current">
         <component_ref component="calcium_background_current"/>
         <component_ref component="potassium_background_current"/>
      </component_ref>
   </group>
   <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_ref component="time_dependent_rectifier_potassium_current">
         <component_ref component="time_dependent_rectifier_potassium_current_x_gate"/>
      </component_ref>
      <component_ref component="transient_outward_current">
         <component_ref component="transient_outward_current_s_gate"/>
         <component_ref component="transient_outward_current_r_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
         <component_ref component="L_type_Ca_channel_d_gate"/>
         <component_ref component="L_type_Ca_channel_f_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="environment" component_2="membrane_potential"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="fast_sodium_current"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="time_dependent_rectifier_potassium_current"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="transient_outward_current"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="L_type_Ca_channel"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="CaMKII_factor"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="RyR"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="calmodulin"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="troponin"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="SR_calcium_concentration"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="Force"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="membrane_potential"/>
      <map_variables variable_1="Cm" variable_2="Cm"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="reversal_potentials"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="time_dependent_rectifier_potassium_current"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="time_independent_potassium_current"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="L_type_Ca_channel"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="sodium_calcium_exchanger"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="Na_o" variable_2="Na_o"/>
      <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="sodium_potassium_pump"/>
      <map_variables variable_1="R" variable_2="R"/>
      <map_variables variable_1="T" variable_2="T"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="K_o" variable_2="K_o"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="v_SR" variable_2="v_SR"/>
      <map_variables variable_1="v_i" variable_2="v_i"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="SR_calcium_concentration"/>
      <map_variables variable_1="v_SR" variable_2="v_SR"/>
      <map_variables variable_1="v_i" variable_2="v_i"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="v_i" variable_2="v_i"/>
   </connection>
   <connection>
      <map_components component_1="cell_parameters" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="F" variable_2="F"/>
      <map_variables variable_1="v_i" variable_2="v_i"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="fast_sodium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="time_dependent_rectifier_potassium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_K" variable_2="i_K"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="time_independent_potassium_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="transient_outward_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_to" variable_2="i_to"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="L_type_Ca_channel"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Ca_L" variable_2="i_Ca_L"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="sodium_calcium_exchanger"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="sodium_potassium_pump"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="calcium_background_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_b_Ca" variable_2="i_b_Ca"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="potassium_background_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_b_K" variable_2="i_b_K"/>
   </connection>
   <connection>
      <map_components component_1="membrane_potential" component_2="sodium_background_current"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_b_Na" variable_2="i_b_Na"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="fast_sodium_current"/>
      <map_variables variable_1="E_mh" variable_2="E_mh"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="time_independent_potassium_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="transient_outward_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="calcium_background_current"/>
      <map_variables variable_1="E_Ca" variable_2="E_Ca"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="potassium_background_current"/>
      <map_variables variable_1="E_K" variable_2="E_K"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="sodium_background_current"/>
      <map_variables variable_1="E_Na" variable_2="E_Na"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
   </connection>
   <connection>
      <map_components component_1="reversal_potentials" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_Na" variable_2="i_Na"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_m_gate"/>
      <map_variables variable_1="m" variable_2="m"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_sodium_current" component_2="fast_sodium_current_h_gate"/>
      <map_variables variable_1="h" variable_2="h"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="time_dependent_rectifier_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_K" variable_2="i_K"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
   </connection>
   <connection>
      <map_components component_1="time_dependent_rectifier_potassium_current" component_2="time_dependent_rectifier_potassium_current_x_gate"/>
      <map_variables variable_1="x" variable_2="x"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="time_independent_potassium_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_K1" variable_2="i_K1"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_to" variable_2="i_to"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="transient_outward_current_s_gate"/>
      <map_variables variable_1="s" variable_2="s"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="transient_outward_current" component_2="transient_outward_current_r_gate"/>
      <map_variables variable_1="r" variable_2="r"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="i_Ca_L_Ca" variable_2="i_Ca_L_Ca"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="K_i" variable_2="K_i"/>
      <map_variables variable_1="i_Ca_L_K" variable_2="i_Ca_L_K"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
      <map_variables variable_1="i_Ca_L_Na" variable_2="i_Ca_L_Na"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="RyR"/>
      <map_variables variable_1="i_Ca_L" variable_2="i_Ca_L"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_d_gate"/>
      <map_variables variable_1="d" variable_2="d"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_f_gate"/>
      <map_variables variable_1="f" variable_2="f"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_1="sodium_calcium_exchanger" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="Na_i" variable_2="Na_i"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
   </connection>
   <connection>
      <map_components component_1="sodium_potassium_pump" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_NaK" variable_2="i_NaK"/>
   </connection>
   <connection>
      <map_components component_1="calcium_background_current" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="i_b_Ca" variable_2="i_b_Ca"/>
   </connection>
   <connection>
      <map_components component_1="potassium_background_current" component_2="intracellular_potassium_concentration"/>
      <map_variables variable_1="i_b_K" variable_2="i_b_K"/>
   </connection>
   <connection>
      <map_components component_1="sodium_background_current" component_2="intracellular_sodium_concentration"/>
      <map_variables variable_1="i_b_Na" variable_2="i_b_Na"/>
   </connection>
   <connection>
      <map_components component_1="CaMKII_factor" component_2="calmodulin"/>
      <map_variables variable_1="Cmdn_Ca" variable_2="Cmdn_Ca"/>
   </connection>
   <connection>
      <map_components component_1="CaMKII_factor" component_2="RyR"/>
      <map_variables variable_1="F_CaMK" variable_2="F_CaMK"/>
   </connection>
   <connection>
      <map_components component_1="CaMKII_factor" component_2="SERCA"/>
      <map_variables variable_1="F_CaMK" variable_2="F_CaMK"/>
   </connection>
   <connection>
      <map_components component_1="RyR" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="j_rel" variable_2="j_rel"/>
   </connection>
   <connection>
      <map_components component_1="RyR" component_2="SR_calcium_concentration"/>
      <map_variables variable_1="Ca_SR" variable_2="Ca_SR"/>
      <map_variables variable_1="j_rel" variable_2="j_rel"/>
   </connection>
   <connection>
      <map_components component_1="SERCA" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="j_up" variable_2="j_up"/>
   </connection>
   <connection>
      <map_components component_1="SERCA" component_2="SR_calcium_concentration"/>
      <map_variables variable_1="Ca_SR" variable_2="Ca_SR"/>
      <map_variables variable_1="j_up" variable_2="j_up"/>
   </connection>
   <connection>
      <map_components component_1="calmodulin" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Cmdn_Ca" variable_2="Cmdn_Ca"/>
      <map_variables variable_1="dCmdn_Ca_dtime" variable_2="dCmdn_Ca_dtime"/>
   </connection>
   <connection>
      <map_components component_1="troponin" component_2="intracellular_calcium_concentration"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="Trpn_Ca" variable_2="Trpn_Ca"/>
      <map_variables variable_1="dTrpn_Ca_dtime" variable_2="dTrpn_Ca_dtime"/>
   </connection>
   <connection>
      <map_components component_1="troponin" component_2="Force"/>
      <map_variables variable_1="Trpn_Ca" variable_2="Trpn_Ca"/>
      <map_variables variable_1="Trpn_tot" variable_2="Trpn_tot"/>
      <map_variables variable_1="alpha_trpn" variable_2="alpha_trpn"/>
      <map_variables variable_1="beta_trpn" variable_2="beta_trpn"/>
      <map_variables variable_1="Force_norm" variable_2="Force_norm"/>
   </connection>

<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"><rdf:Description rdf:nodeID="n1"><maximumStepSize xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">0.1</maximumStepSize><pointDensity xmlns="http://www.cellml.org/metadata/simulation/1.0#nonstandard-" xml:lang="en">1000</pointDensity><endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">1</endingValue></rdf:Description><rdf:Description rdf:about="rdf:#04f8d57c-80dd-4af3-9af0-d2ca845b9aaa"><ORG xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#826ed13e-3a92-49aa-a5be-6cef1d7bb8d6"/></ORG><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#90981391-436d-418e-8852-b859fd0ff88f"/></N><EMAIL xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#168f95b1-c557-4910-86ea-13b1d022084e"/></EMAIL></rdf:Description><rdf:Description rdf:about="rdf:#5b6136a0-50ec-442f-9953-de1bc5a21b05"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Philos Transact A Math Phys Eng Sci</title></rdf:Description><rdf:Description rdf:about="rdf:#71d66a3a-219d-4543-aa4a-b38f63d3f3e0"><cmeta:modifier><rdf:Description rdf:about="rdf:#a90f9cc9-eb2f-4842-9959-e80930fcba04"/></cmeta:modifier><rdf:value xml:lang="en">Curation of units and fixed a couple of errors.</rdf:value><modified xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:about="rdf:#122b7ea1-1a96-4d42-b2e3-ad95d0ca91fb"/></modified></rdf:Description><rdf:Description rdf:about="rdf:#7730ac84-575e-4aa6-88d6-6a2ba642875b"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Peter</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Kohl</Family></rdf:Description><rdf:Description rdf:about="rdf:#122b7ea1-1a96-4d42-b2e3-ad95d0ca91fb"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2008-05-08T04:13:19+12:00</W3CDTF></rdf:Description><rdf:Description rdf:about="rdf:#168f95b1-c557-4910-86ea-13b1d022084e"><rdf:type><rdf:Description rdf:about="http://imc.org/vCard/3.0#internet"/></rdf:type><rdf:value xml:lang="en">penny.noble@dpag.ox.ac.uk</rdf:value></rdf:Description><rdf:Description rdf:nodeID="n2"><boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n3"/></boundIntervals></rdf:Description><rdf:Description rdf:about="rdf:#826ed13e-3a92-49aa-a5be-6cef1d7bb8d6"><Orgunit xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en"></Orgunit><Orgname xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">University of Oxford</Orgname></rdf:Description><rdf:Description rdf:about="rdf:#e4675846-84f5-4f18-85f1-58aacc0ddaf6"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Denis</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Noble</Family></rdf:Description><rdf:Description rdf:nodeID="n3"><rdf:first><rdf:Description rdf:nodeID="n1"/></rdf:first><rdf:rest><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/></rdf:rest></rdf:Description><rdf:Description rdf:about="rdf:#7d6d355d-e937-4463-9470-65ed8e2b346c"><rdf:type><rdf:Description rdf:about="http://www.cellml.org/bqs/1.0#Person"/></rdf:type><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#e4675846-84f5-4f18-85f1-58aacc0ddaf6"/></N></rdf:Description><rdf:Description rdf:about="rdf:#642d0f5b-a4f8-4baa-ad5f-65bffc370e3a"><FN xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Penny Noble</FN></rdf:Description><rdf:Description rdf:about="rdf:#a90f9cc9-eb2f-4842-9959-e80930fcba04"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#d1b29776-8a7a-40db-9f84-af7208ea8496"/></N></rdf:Description><rdf:Description rdf:about="#iribe_kohl_noble_2006"><cmeta:comment><rdf:Description rdf:about="rdf:#9b12eab5-0f1b-4cad-9f5d-3615bb1a4f8c"/></cmeta:comment><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#190a1649-9332-418f-b3b3-0914e6fa941a"/></reference><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n4"/></reference><simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n2"/></simulation></rdf:Description><rdf:Description rdf:about="rdf:#7f5fbf85-77fb-49aa-8827-d6ca69ce7f1b"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2008-05-08T00:00:00+00:00</W3CDTF></rdf:Description><rdf:Description rdf:about="rdf:#eeb779cc-ab04-4a63-aa51-f04ae493c585"><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Seq"/></rdf:type><rdf:_1><rdf:Description rdf:about="rdf:#4f389751-d839-4496-982b-3a284cadd91f"/></rdf:_1><rdf:_2><rdf:Description rdf:about="rdf:#e13cad3d-b6a3-4507-8c89-3fc44eae9c9e"/></rdf:_2><rdf:_3><rdf:Description rdf:about="rdf:#7d6d355d-e937-4463-9470-65ed8e2b346c"/></rdf:_3></rdf:Description><rdf:Description rdf:about="rdf:#90981391-436d-418e-8852-b859fd0ff88f"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Penny</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Noble</Family><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en"></Other></rdf:Description><rdf:Description rdf:about="rdf:#190a1649-9332-418f-b3b3-0914e6fa941a"><Pubmed_id xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">16608699</Pubmed_id><JournalArticle xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#8946cd21-9768-4b1c-957e-0fb961441ea5"/></JournalArticle></rdf:Description><rdf:Description rdf:about="rdf:#9b12eab5-0f1b-4cad-9f5d-3615bb1a4f8c"><rdf:value xml:lang="en"></rdf:value><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#dc560eef-7743-40d9-b33a-4d2bcf92be67"/></creator></rdf:Description><rdf:Description rdf:about="rdf:#cc267f5a-4a46-44ae-bbbc-b74d38f47558"><rdf:value xml:lang="en">We hypothesize that slow inactivation of Ca2+/calmodulin-dependent kinase II (CaMKII) and its modulatory effect on sarcoplasmic reticulum (SR) Ca2+ handling are important for various interval-force (I-F) relations, in particular for the beat interval dependency in transient alternans during the decay of post-extrasystolic potentiation. We have developed a mathematical model of a single cardiomyocyte to integrate various I-F relations, including alternans, by incorporating a conceptual CaMKII kinetics model into the SR Ca2+ handling model. Our model integrates I-F relations, such as the beat interval-dependent twitch force duration, restitution and potentiation, positive staircase phenomenon and alternans. We found that CaMKII affects more or less all I-F relations, and it is a key factor for integration of the various I-F relations in our model. Alternans arises, in the model, out of a steep relation between SR Ca2+ load and release, owing to SR load-dependent changes in the releasability of Ca2+ via the ryanodine receptor. Beat interval-dependent CaMKII activity, owing to its kinetic properties and amplifying effect on SR Ca2+ load dependency of Ca2+ release, replicated the beat interval dependency of alternans, as observed experimentally. Additionally, our model enabled reproduction of the effects of various interventions on alternans, such as the slowing or accelerating of Ca2+ release and/or uptake. We conclude that a slow time-dependent factor, represented in the model by CaMKII, is important for the integration of I-F relations, including alternans, and that our model offers a useful tool for further analysis of the roles of integrative Ca2+ handling in myocardial I-F relations.</rdf:value><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#642d0f5b-a4f8-4baa-ad5f-65bffc370e3a"/></creator></rdf:Description><rdf:Description rdf:about="rdf:#e13cad3d-b6a3-4507-8c89-3fc44eae9c9e"><rdf:type><rdf:Description rdf:about="http://www.cellml.org/bqs/1.0#Person"/></rdf:type><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#7730ac84-575e-4aa6-88d6-6a2ba642875b"/></N></rdf:Description><rdf:Description rdf:about="rdf:#4f389751-d839-4496-982b-3a284cadd91f"><rdf:type><rdf:Description rdf:about="http://www.cellml.org/bqs/1.0#Person"/></rdf:type><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#25c5ba58-66cf-419d-83bd-0e55b3900923"/></N></rdf:Description><rdf:Description rdf:about="rdf:#8946cd21-9768-4b1c-957e-0fb961441ea5"><volume xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">364</volume><issued xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:about="rdf:#90e5579c-142f-48f7-bdac-36b41b2ab739"/></issued><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#eeb779cc-ab04-4a63-aa51-f04ae493c585"/></creator><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Modulatory effect of calmodulin-dependent kinase II (CaMKII) on sarcoplasmic reticulum Ca2+ handling and interval-force relations: a modelling study.</title><first_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">1107</first_page><last_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">1133</last_page><Journal xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#5b6136a0-50ec-442f-9953-de1bc5a21b05"/></Journal></rdf:Description><rdf:Description rdf:about=""><cmeta:comment><rdf:Description rdf:about="rdf:#cc267f5a-4a46-44ae-bbbc-b74d38f47558"/></cmeta:comment><created xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:about="rdf:#7f5fbf85-77fb-49aa-8827-d6ca69ce7f1b"/></created><cmeta:modification><rdf:Description rdf:about="rdf:#71d66a3a-219d-4543-aa4a-b38f63d3f3e0"/></cmeta:modification><publisher xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en"></publisher><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#04f8d57c-80dd-4af3-9af0-d2ca845b9aaa"/></creator></rdf:Description><rdf:Description rdf:about="rdf:#25c5ba58-66cf-419d-83bd-0e55b3900923"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Gentaro</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Iribe</Family></rdf:Description><rdf:Description rdf:about="rdf:#90e5579c-142f-48f7-bdac-36b41b2ab739"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2006-05-00 00:00</W3CDTF></rdf:Description><rdf:Description rdf:about="rdf:#dc560eef-7743-40d9-b33a-4d2bcf92be67"><FN xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en"></FN></rdf:Description><rdf:Description rdf:nodeID="n5"><subject_type xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">keyword</subject_type><rdf:value><rdf:Description rdf:nodeID="n6"/></rdf:value></rdf:Description><rdf:Description rdf:about="rdf:#d1b29776-8a7a-40db-9f84-af7208ea8496"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Penny</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Noble</Family></rdf:Description><rdf:Description rdf:nodeID="n6"><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Bag"/></rdf:type><rdf:_1 xml:lang="en">calcium dynamics</rdf:_1><rdf:_2 xml:lang="en">camkii</rdf:_2><rdf:_3 xml:lang="en">electrophysiology</rdf:_3><rdf:_4 xml:lang="en">excitation-contraction coupling</rdf:_4></rdf:Description><rdf:Description rdf:nodeID="n4"><subject xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n5"/></subject></rdf:Description></rdf:RDF></model>