- Author:
- pmr2.import <nobody@models.cellml.org>
- Date:
- 2008-07-15 03:14:42+12:00
- Desc:
- committing version01 of schneider_shimayoshi_amano_matsuda_2006
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/schneider_shimayoshi_amano_matsuda_2006/rawfile/3273fd2e73032f00f7a81bb38e82127df473b141/schneider_shimayoshi_amano_matsuda_2006.cellml
<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 14/07/2008 at 2:41:35 at p.m. using:
COR (0.9.31.979)
Copyright 2002-2008 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk
CellML 1.0 was used to generate this model
http://www.CellML.org/
--><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#" name="schneider_shimayoshi_amano_matsuda_2006" cmeta:id="schneider_shimayoshi_amano_matsuda_2006">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Mechanism of the Frank-Starling law-- A simulation study with a novel cardiac muscle contraction model that includes titin and tropomyosin</title>
<author>
<firstname>Geoffrey</firstname>
<surname>Nunns</surname>
<affiliation>
<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This model is known to run in both PCEnv and COR. The current model is based on the M1 model from the paper, but other variants are possible via simple modifications. In particular, to recreate the published results found in most of the paper, Ca(2+) concentration was modeled as an explicit function of time, and graphed on a log scale.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
Abstract: A stretch-induced increase of active tension is one of the most important properties of the heart, known as the Frank-Starling law. Although a variation of myofilament Ca(2+) sensitivity with sarcomere length (SL) change was found to be involved, the underlying molecular mechanisms are not fully clarified. Some recent experimental studies indicate that a reduction of the lattice spacing between thin and thick filaments, through the increase of passive tension caused by the sarcomeric protein titin with an increase in SL within the physiological range, promotes formation of force-generating crossbridges (Xbs). However, the mechanism by which the Xb concentration determines the degree of cooperativity for a given SL has so far evaded experimental elucidation. In this simulation study, a novel, rather simple molecular-based cardiac contraction model, appropriate for integration into a ventricular cell model, was designed, being the first model to introduce experimental data on titin-based radial tension to account for the SL-dependent modulation of the interfilament lattice spacing and to include a conformational change of troponin I (TnI). Simulation results for the isometric twitch contraction time course, the length-tension and the force-[Ca(2+)] relationships are comparable to experimental data. A complete potential Frank-Starling mechanism was analyzed by this simulation study. The SL-dependent modulation of the myosin binding rate through titin's passive tension determines the Xb concentration which then alters the degree of positive cooperativity affecting the rate of the TnI conformation change and causing the Hill coefficient to be SL-dependent.</para>
<para>
The complete original paper reference is cited below:
</para>
<para>
<ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6WK6-4KGG4FG-1&_user=140507&_rdoc=1&_fmt=&_orig=search&_sort=d&view=c&_acct=C000011498&_version=1&_urlVersion=0&_userid=140507&md5=5c8dba4535cc3667fcef1abaf8b54a9f">Mechanism of the Frank-Starling law-- A simulation study with a novel cardiac muscle contraction model that includes titin and tropomyosin</ulink>, Natalie S. Schneider, Takao Shimayoshi, Akira Amano, Tetsuya Matsuda, 2006, <ulink url="http://www.usouthal.edu/ishr/journal.htm">
<emphasis>Journal of Molecular and Cellular Cardiology</emphasis>
</ulink>, 41, 522-536. (<ulink url="http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6WK6-4KGG4FG-1-H&_cdi=6898&_user=140507&_orig=search&_coverDate=09%2F30%2F2006&_sk=999589996&view=c&wchp=dGLbVlz-zSkWA&md5=ee86e84e4d20580a0179e4d4b2cae871&ie=/sdarticle.pdf">PDF</ulink> versions of the article are available on the Science Direct website.) <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=16860336&dopt=Abstract">PubMed ID: 16860336 </ulink>
</para>
</sect1>
</article>
</documentation>
<units name="millisecond">
<unit units="second" prefix="milli"/>
</units>
<units name="mM">
<unit units="mole" prefix="milli"/>
<unit units="liter" exponent="-1"/>
</units>
<units name="mM_per_msec">
<unit units="mM"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="per_msec">
<unit units="millisecond" exponent="-1"/>
</units>
<units name="per_mM_per_msec">
<unit units="mM" exponent="-1"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="mN_per_mm2">
<unit units="newton" prefix="milli"/>
<unit units="meter" prefix="milli" exponent="-2"/>
</units>
<units name="mmeter">
<unit units="meter" prefix="milli"/>
</units>
<units name="mN_per_mm2_per_mM">
<unit units="mN_per_mm2"/>
<unit units="mM" exponent="-1"/>
</units>
<component name="environment">
<variable units="millisecond" public_interface="out" name="time"/>
</component>
<component name="parameters">
<variable units="mN_per_mm2_per_mM" public_interface="out" name="A_iso" initial_value="9000"/>
<variable units="per_msec" public_interface="out" name="k_off" initial_value="0.2"/>
<variable units="per_mM_per_msec" public_interface="out" name="k_on" initial_value="17.3"/>
<variable units="per_msec" public_interface="out" name="k_onI" initial_value="0.2"/>
<variable units="per_msec" public_interface="out" name="k_offi" initial_value="0.075"/>
<variable units="dimensionless" public_interface="out" name="k_OFF" initial_value="-0.32"/>
<variable units="dimensionless" public_interface="out" name="k_xboff" initial_value="-1.37"/>
<variable units="per_msec" public_interface="out" name="k_tmoff" initial_value="0.067"/>
<variable units="per_msec" public_interface="out" name="k_tmonc" initial_value="0.014"/>
<variable units="dimensionless" public_interface="out" name="k_tmRU" initial_value="10"/>
<variable units="dimensionless" public_interface="out" name="k_tmxb" initial_value="-1.86"/>
<variable units="per_mM_per_msec" public_interface="out" name="k_12" initial_value="2"/>
<variable units="per_msec" public_interface="out" name="k_minus12" initial_value="0.6"/>
<variable units="dimensionless" public_interface="out" name="k_Titin" initial_value="0.2"/>
<variable units="per_msec" public_interface="out" name="k_3c" initial_value="0.025"/>
<variable units="dimensionless" public_interface="out" name="k_3f" initial_value="50"/>
<variable units="dimensionless" public_interface="out" name="k_3xb" initial_value="-1.23"/>
<variable units="per_msec" public_interface="out" name="k_minus3" initial_value="0.008"/>
<variable units="per_msec" public_interface="out" name="k_4" initial_value="0.077"/>
<variable units="per_msec" public_interface="out" name="k_minus4" initial_value="0.001"/>
<variable units="per_msec" public_interface="out" name="k_5" initial_value="0.03723"/>
</component>
<component name="active_force">
<variable units="mN_per_mm2" name="F_b"/>
<variable units="mN_per_mm2_per_mM" public_interface="in" name="A_iso"/>
<variable units="mM" public_interface="in" name="RUA_MADPPi"/>
<variable units="mM" public_interface="in" name="RUA_MADP"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>F_b</ci>
<apply>
<times/>
<ci>A_iso</ci>
<apply>
<plus/>
<ci>RUA_MADPPi</ci>
<ci>RUA_MADP</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Ca_binding_rate">
<variable units="mM_per_msec" public_interface="out" name="Q_CaB"/>
<variable units="mM" public_interface="in" name="RUNA"/>
<variable units="mM" public_interface="in" name="Ca"/>
<variable units="mM" public_interface="in" name="RUTCa_off"/>
<variable units="per_mM_per_msec" public_interface="in" name="k_on"/>
<variable units="per_msec" public_interface="in" name="k_off"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Q_CaB</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_on</ci>
<ci>RUNA</ci>
<ci>Ca</ci>
</apply>
<apply>
<times/>
<ci>k_off</ci>
<ci>RUTCa_off</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="TnI_conf_change_rate">
<variable units="mM_per_msec" public_interface="out" name="Q_TCaA"/>
<variable units="per_msec" name="k_offI"/>
<variable units="per_msec" public_interface="in" name="k_onI"/>
<variable units="per_msec" public_interface="in" name="k_offi"/>
<variable units="dimensionless" public_interface="in" name="k_OFF"/>
<variable units="dimensionless" public_interface="in" name="k_xboff"/>
<variable units="mM" public_interface="in" name="RUTCa_off"/>
<variable units="mM" public_interface="in" name="RUTCa_on"/>
<variable units="mM" public_interface="in" name="RUA"/>
<variable units="mM" public_interface="in" name="RU_total"/>
<variable units="mM" public_interface="in" name="RUA_MADPPi"/>
<variable units="mM" public_interface="in" name="RUA_MADP"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Q_TCaA</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_onI</ci>
<ci>RUTCa_off</ci>
</apply>
<apply>
<times/>
<ci>k_offI</ci>
<ci>RUTCa_on</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_offI</ci>
<apply>
<times/>
<ci>k_offi</ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>k_OFF</ci>
<apply>
<divide/>
<ci>RUA</ci>
<ci>RU_total</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>k_xboff</ci>
<apply>
<divide/>
<apply>
<plus/>
<ci>RUA_MADPPi</ci>
<ci>RUA_MADP</ci>
</apply>
<ci>RU_total</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">4.4</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Tm_conf_change_rate">
<variable units="mM_per_msec" public_interface="out" name="Q_TMA"/>
<variable units="per_msec" name="k_tmon"/>
<variable units="per_msec" public_interface="in" name="k_tmoff"/>
<variable units="dimensionless" public_interface="in" name="k_tmRU"/>
<variable units="per_msec" public_interface="in" name="k_tmonc"/>
<variable units="dimensionless" public_interface="in" name="k_tmxb"/>
<variable units="mM" public_interface="in" name="RUTM_on"/>
<variable units="mM" public_interface="in" name="RUTCa_on"/>
<variable units="mM" public_interface="in" name="RU_total"/>
<variable units="mM" public_interface="in" name="RUA_MADPPi"/>
<variable units="mM" public_interface="in" name="RUA_MADP"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Q_TMA</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_tmon</ci>
<ci>RUTCa_on</ci>
</apply>
<apply>
<times/>
<ci>k_tmoff</ci>
<ci>RUTM_on</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_tmon</ci>
<apply>
<times/>
<ci>k_tmonc</ci>
<apply>
<divide/>
<ci>RUTCa_on</ci>
<ci>RU_total</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>k_tmRU</ci>
<apply>
<divide/>
<ci>RUTM_on</ci>
<ci>RU_total</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>k_tmxb</ci>
<apply>
<divide/>
<apply>
<plus/>
<ci>RUA_MADPPi</ci>
<ci>RUA_MADP</ci>
</apply>
<ci>RU_total</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="myosin_binding_rate">
<variable units="mM_per_msec" public_interface="out" name="Q_MB"/>
<variable units="dimensionless" private_interface="in" name="KTitin"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="k_Titin"/>
<variable units="mmeter" public_interface="in" private_interface="out" name="L"/>
<variable units="mmeter" public_interface="in" private_interface="out" name="L_0"/>
<variable units="per_mM_per_msec" public_interface="in" name="k_12"/>
<variable units="per_msec" public_interface="in" name="k_minus12"/>
<variable units="mM" public_interface="in" name="MADPPi"/>
<variable units="mM" public_interface="in" name="RUTM_on_eff"/>
<variable units="mM" public_interface="in" name="RUAMADPPi"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Q_MB</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_12</ci>
<ci>KTitin</ci>
<ci>MADPPi</ci>
<ci>RUTM_on_eff</ci>
</apply>
<apply>
<times/>
<ci>k_minus12</ci>
<ci>RUAMADPPi</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="KTitin">
<variable units="dimensionless" public_interface="out" name="KTitin"/>
<variable units="dimensionless" name="f_titinNormal"/>
<variable units="dimensionless" name="f_titinDamaged"/>
<variable units="dimensionless" public_interface="in" name="k_Titin"/>
<variable units="mmeter" public_interface="in" name="L"/>
<variable units="mmeter" public_interface="in" name="L_0"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>f_titinNormal</ci>
<apply>
<times/>
<ci>k_Titin</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">9.9663</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<divide/>
<ci>L</ci>
<ci>L_0</ci>
</apply>
<cn cellml:units="dimensionless">1.0639</cn>
</apply>
<cn cellml:units="dimensionless">0.0696</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>f_titinDamaged</ci>
<apply>
<times/>
<ci>k_Titin</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">31</cn>
</apply>
<apply>
<divide/>
<ci>L</ci>
<ci>L_0</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">40</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>KTitin</ci>
<piecewise>
<piece>
<ci>f_titinNormal</ci>
<apply>
<leq/>
<ci>L</ci>
<cn cellml:units="mmeter">1.1</cn>
</apply>
</piece>
<otherwise>
<ci>f_titinDamaged</ci>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="Force_Generating_rate">
<variable units="mM_per_msec" public_interface="out" name="Q_Fgen"/>
<variable units="per_msec" name="k_3"/>
<variable units="per_msec" public_interface="in" name="k_3c"/>
<variable units="dimensionless" public_interface="in" name="k_3f"/>
<variable units="dimensionless" public_interface="in" name="k_3xb"/>
<variable units="per_msec" public_interface="in" name="k_minus3"/>
<variable units="mM" public_interface="in" name="RUAMADPPi"/>
<variable units="mM" public_interface="in" name="RUA_MADPPi"/>
<variable units="mM" public_interface="in" name="RUTCa_on_eff"/>
<variable units="mM" public_interface="in" name="RUTM_on_eff"/>
<variable units="mM" public_interface="in" name="RUA_MADP"/>
<variable units="mM" public_interface="in" name="RU_total"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Q_Fgen</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_3</ci>
<ci>RUAMADPPi</ci>
</apply>
<apply>
<times/>
<ci>k_minus3</ci>
<ci>RUA_MADPPi</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_3</ci>
<apply>
<times/>
<ci>k_3c</ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>k_3f</ci>
<apply>
<plus/>
<ci>RUTCa_on_eff</ci>
<ci>RUTM_on_eff</ci>
<ci>RUAMADPPi</ci>
</apply>
</apply>
<ci>RU_total</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>k_3xb</ci>
<apply>
<plus/>
<ci>RUA_MADPPi</ci>
<ci>RUA_MADP</ci>
</apply>
</apply>
<ci>RU_total</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Pi_release_rate">
<variable units="mM_per_msec" public_interface="out" name="Q_PiR"/>
<variable units="per_msec" public_interface="in" name="k_4"/>
<variable units="per_msec" public_interface="in" name="k_minus4"/>
<variable units="mM" public_interface="in" name="RUA_MADPPi"/>
<variable units="mM" public_interface="in" name="RUA_MADP"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Q_PiR</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_4</ci>
<ci>RUA_MADPPi</ci>
</apply>
<apply>
<times/>
<ci>k_minus4</ci>
<ci>RUA_MADP</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="ADP_release_rate">
<variable units="mM_per_msec" public_interface="out" name="Q_ADPR"/>
<variable units="per_msec" public_interface="in" name="k_5"/>
<variable units="mM" public_interface="in" name="RUA_MADP"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Q_ADPR</ci>
<apply>
<times/>
<ci>k_5</ci>
<ci>RUA_MADP</ci>
</apply>
</apply>
</math>
</component>
<component name="RU_species">
<variable units="mM" public_interface="out" name="RUTCa_off" initial_value="0"/>
<variable units="mM" public_interface="out" name="RUTCa_on" initial_value="0"/>
<variable units="mM" public_interface="out" name="RUTM_on" initial_value="0"/>
<variable units="mM" public_interface="out" name="RUAMADPPi" initial_value="0"/>
<variable units="mM" public_interface="out" name="RUA_MADPPi" initial_value="0"/>
<variable units="mM" public_interface="out" name="RUA_MADP" initial_value="0"/>
<variable units="mM" public_interface="out" name="Ca" initial_value="0.01"/>
<variable units="mM" public_interface="out" name="RUNA"/>
<variable units="mM" public_interface="out" name="RUA"/>
<variable units="mM" public_interface="out" name="RUTM_on_eff"/>
<variable units="mM" public_interface="out" name="RUTCa_on_eff"/>
<variable units="mM" public_interface="out" name="RU_total" initial_value="0.0726"/>
<variable units="mM" public_interface="out" name="MADPPi" initial_value="0.1375"/>
<variable units="dimensionless" private_interface="in" name="alpha"/>
<variable units="mmeter" public_interface="in" private_interface="out" name="L"/>
<variable units="mmeter" public_interface="in" private_interface="out" name="L_0"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="mM_per_msec" public_interface="in" name="Q_CaB"/>
<variable units="mM_per_msec" public_interface="in" name="Q_TCaA"/>
<variable units="mM_per_msec" public_interface="in" name="Q_TMA"/>
<variable units="mM_per_msec" public_interface="in" name="Q_MB"/>
<variable units="mM_per_msec" public_interface="in" name="Q_ADPR"/>
<variable units="mM_per_msec" public_interface="in" name="Q_Fgen"/>
<variable units="mM_per_msec" public_interface="in" name="Q_PiR"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>RUTCa_off</ci>
</apply>
<apply>
<minus/>
<ci>Q_CaB</ci>
<ci>Q_TCaA</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>RUTCa_on</ci>
</apply>
<apply>
<minus/>
<ci>Q_TCaA</ci>
<ci>Q_TMA</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>RUTM_on</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>Q_TMA</ci>
<ci>Q_MB</ci>
</apply>
<ci>Q_ADPR</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>RUAMADPPi</ci>
</apply>
<apply>
<minus/>
<ci>Q_MB</ci>
<ci>Q_Fgen</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>RUA_MADPPi</ci>
</apply>
<apply>
<minus/>
<ci>Q_Fgen</ci>
<ci>Q_PiR</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>RUA_MADP</ci>
</apply>
<apply>
<minus/>
<ci>Q_PiR</ci>
<ci>Q_ADPR</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca</ci>
</apply>
<apply>
<minus/>
<ci>Q_CaB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>RUNA</ci>
<apply>
<minus/>
<ci>RU_total</ci>
<apply>
<plus/>
<ci>RUTCa_off</ci>
<ci>RUTCa_on</ci>
<ci>RUTM_on</ci>
<ci>RUAMADPPi</ci>
<ci>RUA_MADPPi</ci>
<ci>RUA_MADP</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>RUA</ci>
<apply>
<plus/>
<ci>RUTCa_on</ci>
<ci>RUTM_on</ci>
<ci>RUAMADPPi</ci>
<ci>RUA_MADPPi</ci>
<ci>RUA_MADP</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>RUTM_on_eff</ci>
<apply>
<times/>
<ci>alpha</ci>
<ci>RUTM_on</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>RUTCa_on_eff</ci>
<apply>
<times/>
<ci>alpha</ci>
<ci>RUTCa_on</ci>
</apply>
</apply>
</math>
</component>
<component name="alpha">
<variable units="dimensionless" public_interface="out" name="alpha"/>
<variable units="mmeter" public_interface="in" name="L"/>
<variable units="mmeter" public_interface="in" name="L_0"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha</ci>
<piecewise>
<piece>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1.5</cn>
<apply>
<divide/>
<ci>L</ci>
<ci>L_0</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">0.5</cn>
</apply>
<apply>
<lt/>
<ci>L</ci>
<cn cellml:units="mmeter">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">1</cn>
<apply>
<and/>
<apply>
<leq/>
<cn cellml:units="mmeter">1</cn>
<ci>L</ci>
</apply>
<apply>
<lt/>
<ci>L</ci>
<cn cellml:units="mmeter">1.1</cn>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1.6</cn>
</apply>
<apply>
<divide/>
<ci>L</ci>
<ci>L_0</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2.76</cn>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="length">
<variable units="mmeter" public_interface="out" name="L"/>
<variable units="mmeter" public_interface="out" name="L_0" initial_value="0.001"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>L</ci>
<ci>L_0</ci>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="myosin_binding_rate">
<component_ref component="KTitin"/>
</component_ref>
<component_ref component="RU_species">
<component_ref component="alpha"/>
</component_ref>
</group>
<connection>
<map_components component_2="RU_species" component_1="environment"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="active_force" component_1="parameters"/>
<map_variables variable_2="A_iso" variable_1="A_iso"/>
</connection>
<connection>
<map_components component_2="Ca_binding_rate" component_1="parameters"/>
<map_variables variable_2="k_on" variable_1="k_on"/>
<map_variables variable_2="k_off" variable_1="k_off"/>
</connection>
<connection>
<map_components component_2="TnI_conf_change_rate" component_1="parameters"/>
<map_variables variable_2="k_onI" variable_1="k_onI"/>
<map_variables variable_2="k_offi" variable_1="k_offi"/>
<map_variables variable_2="k_OFF" variable_1="k_OFF"/>
<map_variables variable_2="k_xboff" variable_1="k_xboff"/>
</connection>
<connection>
<map_components component_2="Tm_conf_change_rate" component_1="parameters"/>
<map_variables variable_2="k_tmoff" variable_1="k_tmoff"/>
<map_variables variable_2="k_tmRU" variable_1="k_tmRU"/>
<map_variables variable_2="k_tmonc" variable_1="k_tmonc"/>
<map_variables variable_2="k_tmxb" variable_1="k_tmxb"/>
</connection>
<connection>
<map_components component_2="myosin_binding_rate" component_1="parameters"/>
<map_variables variable_2="k_12" variable_1="k_12"/>
<map_variables variable_2="k_minus12" variable_1="k_minus12"/>
<map_variables variable_2="k_Titin" variable_1="k_Titin"/>
</connection>
<connection>
<map_components component_2="Force_Generating_rate" component_1="parameters"/>
<map_variables variable_2="k_3c" variable_1="k_3c"/>
<map_variables variable_2="k_3f" variable_1="k_3f"/>
<map_variables variable_2="k_3xb" variable_1="k_3xb"/>
<map_variables variable_2="k_minus3" variable_1="k_minus3"/>
</connection>
<connection>
<map_components component_2="Pi_release_rate" component_1="parameters"/>
<map_variables variable_2="k_4" variable_1="k_4"/>
<map_variables variable_2="k_minus4" variable_1="k_minus4"/>
</connection>
<connection>
<map_components component_2="ADP_release_rate" component_1="parameters"/>
<map_variables variable_2="k_5" variable_1="k_5"/>
</connection>
<connection>
<map_components component_2="RU_species" component_1="active_force"/>
<map_variables variable_2="RUA_MADPPi" variable_1="RUA_MADPPi"/>
<map_variables variable_2="RUA_MADP" variable_1="RUA_MADP"/>
</connection>
<connection>
<map_components component_2="RU_species" component_1="Ca_binding_rate"/>
<map_variables variable_2="Q_CaB" variable_1="Q_CaB"/>
<map_variables variable_2="RUNA" variable_1="RUNA"/>
<map_variables variable_2="Ca" variable_1="Ca"/>
<map_variables variable_2="RUTCa_off" variable_1="RUTCa_off"/>
</connection>
<connection>
<map_components component_2="RU_species" component_1="TnI_conf_change_rate"/>
<map_variables variable_2="Q_TCaA" variable_1="Q_TCaA"/>
<map_variables variable_2="RUTCa_off" variable_1="RUTCa_off"/>
<map_variables variable_2="RUTCa_on" variable_1="RUTCa_on"/>
<map_variables variable_2="RUA" variable_1="RUA"/>
<map_variables variable_2="RU_total" variable_1="RU_total"/>
<map_variables variable_2="RUA_MADPPi" variable_1="RUA_MADPPi"/>
<map_variables variable_2="RUA_MADP" variable_1="RUA_MADP"/>
</connection>
<connection>
<map_components component_2="RU_species" component_1="Tm_conf_change_rate"/>
<map_variables variable_2="Q_TMA" variable_1="Q_TMA"/>
<map_variables variable_2="RUTM_on" variable_1="RUTM_on"/>
<map_variables variable_2="RUTCa_on" variable_1="RUTCa_on"/>
<map_variables variable_2="RU_total" variable_1="RU_total"/>
<map_variables variable_2="RUA_MADPPi" variable_1="RUA_MADPPi"/>
<map_variables variable_2="RUA_MADP" variable_1="RUA_MADP"/>
</connection>
<connection>
<map_components component_2="RU_species" component_1="myosin_binding_rate"/>
<map_variables variable_2="Q_MB" variable_1="Q_MB"/>
<map_variables variable_2="RUTM_on_eff" variable_1="RUTM_on_eff"/>
<map_variables variable_2="RUAMADPPi" variable_1="RUAMADPPi"/>
<map_variables variable_2="MADPPi" variable_1="MADPPi"/>
</connection>
<connection>
<map_components component_2="KTitin" component_1="myosin_binding_rate"/>
<map_variables variable_2="KTitin" variable_1="KTitin"/>
<map_variables variable_2="k_Titin" variable_1="k_Titin"/>
<map_variables variable_2="L" variable_1="L"/>
<map_variables variable_2="L_0" variable_1="L_0"/>
</connection>
<connection>
<map_components component_2="RU_species" component_1="Force_Generating_rate"/>
<map_variables variable_2="Q_Fgen" variable_1="Q_Fgen"/>
<map_variables variable_2="RUTCa_on_eff" variable_1="RUTCa_on_eff"/>
<map_variables variable_2="RUTM_on_eff" variable_1="RUTM_on_eff"/>
<map_variables variable_2="RUAMADPPi" variable_1="RUAMADPPi"/>
<map_variables variable_2="RU_total" variable_1="RU_total"/>
<map_variables variable_2="RUA_MADPPi" variable_1="RUA_MADPPi"/>
<map_variables variable_2="RUA_MADP" variable_1="RUA_MADP"/>
</connection>
<connection>
<map_components component_2="RU_species" component_1="Pi_release_rate"/>
<map_variables variable_2="Q_PiR" variable_1="Q_PiR"/>
<map_variables variable_2="RUA_MADPPi" variable_1="RUA_MADPPi"/>
<map_variables variable_2="RUA_MADP" variable_1="RUA_MADP"/>
</connection>
<connection>
<map_components component_2="RU_species" component_1="ADP_release_rate"/>
<map_variables variable_2="Q_ADPR" variable_1="Q_ADPR"/>
<map_variables variable_2="RUA_MADP" variable_1="RUA_MADP"/>
</connection>
<connection>
<map_components component_2="alpha" component_1="RU_species"/>
<map_variables variable_2="alpha" variable_1="alpha"/>
<map_variables variable_2="L" variable_1="L"/>
<map_variables variable_2="L_0" variable_1="L_0"/>
</connection>
<connection>
<map_components component_2="length" component_1="myosin_binding_rate"/>
<map_variables variable_2="L" variable_1="L"/>
<map_variables variable_2="L_0" variable_1="L_0"/>
</connection>
<connection>
<map_components component_2="length" component_1="RU_species"/>
<map_variables variable_2="L" variable_1="L"/>
<map_variables variable_2="L_0" variable_1="L_0"/>
</connection>
<rdf:RDF>
<rdf:Seq rdf:about="rdf:#55ce0f30-8d03-4881-a10d-a92b4903d87b">
<rdf:li rdf:resource="rdf:#9f62392f-84ea-4eb4-b556-0b5de13767ed"/>
<rdf:li rdf:resource="rdf:#565d3862-1898-4324-9628-08e1c4ec5e01"/>
<rdf:li rdf:resource="rdf:#56571cf7-6691-434b-a6b6-9facc5261831"/>
<rdf:li rdf:resource="rdf:#e8c058c7-bca8-4828-acf7-d8b2ac78c43a"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#f8a02b07-d159-4957-8f17-d54c4d742c99">
<dc:creator rdf:resource="rdf:#f98c9f15-3400-43af-bde6-d184e5a24063"/>
<rdf:value>This model is known to run in both PCEnv and COR. The current model is based on the M1 model from the paper, but other variants are possible via simple modifications. In particular, to recreate the published results found in most of the paper, Ca(2+) concentration was modeled as an explicit function of time, and graphed on a log scale.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#64f78a04-77c1-4300-a207-29b21b1683f1">
<vCard:Given>Natalie</vCard:Given>
<vCard:Family>Schneider</vCard:Family>
<vCard:Other>S</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher/>
<cmeta:comment rdf:resource="rdf:#f8a02b07-d159-4957-8f17-d54c4d742c99"/>
<dcterms:created rdf:resource="rdf:#61f40e62-74c8-434d-b2b1-ad4f37816b9e"/>
<dc:creator rdf:resource="rdf:#0dd5c104-040e-450a-9a58-a8098044f30f"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fabb1d98-babc-469e-8e65-65b0ff514d9b">
<dc:creator rdf:resource="rdf:#55ce0f30-8d03-4881-a10d-a92b4903d87b"/>
<dc:title>Mechanism of the Frank-Starling law-- A simulation study with a novel cardiac muscle contraction model that includes titin and tropomyosin</dc:title>
<bqs:volume>41</bqs:volume>
<bqs:first_page>522</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#0e8c7f06-364c-4397-be49-7bb93c07e435"/>
<dcterms:issued rdf:resource="rdf:#32a5cfd9-72df-482d-9fe4-3e9d7b379052"/>
<bqs:last_page>536</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#56571cf7-6691-434b-a6b6-9facc5261831">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#fa99492e-f80a-41a1-b4d8-7a3b09cf35a0"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#565d3862-1898-4324-9628-08e1c4ec5e01">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#c1850d66-351f-47b1-a6b0-c8696e29b30d"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fd691468-c29d-49b9-a1c3-b02511eb2716">
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
<rdf:value>gnunns1@jhem.jhu.edu</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fa99492e-f80a-41a1-b4d8-7a3b09cf35a0">
<vCard:Given>Akira</vCard:Given>
<vCard:Family>Amano</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#453b3e52-b83b-43b4-ba8f-0d9766af26b0">
<vCard:Orgname>Auckland Bioengineering Institute</vCard:Orgname>
<vCard:Orgunit>CellML Team</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9f62392f-84ea-4eb4-b556-0b5de13767ed">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#64f78a04-77c1-4300-a207-29b21b1683f1"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1a7784ed-f764-4947-aa9c-a811669fc9f1">
<bqs:Pubmed_id>16860336</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#fabb1d98-babc-469e-8e65-65b0ff514d9b"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c38b8a61-a4c4-49c8-b319-af3b16afad61">
<vCard:Given>Tetsuya</vCard:Given>
<vCard:Family>Matsuda</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0dd5c104-040e-450a-9a58-a8098044f30f">
<vCard:ORG rdf:resource="rdf:#453b3e52-b83b-43b4-ba8f-0d9766af26b0"/>
<vCard:EMAIL rdf:resource="rdf:#fd691468-c29d-49b9-a1c3-b02511eb2716"/>
<vCard:N rdf:resource="rdf:#75d46cc5-875b-4176-8e48-678ee3fd513d"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c1850d66-351f-47b1-a6b0-c8696e29b30d">
<vCard:Given>Takao</vCard:Given>
<vCard:Family>Shimayoshi</vCard:Family>
<vCard:Other/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f98c9f15-3400-43af-bde6-d184e5a24063">
<vCard:FN>Geoff Nunns</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e8c058c7-bca8-4828-acf7-d8b2ac78c43a">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#c38b8a61-a4c4-49c8-b319-af3b16afad61"/>
</rdf:Description>
<rdf:Description rdf:about="#schneider_shimayoshi_amano_matsuda_2006">
<bqs:reference rdf:resource="rdf:#1a7784ed-f764-4947-aa9c-a811669fc9f1"/>
<cmeta:comment rdf:resource="rdf:#a4f7fba2-d7ab-42b7-8666-586bb10d48c3"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#61f40e62-74c8-434d-b2b1-ad4f37816b9e">
<dcterms:W3CDTF>2008-07-14T00:00:00+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#954938f8-595b-406d-b80f-6e6769c42615">
<vCard:FN>Geoff Nunns</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a4f7fba2-d7ab-42b7-8666-586bb10d48c3">
<dc:creator rdf:resource="rdf:#954938f8-595b-406d-b80f-6e6769c42615"/>
<rdf:value>This model is known to run in both PCEnv and COR. The current model is based on the M1 model from the paper, but other variants are possible via simple modifications. In particular, to recreate the published results found in most of the paper, Ca(2+) concentration was modeled as an explicit function of time, and graphed on a log scale.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#32a5cfd9-72df-482d-9fe4-3e9d7b379052">
<dcterms:W3CDTF>2006-06-00 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0e8c7f06-364c-4397-be49-7bb93c07e435">
<dc:title>Journal of Molecular and Cellular Cardiology</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#75d46cc5-875b-4176-8e48-678ee3fd513d">
<vCard:Given>Geoffrey</vCard:Given>
<vCard:Family>Nunns</vCard:Family>
<vCard:Other>Rogan</vCard:Other>
</rdf:Description>
</rdf:RDF>
</model>
