Location: Rice, Winslow, Jafri, 1999 @ 7a1de088ee46 / rice_winslow_jafri_1999.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2006-07-09 07:52:38+12:00
Desc:
committing version01 of rice_winslow_jafri_1999
Permanent Source URI:
http://models.cellml.org/workspace/rice_winslow_jafri_1999/rawfile/7a1de088ee46241e1ee4a8a93bd5254c442eda51/rice_winslow_jafri_1999.cellml

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

CREATED :  29th June 2002

LAST MODIFIED : 11th September 2003

AUTHOR :  Catherine Lloyd
          The 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/01/2002 CellML Metadata 1.0 Specification.

DESCRIPTION :  This file contains a CellML description of Rice et al's 1999 model of the functional Ca2+ release unit.

CHANGES:  
  18/07/2002 - CML - Added more metadata.
  09/04/2003 - AAC - Added publication date information.  
  11/09/2003 - AAC - Fixed "J_RyR_calculation".  
--><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:pathway_editor="http://www.physiome.com/pathway_editor/1.0#" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" pathway_editor:rendering_config_file="RJW_model_1999_CellMLrender.xml" cmeta:id="rice_winslow_jafri_1999_version01" name="rice_winslow_jafri_1999_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modelling the Functional Ca2+ Release Unit</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 is the original unchecked version of the model imported from the previous
            CellML model repository, 24-Jan-2006.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
Calcium-induced calcium release (CICR) in cardiac muscle displays both gradedness and high gain.  It is graded because release of Ca<superscript>2+</superscript> from the sarcoplasmic reticulum (SR) is proportional to the influx of trigger Ca<superscript>2+</superscript>, and high gain refers to the observation that Ca<superscript>2+</superscript> release from the SR is much greater than the trigger influx.  Inherent to such high-gain systems is positive feedback, which causes an all-or-nothing response rather than graded Ca<superscript>2+</superscript> release.  Such behaviour is predicted by all mathematical models in which sarcolemmal Ca<superscript>2+</superscript> influx and SR Ca<superscript>2+</superscript> release are directed into a single compartment, or common pool model.  
</para>

<para>
The previous (1998) model of cardiac calcium handling, (see <ulink url="${HTML_EXMPL_JRW_MODEL}">Jafri-Rice-Winslow Ventricular Model, 1998</ulink>) by the authors was an example of a common pool model.  While it was able to simulate frequency-dependent aspects of cardiac Ca<superscript>2+</superscript> cycling and high gain, it also displayed an all-or-none rather than a graded Ca<superscript>2+</superscript> release.  By contrast, their 1999 model (see <xref linkend="fig_cell_diagram"/> below) is based on functional release units (FRUs) formulated from one sarcolemmal dihydropyridine receptor (DHPR), eight ryanodine receptor (RyR) channels, and the volume (subspace) surrounding them.  It is assumed that a variable number of independent FRUs can be recruited to produce graded Ca<superscript>2+</superscript> release.  Voltage-gated DHPR channels open adjacent RyRs via a localised increase in [Ca<superscript>2+</superscript>].
</para>

<para>
Results obtained with the stochastic model suggest that high gain and gradedness can occur by the recruitment of independent FRUs without requiring spatial [Ca<superscript>2+</superscript>] gradients within a functional unit, or coupling between adjacent units.
</para>

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

<para>
<ulink url="http://www.biophysj.org/cgi/content/abstract/77/4/1871">Modeling Gain and Gradedness of Ca<superscript>2+</superscript> Release in the Functional Unit of the Cardiac Diadic Space</ulink> John J. Rice, Saleet Jafri and Raimond L. Winslow, 1999, <ulink url="http://www.biophysj.org/">
            <emphasis>Biophysical Journal</emphasis>
          </ulink>, 77, 1871-1884.  (<ulink url="http://www.biophysj.org/cgi/content/full/77/4/1871">Full text</ulink> and <ulink url="http://www.biophysj.org/cgi/reprint/77/4/1871.pdf">PDF</ulink> versions of the article are available for Journal Members on the Biophysical Journal website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=10512809&amp;dopt=Abstract">PubMed ID: 10512809</ulink> 
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell_diagram</title>
    </objectinfo>
    <imagedata fileref="rice_1999.png"/>
  </imageobject>
</mediaobject>
<caption>A) Schematic diagram of the functional release unit model.  B) State diagram of the DHPR Channel.  C) State diagram of the RyR channel.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
  
  
  <!--
    Below, we define some additional units for association with variables and
    constants within the model. The identifiers are fairly self-explanatory.
  -->
  
  <units name="millisecond">
    <unit units="second" prefix="milli"/>
  </units>
  
  <units name="cm_per_millisecond">
    <unit units="metre" prefix="centi"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="millimolar">
    <unit units="mole" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="micrometre">
    <unit units="metre" prefix="micro"/>
  </units>
  
  <units name="per_micromolar">
    <unit units="micromolar" exponent="-1"/>
  </units>
  
  <units name="flux">
    <unit units="micromolar"/>
    <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="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>
  
  <units name="microA_per_microF">
    <unit units="ampere" prefix="micro"/>
    <unit units="farad" prefix="micro" exponent="-1"/>
  </units>
  
  <units name="joule_per_mole_kelvin">
    <unit units="joule"/>
    <unit units="mole" exponent="-1"/>
    <unit units="kelvin" exponent="-1"/>
  </units>
  
  <units name="coulomb_per_mole">
    <unit units="coulomb"/>
    <unit units="mole" exponent="-1"/>
  </units>
  
  <units name="micro_litre">
    <unit units="litre" prefix="micro"/>
  </units>
  
  <!--
    The "environment" component is used to declare variables that are used by
    all or most of the other components, in this case just "time".
  -->
  
  <component name="environment">
    <variable units="second" public_interface="out" name="time"/>
  </component>
  
  <component name="sarcolemma">
    <variable units="millivolt" public_interface="out" name="V" initial_value="-84.624"/>
    <variable units="joule_per_mole_kelvin" public_interface="out" name="R" initial_value="8.314"/>
    <variable units="kelvin" public_interface="out" name="T" initial_value="310.0"/>
    <variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96500.0"/>
    <variable units="millimolar" public_interface="out" name="Ca_o" initial_value="1.8"/>
  </component>
  
  <component name="subspace">
    <variable units="micro_litre" public_interface="out" name="V_SS" initial_value="2.03E-12"/>
    <variable units="flux" public_interface="out" name="J_RyR"/>
    <variable units="micromolar" public_interface="out" name="Ca_SS"/>
   
    <variable units="flux" name="J_DHPR"/>
    <variable units="dimensionless" name="DHPR_open"/>
    <variable units="dimensionless" name="RyR_open"/>
    <variable units="microA_per_microF" name="I_DHPR"/>
    <variable units="first_order_rate_constant" name="J_RyR_perm" initial_value="3960.0"/>
    <variable units="micromolar" name="beta_SS"/>
    <variable units="micromolar" name="B_SR" initial_value="47.0"/>
    <variable units="micromolar" name="K_BSR" initial_value="0.87"/>
    <variable units="micromolar" name="B_SL" initial_value="1124.0"/>
    <variable units="micromolar" name="K_BSL" initial_value="8.7"/>
    <variable units="cm_per_millisecond" name="P_Ca" initial_value="1.51E-11"/>
    <variable units="flux" name="J_xfer"/>
    <variable units="millisecond" name="tau_xfer" initial_value="0.0007"/>
    <variable units="micromolar" name="Ca_myo" initial_value="0.1"/>
    <variable units="dimensionless" name="i" initial_value="0.0"/>
    <variable units="micro_litre" name="V_myo" initial_value="0.1"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="joule_per_mole_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="Ca_o"/>
    <variable units="millimolar" public_interface="in" name="O"/>
    <variable units="millimolar" public_interface="in" name="O_Ca"/>
    <variable units="millimolar" public_interface="in" name="P_O3"/>
    <variable units="millimolar" public_interface="in" name="P_O4"/>
    <variable units="micromolar" public_interface="in" name="Ca_JSR"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="J_xfer_calculation">
        <eq/>
        <ci> J_xfer </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> tau_xfer </ci>
          </apply>
          <apply>
            <minus/>
            <ci> Ca_SS </ci>
            <ci> Ca_myo </ci>
          </apply>
        </apply>
      </apply>
     
      <apply id="J_DHPR_calculation">
        <eq/>
        <ci> J_DHPR </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> I_DHPR </ci>
            <ci> DHPR_open </ci>
          </apply>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> F </ci>
            <ci> V_SS </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="I_DHPR_calculation">
        <eq/>
        <ci> I_DHPR </ci>
        <apply>
          <times/>
          <ci> P_Ca </ci>
          <cn cellml:units="dimensionless"> 4.0 </cn>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> V </ci>
              <apply>
                <power/>
                <ci> F </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </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>
                    <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>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.341 </cn>
                <ci> Ca_o </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <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>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="DHPR_open_calculation">
        <eq/>
        <ci> DHPR_open </ci>
        <piecewise>
          <piece>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <or/>
              <apply>
                <eq/>
                <ci> i </ci>
                <ci> O </ci>
              </apply>
              <apply>
                <eq/>
                <ci> i </ci>
                <ci> O_Ca </ci>
              </apply> 
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="dimensionless"> 0.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
      
      <apply id="J_RyR_calculation">
        <eq/>
        <ci> J_RyR </ci>
        <apply>
          <sum/>
          <bvar>
            <ci> i </ci>
          </bvar>
          <lowlimit>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </lowlimit>
          <uplimit>
            <cn cellml:units="dimensionless"> 8.0 </cn>
          </uplimit>
          <apply>
            <times/>
            <ci> J_RyR_perm </ci>
            <ci> RyR_open </ci>
            <apply>
              <minus/>
              <ci> Ca_JSR </ci>
              <ci> Ca_SS </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="RyR_open_calculation">
        <eq/>
        <ci> RyR_open </ci>
        <piecewise>
          <piece>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <or/>
              <apply>
                <eq/>
                <ci> i </ci>
                <ci> P_O3 </ci>
              </apply>
              <apply>
                <eq/>
                <ci> i </ci>
                <ci> P_O4 </ci>
              </apply> 
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="dimensionless"> 0.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
      
      <apply id="calcium_subspace_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_SS </ci>
        </apply>
        <apply>
          <times/>
          <ci> beta_SS </ci>
          <apply>
            <plus/>
            <ci> J_DHPR </ci>
            <ci> J_RyR </ci>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> V_myo </ci>
                <ci> V_SS </ci>
              </apply>
              <ci> J_xfer </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_SS_calculation">
        <eq/>
        <ci> beta_SS </ci>
        <apply>
          <power/>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> B_SR </ci>
                <ci> K_BSR </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <plus/>
                  <ci> K_BSR </ci>
                  <ci> Ca_SS </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> B_SL </ci>
                <ci> K_BSL </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <plus/>
                  <ci> K_BSL </ci>
                  <ci> Ca_SS </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless"> -1.0 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sarcoplasmic_reticulum">
    <variable units="micromolar" public_interface="out" name="Ca_JSR"/>
    
    <variable units="micro_litre" name="V_JSR" initial_value="2.03E-12"/>
    <variable units="millimolar" name="beta_JSR"/>
    <variable units="millimolar" name="CSQN_tot" initial_value="15.0"/>
    <variable units="millimolar" name="K_CSQN" initial_value="0.8"/>
    <variable units="flux" name="J_tr"/>
    <variable units="millisecond" name="tau_tr" initial_value="0.005"/>
    <variable units="micromolar" name="Ca_NSR" initial_value="800.0"/>
      
    <variable units="flux" public_interface="in" name="J_RyR"/>
    <variable units="micro_litre" public_interface="in" name="V_SS"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="calcium_JSR_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_JSR </ci>
        </apply>
        <apply>
          <times/>
          <ci> beta_JSR </ci>
          <apply>
            <minus/>
            <ci> J_tr </ci>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> V_SS </ci>
                <ci> V_JSR </ci>
              </apply>
              <ci> J_RyR </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_JSR_calculation">
        <eq/>
        <ci> beta_JSR </ci>
        <apply>
          <power/>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <apply>
                <times/>
                <ci> CSQN_tot </ci>
                <ci> K_CSQN </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <plus/>
                  <ci> K_CSQN </ci>
                  <ci> Ca_JSR </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless"> -1.0 </cn>
        </apply>
      </apply>
      
      <apply id="J_tr_calculation">
        <eq/>
        <ci> J_tr </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> tau_tr </ci>
          </apply>
          <apply>
            <minus/>
            <ci> Ca_NSR </ci>
            <ci> Ca_JSR </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
 
  <!--
    The following components describe all the channel states involved in the 
    model. 
  -->
   
   <component cmeta:id="C0" name="C0">
      <variable units="micromolar" public_interface="out" name="C0" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C0_rxn10"/>
      <variable units="flux" public_interface="in" name="delta_C0_rxn3"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C0</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C0_rxn10</ci>
               <ci>delta_C0_rxn3</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component cmeta:id="C1" name="C1">
      <variable units="micromolar" public_interface="out" name="C1" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C1_rxn3"/>
      <variable units="flux" public_interface="in" name="delta_C1_rxn11"/>
      <variable units="flux" public_interface="in" name="delta_C1_rxn2"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C1</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C1_rxn3</ci>
               <ci>delta_C1_rxn11</ci>
               <ci>delta_C1_rxn2</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component cmeta:id="C3" name="C3">
      <variable units="micromolar" public_interface="out" name="C3" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C3_rxn4"/>
      <variable units="flux" public_interface="in" name="delta_C3_rxn13"/>
      <variable units="flux" public_interface="in" name="delta_C3_rxn1"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C3</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C3_rxn4</ci>
               <ci>delta_C3_rxn13</ci>
               <ci>delta_C3_rxn1</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component cmeta:id="C4" name="C4">
      <variable units="micromolar" public_interface="out" name="C4" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C4_rxn1"/>
      <variable units="flux" public_interface="in" name="delta_C4_rxn14"/>
      <variable units="flux" public_interface="in" name="delta_C4_rxn0"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C4</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C4_rxn1</ci>
               <ci>delta_C4_rxn14</ci>
               <ci>delta_C4_rxn0</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component cmeta:id="O" name="O">
      <variable units="micromolar" public_interface="out" name="O" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_O_rxn0"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O</ci>
            </apply>
            <ci>delta_O_rxn0</ci>
         </apply>
      </math>
   </component>
  
   <component cmeta:id="C_Ca0" name="C_Ca0">
      <variable units="micromolar" public_interface="out" name="C_Ca0" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca0_rxn9"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca0_rxn10"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_Ca0</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C_Ca0_rxn9</ci>
               <ci>delta_C_Ca0_rxn10</ci>
            </apply>
         </apply>
      </math>
   </component>
  
   <component cmeta:id="C_Ca1" name="C_Ca1">
      <variable units="micromolar" public_interface="out" name="C_Ca1" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca1_rxn9"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca1_rxn8"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca1_rxn11"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_Ca1</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C_Ca1_rxn9</ci>
               <ci>delta_C_Ca1_rxn8</ci>
               <ci>delta_C_Ca1_rxn11</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component cmeta:id="C_Ca2" name="C_Ca2">
      <variable units="micromolar" public_interface="out" name="C_Ca2" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca2_rxn8"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca2_rxn7"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca2_rxn12"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_Ca2</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C_Ca2_rxn8</ci>
               <ci>delta_C_Ca2_rxn7</ci>
               <ci>delta_C_Ca2_rxn12</ci>
            </apply>
         </apply>
      </math>
   </component>
  
   <component cmeta:id="C_Ca3" name="C_Ca3">
      <variable units="micromolar" public_interface="out" name="C_Ca3" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca3_rxn7"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca3_rxn6"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca3_rxn13"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_Ca3</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C_Ca3_rxn7</ci>
               <ci>delta_C_Ca3_rxn6</ci>
               <ci>delta_C_Ca3_rxn13</ci>
            </apply>
         </apply>
      </math>
   </component>
  
   <component cmeta:id="C_Ca4" name="C_Ca4">
      <variable units="micromolar" public_interface="out" name="C_Ca4" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca4_rxn6"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca4_rxn5"/>
      <variable units="flux" public_interface="in" name="delta_C_Ca4_rxn14"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C_Ca4</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C_Ca4_rxn6</ci>
               <ci>delta_C_Ca4_rxn5</ci>
               <ci>delta_C_Ca4_rxn14</ci>
            </apply>
         </apply>
      </math>
   </component>
  
   <component cmeta:id="O_Ca" name="O_Ca">
      <variable units="micromolar" public_interface="out" name="O_Ca" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_O_Ca_rxn5"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>O_Ca</ci>
            </apply>
            <ci>delta_O_Ca_rxn5</ci>
         </apply>
      </math>
   </component>
  
   <component cmeta:id="C2" name="C2">
      <variable units="micromolar" public_interface="out" name="C2" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_C2_rxn2"/>
      <variable units="flux" public_interface="in" name="delta_C2_rxn12"/>
      <variable units="flux" public_interface="in" name="delta_C2_rxn4"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C2</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_C2_rxn2</ci>
               <ci>delta_C2_rxn12</ci>
               <ci>delta_C2_rxn4</ci>
            </apply>
         </apply>
      </math>
   </component>
  
   <component cmeta:id="P_O3" name="P_O3">
      <variable units="micromolar" public_interface="out" name="P_O3" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_P_O3_rxn16"/>
      <variable units="flux" public_interface="in" name="delta_P_O3_rxn17"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_O3</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_P_O3_rxn16</ci>
               <ci>delta_P_O3_rxn17</ci>
            </apply>
         </apply>
      </math>
   </component>
 
   <component cmeta:id="P_O4" name="P_O4">
      <variable units="micromolar" public_interface="out" name="P_O4" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_P_O4_rxn17"/>
      <variable units="flux" public_interface="in" name="delta_P_O4_rxn18"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_O4</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_P_O4_rxn17</ci>
               <ci>delta_P_O4_rxn18</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component cmeta:id="P_C2" name="P_C2">
      <variable units="micromolar" public_interface="out" name="P_C2" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_P_C2_rxn15"/>
      <variable units="flux" public_interface="in" name="delta_P_C2_rxn16"/>
      <variable units="flux" public_interface="in" name="delta_P_C2_rxn19"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_C2</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_P_C2_rxn15</ci>
               <ci>delta_P_C2_rxn16</ci>
               <ci>delta_P_C2_rxn19</ci>
            </apply>
         </apply>
      </math>
   </component>
  
   <component cmeta:id="P_C5" name="P_C5">
      <variable units="micromolar" public_interface="out" name="P_C5" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_P_C5_rxn18"/>
      <variable units="flux" public_interface="in" name="delta_P_C5_rxn19"/>
      <variable units="flux" public_interface="in" name="delta_P_C5_rxn20"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_C5</ci>
            </apply>
            <apply>
               <plus/>
               <ci>delta_P_C5_rxn18</ci>
               <ci>delta_P_C5_rxn19</ci>
               <ci>delta_P_C5_rxn20</ci>
            </apply>
         </apply>
      </math>
   </component>
 
   <component cmeta:id="P_C1" name="P_C1">
      <variable units="micromolar" public_interface="out" name="P_C1" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_P_C1_rxn15"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_C1</ci>
            </apply>
            <ci>delta_P_C1_rxn15</ci>
         </apply>
      </math>
   </component>
   
   <component cmeta:id="P_C6" name="P_C6">
      <variable units="micromolar" public_interface="out" name="P_C6" initial_value="1.0"/>
      <variable units="flux" public_interface="in" name="delta_P_C6_rxn20"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_C6</ci>
            </apply>
            <ci>delta_P_C6_rxn20</ci>
         </apply>
      </math>
   </component>
   
   <!--
    The following components describe the channel state changes of the model.
   -->
   <component name="reaction_rate_constants">
     <variable units="first_order_rate_constant" public_interface="out" name="alpha" initial_value="1.0"/>
     <variable units="first_order_rate_constant" public_interface="out" name="gamma" initial_value="1.0"/>
     <variable units="first_order_rate_constant" public_interface="out" name="beta" initial_value="1.0"/>
     <variable units="first_order_rate_constant" public_interface="out" name="alpha_Ca" initial_value="1.0"/>
     <variable units="first_order_rate_constant" public_interface="out" name="beta_Ca" initial_value="1.0"/>
     <variable units="first_order_rate_constant" public_interface="out" name="omega" initial_value="10.0"/>
     <variable units="first_order_rate_constant" public_interface="out" name="a" initial_value="2.0"/>
     <variable units="first_order_rate_constant" public_interface="out" name="b" initial_value="2.0"/>
               
     <variable units="millivolt" public_interface="in" name="V"/>
     <variable units="micromolar" public_interface="in" name="Ca_SS"/>
    
     <math xmlns="http://www.w3.org/1998/Math/MathML">
       <apply id="gamma_calculation">
        <eq/>
         <ci> gamma </ci>
         <apply>
          <times/>
           <ci> Ca_SS </ci>
           <cn cellml:units="second_order_rate_constant"> 1.0 </cn>
         </apply>
       </apply>
      
      <apply id="alpha_calculation">
        <eq/>
        <ci> alpha </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 0.4 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <plus/>
                <ci> V </ci>
                <cn cellml:units="millivolt"> 6.0 </cn>
              </apply>
              <cn cellml:units="millivolt"> 25.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_calculation">
        <eq/>
        <ci> beta </ci>
        <apply>
          <times/>
          <cn cellml:units="first_order_rate_constant"> 0.05 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 6.0 </cn>
                </apply>
              </apply>
              <cn cellml:units="millivolt"> 29.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="alpha_Ca_calculation">
        <eq/>
        <ci> alpha_Ca </ci>
        <apply>
          <times/>
          <ci> alpha </ci>
          <ci> a </ci>
        </apply>
      </apply>
      
      <apply id="beta_Ca_calculation">
        <eq/>
        <ci> beta_Ca </ci>
        <apply>
          <divide/>
          <ci> beta </ci>
          <ci> b </ci>
        </apply>
      </apply>
    </math>     
  </component> 
  
   <component name="reaction0">
      <variable units="micromolar" public_interface="in" name="C4"/>
      <variable units="micromolar" public_interface="in" name="O"/>
      <variable units="flux" public_interface="out" name="delta_C4_rxn0"/>
      <variable units="flux" public_interface="out" name="delta_O_rxn0"/>
      <variable units="first_order_rate_constant" name="k0" initial_value="300.0"/>
      <variable units="first_order_rate_constant" name="k0_" initial_value="2000.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C4">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C4_rxn0" role="reactant"/>
         </variable_ref>
         <variable_ref variable="O">
            <role stoichiometry="1" direction="forward" delta_variable="delta_O_rxn0" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k0</ci>
                           <ci>C4</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k0_</ci>
                              <ci>O</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
   </component>
   
   <component name="reaction1">
      <variable units="micromolar" public_interface="in" name="C3"/>
      <variable units="micromolar" public_interface="in" name="C4"/>
      <variable units="flux" public_interface="out" name="delta_C3_rxn1"/>
      <variable units="flux" public_interface="out" name="delta_C4_rxn1"/>
      <variable units="first_order_rate_constant" public_interface="in" name="alpha"/>
      <variable units="first_order_rate_constant" public_interface="in" name="beta"/>
      <variable units="first_order_rate_constant" name="k1" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k1_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C3">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C3_rxn1" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C4">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C4_rxn1" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k1</ci>
                           <ci>C3</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k1_</ci>
                              <ci>C4</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k1_calculation">
        <eq/>
          <ci> k1 </ci>
          <ci> alpha </ci>
        </apply>
        
        <apply id="k1__calculation">
        <eq/>
          <ci> k1_ </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 4.0 </cn>
            <ci> beta </ci>
          </apply>
        </apply>
      </math>
   </component>
  
   <component name="reaction2">
      <variable units="micromolar" public_interface="in" name="C1"/>
      <variable units="micromolar" public_interface="in" name="C2"/>
      <variable units="flux" public_interface="out" name="delta_C1_rxn2"/>
      <variable units="flux" public_interface="out" name="delta_C2_rxn2"/>
      <variable units="first_order_rate_constant" public_interface="in" name="alpha"/>
      <variable units="first_order_rate_constant" public_interface="in" name="beta"/>
      <variable units="first_order_rate_constant" name="k2" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k2_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C1">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C1_rxn2" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C2">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C2_rxn2" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k2</ci>
                           <ci>C1</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k2_</ci>
                              <ci>C2</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k2_calculation">
        <eq/>
          <ci> k2 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> alpha </ci>
          </apply>
        </apply>
        
        <apply id="k2__calculation">
        <eq/>
          <ci> k2_ </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> beta </ci>
          </apply>
        </apply>
      </math>
   </component>
 
   <component name="reaction3">
      <variable units="micromolar" public_interface="in" name="C0"/>
      <variable units="micromolar" public_interface="in" name="C1"/>
      <variable units="flux" public_interface="out" name="delta_C0_rxn3"/>
      <variable units="flux" public_interface="out" name="delta_C1_rxn3"/>
      <variable units="first_order_rate_constant" public_interface="in" name="alpha"/>
      <variable units="first_order_rate_constant" public_interface="in" name="beta"/>
      <variable units="first_order_rate_constant" name="k3" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k3_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C0">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C0_rxn3" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C1">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C1_rxn3" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k3</ci>
                           <ci>C0</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k3_</ci>
                              <ci>C1</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k3_calculation">
        <eq/>
          <ci> k3 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 4.0 </cn>
            <ci> alpha </ci>
          </apply>
        </apply>
        
        <apply id="k3__calculation">
        <eq/>
          <ci> k3_ </ci>
          <ci> beta </ci>
        </apply>
      </math>
   </component>
  
   <component name="reaction4">
      <variable units="micromolar" public_interface="in" name="C2"/>
      <variable units="micromolar" public_interface="in" name="C3"/>
      <variable units="flux" public_interface="out" name="delta_C2_rxn4"/>
      <variable units="flux" public_interface="out" name="delta_C3_rxn4"/>
      <variable units="first_order_rate_constant" public_interface="in" name="alpha"/>
      <variable units="first_order_rate_constant" public_interface="in" name="beta"/>
      <variable units="first_order_rate_constant" name="k4" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k4_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C2">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C2_rxn4" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C3">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C3_rxn4" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k4</ci>
                           <ci>C2</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k4_</ci>
                              <ci>C3</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k4_calculation">
        <eq/>
          <ci> k4 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> alpha </ci>
          </apply>
        </apply>
        
        <apply id="k4__calculation">
        <eq/>
          <ci> k4_ </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> beta </ci>
          </apply>
        </apply>
      </math>
   </component>
  
   <component name="reaction5">
      <variable units="micromolar" public_interface="in" name="C_Ca4"/>
      <variable units="micromolar" public_interface="in" name="O_Ca"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca4_rxn5"/>
      <variable units="flux" public_interface="out" name="delta_O_Ca_rxn5"/>
      <variable units="first_order_rate_constant" name="k5" initial_value="5.0"/>
      <variable units="first_order_rate_constant" name="k5_" initial_value="7000.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca4">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca4_rxn5" role="reactant"/>
         </variable_ref>
         <variable_ref variable="O_Ca">
            <role stoichiometry="1" direction="forward" delta_variable="delta_O_Ca_rxn5" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k5</ci>
                           <ci>C_Ca4</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k5_</ci>
                              <ci>O_Ca</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
   </component>
   
   <component name="reaction6">
      <variable units="micromolar" public_interface="in" name="C_Ca3"/>
      <variable units="micromolar" public_interface="in" name="C_Ca4"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca3_rxn6"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca4_rxn6"/>
      <variable units="first_order_rate_constant" public_interface="in" name="alpha_Ca"/>
      <variable units="first_order_rate_constant" public_interface="in" name="beta_Ca"/>
      <variable units="first_order_rate_constant" name="k6" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k6_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca3">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca3_rxn6" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C_Ca4">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca4_rxn6" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k6</ci>
                           <ci>C_Ca3</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k6_</ci>
                              <ci>C_Ca4</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k6_calculation">
        <eq/>
          <ci> k6 </ci>
          <ci> alpha_Ca </ci>
        </apply>
        
        <apply id="k6__calculation">
        <eq/>
          <ci> k6_ </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 4.0 </cn>
            <ci> beta_Ca </ci>
          </apply>
        </apply>
      </math>
   </component>
 
   <component name="reaction7">
      <variable units="micromolar" public_interface="in" name="C_Ca2"/>
      <variable units="micromolar" public_interface="in" name="C_Ca3"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca2_rxn7"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca3_rxn7"/>
      <variable units="first_order_rate_constant" public_interface="in" name="alpha_Ca"/>
      <variable units="first_order_rate_constant" public_interface="in" name="beta_Ca"/>
      <variable units="first_order_rate_constant" name="k7" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k7_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca2">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca2_rxn7" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C_Ca3">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca3_rxn7" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k7</ci>
                           <ci>C_Ca2</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k7_</ci>
                              <ci>C_Ca3</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k7_calculation">
        <eq/>
          <ci> k7 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> alpha_Ca </ci>
          </apply>
        </apply>
        
        <apply id="k7__calculation">
        <eq/>
          <ci> k7_ </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> beta_Ca </ci>
          </apply>
        </apply>
      </math>
   </component>
   
   <component name="reaction8">
      <variable units="micromolar" public_interface="in" name="C_Ca1"/>
      <variable units="micromolar" public_interface="in" name="C_Ca2"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca1_rxn8"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca2_rxn8"/>
      <variable units="first_order_rate_constant" public_interface="in" name="alpha_Ca"/>
      <variable units="first_order_rate_constant" public_interface="in" name="beta_Ca"/>
      <variable units="first_order_rate_constant" name="k8" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k8_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca1">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca1_rxn8" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C_Ca2">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca2_rxn8" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k8</ci>
                           <ci>C_Ca1</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k8_</ci>
                              <ci>C_Ca2</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k8_calculation">
        <eq/>
          <ci> k8 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> alpha_Ca </ci>
          </apply>
        </apply>
        
        <apply id="k8__calculation">
        <eq/>
          <ci> k8_ </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> beta_Ca </ci>
          </apply>
        </apply>
      </math>
   </component>
 
   <component name="reaction9">
      <variable units="micromolar" public_interface="in" name="C_Ca0"/>
      <variable units="micromolar" public_interface="in" name="C_Ca1"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca0_rxn9"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca1_rxn9"/>
      <variable units="first_order_rate_constant" public_interface="in" name="alpha_Ca"/>
      <variable units="first_order_rate_constant" public_interface="in" name="beta_Ca"/>
      <variable units="first_order_rate_constant" name="k9" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k9_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca0">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca0_rxn9" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C_Ca1">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca1_rxn9" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k9</ci>
                           <ci>C_Ca0</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k9_</ci>
                              <ci>C_Ca1</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k9_calculation">
        <eq/>
          <ci> k9 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 4.0 </cn>
            <ci> alpha_Ca </ci>
          </apply>
        </apply>
        
        <apply id="k9__calculation">
        <eq/>
          <ci> k9_ </ci>
          <ci> beta_Ca </ci>
        </apply>
      </math>
   </component>
   
   <component name="reaction10">
      <variable units="micromolar" public_interface="in" name="C_Ca0"/>
      <variable units="micromolar" public_interface="in" name="C0"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca0_rxn10"/>
      <variable units="flux" public_interface="out" name="delta_C0_rxn10"/>
      <variable units="first_order_rate_constant" public_interface="in" name="omega"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gamma"/>
      <variable units="first_order_rate_constant" name="k10" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k10_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca0">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca0_rxn10" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C0">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C0_rxn10" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k10</ci>
                           <ci>C_Ca0</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k10_</ci>
                              <ci>C0</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k10_calculation">
        <eq/>
          <ci> k10 </ci>
          <ci> omega </ci>
        </apply>
        
        <apply id="k10__calculation">
        <eq/>
          <ci> k10_ </ci>
          <ci> gamma </ci>
        </apply>
      </math>
   </component>
   
   <component name="reaction11">
      <variable units="micromolar" public_interface="in" name="C_Ca1"/>
      <variable units="micromolar" public_interface="in" name="C1"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca1_rxn11"/>
      <variable units="flux" public_interface="out" name="delta_C1_rxn11"/>
      <variable units="first_order_rate_constant" public_interface="in" name="omega"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gamma"/>
      <variable units="first_order_rate_constant" public_interface="in" name="a"/>
      <variable units="first_order_rate_constant" public_interface="in" name="b"/>
      <variable units="first_order_rate_constant" name="k11" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k11_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca1">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca1_rxn11" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C1">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C1_rxn11" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k11</ci>
                           <ci>C_Ca1</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k11_</ci>
                              <ci>C1</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k11_calculation">
        <eq/>
          <ci> k11 </ci>
          <apply>
          <divide/>
            <ci> omega </ci>
            <ci> b </ci>
          </apply>
        </apply>
        
        <apply id="k11__calculation">
        <eq/>
          <ci> k11_ </ci>
          <apply>
          <times/>
            <ci> gamma </ci>
            <ci> a </ci>
          </apply>
        </apply>
      </math>
   </component>
  
   <component name="reaction12">
      <variable units="micromolar" public_interface="in" name="C_Ca2"/>
      <variable units="micromolar" public_interface="in" name="C2"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca2_rxn12"/>
      <variable units="flux" public_interface="out" name="delta_C2_rxn12"/>
      <variable units="first_order_rate_constant" public_interface="in" name="omega"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gamma"/>
      <variable units="first_order_rate_constant" public_interface="in" name="a"/>
      <variable units="first_order_rate_constant" public_interface="in" name="b"/>
      <variable units="first_order_rate_constant" name="k12" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k12_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca2">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca2_rxn12" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C2">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C2_rxn12" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k12</ci>
                           <ci>C_Ca2</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k12_</ci>
                              <ci>C2</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k12_calculation">
        <eq/>
          <ci> k12 </ci>
          <apply>
          <divide/>
            <ci> omega </ci>
            <apply>
            <power/>
              <ci> b </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
        </apply>
        
        <apply id="k12__calculation">
        <eq/>
          <ci> k12_ </ci>
          <apply>
          <times/>
            <ci> gamma </ci>
            <apply>
            <power/>
              <ci> a </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
        </apply>
      </math>
   </component>
  
   <component name="reaction13">
      <variable units="micromolar" public_interface="in" name="C_Ca3"/>
      <variable units="micromolar" public_interface="in" name="C3"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca3_rxn13"/>
      <variable units="flux" public_interface="out" name="delta_C3_rxn13"/>
      <variable units="first_order_rate_constant" public_interface="in" name="omega"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gamma"/>
      <variable units="first_order_rate_constant" public_interface="in" name="a"/>
      <variable units="first_order_rate_constant" public_interface="in" name="b"/>
      <variable units="first_order_rate_constant" name="k13" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k13_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca3">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca3_rxn13" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C3">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C3_rxn13" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k13</ci>
                           <ci>C_Ca3</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k13_</ci>
                              <ci>C3</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k13_calculation">
        <eq/>
          <ci> k13 </ci>
          <apply>
          <divide/>
            <ci> omega </ci>
            <apply>
            <power/>
              <ci> b </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
        </apply>
        
        <apply id="k13__calculation">
        <eq/>
          <ci> k13_ </ci>
          <apply>
          <times/>
            <ci> gamma </ci>
            <apply>
            <power/>
              <ci> a </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
        </apply>
      </math>
   </component>
   
   <component name="reaction14">
      <variable units="micromolar" public_interface="in" name="C_Ca4"/>
      <variable units="micromolar" public_interface="in" name="C4"/>
      <variable units="flux" public_interface="out" name="delta_C_Ca4_rxn14"/>
      <variable units="flux" public_interface="out" name="delta_C4_rxn14"/>
      <variable units="first_order_rate_constant" public_interface="in" name="omega"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gamma"/>
      <variable units="first_order_rate_constant" public_interface="in" name="a"/>
      <variable units="first_order_rate_constant" public_interface="in" name="b"/>
      <variable units="first_order_rate_constant" name="k14" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k14_" initial_value="1.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="C_Ca4">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C_Ca4_rxn14" role="reactant"/>
         </variable_ref>
         <variable_ref variable="C4">
            <role stoichiometry="1" direction="forward" delta_variable="delta_C4_rxn14" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k14</ci>
                           <ci>C_Ca4</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k14_</ci>
                              <ci>C4</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k14_calculation">
        <eq/>
          <ci> k14 </ci>
          <apply>
          <divide/>
            <ci> omega </ci>
            <apply>
            <power/>
              <ci> b </ci>
              <cn cellml:units="dimensionless"> 4.0 </cn>
            </apply>
          </apply>
        </apply>
        
        <apply id="k14__calculation">
        <eq/>
          <ci> k14_ </ci>
          <apply>
          <times/>
            <ci> gamma </ci>
            <apply>
            <power/>
              <ci> a </ci>
              <cn cellml:units="dimensionless"> 4.0 </cn>
            </apply>
          </apply>
        </apply>
      </math>
   </component>
  
   <component name="reaction15">
      <variable units="micromolar" public_interface="in" name="P_C1"/>
      <variable units="micromolar" public_interface="in" name="P_C2"/>
      <variable units="flux" public_interface="out" name="delta_P_C1_rxn15"/>
      <variable units="flux" public_interface="out" name="delta_P_C2_rxn15"/>
      <variable units="micromolar" public_interface="in" name="Ca_SS"/>
      <variable units="first_order_rate_constant" name="k15" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k15_" initial_value="250000.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="P_C1">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_C1_rxn15" role="reactant"/>
         </variable_ref>
         <variable_ref variable="P_C2">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_C2_rxn15" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k15</ci>
                           <ci>P_C1</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k15_</ci>
                              <ci>P_C2</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k15_calculation">
        <eq/>
          <ci> k15 </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <cn cellml:units="first_order_rate_constant"> 3000000.0 </cn>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
            <apply>
            <plus/>
              <apply>
              <power/>
                <cn cellml:units="micromolar"> 2.68 </cn>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>
   </component>
   
   <component name="reaction16">
      <variable units="micromolar" public_interface="in" name="P_C2"/>
      <variable units="micromolar" public_interface="in" name="P_O3"/>
      <variable units="flux" public_interface="out" name="delta_P_C2_rxn16"/>
      <variable units="flux" public_interface="out" name="delta_P_O3_rxn16"/>
      <variable units="micromolar" public_interface="in" name="Ca_SS"/>
      <variable units="first_order_rate_constant" name="k16" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k16_" initial_value="9600.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="P_C2">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_C2_rxn16" role="reactant"/>
         </variable_ref>
         <variable_ref variable="P_O3">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_O3_rxn16" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k16</ci>
                           <ci>P_C2</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k16_</ci>
                              <ci>P_O3</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k16_calculation">
        <eq/>
          <ci> k16 </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <cn cellml:units="first_order_rate_constant"> 30000000.0 </cn>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
            <apply>
            <plus/>
              <apply>
              <power/>
                <cn cellml:units="micromolar"> 2.68 </cn>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>
   </component>
 
   <component name="reaction17">
      <variable units="micromolar" public_interface="in" name="P_O3"/>
      <variable units="micromolar" public_interface="in" name="P_O4"/>
      <variable units="flux" public_interface="out" name="delta_P_O3_rxn17"/>
      <variable units="flux" public_interface="out" name="delta_P_O4_rxn17"/>
      <variable units="micromolar" public_interface="in" name="Ca_SS"/>
      <variable units="first_order_rate_constant" name="k17" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k17_" initial_value="13000.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="P_O3">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_O3_rxn17" role="reactant"/>
         </variable_ref>
         <variable_ref variable="P_O4">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_O4_rxn17" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k17</ci>
                           <ci>P_O3</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k17_</ci>
                              <ci>P_O4</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k17_calculation">
        <eq/>
          <ci> k17 </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <cn cellml:units="first_order_rate_constant"> 3000000.0 </cn>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
            <apply>
            <plus/>
              <apply>
              <power/>
                <cn cellml:units="micromolar"> 2.68 </cn>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>
   </component>
   
   <component name="reaction18">
      <variable units="micromolar" public_interface="in" name="P_O4"/>
      <variable units="micromolar" public_interface="in" name="P_C5"/>
      <variable units="flux" public_interface="out" name="delta_P_O4_rxn18"/>
      <variable units="flux" public_interface="out" name="delta_P_C5_rxn18"/>
      <variable units="micromolar" public_interface="in" name="Ca_SS"/>
      <variable units="first_order_rate_constant" name="k18" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k18_" initial_value="66.67"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="P_O4">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_O4_rxn18" role="reactant"/>
         </variable_ref>
         <variable_ref variable="P_C5">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_C5_rxn18" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k18</ci>
                           <ci>P_O4</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k18_</ci>
                              <ci>P_C5</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k18_calculation">
        <eq/>
          <ci> k18 </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <cn cellml:units="first_order_rate_constant"> 198.0 </cn>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
            <apply>
            <plus/>
              <apply>
              <power/>
                <cn cellml:units="micromolar"> 2.68 </cn>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>
   </component>
  
   <component name="reaction19">
      <variable units="micromolar" public_interface="in" name="P_C2"/>
      <variable units="micromolar" public_interface="in" name="P_C5"/>
      <variable units="flux" public_interface="out" name="delta_P_C2_rxn19"/>
      <variable units="flux" public_interface="out" name="delta_P_C5_rxn19"/>
      <variable units="micromolar" public_interface="in" name="Ca_SS"/>
      <variable units="first_order_rate_constant" name="k19" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k19_" initial_value="1.235"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="P_C2">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_C2_rxn19" role="reactant"/>
         </variable_ref>
         <variable_ref variable="P_C5">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_C5_rxn19" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k19</ci>
                           <ci>P_C2</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k19_</ci>
                              <ci>P_C5</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k19_calculation">
        <eq/>
          <ci> k19 </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <cn cellml:units="first_order_rate_constant"> 300000.0 </cn>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
            <apply>
            <plus/>
              <apply>
              <power/>
                <cn cellml:units="micromolar"> 2.68 </cn>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>
   </component>
  
   <component name="reaction20">
      <variable units="micromolar" public_interface="in" name="P_C5"/>
      <variable units="micromolar" public_interface="in" name="P_C6"/>
      <variable units="flux" public_interface="out" name="delta_P_C5_rxn20"/>
      <variable units="flux" public_interface="out" name="delta_P_C6_rxn20"/>
      <variable units="micromolar" public_interface="in" name="Ca_SS"/>
      <variable units="first_order_rate_constant" name="k20" initial_value="1.0"/>
      <variable units="first_order_rate_constant" name="k20_" initial_value="3000000.0"/>
      <variable units="flux" name="rate"/>
      <reaction reversible="yes">
         <variable_ref variable="P_C5">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_C5_rxn20" role="reactant"/>
         </variable_ref>
         <variable_ref variable="P_C6">
            <role stoichiometry="1" direction="forward" delta_variable="delta_P_C6_rxn20" role="product"/>
         </variable_ref>
         <variable_ref variable="rate">
            <role role="rate">
               <math xmlns="http://www.w3.org/1998/Math/MathML">
                  <apply>
                     <eq/>
                     <ci>rate</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>k20</ci>
                           <ci>P_C5</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k20_</ci>
                              <ci>P_C6</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </math>
            </role>
         </variable_ref>
      </reaction>
      
      <math xmlns="http://www.w3.org/1998/Math/MathML">
        <apply id="k20_calculation">
        <eq/>
          <ci> k20 </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <cn cellml:units="first_order_rate_constant"> 3000000.0 </cn>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
            <apply>
            <plus/>
              <apply>
              <power/>
                <cn cellml:units="micromolar"> 2.68 </cn>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>
   </component>
   
   <connection>
      <map_components component_2="environment" component_1="subspace"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   
   <connection>
      <map_components component_2="environment" component_1="sarcoplasmic_reticulum"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   
   <connection>
    <map_components component_2="sarcolemma" component_1="subspace"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="Ca_o" variable_1="Ca_o"/>
  </connection>
  
  <connection>
    <map_components component_2="sarcoplasmic_reticulum" component_1="subspace"/>
    <map_variables variable_2="J_RyR" variable_1="J_RyR"/>
    <map_variables variable_2="V_SS" variable_1="V_SS"/>
    <map_variables variable_2="Ca_JSR" variable_1="Ca_JSR"/>
  </connection>
  
   <connection>
      <map_components component_2="reaction10" component_1="C0"/>
      <map_variables variable_2="C0" variable_1="C0"/>
      <map_variables variable_2="delta_C0_rxn10" variable_1="delta_C0_rxn10"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction3" component_1="C0"/>
      <map_variables variable_2="C0" variable_1="C0"/>
      <map_variables variable_2="delta_C0_rxn3" variable_1="delta_C0_rxn3"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction3" component_1="C1"/>
      <map_variables variable_2="C1" variable_1="C1"/>
      <map_variables variable_2="delta_C1_rxn3" variable_1="delta_C1_rxn3"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction11" component_1="C1"/>
      <map_variables variable_2="C1" variable_1="C1"/>
      <map_variables variable_2="delta_C1_rxn11" variable_1="delta_C1_rxn11"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction2" component_1="C1"/>
      <map_variables variable_2="C1" variable_1="C1"/>
      <map_variables variable_2="delta_C1_rxn2" variable_1="delta_C1_rxn2"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction4" component_1="C3"/>
      <map_variables variable_2="C3" variable_1="C3"/>
      <map_variables variable_2="delta_C3_rxn4" variable_1="delta_C3_rxn4"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction13" component_1="C3"/>
      <map_variables variable_2="C3" variable_1="C3"/>
      <map_variables variable_2="delta_C3_rxn13" variable_1="delta_C3_rxn13"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction1" component_1="C3"/>
      <map_variables variable_2="C3" variable_1="C3"/>
      <map_variables variable_2="delta_C3_rxn1" variable_1="delta_C3_rxn1"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction1" component_1="C4"/>
      <map_variables variable_2="C4" variable_1="C4"/>
      <map_variables variable_2="delta_C4_rxn1" variable_1="delta_C4_rxn1"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction14" component_1="C4"/>
      <map_variables variable_2="C4" variable_1="C4"/>
      <map_variables variable_2="delta_C4_rxn14" variable_1="delta_C4_rxn14"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction0" component_1="C4"/>
      <map_variables variable_2="C4" variable_1="C4"/>
      <map_variables variable_2="delta_C4_rxn0" variable_1="delta_C4_rxn0"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction0" component_1="O"/>
      <map_variables variable_2="O" variable_1="O"/>
      <map_variables variable_2="delta_O_rxn0" variable_1="delta_O_rxn0"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction9" component_1="C_Ca0"/>
      <map_variables variable_2="C_Ca0" variable_1="C_Ca0"/>
      <map_variables variable_2="delta_C_Ca0_rxn9" variable_1="delta_C_Ca0_rxn9"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction10" component_1="C_Ca0"/>
      <map_variables variable_2="C_Ca0" variable_1="C_Ca0"/>
      <map_variables variable_2="delta_C_Ca0_rxn10" variable_1="delta_C_Ca0_rxn10"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction9" component_1="C_Ca1"/>
      <map_variables variable_2="C_Ca1" variable_1="C_Ca1"/>
      <map_variables variable_2="delta_C_Ca1_rxn9" variable_1="delta_C_Ca1_rxn9"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction8" component_1="C_Ca1"/>
      <map_variables variable_2="C_Ca1" variable_1="C_Ca1"/>
      <map_variables variable_2="delta_C_Ca1_rxn8" variable_1="delta_C_Ca1_rxn8"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction11" component_1="C_Ca1"/>
      <map_variables variable_2="C_Ca1" variable_1="C_Ca1"/>
      <map_variables variable_2="delta_C_Ca1_rxn11" variable_1="delta_C_Ca1_rxn11"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction8" component_1="C_Ca2"/>
      <map_variables variable_2="C_Ca2" variable_1="C_Ca2"/>
      <map_variables variable_2="delta_C_Ca2_rxn8" variable_1="delta_C_Ca2_rxn8"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction7" component_1="C_Ca2"/>
      <map_variables variable_2="C_Ca2" variable_1="C_Ca2"/>
      <map_variables variable_2="delta_C_Ca2_rxn7" variable_1="delta_C_Ca2_rxn7"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction12" component_1="C_Ca2"/>
      <map_variables variable_2="C_Ca2" variable_1="C_Ca2"/>
      <map_variables variable_2="delta_C_Ca2_rxn12" variable_1="delta_C_Ca2_rxn12"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction7" component_1="C_Ca3"/>
      <map_variables variable_2="C_Ca3" variable_1="C_Ca3"/>
      <map_variables variable_2="delta_C_Ca3_rxn7" variable_1="delta_C_Ca3_rxn7"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction6" component_1="C_Ca3"/>
      <map_variables variable_2="C_Ca3" variable_1="C_Ca3"/>
      <map_variables variable_2="delta_C_Ca3_rxn6" variable_1="delta_C_Ca3_rxn6"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction13" component_1="C_Ca3"/>
      <map_variables variable_2="C_Ca3" variable_1="C_Ca3"/>
      <map_variables variable_2="delta_C_Ca3_rxn13" variable_1="delta_C_Ca3_rxn13"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction6" component_1="C_Ca4"/>
      <map_variables variable_2="C_Ca4" variable_1="C_Ca4"/>
      <map_variables variable_2="delta_C_Ca4_rxn6" variable_1="delta_C_Ca4_rxn6"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction5" component_1="C_Ca4"/>
      <map_variables variable_2="C_Ca4" variable_1="C_Ca4"/>
      <map_variables variable_2="delta_C_Ca4_rxn5" variable_1="delta_C_Ca4_rxn5"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction14" component_1="C_Ca4"/>
      <map_variables variable_2="C_Ca4" variable_1="C_Ca4"/>
      <map_variables variable_2="delta_C_Ca4_rxn14" variable_1="delta_C_Ca4_rxn14"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction5" component_1="O_Ca"/>
      <map_variables variable_2="O_Ca" variable_1="O_Ca"/>
      <map_variables variable_2="delta_O_Ca_rxn5" variable_1="delta_O_Ca_rxn5"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction2" component_1="C2"/>
      <map_variables variable_2="C2" variable_1="C2"/>
      <map_variables variable_2="delta_C2_rxn2" variable_1="delta_C2_rxn2"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction12" component_1="C2"/>
      <map_variables variable_2="C2" variable_1="C2"/>
      <map_variables variable_2="delta_C2_rxn12" variable_1="delta_C2_rxn12"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction4" component_1="C2"/>
      <map_variables variable_2="C2" variable_1="C2"/>
      <map_variables variable_2="delta_C2_rxn4" variable_1="delta_C2_rxn4"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction16" component_1="P_O3"/>
      <map_variables variable_2="P_O3" variable_1="P_O3"/>
      <map_variables variable_2="delta_P_O3_rxn16" variable_1="delta_P_O3_rxn16"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction17" component_1="P_O3"/>
      <map_variables variable_2="P_O3" variable_1="P_O3"/>
      <map_variables variable_2="delta_P_O3_rxn17" variable_1="delta_P_O3_rxn17"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction17" component_1="P_O4"/>
      <map_variables variable_2="P_O4" variable_1="P_O4"/>
      <map_variables variable_2="delta_P_O4_rxn17" variable_1="delta_P_O4_rxn17"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction18" component_1="P_O4"/>
      <map_variables variable_2="P_O4" variable_1="P_O4"/>
      <map_variables variable_2="delta_P_O4_rxn18" variable_1="delta_P_O4_rxn18"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction15" component_1="P_C2"/>
      <map_variables variable_2="P_C2" variable_1="P_C2"/>
      <map_variables variable_2="delta_P_C2_rxn15" variable_1="delta_P_C2_rxn15"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction16" component_1="P_C2"/>
      <map_variables variable_2="P_C2" variable_1="P_C2"/>
      <map_variables variable_2="delta_P_C2_rxn16" variable_1="delta_P_C2_rxn16"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction19" component_1="P_C2"/>
      <map_variables variable_2="P_C2" variable_1="P_C2"/>
      <map_variables variable_2="delta_P_C2_rxn19" variable_1="delta_P_C2_rxn19"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction18" component_1="P_C5"/>
      <map_variables variable_2="P_C5" variable_1="P_C5"/>
      <map_variables variable_2="delta_P_C5_rxn18" variable_1="delta_P_C5_rxn18"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction19" component_1="P_C5"/>
      <map_variables variable_2="P_C5" variable_1="P_C5"/>
      <map_variables variable_2="delta_P_C5_rxn19" variable_1="delta_P_C5_rxn19"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction20" component_1="P_C5"/>
      <map_variables variable_2="P_C5" variable_1="P_C5"/>
      <map_variables variable_2="delta_P_C5_rxn20" variable_1="delta_P_C5_rxn20"/>
   </connection>
  
   <connection>
      <map_components component_2="reaction15" component_1="P_C1"/>
      <map_variables variable_2="P_C1" variable_1="P_C1"/>
      <map_variables variable_2="delta_P_C1_rxn15" variable_1="delta_P_C1_rxn15"/>
   </connection>
 
   <connection>
      <map_components component_2="reaction20" component_1="P_C6"/>
      <map_variables variable_2="P_C6" variable_1="P_C6"/>
      <map_variables variable_2="delta_P_C6_rxn20" variable_1="delta_P_C6_rxn20"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="C0"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="C1"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="C3"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="C4"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="O"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="C_Ca0"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="C_Ca1"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="C_Ca2"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="C_Ca3"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="C_Ca4"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="O_Ca"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="C2"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="P_O3"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="P_O4"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="P_C2"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="P_C5"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
 
   <connection>
      <map_components component_2="environment" component_1="P_C1"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="P_C6"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   
   <connection>
      <map_components component_2="sarcolemma" component_1="reaction_rate_constants"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="reaction_rate_constants"/>
      <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction1" component_1="reaction_rate_constants"/>
      <map_variables variable_2="alpha" variable_1="alpha"/>
      <map_variables variable_2="beta" variable_1="beta"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction2" component_1="reaction_rate_constants"/>
      <map_variables variable_2="alpha" variable_1="alpha"/>
      <map_variables variable_2="beta" variable_1="beta"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction3" component_1="reaction_rate_constants"/>
      <map_variables variable_2="alpha" variable_1="alpha"/>
      <map_variables variable_2="beta" variable_1="beta"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction4" component_1="reaction_rate_constants"/>
      <map_variables variable_2="alpha" variable_1="alpha"/>
      <map_variables variable_2="beta" variable_1="beta"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction6" component_1="reaction_rate_constants"/>
      <map_variables variable_2="alpha_Ca" variable_1="alpha_Ca"/>
      <map_variables variable_2="beta_Ca" variable_1="beta_Ca"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction7" component_1="reaction_rate_constants"/>
      <map_variables variable_2="alpha_Ca" variable_1="alpha_Ca"/>
      <map_variables variable_2="beta_Ca" variable_1="beta_Ca"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction8" component_1="reaction_rate_constants"/>
      <map_variables variable_2="alpha_Ca" variable_1="alpha_Ca"/>
      <map_variables variable_2="beta_Ca" variable_1="beta_Ca"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction9" component_1="reaction_rate_constants"/>
      <map_variables variable_2="alpha_Ca" variable_1="alpha_Ca"/>
      <map_variables variable_2="beta_Ca" variable_1="beta_Ca"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction10" component_1="reaction_rate_constants"/>
      <map_variables variable_2="omega" variable_1="omega"/>
      <map_variables variable_2="gamma" variable_1="gamma"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction11" component_1="reaction_rate_constants"/>
      <map_variables variable_2="omega" variable_1="omega"/>
      <map_variables variable_2="gamma" variable_1="gamma"/>
      <map_variables variable_2="a" variable_1="a"/>
      <map_variables variable_2="b" variable_1="b"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction12" component_1="reaction_rate_constants"/>
      <map_variables variable_2="omega" variable_1="omega"/>
      <map_variables variable_2="gamma" variable_1="gamma"/>
      <map_variables variable_2="a" variable_1="a"/>
      <map_variables variable_2="b" variable_1="b"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction13" component_1="reaction_rate_constants"/>
      <map_variables variable_2="omega" variable_1="omega"/>
      <map_variables variable_2="gamma" variable_1="gamma"/>
      <map_variables variable_2="a" variable_1="a"/>
      <map_variables variable_2="b" variable_1="b"/>
   </connection>
   
   <connection>
      <map_components component_2="reaction14" component_1="reaction_rate_constants"/>
      <map_variables variable_2="omega" variable_1="omega"/>
      <map_variables variable_2="gamma" variable_1="gamma"/>
      <map_variables variable_2="a" variable_1="a"/>
      <map_variables variable_2="b" variable_1="b"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="reaction15"/>
      <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="reaction16"/>
      <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="reaction17"/>
      <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="reaction18"/>
      <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="reaction19"/>
      <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="reaction20"/>
      <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="O"/>
      <map_variables variable_2="O" variable_1="O"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="O_Ca"/>
      <map_variables variable_2="O_Ca" variable_1="O_Ca"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="P_O3"/>
      <map_variables variable_2="P_O3" variable_1="P_O3"/>
   </connection>
   
   <connection>
      <map_components component_2="subspace" component_1="P_O4"/>
      <map_variables variable_2="P_O4" variable_1="P_O4"/>
   </connection>



<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#ec09e0d2-a605-4d35-9afa-d0c5302c2d43">
    <rdf:li>Cardiac Myocyte</rdf:li>
    <rdf:li>calcium dynamics</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#8aed1500-d219-4d8f-9cfa-d9ad9f15a643">
    <rdf:li rdf:resource="rdf:#6c5ee79c-adbe-45cd-8de8-c736477ac9d7"/>
    <rdf:li rdf:resource="rdf:#0361025a-e791-4ce8-b080-859fc0f0039b"/>
    <rdf:li rdf:resource="rdf:#4514839d-e74e-4434-8c92-fb716dd38414"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#2be31a85-b778-4ef7-8bfb-1e92a7f2fc0b">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6c5ee79c-adbe-45cd-8de8-c736477ac9d7">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#faa5ede7-e2bd-48d2-8a49-6c4d5f8220fc"/>
  </rdf:Description>
  <rdf:Description rdf:about="#rice_winslow_jafri_1999_version01">
    <dc:title>
        The Rice et al. 1999 model of the functional Ca2+ release unit. 
      </dc:title>
    <cmeta:bio_entity>Cardiac Myocyte</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#032205d2-1b7f-4cec-aaf2-527096eaffa8"/>
    <bqs:reference rdf:resource="rdf:#7c4fe4b5-a46f-40e4-a16d-dc5fc8e36755"/>
    <bqs:reference rdf:resource="rdf:#93ec15fe-b7da-4b05-bde9-35ab06c34130"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#499d0c2d-3409-4804-a735-b09139722b44">
    <vCard:N rdf:resource="rdf:#5224a79e-8a89-4f61-8193-78cba8e73246"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a22e65f7-7fe0-45f7-aab8-6be4b3054c37">
    <dcterms:W3CDTF>2002-06-29</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cd68b687-4af5-4fae-a412-dcdbfd59cd14">
    <dc:creator rdf:resource="rdf:#8aed1500-d219-4d8f-9cfa-d9ad9f15a643"/>
    <dc:title>
            Modeling Gain and Gradedness of Ca2+ Release in the Functional Unit              of the Cardiac Diadic Space
          </dc:title>
    <bqs:volume>77</bqs:volume>
    <bqs:first_page>1871</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#53695025-8bd8-4f0b-ae20-1d62522f0e02"/>
    <dcterms:issued rdf:resource="rdf:#4fc02258-4d9a-4fdb-9d64-fe354cb51f39"/>
    <bqs:last_page>1884</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#53695025-8bd8-4f0b-ae20-1d62522f0e02">
    <dc:title>Biophysical Journal</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#19be5ad1-95e7-41b0-8060-70a5c2976011">
    <vCard:Given>Raimond</vCard:Given>
    <vCard:Family>Winslow</vCard:Family>
    <vCard:Other>L</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c7a0ef07-f123-4b77-88e2-e2df1a240b0a">
    <dcterms:modified rdf:resource="rdf:#4dac4b00-40a6-4882-b3e4-303d5273edbd"/>
    <rdf:value>
          Fixed "J_RyR_calculation".
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#499d0c2d-3409-4804-a735-b09139722b44"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#032205d2-1b7f-4cec-aaf2-527096eaffa8">
    <dc:creator rdf:resource="rdf:#2be31a85-b778-4ef7-8bfb-1e92a7f2fc0b"/>
    <rdf:value>
          This is the CellML description of Rice et al's model of the functional           Ca2+ release unit.
        </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4514839d-e74e-4434-8c92-fb716dd38414">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#4df286f0-2e10-40e9-bc52-6fe9b37092ff"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f4d83975-8097-4aba-80ac-6d6961d5b952">
    <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:#21a3ef0f-bdda-442d-abad-f3a19b0f54f4">
    <vCard:ORG rdf:resource="rdf:#73b45ddd-5b85-44c2-9ec6-96df914851b0"/>
    <vCard:EMAIL rdf:resource="rdf:#f4d83975-8097-4aba-80ac-6d6961d5b952"/>
    <vCard:N rdf:resource="rdf:#5c2df81a-5274-44f2-ac26-ff815319aece"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#93ec15fe-b7da-4b05-bde9-35ab06c34130">
    <bqs:Pubmed_id>10512809</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#cd68b687-4af5-4fae-a412-dcdbfd59cd14"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0361025a-e791-4ce8-b080-859fc0f0039b">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#19be5ad1-95e7-41b0-8060-70a5c2976011"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4fc02258-4d9a-4fdb-9d64-fe354cb51f39">
    <dcterms:W3CDTF>1999-10</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Bioengineering Institute
      </dc:publisher>
    <cmeta:modification rdf:resource="rdf:#c7a0ef07-f123-4b77-88e2-e2df1a240b0a"/>
    <dcterms:created rdf:resource="rdf:#a22e65f7-7fe0-45f7-aab8-6be4b3054c37"/>
    <dc:creator rdf:resource="rdf:#21a3ef0f-bdda-442d-abad-f3a19b0f54f4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4df286f0-2e10-40e9-bc52-6fe9b37092ff">
    <vCard:Given>M</vCard:Given>
    <vCard:Family>Jafri</vCard:Family>
    <vCard:Other>Saleet</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7c4fe4b5-a46f-40e4-a16d-dc5fc8e36755">
    <dc:subject rdf:resource="rdf:#263865bc-7ee9-4520-83cd-77dd03b67ada"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5224a79e-8a89-4f61-8193-78cba8e73246">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4dac4b00-40a6-4882-b3e4-303d5273edbd">
    <dcterms:W3CDTF>2003-09-11</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#73b45ddd-5b85-44c2-9ec6-96df914851b0">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c2df81a-5274-44f2-ac26-ff815319aece">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#faa5ede7-e2bd-48d2-8a49-6c4d5f8220fc">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Rice</vCard:Family>
    <vCard:Other>Jeremy</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#263865bc-7ee9-4520-83cd-77dd03b67ada">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#ec09e0d2-a605-4d35-9afa-d0c5302c2d43"/>
  </rdf:Description>
</rdf:RDF>
</model>