Location: Saucerman, Brunton, Michailova, Mcculloch, 2003 @ e5422b92bc7c / saucerman_brunton_michailova_mcculloch_2003.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2007-05-23 05:31:15+12:00
Desc:
committing version02 of saucerman_brunton_michailova_mcculloch_2003
Permanent Source URI:
https://models.physiomeproject.org/workspace/saucerman_brunton_michailova_mcculloch_2003/rawfile/e5422b92bc7cca4e36c92f8866d76b74f22662c3/saucerman_brunton_michailova_mcculloch_2003.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  saucerman_model_2003.xml

CREATED :  18th November 2003

LAST MODIFIED : 18th November 2003

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/1/02 CellML Metadata 1.0 Specification.

DESCRIPTION :  This file contains a CellML description of Saucerman et al.'s 2003 mathematical model of beta-adrenergic control of cardiac myocyte contractility. 

CHANGES:  
  
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="saucerman_model_2003" name="saucerman_2003_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modelling Beta-adrenergic Control of Cardiac Myocyte Contractility in Silico</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
        This CellML version of the model can not be run in either PCEnv or COR at the moment due to the presence of circular arguments.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
In cardiac myocytes, the beta-adrenergic signalling network is triggered in response to norepinephrine and epinephrine binding to the Gs-protein coupled beta-adrenergic receptor.  Receptor-ligand binding activates the coupled Gs-protein, which in turn activates adenylate cyclase. ATP is converted into the secondary messenger cyclic AMP (cAMP).  cAMP promotes dissociation of the protein kinase A (PKA) holoenzyme, whose catalytic subunits go on to phosphorylate a wide range of target proteins, two of which are the L-type calcium channel and phospholamban - which both play essential roles in the regulation of calcium dynamics and transport.  In healthy cardiac myocytes, the end result of this signal transduction pathway is to provide coordinated control of contractility, metabolism and gene regulation.  However, altered beta-adrenergic signalling may also play a role in the progression of heart failure.
</para> 

<para>
Molecular components of this signalling pathway have been studied in detail as potential therapeutic targets in heart failure.  However, the situation is complex, with intracellular compartmentation and functional integration of pathways being suggested as playing an important role in the signalling outcome.  Systems level mathematical modelling will help make sense of large quantities of experimental data.  Since the Hodgkin-Huxley Squid Axon Model was published in 1952, many ionic models of cell electrophysiology have been developed.  Recently systems modelling has been identified as a potential method for improving understanding of signaling networks and the response to genetic and pharmaceutical perturbations. 
</para>

<para>
While to date portions of the mammalian beta-adrenergic signalling network have been described by mathematical models of neuron and HEK-293 cells, no analysis has modelled and validated an entire pathway from ligand to effectors such that neuro-hormonal regulation of cell physiology can be predicted.  In the publication described here, Jeffrey Saucerman <emphasis>et al.</emphasis> have developed a systems biology model to investigate the molecular mechanisms underlying the control of the beta-adrenergic signalling network over cardiac myocyte contractility (see figure below).  The model was validated using a wide range of experimental data, and model simulations were used to investigate quantitatively the effects of specific molecular perturbations (with potential for pharmaceutical testing).  By performing systems analysis the effects of molecular perturbations in the beta-adrenergic signalling network may be understood within the context of integrative physiology.  
</para>

<para>
The model of the signal transduction pathway was embedded within an extension of the Luo-Rudy Ventricular Model II (dynamic), 1994, which was modified for the rabbit ventricular myocyte (Puglisi-Bers Rabbit Ventricular Myocyte Model, 2001).  The model was further adapted by inserting the equations for the L-type calcium channel from the Jafri-Rice-Winslow Ventricular Model, 1998, and also the steady-state potassium currents from the Pandit <emphasis>et al.</emphasis> Adult Rat Left Ventricular Myocyte Model, 2001. 
</para>  

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

<para>
<ulink url="http://www.jbc.org/cgi/content/abstract/278/48/47997">Modeling beta-adrenergic control of cardiac myocyte contractility <emphasis>in silico</emphasis>
          </ulink>, Jeffrey J. Saucerman, Laurence L. Brunton, Anushka P. Michailova, and Andrew D. McCulloch, 2003, <ulink url="http://www.jbc.org/">
            <emphasis>Journal of Biological Chemistry</emphasis>
          </ulink>, 48, 47997-48003.  (<ulink url="http://www.jbc.org/cgi/content/full/278/48/47997">Full text (HTML)</ulink> and <ulink url="http://www.jbc.org/cgi/reprint/278/48/47997.pdf">PDF</ulink> versions of the article are available on the <emphasis>Journal of Biological Chemistry</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=12972422&amp;dopt=Abstract">PubMed ID: 12972422</ulink>
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram</title>
    </objectinfo>
    <imagedata fileref="saucerman_mcculloch_2003.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the integrated model components, including the beta-adrenergic network, calcium handling, and the electrophysiology of the rat ventricular myocyte.</caption>
</informalfigure>

<informalfigure float="0" id="fig_pathway_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>pathway diagram</title>
    </objectinfo>
    <imagedata fileref="saucerman_2003.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram highlighting the signalling network component of the model.</caption>
</informalfigure>

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



  
  
  
  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="millimolar">
    <unit units="mole" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>
  
  <units name="microlitre">
    <unit units="litre" prefix="micro"/>
  </units>
  
  <units name="millisecond">
    <unit units="second" prefix="milli"/>
  </units>
  
  <units name="flux">
    <unit units="micromolar"/>
    <unit units="second" exponent="-1"/>
  </units>

  <units name="millimolar_per_second">
    <unit units="millimolar"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="second_order_rate_constant">
    <unit units="micromolar" exponent="-1"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="per_millisecond">
    <unit units="second" prefix="milli" exponent="-1"/>
  </units>

  <units name="per_millivolt">
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   
   <units name="per_millivolt_second">
      <unit units="volt" prefix="milli" exponent="-1"/>
	  <unit units="second" exponent="-1"/>
   </units>
   
   <units name="cm_per_second">
      <unit units="metre" prefix="centi"/>
	  <unit units="second" exponent="-1"/>
   </units>
   
   <units name="cm2">
      <unit units="metre" prefix="centi" exponent="2"/>
   </units>
   
   <units name="milliS_per_microF">
      <unit units="siemens" prefix="milli"/>
	  <unit units="farad" prefix="micro" exponent="-1"/>
   </units>
   
   <units name="microA_per_microF">
      <unit units="ampere" prefix="micro"/>
      <unit units="farad" prefix="micro"/>
   </units>
   
   <units name="millimolar3_per_second">
      <unit units="millimolar" exponent="-3"/>
      <unit units="second" exponent="-1"/>
   </units>
   
   <units name="millimolar4_per_second">
      <unit units="millimolar" exponent="-4"/>
      <unit units="second" exponent="-1"/>
   </units>
   
   <units name="millimolar_4">
      <unit units="millimolar" exponent="-4"/>
   </units>
   
   <units name="joule_per_kilomole_kelvin">
      <unit units="joule"/>
      <unit units="mole" prefix="kilo" exponent="-1"/>
      <unit units="kelvin" exponent="-1"/>
   </units>
   
   <units name="coulomb_per_mole">
      <unit units="coulomb" exponent="-1"/>
      <unit units="mole"/>
   </units>
   
   <units name="pico_litre">
      <unit units="litre" prefix="pico"/>
   </units>
  
  <component name="environment">
    <variable units="second" public_interface="out" name="time"/>
  </component>

<!-- Signal transduction pathway component of the model from the Saucerman et al. 2003 paper -->
  
  <component name="beta_adrenergic_receptor_module"> 
    <variable units="micromolar" public_interface="out" name="LRG"/>
    <variable units="micromolar" public_interface="out" name="RG"/>
    
    <variable units="micromolar" name="LR"/>
    <variable units="micromolar" name="beta1_AR" initial_value="0.01205"/>
    <variable units="micromolar" name="Gs" initial_value="3.829"/>
    <variable units="micromolar" name="beta1_ARact" initial_value="0.01205"/>
    <variable units="micromolar" name="beta1_AR_S464" initial_value="0.0"/>
    <variable units="micromolar" name="beta1_AR_S301" initial_value="1.154E-3"/>
    <variable units="micromolar" name="L" initial_value="10.0"/>
    <variable units="micromolar" name="Ltotmax" initial_value="1.0"/>
    <variable units="micromolar" name="KL" initial_value="0.285"/>
    <variable units="micromolar" name="KR" initial_value="0.062"/>
    <variable units="micromolar" name="KC" initial_value="33.0"/>
    <variable units="micromolar" name="Gstot" initial_value="3.83"/>
    <variable units="first_order_rate_constant" name="k_betaARK_plus" initial_value="1.1E-3"/>
    <variable units="first_order_rate_constant" name="k_betaARK_minus" initial_value="2.2E-3"/>
    <variable units="second_order_rate_constant" name="k_PKA_plus" initial_value="3.6E-3"/>
    <variable units="first_order_rate_constant" name="k_PKA_minus" initial_value="2.2E-3"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="micromolar" public_interface="in" name="Gs_beta_gamma"/>
    <variable units="micromolar" public_interface="in" name="PKACI"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="1"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> L </ci> 
        </apply> 
        <apply>
          <minus/>
          <ci> Ltotmax </ci>
          <apply>
            <plus/>
            <ci> LR </ci>
            <ci> LRG </ci>
            <ci> L </ci> 
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>LR</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> L </ci>
            <ci> beta1_AR </ci>
          </apply>
          <ci> KL </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>LRG</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> L </ci>
            <ci> beta1_AR </ci>
            <ci> Gs </ci>
          </apply>
          <apply>
            <times/>
            <ci> KL </ci>
            <ci> KR </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>RG</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> beta1_AR </ci>
            <ci> Gs </ci>
          </apply>
          <ci> KC </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> Gs </ci> 
        </apply> 
        <apply>
          <minus/>
          <ci> Gstot </ci>
          <apply>
            <plus/>
            <ci> RG </ci>
            <ci> LRG </ci>
            <ci> Gs </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> beta1_AR </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> beta1_ARact </ci>
          <apply>
            <plus/>
            <ci> LR </ci>
            <ci> LRG </ci>
            <ci> RG </ci>
            <ci> beta1_AR </ci> 
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>beta1_ARact</ci> 
        </apply> 
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k_betaARK_minus </ci>
              <ci> beta1_AR_S464 </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_betaARK_plus </ci>
              <apply>
                <plus/>
                <ci> LR </ci>
                <ci> LRG </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k_PKA_minus </ci>
              <ci> beta1_AR_S301 </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_PKA_plus </ci>
              <ci> PKACI </ci>
              <ci> beta1_ARact </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>beta1_AR_S464</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_betaARK_plus </ci>
            <apply>
              <plus/>
              <ci> LR </ci>
              <ci> LRG </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_betaARK_minus </ci>
            <ci> beta1_AR_S464 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>beta1_AR_S301</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_PKA_plus </ci>
            <ci> PKACI </ci>
            <ci> beta1_ARact </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_PKA_minus </ci>
            <ci> beta1_AR_S301 </ci>
          </apply>
        </apply>
      </apply>        
    </math>
  </component>
  
  <component name="Gs_activation_module"> 
    <variable units="micromolar" public_interface="out" name="Gs_beta_gamma" initial_value="0.02569"/>
    <variable units="micromolar" public_interface="out" name="Gs_alpha_GTPtot" initial_value="0.02505"/>
   
    <variable units="micromolar" name="Gs_alpha_GDP" initial_value="6.44E-4"/>
    <variable units="first_order_rate_constant" name="k_gact" initial_value="16.0"/>
    <variable units="first_order_rate_constant" name="k_hyd" initial_value="0.8"/>
    <variable units="second_order_rate_constant" name="k_reassoc" initial_value="1.21E3"/>
    
    <variable units="micromolar" public_interface="in" name="RG"/>
    <variable units="micromolar" public_interface="in" name="LRG"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="2">
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Gs_alpha_GTPtot</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_gact </ci>
            <apply>
              <plus/>
              <ci> RG </ci>
              <ci> LRG </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_hyd </ci>
            <ci> Gs_alpha_GTPtot </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Gs_beta_gamma</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_gact </ci>
            <apply>
              <plus/>
              <ci> RG </ci>
              <ci> LRG </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_reassoc </ci>
            <ci> Gs_alpha_GDP </ci>
            <ci> Gs_beta_gamma </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Gs_alpha_GDP</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_hyd </ci>
            <ci> Gs_alpha_GTPtot </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_reassoc </ci>
            <ci> Gs_alpha_GDP </ci>
            <ci> Gs_beta_gamma </ci>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>
  
  <component name="cyclic_AMP_metabolism_module">
    <variable units="micromolar" public_interface="out" name="cAMPtot" initial_value="0.8453"/>
    
    <variable units="micromolar" name="AC" initial_value="0.04706"/>
    <variable units="micromolar" name="Gs_alpha_GTP" initial_value="0.02241"/>
    <variable units="micromolar" name="Gs_alpha_GTP_AC"/>
    <variable units="micromolar" name="PDE" initial_value="0.0389"/>
    <!-- <variable units="micromolar" name="PDE_ACT"/> -->
    <variable units="micromolar" name="ACtot" initial_value="49.7E-3"/>
    <variable units="micromolar" name="ATP" initial_value="5.0E3"/>
    <variable units="micromolar" name="PDEtot" initial_value="38.9E-3"/>
    <variable units="micromolar" name="PDEinhib"/>
    <variable units="micromolar" name="IBMXtot" initial_value="0.0"/>
    <variable units="micromolar" name="IBMX" initial_value="0.0"/>
    <variable units="micromolar" name="fsktot" initial_value="0.0"/>
    <variable units="micromolar" name="fsk" initial_value="0.0"/>
    <variable units="micromolar" name="fsk_AC"/>
    <variable units="first_order_rate_constant" name="kAC_basal" initial_value="0.2"/>
    <variable units="first_order_rate_constant" name="kAC_fsk" initial_value="7.3"/>
    <variable units="first_order_rate_constant" name="k_PDE" initial_value="5.0"/>
    <variable units="first_order_rate_constant" name="kAC_Gs_alpha" initial_value="8.5"/>
    <variable units="micromolar" name="Km_basal" initial_value="1.03E3"/>
    <variable units="micromolar" name="Km_PDE" initial_value="1.3"/>
    <variable units="micromolar" name="Km_fsk" initial_value="860.0"/>
    <variable units="micromolar" name="K_fsk" initial_value="44.0"/>
    <variable units="micromolar" name="Km_Gs_alpha_GTP" initial_value="0.4"/>
    <variable units="micromolar" name="K_Gs_alpha" initial_value="0.4"/>
    <variable units="micromolar" name="K_IBMX" initial_value="30.0"/>
    
    <variable units="micromolar" public_interface="in" name="Gs_alpha_GTPtot"/>
    <variable units="second" public_interface="in" name="time"/>
    <variable units="micromolar" public_interface="in" name="cAMP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="3"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> Gs_alpha_GTP </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> Gs_alpha_GTPtot </ci>
          <apply>
            <plus/>
            <ci> Gs_alpha_GTP_AC </ci>
            <ci> Gs_alpha_GTP </ci> 
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> fsk </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> fsktot </ci>
          <apply>
            <plus/>
            <ci> fsk_AC </ci>
            <ci> fsk </ci> 
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> AC </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> ACtot </ci>
          <apply>
            <plus/>
            <ci> Gs_alpha_GTP_AC </ci>
            <ci> AC </ci> 
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> PDE </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> PDEtot </ci>
          <apply>
            <plus/>
            <ci> PDEinhib </ci>
            <ci> PDE </ci> 
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> IBMX </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> IBMXtot </ci>
          <apply>
            <plus/>
            <ci> PDEinhib </ci>
            <ci> IBMX </ci> 
          </apply>
        </apply>
      </apply>
<!--
      <apply><eq/> 
        <apply><diff/> 
          <bvar><ci> time </ci></bvar> 
          <ci> PDE_ACT </ci> 
        </apply>
        <apply><divide/>
          <apply><times/>
            <ci> k_PDE </ci>
            <ci> PDE </ci>
            <ci> cAMP </ci>
          </apply>
          <apply><plus/>
            <ci> Km_PDE </ci>
            <ci> cAMP </ci> 
          </apply>
        </apply>
      </apply>
	  
	  -->
      
      <apply>
        <eq/> 
        <ci>Gs_alpha_GTP_AC</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> Gs_alpha_GTP </ci>
            <ci> AC </ci>
          </apply>
          <ci> K_Gs_alpha </ci>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci>fsk_AC</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> fsk </ci>
            <ci> AC </ci>
          </apply>
          <ci> K_fsk </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>PDEinhib</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> PDE </ci>
            <ci> IBMX </ci>
          </apply>
          <ci> K_IBMX </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>cAMPtot</ci> 
        </apply>
        <apply>
          <minus/> 
          <apply>
            <plus/>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> kAC_basal </ci>
                <ci> AC </ci>
                <ci> ATP </ci>
              </apply>
              <apply>
                <plus/>
                <ci> Km_basal </ci>
                <ci> ATP </ci>
              </apply>
            </apply>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> kAC_Gs_alpha </ci>
                <ci> Gs_alpha_GTP_AC </ci>
                <ci> ATP </ci>
              </apply>
              <apply>
                <plus/>
                <ci> Km_Gs_alpha_GTP </ci>
                <ci> ATP </ci>
              </apply>
            </apply>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> kAC_fsk </ci>
                <ci> fsk_AC </ci>
                <ci> ATP </ci>
              </apply>
              <apply>
                <plus/>
                <ci> Km_fsk </ci>
                <ci> ATP </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> k_PDE </ci>
              <ci> PDE </ci>
              <ci> cAMP </ci>
            </apply>
            <apply>
              <plus/>
              <ci> Km_PDE </ci>
              <ci> cAMP </ci>
            </apply>
          </apply>
        </apply>
      </apply>        
    </math>
  </component>
  
  <component name="PKA_activation_module">
    <variable units="micromolar" public_interface="out" name="PKACI"/>
    <variable units="micromolar" public_interface="out" name="PKACII"/>
    <variable units="micromolar" public_interface="out" name="cAMP"/>
    
    <variable units="micromolar" name="ARCI"/>
    <variable units="micromolar" name="A2RCI"/>
    <variable units="micromolar" name="A2RCII"/>
    <variable units="micromolar" name="A2RI"/>
    <variable units="micromolar" name="A2RII"/>
    <variable units="micromolar" name="ARCII"/>
    <variable units="micromolar" name="RCI"/>
    <variable units="micromolar" name="PKACI_PKI"/>
    <variable units="micromolar" name="PKAtemp"/>
    <variable units="micromolar" name="PKI"/>
    <variable units="micromolar" name="RCII"/>
    <variable units="micromolar" name="PKACII_PKI"/>
    <variable units="micromolar" name="Ki_pki" initial_value="0.2E-3"/>
    <variable units="micromolar" name="PKAItot" initial_value="0.59"/>
    <variable units="micromolar" name="PKAIItot" initial_value="0.025"/>
    <variable units="micromolar" name="PKItot" initial_value="0.18"/>
    <variable units="micromolar" name="KA" initial_value="9.14"/>
    <variable units="micromolar" name="KB" initial_value="1.64"/>
    <variable units="micromolar" name="KD" initial_value="4.375"/>
    <variable units="micromolar" name="KPKI" initial_value="2.0E-4"/>
    
    <variable units="micromolar" public_interface="in" name="cAMPtot"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="4"> 
      <apply>
        <eq/> 
        <ci>cAMP</ci> 
        <apply>
          <minus/>
          <ci> cAMPtot </ci>
          <apply>
            <plus/>
            <apply>
              <plus/>
              <ci> ARCI </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> A2RCI </ci>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> A2RI </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <ci> ARCII </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> A2RCII </ci>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> A2RII </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> PKAtemp </ci> 
        <apply>
          <plus/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> KA </ci>
              <ci> KB </ci>
            </apply>
            <ci> KD </ci>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> KA </ci>
              <ci> cAMP </ci>
            </apply>
            <ci> KD </ci>
          </apply>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> cAMP </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <ci> KD </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PKACI </ci> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> PKAItot </ci>
            <apply>
              <power/>
              <ci> cAMP </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> PKAtemp </ci>
                <ci> PKACI </ci> 
              </apply>
              <apply>
                <power/>
                <ci> cAMP </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> PKACI </ci>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> PKI </ci>
                <ci> Ki_pki </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> PKACII </ci> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> PKAIItot </ci>
            <apply>
              <power/>
              <ci> cAMP </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> PKAtemp </ci>
                <ci> PKACII </ci> 
              </apply>
              <apply>
                <power/>
                <ci> cAMP </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> PKACII </ci>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> PKI </ci>
                <ci> Ki_pki </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>RCI</ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> KA </ci>
              <ci> KB </ci>
            </apply>
            <apply>
              <power/>
              <ci> cAMP </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> PKACI </ci>
              <ci> KD </ci>
            </apply>
            <ci> PKACI </ci>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> PKI </ci>
                <ci> Ki_pki </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>ARCI</ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <ci> KA </ci>
            <ci> cAMP </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> PKACI </ci>
              <ci> KD </ci>
            </apply>
            <ci> PKACI </ci>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> PKI </ci>
                <ci> Ki_pki </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci>ARCII</ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <ci> KA </ci>
            <ci> cAMP </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> PKACII </ci>
              <ci> KD </ci>
            </apply>
            <ci> PKACII </ci>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> PKI </ci>
                <ci> Ki_pki </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>A2RCI</ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <ci> PKACI </ci>
            <ci> KD </ci>
          </apply>
          <ci> PKACI </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> PKI </ci>
              <ci> Ki_pki </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>A2RI</ci> 
        <apply>
          <times/>
          <ci> PKACI </ci> 
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> PKI </ci>
              <ci> Ki_pki </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>RCII</ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> KA </ci>
              <ci> KB </ci>
            </apply>
            <apply>
              <power/>
              <ci> cAMP </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> PKACII </ci>
              <ci> KD </ci>
            </apply>
            <ci> PKACII </ci>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> PKI </ci>
                <ci> Ki_pki </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>A2RCII</ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <ci> PKACII </ci>
            <ci> KD </ci>
          </apply>
          <ci> PKACII </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> PKI </ci>
              <ci> Ki_pki </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>A2RII</ci> 
        <apply>
          <times/>
          <ci> PKACII </ci> 
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> PKI </ci>
              <ci> Ki_pki </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>PKI</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> KPKI </ci>
            <ci> PKItot </ci>
          </apply>
          <apply>
            <plus/>
            <ci> KPKI </ci>
            <ci> PKACI </ci> 
            <ci> PKACII </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>PKACI_PKI</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> PKACI </ci>
            <ci> PKItot </ci>
          </apply>
          <apply>
            <plus/>
            <ci> KPKI </ci>
            <ci> PKACI </ci> 
            <ci> PKACII </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>PKACII_PKI</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> PKACI </ci>
            <ci> PKItot </ci>
          </apply>
          <apply>
            <plus/>
            <ci> KPKI </ci>
            <ci> PKACI </ci> 
            <ci> PKACII </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="phospholamban_regulation_module">
    <variable units="dimensionless" public_interface="out" name="fracPLBp"/>
    <variable units="dimensionless" public_interface="out" name="fracPLB"/>
    <variable units="dimensionless" public_interface="out" name="fracPLBo" initial_value="0.9613"/>
 
    <variable units="micromolar" name="PLBp" initial_value="4.105"/>
    <variable units="micromolar" name="PLB"/>
    <variable units="micromolar" name="Inhib1ptot" initial_value="0.0526"/>
    <variable units="micromolar" name="Inhib1" initial_value="0.2474"/>
    <variable units="micromolar" name="Inhib1p" initial_value="6.3E-5"/>
    <variable units="micromolar" name="PP1_Inhib1p"/>
    <variable units="micromolar" name="PP1" initial_value="0.838"/>
    <variable units="micromolar" name="PP1tot" initial_value="0.89"/>
    <variable units="micromolar" name="PLBtot" initial_value="106.0"/>
    <variable units="micromolar" name="Inhib1tot" initial_value="0.3"/>
    <variable units="dimensionless" name="epsilon" initial_value="10.0"/>
    <variable units="first_order_rate_constant" name="kPKA_PLB" initial_value="54.0"/>
    <variable units="micromolar" name="KmPKA_PLB" initial_value="21.0"/>
    <variable units="first_order_rate_constant" name="kPKA_Inhib1" initial_value="60.0"/>
    <variable units="first_order_rate_constant" name="kPP1_PLB" initial_value="8.5"/>
    <variable units="micromolar" name="KmPP1_PLB" initial_value="7.0"/>
    <variable units="micromolar" name="KmPKA_Inhib1" initial_value="1.0"/>
    <variable units="flux" name="VmaxPP2A_Inhib1" initial_value="14.0"/>
    <variable units="micromolar" name="KmPP2A_Inhib1" initial_value="1.0"/>
    <variable units="micromolar" name="KInhib1" initial_value="1.0E-3"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="micromolar" public_interface="in" name="PKACI"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="5"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>PP1</ci> 
        </apply>  
        <apply>
          <minus/>
          <ci> PP1tot </ci>
          <apply>
            <plus/>
            <ci> PP1_Inhib1p </ci>
            <ci> PP1 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci>PP1_Inhib1p</ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> PP1 </ci>
            <ci> Inhib1p </ci>
          </apply>
          <ci> KInhib1 </ci>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Inhib1p</ci> 
        </apply>  
        <apply>
          <minus/>
          <ci> Inhib1ptot </ci>
          <apply>
            <plus/>
            <ci> PP1_Inhib1p </ci>
            <ci> Inhib1p </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>PLBp</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> kPKA_PLB </ci>
              <ci> PKACI </ci>
              <ci> PLB </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_PLB </ci>
              <ci> PLB </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> kPP1_PLB </ci>
              <ci> PP1 </ci>
              <ci> PLBp </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP1_PLB </ci>
              <ci> PLBp </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <ci>PLB</ci> 
        <apply>
          <minus/>
          <ci> PLBtot </ci>
          <ci> PLBp </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Inhib1ptot</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> kPKA_Inhib1 </ci>
              <ci> PKACI </ci>
              <ci> Inhib1ptot </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_Inhib1 </ci>
              <ci> Inhib1ptot </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> VmaxPP2A_Inhib1 </ci>
              <ci> Inhib1 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_Inhib1 </ci>
              <ci> Inhib1 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Inhib1</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> VmaxPP2A_Inhib1 </ci>
              <ci> Inhib1 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_Inhib1 </ci>
              <ci> Inhib1 </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> kPKA_Inhib1 </ci>
              <ci> PKACI </ci>
              <ci> Inhib1ptot </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_Inhib1 </ci>
              <ci> Inhib1ptot </ci>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> fracPLBp </ci> 
        <apply>
          <divide/>
          <ci> PLBp </ci>
          <ci> PLBtot </ci>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> fracPLB </ci> 
        <apply>
          <divide/>
          <ci> PLB </ci>
          <ci> PLBtot </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="L_type_Ca_channel_regulation_module">
    <variable units="dimensionless" public_interface="out" name="fracLCCap"/>
    <variable units="dimensionless" public_interface="out" name="fracLCCbp"/>
    <variable units="dimensionless" public_interface="out" name="fracLCCapo" initial_value="0.2041"/>
    <variable units="dimensionless" public_interface="out" name="fracLCCbpo" initial_value="0.2336"/>
 
    <variable units="micromolar" name="LCC_beta2_p" initial_value="1.27575E-4"/>
    <variable units="micromolar" name="LCC_beta2" initial_value="1.46025E-3"/>
    <!--
	<variable units="micromolar" name="P1"/>
    <variable units="micromolar" name="P2"/>
    <variable units="micromolar" name="P3"/>
    <variable units="micromolar" name="P4"/>
    <variable units="micromolar" name="P5"/>
    <variable units="micromolar" name="P6"/>
    <variable units="micromolar" name="P7"/>
    <variable units="micromolar" name="P8"/>
    <variable units="micromolar" name="P9"/>
    <variable units="micromolar" name="P10"/>
    <variable units="micromolar" name="P11"/>
    <variable units="micromolar" name="P12"/>
    <variable units="micromolar" name="P13"/>
    <variable units="micromolar" name="P14"/>
    <variable units="micromolar" name="P15"/>
    <variable units="micromolar" name="P16"/>
    <variable units="micromolar" name="P17"/>
    <variable units="micromolar" name="P18"/>
    <variable units="micromolar" name="P19"/>
    <variable units="micromolar" name="P20"/>
    <variable units="micromolar" name="P21"/>
    <variable units="micromolar" name="P22"/>
    <variable units="micromolar" name="P23"/>
    <variable units="micromolar" name="P24"/>
    <variable units="micromolar" name="PHOSPH1"/>
    <variable units="micromolar" name="PHOSPH2"/>
    <variable units="micromolar" name="PHOSPH3"/>
    <variable units="micromolar" name="PHOSPH4"/>
    <variable units="micromolar" name="PHOSPH5"/>
    <variable units="micromolar" name="PHOSPH6"/>
    <variable units="micromolar" name="PHOSPH7"/>
    <variable units="micromolar" name="PHOSPH8"/>
    <variable units="micromolar" name="PHOSPH9"/>
    <variable units="micromolar" name="PHOSPH10"/>
    <variable units="micromolar" name="PHOSPH11"/>
    <variable units="micromolar" name="PHOSPH12"/>
	-->
    <variable units="micromolar" name="LCCtot" initial_value="25.0E-3"/>
    <variable units="micromolar" name="PP1_LCC" initial_value="25.0E-3"/>
    <variable units="micromolar" name="PP2A_LCC" initial_value="25E-3"/>
    <variable units="dimensionless" name="epsilon" initial_value="10.0"/>
    <variable units="first_order_rate_constant" name="kPKA_LCC" initial_value="54.0"/>
    <variable units="micromolar" name="KmPKA_LCC" initial_value="21.0"/>
    <variable units="micromolar" name="KmPP1_LCC" initial_value="3.0"/>
    <variable units="first_order_rate_constant" name="kPP1_LCC" initial_value="8.52"/>
    <variable units="first_order_rate_constant" name="kPP2A_LCC" initial_value="10.1"/>
    <variable units="micromolar" name="KmPP2A_LCC" initial_value="3.0"/>
    <variable units="first_order_rate_constant" name="omega_LCC" initial_value="10.0"/>
    <variable units="first_order_rate_constant" name="gamma0_LCC" initial_value="2.8E3"/>
    <variable units="first_order_rate_constant" name="fn_LCC" initial_value="200.0"/>
    <variable units="first_order_rate_constant" name="gn_LCC" initial_value="2.0E3"/>
    <variable units="first_order_rate_constant" name="fc_LCC" initial_value="5.0"/>
    <variable units="first_order_rate_constant" name="gc_LCC" initial_value="7.0E3"/>
    <variable units="first_order_rate_constant" name="fp_LCC" initial_value="800.0"/>
    <variable units="first_order_rate_constant" name="gp_LCC" initial_value="2.0E3"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="micromolar" public_interface="in" name="PKACII"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="6"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>LCC_beta2_p</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <apply>
                <minus/>
                <ci> LCCtot </ci>
                <ci> LCC_beta2_p </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <apply>
                  <minus/>
                  <ci> LCCtot </ci>
                  <ci> LCC_beta2_p </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> LCC_beta2_p </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCC_beta2_p </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>LCC_beta2</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <apply>
                <minus/>
                <ci> LCCtot </ci>
                <ci> LCC_beta2 </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <apply>
                  <minus/>
                  <ci> LCCtot </ci>
                  <ci> LCC_beta2 </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP1_LCC </ci>
              <ci> PP1_LCC </ci>
              <ci> LCC_beta2 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP1_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCC_beta2 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> fracLCCbp </ci> 
        <apply>
          <divide/>
          <ci> LCC_beta2 </ci>
          <ci> LCCtot </ci>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> fracLCCap </ci> 
        <apply>
          <divide/>
          <ci> LCC_beta2_p </ci>
          <ci> LCCtot </ci>
        </apply>
      </apply>

	  
<!-- Mode Normal -->
      <!--
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P1</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> beta_LCC </ci>
              <ci> P2 </ci>
            </apply>
            <apply>
              <times/>
              <ci> omega_LCC </ci>
              <ci> P7 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                  <ci> alpha_LCC </ci>
                </apply>
                <ci> gamma_LCC </ci>
              </apply>
              <ci> P1 </ci>
            </apply>
            <ci> PHOSPH1 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P2</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P1 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> beta_LCC </ci>
              <ci> P3 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> P8 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <ci> beta_LCC </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                  <ci> alpha_LCC </ci>
                </apply>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> gamma_LCC </ci>
                </apply>
              </apply>
              <ci> P2 </ci>
            </apply>
            <ci> PHOSPH2 </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P3</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P2 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> beta_LCC </ci>
              <ci> P4 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
              <ci> P9 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> beta_LCC </ci>
                </apply>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> alpha_LCC </ci>
                </apply>
                <apply>
                  <times/>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                  </apply>
                  <ci> gamma_LCC </ci>
                </apply>
              </apply>
              <ci> P3 </ci>
            </apply>
            <ci> PHOSPH3 </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P4</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P3 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <ci> beta_LCC </ci>
              <ci> P5 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
              </apply>
              <ci> P10 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                  <ci> beta_LCC </ci>
                </apply>
                <ci> alpha_LCC </ci>
                <apply>
                  <times/>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <cn cellml:units="dimensionless"> 3.0 </cn>
                  </apply>
                  <ci> gamma_LCC </ci>
                </apply>
              </apply>
              <ci> P4 </ci>
            </apply>
            <ci> PHOSPH4 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P5</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> alpha_LCC </ci>
              <ci> P4 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gn_LCC </ci>
              <ci> P6 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                </apply>
              </apply>
              <ci> P11 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                  <ci> beta_LCC </ci>
                </apply>
                <ci> fn_LCC </ci>
                <apply>
                  <times/>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <cn cellml:units="dimensionless"> 4.0 </cn>
                  </apply>
                  <ci> gamma_LCC </ci>
                </apply>
              </apply>
              <ci> P5 </ci>
            </apply>
            <ci> PHOSPH5 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P6</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> fc_LCC </ci>
            <ci> P5 </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> gc_LCC </ci>
              <ci> P6 </ci>
            </apply>
            <ci> PHOSPH6 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P7</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> beta_LCC </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> P8 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gamma_LCC </ci>
              <ci> P1 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> alpha_LCC </ci>
                </apply>
                <ci> omega_LCC </ci>
              </apply>
              <ci> P7 </ci>
            </apply>
            <ci> PHOSPH7 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P8</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P7 </ci>
            </apply>
            <apply>
              <times/>
              <ci> beta_LCC </ci>
              <ci> P9 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> gamma_LCC </ci>
              <ci> P2 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <divide/>
                  <ci> beta_LCC </ci>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> alpha_LCC </ci>
                </apply>
                <apply>
                  <divide/>
                  <ci> omega_LCC </ci>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
              <ci> P8 </ci>
            </apply>
            <ci> PHOSPH8 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P9</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P8 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                  <ci> beta_LCC </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> P10 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <power/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> gamma_LCC </ci>
              <ci> P3 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <ci> beta_LCC </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> alpha_LCC </ci>
                </apply>
                <apply>
                  <divide/>
                  <ci> omega_LCC </ci>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                  </apply>
                </apply>
              </apply>
              <ci> P9 </ci>
            </apply>
            <ci> PHOSPH9 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P10</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P9 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> beta_LCC </ci>
              <ci> P11 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <power/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <ci> gamma_LCC </ci>
              <ci> P4 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 3.0 </cn>
                    <ci> beta_LCC </ci>
                  </apply>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> alpha_LCC </ci>
                </apply>
                <apply>
                  <divide/>
                  <ci> omega_LCC </ci>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <cn cellml:units="dimensionless"> 3.0 </cn>
                  </apply>
                </apply>
              </apply>
              <ci> P10 </ci>
            </apply>
            <ci> PHOSPH10 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P11</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P10 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gc_LCC </ci>
              <ci> P12 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <power/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <ci> gamma_LCC </ci>
              <ci> P5 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 4.0 </cn>
                    <ci> beta_LCC </ci>
                  </apply>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
                <ci> fc_LCC </ci>
                <apply>
                  <divide/>
                  <ci> omega_LCC </ci>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <cn cellml:units="dimensionless"> 4.0 </cn>
                  </apply>
                </apply>
              </apply>
              <ci> P11 </ci>
            </apply>
            <ci> PHOSPH11 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P12</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> fc_LCC </ci>
            <ci> P11 </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> gc_LCC </ci>
              <ci> P12 </ci>
            </apply>
            <ci> PHOSPH12 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P13</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> beta_LCC </ci>
              <ci> P14 </ci>
            </apply>
            <apply>
              <times/>
              <ci> omega_LCC </ci>
              <ci> P19 </ci>
            </apply>
            <ci> PHOSPH1 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 4.0 </cn>
                <ci> alpha_LCC </ci>
              </apply>
              <ci> gamma_LCC </ci>
            </apply>
            <ci> P13 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P14</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P13 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> beta_LCC </ci>
              <ci> P15 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> P20 </ci>
            </apply>
            <ci> PHOSPH2 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <ci> beta_LCC </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 3.0 </cn>
                <ci> alpha_LCC </ci>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> gamma_LCC </ci>
              </apply>
            </apply>
            <ci> P14 </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P15</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P14 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> beta_LCC </ci>
              <ci> P16 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
              <ci> P21 </ci>
            </apply>
            <ci> PHOSPH3 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> beta_LCC </ci>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> alpha_LCC </ci>
              </apply>
              <apply>
                <times/>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
                <ci> gamma_LCC </ci>
              </apply>
            </apply>
            <ci> P15 </ci>
          </apply>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P16</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P15 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <ci> beta_LCC </ci>
              <ci> P17 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
              </apply>
              <ci> P22 </ci>
            </apply>
            <ci> PHOSPH4 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 3.0 </cn>
                <ci> beta_LCC </ci>
              </apply>
              <ci> alpha_LCC </ci>
              <apply>
                <times/>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
                <ci> gamma_LCC </ci>
              </apply>
            </apply>
            <ci> P16 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P17</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> alpha_LCC </ci>
              <ci> P16 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gn_LCC </ci>
              <ci> P18 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                </apply>
              </apply>
              <ci> P23 </ci>
            </apply>
            <ci> PHOSPH5 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 4.0 </cn>
                <ci> beta_LCC </ci>
              </apply>
              <ci> fn_LCC </ci>
              <apply>
                <times/>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                </apply>
                <ci> gamma_LCC </ci>
              </apply>
            </apply>
            <ci> P17 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P18</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> fc_LCC </ci>
              <ci> P17 </ci>
            </apply>
            <ci> PHOSPH6 </ci>
          </apply>
          <apply>
            <times/>
            <ci> gc_LCC </ci>
            <ci> P18 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P19</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> beta_LCC </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <ci> P20 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gamma_LCC </ci>
              <ci> P13 </ci>
            </apply>
            <ci> PHOSPH7 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 4.0 </cn>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> alpha_LCC </ci>
              </apply>
              <ci> omega_LCC </ci>
            </apply>
            <ci> P19 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P20</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 4.0 </cn>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P19 </ci>
            </apply>
            <apply>
              <times/>
              <ci> beta_LCC </ci>
              <ci> P21 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> gamma_LCC </ci>
              <ci> P14 </ci>
            </apply>
            <ci> PHOSPH8 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <divide/>
                <ci> beta_LCC </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 3.0 </cn>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> alpha_LCC </ci>
              </apply>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <ci> P20 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P21</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P18 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                  <ci> beta_LCC </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> P20 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <power/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> gamma_LCC </ci>
              <ci> P15 </ci>
            </apply>
            <ci> PHOSPH9 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <ci> beta_LCC </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> alpha_LCC </ci>
              </apply>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> P21 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P22</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P17 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> beta_LCC </ci>
              <ci> P19 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <power/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <ci> gamma_LCC </ci>
              <ci> P16 </ci>
            </apply>
            <ci> PHOSPH10 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                  <ci> beta_LCC </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> alpha_LCC </ci>
              </apply>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> P22 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P23</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> alpha_LCC </ci>
              <ci> P16 </ci>
            </apply>
            <apply>
              <times/>
              <ci> gc_LCC </ci>
              <ci> P24 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <power/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <ci> gamma_LCC </ci>
              <ci> P17 </ci>
            </apply>
            <ci> PHOSPH11 </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                  <ci> beta_LCC </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> fc_LCC </ci>
              <apply>
                <divide/>
                <ci> omega_LCC </ci>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                </apply>
              </apply>
            </apply>
            <ci> P23 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>P24</ci> 
        </apply> 
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> fc_LCC </ci>
              <ci> P23 </ci>
            </apply>
            <ci> PHOSPH12 </ci>
          </apply>
          <apply>
            <times/>
            <ci> gc_LCC </ci>
            <ci> P24 </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH1 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P1 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P1 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P13 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P13 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH2 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P2 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P2 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P14 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P14 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH3 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P3 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P3 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P15 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P15 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH4 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P4 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P4 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P16 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P16 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH5 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P5 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P5 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P17 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P17 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH6 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P6 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P6 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P18 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P18 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH7 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P7 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P7 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P19 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P19 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH8 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P8 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P8 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P20 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P20 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH9 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P9 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P9 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P21 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P21 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH10 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P10 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P10 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P22 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P22 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH11 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P11 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P11 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P23 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P23 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> PHOSPH12 </ci>
        <apply>
          <minus/> 
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPKA_LCC </ci>
              <ci> PKACII </ci>
              <ci> P12 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPKA_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P12 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <ci> kPP2A_LCC </ci>
              <ci> PP2A_LCC </ci>
              <ci> P24 </ci>
            </apply>
            <apply>
              <plus/>
              <ci> KmPP2A_LCC </ci>
              <apply>
                <times/>
                <ci> epsilon </ci>
                <ci> LCCtot </ci>
                <ci> P24 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
	  
	  -->
     
    </math>
  </component>

<!-- EP component of the model from the Saucerman et al. 2003 paper -->

   <component name="membrane">
      <variable units="millivolt" public_interface="out" name="V" initial_value="-85.66"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="out" name="R" initial_value="8314"/>
      <variable units="kelvin" public_interface="out" name="T" initial_value="310"/>
      <variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96485"/>
      
      <variable units="second" public_interface="in" name="time"/>
      <variable units="microA_per_microF" public_interface="in" name="I_Na"/>
      <variable units="microA_per_microF" public_interface="in" name="I_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="I_CaK"/>
      <variable units="microA_per_microF" public_interface="in" name="I_NaCa"/>
      <variable units="microA_per_microF" public_interface="in" name="I_Ki"/>
      <variable units="microA_per_microF" public_interface="in" name="I_Kp"/>
      <variable units="microA_per_microF" public_interface="in" name="I_p_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="I_Na_b"/>
      <variable units="microA_per_microF" public_interface="in" name="I_Ca_b"/>
      <variable units="microA_per_microF" public_interface="in" name="I_NaK"/>
      <variable units="microA_per_microF" public_interface="in" name="I_to"/>
      <variable units="microA_per_microF" public_interface="in" name="I_ss"/>
      
      <variable units="microA_per_microF" name="I_app"/>
      <variable units="millivolt" name="R_clamp" initial_value="0.2"/>
      <variable units="millivolt" name="V_clamp"/>
      <variable units="millivolt" name="V_test" initial_value="-10.0"/>
      <variable units="millivolt" name="V_hold" initial_value="-40.0"/>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="7">
        <apply>
        <eq/>
          <ci> I_app </ci>
          <apply>
          <divide/>
            <apply>
            <minus/>
              <ci> V_clamp </ci>
              <ci> V </ci>
            </apply>
            <ci> R_clamp </ci>
          </apply>
        </apply>
      
        <apply>
        <eq/>
          <ci> V_clamp </ci>
          <piecewise>
            <piece>
              <ci> V_test </ci>
              <apply>
              <and/>
                <apply>
                <gt/>
                  <ci>time</ci>
                  <cn cellml:units="millisecond">59.1</cn>
                </apply>
                <apply>
                <leq/>
                  <ci>time</ci>
                  <cn cellml:units="millisecond">59.5</cn>
                </apply>
              </apply>
            </piece>
            <otherwise>
              <ci> V_hold </ci>
            </otherwise>
          </piecewise>
        </apply>
         
        <apply>
        <eq/>
            <apply>
          <diff/>
               <bvar>
            <ci>time</ci>
          </bvar>
               <ci>V</ci>
            </apply>
            <apply>
          <times/>
               <cn cellml:units="dimensionless">-1E3</cn>
               <apply>
            <plus/>
                  <ci>I_Na</ci>
                  <ci>I_Na_b</ci>
                  <ci>I_NaCa</ci>
                  <ci>I_NaK</ci>
                  <ci>I_to</ci>
                  <ci>I_ss</ci>
                  <ci>I_Ki</ci>
                  <ci>I_Kp</ci>
                  <ci>I_CaK</ci>
                  <ci>I_Ca</ci>
                  <ci>I_Ca_b</ci>
                  <ci>I_p_Ca</ci>
                  <apply>
              <minus/>
                    <ci>I_app</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>

<!-- the following descriptions of the transient outward and steady state potassium currents have been taken from the Pandit et al. 2001 model. -->  

   <component name="Ca_independent_transient_outward_K_current">
      <variable units="microA_per_microF" public_interface="out" name="I_to"/>
      
      <variable units="milliS_per_microF" name="g_to" initial_value="0.35"/>
      <variable units="dimensionless" name="a" initial_value="0.886"/>
      <variable units="dimensionless" name="b" initial_value="0.114"/>
      
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      
      <variable units="dimensionless" private_interface="in" name="r"/>
      <variable units="dimensionless" private_interface="in" name="s"/>
      <variable units="dimensionless" private_interface="in" name="s_slow"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="8">
         <apply>
            <eq/>
            <ci>I_to</ci>
            <apply>
               <times/>
               <ci>g_to</ci>
               <ci>r</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>a</ci>
                     <ci>s</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>b</ci>
                     <ci>s_slow</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
  
   <component name="Ca_independent_transient_outward_K_current_r_gate">
      <variable units="dimensionless" public_interface="out" name="r" initial_value="1.4E-3"/>
      <variable units="second" name="tau_r"/>
      <variable units="dimensionless" name="r_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="9">
         <apply>
            <eq/>
            <ci>r_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">10.6</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">11.42</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
         <apply>
            <eq/>
            <ci>tau_r</ci>
            <apply>
               <divide/>
               <cn cellml:units="second">1</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">45.16</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.03577</cn>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">50</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">98.9</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">0.1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">38</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>r_infinity</ci>
                  <ci>r</ci>
               </apply>
               <ci>tau_r</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="Ca_independent_transient_outward_K_current_s_gate">
      <variable units="dimensionless" public_interface="out" name="s" initial_value="1.0"/>
      <variable units="second" name="tau_s"/>
      <variable units="dimensionless" name="s_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="10">
         <apply>
            <eq/>
            <ci>s_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">45.3</cn>
                        </apply>
                        <cn cellml:units="millivolt">6.8841</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_s</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">0.35</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">70</cn>
                              </apply>
                              <cn cellml:units="millivolt">15</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.035</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>s_infinity</ci>
                  <ci>s</ci>
               </apply>
               <ci>tau_s</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_independent_transient_outward_K_current_s_slow_gate">
      <variable units="dimensionless" public_interface="out" name="s_slow" initial_value="0.613"/>
      <variable units="second" name="tau_s_slow"/>
      <variable units="dimensionless" name="s_slow_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="11">
         <apply>
            <eq/>
            <ci>s_slow_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">45.3</cn>
                        </apply>
                        <cn cellml:units="millivolt">6.8841</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_s_slow</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <cn cellml:units="second">3.7</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">70</cn>
                              </apply>
                              <cn cellml:units="millivolt">30</cn>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="second">0.035</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s_slow</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>s_slow_infinity</ci>
                  <ci>s_slow</ci>
               </apply>
               <ci>tau_s_slow</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="steady_state_outward_K_current">
      <variable units="microA_per_microF" public_interface="out" name="I_ss"/>
      <variable units="milliS_per_microF" name="g_ss" initial_value="0.07"/>
      
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="dimensionless" private_interface="in" name="r_ss"/>
      <variable units="dimensionless" private_interface="in" name="s_ss"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="12">
         <apply>
            <eq/>
            <ci>I_ss</ci>
            <apply>
               <times/>
               <ci>g_ss</ci>
               <ci>r_ss</ci>
               <ci>s_ss</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="steady_state_outward_K_current_r_ss_gate">
      <variable units="dimensionless" public_interface="out" name="r_ss" initial_value="198E-3"/>
      <variable units="second" name="tau_r_ss"/>
      <variable units="dimensionless" name="r_ss_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>r_ss_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">11.5</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">11.82</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_r_ss</ci>
            <apply>
               <divide/>
               <cn cellml:units="second">10</cn>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">45.16</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">0.03577</cn>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">50</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">98.9</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">0.1</cn>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">38</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>r_ss</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>r_ss_infinity</ci>
                  <ci>r_ss</ci>
               </apply>
               <ci>tau_r_ss</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="steady_state_outward_K_current_s_ss_gate">
      <variable units="dimensionless" public_interface="out" name="s_ss" initial_value="0.43"/>
      <variable units="second" name="tau_s_ss"/>
      <variable units="dimensionless" name="s_ss_infinity"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="13">
         <apply>
            <eq/>
            <ci>s_ss_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <cn cellml:units="millivolt">87.5</cn>
                        </apply>
                        <cn cellml:units="millivolt">10.3</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_s_ss</ci>
            <cn cellml:units="second">2.1</cn>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>s_ss</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>s_ss_infinity</ci>
                  <ci>s_ss</ci>
               </apply>
               <ci>tau_s_ss</ci>
            </apply>
         </apply>
      </math>
   </component>

<!-- the following components have been taken from an adapted version of the Luo-Rudy model published by Puglisi and Bers 2001. -->

   <component name="fast_sodium_current">
      <variable units="microA_per_microF" public_interface="out" name="I_Na"/>
      <variable units="millivolt" public_interface="out" name="E_Na"/>
      
      <variable units="milliS_per_microF" name="g_Na" initial_value="8.0"/>
      
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      
      <variable units="dimensionless" private_interface="in" name="m"/>
      <variable units="dimensionless" private_interface="in" name="h"/>
      <variable units="dimensionless" private_interface="in" name="j"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="14">
         <apply>
            <eq/>
            <ci>E_Na</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <log/>
                  <apply>
                     <divide/>
                     <ci>Nao</ci>
                     <ci>Nai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         
         <apply>
            <eq/>
            <ci>I_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <ci>j</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable units="dimensionless" public_interface="out" name="m" initial_value="1.4E-3"/>
      <variable units="first_order_rate_constant" name="alpha_m"/>
      <variable units="first_order_rate_constant" name="beta_m"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="15">
         <apply>
            <eq/>
            <ci>alpha_m</ci>
            <apply>
          <divide/>
              <apply>
            <times/>
                <cn cellml:units="per_millivolt_second">0.32</cn>
                <apply>
              <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">47.13</cn>
                </apply>
              </apply>
              <apply>
            <minus/>
                <cn cellml:units="dimensionless">1</cn>
                <apply>
              <exp/>
                  <apply>
                <times/>
                    <apply>
                  <minus/>
                      <cn cellml:units="per_millivolt">0.1</cn>
                    </apply>
                    <apply>
                  <plus/>
                      <ci>V</ci>
                      <cn cellml:units="millivolt">47.13</cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
         
          <apply>
            <eq/>
            <ci>beta_m</ci>
            <apply>
               <times/>
               <cn cellml:units="first_order_rate_constant">0.08</cn>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                     </apply>
                     <cn cellml:units="millivolt">11</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
          <times/>
              <cn cellml:units="dimensionless">1E3</cn>
              <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_m</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>m</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_m</ci>
                  <ci>m</ci>
               </apply>
              </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_h_gate">
      <variable units="dimensionless" public_interface="out" name="h" initial_value="0.99"/>
      <variable units="first_order_rate_constant" name="alpha_h"/>
      <variable units="first_order_rate_constant" name="beta_h"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="16">
         <apply>
            <eq/>
            <ci>alpha_h</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <cn cellml:units="first_order_rate_constant">0.135</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <cn cellml:units="millivolt">80</cn>
                              <ci>V</ci>
                           </apply>
                           <apply>
                              <minus/>
                              <cn cellml:units="millivolt">6.8</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="first_order_rate_constant">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>beta_h</ci>
            <piecewise>
               <piece>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="first_order_rate_constant">3.56</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.079</cn>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="first_order_rate_constant">3.1E5</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt">0.35</cn>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <cn cellml:units="first_order_rate_constant">1.0</cn>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">0.13</cn>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <plus/>
                                    <ci>V</ci>
                                    <cn cellml:units="millivolt">10.66</cn>
                                 </apply>
                                 <apply>
                                    <minus/>
                                    <cn cellml:units="millivolt">11.1</cn>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
          <times/>
              <cn cellml:units="dimensionless">1E3</cn>
              <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_h</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>h</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_h</ci>
                  <ci>h</ci>
               </apply>
              </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_j_gate">
      <variable units="dimensionless" public_interface="out" name="j" initial_value="0.99"/>
      <variable units="first_order_rate_constant" name="alpha_j"/>
      <variable units="first_order_rate_constant" name="beta_j"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="17">
         <apply>
            <eq/>
            <ci>alpha_j</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt_second">1.2714E5</cn>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="per_millivolt">0.2444</cn>
                                    <ci>V</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="per_millivolt_second">3.474e-5</cn>
                              <apply>
                                 <exp/>
                                 <apply>
                                    <times/>
                                    <apply>
                                       <minus/>
                                       <cn cellml:units="per_millivolt">0.04391</cn>
                                    </apply>
                                    <ci>V</ci>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">37.78</cn>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <times/>
                           <cn cellml:units="per_millivolt">0.311</cn>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <cn cellml:units="millivolt">79.23</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="first_order_rate_constant">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         
         <apply>
            <eq/>
            <ci>beta_j</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="first_order_rate_constant">0.1212</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.01052</cn>
                              </apply>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1378</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">40.14</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>V</ci>
                     <apply>
                        <minus/>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <cn cellml:units="first_order_rate_constant">0.3</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">2.535e-7</cn>
                              </apply>
                              <ci>V</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="per_millivolt">0.1</cn>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <cn cellml:units="millivolt">32</cn>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>j</ci>
            </apply>
            <apply>
          <times/>
              <cn cellml:units="dimensionless">1E3</cn>
              <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha_j</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>j</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_j</ci>
                  <ci>j</ci>
               </apply>
              </apply>
            </apply>
         </apply>
      </math>
   </component>

 <!-- the following description of the L-type calcium channel has been adapted from the Jafri, Rice, and Winslow 1998 model. -->  
   
   <component name="L_type_Ca_channel">
    <variable units="microA_per_microF" public_interface="out" name="I_Ca"/>
    <variable units="microA_per_microF" public_interface="out" name="I_CaK"/>
    
    <variable units="first_order_rate_constant" private_interface="out" name="alpha"/>
    <variable units="first_order_rate_constant" private_interface="out" name="beta"/>
    <variable units="first_order_rate_constant" private_interface="out" name="gamma"/>
    
    <variable units="cm_per_second" name="P_Ca" initial_value="1.7469E-8"/>
    <variable units="cm_per_second" name="P_K" initial_value="3.234E-11"/>
    <variable units="microA_per_microF" name="I_Ca05" initial_value="-0.458"/>
    <variable units="microA_per_microF" name="ibarca"/>
    <variable units="microA_per_microF" name="ibark"/>
    <variable units="second_order_rate_constant" name="gammao" initial_value="5187.5"/>
    <variable units="dimensionless" name="Nlcc" initial_value="3E5"/>
    <variable units="dimensionless" name="favail"/>  
    
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="dimensionless" public_interface="in" name="fracLCCbp"/>
    <variable units="dimensionless" public_interface="in" name="fracLCCbpo"/>
    <variable units="dimensionless" public_interface="in" private_interface="out" name="fracLCCap"/>
    <variable units="dimensionless" public_interface="in" private_interface="out" name="fracLCCapo"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
	<variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
	<variable units="millimolar" public_interface="in" name="Ki"/>
    <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    
    <variable units="dimensionless" private_interface="in" name="y"/>
    <variable units="dimensionless" private_interface="in" name="v"/>
    <variable units="dimensionless" private_interface="in" name="z"/>
	<variable units="dimensionless" private_interface="in" name="x"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="18">    
      <apply id="I_Ca_calculation">
        <eq/>
        <ci> I_Ca </ci>
        <apply>
          <times/>
          <ci> ibarca </ci>
          <ci> Nlcc </ci>
          <ci> favail </ci>
          <apply>
            <power/>
            <ci> v </ci>
            <cn cellml:units="dimensionless"> 4.0 </cn>
          </apply>
          <ci> x </ci>
          <ci> y </ci>
          <ci> z </ci>
        </apply>
      </apply>

      <apply id="I_CaK_calculation">
        <eq/>
        <ci> I_CaK </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <ci> ibark </ci>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> I_Ca </ci>
                <ci> I_Ca05 </ci>
              </apply>
            </apply>
          </apply>
          <ci> Nlcc </ci>
          <ci> favail </ci>
          <apply>
            <power/>
            <ci> v </ci>
            <cn cellml:units="dimensionless"> 4.0 </cn>
          </apply>
          <ci> x </ci>
          <ci> y </ci>
          <ci> z </ci>
        </apply>
      </apply>
      
      <apply id="ibarca_calculation">
        <eq/>
        <ci> ibarca </ci>
        <apply>
          <times/>
          <ci> P_Ca </ci>
          <cn cellml:units="dimensionless"> 4.0 </cn>
          <apply>
            <divide/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> V </ci>
              <ci> F </ci>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.001 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> V </ci>
                    <apply>
                      <divide/>
                      <ci> F </ci>
                      <apply>
                        <times/>
                        <ci> R </ci>
                        <ci> T </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.341 </cn>
                <ci> Cao </ci>
              </apply>
            </apply> 
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> V </ci>
                  <apply>
                    <divide/>
                    <ci> F </ci>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply id="ibark_calculation">
        <eq/>
        <ci> ibark </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> P_K </ci>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> V </ci>
                <ci> F </ci>
              </apply>
              <apply>
                <times/>
                <ci> R </ci>
                <ci> T </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> Ki </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <ci> Ko </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> V </ci>
                  <ci> F </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> R </ci>
                  <ci> T </ci>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>

      <apply id="alpha_calculation">
        <eq/>
        <ci> alpha </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 400.0 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <plus/>
                <ci> V </ci>
                <cn cellml:units="millivolt"> 2.0 </cn>
              </apply>
              <cn cellml:units="millivolt"> 10.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_calculation">
        <eq/>
        <ci> beta </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 50.0 </cn>
          <apply>
            <exp/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> -1.0 </cn>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 2.0 </cn>
                </apply>
                <cn cellml:units="millivolt"> 13.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply id="gamma_calculation">
        <eq/>
        <ci> gamma </ci>
        <apply>
          <times/>
          <ci> gammao </ci>
          <ci> Cai </ci>
        </apply>
      </apply>

      <apply id="favail_calculation">
        <eq/>
        <ci> favail </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.5 </cn>
          <apply>
            <divide/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.4 </cn>
              <ci> fracLCCbp </ci>
            </apply>
            <apply>
              <plus/>
              <ci> fracLCCbpo </ci>
              <cn cellml:units="dimensionless"> 0.60 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="L_type_Ca_channel_v_gate">
    <variable units="dimensionless" public_interface="out" name="v" initial_value="0.0"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="alpha"/>
    <variable units="first_order_rate_constant" public_interface="in" name="beta"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="19">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> v </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> v </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta </ci>
            <ci> v </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="L_type_Ca_channel_w_gate">
    <variable units="dimensionless" public_interface="out" name="w" initial_value="0.0"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="alpha"/>
    <variable units="first_order_rate_constant" public_interface="in" name="beta"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="20">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> w </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> alpha </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> w </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> beta </ci>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> w </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="L_type_Ca_channel_x_gate">
    <variable units="dimensionless" public_interface="out" name="x" initial_value="0.13"/>
    <variable units="first_order_rate_constant" public_interface="out" name="g" initial_value="2E3"/>
    <variable units="first_order_rate_constant" public_interface="out" name="flcc"/>
    
    <variable units="first_order_rate_constant" name="f" initial_value="300.0"/>
    
    <variable units="dimensionless" public_interface="in" name="fracLCCap"/>
    <variable units="dimensionless" public_interface="in" name="fracLCCapo"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="21">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> x </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> flcc </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> x </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> g </ci>
            <ci> x </ci>
          </apply>
        </apply>
      </apply>

      <apply id="flcc_calculation">
        <eq/>
        <ci> flcc </ci>
        <apply>
          <times/>
          <ci> f </ci>
          <apply>
            <divide/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.375 </cn>
              <ci> fracLCCap </ci>
            </apply>
            <apply>
              <plus/>
              <ci> fracLCCapo </ci>
              <cn cellml:units="dimensionless"> 0.625 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
      
  <component name="L_type_Ca_channel_y_gate">
    <variable units="dimensionless" public_interface="out" name="y" initial_value="0.96"/>
    
    <variable units="dimensionless" name="y_infinity"/>
    <variable units="dimensionless" name="tau_y"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="22">
      <apply id="y_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> y </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> y_infinity </ci>
            <ci> y </ci>
          </apply>
          <ci> tau_y </ci>
        </apply>
      </apply>
      
      <apply id="y_infinity_calculation">
        <eq/>
        <ci> y_infinity </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 55.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 7.5 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 0.1 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <apply>
                      <minus/>
                      <ci> V </ci>
                    </apply>
                    <cn cellml:units="millivolt"> 21.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 6.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_y_calculation">
        <eq/>
        <ci> tau_y </ci>
        <apply>
          <plus/>
          <cn cellml:units="dimensionless"> 0.02 </cn>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 0.3 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 30.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 9.5 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="L_type_Ca_channel_z_gate">
    <variable units="dimensionless" public_interface="out" name="z" initial_value="0.92"/>
    
    <variable units="first_order_rate_constant" name="omega" initial_value="10.0"/>
    <variable units="first_order_rate_constant" name="v_omega"/>
    <variable units="second_order_rate_constant" name="v_gamma"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="gamma"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g"/>
    <variable units="dimensionless" private_interface="in" name="v"/>
    <variable units="first_order_rate_constant" public_interface="in" name="flcc"/>
    <variable units="dimensionless" public_interface="in" name="w"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="23">
      <apply id="z_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> z </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> v_omega </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> z </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> v_gamma </ci>
            <ci> z </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="v_omega_calculation">
        <eq/>
        <ci> v_omega </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> omega </ci>
            <apply>
              <power/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> w </ci>
              </apply>
              <cn cellml:units="dimensionless"> 4.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.5 </cn>
            <ci> w </ci>
            <apply>
              <power/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> w </ci>
              </apply>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.25 </cn>
            <apply>
              <power/>
              <ci> w </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> w </ci>
            </apply>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.125 </cn>
            <apply>
              <power/>
              <ci> w </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> w </ci>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.0625 </cn>
            <apply>
              <power/>
              <ci> w </ci>
              <cn cellml:units="dimensionless"> 4.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="v_gamma_calculation">
        <eq/>
        <ci> v_gamma </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> gamma </ci>
            <apply>
              <power/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> v </ci>
              </apply>
              <cn cellml:units="dimensionless"> 4.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> v </ci>
            <apply>
              <power/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> v </ci>
              </apply>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 4.0 </cn>
            <apply>
              <power/>
              <ci> v </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> v </ci>
            </apply>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 8.0 </cn>
            <apply>
              <power/>
              <ci> v </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> v </ci>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 16.0 </cn>
            <apply>
              <power/>
              <ci> v </ci>
              <cn cellml:units="dimensionless"> 4.0 </cn>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> flcc </ci>
                <ci> g </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <!-- back to the Puglisi and Bers 2001 model. -->
   
  <component name="time_independent_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="I_Ki"/>
      <variable units="millivolt" public_interface="out" private_interface="out" name="E_K"/>
      <variable units="milliS_per_microF" name="g_Ki" initial_value="0.24"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
	  <variable units="millimolar" public_interface="in" name="Ki"/>
	  <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="dimensionless" private_interface="in" name="Ki_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="24">
	    <apply>
            <eq/>
            <ci>E_K</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <log/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <ci>Ki</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
	  
         <apply>
        <eq/>
           <ci>I_Ki</ci>
           <apply>
          <times/>
             <ci>g_Ki</ci>
             <apply>
            <root/>
               <apply>
              <divide/>
                 <ci>Ko</ci>
                 <cn cellml:units="millimolar"> 5.4 </cn>
               </apply>
             </apply>
             <ci>Ki_infinity</ci>
             <apply>
            <minus/>
               <ci>V</ci>
               <ci>E_K</ci>
             </apply>
           </apply>
         </apply>
      </math>
   </component>
   
   <component name="time_independent_potassium_current_Ki_gate">
      <variable units="dimensionless" public_interface="out" name="Ki_infinity"/>
      <variable units="first_order_rate_constant" name="alpha_Ki"/>
      <variable units="first_order_rate_constant" name="beta_Ki"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="25">
         <apply>
            <eq/>
            <ci>alpha_Ki</ci>
            <apply>
               <divide/>
               <cn cellml:units="first_order_rate_constant">1.02</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <cn cellml:units="per_millivolt">0.2385</cn>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <ci>E_K</ci>
                           </apply>
                           <cn cellml:units="millivolt">59.215</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta_Ki</ci>
            <apply>
               <divide/>
               <apply>
            <plus/>
                  <apply>
              <times/>
                     <cn cellml:units="first_order_rate_constant">0.49124</cn>
                     <apply>
                <exp/>
                        <apply>
                  <times/>
                           <cn cellml:units="dimensionless">0.08032</cn>
                           <apply>
                    <plus/>
                              <apply>
                      <minus/>
                                 <ci>V</ci>
                                 <ci>E_K</ci>
                              </apply>
                              <cn cellml:units="millivolt">5.476</cn>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
              <exp/>
                     <apply>
                <times/>
                        <cn cellml:units="dimensionless">0.06175</cn>
                        <apply>
                  <plus/>
                           <apply>
                    <minus/>
                              <ci>V</ci>
                              <ci>E_K</ci>
                           </apply>
                           <cn cellml:units="millivolt">-594.31</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="per_millivolt">0.5143</cn>
                        </apply>
                        <apply>
                           <plus/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <ci>E_K</ci>
                           </apply>
                           <cn cellml:units="millivolt">4.753</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Ki_infinity</ci>
            <apply>
               <divide/>
               <ci>alpha_Ki</ci>
               <apply>
                  <plus/>
                  <ci>alpha_Ki</ci>
                  <ci>beta_Ki</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="plateau_potassium_current">
      <variable units="microA_per_microF" public_interface="out" name="I_Kp"/>
      <variable units="milliS_per_microF" name="g_Kp" initial_value="0.008"/>
      <variable units="dimensionless" name="Kp"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="26">
         <apply>
            <eq/>
            <ci>Kp</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <cn cellml:units="millivolt">7.488</cn>
                           <ci>V</ci>
                        </apply>
                        <cn cellml:units="millivolt">5.98</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_Kp</ci>
            <apply>
               <times/>
               <ci>g_Kp</ci>
               <ci>Kp</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="sarcolemmal_calcium_pump">
      <variable units="microA_per_microF" public_interface="out" name="I_p_Ca"/>
      <variable units="millimolar" name="K_mpCa" initial_value="0.5E-3"/>
      <variable units="microA_per_microF" name="I_pCa" initial_value="1.15"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="27">
         <apply>
            <eq/>
            <ci>I_p_Ca</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>I_pCa</ci>
                  <ci>Cai</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_mpCa</ci>
                  <ci>Cai</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="sodium_background_current">
      <variable units="microA_per_microF" public_interface="out" name="I_Na_b"/>
      <variable units="microA_per_microF" name="g_Nab" initial_value="1.18E-3"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="28">
         <apply>
            <eq/>
            <ci>I_Na_b</ci>
            <apply>
               <times/>
               <ci>g_Nab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="calcium_background_current">
      <variable units="microA_per_microF" public_interface="out" name="I_Ca_b"/>
      <variable units="microA_per_microF" name="g_Cab" initial_value="2.8E-3"/>
      <variable units="millivolt" public_interface="out" name="E_Ca"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <variable units="millimolar" public_interface="in" name="Cao"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="29">
         <apply>
            <eq/>
            <ci>E_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>F</ci>
                  </apply>
               </apply>
               <apply>
                  <log/>
                  <apply>
                     <divide/>
                     <ci>Cao</ci>
                     <ci>Cai</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_Ca_b</ci>
            <apply>
               <times/>
               <ci>g_Cab</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="sodium_potassium_pump">
      <variable units="microA_per_microF" public_interface="out" name="I_NaK"/>
      
      <variable units="microA_per_microF" name="i_NaK" initial_value="1.1"/>
      <variable units="dimensionless" name="f_NaK"/>
      <variable units="millimolar" name="K_mNai" initial_value="10"/>
      <variable units="millimolar" name="K_mKo" initial_value="1.5"/>
      <variable units="dimensionless" name="sigma"/>
      
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
      <variable units="kelvin" public_interface="in" name="T"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="millimolar" public_interface="in" name="Nai"/>
      <variable units="millimolar" public_interface="in" name="Nao"/>
      <variable units="millimolar" public_interface="in" name="Ko"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="30">
         <apply>
            <eq/>
            <ci>sigma</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <cn cellml:units="dimensionless">7</cn>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <ci>Nao</ci>
                        <cn cellml:units="dimensionless">67.3</cn>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>f_NaK</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.1245</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">0.1</cn>
                              </apply>
                              <ci>V</ci>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.0365</cn>
                     <ci>sigma</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                              </apply>
                              <ci>F</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>R</ci>
                              <ci>T</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>I_NaK</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>i_NaK</ci>
                        <ci>f_NaK</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <ci>K_mNai</ci>
                              <ci>Nai</ci>
                           </apply>
                           <cn cellml:units="dimensionless">1.5</cn>
                        </apply>
                     </apply>
                  </apply>
                  <ci>Ko</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Ko</ci>
                  <ci>K_mKo</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="Na_Ca_exchanger">
    <variable units="microA_per_microF" public_interface="out" name="I_NaCa"/>
    
    <variable units="microA_per_microF" name="K_NaCa" initial_value="1483.0"/>
    <variable units="millimolar" name="K_mNa" initial_value="87.5"/>
    <variable units="dimensionless" name="K_sat" initial_value="0.1"/>
    <variable units="dimensionless" name="eta" initial_value="0.35"/>
    <variable units="millimolar" name="K_mCa" initial_value="1.38"/>
   
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="31">
      <apply id="Na_Ca_exchanger">
        <eq/>
        <ci> I_NaCa </ci>
        <apply>
          <times/> 
          <ci> K_NaCa </ci>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <apply>
                <power/>
                <ci> K_mNa </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <apply>
                <power/>
                <ci> Nao </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <ci> K_mCa </ci>
              <ci> Cao </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <ci> K_sat </ci>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <ci> eta </ci>
                      <cn cellml:units="dimensionless"> 1.0 </cn>
                    </apply>
                    <ci> V </ci>
                    <apply>
                      <divide/>
                      <ci> F </ci>
                      <apply>
                        <times/>
                        <ci> R </ci>
                        <ci> T </ci>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <ci> eta </ci>
                  <ci> V </ci>
                  <apply>
                    <divide/>
                    <ci> F </ci>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> Nai </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <ci> Cao </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <ci> eta </ci>
                    <cn cellml:units="dimensionless"> 1.0 </cn>
                  </apply>
                  <ci> V </ci>
                  <apply>
                    <divide/>
                    <ci> F </ci>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> Nao </ci>
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <ci> Cai </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

   <component name="calcium_dynamics">
      <variable units="flux" public_interface="out" name="I_rel"/>
      <variable units="flux" public_interface="out" name="I_up"/>
      <variable units="flux" public_interface="out" name="I_leak"/>
      <variable units="flux" name="I_tr"/>
      <variable units="millimolar" name="Ca_JSR" initial_value="1.92"/>
      <variable units="millimolar" name="Ca_NSR" initial_value="1.92"/>
      <variable units="dimensionless" public_interface="out" name="B_myo"/>
      
      <variable units="first_order_rate_constant" name="G_rel"/>
      <variable units="flux" name="G_rel_max" initial_value="60E3"/>
      <variable units="first_order_rate_constant" name="G_rel_overload" initial_value="4000"/>
      <variable units="second" name="tau_tr" initial_value="5.7E-4"/>
      <variable units="millimolar" name="K_mrel" initial_value="0.8E-3"/>
      <variable units="millimolar" name="delta_Ca_ith" initial_value="0.18E-3"/>
      <variable units="millimolar" name="CSQN_max" initial_value="15"/>
      <variable units="millimolar" name="K_mCSQN" initial_value="0.8"/>
      <variable units="millimolar" name="K_mup" initial_value="3E-4"/>
      <variable units="dimensionless" name="Km_up"/>
      <variable units="flux" name="i_up" initial_value="4.7"/>
      <variable units="millimolar" name="Ca_NSR_max" initial_value="15"/>
      <variable units="millimolar" name="Km_Tn" initial_value="0.5128E-3"/>
      <variable units="millimolar" name="Km_CMDN" initial_value="2.38E-3"/>
      <variable units="millimolar" name="Km_INDO" initial_value="8.44E-4"/>
      <variable units="millimolar" name="Tn_max" initial_value="0.07"/>
      <variable units="millimolar" name="CMDN_max" initial_value="0.05"/>
      <variable units="millimolar" name="INDO_max" initial_value="0.07"/>
      <variable units="dimensionless" name="b_CMDN"/>
      <variable units="dimensionless" name="b_Tn"/>
      <variable units="dimensionless" name="b_INDO"/>
      <variable units="dimensionless" name="ryr_on"/>
      <variable units="dimensionless" name="ryr_off"/>
      <variable units="second" name="trel"/>
      <variable units="second" name="y44" initial_value="0.898"/>
      <variable units="micromolar" name="SR_content"/>
      <variable units="second" name="tau_on" initial_value="2E-3"/>
      <variable units="second" name="tau_off" initial_value="2E-3"/>
      
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="coulomb_per_mole" public_interface="in" name="F"/>
      <variable units="microlitre" public_interface="in" name="V_myo"/>
      <variable units="microlitre" public_interface="in" name="V_JSR"/>
      <variable units="microlitre" public_interface="in" name="V_NSR"/>
      <variable units="cm2" public_interface="in" name="A_cap"/>
      <variable units="microA_per_microF" public_interface="in" name="I_NaCa"/>
      <variable units="microA_per_microF" public_interface="in" name="I_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="I_p_Ca"/>
      <variable units="microA_per_microF" public_interface="in" name="I_Ca_b"/>
      <variable units="millimolar" public_interface="in" name="Cai"/>
      <variable units="millimolar" public_interface="in" name="Cao"/>
      <variable units="dimensionless" public_interface="in" name="fracPLB"/>
      <variable units="dimensionless" public_interface="in" name="fracPLBo"/>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="32">
         <apply>
        <eq/>
            <ci>ryr_on</ci>
            <apply>
          <minus/>
               <cn cellml:units="dimensionless"> 2.0 </cn>
               <apply>
            <exp/>
                  <apply>
              <divide/>
                     <apply>
                <minus/>
                        <ci>trel</ci>
                     </apply>
                     <ci>tau_on</ci>
                  </apply>
               </apply>
            </apply>
         </apply>

         <apply>
        <eq/>
            <ci>ryr_off</ci>
            <apply>
          <exp/>
               <apply>
            <divide/>
                  <apply>
              <minus/>
                     <ci>trel</ci>
                  </apply>
                  <ci>tau_off</ci>
               </apply>
            </apply>
         </apply>

         <apply>
        <eq/>
            <ci>I_rel</ci>
            <apply>
          <times/>
               <ci>G_rel</ci>
               <ci>ryr_on</ci>
               <ci>ryr_off</ci>
               <apply>
            <minus/>
                  <ci>Ca_JSR</ci>
                  <ci>Cai</ci>
               </apply>
            </apply>
         </apply>
		 
		 <apply>
        <eq/>
            <ci>trel</ci>
            <apply>
          <plus/>
               <ci>y44</ci>
               <cn cellml:units="second"> 2E-3 </cn>
            </apply>
         </apply>

      <apply>
        <eq/>
        <ci> G_rel </ci>
        <apply>
          <divide/>
          <ci> G_rel_max </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <cn cellml:units="microA_per_microF"> 5.0 </cn>
                    <ci>I_Ca</ci>
                    <ci>I_Ca_b</ci>
                    <ci>I_p_Ca</ci>
                  </apply>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless">2</cn>
                    <ci>I_NaCa</ci>
                  </apply>
                </apply>
                <cn cellml:units="microA_per_microF"> 0.9 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
         
         <apply>
        <eq/>
            <ci>I_up</ci>
            <apply>
          <divide/>
               <apply>
            <times/>
                  <ci>i_up</ci>
                  <apply>
              <power/>
                     <ci>Cai</ci>
                     <cn cellml:units="dimensionless"> 2.0 </cn>
                  </apply>
               </apply>
               <apply>
            <plus/>
                  <apply>
              <power/>
                     <ci>K_mup</ci>
                     <cn cellml:units="dimensionless"> 2.0 </cn>
                  </apply>
                  <apply>
              <power/>
                     <ci>Cai</ci>
                     <cn cellml:units="dimensionless"> 2.0 </cn>
                  </apply>
               </apply>
            </apply>
         </apply>

         <apply>
        <eq/>
            <ci>Km_up</ci>
            <apply>
          <times/>
               <ci>K_mup</ci>
               <apply>
            <divide/>
                  <apply>
              <plus/>
                     <cn cellml:units="dimensionless"> 1.0 </cn>
                     <apply>
                <times/>
                        <cn cellml:units="dimensionless"> 2.0 </cn>
                        <ci>fracPLB</ci>
                     </apply>
                  </apply>
                  <apply>
              <plus/>
                     <cn cellml:units="dimensionless"> 1.0 </cn>
                     <apply>
                <times/>
                        <cn cellml:units="dimensionless"> 2.0 </cn>
                        <ci>fracPLBo</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
         <apply>
        <eq/>
            <ci>I_leak</ci>
            <apply>
          <divide/>
               <apply>
            <times/>
                  <ci>i_up</ci>
                  <ci>Ca_NSR_max</ci>
               </apply>
               <ci>Ca_NSR</ci>
            </apply>
         </apply>
         
         <apply>
        <eq/>
            <ci>I_tr</ci>
            <apply>
          <divide/>
               <apply>
            <minus/>
                  <ci>Ca_NSR</ci>
                  <ci>Ca_JSR</ci>
               </apply>
               <ci>tau_tr</ci>
            </apply>
         </apply>
         
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_JSR</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>CSQN_max</ci>
                           <ci>K_mCSQN</ci>
                        </apply>
                        <apply>
                           <power/>
                           <apply>
                              <plus/>
                              <ci>K_mCSQN</ci>
                              <ci>Ca_JSR</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>I_tr</ci>
                  <ci>I_rel</ci>
               </apply>
            </apply>
         </apply>
         
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_NSR</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>I_tr</ci>
                        </apply>
                        <ci>V_JSR</ci>
                     </apply>
                     <ci>V_NSR</ci>
                  </apply>
                  <ci>I_leak</ci>
               </apply>
               <ci>I_up</ci>
            </apply>
         </apply>

         <apply>
        <eq/>
            <ci>SR_content</ci>
            <apply>
          <times/>
               <cn cellml:units="dimensionless">1E3</cn>
               <apply>
            <plus/>
                  <apply>
              <times/>
                     <apply>
                <plus/>
                        <ci>Ca_JSR</ci>
                        <apply>
                  <divide/>
                           <ci>Ca_JSR</ci>
                           <apply>
                    <divide/>
                              <cn cellml:units="dimensionless">1</cn>
                              <apply>
                      <plus/>
                                 <cn cellml:units="dimensionless">1</cn>
                                 <apply>
                        <divide/>
                                    <apply>
                          <times/>
                                       <ci>CSQN_max</ci>
                                       <ci>K_mCSQN</ci>
                                    </apply>
                                    <apply>
                          <power/>
                                       <apply>
                            <plus/>
                                          <ci>K_mCSQN</ci>
                                          <ci>Ca_JSR</ci>
                                       </apply>
                                       <cn cellml:units="dimensionless">2</cn>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                <divide/>
                        <ci>V_JSR</ci>
                        <ci>V_myo</ci>
                     </apply>
                  </apply>
                  <apply>
              <times/>
                     <ci>Ca_NSR</ci>
                     <apply>
                <divide/>
                        <ci>V_NSR</ci>
                        <ci>V_myo</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>

      <apply>
        <eq/> 
        <ci> B_myo </ci> 
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> b_CMDN </ci>
            <ci> b_Tn </ci>
            <ci> b_Tn </ci>
            <ci> b_INDO </ci>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> b_CMDN </ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> CMDN_max </ci>
            <ci> Km_CMDN </ci>
          </apply>
          <apply>
            <power/>
            <apply>
              <plus/>
              <ci> Km_CMDN </ci>
              <ci> Cai </ci>
            </apply>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
        </apply>
       </apply>

       <apply>
        <eq/> 
        <ci> b_Tn </ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> Tn_max </ci>
            <ci> Km_Tn </ci>
          </apply>
          <apply>
            <power/>
            <apply>
              <plus/>
              <ci> Km_Tn </ci>
              <ci> Cai </ci>
            </apply>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
        </apply>
       </apply>

       <apply>
        <eq/> 
        <ci> b_INDO </ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> INDO_max </ci>
            <ci> Km_INDO </ci>
          </apply>
          <apply>
            <power/>
            <apply>
              <plus/>
              <ci> Km_INDO </ci>
              <ci> Cai </ci>
            </apply>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
        </apply>
       </apply>

      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci> y44 </ci>
        </apply>
        <piecewise>
          <piece>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1 </cn>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1E4 </cn>
                <ci> y44 </ci>
              </apply>
            </apply>
            <apply>
              <gt/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 30E3 </cn>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="second"> 1.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
     </math>
   </component>

   <component name="ionic_concentrations">
    <variable units="millimolar" public_interface="out" name="Nai" initial_value="16.0"/>
    <variable units="millimolar" public_interface="out" name="Nao" initial_value="140.0"/>
    <variable units="millimolar" public_interface="out" name="Cai" initial_value="1.58E-4"/>
    <variable units="millimolar" public_interface="out" name="Cao" initial_value="1.8"/>
    <variable units="millimolar" public_interface="out" name="Ki" initial_value="145.0"/>
    <variable units="millimolar" public_interface="out" name="Ko" initial_value="5.4"/>
    
    <variable units="dimensionless" public_interface="in" name="B_myo"/>
    <variable units="second" public_interface="in" name="time"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    <variable units="microA_per_microF" public_interface="in" name="I_Na"/>       
    <variable units="microA_per_microF" public_interface="in" name="I_Na_b"/>
    <variable units="microA_per_microF" public_interface="in" name="I_NaCa"/>
    <variable units="microA_per_microF" public_interface="in" name="I_NaK"/>
    <variable units="microA_per_microF" public_interface="in" name="I_Ca"/>
    <variable units="microA_per_microF" public_interface="in" name="I_CaK"/>
    <variable units="microA_per_microF" public_interface="in" name="I_p_Ca"/>
    <variable units="microA_per_microF" public_interface="in" name="I_Ca_b"/>
    <variable units="microA_per_microF" public_interface="in" name="I_ss"/>
    <variable units="microA_per_microF" public_interface="in" name="I_Ki"/>
    <variable units="microA_per_microF" public_interface="in" name="I_Kp"/>
    <variable units="microA_per_microF" public_interface="in" name="I_to"/>       
    <variable units="flux" public_interface="in" name="I_leak"/>   
    <variable units="flux" public_interface="in" name="I_up"/>
	<variable units="flux" public_interface="in" name="I_rel"/>
    <variable units="microlitre" public_interface="in" name="V_myo"/>
    <variable units="microlitre" public_interface="in" name="V_JSR"/>
    <variable units="microlitre" public_interface="in" name="V_NSR"/>
    <variable units="cm2" public_interface="in" name="A_cap"/>
              
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="33">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Nai</ci>
        </apply>
        <apply>
          <times/>
          <cn cellml:units="dimensionless">-1.0E3</cn>
          <apply>
            <plus/>
            <ci>I_Na</ci>
            <ci>I_Na_b</ci>
            <apply>
              <times/>
              <ci>I_NaCa</ci>
              <cn cellml:units="dimensionless">3</cn>
            </apply>
            <apply>
              <times/>
              <ci>I_NaK</ci>
              <cn cellml:units="dimensionless">3</cn>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci>A_cap</ci>
            <apply>
              <times/>
              <ci>V_myo</ci>
              <ci>F</ci>
            </apply>
          </apply>
        </apply>
      </apply>
         
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Ki</ci>
        </apply>
        <apply>
          <times/>
          <cn cellml:units="dimensionless">-1.0E3</cn>
          <apply>
            <plus/>
            <ci>I_to</ci>
            <ci>I_ss</ci>
            <ci>I_Ki</ci>
            <ci>I_Kp</ci>
            <ci>I_CaK</ci>
            <ci>I_Kp</ci>
            <ci>I_to</ci>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci>I_NaK</ci>
                <cn cellml:units="dimensionless">2</cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci>A_cap</ci>
            <apply>
              <times/>
              <ci>V_myo</ci>
              <ci>F</ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Cai</ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless">1E3</cn>
              <apply>
                <minus/>
                <apply>
                  <plus/>
                  <ci>I_Ca</ci>
                  <ci>I_Ca_b</ci>
                </apply>
                <apply>
                  <plus/>
                  <ci>I_p_Ca</ci>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless">2</cn>
                    <ci>I_NaCa</ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <divide/>
                <ci>A_cap</ci>
                <apply>
                  <times/>
                  <ci>V_myo</ci>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>F</ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <ci>I_up</ci>
                <ci>I_leak</ci>
              </apply>
              <apply>
                <divide/>
                <ci>V_NSR</ci>
                <ci>V_myo</ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <ci>B_myo</ci>
            <apply>
              <times/>
              <ci>I_rel</ci>
              <apply>
                <divide/>
                <ci>V_JSR</ci>
                <ci>V_myo</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="model_constants">
    <variable units="cm2" public_interface="out" name="A_cap" initial_value="1.534E-4"/>
    <variable units="microlitre" public_interface="out" name="V_myo" initial_value="20.8E-6"/>
    <variable units="microlitre" public_interface="out" name="V_JSR" initial_value="9.3E-8"/>
    <variable units="microlitre" public_interface="out" name="V_NSR" initial_value="9.88E-7"/>  
  </component>
  
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="fast_sodium_current">
            <component_ref component="fast_sodium_current_m_gate"/>
            <component_ref component="fast_sodium_current_h_gate"/>
            <component_ref component="fast_sodium_current_j_gate"/>
         </component_ref>
         <component_ref component="L_type_Ca_channel">
            <component_ref component="L_type_Ca_channel_v_gate"/>
            <component_ref component="L_type_Ca_channel_w_gate"/>
            <component_ref component="L_type_Ca_channel_x_gate"/>
            <component_ref component="L_type_Ca_channel_y_gate"/>
            <component_ref component="L_type_Ca_channel_z_gate"/>
         </component_ref>
         <component_ref component="time_independent_potassium_current">
            <component_ref component="time_independent_potassium_current_Ki_gate"/>
         </component_ref>
         <component_ref component="Ca_independent_transient_outward_K_current">
            <component_ref component="Ca_independent_transient_outward_K_current_r_gate"/>
            <component_ref component="Ca_independent_transient_outward_K_current_s_gate"/>
            <component_ref component="Ca_independent_transient_outward_K_current_s_slow_gate"/>
         </component_ref>
         <component_ref component="steady_state_outward_K_current">
            <component_ref component="steady_state_outward_K_current_r_ss_gate"/>
            <component_ref component="steady_state_outward_K_current_s_ss_gate"/>
         </component_ref>
         <component_ref component="Na_Ca_exchanger"/>
         <component_ref component="plateau_potassium_current"/>
         <component_ref component="sarcolemmal_calcium_pump"/>
         <component_ref component="sodium_background_current"/>
         <component_ref component="calcium_background_current"/>
         <component_ref component="sodium_potassium_pump"/>
         <component_ref component="ionic_concentrations"/>
         <component_ref component="calcium_dynamics"/>
      </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="fast_sodium_current_j_gate"/>
      </component_ref>
      <component_ref component="L_type_Ca_channel">
         <component_ref component="L_type_Ca_channel_v_gate"/>
         <component_ref component="L_type_Ca_channel_w_gate"/>
         <component_ref component="L_type_Ca_channel_x_gate"/>
         <component_ref component="L_type_Ca_channel_y_gate"/>
         <component_ref component="L_type_Ca_channel_z_gate"/>
      </component_ref>
      <component_ref component="time_independent_potassium_current">
         <component_ref component="time_independent_potassium_current_Ki_gate"/>
      </component_ref>
      <component_ref component="Ca_independent_transient_outward_K_current">
         <component_ref component="Ca_independent_transient_outward_K_current_r_gate"/>
         <component_ref component="Ca_independent_transient_outward_K_current_s_gate"/>
         <component_ref component="Ca_independent_transient_outward_K_current_s_slow_gate"/>
      </component_ref>
      <component_ref component="steady_state_outward_K_current">
         <component_ref component="steady_state_outward_K_current_r_ss_gate"/>
         <component_ref component="steady_state_outward_K_current_s_ss_gate"/>
      </component_ref>
   </group>

<!-- EP model connections -->
  
   <connection>
      <map_components component_2="environment" component_1="membrane"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="steady_state_outward_K_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="fast_sodium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="ionic_concentrations"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="calcium_dynamics"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_Na" variable_1="I_Na"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="Ca_independent_transient_outward_K_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_to" variable_1="I_to"/>
   </connection>
   <connection>
      <map_components component_2="steady_state_outward_K_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_ss" variable_1="I_ss"/>
   </connection>
   <connection>
      <map_components component_2="L_type_Ca_channel" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_Ca" variable_1="I_Ca"/>
      <map_variables variable_2="I_CaK" variable_1="I_CaK"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="Na_Ca_exchanger" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_NaCa" variable_1="I_NaCa"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_Ki" variable_1="I_Ki"/>
	  <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="plateau_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_Kp" variable_1="I_Kp"/>
   </connection>
   <connection>
      <map_components component_2="sarcolemmal_calcium_pump" component_1="membrane"/>
      <map_variables variable_2="I_p_Ca" variable_1="I_p_Ca"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_Na_b" variable_1="I_Na_b"/>
   </connection>
   <connection>
      <map_components component_2="calcium_background_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_Ca_b" variable_1="I_Ca_b"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="sodium_potassium_pump" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="I_NaK" variable_1="I_NaK"/>
      <map_variables variable_2="R" variable_1="R"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="membrane"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="membrane"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="I_to" variable_1="I_to"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="steady_state_outward_K_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="steady_state_outward_K_current"/>
      <map_variables variable_2="I_ss" variable_1="I_ss"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="fast_sodium_current"/>
      <map_variables variable_2="I_Na" variable_1="I_Na"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
   </connection>
   <connection>
      <map_components component_2="sodium_background_current" component_1="fast_sodium_current"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
	  <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
	  <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="I_Ca" variable_1="I_Ca"/>
      <map_variables variable_2="I_CaK" variable_1="I_CaK"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="I_Ca" variable_1="I_Ca"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="Na_Ca_exchanger"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
      <map_variables variable_2="I_NaCa" variable_1="I_NaCa"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="Na_Ca_exchanger"/>
      <map_variables variable_2="I_NaCa" variable_1="I_NaCa"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
	  <map_variables variable_2="Ki" variable_1="Ki"/>
      <map_variables variable_2="I_Ki" variable_1="I_Ki"/>
   </connection>
   <connection>
      <map_components component_2="time_independent_potassium_current" component_1="plateau_potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="plateau_potassium_current"/>
      <map_variables variable_2="I_Kp" variable_1="I_Kp"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="sarcolemmal_calcium_pump"/>
      <map_variables variable_2="I_p_Ca" variable_1="I_p_Ca"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="sarcolemmal_calcium_pump"/>
      <map_variables variable_2="I_p_Ca" variable_1="I_p_Ca"/>
	  <map_variables variable_2="Cai" variable_1="Cai"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="sodium_background_current"/>
      <map_variables variable_2="I_Na_b" variable_1="I_Na_b"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dynamics" component_1="calcium_background_current"/>
      <map_variables variable_2="I_Ca_b" variable_1="I_Ca_b"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="calcium_background_current"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
      <map_variables variable_2="I_Ca_b" variable_1="I_Ca_b"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="sodium_potassium_pump"/>
      <map_variables variable_2="Nai" variable_1="Nai"/>
      <map_variables variable_2="Nao" variable_1="Nao"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="I_NaK" variable_1="I_NaK"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_m_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="m" variable_1="m"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_h_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="h" variable_1="h"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_j_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="j" variable_1="j"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   
   <connection>
      <map_components component_2="L_type_Ca_channel_y_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="y" variable_1="y"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>

   <connection>
      <map_components component_2="L_type_Ca_channel_v_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="v" variable_1="v"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="alpha" variable_1="alpha"/>
      <map_variables variable_2="beta" variable_1="beta"/>
   </connection>

   <connection>
      <map_components component_2="L_type_Ca_channel_w_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="alpha" variable_1="alpha"/>
      <map_variables variable_2="beta" variable_1="beta"/>
   </connection>
  
   <connection>
      <map_components component_2="L_type_Ca_channel_x_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="x" variable_1="x"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="fracLCCap" variable_1="fracLCCap"/>
      <map_variables variable_2="fracLCCapo" variable_1="fracLCCapo"/>
   </connection>

   <connection>
      <map_components component_2="L_type_Ca_channel_z_gate" component_1="L_type_Ca_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="z" variable_1="z"/>
	  <map_variables variable_2="gamma" variable_1="gamma"/>
   </connection>

   <connection>
      <map_components component_2="L_type_Ca_channel_x_gate" component_1="L_type_Ca_channel_z_gate"/>
      <map_variables variable_2="g" variable_1="g"/>
	  <map_variables variable_2="flcc" variable_1="flcc"/>
   </connection>
   
   <connection>
      <map_components component_2="L_type_Ca_channel_w_gate" component_1="L_type_Ca_channel_z_gate"/>
      <map_variables variable_2="w" variable_1="w"/>
   </connection>
   
   <connection>
      <map_components component_2="L_type_Ca_channel_v_gate" component_1="L_type_Ca_channel_z_gate"/>
      <map_variables variable_2="v" variable_1="v"/>
   </connection>

   <connection>
      <map_components component_2="time_independent_potassium_current_Ki_gate" component_1="time_independent_potassium_current"/>
      <map_variables variable_2="Ki_infinity" variable_1="Ki_infinity"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="Ca_independent_transient_outward_K_current_r_gate" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="r" variable_1="r"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="Ca_independent_transient_outward_K_current_s_gate" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="s" variable_1="s"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="Ca_independent_transient_outward_K_current_s_slow_gate" component_1="Ca_independent_transient_outward_K_current"/>
      <map_variables variable_2="s_slow" variable_1="s_slow"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="steady_state_outward_K_current_r_ss_gate" component_1="steady_state_outward_K_current"/>
      <map_variables variable_2="r_ss" variable_1="r_ss"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="steady_state_outward_K_current_s_ss_gate" component_1="steady_state_outward_K_current"/>
      <map_variables variable_2="s_ss" variable_1="s_ss"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="ionic_concentrations" component_1="calcium_dynamics"/>
      <map_variables variable_2="I_leak" variable_1="I_leak"/>
      <map_variables variable_2="I_up" variable_1="I_up"/>
	  <map_variables variable_2="I_rel" variable_1="I_rel"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
      <map_variables variable_2="Cao" variable_1="Cao"/>
      <map_variables variable_2="B_myo" variable_1="B_myo"/>
   </connection>

   <connection>
      <map_components component_2="model_constants" component_1="calcium_dynamics"/>
      <map_variables variable_2="A_cap" variable_1="A_cap"/>
      <map_variables variable_2="V_NSR" variable_1="V_NSR"/>
      <map_variables variable_2="V_JSR" variable_1="V_JSR"/>
      <map_variables variable_2="V_myo" variable_1="V_myo"/>
   </connection>

   <connection>
      <map_components component_2="model_constants" component_1="ionic_concentrations"/>
      <map_variables variable_2="A_cap" variable_1="A_cap"/>
      <map_variables variable_2="V_NSR" variable_1="V_NSR"/>
      <map_variables variable_2="V_JSR" variable_1="V_JSR"/>
      <map_variables variable_2="V_myo" variable_1="V_myo"/>
   </connection>

  
<!-- Signalling model connections -->

  <connection>
    <map_components component_2="environment" component_1="beta_adrenergic_receptor_module"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Gs_activation_module"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="cyclic_AMP_metabolism_module"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="PKA_activation_module"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="phospholamban_regulation_module"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="L_type_Ca_channel_regulation_module"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="Gs_activation_module" component_1="beta_adrenergic_receptor_module"/>
    <map_variables variable_2="Gs_beta_gamma" variable_1="Gs_beta_gamma"/>
    <map_variables variable_2="LRG" variable_1="LRG"/>
    <map_variables variable_2="RG" variable_1="RG"/>
  </connection>
  
  <connection>
    <map_components component_2="PKA_activation_module" component_1="beta_adrenergic_receptor_module"/>
    <map_variables variable_2="PKACI" variable_1="PKACI"/>
  </connection>
  
  <connection>
    <map_components component_2="Gs_activation_module" component_1="cyclic_AMP_metabolism_module"/>
    <map_variables variable_2="Gs_alpha_GTPtot" variable_1="Gs_alpha_GTPtot"/>
  </connection>
  
  <connection>
    <map_components component_2="PKA_activation_module" component_1="cyclic_AMP_metabolism_module"/>
    <map_variables variable_2="cAMP" variable_1="cAMP"/>
    <map_variables variable_2="cAMPtot" variable_1="cAMPtot"/>
  </connection>
  
  <connection>
    <map_components component_2="phospholamban_regulation_module" component_1="PKA_activation_module"/>
    <map_variables variable_2="PKACI" variable_1="PKACI"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_channel_regulation_module" component_1="PKA_activation_module"/>
    <map_variables variable_2="PKACII" variable_1="PKACII"/>
  </connection>

<!-- connecting the two models -->

  <connection>
    <map_components component_2="L_type_Ca_channel_regulation_module" component_1="L_type_Ca_channel"/>
    <map_variables variable_2="fracLCCbp" variable_1="fracLCCbp"/>
    <map_variables variable_2="fracLCCbpo" variable_1="fracLCCbpo"/>
    <map_variables variable_2="fracLCCap" variable_1="fracLCCap"/>
    <map_variables variable_2="fracLCCapo" variable_1="fracLCCapo"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_dynamics" component_1="phospholamban_regulation_module"/>
    <map_variables variable_2="fracPLB" variable_1="fracPLB"/>
    <map_variables variable_2="fracPLBo" variable_1="fracPLBo"/>
  </connection>
  


<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <rdf:Bag rdf:about="rdf:#a8556578-c204-404c-b0c9-2cfca30cb5f7">
    <rdf:li>Cardiac Myocyte</rdf:li>
    <rdf:li>metabolism</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>mechanics</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#ad72bc9a-a2c0-48d8-9ed5-a63e549fef6c">
    <rdf:li rdf:resource="rdf:#be31d7e4-9620-403d-bffb-cc0554221966"/>
    <rdf:li rdf:resource="rdf:#53f1d965-ecf1-4cb5-9cc6-e720ae5578d1"/>
    <rdf:li rdf:resource="rdf:#e419231f-fc6c-4175-b366-6494523e15a3"/>
    <rdf:li rdf:resource="rdf:#3cec4432-f09e-492b-980d-41c3f325f8cd"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#3cec4432-f09e-492b-980d-41c3f325f8cd">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#470c1dbb-3c39-4107-ba57-b0135c5eae1a"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#f78bffe1-cd0f-4ec3-86c9-3f36f876f6cd"/>
    <dcterms:created rdf:resource="rdf:#49f4d8ee-41b0-4032-9264-853069efd766"/>
    <dc:creator rdf:resource="rdf:#a55f9784-78de-4dfe-a483-4e05b2c0c605"/>
    <cmeta:modification rdf:resource="rdf:#f7c8a1c3-6c98-41bf-85e2-1de40485d219"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#53f1d965-ecf1-4cb5-9cc6-e720ae5578d1">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#81c1dfab-72c7-40d0-b187-ccf727ed171d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9db2ab2d-13dd-42aa-b8ef-adb3f0a0082e">
    <dc:subject rdf:resource="rdf:#219fdb75-145a-498a-bf98-ec6a2091d1db"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0de86d89-0ae0-44d6-8038-df99a786eded">
    <vCard:Given>Anushka</vCard:Given>
    <vCard:Family>Michailova</vCard:Family>
    <vCard:Other>P</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f7c8a1c3-6c98-41bf-85e2-1de40485d219">
    <dcterms:modified rdf:resource="rdf:#5c4e73c9-cb5e-4e08-88f9-c26ba6fe32a1"/>
    <rdf:value>The model has been drastically changed.  An electrophysiological excitable cell component has been added, and the author's source code was used to add initial values etc.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#0b6603ba-a263-405e-9372-eddfb8b986b4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#470c1dbb-3c39-4107-ba57-b0135c5eae1a">
    <vCard:Given>Andrew</vCard:Given>
    <vCard:Family>McCulloch</vCard:Family>
    <vCard:Other>D</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#be31d7e4-9620-403d-bffb-cc0554221966">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#481509e2-a814-4d47-8605-2db1a7cb08bf"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c4e73c9-cb5e-4e08-88f9-c26ba6fe32a1">
    <dcterms:W3CDTF>2007-09-21T11:43:40+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#481509e2-a814-4d47-8605-2db1a7cb08bf">
    <vCard:Given>Jeffrey</vCard:Given>
    <vCard:Family>Saucerman</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0efe1568-9510-478e-97b5-d33d7d62e77d">
    <dcterms:W3CDTF>2003-09-12</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#219fdb75-145a-498a-bf98-ec6a2091d1db">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#a8556578-c204-404c-b0c9-2cfca30cb5f7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ec188067-d335-4afc-b707-1f3c17a2fe08">
    <dc:title>Journal of Biological Chemistry</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f78bffe1-cd0f-4ec3-86c9-3f36f876f6cd">
    <dc:creator rdf:resource="rdf:#ef1fa92d-7447-4331-9ed0-e1e1d093cbdb"/>
    <rdf:value>Although much closer now to the published model, this CellML version of the model still can't be run in PCEnv or COR due to the presence of circular arguments.  It's hoped that this problem will be addressed by the further development of the tools.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#saucerman_model_2003">
    <dc:title>
        Saucerman et al.'s 2003 mathematical model of beta-adrenergic control of cardiac myocyte contractility.
      </dc:title>
    <cmeta:bio_entity>Cardiac Myocyte</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#bc9da357-61a4-490f-98fa-682c8648607a"/>
    <bqs:reference rdf:resource="rdf:#9db2ab2d-13dd-42aa-b8ef-adb3f0a0082e"/>
    <bqs:reference rdf:resource="rdf:#ccc52e45-21b5-4448-84e5-f5a29a280b1a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#98fa7ace-c01f-44d2-9b9d-4874eafd1744">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2d548f1a-f7ac-4be9-9a87-c694dd1435b5">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a55f9784-78de-4dfe-a483-4e05b2c0c605">
    <vCard:ORG rdf:resource="rdf:#2d548f1a-f7ac-4be9-9a87-c694dd1435b5"/>
    <vCard:EMAIL rdf:resource="rdf:#23ff5da7-0c81-4389-8ef9-8d71fcc139ec"/>
    <vCard:N rdf:resource="rdf:#415a7e10-3261-4cb8-b735-35883e7b7f43"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ccc52e45-21b5-4448-84e5-f5a29a280b1a">
    <bqs:Pubmed_id>12972422</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#3b339aac-040a-47a4-be05-c0ba8b381cb9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4f762dc2-c157-4e7f-8ec6-d38089f38348">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ef1fa92d-7447-4331-9ed0-e1e1d093cbdb">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#23ff5da7-0c81-4389-8ef9-8d71fcc139ec">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bc9da357-61a4-490f-98fa-682c8648607a">
    <dc:creator rdf:resource="rdf:#4f762dc2-c157-4e7f-8ec6-d38089f38348"/>
    <rdf:value>This is the CellML description of Saucerman et al.'s 2003 mathematical model of beta-adrenergic control of cardiac myocyte contractility.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#81c1dfab-72c7-40d0-b187-ccf727ed171d">
    <vCard:Given>Laurence</vCard:Given>
    <vCard:Family>Brunton</vCard:Family>
    <vCard:Other>L</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0b6603ba-a263-405e-9372-eddfb8b986b4">
    <vCard:N rdf:resource="rdf:#98fa7ace-c01f-44d2-9b9d-4874eafd1744"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#415a7e10-3261-4cb8-b735-35883e7b7f43">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3b339aac-040a-47a4-be05-c0ba8b381cb9">
    <dc:creator rdf:resource="rdf:#ad72bc9a-a2c0-48d8-9ed5-a63e549fef6c"/>
    <dc:title>Modeling beta-adrenergic control of cardiac myocyte contractility in silico.</dc:title>
    <bqs:volume>278</bqs:volume>
    <bqs:first_page>47997</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#ec188067-d335-4afc-b707-1f3c17a2fe08"/>
    <dcterms:issued rdf:resource="rdf:#0efe1568-9510-478e-97b5-d33d7d62e77d"/>
    <bqs:last_page>48003</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e419231f-fc6c-4175-b366-6494523e15a3">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#0de86d89-0ae0-44d6-8038-df99a786eded"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#49f4d8ee-41b0-4032-9264-853069efd766">
    <dcterms:W3CDTF>2003-11-18T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
</rdf:RDF>
</model>