Location: Sobie, Dilly, Dossantoscruz, Lederer, Jafri, 2002 @ 19a4c1b7442f / sobie_dilly_dossantoscruz_lederer_jafri_2002.cellml

Author:
Geoff Nunns <gnunns1@jhu.edu>
Date:
2010-05-12 16:33:35+12:00
Desc:
Edited model status and simplified structure to include only abstract.
Permanent Source URI:
https://models.physiomeproject.org/workspace/sobie_dilly_dossantoscruz_lederer_jafri_2002/rawfile/19a4c1b7442f814e91d31410d752d8b5512d2c03/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 model is a faithful representation of the published equations, but it is not compatible with any CellML tools. It contains the sum MathML element, which is not supported by PCEnv, OpenCell or COR, rendering further curation useless. Inspection of the model suggests that there are significant errors beyond the presence of the sum operator that would render it unsuitable for description in CellML. 
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: A Ca(2+) spark arises when a cluster of sarcoplasmic reticulum (SR) channels (ryanodine receptors or RyRs) opens to release calcium in a locally regenerative manner. Normally triggered by Ca(2+) influx across the sarcolemmal or transverse tubule membrane neighboring the cluster, the Ca(2+) spark has been shown to be the elementary Ca(2+) signaling event of excitation-contraction coupling in heart muscle. However, the question of how the Ca(2+) spark terminates remains a central, unresolved issue. Here we present a new model, "sticky cluster," of SR Ca(2+) release that simulates Ca(2+) spark behavior and enables robust Ca(2+) spark termination. Two newly documented features of RyR behavior have been incorporated in this otherwise simple model: "coupled gating" and an opening rate that depends on SR lumenal [Ca(2+)]. Using a Monte Carlo method, local Ca(2+)-induced Ca(2+) release from clusters containing between 10 and 100 RyRs is modeled. After release is triggered, Ca(2+) flux from RyRs diffuses into the cytosol and binds to intracellular buffers and the fluorescent Ca(2+) indicator fluo-3 to produce the model Ca(2+) spark. Ca(2+) sparks generated by the sticky cluster model resemble those observed experimentally, and Ca(2+) spark duration and amplitude are largely insensitive to the number of RyRs in a cluster. As expected from heart cell investigation, the spontaneous Ca(2+) spark rate in the model increases with elevated cytosolic or SR lumenal [Ca(2+)]. Furthermore, reduction of RyR coupling leads to prolonged model Ca(2+) sparks just as treatment with FK506 lengthens Ca(2+) sparks in heart cells. This new model of Ca(2+) spark behavior provides a "proof of principle" test of a new hypothesis for Ca(2+) spark termination and reproduces critical features of Ca(2+) sparks observed experimentally.
</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>