Location: semantic-annotation @ 6ae217019bf5 / niederer_hunter_smith_2006.cellml

Author:
Dewan Sarwar <sarwarcse@gmail.com>
Date:
2019-02-07 13:16:58+13:00
Desc:
change PR id for warren model
Permanent Source URI:
https://models.physiomeproject.org/workspace/267/rawfile/6ae217019bf59fc1e22e2dc5235751c6fcf6be24/niederer_hunter_smith_2006.cellml

<?xml version="1.0" encoding="UTF-8"?>
<model xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:semsim="http://www.bhi.washington.edu/SemSim#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="niederer_2006" cmeta:id="niederer_2006">
  <documentation xmlns="http://cellml.org/tmp-documentation">
    <article>
      <section id="sec_status">
        <title>Model Status</title>
        <para>This model is known to run in OpenCell and COR to reproduce the output shown in the publication. The units have been checked and they are consistent.</para>
      </section>
      <sect1 id="sec_structure">
        <title>Model Structure</title>
        <para>Abstract: The determinants of relaxation in cardiac muscle are poorly understood, yet compromised relaxation accompanies various pathologies and impaired pump function. In this study, we develop a model of active contraction to elucidate the relative importance of the [Ca2+]i transient magnitude, the unbinding of Ca2+ from troponin C (TnC), and the length-dependence of tension and Ca2+ sensitivity on relaxation. Using the framework proposed by one of our researchers, we extensively reviewed experimental literature, to quantitatively characterize the binding of Ca2+ to TnC, the kinetics of tropomyosin, the availability of binding sites, and the kinetics of crossbridge binding after perturbations in sarcomere length. Model parameters were determined from multiple experimental results and modalities (skinned and intact preparations) and model results were validated against data from length step, caged Ca2+, isometric twitches, and the half-time to relaxation with increasing sarcomere length experiments. A factorial analysis found that the [Ca2+]i transient and the unbinding of Ca2+ from TnC were the primary determinants of relaxation, with a fivefold greater effect than that of length-dependent maximum tension and twice the effect of tension-dependent binding of Ca2+ to TnC and length-dependent Ca2+ sensitivity. The affects of the [Ca2+]i transient and the unbinding rate of Ca2+ from TnC were tightly coupled with the effect of increasing either factor, depending on the reference [Ca2+]i transient and unbinding rate.</para>
        <informalfigure float="0" id="fig_reaction_diagram">
          <mediaobject>
            <imageobject>
              <objectinfo>
                <title>model diagram</title>
              </objectinfo>
              <imagedata fileref="niederer_2006.png" />
            </imageobject>
          </mediaobject>
          <caption>Schematic diagram depicting the relationships of the active contraction framework proposed by Hunter et al. (11). The model is driven by SL and sarcomere velocity, and intracellular [Ca2+]i. Inputs are in bold, algebraic length dependencies are in italics, processes described by differential equations are standard font.</caption>
        </informalfigure>
        <para>The original publication reference is cited below:</para>
        <para>
          A Quantitative Analysis of Cardiac Myocyte Relaxation: A Simulation Study, Steven Niederer, Peter Hunter, Nicholas Smith, 2006
          <emphasis>Biophysical Journal</emphasis>
          , 90 1697-1722
          <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/16339881?ordinalpos=1&amp;itool=EntrezSystem2.PEntrez.Pubmed.Pubmed_ResultsPanel.Pubmed_RVDocSum">PubMed ID: 16339881</ulink>
        </para>
      </sect1>
    </article>
  </documentation>
  <units name="ms">
    <unit prefix="milli" units="second" />
  </units>
  <units name="mM_per_ms">
    <unit exponent="-1.0" prefix="milli" units="second" />
    <unit units="mM" />
  </units>
  <units name="per_mM_per_ms">
    <unit exponent="-1.0" units="mM" />
    <unit exponent="-1.0" prefix="milli" units="second" />
  </units>
  <units name="mM">
    <unit prefix="milli" units="mole" />
    <unit exponent="-1.0" units="litre" />
  </units>
  <units name="per_ms">
    <unit exponent="-1.0" prefix="milli" units="second" />
  </units>
  <units name="N_per_mm2">
    <unit exponent="-2.0" prefix="milli" units="metre" />
    <unit units="newton" />
  </units>
  <component name="troponin">
    <variable cmeta:id="troponin.gamma_trpn" name="gamma_trpn" public_interface="out" units="dimensionless" />
    <variable cmeta:id="troponin.k_on" name="k_on" public_interface="out" units="per_mM_per_ms" />
    <variable cmeta:id="troponin.Ca_TRPN_Max" name="Ca_TRPN_Max" public_interface="out" units="mM" />
    <variable cmeta:id="troponin.TRPN" name="TRPN" public_interface="in" units="mM" />
    <variable cmeta:id="troponin.k_Ref_off" name="k_Ref_off" public_interface="out" units="per_ms" />
    <variable cmeta:id="troponin.k_off" name="k_off" units="per_ms" />
    <variable cmeta:id="troponin.J_TRPN" name="J_TRPN" public_interface="out" units="mM_per_ms" />
    <variable cmeta:id="troponin.Ca_i" name="Ca_i" public_interface="in" units="mM" />
    <variable cmeta:id="troponin.Tension" name="Tension" public_interface="in" units="N_per_mm2" />
    <variable cmeta:id="troponin.T_ref" name="T_ref" public_interface="in" units="N_per_mm2" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>k_off</ci>
        <piecewise>
          <piece>
            <apply>
              <times />
              <ci>k_Ref_off</ci>
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <apply>
                  <divide />
                  <ci>Tension</ci>
                  <apply>
                    <times />
                    <ci>gamma_trpn</ci>
                    <ci>T_ref</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <gt />
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <apply>
                  <divide />
                  <ci>Tension</ci>
                  <apply>
                    <times />
                    <ci>gamma_trpn</ci>
                    <ci>T_ref</ci>
                  </apply>
                </apply>
              </apply>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <times />
              <ci>k_Ref_off</ci>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.1</cn>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <ci>J_TRPN</ci>
        <apply>
          <minus />
          <apply>
            <times />
            <apply>
              <minus />
              <ci>Ca_TRPN_Max</ci>
              <ci>TRPN</ci>
            </apply>
            <ci>k_off</ci>
          </apply>
          <apply>
            <times />
            <ci>Ca_i</ci>
            <ci>TRPN</ci>
            <ci>k_on</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="isometric_tension">
    <variable cmeta:id="isometric_tension.overlap" name="overlap" public_interface="in" units="dimensionless" />
    <variable cmeta:id="isometric_tension.T_Base" name="T_Base" public_interface="in" units="N_per_mm2" />
    <variable cmeta:id="isometric_tension_T0" name="T_0" public_interface="out" units="N_per_mm2" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>T_0</ci>
        <apply>
          <times />
          <ci>T_Base</ci>
          <ci>overlap</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="environment">
    <variable cmeta:id="environment_time" name="time" public_interface="out" units="ms" />
  </component>
  <component name="Cross_Bridges">
    <variable cmeta:id="Cross_Bridges.Q_1" name="Q_1" units="dimensionless" />
    <variable cmeta:id="Cross_Bridges.a" name="a" units="dimensionless" />
    <variable cmeta:id="Cross_Bridges.Q_3" name="Q_3" units="dimensionless" />
    <variable cmeta:id="Cross_Bridges.T_0" name="T_0" public_interface="in" units="N_per_mm2" />
    <variable cmeta:id="Cross_Bridges.alpha_2" name="alpha_2" units="per_ms" />
    <variable cmeta:id="Cross_Bridges.Q_2" name="Q_2" units="dimensionless" />
    <variable cmeta:id="Cross_Bridges_Tension" name="Tension" public_interface="out" units="N_per_mm2" />
    <variable cmeta:id="Cross_Bridges.alpha_1" name="alpha_1" units="per_ms" />
    <variable cmeta:id="Cross_Bridges.alpha_3" name="alpha_3" units="per_ms" />
    <variable cmeta:id="Cross_Bridges.dExtensionRatiodt" name="dExtensionRatiodt" public_interface="in" units="per_ms" />
    <variable cmeta:id="Cross_Bridges.time" name="time" public_interface="in" units="ms" />
    <variable cmeta:id="Cross_Bridges.A_1" name="A_1" units="dimensionless" />
    <variable cmeta:id="Cross_Bridges.Q" name="Q" units="dimensionless" />
    <variable cmeta:id="Cross_Bridges.A_2" name="A_2" units="dimensionless" />
    <variable cmeta:id="Cross_Bridges.A_3" name="A_3" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>Q</ci>
        <apply>
          <plus />
          <ci>Q_1</ci>
          <ci>Q_2</ci>
          <ci>Q_3</ci>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>Tension</ci>
        <piecewise>
          <piece>
            <apply>
              <divide />
              <apply>
                <times />
                <ci>T_0</ci>
                <apply>
                  <plus />
                  <apply>
                    <times />
                    <ci>a</ci>
                    <ci>Q</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                </apply>
              </apply>
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <ci>Q</ci>
              </apply>
            </apply>
            <apply>
              <lt />
              <ci>Q</ci>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide />
              <apply>
                <times />
                <ci>T_0</ci>
                <apply>
                  <plus />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                    <times />
                    <apply>
                      <plus />
                      <ci>a</ci>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                    </apply>
                    <ci>Q</ci>
                  </apply>
                </apply>
              </apply>
              <apply>
                <plus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <ci>Q</ci>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Q_1</ci>
        </apply>
        <apply>
          <minus />
          <apply>
            <times />
            <ci>A_1</ci>
            <ci>dExtensionRatiodt</ci>
          </apply>
          <apply>
            <times />
            <ci>alpha_1</ci>
            <ci>Q_1</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Q_2</ci>
        </apply>
        <apply>
          <minus />
          <apply>
            <times />
            <ci>A_2</ci>
            <ci>dExtensionRatiodt</ci>
          </apply>
          <apply>
            <times />
            <ci>alpha_2</ci>
            <ci>Q_2</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Q_3</ci>
        </apply>
        <apply>
          <minus />
          <apply>
            <times />
            <ci>A_3</ci>
            <ci>dExtensionRatiodt</ci>
          </apply>
          <apply>
            <times />
            <ci>alpha_3</ci>
            <ci>Q_3</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="thinfilaments" />
  <component name="filament_overlap">
    <variable cmeta:id="filament_overlap.overlap" name="overlap" public_interface="out" units="dimensionless" />
    <variable cmeta:id="filament_overlap.lambda" name="lambda" public_interface="in" units="dimensionless" />
    <variable cmeta:id="filament_overlap.beta_0" name="beta_0" public_interface="out" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>overlap</ci>
        <apply>
          <plus />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
          <apply>
            <times />
            <ci>beta_0</ci>
            <apply>
              <minus />
              <ci>lambda</ci>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="tropomyosin">
    <variable cmeta:id="tropomyosin.n_Hill" name="n_Hill" units="dimensionless" />
    <variable cmeta:id="tropomyosin.Ca_b" name="Ca_b" public_interface="in" units="mM" />
    <variable cmeta:id="tropomyosin.z_p" name="z_p" units="dimensionless" />
    <variable cmeta:id="tropomyosin.alpha_r2" name="alpha_r2" units="per_ms" />
    <variable cmeta:id="tropomyosin.alpha_r1" name="alpha_r1" units="per_ms" />
    <variable cmeta:id="tropomyosin.Ca_TRPN_Max" name="Ca_TRPN_Max" public_interface="in" units="mM" />
    <variable cmeta:id="tropomyosin.lambda" name="lambda" public_interface="in" units="dimensionless" />
    <variable cmeta:id="tropomyosin.alpha_0" name="alpha_0" units="per_ms" />
    <variable cmeta:id="tropomyosin.Ca_50ref" name="Ca_50ref" units="mM" />
    <variable cmeta:id="tropomyosin.k_on" name="k_on" public_interface="in" units="per_mM_per_ms" />
    <variable cmeta:id="tropomyosin.beta_0" name="beta_0" public_interface="in" units="dimensionless" />
    <variable cmeta:id="tropomyosin.beta_1" name="beta_1" units="dimensionless" />
    <variable cmeta:id="tropomyosin.K_2" name="K_2" units="per_ms" />
    <variable cmeta:id="tropomyosin.alpha_Tm" name="alpha_Tm" units="per_ms" />
    <variable cmeta:id="tropomyosin.K_1" name="K_1" units="per_ms" />
    <variable cmeta:id="tropomyosin.z_max" name="z_max" public_interface="out" units="dimensionless" />
    <variable cmeta:id="tropomyosin.K_z" name="K_z" units="dimensionless" />
    <variable cmeta:id="tropomyosin.gamma_trpn" name="gamma_trpn" public_interface="in" units="dimensionless" />
    <variable cmeta:id="tropomyosin.n_Rel" name="n_Rel" units="dimensionless" />
    <variable cmeta:id="tropomyosin.Ca_TRPN_50" name="Ca_TRPN_50" units="mM" />
    <variable cmeta:id="tropomyosin.time" name="time" public_interface="in" units="ms" />
    <variable cmeta:id="tropomyosin.beta_Tm" name="beta_Tm" units="per_ms" />
    <variable cmeta:id="tropomyosin.k_Ref_off" name="k_Ref_off" public_interface="in" units="per_ms" />
    <variable cmeta:id="tropomyosin.Ca_50" name="Ca_50" units="mM" />
    <variable cmeta:id="tropomyosin.z" name="z" public_interface="out" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>K_2</ci>
        <apply>
          <times />
          <apply>
            <divide />
            <apply>
              <times />
              <ci>alpha_r2</ci>
              <apply>
                <power />
                <ci>z_p</ci>
                <ci>n_Rel</ci>
              </apply>
            </apply>
            <apply>
              <plus />
              <apply>
                <power />
                <ci>z_p</ci>
                <ci>n_Rel</ci>
              </apply>
              <apply>
                <power />
                <ci>K_z</ci>
                <ci>n_Rel</ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            <apply>
              <divide />
              <apply>
                <times />
                <ci>n_Rel</ci>
                <apply>
                  <power />
                  <ci>K_z</ci>
                  <ci>n_Rel</ci>
                </apply>
              </apply>
              <apply>
                <plus />
                <apply>
                  <power />
                  <ci>z_p</ci>
                  <ci>n_Rel</ci>
                </apply>
                <apply>
                  <power />
                  <ci>K_z</ci>
                  <ci>n_Rel</ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>K_1</ci>
        <apply>
          <divide />
          <apply>
            <times />
            <ci>alpha_r2</ci>
            <apply>
              <power />
              <ci>z_p</ci>
              <apply>
                <minus />
                <ci>n_Rel</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
              </apply>
            </apply>
            <ci>n_Rel</ci>
            <apply>
              <power />
              <ci>K_z</ci>
              <ci>n_Rel</ci>
            </apply>
          </apply>
          <apply>
            <power />
            <apply>
              <plus />
              <apply>
                <power />
                <ci>z_p</ci>
                <ci>n_Rel</ci>
              </apply>
              <apply>
                <power />
                <ci>K_z</ci>
                <ci>n_Rel</ci>
              </apply>
            </apply>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>z_max</ci>
        <apply>
          <divide />
          <apply>
            <minus />
            <apply>
              <divide />
              <ci>alpha_0</ci>
              <apply>
                <power />
                <apply>
                  <divide />
                  <ci>Ca_TRPN_50</ci>
                  <ci>Ca_TRPN_Max</ci>
                </apply>
                <ci>n_Hill</ci>
              </apply>
            </apply>
            <ci>K_2</ci>
          </apply>
          <apply>
            <plus />
            <ci>alpha_r1</ci>
            <ci>K_1</ci>
            <apply>
              <divide />
              <ci>alpha_0</ci>
              <apply>
                <power />
                <apply>
                  <divide />
                  <ci>Ca_TRPN_50</ci>
                  <ci>Ca_TRPN_Max</ci>
                </apply>
                <ci>n_Hill</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>Ca_50</ci>
        <apply>
          <times />
          <ci>Ca_50ref</ci>
          <apply>
            <plus />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            <apply>
              <times />
              <ci>beta_1</ci>
              <apply>
                <minus />
                <ci>lambda</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>Ca_TRPN_50</ci>
        <apply>
          <divide />
          <apply>
            <times />
            <ci>Ca_50</ci>
            <ci>Ca_TRPN_Max</ci>
          </apply>
          <apply>
            <plus />
            <ci>Ca_50</ci>
            <apply>
              <times />
              <apply>
                <divide />
                <ci>k_Ref_off</ci>
                <ci>k_on</ci>
              </apply>
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <apply>
                  <divide />
                  <apply>
                    <times />
                    <apply>
                      <plus />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                      <apply>
                        <times />
                        <ci>beta_0</ci>
                        <apply>
                          <minus />
                          <ci>lambda</ci>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        </apply>
                      </apply>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.5</cn>
                  </apply>
                  <ci>gamma_trpn</ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>alpha_Tm</ci>
        <apply>
          <times />
          <ci>alpha_0</ci>
          <apply>
            <power />
            <apply>
              <divide />
              <ci>Ca_b</ci>
              <ci>Ca_TRPN_50</ci>
            </apply>
            <ci>n_Hill</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>beta_Tm</ci>
        <apply>
          <plus />
          <ci>alpha_r1</ci>
          <apply>
            <divide />
            <apply>
              <times />
              <ci>alpha_r2</ci>
              <apply>
                <power />
                <ci>z</ci>
                <apply>
                  <minus />
                  <ci>n_Rel</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <plus />
              <apply>
                <power />
                <ci>z</ci>
                <ci>n_Rel</ci>
              </apply>
              <apply>
                <power />
                <ci>K_z</ci>
                <ci>n_Rel</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>z</ci>
        </apply>
        <apply>
          <minus />
          <apply>
            <times />
            <ci>alpha_Tm</ci>
            <apply>
              <minus />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
              <ci>z</ci>
            </apply>
          </apply>
          <apply>
            <times />
            <ci>beta_Tm</ci>
            <ci>z</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="intracellular_ion_concentrations">
    <variable cmeta:id="intracellular_ion_concentrations.Ca_TRPN_Max" name="Ca_TRPN_Max" public_interface="in" units="mM" />
    <variable cmeta:id="intracellular_ion_concentrations.J_TRPN" name="J_TRPN" public_interface="in" units="mM_per_ms" />
    <variable cmeta:id="intracellular_ion_concentrations.Ca_b" name="Ca_b" public_interface="out" units="mM" />
    <variable cmeta:id="intracellular_ion_concentrations.TRPN" name="TRPN" public_interface="out" units="mM" />
    <variable cmeta:id="intracellular_ion_concentrations.time" name="time" public_interface="in" units="ms" />
    <variable cmeta:id="intracellular_ion_concentrations_Ca_i" name="Ca_i" public_interface="out" units="mM" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>Ca_b</ci>
        <apply>
          <minus />
          <ci>Ca_TRPN_Max</ci>
          <ci>TRPN</ci>
        </apply>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>TRPN</ci>
        </apply>
        <ci>J_TRPN</ci>
      </apply>
      <apply>
        <eq />
        <ci>Ca_i</ci>
        <piecewise>
          <piece>
            <apply>
              <times />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1000</cn>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                1.8433
                <sep />
                -7
              </cn>
            </apply>
            <apply>
              <lt />
              <ci>time</ci>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1</cn>
            </apply>
          </piece>
          <piece>
            <apply>
              <times />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1000</cn>
              <apply>
                <minus />
                <apply>
                  <plus />
                  <apply>
                    <minus />
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM">1.055</cn>
                      <apply>
                        <power />
                        <apply>
                          <divide />
                          <ci>time</ci>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                      </apply>
                    </apply>
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM">0.03507</cn>
                      <apply>
                        <power />
                        <apply>
                          <divide />
                          <ci>time</ci>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                      </apply>
                    </apply>
                  </apply>
                  <apply>
                    <divide />
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM">0.0003992</cn>
                      <ci>time</ci>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                  </apply>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                  1.356
                  <sep />
                  -6
                </cn>
              </apply>
            </apply>
            <apply>
              <and />
              <apply>
                <geq />
                <ci>time</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">10</cn>
              </apply>
              <apply>
                <lt />
                <ci>time</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">15</cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <times />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1000</cn>
              <apply>
                <minus />
                <apply>
                  <plus />
                  <apply>
                    <minus />
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM">0.014</cn>
                      <apply>
                        <power />
                        <apply>
                          <divide />
                          <ci>time</ci>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                      </apply>
                    </apply>
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM">0.002555</cn>
                      <apply>
                        <power />
                        <apply>
                          <divide />
                          <ci>time</ci>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                      </apply>
                    </apply>
                  </apply>
                  <apply>
                    <divide />
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM">0.0001494</cn>
                      <ci>time</ci>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                  </apply>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                  1.428
                  <sep />
                  -6
                </cn>
              </apply>
            </apply>
            <apply>
              <and />
              <apply>
                <geq />
                <ci>time</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">15</cn>
              </apply>
              <apply>
                <lt />
                <ci>time</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">55</cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <times />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1000</cn>
              <apply>
                <plus />
                <apply>
                  <minus />
                  <apply>
                    <minus />
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                        1.739
                        <sep />
                        -5
                      </cn>
                      <apply>
                        <power />
                        <apply>
                          <divide />
                          <ci>time</ci>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                      </apply>
                    </apply>
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                        3.209
                        <sep />
                        -6
                      </cn>
                      <apply>
                        <power />
                        <apply>
                          <divide />
                          <ci>time</ci>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                      </apply>
                    </apply>
                  </apply>
                  <apply>
                    <divide />
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                        5.689
                        <sep />
                        -6
                      </cn>
                      <ci>time</ci>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                  </apply>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                  1.719
                  <sep />
                  -6
                </cn>
              </apply>
            </apply>
            <apply>
              <and />
              <apply>
                <geq />
                <ci>time</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">55</cn>
              </apply>
              <apply>
                <lt />
                <ci>time</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">250</cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <times />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1000</cn>
              <apply>
                <plus />
                <apply>
                  <minus />
                  <apply>
                    <plus />
                    <apply>
                      <minus />
                      <apply>
                        <times />
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM">0.0001321</cn>
                        <apply>
                          <power />
                          <apply>
                            <divide />
                            <ci>time</ci>
                            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                          </apply>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">4</cn>
                        </apply>
                      </apply>
                      <apply>
                        <times />
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM">0.0002197</cn>
                        <apply>
                          <power />
                          <apply>
                            <divide />
                            <ci>time</ci>
                            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                          </apply>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
                        </apply>
                      </apply>
                    </apply>
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM">0.0001374</cn>
                      <apply>
                        <power />
                        <apply>
                          <divide />
                          <ci>time</ci>
                          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                      </apply>
                    </apply>
                  </apply>
                  <apply>
                    <divide />
                    <apply>
                      <times />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                        3.895
                        <sep />
                        -5
                      </cn>
                      <ci>time</ci>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
                  </apply>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                  4.441
                  <sep />
                  -6
                </cn>
              </apply>
            </apply>
            <apply>
              <and />
              <apply>
                <geq />
                <ci>time</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">250</cn>
              </apply>
              <apply>
                <lt />
                <ci>time</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">490</cn>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <times />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1000</cn>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" type="e-notation" cellml:units="mM">
                1.2148
                <sep />
                -7
              </cn>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
    </math>
  </component>
  <component name="length_independent_tension">
    <variable cmeta:id="length_independent_tension.T_Base" name="T_Base" public_interface="out" units="N_per_mm2" />
    <variable cmeta:id="length_independent_tension.z_max" name="z_max" public_interface="in" units="dimensionless" />
    <variable cmeta:id="length_independent_tension.T_ref" name="T_ref" public_interface="out" units="N_per_mm2" />
    <variable cmeta:id="length_independent_tension.z" name="z" public_interface="in" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>T_Base</ci>
        <apply>
          <divide />
          <apply>
            <times />
            <ci>T_ref</ci>
            <ci>z</ci>
          </apply>
          <ci>z_max</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="Myofilaments">
    <variable cmeta:id="Myofilaments.lambda" name="lambda" public_interface="out" units="dimensionless" />
    <variable cmeta:id="Myofilaments.dExtensionRatiodt" name="dExtensionRatiodt" public_interface="out" units="per_ms" />
    <variable cmeta:id="Myofilaments.lambda_prev" name="lambda_prev" public_interface="out" units="dimensionless" />
    <variable cmeta:id="Myofilaments.ExtensionRatio" name="ExtensionRatio" public_interface="out" units="dimensionless" />
    <variable cmeta:id="Myofilaments.time" name="time" public_interface="in" units="ms" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>ExtensionRatio</ci>
        <piecewise>
          <piece>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            <apply>
              <gt />
              <ci>time</ci>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">1000</cn>
            </apply>
          </piece>
          <otherwise>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <ci>lambda_prev</ci>
        <ci>ExtensionRatio</ci>
      </apply>
      <apply>
        <eq />
        <ci>dExtensionRatiodt</ci>
        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_ms">0</cn>
      </apply>
      <apply>
        <eq />
        <ci>lambda</ci>
        <piecewise>
          <piece>
            <ci>ExtensionRatio</ci>
            <apply>
              <and />
              <apply>
                <gt />
                <ci>ExtensionRatio</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.8</cn>
              </apply>
              <apply>
                <leq />
                <ci>ExtensionRatio</ci>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.15</cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.15</cn>
            <apply>
              <gt />
              <ci>ExtensionRatio</ci>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.15</cn>
            </apply>
          </piece>
          <otherwise>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.8</cn>
          </otherwise>
        </piecewise>
      </apply>
    </math>
  </component>
  <connection>
    <map_components component_1="length_independent_tension" component_2="troponin" />
    <map_variables variable_1="T_ref" variable_2="T_ref" />
  </connection>
  <connection>
    <map_components component_1="intracellular_ion_concentrations" component_2="troponin" />
    <map_variables variable_1="TRPN" variable_2="TRPN" />
    <map_variables variable_1="Ca_i" variable_2="Ca_i" />
  </connection>
  <connection>
    <map_components component_1="intracellular_ion_concentrations" component_2="tropomyosin" />
    <map_variables variable_1="Ca_b" variable_2="Ca_b" />
  </connection>
  <connection>
    <map_components component_1="troponin" component_2="tropomyosin" />
    <map_variables variable_1="k_Ref_off" variable_2="k_Ref_off" />
    <map_variables variable_1="gamma_trpn" variable_2="gamma_trpn" />
    <map_variables variable_1="Ca_TRPN_Max" variable_2="Ca_TRPN_Max" />
    <map_variables variable_1="k_on" variable_2="k_on" />
  </connection>
  <connection>
    <map_components component_1="Myofilaments" component_2="Cross_Bridges" />
    <map_variables variable_1="dExtensionRatiodt" variable_2="dExtensionRatiodt" />
  </connection>
  <connection>
    <map_components component_1="Myofilaments" component_2="filament_overlap" />
    <map_variables variable_1="lambda" variable_2="lambda" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="intracellular_ion_concentrations" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="tropomyosin" component_2="length_independent_tension" />
    <map_variables variable_1="z_max" variable_2="z_max" />
    <map_variables variable_1="z" variable_2="z" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="Cross_Bridges" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="filament_overlap" component_2="isometric_tension" />
    <map_variables variable_1="overlap" variable_2="overlap" />
  </connection>
  <connection>
    <map_components component_1="Myofilaments" component_2="tropomyosin" />
    <map_variables variable_1="lambda" variable_2="lambda" />
  </connection>
  <connection>
    <map_components component_1="isometric_tension" component_2="Cross_Bridges" />
    <map_variables variable_1="T_0" variable_2="T_0" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="tropomyosin" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="Cross_Bridges" component_2="troponin" />
    <map_variables variable_1="Tension" variable_2="Tension" />
  </connection>
  <connection>
    <map_components component_1="filament_overlap" component_2="tropomyosin" />
    <map_variables variable_1="beta_0" variable_2="beta_0" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="Myofilaments" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="troponin" component_2="intracellular_ion_concentrations" />
    <map_variables variable_1="Ca_TRPN_Max" variable_2="Ca_TRPN_Max" />
    <map_variables variable_1="J_TRPN" variable_2="J_TRPN" />
  </connection>
  <connection>
    <map_components component_1="length_independent_tension" component_2="isometric_tension" />
    <map_variables variable_1="T_Base" variable_2="T_Base" />
  </connection>
  <rdf:RDF xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:bqbiol="http://biomodels.net/biology-qualifiers/" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" xmlns:ro="http://www.obofoundry.org/ro/ro.owl#" xmlns:opb="http://bhi.washington.edu/OPB#" xmlns:j.0="http://imc.org/vCard/3.0#">
    <rdf:Description rdf:about="#troponin.k_on">
      <dcterms:description>The binding rate of calcium (2+) to troponin C site II</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.alpha_r2">
      <dcterms:description>Fast relaxation rate parameter for computing z</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.dExtensionRatiodt">
      <dcterms:description>Change in ExtensionRatio with respect to time</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#intracellular_ion_concentrations.time">
      <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01023" />
      <dcterms:description>Temporal solution domain</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Myofilaments.dExtensionRatiodt">
      <dcterms:description>Change in ExtensionRatio with respect to time</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="http://www.cellml.org/cellml/1.0#niederer_2006">
      <ns7:simulation>
        <rdf:Description rdf:about="rdf:#$+WTRo2">
          <ns7:boundIntervals>
            <rdf:Description rdf:about="rdf:#$0XTRo2">
              <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil" />
              <rdf:first rdf:resource="rdf:#$1XTRo2" ns7:endingValue="1000" ns7:nonstandard-pointDensity="10000" ns7:maximumStepSize="0.1" />
            </rdf:Description>
          </ns7:boundIntervals>
        </rdf:Description>
      </ns7:simulation>
      <bqs:reference>
        <rdf:Description rdf:about="rdf:#4f08decc-3fc3-49f4-9dba-6cfaa15a3c57" bqs:Pubmed_id="16339881">
          <bqs:JournalArticle>
            <rdf:Description rdf:about="rdf:#52c457d5-26a5-4a89-89a2-aec7970789e1" bqs:volume="90(5)" bqs:first_page="1697" bqs:last_page="1722">
              <dc:creator>
                <rdf:Seq rdf:about="rdf:#887c750e-c9a6-4af3-bf5c-e07638d57973">
                  <rdf:li>
                    <bqs:Person rdf:about="rdf:#cd4567e9-796d-415f-8463-bfcc4af7136a">
                      <vCard:N rdf:resource="rdf:#4749a01a-3930-41a0-b116-99b69afa57f2" vCard:Given="S" vCard:Family="Niederer" vCard:Other="A" />
                    </bqs:Person>
                  </rdf:li>
                  <rdf:li>
                    <bqs:Person rdf:about="rdf:#f475ac2a-312b-4d06-9917-68437cf6c8da">
                      <vCard:N rdf:resource="rdf:#5b2dd97e-1afe-4815-8e78-2542050b7b7d" vCard:Given="P" vCard:Family="Hunter" vCard:Other="J" />
                    </bqs:Person>
                  </rdf:li>
                  <rdf:li>
                    <bqs:Person rdf:about="rdf:#8fbfc9a4-eb83-48fe-8882-d7ca7e25046b">
                      <vCard:N rdf:resource="rdf:#dc21789e-b06b-41ab-af35-2d6feb9c5dbd" vCard:Given="N" vCard:Family="Smith" vCard:Other="P" />
                    </bqs:Person>
                  </rdf:li>
                </rdf:Seq>
              </dc:creator>
              <dc:title>A quantitative analysis of cardiac myocyte relaxation: a simulation study</dc:title>
              <bqs:Journal>
                <rdf:Description rdf:about="rdf:#4be66846-659c-4208-a256-233bf4cfad34">
                  <dc:title>Biophysical Journal</dc:title>
                </rdf:Description>
              </bqs:Journal>
              <dcterms:issued>
                <rdf:Description rdf:about="rdf:#ad018539-ea8c-4542-954e-02753e2863fc">
                  <dcterms:W3CDTF>2006-03-01 00:00</dcterms:W3CDTF>
                </rdf:Description>
              </dcterms:issued>
            </rdf:Description>
          </bqs:JournalArticle>
        </rdf:Description>
      </bqs:reference>
      <bqs:reference rdf:parseType="Resource">
        <dc:subject rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>cardiac myocyte</rdf:li>
              <rdf:li>endocrine</rdf:li>
              <rdf:li>electrophysiology</rdf:li>
              <rdf:li>relaxation</rdf:li>
              <rdf:li>excitation-contraction coupling</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>
      <cmeta:comment>
        <rdf:Description rdf:about="rdf:#64bb668b-6f2c-493f-ae16-e2071198e2dd">
          <dc:creator>
            <rdf:Description rdf:about="rdf:#982f649b-9600-4031-b4c1-3d8308448733">
              <vCard:FN>James Lawson</vCard:FN>
            </rdf:Description>
          </dc:creator>
          <rdf:value>This is a CellML description of Niederer et al.'s computational model of cardiac myocyte relaxation as described in the 2006 paper "A quantitative analysis of cardiac myocyte relaxation: a simulation study"</rdf:value>
        </rdf:Description>
      </cmeta:comment>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.alpha_1">
      <dcterms:description>Exponential rate constant for process 1 of the crossbridge kinetics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.Ca_TRPN_Max">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_7">
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>The maximum concentration of calcium ions that can bind to site II on troponin C</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.time">
      <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01023" />
      <dcterms:description>Temporal solution domain</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Myofilaments.time">
      <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01023" />
      <dcterms:description>Temporal solution domain</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#troponin.J_TRPN">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_2">
          <semsim:physicalPropertyOf>
            <rdf:Description rdf:about="#process_0">
              <semsim:hasSinkParticipant>
                <rdf:Description rdf:about="#sink_1" semsim:hasMultiplier="1.0">
                  <semsim:hasPhysicalEntityReference>
                    <rdf:Description rdf:about="#entity_4">
                      <ro:part_of>
                        <rdf:Description rdf:about="#entity_1">
                          <ro:part_of>
                            <rdf:Description rdf:about="#entity_2">
                              <semsim:hasPhysicalDefinition rdf:resource="http://purl.obolibrary.org/obo/FMA_14067" />
                            </rdf:Description>
                          </ro:part_of>
                          <semsim:hasPhysicalDefinition rdf:resource="http://purl.obolibrary.org/obo/FMA_66836" />
                        </rdf:Description>
                      </ro:part_of>
                      <semsim:hasPhysicalDefinition rdf:resource="http://purl.obolibrary.org/obo/CHEBI_29108" />
                    </rdf:Description>
                  </semsim:hasPhysicalEntityReference>
                </rdf:Description>
              </semsim:hasSinkParticipant>
              <semsim:hasSinkParticipant>
                <rdf:Description rdf:about="#sink_0" semsim:hasMultiplier="1.0">
                  <semsim:hasPhysicalEntityReference>
                    <rdf:Description rdf:about="#entity_0">
                      <ro:part_of rdf:resource="#entity_1" />
                      <semsim:hasPhysicalDefinition rdf:resource="http://purl.obolibrary.org/obo/PR_P63316" />
                    </rdf:Description>
                  </semsim:hasPhysicalEntityReference>
                </rdf:Description>
              </semsim:hasSinkParticipant>
              <semsim:hasSourceParticipant>
                <rdf:Description rdf:about="#source_0" semsim:hasMultiplier="1.0">
                  <semsim:hasPhysicalEntityReference>
                    <rdf:Description rdf:about="#entity_3">
                      <ro:part_of rdf:resource="#entity_1" />
                      <dcterms:description>Calcium (2+) ion bound to site II of troponin C</dcterms:description>
                      <semsim:name>Calcium (2+) bound to troponin C</semsim:name>
                    </rdf:Description>
                  </semsim:hasPhysicalEntityReference>
                </rdf:Description>
              </semsim:hasSourceParticipant>
              <dcterms:description>In this model, the formation of troponin C comes from the dissociation of the calcium (2+) / troponin C complex</dcterms:description>
              <semsim:name>Dissociation of troponin C bound to calcium</semsim:name>
            </rdf:Description>
          </semsim:physicalPropertyOf>
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00593" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Intracellular troponin flux</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#troponin.Ca_i">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_3">
          <semsim:physicalPropertyOf rdf:resource="#entity_4" />
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Intracellular calcium (2+) concentration</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.z_max">
      <dcterms:description>the maximum fraction of available actin sites at a given sarcomere length</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#intracellular_ion_concentrations_Ca_i">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_15">
          <semsim:physicalPropertyOf rdf:resource="#entity_4" />
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Intracellular calcium (2+) concentration</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Myofilaments.ExtensionRatio">
      <dcterms:description>Sarcomere strain</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.alpha_2">
      <dcterms:description>Exponential rate constant for process 2 of the crossbridge kinetics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.T_0">
      <dcterms:description>The isometric tension at a given intracellular calcium ion concentration and sarcomere length</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges_Tension">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_5">
          <semsim:physicalPropertyOf rdf:resource="#entity_2" />
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01097" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Active tension of cardiac myocyte</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.Q_1">
      <dcterms:description>Exponential process 1 of the crossbridge kinetics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.alpha_3">
      <dcterms:description>Exponential rate constant for process 3 of the crossbridge kinetics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.gamma_trpn">
      <dcterms:description>A measure of the effect of tension on the unbinding rate of calcium (2+) from troponin C</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.A_1">
      <dcterms:description>Weighting coefficient for process 1 of the crossbridge kinetics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#intracellular_ion_concentrations.Ca_b">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_13">
          <semsim:physicalPropertyOf rdf:resource="#entity_3" />
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Concentration of calcium (2+) bound to troponin C</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#troponin.T_ref">
      <dcterms:description>Maximum tension at resting sarcomere length</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#filament_overlap.beta_0">
      <dcterms:description>Scaling factor (slope) for the relationship between lambda and T_0</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.alpha_r1">
      <dcterms:description>Slow relaxation rate parameter for computing z</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.k_on">
      <dcterms:description>The binding rate of calcium (2+) to troponin C site II</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.k_Ref_off">
      <dcterms:description>Rate of calcium (2+) unbinding from Troponin C in the absence of tension</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.Q_2">
      <dcterms:description>Exponential process 2 of the crossbridge kinetics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.alpha_0">
      <dcterms:description>The monoexponential activation rate seen in caged calcium (2+) experiments. Used to compute z.</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.z">
      <dcterms:description>The fraction of actin sites available for crossbridge binding as determined by thin filament kinetics</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.A_2">
      <dcterms:description>Weighting coefficient for process 2 of the crossbridge kinetics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#troponin.k_off">
      <dcterms:description>The unbinding rate of calcium (2+) from troponin C site II</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.Ca_TRPN_50">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_9">
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Calcium (2+) bound to TnC at half-activation of maximum tension</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.n_Rel">
      <dcterms:description>Fast relaxation rate parameter for computing z</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#environment_time">
      <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01023" />
      <dcterms:description>Temporal solution domain</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.Q_3">
      <dcterms:description>Exponential process 3 of the crossbridge kinetics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.K_2">
      <dcterms:description>Linearization term for computing z_max</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.alpha_Tm">
      <dcterms:description>Intermediate term for first component of varaible z's ODE (equation 9 in reference paper)</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#intracellular_ion_concentrations.TRPN">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_14">
          <semsim:physicalPropertyOf rdf:resource="#entity_0" />
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Intracellular troponin concentration</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#length_independent_tension.z_max">
      <dcterms:description>The fraction of actin sites available at full activation</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.time">
      <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01023" />
      <dcterms:description>Temporal solution domain</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.A_3">
      <dcterms:description>Weighting coefficient for process 3 of the crossbridge kinetics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="http://www.cellml.org/cellml/1.0">
      <dc:publisher>Auckland Bioengineering Institute</dc:publisher>
      <cmeta:comment>
        <rdf:Description rdf:about="rdf:#22da1e29-04c8-420d-936b-bd8b2e486104">
          <dc:creator>
            <rdf:Description rdf:about="rdf:#314dbed2-8642-4da9-80fd-945b5d1ec7ea">
              <vCard:FN>James Lawson</vCard:FN>
            </rdf:Description>
          </dc:creator>
          <rdf:value>This model can be simulated in COR and PCEnv and is able to reproduce the results published in Niederer et al.'s 2006 paper.</rdf:value>
        </rdf:Description>
      </cmeta:comment>
      <dcterms:created rdf:resource="rdf:#47b74e5d-97ed-4aaa-a687-07fd8f195b42" dcterms:W3CDTF="2008-01-07T00:00:00+00:00" />
      <dc:creator>
        <rdf:Description rdf:about="rdf:#e137d61c-82dc-4189-b9f0-29b66171c2e1">
          <vCard:ORG>
            <rdf:Description rdf:about="rdf:#ef3cf40f-3b42-4ed0-993c-730bdf918afc">
              <vCard:Orgname>The University of Auckland</vCard:Orgname>
              <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
            </rdf:Description>
          </vCard:ORG>
          <vCard:EMAIL rdf:resource="rdf:#53e5d308-54da-4612-80cc-b8391fe05593" rdf:value="j.terkildsen@auckland.ac.nz" rdf:type="http://imc.org/vCard/3.0#internet" />
          <vCard:N rdf:resource="rdf:#02fbf752-3e20-4809-818c-bd7725477072" vCard:Given="Jonna" vCard:Family="Terkildsen" vCard:Other="" />
        </rdf:Description>
      </dc:creator>
      <cmeta:modification>
        <rdf:Description rdf:about="rdf:#09f3904b-dd66-4a28-abac-7bb5cd351965">
          <dcterms:modified rdf:resource="rdf:#a4cc9a49-f551-4ae0-8ce6-2c5513162628" dcterms:W3CDTF="2008-08-05T10:32:30+12:00" />
          <rdf:value>added cmeta id to T_0 variable</rdf:value>
          <cmeta:modifier>
            <rdf:Description rdf:about="rdf:#7f8791df-4ef7-4dcc-963d-9a0480c65a43">
              <vCard:N rdf:resource="rdf:#5558142b-a43a-435c-aa8d-7da64599a4cc" vCard:Given="James" vCard:Family="Lawson" vCard:Other="Richard" />
            </rdf:Description>
          </cmeta:modifier>
        </rdf:Description>
      </cmeta:modification>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.lambda">
      <dcterms:description>Sarcomere strain?</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.z_p">
      <dcterms:description>Point around which the non-linear component of equation 9 in the paper (solving for dz/dt) is linearized to maintain an explicit formulation for the model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Myofilaments.lambda">
      <dcterms:description>Sarcomere strain, constrained to values between 0.8 and 1.15</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.Ca_b">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_6">
          <semsim:physicalPropertyOf rdf:resource="#entity_3" />
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Concentration of calcium (2+) bound to troponin C</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.beta_0">
      <dcterms:description>Scaling factor (slope) for the relationship between lambda and T_0</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#troponin.Ca_TRPN_Max">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_0">
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>The maximum concentration of calcium ions that can bind to site II on troponin C</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#intracellular_ion_concentrations.J_TRPN">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_12">
          <semsim:physicalPropertyOf rdf:resource="#process_0" />
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00593" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Intracellular troponin flux</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#length_independent_tension.z">
      <dcterms:description>The fraction of actin sites available for crossbridge binding</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#length_independent_tension.T_ref">
      <dcterms:description>Maximum tension at resting sarcomere length</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.Ca_50">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_10">
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Concentration of free calcium (2+) at half-activation of maximum tension</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#intracellular_ion_concentrations.Ca_TRPN_Max">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_11">
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>The maximum concentration of calcium ions that can bind to site II on troponin C</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#isometric_tension.T_Base">
      <dcterms:description>Isometric tension of cell</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Myofilaments.lambda_prev">
      <dcterms:description>Mystery variable mathematically equal to sarcomere strain</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#filament_overlap.lambda">
      <dcterms:description>Sarcomere strain?</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#filament_overlap.overlap">
      <dcterms:description>Sarcomere overlap?</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.K_z">
      <dcterms:description>Fast relaxation rate parameter for computing z</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.beta_1">
      <dcterms:description>Scaling parameter for Ca_50 vs. lambda dependency</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#troponin.k_Ref_off">
      <dcterms:description>Rate of calcium (2+) unbinding from Troponin C in the absence of tension</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.n_Hill">
      <dcterms:description>Analogous to the Hill coefficient in the steady-state force-calcium (2+) curve (F-pCa) and provides a phenomenological rep- resentation of the high cooperativity, due to nearest-neighbor interactions between tropomyosin and/or crossbridges, seen in the activation of tension in cardiac muscle</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.Ca_50ref">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_8">
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Offset parameter for Ca_50 vs. lambda dependency</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#troponin.TRPN">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_1">
          <semsim:physicalPropertyOf rdf:resource="#entity_0" />
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Intracellular troponin concentration</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.K_1">
      <dcterms:description>Linearization term for computing z_max</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.a">
      <dcterms:description>Curve-shaping constant for relationship between relative sarcomere velocity and relative tension</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#tropomyosin.beta_Tm">
      <dcterms:description>Intermediate term for second component of variable z's ODE (equation 9 in reference paper)</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#troponin.gamma_trpn">
      <dcterms:description>A measure of the effect of tension on the unbinding rate of calcium (2+) from troponin C</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#length_independent_tension.T_Base">
      <dcterms:description>Isometric tension of cell</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#troponin.Tension">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_4">
          <semsim:physicalPropertyOf rdf:resource="#entity_2" />
          <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/opb/OPB_01097" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>Active tension of cardiac myocyte</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Cross_Bridges.Q">
      <dcterms:description>Sum of the three exponential processes in the crossbridge dynamics fading memory model</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#isometric_tension_T0">
      <dcterms:description>The isometric tension at a given intracellular calcium ion concentration and sarcomere length</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#niederer_hunter_smith_2006" semsim:modelId="niederer_2006" semsim:modelName="niederer_2006" />
    <rdf:Description rdf:about="#niederer_hunter_smith_2006" semsim:modelName="niederer_hunter_smith_2006" semsim:unknown="niederer_hunter_smith_2006">
      <ro:modelOf rdf:resource="http://purl.obolibrary.org/obo/PR_P63316" />
    </rdf:Description>
  </rdf:RDF>
</model>