- Author:
- Catherine Lloyd <c.lloyd@auckland.ac.nz>
- Date:
- 2010-04-01 11:18:54+13:00
- Desc:
- Fixed e-notation and updated curation comment.
- Permanent Source URI:
- http://models.cellml.org/workspace/niederer_hunter_smith_2006/rawfile/59a657c9a78c16e4f8f84a2a8e6619fd7c021870/niederer_hunter_smith_2006.cellml
<?xml version='1.0' encoding='utf-8'?>
<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#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="niederer_2006" name="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&itool=EntrezSystem2.PEntrez.Pubmed.Pubmed_ResultsPanel.Pubmed_RVDocSum">PubMed ID: 16339881</ulink>
</para>
</sect1>
</article>
</documentation>
<units name="per_ms">
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<units name="ms">
<unit units="second" prefix="milli"/>
</units>
<units name="mM">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="N_per_mm2">
<unit units="newton"/>
<unit units="metre" prefix="milli" exponent="-2"/>
</units>
<units name="per_mM_per_ms">
<unit units="mM" exponent="-1"/>
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<units name="mM_per_ms">
<unit units="mM"/>
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<component name="environment">
<variable units="ms" public_interface="out" cmeta:id="environment_time" name="time"/>
</component>
<component name="intracellular_ion_concentrations">
<variable units="ms" public_interface="in" name="time"/>
<variable units="mM" public_interface="out" cmeta:id="intracellular_ion_concentrations_Ca_i" name="Ca_i"/>
<variable units="mM" public_interface="out" name="Ca_b"/>
<variable units="mM" public_interface="out" name="TRPN" initial_value="0.067593139865"/>
<variable units="mM" public_interface="in" name="Ca_TRPN_Max"/>
<variable units="mM_per_ms" public_interface="in" name="J_TRPN"/>
<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 cellml:units="dimensionless">1000</cn>
<cn type="e-notation" cellml:units="mM">1.8433<sep/>-7</cn>
</apply>
<apply>
<lt/>
<ci>time</ci>
<cn cellml:units="ms">1</cn>
</apply>
</piece>
<piece>
<apply>
<times/>
<cn cellml:units="dimensionless">1000</cn>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="mM">1.055</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="mM">0.03507</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="mM">0.0003992</cn>
<ci>time</ci>
</apply>
<cn cellml:units="ms">1000</cn>
</apply>
</apply>
<cn type="e-notation" cellml:units="mM">1.356<sep/>-6</cn>
</apply>
</apply>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<cn cellml:units="ms">10</cn>
</apply>
<apply>
<lt/>
<ci>time</ci>
<cn cellml:units="ms">15</cn>
</apply>
</apply>
</piece>
<piece>
<apply>
<times/>
<cn cellml:units="dimensionless">1000</cn>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="mM">0.014</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="mM">0.002555</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="mM">0.0001494</cn>
<ci>time</ci>
</apply>
<cn cellml:units="ms">1000</cn>
</apply>
</apply>
<cn type="e-notation" cellml:units="mM">1.428<sep/>-6</cn>
</apply>
</apply>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<cn cellml:units="ms">15</cn>
</apply>
<apply>
<lt/>
<ci>time</ci>
<cn cellml:units="ms">55</cn>
</apply>
</apply>
</piece>
<piece>
<apply>
<times/>
<cn cellml:units="dimensionless">1000</cn>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<cn type="e-notation" cellml:units="mM">1.739<sep/>-5</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
<apply>
<times/>
<cn type="e-notation" cellml:units="mM">3.209<sep/>-6</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn type="e-notation" cellml:units="mM">5.689<sep/>-6</cn>
<ci>time</ci>
</apply>
<cn cellml:units="ms">1000</cn>
</apply>
</apply>
<cn type="e-notation" cellml:units="mM">1.719<sep/>-6</cn>
</apply>
</apply>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<cn cellml:units="ms">55</cn>
</apply>
<apply>
<lt/>
<ci>time</ci>
<cn cellml:units="ms">250</cn>
</apply>
</apply>
</piece>
<piece>
<apply>
<times/>
<cn cellml:units="dimensionless">1000</cn>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="mM">0.0001321</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="mM">0.0002197</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="mM">0.0001374</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn type="e-notation" cellml:units="mM">3.895<sep/>-5</cn>
<ci>time</ci>
</apply>
<cn cellml:units="ms">1000</cn>
</apply>
</apply>
<cn type="e-notation" cellml:units="mM">4.441<sep/>-6</cn>
</apply>
</apply>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<cn cellml:units="ms">250</cn>
</apply>
<apply>
<lt/>
<ci>time</ci>
<cn cellml:units="ms">490</cn>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<cn cellml:units="dimensionless">1000</cn>
<cn type="e-notation" cellml:units="mM">1.2148<sep/>-7</cn>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="thinfilaments"/>
<component name="tropomyosin">
<variable units="dimensionless" public_interface="out" name="z" initial_value="0.014417937837"/>
<variable units="dimensionless" public_interface="out" name="z_max"/>
<variable units="ms" public_interface="in" name="time"/>
<variable units="per_mM_per_ms" public_interface="in" name="k_on"/>
<variable units="per_ms" public_interface="in" name="k_Ref_off"/>
<variable units="dimensionless" public_interface="in" name="gamma_trpn"/>
<variable units="mM" public_interface="in" name="Ca_TRPN_Max"/>
<variable units="per_ms" name="alpha_0" initial_value="8e-3"/>
<variable units="per_ms" name="alpha_r1" initial_value="2e-3"/>
<variable units="per_ms" name="alpha_r2" initial_value="1.75e-3"/>
<variable units="dimensionless" name="n_Rel" initial_value="3"/>
<variable units="dimensionless" name="K_z" initial_value="0.15"/>
<variable units="dimensionless" name="n_Hill" initial_value="3"/>
<variable units="mM" name="Ca_50ref" initial_value="1.05e-3"/>
<variable units="dimensionless" name="z_p" initial_value="0.85"/>
<variable units="dimensionless" name="beta_1" initial_value="-4"/>
<variable units="mM" name="Ca_50"/>
<variable units="mM" name="Ca_TRPN_50"/>
<variable units="per_ms" name="K_2"/>
<variable units="per_ms" name="K_1"/>
<variable units="per_ms" name="alpha_Tm"/>
<variable units="per_ms" name="beta_Tm"/>
<variable units="dimensionless" public_interface="in" name="beta_0"/>
<variable units="mM" public_interface="in" name="Ca_b"/>
<variable units="dimensionless" public_interface="in" name="lambda"/>
<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 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 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 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 cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>beta_1</ci>
<apply>
<minus/>
<ci>lambda</ci>
<cn 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 cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>beta_0</ci>
<apply>
<minus/>
<ci>lambda</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<cn 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 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 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="troponin">
<variable units="mM_per_ms" public_interface="out" name="J_TRPN"/>
<variable units="mM" public_interface="out" name="Ca_TRPN_Max" initial_value="70e-3"/>
<variable units="mM" public_interface="in" name="TRPN"/>
<variable units="mM" public_interface="in" name="Ca_i"/>
<variable units="per_ms" name="k_off"/>
<variable units="per_mM_per_ms" public_interface="out" name="k_on" initial_value="100"/>
<variable units="per_ms" public_interface="out" name="k_Ref_off" initial_value="0.2"/>
<variable units="dimensionless" public_interface="out" name="gamma_trpn" initial_value="2"/>
<variable units="N_per_mm2" public_interface="in" name="Tension"/>
<variable units="N_per_mm2" public_interface="in" name="T_ref"/>
<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 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 cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Tension</ci>
<apply>
<times/>
<ci>gamma_trpn</ci>
<ci>T_ref</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">0.1</cn>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<ci>k_Ref_off</ci>
<cn 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="Myofilaments">
<variable units="ms" public_interface="in" name="time"/>
<variable units="dimensionless" public_interface="out" name="lambda"/>
<variable units="dimensionless" public_interface="out" name="ExtensionRatio"/>
<variable units="per_ms" public_interface="out" name="dExtensionRatiodt"/>
<variable units="dimensionless" public_interface="out" name="lambda_prev"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>ExtensionRatio</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">1</cn>
<apply>
<gt/>
<ci>time</ci>
<cn cellml:units="ms">1000</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>lambda_prev</ci>
<ci>ExtensionRatio</ci>
</apply>
<apply>
<eq/>
<ci>dExtensionRatiodt</ci>
<cn 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 cellml:units="dimensionless">0.8</cn>
</apply>
<apply>
<leq/>
<ci>ExtensionRatio</ci>
<cn cellml:units="dimensionless">1.15</cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">1.15</cn>
<apply>
<gt/>
<ci>ExtensionRatio</ci>
<cn cellml:units="dimensionless">1.15</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0.8</cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="filament_overlap">
<variable units="dimensionless" public_interface="in" name="lambda"/>
<variable units="dimensionless" public_interface="out" name="overlap"/>
<variable units="dimensionless" public_interface="out" name="beta_0" initial_value="4.9"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>overlap</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>beta_0</ci>
<apply>
<minus/>
<ci>lambda</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="length_independent_tension">
<variable units="N_per_mm2" public_interface="out" name="T_ref" initial_value="56.2"/>
<variable units="N_per_mm2" public_interface="out" name="T_Base"/>
<variable units="dimensionless" public_interface="in" name="z"/>
<variable units="dimensionless" public_interface="in" name="z_max"/>
<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="isometric_tension">
<variable units="N_per_mm2" public_interface="out" cmeta:id="isometric_tension_T0" name="T_0"/>
<variable units="N_per_mm2" public_interface="in" name="T_Base"/>
<variable units="dimensionless" public_interface="in" name="overlap"/>
<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="Cross_Bridges">
<variable units="ms" public_interface="in" name="time"/>
<variable units="per_ms" public_interface="in" name="dExtensionRatiodt"/>
<variable units="dimensionless" name="Q"/>
<variable units="dimensionless" name="a" initial_value="0.35"/>
<variable units="dimensionless" name="Q_1" initial_value="0"/>
<variable units="dimensionless" name="Q_2" initial_value="0"/>
<variable units="dimensionless" name="Q_3" initial_value="0"/>
<variable units="dimensionless" name="A_1" initial_value="-29"/>
<variable units="dimensionless" name="A_2" initial_value="138"/>
<variable units="dimensionless" name="A_3" initial_value="129"/>
<variable units="per_ms" name="alpha_1" initial_value="0.03"/>
<variable units="per_ms" name="alpha_2" initial_value="0.13"/>
<variable units="per_ms" name="alpha_3" initial_value="0.625"/>
<variable units="N_per_mm2" public_interface="in" name="T_0"/>
<variable units="N_per_mm2" public_interface="out" cmeta:id="Cross_Bridges_Tension" name="Tension"/>
<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 cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>Q</ci>
</apply>
</apply>
<apply>
<lt/>
<ci>Q</ci>
<cn cellml:units="dimensionless">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<ci>T_0</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<plus/>
<ci>a</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>Q</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn 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>
<connection>
<map_components component_2="environment" component_1="intracellular_ion_concentrations"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="tropomyosin"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="Myofilaments"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="Cross_Bridges"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="troponin"/>
<map_variables variable_2="Ca_TRPN_Max" variable_1="Ca_TRPN_Max"/>
<map_variables variable_2="TRPN" variable_1="TRPN"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
<map_variables variable_2="J_TRPN" variable_1="J_TRPN"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="tropomyosin"/>
<map_variables variable_2="Ca_b" variable_1="Ca_b"/>
</connection>
<connection>
<map_components component_2="tropomyosin" component_1="length_independent_tension"/>
<map_variables variable_2="z" variable_1="z"/>
<map_variables variable_2="z_max" variable_1="z_max"/>
</connection>
<connection>
<map_components component_2="Cross_Bridges" component_1="troponin"/>
<map_variables variable_2="Tension" variable_1="Tension"/>
</connection>
<connection>
<map_components component_2="length_independent_tension" component_1="troponin"/>
<map_variables variable_2="T_ref" variable_1="T_ref"/>
</connection>
<connection>
<map_components component_2="troponin" component_1="tropomyosin"/>
<map_variables variable_2="k_on" variable_1="k_on"/>
<map_variables variable_2="k_Ref_off" variable_1="k_Ref_off"/>
<map_variables variable_2="gamma_trpn" variable_1="gamma_trpn"/>
<map_variables variable_2="Ca_TRPN_Max" variable_1="Ca_TRPN_Max"/>
</connection>
<connection>
<map_components component_2="filament_overlap" component_1="tropomyosin"/>
<map_variables variable_2="beta_0" variable_1="beta_0"/>
</connection>
<connection>
<map_components component_2="Myofilaments" component_1="tropomyosin"/>
<map_variables variable_2="lambda" variable_1="lambda"/>
</connection>
<connection>
<map_components component_2="Myofilaments" component_1="filament_overlap"/>
<map_variables variable_2="lambda" variable_1="lambda"/>
</connection>
<connection>
<map_components component_2="length_independent_tension" component_1="isometric_tension"/>
<map_variables variable_2="T_Base" variable_1="T_Base"/>
</connection>
<connection>
<map_components component_2="filament_overlap" component_1="isometric_tension"/>
<map_variables variable_2="overlap" variable_1="overlap"/>
</connection>
<connection>
<map_components component_2="Myofilaments" component_1="Cross_Bridges"/>
<map_variables variable_2="dExtensionRatiodt" variable_1="dExtensionRatiodt"/>
</connection>
<connection>
<map_components component_2="isometric_tension" component_1="Cross_Bridges"/>
<map_variables variable_2="T_0" variable_1="T_0"/>
</connection>
<rdf:RDF>
<rdf:Seq rdf:about="rdf:#887c750e-c9a6-4af3-bf5c-e07638d57973">
<rdf:li rdf:resource="rdf:#cd4567e9-796d-415f-8463-bfcc4af7136a"/>
<rdf:li rdf:resource="rdf:#f475ac2a-312b-4d06-9917-68437cf6c8da"/>
<rdf:li rdf:resource="rdf:#8fbfc9a4-eb83-48fe-8882-d7ca7e25046b"/>
</rdf:Seq>
<rdf:Description rdf:about="">
<dc:publisher>Auckland Bioengineering Institute</dc:publisher>
<cmeta:comment rdf:resource="rdf:#22da1e29-04c8-420d-936b-bd8b2e486104"/>
<dcterms:created rdf:resource="rdf:#47b74e5d-97ed-4aaa-a687-07fd8f195b42"/>
<dc:creator rdf:resource="rdf:#e137d61c-82dc-4189-b9f0-29b66171c2e1"/>
<cmeta:modification rdf:resource="rdf:#09f3904b-dd66-4a28-abac-7bb5cd351965"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8fbfc9a4-eb83-48fe-8882-d7ca7e25046b">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#dc21789e-b06b-41ab-af35-2d6feb9c5dbd"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#53e5d308-54da-4612-80cc-b8391fe05593">
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
<rdf:value>j.terkildsen@auckland.ac.nz</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#64bb668b-6f2c-493f-ae16-e2071198e2dd">
<dc:creator rdf:resource="rdf:#982f649b-9600-4031-b4c1-3d8308448733"/>
<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>
<rdf:Description rdf:about="rdf:#dc21789e-b06b-41ab-af35-2d6feb9c5dbd">
<vCard:Given>N</vCard:Given>
<vCard:Family>Smith</vCard:Family>
<vCard:Other>P</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5b2dd97e-1afe-4815-8e78-2542050b7b7d">
<vCard:Given>P</vCard:Given>
<vCard:Family>Hunter</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#09f3904b-dd66-4a28-abac-7bb5cd351965">
<dcterms:modified rdf:resource="rdf:#a4cc9a49-f551-4ae0-8ce6-2c5513162628"/>
<rdf:value>added cmeta id to T_0 variable</rdf:value>
<cmeta:modifier rdf:resource="rdf:#7f8791df-4ef7-4dcc-963d-9a0480c65a43"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#314dbed2-8642-4da9-80fd-945b5d1ec7ea">
<vCard:FN>James Lawson</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#982f649b-9600-4031-b4c1-3d8308448733">
<vCard:FN>James Lawson</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#52c457d5-26a5-4a89-89a2-aec7970789e1">
<dc:creator rdf:resource="rdf:#887c750e-c9a6-4af3-bf5c-e07638d57973"/>
<dc:title>A quantitative analysis of cardiac myocyte relaxation: a simulation study</dc:title>
<bqs:volume>90(5)</bqs:volume>
<bqs:first_page>1697</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#4be66846-659c-4208-a256-233bf4cfad34"/>
<dcterms:issued rdf:resource="rdf:#ad018539-ea8c-4542-954e-02753e2863fc"/>
<bqs:last_page>1722</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$1XTRo2">
<ns7:endingValue>1000</ns7:endingValue>
<ns7:nonstandard-pointDensity>10000</ns7:nonstandard-pointDensity>
<ns7:maximumStepSize>0.1</ns7:maximumStepSize>
</rdf:Description>
<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>
<rdf:Description rdf:about="rdf:#02fbf752-3e20-4809-818c-bd7725477072">
<vCard:Given>Jonna</vCard:Given>
<vCard:Family>Terkildsen</vCard:Family>
<vCard:Other/>
</rdf:Description>
<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"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$+WTRo2">
<ns7:boundIntervals rdf:resource="rdf:#$0XTRo2"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7f8791df-4ef7-4dcc-963d-9a0480c65a43">
<vCard:N rdf:resource="rdf:#5558142b-a43a-435c-aa8d-7da64599a4cc"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4f08decc-3fc3-49f4-9dba-6cfaa15a3c57">
<bqs:Pubmed_id>16339881</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#52c457d5-26a5-4a89-89a2-aec7970789e1"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f475ac2a-312b-4d06-9917-68437cf6c8da">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#5b2dd97e-1afe-4815-8e78-2542050b7b7d"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4749a01a-3930-41a0-b116-99b69afa57f2">
<vCard:Given>S</vCard:Given>
<vCard:Family>Niederer</vCard:Family>
<vCard:Other>A</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5558142b-a43a-435c-aa8d-7da64599a4cc">
<vCard:Given>James</vCard:Given>
<vCard:Family>Lawson</vCard:Family>
<vCard:Other>Richard</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#cd4567e9-796d-415f-8463-bfcc4af7136a">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#4749a01a-3930-41a0-b116-99b69afa57f2"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4be66846-659c-4208-a256-233bf4cfad34">
<dc:title>Biophysical Journal</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a4cc9a49-f551-4ae0-8ce6-2c5513162628">
<dcterms:W3CDTF>2008-08-05T10:32:30+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="#niederer_2006">
<ns7:simulation rdf:resource="rdf:#$+WTRo2"/>
<bqs:reference rdf:resource="rdf:#4f08decc-3fc3-49f4-9dba-6cfaa15a3c57"/>
<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:resource="rdf:#64bb668b-6f2c-493f-ae16-e2071198e2dd"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#22da1e29-04c8-420d-936b-bd8b2e486104">
<dc:creator rdf:resource="rdf:#314dbed2-8642-4da9-80fd-945b5d1ec7ea"/>
<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>
<rdf:Description rdf:about="rdf:#47b74e5d-97ed-4aaa-a687-07fd8f195b42">
<dcterms:W3CDTF>2008-01-07T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ad018539-ea8c-4542-954e-02753e2863fc">
<dcterms:W3CDTF>2006-03-01 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e137d61c-82dc-4189-b9f0-29b66171c2e1">
<vCard:ORG rdf:resource="rdf:#ef3cf40f-3b42-4ed0-993c-730bdf918afc"/>
<vCard:EMAIL rdf:resource="rdf:#53e5d308-54da-4612-80cc-b8391fe05593"/>
<vCard:N rdf:resource="rdf:#02fbf752-3e20-4809-818c-bd7725477072"/>
</rdf:Description>
</rdf:RDF>
</model>
