Location: Rice, Winslow, Hunter, 1999 @ 1403ca54f018 / rice_winslow_hunter_1999_d.cellml

Author:
Dougal Cowan <devnull@localhost>
Date:
2011-04-15 08:54:31+12:00
Desc:
fixed the C model K_Ca initial value (to 0.5) noted the model number for each model
Permanent Source URI:
https://models.physiomeproject.org/workspace/rice_winslow_hunter_1999/rawfile/1403ca54f01826624e0c0deb39344c4ef20f0203/rice_winslow_hunter_1999_d.cellml

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

CREATED :  24th May 2007

LAST MODIFIED : 24th May 2007

AUTHOR :  Catherine Lloyd
          The Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.1 Specification.

DESCRIPTION :  This file contains a CellML description of Rice et al's 1999 1st
               model of isometric force generation in cardiac myofilaments.

CHANGES:  
  
--><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:pathway_editor="http://www.physiome.com/pathway_editor/1.0#" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" pathway_editor:rendering_config_file="rice_model5_1999_CellMLrender.xml" cmeta:id="rice_model4_1999_version02" name="rice_winslow_hunter_1999_version05">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Cooperative Mechanisms in Cardiac Muscle (model 4)</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
In cardiac muscle, steady-state force-Ca<superscript>2+</superscript> (F-Ca) relations exhibit more cooperativity than that predicted by the single Ca<superscript>2+</superscript> binding site on troponin.  The exact mechanisms underlying this high cooperativity are unknown.  In their 1999 paper, J. Jeremy Rice, Raimond L. Winslow and William C. Hunter present five potential models for force generation in cardiac muscle (see the figure below).  These models were constructed by assuming different subsets of three possible cooperative mechanisms: 
<itemizedlist>
  <listitem>
              <para>
                <emphasis role="bold">Cooperative mechanism 1</emphasis> 
 is based on the theory that cross bridge formation between actin and myosin increases the affinity of troponin for Ca<superscript>2+</superscript>.</para>
            </listitem>
  <listitem>
              <para>
                <emphasis role="bold">Cooperative mechanism 2</emphasis> 
 assumes that the binding of a cross bridge increases the rate of formation of neighbouring cross bridges and that multiple cross bridges can actin activation even in the absence of Ca<superscript>2+</superscript>.</para>
            </listitem>
  <listitem>
              <para>
                <emphasis role="bold">Cooperative mechanism 2</emphasis> 
 simulates end-to-end interactions between adjacent troponin and tropomyosin.</para>
            </listitem>
</itemizedlist>   
</para>

<para>
Comparison of putative cooperative mechanisms in cardiac muscle: length dependence and dynamic responses,J. Jeremy Rice, Raimond L. Winslow and William C. Hunter, 1999, <emphasis>American Journal of Physiology</emphasis>, 276, H1734-H1754. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=10330260&amp;dopt=Abstract">PubMed ID: 10330260</ulink>
</para>

<informalfigure float="0" id="fig_reaction_diagrams">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>reaction_diagrams</title>
    </objectinfo>
    <imagedata fileref="rice_hunter_1999.png"/>
  </imageobject>
</mediaobject>
<caption>State diagrams for the five models of isometric force generation in cardiac muscle. T represents tropomyosin, TCa is Ca<superscript>2+</superscript> bound tropomyosin, N0, N1, P0 and P1 are the non-permissive and permissive tropomyosin states.</caption>
</informalfigure>

<para>
All the models are similar in that they are structured around a functional unit of troponin, tropomyosin and actin.  Tropomyosin can exist in four states, two permissive or two non-permissive (referring to whether or not the actin sites are available for binding to myosin and hence cross bridge formation).  Depending on the model, one or more cross bridges exist, and these are either weakly-bound (non-force generating) or strongly bound (force generating).
</para>

<para>
The paper (cited below) tests the behaviours of the five models of force generation in cardiac myocytes.  The first two models provide a baseline of performance for comparison.  Models 3 to 5 are developed to incorporate more cooperative mechanisms.  From the results of these simulations, which were compared to and consistent with experimental data, it is hypothesised that multiple mechanisms of cooperativity may coexist and contribute to the responses of cardiac muscle.
</para>

</sect1>
</article>
</documentation>
  
  
  
  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="micrometre">
    <unit units="metre" prefix="micro"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="second_order_rate_constant">
    <unit units="micromolar" exponent="-1"/>
    <unit units="second" exponent="-1"/>
  </units>
  
 
  <component name="environment">
    <variable units="second" public_interface="out" name="time"/>
  </component>
  
   
  <component cmeta:id="T" name="T">
    <variable units="dimensionless" public_interface="out" name="T"/>
    
    <variable units="dimensionless" public_interface="in" name="TCa"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> T </ci> 
        <apply>
          <minus/>
          <cn cellml:units="dimensionless"> 1.0 </cn> 
          <ci> TCa </ci>
        </apply>  
      </apply> 
    </math> 
  </component>

  <component cmeta:id="TCa" name="TCa">
    <variable units="dimensionless" public_interface="out" name="TCa"/>
    
    <variable units="second_order_rate_constant" public_interface="in" name="kon"/>
    <variable units="first_order_rate_constant" public_interface="in" name="koff"/>
    <variable units="dimensionless" public_interface="in" name="T"/>
    <variable units="micromolar" public_interface="in" name="Ca"/>  
    <variable units="second" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> TCa </ci> 
        </apply>
        <apply>
          <minus/> 
          <apply>
            <times/>
            <ci> kon </ci>
            <ci> Ca </ci> 
            <ci> T </ci>
          </apply>
          <apply>
            <times/>  
            <ci> koff </ci>
            <ci> TCa </ci>
          </apply>
        </apply>  
      </apply> 
    </math> 
  </component>

  <component cmeta:id="N0" name="N0">
    <variable units="dimensionless" public_interface="out" name="N0"/>
     
    <variable units="dimensionless" public_interface="in" name="P0"/>
    <variable units="dimensionless" public_interface="in" name="N1"/>
    <variable units="dimensionless" public_interface="in" name="P1"/>  
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> N0 </ci> 
        <apply>
          <minus/>
          <cn cellml:units="dimensionless"> 1.0 </cn> 
          <apply>
            <plus/>
            <ci> P0 </ci>
            <ci> N1 </ci>
            <ci> P1 </ci>
          </apply>
        </apply>  
      </apply> 
    </math> 
  </component>
  
  <component cmeta:id="N1" name="N1">
    <variable units="dimensionless" public_interface="out" name="N1"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="k1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k1_"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g10"/>
    <variable units="dimensionless" public_interface="in" name="N0"/>
    <variable units="dimensionless" public_interface="in" name="P1"/>
    <variable units="second" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> N1 </ci> 
        </apply>
        <apply>
          <minus/> 
          <apply>
            <times/>
            <ci> k1_ </ci>
            <ci> P1 </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>  
              <ci> g10 </ci>
              <ci> N1 </ci> 
            </apply>
            <apply>
              <times/>  
              <ci> k1 </ci>
              <ci> N1 </ci>
            </apply>
          </apply>
        </apply>  
      </apply> 
    </math> 
  </component>
   
  <component cmeta:id="P0" name="P0">
    <variable units="dimensionless" public_interface="out" name="P0"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="k1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k1_"/>
    <variable units="first_order_rate_constant" public_interface="in" name="f01"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g10"/>
    <variable units="dimensionless" public_interface="in" name="N0"/>
    <variable units="dimensionless" public_interface="in" name="P1"/>
    <variable units="second" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> P0 </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/> 
            <apply>
              <times/>
              <ci> k1 </ci>
              <ci> N0 </ci>
            </apply>
            <apply>
              <times/>
              <ci> g10 </ci>
              <ci> P1 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>  
              <ci> k1_ </ci>
              <ci> P0 </ci> 
            </apply>
            <apply>
              <times/>  
              <ci> f01 </ci>
              <ci> P0 </ci>
            </apply>
          </apply>
        </apply>  
      </apply> 
    </math> 
  </component>
  
  <component cmeta:id="P1" name="P1">
    <variable units="dimensionless" public_interface="out" name="P1"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="k1"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k1_"/>
    <variable units="first_order_rate_constant" public_interface="in" name="f01"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g10"/>
    <variable units="first_order_rate_constant" public_interface="in" name="f12"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g21"/>
    <variable units="dimensionless" public_interface="in" name="N1"/>
    <variable units="dimensionless" public_interface="in" name="P2"/>
    <variable units="dimensionless" public_interface="in" name="P0"/>
    <variable units="second" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> P1 </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/> 
            <apply>
              <times/>
              <ci> k1 </ci>
              <ci> N1 </ci>
            </apply>
            <apply>
              <times/>
              <ci> f01 </ci>
              <ci> P0 </ci>
            </apply>
            <apply>
              <times/>
              <ci> g21 </ci>
              <ci> P2 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>  
              <ci> k1_ </ci>
              <ci> P1 </ci> 
            </apply>
            <apply>
              <times/>  
              <ci> g10 </ci>
              <ci> P1 </ci>
            </apply>
            <apply>
              <times/>  
              <ci> f12 </ci>
              <ci> P1 </ci>
            </apply>
          </apply>
        </apply>  
      </apply> 
    </math> 
  </component>

  <component cmeta:id="P2" name="P2">
    <variable units="dimensionless" public_interface="out" name="P2"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="f23"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g32"/>
    <variable units="first_order_rate_constant" public_interface="in" name="f12"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g21"/>
    <variable units="dimensionless" public_interface="in" name="P1"/>
    <variable units="dimensionless" public_interface="in" name="P3"/>
    <variable units="second" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> P2 </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/> 
            <apply>
              <times/>
              <ci> f12 </ci>
              <ci> P1 </ci>
            </apply>
            <apply>
              <times/>
              <ci> g32 </ci>
              <ci> P3 </ci>
            </apply>
          </apply>
          <apply>
            <plus/> 
            <apply>
              <times/>
              <ci> f23 </ci>
              <ci> P2 </ci>
            </apply>
            <apply>
              <times/>
              <ci> g21 </ci>
              <ci> P2 </ci>
            </apply>
          </apply>
        </apply>  
      </apply> 
    </math> 
  </component>

  <component cmeta:id="P3" name="P3">
    <variable units="dimensionless" public_interface="out" name="P3"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="f23"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g32"/>
    <variable units="dimensionless" public_interface="in" name="P2"/>
    <variable units="second" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> P3 </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> f23 </ci>
            <ci> P2 </ci>
          </apply>
          <apply>
            <times/>
            <ci> g32 </ci>
            <ci> P3 </ci>
          </apply>
        </apply>  
      </apply> 
    </math> 
  </component>
   
  <component name="constants">
      <variable units="micromolar" public_interface="out" name="Ca" initial_value="1.0"/>
      <variable units="second_order_rate_constant" public_interface="out" name="kon" initial_value="20.0"/>
      <variable units="first_order_rate_constant" public_interface="out" name="koff" initial_value="40.0"/>
      <variable units="first_order_rate_constant" public_interface="out" name="k1"/>
      <variable units="first_order_rate_constant" public_interface="out" name="k1_" initial_value="45.0"/>
      <variable units="first_order_rate_constant" public_interface="out" name="f01"/>
      <variable units="first_order_rate_constant" public_interface="out" name="g10"/>
      <variable units="first_order_rate_constant" public_interface="out" name="f12"/>
      <variable units="first_order_rate_constant" public_interface="out" name="g21"/>
      <variable units="first_order_rate_constant" public_interface="out" name="f23"/>
      <variable units="first_order_rate_constant" public_interface="out" name="g32"/>
    
      <variable units="first_order_rate_constant" name="f" initial_value="10.0"/>
      <variable units="first_order_rate_constant" name="g" initial_value="30.0"/>
      <variable units="first_order_rate_constant" name="gSL"/>
      <variable units="micromolar" name="K_1_2"/>
      <variable units="micromolar" name="K_Ca" initial_value="1.47"/>
      <variable units="dimensionless" name="N"/>
      <variable units="micrometre" name="SL" initial_value="1.7"/>
      <variable units="micrometre" name="SL_norm"/>
     
      <variable units="dimensionless" public_interface="in" name="TCa"/>

      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1">
        <apply id="f01_calculation">
        <eq/>
          <ci> f01 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> f </ci>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1a">  
      <apply id="f12_calculation">
        <eq/>
          <ci> f12 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 10.0 </cn>
            <ci> f </ci>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1b">  
        <apply id="f23_calculation">
        <eq/>
          <ci> f23 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 7.0 </cn>
            <ci> f </ci>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1c">  
        <apply id="gSL_calculation">
        <eq/>
          <ci> gSL </ci>
          <apply>
          <times/>
            <ci> g </ci>
            <apply>
            <minus/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <apply>
              <power/>
                <ci> SL_norm </ci>
                <cn cellml:units="dimensionless"> 1.6 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1d">  
        <apply id="g10_calculation">
        <eq/>
          <ci> g10 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> gSL </ci>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1e">  
        <apply id="g21_calculation">
        <eq/>
          <ci> g21 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <ci> gSL </ci>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1f">    
        <apply id="g32_calculation">
        <eq/>
          <ci> g32 </ci>
          <apply>
          <times/>
            <cn cellml:units="dimensionless"> 3.0 </cn>
            <ci> gSL </ci>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1g">  
        <apply id="k1_calculation">
        <eq/>
          <ci> k1 </ci>
          <apply>
          <times/>
            <ci> k1_ </ci>
            <apply>
            <power/>
              <apply>
              <divide/>
                <ci> TCa </ci>
                <ci> K_1_2 </ci>
              </apply>
              <ci> N </ci>
            </apply>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1h">  
        <apply id="K_1_2_calculation">
        <eq/>
          <ci> K_1_2 </ci>
          <apply>
          <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
            <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
              <divide/>
                <ci> K_Ca </ci>
                <apply>
                <minus/>
                  <cn cellml:units="micromolar"> 1.5 </cn>
                  <apply>
                  <times/>
                    <ci> SL_norm </ci>
                    <cn cellml:units="micromolar"> 1.0 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1i">       
        <apply id="N_calculation">
        <eq/>
          <ci> N </ci>
          <apply>
          <plus/>
            <cn cellml:units="dimensionless"> 5.0 </cn>
            <apply>
            <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> SL_norm </ci>
            </apply>
          </apply>
        </apply>
      </math>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1j">          
        <apply id="SL_norm_calculation">
        <eq/>
          <ci> SL_norm </ci>
          <apply>
          <divide/>
            <apply>
            <minus/>
              <ci> SL </ci>
              <cn cellml:units="micrometre"> 1.7 </cn>
            </apply>
            <apply>
            <minus/>
              <cn cellml:units="micrometre"> 2.3 </cn>
              <cn cellml:units="micrometre"> 1.7 </cn>
            </apply>
          </apply>
        </apply>
      </math>        
   </component>
   
   <component name="force">
    <variable units="dimensionless" name="F"/>
    <variable units="dimensionless" name="alpha" initial_value="1.0"/>
    <variable units="dimensionless" name="Fmax"/>
    <variable units="dimensionless" name="P1_max"/>
    <variable units="dimensionless" name="P2_max"/>
    <variable units="dimensionless" name="P3_max"/>
    
    <variable units="dimensionless" public_interface="in" name="P1"/>
    <variable units="dimensionless" public_interface="in" name="N1"/>
    <variable units="dimensionless" public_interface="in" name="P2"/>
    <variable units="dimensionless" public_interface="in" name="P3"/>
    <variable units="first_order_rate_constant" public_interface="in" name="f01"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g10"/>
    <variable units="first_order_rate_constant" public_interface="in" name="f12"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g21"/>
    <variable units="first_order_rate_constant" public_interface="in" name="f23"/>
    <variable units="first_order_rate_constant" public_interface="in" name="g32"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="F_calculation">
        <eq/>
          <ci> F </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <ci> alpha </ci>
              <apply>
              <plus/>
                <ci> P1 </ci>
                <ci> N1 </ci>
                <apply>
                <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> P2 </ci>
                </apply>
                <apply>
                <times/>
                  <cn cellml:units="dimensionless"> 3.0 </cn>
                  <ci> P3 </ci>
                </apply>
              </apply>
            </apply>
            <ci> Fmax </ci>
          </apply>
        </apply>
        
        <apply id="Fmax_calculation">
        <eq/>
          <ci> Fmax </ci>
          <apply>
          <plus/>
            <ci> P1_max </ci>
            <apply>
            <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> P2_max </ci>
            </apply>
            <apply>
            <times/>
              <cn cellml:units="dimensionless"> 3.0 </cn>
              <ci> P3_max </ci>
            </apply>
          </apply>
        </apply>
        
        <apply id="P1_max_calculation">
        <eq/>
          <ci> P1_max </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <ci> f01 </ci>
              <ci> g21 </ci>
              <ci> g32 </ci>
            </apply>
            <apply>
            <plus/>
              <apply>
              <times/>
                <ci> g10 </ci>
                <ci> g21 </ci>
                <ci> g32 </ci>
              </apply>
              <apply>
              <times/>
                <ci> f01 </ci>
                <ci> g21 </ci>
                <ci> g32 </ci>
              </apply>
              <apply>
              <times/>
                <ci> f01 </ci>
                <ci> f12 </ci>
                <ci> g32 </ci>
              </apply>
              <apply>
              <times/>
                <ci> f01 </ci>
                <ci> f12 </ci>
                <ci> f23 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
        
        <apply id="P2_max_calculation">
        <eq/>
          <ci> P2_max </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <ci> f01 </ci>
              <ci> f12 </ci>
              <ci> g32 </ci>
            </apply>
            <apply>
            <plus/>
              <apply>
              <times/>
                <ci> g10 </ci>
                <ci> g21 </ci>
                <ci> g32 </ci>
              </apply>
              <apply>
              <times/>
                <ci> f01 </ci>
                <ci> g21 </ci>
                <ci> g32 </ci>
              </apply>
              <apply>
              <times/>
                <ci> f01 </ci>
                <ci> f12 </ci>
                <ci> g32 </ci>
              </apply>
              <apply>
              <times/>
                <ci> f01 </ci>
                <ci> f12 </ci>
                <ci> f23 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
        
        <apply id="P3_max_calculation">
        <eq/>
          <ci> P3_max </ci>
          <apply>
          <divide/>
            <apply>
            <times/>
              <ci> f01 </ci>
              <ci> f12 </ci>
              <ci> f23 </ci>
            </apply>
            <apply>
            <plus/>
              <apply>
              <times/>
                <ci> g10 </ci>
                <ci> g21 </ci>
                <ci> g32 </ci>
              </apply>
              <apply>
              <times/>
                <ci> f01 </ci>
                <ci> g21 </ci>
                <ci> g32 </ci>
              </apply>
              <apply>
              <times/>
                <ci> f01 </ci>
                <ci> f12 </ci>
                <ci> g32 </ci>
              </apply>
              <apply>
              <times/>
                <ci> f01 </ci>
                <ci> f12 </ci>
                <ci> f23 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>          
   </component>
  
   
   <connection>
      <map_components component_2="environment" component_1="N1"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="P0"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>

   <connection>
      <map_components component_2="environment" component_1="TCa"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="P1"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
  
   <connection>
      <map_components component_2="environment" component_1="P2"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>

   <connection>
      <map_components component_2="environment" component_1="P3"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>

   <connection>
      <map_components component_2="TCa" component_1="T"/>
      <map_variables variable_2="T" variable_1="T"/>
      <map_variables variable_2="TCa" variable_1="TCa"/>
   </connection>

   <connection>
      <map_components component_2="P0" component_1="N0"/>
      <map_variables variable_2="N0" variable_1="N0"/>
      <map_variables variable_2="P0" variable_1="P0"/>
   </connection>
  
   <connection>
      <map_components component_2="N1" component_1="N0"/>
      <map_variables variable_2="N0" variable_1="N0"/>
      <map_variables variable_2="N1" variable_1="N1"/>
   </connection>

   <connection>
      <map_components component_2="P1" component_1="N0"/>
      <map_variables variable_2="P1" variable_1="P1"/>
   </connection>
  
   <connection>
      <map_components component_2="P1" component_1="N1"/>
      <map_variables variable_2="N1" variable_1="N1"/>
      <map_variables variable_2="P1" variable_1="P1"/>
   </connection>

   <connection>
      <map_components component_2="P0" component_1="P1"/>
      <map_variables variable_2="P1" variable_1="P1"/>
      <map_variables variable_2="P0" variable_1="P0"/>
   </connection>

   <connection>
      <map_components component_2="P2" component_1="P1"/>
      <map_variables variable_2="P1" variable_1="P1"/>
      <map_variables variable_2="P2" variable_1="P2"/>
   </connection>

   <connection>
      <map_components component_2="P2" component_1="P3"/>
      <map_variables variable_2="P3" variable_1="P3"/>
      <map_variables variable_2="P2" variable_1="P2"/>
   </connection>

   <connection>
      <map_components component_2="constants" component_1="TCa"/>
      <map_variables variable_2="kon" variable_1="kon"/>
      <map_variables variable_2="koff" variable_1="koff"/>
      <map_variables variable_2="TCa" variable_1="TCa"/>
      <map_variables variable_2="Ca" variable_1="Ca"/>
   </connection>
   
   <connection>
      <map_components component_2="constants" component_1="N1"/>
      <map_variables variable_2="k1" variable_1="k1"/>
      <map_variables variable_2="k1_" variable_1="k1_"/>
      <map_variables variable_2="g10" variable_1="g10"/>
   </connection>

   <connection>
      <map_components component_2="constants" component_1="P0"/>
      <map_variables variable_2="k1" variable_1="k1"/>
      <map_variables variable_2="k1_" variable_1="k1_"/>
      <map_variables variable_2="f01" variable_1="f01"/>
      <map_variables variable_2="g10" variable_1="g10"/>
   </connection>

   <connection>
      <map_components component_2="constants" component_1="P1"/>
      <map_variables variable_2="k1" variable_1="k1"/>
      <map_variables variable_2="k1_" variable_1="k1_"/>
      <map_variables variable_2="f01" variable_1="f01"/>
      <map_variables variable_2="g10" variable_1="g10"/>
      <map_variables variable_2="f12" variable_1="f12"/>
      <map_variables variable_2="g21" variable_1="g21"/>
   </connection>

   <connection>
      <map_components component_2="constants" component_1="P2"/>
      <map_variables variable_2="f23" variable_1="f23"/>
      <map_variables variable_2="g32" variable_1="g32"/>
      <map_variables variable_2="f12" variable_1="f12"/>
      <map_variables variable_2="g21" variable_1="g21"/>
   </connection>

   <connection>
      <map_components component_2="constants" component_1="P3"/>
      <map_variables variable_2="f23" variable_1="f23"/>
      <map_variables variable_2="g32" variable_1="g32"/>
   </connection>

   <connection>
      <map_components component_2="P1" component_1="force"/>
      <map_variables variable_2="P1" variable_1="P1"/>
   </connection>

   <connection>
      <map_components component_2="N1" component_1="force"/>
      <map_variables variable_2="N1" variable_1="N1"/>
   </connection>

   <connection>
      <map_components component_2="P2" component_1="force"/>
      <map_variables variable_2="P2" variable_1="P2"/>
   </connection>

   <connection>
      <map_components component_2="P3" component_1="force"/>
      <map_variables variable_2="P3" variable_1="P3"/>
   </connection>

   <connection>
      <map_components component_2="constants" component_1="force"/>
      <map_variables variable_2="f01" variable_1="f01"/>
      <map_variables variable_2="g10" variable_1="g10"/>
      <map_variables variable_2="f12" variable_1="f12"/>
      <map_variables variable_2="g21" variable_1="g21"/>
      <map_variables variable_2="f23" variable_1="f23"/>
      <map_variables variable_2="g32" variable_1="g32"/>
   </connection>

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#419c6234-a164-4559-9365-2d05f6f4cc48">
    <rdf:li>Cardiac Myocyte</rdf:li>
    <rdf:li>cooperative mechanisms</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>myofilament mechanics</rdf:li>
    <rdf:li>cardiac</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#8fd8cdd0-cee5-4505-bb3d-84337cc74c04">
    <rdf:li rdf:resource="rdf:#898dbb1d-83fa-494c-8886-3a8bcb8c9670"/>
    <rdf:li rdf:resource="rdf:#f8c58b06-6921-4093-ab9b-c20dee87a644"/>
    <rdf:li rdf:resource="rdf:#8d20d160-1116-4b99-b303-12934f83d569"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#cbd25d97-3e4a-4f10-94d5-b4b00d9ffc12">
    <dc:creator rdf:resource="rdf:#8692b76e-e126-41b3-a769-e1ca78835e87"/>
    <rdf:value>This is the CellML description of Rice et al's 4th model of isometric &#13;
          force generation in cardiac myofilaments taken from their 1999 paper.&#13;
          All five mathematical models published in this paper were constructed &#13;
          by assuming different subsets of three putative cooperative &#13;
          mechanisms.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#86a1aa8d-912f-4b65-bfea-7895e031c0ae">
    <bqs:Pubmed_id>10330260</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#279ac64f-629a-41a5-a6d3-bdb76ac33544"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#72cee108-c761-4be2-b5c9-4b8798862d2a"/>
    <dcterms:created rdf:resource="rdf:#441cdfa4-85db-4dcd-a0b5-927519fb2632"/>
    <dc:creator rdf:resource="rdf:#eba30a70-9384-481d-b6f1-33be64c3f657"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ee0b9aff-3952-4e29-b9e4-749af1f5fc2d">
    <vCard:Given>William</vCard:Given>
    <vCard:Family>Hunter</vCard:Family>
    <vCard:Other>C</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f8c58b06-6921-4093-ab9b-c20dee87a644">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#a3c5664e-1645-4140-a07d-d8acd8e4e6ec"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#279ac64f-629a-41a5-a6d3-bdb76ac33544">
    <dc:creator rdf:resource="rdf:#8fd8cdd0-cee5-4505-bb3d-84337cc74c04"/>
    <dc:title>Comparison of putative cooperative mechanisms in cardiac muscle: length dependence and dynamic responses (model 4)</dc:title>
    <bqs:volume>276</bqs:volume>
    <bqs:first_page>H1734</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#943c0f23-12f9-4e06-ac28-00dadd8359b6"/>
    <dcterms:issued rdf:resource="rdf:#0fc41fa4-63ec-4e6b-8fbe-9c459d9cb51b"/>
    <bqs:last_page>H1754</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a3c5664e-1645-4140-a07d-d8acd8e4e6ec">
    <vCard:Given>Raimond</vCard:Given>
    <vCard:Family>Winslow</vCard:Family>
    <vCard:Other>L</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8d20d160-1116-4b99-b303-12934f83d569">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#ee0b9aff-3952-4e29-b9e4-749af1f5fc2d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#154f39f0-4888-4510-9f63-73e6575bf031">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#eba30a70-9384-481d-b6f1-33be64c3f657">
    <vCard:ORG rdf:resource="rdf:#91a86ecc-2fbd-4dce-96d9-8dd70a8e5584"/>
    <vCard:EMAIL rdf:resource="rdf:#c1d7d62d-7c75-4041-9fbb-c5c260bc0e5d"/>
    <vCard:N rdf:resource="rdf:#154f39f0-4888-4510-9f63-73e6575bf031"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6eccd62d-5848-430f-852b-206117e10870">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Rice</vCard:Family>
    <vCard:Other>Jeremy</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#943c0f23-12f9-4e06-ac28-00dadd8359b6">
    <dc:title>American Journal of Physiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#441cdfa4-85db-4dcd-a0b5-927519fb2632">
    <dcterms:W3CDTF>2007-05-24T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#91a86ecc-2fbd-4dce-96d9-8dd70a8e5584">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0fc41fa4-63ec-4e6b-8fbe-9c459d9cb51b">
    <dcterms:W3CDTF>1999-05</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#146dcc9c-5d72-4a98-95d4-3abfd6f4787f">
    <dc:subject rdf:resource="rdf:#d2aa11e6-a8ef-4592-a956-fa23943824d8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fc199485-29ad-40fa-895d-5e07292aea6c">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c1d7d62d-7c75-4041-9fbb-c5c260bc0e5d">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#rice_model4_1999_version02">
    <dc:title>
        The Rice et al. 1999 4th model of isometric force generation in cardiac 
        myofilaments</dc:title>
    <cmeta:bio_entity>Cardiac Myocyte</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#cbd25d97-3e4a-4f10-94d5-b4b00d9ffc12"/>
    <bqs:reference rdf:resource="rdf:#146dcc9c-5d72-4a98-95d4-3abfd6f4787f"/>
    <bqs:reference rdf:resource="rdf:#86a1aa8d-912f-4b65-bfea-7895e031c0ae"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#72cee108-c761-4be2-b5c9-4b8798862d2a">
    <dc:creator rdf:resource="rdf:#fc199485-29ad-40fa-895d-5e07292aea6c"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#898dbb1d-83fa-494c-8886-3a8bcb8c9670">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#6eccd62d-5848-430f-852b-206117e10870"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d2aa11e6-a8ef-4592-a956-fa23943824d8">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#419c6234-a164-4559-9365-2d05f6f4cc48"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8692b76e-e126-41b3-a769-e1ca78835e87">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
</rdf:RDF>
</model>