Location: Sobie, Dilly, Dossantoscruz, Lederer, Jafri, 2002 @ 5c4a2e4d95b5 / sobie_dilly_dossantoscruz_lederer_jafri_2002.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2009-12-14 16:28:07+13:00
Desc:
Added images in ai and svg format
Permanent Source URI:
https://models.physiomeproject.org/workspace/sobie_dilly_dossantoscruz_lederer_jafri_2002/rawfile/5c4a2e4d95b5b9403a6dcc557f8c6d2e40aaeaf1/sobie_dilly_dossantoscruz_lederer_jafri_2002.cellml

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

CREATED :  26th June 2002

LAST MODIFIED : 9th April 2003

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

DESCRIPTION :  This file contains a CellML description of Sobie et al's
2002 model of calcium-induced calcium release in cardiac muscle.

CHANGES:  
  18/07/2002 - CML - Added more metadata.
  09/04/2003 - AAC - Added publication date information.  
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="sobie_dilly_dossantoscruz_lederer_jafri_2002_version01" name="sobie_dilly_dossantoscruz_lederer_jafri_2002_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Mechanism of Cardiac Ca2+ Spark Termination</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>
The activation of Ca<superscript>2+</superscript> sparks is an essential step in cardiac excitation-contraction coupling.  A Ca<superscript>2+</superscript> spark arises when voltage-gated, L-type Ca<superscript>2+</superscript> channels in the sarcolemmal (SL) membrane open.  Ca<superscript>2+</superscript> diffuses into the cell, causing the concentration of intracellular calcium [Ca<superscript>2+</superscript>]<subscript>i</subscript> to rise, especially in the region close to the channel in the sub-space between the sarcolemmal and the sarcoplasmic reticulum (SR) membranes (see the figure below).  This local increase in [Ca<superscript>2+</superscript>]<subscript>i</subscript> triggers the opening of SR Ca<superscript>2+</superscript>-release channels (ryanodine receptors, or RyRs) located in the sub-space.  Ca<superscript>2+</superscript> from the intracellular SR store is released, causing a large further increase in [Ca<superscript>2+</superscript>]<subscript>i</subscript>, which is recognised as a Ca<superscript>2+</superscript> spark.  Ca<superscript>2+</superscript> spark activation displays both high gain - a small amount of trigger Ca<superscript>2+</superscript> activates a large amount of released Ca<superscript>2+</superscript> - and positive feedback - release of Ca<superscript>2+</superscript> raises [Ca<superscript>2+</superscript>]<subscript>i</subscript> and tends to trigger more Ca<superscript>2+</superscript> release.
</para>

<para>
Despite the extensive knowledge outlined above, there is still uncertainty about the mechanisms underlying the termination of Ca<superscript>2+</superscript> sparks.  In their 2002 paper, Eric A. Sobie, Keith W. Dilly, Jader dos Santos Cruz, W. Jonathan Lederer, and M. Saleet Jafri propose a new hypothesis called <emphasis>sticky cluster</emphasis>, to describe the behaviour of Ca<superscript>2+</superscript> sparks in the heart.  Their mathematical model focuses on how Ca<superscript>2+</superscript> sparks are terminated, and how the variable number of RyRs in a cluster may influence spark behaviour.  They suggest that spark termination is dependent on two factors: the effect of [Ca<superscript>2+</superscript>]<subscript>lumen</subscript> on RyR gating, and and coupled gating between RyRs in a cluster.  
</para>

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

<para>
Termination of Cardiac Ca<superscript>2+</superscript> Sparks: An Investigative Mathematical Model of Calcium-Induced Calcium Release, Eric A. Sobie, Keith W. Dilly, Jader dos Santos Cruz, W. Jonathan Lederer, and M. Saleet Jafri, 2002, <emphasis>Biophysical Journal</emphasis>, 83, 59-78. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=12080100&amp;dopt=Abstract">PubMed ID: 12080100</ulink> 
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell_diagram</title>
    </objectinfo>
    <imagedata fileref="sobie_2002.png"/>
  </imageobject>
</mediaobject>
<caption>Layout of the model elements.</caption>
</informalfigure>

<para>
The model is separated out into two main components: 1) Mathematical description of SR Ca<superscript>2+</superscript> release through an RyR cluster at the extracellular-sarcolemmal junction.  2) Mathematical description of the Ca<superscript>2+</superscript> spark produced by the Ca<superscript>2+</superscript> efflux.  This component of the model is largely based on a previously published model of Ca<superscript>2+</superscript> diffusion, buffering and imaging processes that produce the Ca<superscript>2+</superscript> spark (Smith <emphasis>et al.</emphasis>, 1998).  
</para>


</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="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="nanomolar">
    <unit units="mole" prefix="nano"/>
    <unit units="litre" 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="third_order_rate_constant">
    <unit units="micromolar" exponent="-1"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="picoA">
    <unit units="ampere" prefix="pico"/>
  </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="subspace">
    <variable units="millimolar" public_interface="out" name="Ca_SS"/>
    <variable units="flux" public_interface="out" name="J_RyR"/>
    <variable units="flux" public_interface="out" name="J_buf"/>
    <variable units="micro_litre" public_interface="out" name="V_SS" initial_value="1.0E-13"/>
    
    <variable units="micromolar" name="Ca_myo" initial_value="0.1"/>
    <variable units="coulomb_per_mole" name="F" initial_value="96480.0"/>
    <variable units="flux" name="J_release"/>
    <variable units="flux" name="J_DHPR"/>
    <variable units="flux" name="J_efflux"/>
    <variable units="dimensionless" name="i" initial_value="0.0"/>
    <variable units="dimensionless" name="N" initial_value="50.0"/>
    <variable units="picoA" name="I_DHPR" initial_value="-0.5"/>
    <variable units="first_order_rate_constant" name="D_RyR" initial_value="4000.0"/>
    <variable units="micromolar" name="Ca_B_SR_tot" initial_value="47.0"/>
    <variable units="micromolar" name="Ca_B_CaM_tot" initial_value="24.0"/>
    <variable units="micromolar" name="Ca_B_SL_tot" initial_value="1124.0"/>
    <variable units="micromolar" name="Ca_B_trpnC_tot" initial_value="70.0"/>
    <variable units="millimolar" name="Ca_B_SR"/>
    <variable units="micromolar" name="Ca_B_CaM"/>
    <variable units="micromolar" name="Ca_B_SL"/>
    <variable units="micromolar" name="Ca_B_trpnC"/>
    <variable units="second_order_rate_constant" name="k_on_SR" initial_value="115.0"/>
    <variable units="first_order_rate_constant" name="k_off_SR" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="k_on_CaM" initial_value="100.0"/>
    <variable units="first_order_rate_constant" name="k_off_CaM" initial_value="38.0"/>
    <variable units="second_order_rate_constant" name="k_on_trpnC" initial_value="39.0"/>
    <variable units="first_order_rate_constant" name="k_off_trpnC" initial_value="20.0"/>
    <variable units="second_order_rate_constant" name="k_on_SL" initial_value="115.0"/>
    <variable units="first_order_rate_constant" name="k_off_SL" initial_value="1000.0"/>
    <variable units="second" name="tau_efflux" initial_value="7.0E-7"/>
    <variable units="dimensionless" name="RyR_open" initial_value="0.0"/>
   
    <variable units="millimolar" public_interface="in" name="Ca_lumen"/>
    <variable units="second" public_interface="in" name="time"/>
  
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="calcium_subspace_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_SS </ci>
        </apply>
        <apply>
          <plus/>
          <ci> J_release </ci>
          <ci> J_DHPR </ci>
          <ci> J_efflux </ci>
          <ci> J_buf </ci>
        </apply>
      </apply>
      
      <apply id="J_DHPR_calculation">
        <eq/>
        <ci> J_DHPR </ci>
        <apply>
          <minus/>
          <apply>
            <divide/>
            <ci> I_DHPR </ci>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> F </ci>
              <ci> V_SS </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="J_release_calculation">
        <eq/>
        <ci> J_release </ci>
        <apply>
          <sum/>
          <bvar>
            <ci> i </ci>
          </bvar>
          <lowlimit>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </lowlimit>
          <highlimit>
            <ci> N </ci>
          </highlimit>
          <apply>
            <times/>
            <ci> RyR_open </ci>
            <ci> J_RyR </ci>
          </apply>
        </apply>
      </apply>
    
      <apply id="J_RyR_calculation">
        <eq/>
        <ci> J_RyR </ci>
        <apply>
          <times/>
          <ci> D_RyR </ci>
          <apply>
            <minus/>
            <ci> Ca_lumen </ci>
            <ci> Ca_SS </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="RyR_open_calculation">
        <eq/>
        <ci> RyR_open </ci>
        <piecewise>
          <piece>
            <cn cellml:units="dimensionless"> 0.0 </cn>
            <apply>
              <eq/>
              <ci> i </ci> 
              <cn cellml:units="dimensionless"> 0.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
      
      <apply id="Ca_B_SR_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_B_SR </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_on_SR </ci>
            <ci> Ca_SS </ci>
            <ci> Ca_B_SR </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_off_SR </ci>
            <apply>
              <minus/>
              <ci> Ca_B_SR_tot </ci>
              <ci> Ca_B_SR </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Ca_B_trpnC_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_B_trpnC </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_on_trpnC </ci>
            <ci> Ca_SS </ci>
            <ci> Ca_B_trpnC </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_off_trpnC </ci>
            <apply>
              <minus/>
              <ci> Ca_B_trpnC_tot </ci>
              <ci> Ca_B_trpnC </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Ca_B_CaM_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_B_CaM </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_on_CaM </ci>
            <ci> Ca_SS </ci>
            <ci> Ca_B_CaM </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_off_CaM </ci>
            <apply>
              <minus/>
              <ci> Ca_B_CaM_tot </ci>
              <ci> Ca_B_CaM </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Ca_B_SL_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_B_SL </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k_on_SL </ci>
            <ci> Ca_SS </ci>
            <ci> Ca_B_SL </ci>
          </apply>
          <apply>
            <times/>
            <ci> k_off_SL </ci>
            <apply>
              <minus/>
              <ci> Ca_B_SL_tot </ci>
              <ci> Ca_B_SL </ci>
            </apply>
          </apply>
        </apply>
      </apply>
     
      <apply id="J_buf_calculation">
        <eq/>
        <ci> J_buf </ci>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k_on_CaM </ci>
              <ci> Ca_SS </ci>
              <ci> Ca_B_CaM </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_off_CaM </ci>
              <apply>
                <minus/>
                <ci> Ca_B_CaM_tot </ci>
                <ci> Ca_B_CaM </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k_on_trpnC </ci>
              <ci> Ca_SS </ci>
              <ci> Ca_B_trpnC </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_off_trpnC </ci>
              <apply>
                <minus/>
                <ci> Ca_B_trpnC_tot </ci>
                <ci> Ca_B_trpnC </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k_on_SR </ci>
              <ci> Ca_SS </ci>
              <ci> Ca_B_SR </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_off_SR </ci>
              <apply>
                <minus/>
                <ci> Ca_B_SR_tot </ci>
                <ci> Ca_B_SR </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> k_on_SL </ci>
              <ci> Ca_SS </ci>
              <ci> Ca_B_SL </ci>
            </apply>
            <apply>
              <times/>
              <ci> k_off_SL </ci>
              <apply>
                <minus/>
                <ci> Ca_B_SL_tot </ci>
                <ci> Ca_B_SL </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="J_efflux_calculation">
        <eq/>
        <ci> J_efflux </ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <ci> Ca_myo </ci>
            <ci> Ca_SS </ci>
          </apply>
          <ci> tau_efflux </ci>
        </apply>
      </apply>
    </math>
  </component>
      
  <component name="lumen">
    <variable units="millimolar" public_interface="out" name="Ca_lumen"/>
   
    <variable units="flux" name="J_refill"/>
    <variable units="second" name="tau_refill" initial_value="0.01"/>
    <variable units="millimolar" name="Ca_NSR" initial_value="1.0"/>
    <variable units="micro_litre" name="V_JSR" initial_value="1.0E-11"/>
    <variable units="millimolar" name="CSQ_tot" initial_value="10.0"/>
    <variable units="millimolar" name="CSQ"/>
    <variable units="millimolar" name="K_CSQ" initial_value="0.8"/>
    <variable units="dimensionless" name="beta_JSR"/>
    
    <variable units="flux" public_interface="in" name="J_RyR"/>
    <variable units="micro_litre" public_interface="in" name="V_SS"/>
    <variable units="second" public_interface="in" name="time"/>
  
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="calcium_lumen_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_lumen </ci>
        </apply>
        <apply>
          <times/>
          <ci> beta_JSR </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <minus/>
                <ci> J_RyR </ci>
              </apply>
              <apply>
                <divide/>
                <ci> V_SS </ci>
                <ci> V_JSR </ci>
              </apply>
            </apply>
            <ci> J_refill </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="J_refill_calculation">
        <eq/>
        <ci> J_refill </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> Ca_NSR </ci>
            <ci> Ca_lumen </ci>
          </apply>
          <ci> tau_refill </ci>
        </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> CSQ_tot </ci>
                <ci> K_CSQ </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <plus/>
                  <ci> CSQ </ci>
                  <ci> K_CSQ </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="RyR_receptor_gating">
    <variable units="dimensionless" name="k_coop" initial_value="1.0"/>
    <variable units="first_order_rate_constant" name="k_close"/>
    <variable units="first_order_rate_constant" name="k_open"/>
    <variable units="dimensionless" name="CF_close"/>
    <variable units="dimensionless" name="CF_open"/>
    <variable units="dimensionless" name="N_closed"/>
    <variable units="dimensionless" name="N_open"/>
    <variable units="micromolar" name="K_m"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="millimolar" public_interface="in" name="Ca_lumen"/>
    <variable units="millimolar" public_interface="in" name="Ca_SS"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="k_close_calculation">
        <eq/>
        <ci> k_close </ci>
        <apply>
          <times/>
          <ci> CF_close </ci>
          <cn cellml:units="first_order_rate_constant"> 480.0 </cn>
        </apply>
      </apply>
      
      <apply id="k_open_calculation">
        <eq/>
        <ci> k_open </ci>
        <apply>
          <times/>
          <ci> CF_open </ci>
          <cn cellml:units="first_order_rate_constant"> 30000.0 </cn>
            <apply>
            <divide/>
              <apply>
              <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <apply>
              <plus/>
              <apply>
                <power/>
                <ci> Ca_SS </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
              <apply>
                <power/>
                <ci> K_m </ci>
                <cn cellml:units="dimensionless"> 4.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="K_m_calculation">
        <eq/>
        <ci> K_m </ci>
        <apply>
          <minus/>
          <cn cellml:units="micromolar"> 6.0 </cn>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.0024 </cn>
            <ci> Ca_lumen </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="CF_open_calculation">
        <eq/>
        <ci> CF_open </ci>
        <apply>
          <plus/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <divide/>
            <apply>
              <plus/>
              <ci> N_open </ci>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
            <apply>
              <plus/>
              <ci> N_open </ci>
              <ci> N_closed </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="CF_close_calculation">
        <eq/>
        <ci> CF_close </ci>
        <apply>
          <times/>
          <ci> k_coop </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <apply>
                <plus/>
                <ci> N_closed </ci>
                <cn cellml:units="dimensionless"> 1.0 </cn>
              </apply>
              <apply>
                <plus/>
                <ci> N_open </ci>
                <ci> N_closed </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="spark_model">
    <variable units="millimolar" name="Ca"/>
    <variable units="millimolar" name="Ca_F"/>
    <variable units="millimolar" name="Ca_trpnC"/>
    <variable units="millimolar" name="Ca_CaM"/>
    
    <variable units="third_order_rate_constant" name="DC" initial_value="250.0"/>
    <variable units="third_order_rate_constant" name="DF" initial_value="20.0"/>
    <variable units="micromolar" name="F_T" initial_value="50.0"/>
    <variable units="flux" name="J_pump"/>
    <variable units="flux" name="J_leak"/>
    <variable units="flux" name="J_n"/>
    <variable units="flux" name="J_dye"/>
    <variable units="flux" name="Vmax_pump" initial_value="208.0"/>
    <variable units="micromolar" name="Km_pump" initial_value="0.184"/>
    <variable units="dimensionless" name="m" initial_value="3.9"/>
    <variable units="nanomolar" name="C_infinity" initial_value="100.0"/>
    <variable units="second_order_rate_constant" name="k_F_plus" initial_value="80.0"/>
    <variable units="first_order_rate_constant" name="k_F_minus" initial_value="40.0"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="flux" public_interface="in" name="J_RyR"/>
    <variable units="flux" public_interface="in" name="J_buf"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="cytosolic_calcium_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> DC </ci>
            <ci> Ca </ci>
          </apply>
          <ci> J_dye </ci>
          <ci> J_buf </ci>
          <ci> J_pump </ci>
          <ci> J_leak </ci>
          <ci> J_RyR </ci>
        </apply>
      </apply>
      
      <apply id="calcium_bound_dye_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_F </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> DF </ci>
            <ci> Ca_F </ci>
          </apply>
          <ci> J_dye </ci>
        </apply>
      </apply>
      
      <apply id="J_dye_calculation">
        <eq/>
        <ci> J_dye </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <ci> k_F_plus </ci>
            </apply>
            <ci> Ca </ci>
            <apply>
              <minus/>
              <ci> F_T </ci>
              <ci> Ca_F </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k_F_minus </ci>
            <ci> Ca_F </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="J_pump_calculation">
        <eq/>
        <ci> J_pump </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> Vmax_pump </ci>
            <apply>
              <power/>
              <ci> Ca </ci>
              <ci> m </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <ci> Km_pump </ci>
            <apply>
              <power/>
              <ci> Ca </ci>
              <ci> m </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="J_leak_calculation">
        <eq/>
        <ci> J_leak </ci>
        <apply>
          <minus/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> Vmax_pump </ci>
              <apply>
                <power/>
                <ci> C_infinity </ci>
                <ci> m </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <ci> Km_pump </ci>
              <apply>
                <power/>
                <ci> C_infinity </ci>
                <ci> m </ci>
              </apply>
            </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="lumen"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="RyR_receptor_gating"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="spark_model"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="subspace" component_1="spark_model"/>
        <map_variables variable_2="J_RyR" variable_1="J_RyR"/>
        <map_variables variable_2="J_buf" variable_1="J_buf"/>
  </connection>
  
  <connection>
    <map_components component_2="lumen" 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_lumen" variable_1="Ca_lumen"/>
  </connection>
  
  <connection>
    <map_components component_2="RyR_receptor_gating" component_1="subspace"/>
    <map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
  </connection>
  
  <connection>
    <map_components component_2="RyR_receptor_gating" component_1="lumen"/>
    <map_variables variable_2="Ca_lumen" variable_1="Ca_lumen"/>
  </connection>
  

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#0c06092c-0ce5-42c3-afc7-cf8b7d2e9513">
    <rdf:li>cardiac myocyte</rdf:li>
    <rdf:li>ca sparks</rdf:li>
    <rdf:li>calcium dynamics</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>Cardiac Myocyte</rdf:li>
    <rdf:li>cicr</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#3789e49a-1132-45f9-84d1-65ab35477e59">
    <rdf:li rdf:resource="rdf:#4c1b1a92-02d8-43b1-9112-e468b4d0027f"/>
    <rdf:li rdf:resource="rdf:#44a89961-14b2-43e3-8073-84d60f9cc6e0"/>
    <rdf:li rdf:resource="rdf:#feba6782-78d2-4085-94b9-4b1dee23ec69"/>
    <rdf:li rdf:resource="rdf:#38bbe71e-de75-4c82-967a-0e918405a3cf"/>
    <rdf:li rdf:resource="rdf:#d8577172-8dc5-4280-bb7a-7ab8986c8f60"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#acf72b07-487b-4bff-b4c7-51e13bc7d297">
    <vCard:Given>Jonathan</vCard:Given>
    <vCard:Family>Lederer</vCard:Family>
    <vCard:Other>W</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8d70a8ec-84a8-4d1c-a4a4-4a37d2c0b135">
    <vCard:Given>Eric</vCard:Given>
    <vCard:Family>Sobie</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e2eca531-5ee0-4d84-ac33-fceb0f713606">
    <dc:creator rdf:resource="rdf:#17e0ac8f-5f03-4cca-aef6-d02d905bdbd3"/>
    <rdf:value>
          This is the CellML description of Sobie et al's 2002 model of 
          calcium-induced calcium release in cardiac muscle.
        </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4b9801ab-f585-4a1e-bc89-e99613cfda25">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e346ff5f-6593-4856-a8e5-58966cace86b">
    <dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Bioengineering Institute
      </dc:publisher>
    <cmeta:modification rdf:resource="rdf:#05fbcc77-6976-4ad8-afe3-5867477358a7"/>
    <cmeta:modification rdf:resource="rdf:#6fad6fd2-b9ab-46d8-8f61-8db5652757ae"/>
    <dcterms:created rdf:resource="rdf:#92f6e68b-c4e0-4022-8344-24c197a1866e"/>
    <dc:creator rdf:resource="rdf:#91e40a50-ad90-4db7-82c0-aef4d5cc5af9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#44a89961-14b2-43e3-8073-84d60f9cc6e0">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#da58f571-245e-4c70-8cf2-6a3dbe909a17"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#03dbe7ee-bcb8-4991-910f-26740f0de9c9">
    <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:#32f77fea-fe32-44ec-8907-d2463e603973">
    <dcterms:W3CDTF>2002-07</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a2319129-836e-48be-a6cf-eb5b8ed2cbcd">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#da58f571-245e-4c70-8cf2-6a3dbe909a17">
    <vCard:Given>Keith</vCard:Given>
    <vCard:Family>Dilly</vCard:Family>
    <vCard:Other>W</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d8577172-8dc5-4280-bb7a-7ab8986c8f60">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#78776d53-eb09-40e7-98be-ec22d1f0eda4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e0d5fe24-46ce-4f72-8c27-645e50495865">
    <dc:subject rdf:resource="rdf:#9afb2c30-8444-475b-9518-3f4ab04a6110"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4655d490-7a13-4089-bc08-abb7e873fd90">
    <vCard:N rdf:resource="rdf:#a2319129-836e-48be-a6cf-eb5b8ed2cbcd"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#17e0ac8f-5f03-4cca-aef6-d02d905bdbd3">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6fad6fd2-b9ab-46d8-8f61-8db5652757ae">
    <dcterms:modified rdf:resource="rdf:#e346ff5f-6593-4856-a8e5-58966cace86b"/>
    <rdf:value>
          Added publication date information.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#de151fa2-2641-4500-a591-3bd3252d6e7f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#098b087e-0f5b-43bb-bcfa-f782212e09b2">
    <dc:creator rdf:resource="rdf:#3789e49a-1132-45f9-84d1-65ab35477e59"/>
    <dc:title>
            Termination of Cardiac Ca2+ Sparks: An Investigative Mathematical                Model of Calcium-Induced Calcium Release
          </dc:title>
    <bqs:volume>83</bqs:volume>
    <bqs:first_page>59</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#7669bdbb-711b-4098-90f5-04b2edc3ec17"/>
    <dcterms:issued rdf:resource="rdf:#32f77fea-fe32-44ec-8907-d2463e603973"/>
    <bqs:last_page>78</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2b913c3c-3104-4e98-9f70-f0da21dcb35e">
    <dcterms:W3CDTF>2002-07-18</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#78776d53-eb09-40e7-98be-ec22d1f0eda4">
    <vCard:Given>M</vCard:Given>
    <vCard:Family>Jafri</vCard:Family>
    <vCard:Other>Saleet</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1d71a9df-7961-405e-82a3-e122ff9c80fd">
    <bqs:Pubmed_id>12080100</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#098b087e-0f5b-43bb-bcfa-f782212e09b2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#91e40a50-ad90-4db7-82c0-aef4d5cc5af9">
    <vCard:ORG rdf:resource="rdf:#4b9801ab-f585-4a1e-bc89-e99613cfda25"/>
    <vCard:EMAIL rdf:resource="rdf:#03dbe7ee-bcb8-4991-910f-26740f0de9c9"/>
    <vCard:N rdf:resource="rdf:#c6603a7d-df19-48c1-bc3c-9d6d86c0ee17"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#feba6782-78d2-4085-94b9-4b1dee23ec69">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#285abae7-09ec-4704-ad1d-906af9216f48"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a233a745-7217-412f-b478-fd5e2c5acf91">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#05fbcc77-6976-4ad8-afe3-5867477358a7">
    <dcterms:modified rdf:resource="rdf:#2b913c3c-3104-4e98-9f70-f0da21dcb35e"/>
    <rdf:value>
          Added more metadata.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#4655d490-7a13-4089-bc08-abb7e873fd90"/>
  </rdf:Description>
  <rdf:Description rdf:about="#sobie_dilly_dossantoscruz_lederer_jafri_2002_version01">
    <dc:title>
        Sobie et al's 2002 model of calcium-induced calcium release in cardiac 
        muscle.
      </dc:title>
    <cmeta:bio_entity>Cardiac Myocyte</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#e2eca531-5ee0-4d84-ac33-fceb0f713606"/>
    <bqs:reference rdf:resource="rdf:#1d71a9df-7961-405e-82a3-e122ff9c80fd"/>
    <bqs:reference rdf:resource="rdf:#e0d5fe24-46ce-4f72-8c27-645e50495865"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9afb2c30-8444-475b-9518-3f4ab04a6110">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#0c06092c-0ce5-42c3-afc7-cf8b7d2e9513"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#de151fa2-2641-4500-a591-3bd3252d6e7f">
    <vCard:N rdf:resource="rdf:#a233a745-7217-412f-b478-fd5e2c5acf91"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#92f6e68b-c4e0-4022-8344-24c197a1866e">
    <dcterms:W3CDTF>2002-06-26</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#285abae7-09ec-4704-ad1d-906af9216f48">
    <vCard:Given>Jader</vCard:Given>
    <vCard:Family>dos Santos Cruz</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c6603a7d-df19-48c1-bc3c-9d6d86c0ee17">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7669bdbb-711b-4098-90f5-04b2edc3ec17">
    <dc:title>Biophysical Journal</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4c1b1a92-02d8-43b1-9112-e468b4d0027f">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#8d70a8ec-84a8-4d1c-a4a4-4a37d2c0b135"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#38bbe71e-de75-4c82-967a-0e918405a3cf">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#acf72b07-487b-4bff-b4c7-51e13bc7d297"/>
  </rdf:Description>
</rdf:RDF>
</model>