Rendering of the source text

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  rice_model5_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 5th
               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_model5_1999_version02" name="rice_winslow_hunter_1999_version05">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Cooperative Mechanisms in Cardiac Muscle (model 5)</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>
<ulink url="http://ajpheart.physiology.org/cgi/content/abstract/276/5/H1734">Comparison of putative cooperative mechanisms in cardiac muscle: length dependence and dynamic responses</ulink> J. Jeremy Rice, Raimond L. Winslow and William C. Hunter, 1999, <ulink url="http://ajpheart.physiology.org/">
            <emphasis>American Journal of Physiology</emphasis>
          </ulink>, 276, H1734-H1754.  (<ulink url="http://ajpheart.physiology.org/cgi/content/full/276/5/H1734">Full text</ulink> and <ulink url="http://ajpheart.physiology.org/cgi/reprint/276/5/H1734.pdf">PDF</ulink> versions of the article are available for Journal Members on the American Journal of Physiology website.)  <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 name="T" cmeta:id="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 name="TCa" cmeta:id="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 name="N0" cmeta:id="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 name="N1" cmeta:id="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 name="P0" cmeta:id="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 name="P1" cmeta:id="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 name="P2" cmeta:id="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 name="P3" cmeta:id="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_"/>
      <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="koff" initial_value="40.0"/>
      <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"/>
      <variable units="dimensionless" public_interface="in" name="F"/>

      <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.8 </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"> 2.6 </cn>
            <apply>
            <times/>
              <cn cellml:units="dimensionless"> 1.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>
      <math xmlns="http://www.w3.org/1998/Math/MathML" id="1k">          
        <apply id="koff__calculation">
        <eq/>
          <ci> koff_ </ci>
          <apply>
          <times/>
            <ci> koff </ci>
            <apply>
            <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
              <times/>
                <cn cellml:units="dimensionless"> 0.75 </cn>
                <ci> F </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </math>                
   </component>
   
   <component name="force">
    <variable units="dimensionless" public_interface="out" 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"/>
      <map_variables variable_2="F" variable_1="F"/>
   </connection>


<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#b4344afa-e226-4709-9a11-e0fa1cbf07d7">
    <rdf:li>cardiac myocyte</rdf:li>
    <rdf:li>myofilament mechanics</rdf:li>
    <rdf:li>Cardiac Myocyte</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>cooperative mechanisms</rdf:li>
    <rdf:li>cardiac</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#0d6eabe8-dea3-4f69-a621-c296f03732d8">
    <rdf:li rdf:resource="rdf:#8a0581d1-15ad-4557-bd8d-67b664e92c41"/>
    <rdf:li rdf:resource="rdf:#5fd9ed83-1134-4437-9922-33cbaca843fd"/>
    <rdf:li rdf:resource="rdf:#7b4fc7ca-3e9a-4ec7-bdb6-974a1e13edf7"/>
  </rdf:Seq>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#9b109cc4-ec3c-4f2b-bc64-580d54a61090"/>
    <dcterms:created rdf:resource="rdf:#ebd4cf81-cf23-4ede-b115-4089d381dabc"/>
    <dc:creator rdf:resource="rdf:#e7333d38-59fb-4955-a399-a22e524c837c"/>
    <cmeta:modification rdf:resource="rdf:#42c6dc62-3694-45ad-afbc-b94cfba14388"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ebd4cf81-cf23-4ede-b115-4089d381dabc">
    <dcterms:W3CDTF>2007-05-24T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2c600273-2de6-49e4-bb1a-cef87f8a2a0d">
    <dc:creator rdf:resource="rdf:#c26a5939-4cc0-4696-993b-ab9918cfeee8"/>
    <rdf:value>This is the CellML description of Rice et al's 5th 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:#0a03a20b-9c1f-4e61-b7d1-acb641b07946">
    <vCard:Given>Raimond</vCard:Given>
    <vCard:Family>Winslow</vCard:Family>
    <vCard:Other>L</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#abbb3a6c-2c17-4070-975f-bf8d264ba7a7">
    <dc:title>American Journal of Physiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="#rice_model5_1999_version02">
    <dc:title>
        The Rice et al. 1999 5th model of isometric force generation in cardiac 
        myofilaments</dc:title>
    <cmeta:bio_entity>Cardiac Myocyte</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#2c600273-2de6-49e4-bb1a-cef87f8a2a0d"/>
    <bqs:reference rdf:resource="rdf:#0b0156a3-d21c-4ae7-b34b-53b722143a5d"/>
    <bqs:reference rdf:resource="rdf:#ca15e454-1ce1-4357-8d64-541a8cd4d031"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8a0581d1-15ad-4557-bd8d-67b664e92c41">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#24d6156d-e8d8-4f00-bcb3-74efedfebdde"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5fd9ed83-1134-4437-9922-33cbaca843fd">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#0a03a20b-9c1f-4e61-b7d1-acb641b07946"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a92c9211-8ced-40e7-85f6-d1658cfeb44f">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#701b2a27-c4c4-41ba-a5bb-fdc45481add7">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0b0156a3-d21c-4ae7-b34b-53b722143a5d">
    <bqs:Pubmed_id>10330260</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#6a109e1a-b2da-484e-8a62-7215288ae0f5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#42c6dc62-3694-45ad-afbc-b94cfba14388">
    <dcterms:modified rdf:resource="rdf:#161612c8-4837-4c90-947d-5e4da72d9dbf"/>
    <rdf:value>The new version of this model has been re-coded to remove the reaction element and replace it with a simple MathML description of the model reaction kinetics. This is thought to be truer to the original publication, and information regarding the enzyme kinetics etc will later be added to the metadata through use of an ontology. &#13;
&#13;
The model runs in the PCEnv simulator and gives a nice curved output.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#ad5a2403-5d45-40cf-8024-be7a273af5bf"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9b109cc4-ec3c-4f2b-bc64-580d54a61090">
    <dc:creator rdf:resource="rdf:#f443fef9-d85b-4b6a-a0b6-dda7356d0b39"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c26a5939-4cc0-4696-993b-ab9918cfeee8">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ca15e454-1ce1-4357-8d64-541a8cd4d031">
    <dc:subject rdf:resource="rdf:#3d9dd3d4-5ed0-4b8c-9503-c5558ae2311b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e7333d38-59fb-4955-a399-a22e524c837c">
    <vCard:ORG rdf:resource="rdf:#d8a3e810-3dda-4be5-add6-c88a78c4808e"/>
    <vCard:EMAIL rdf:resource="rdf:#a92c9211-8ced-40e7-85f6-d1658cfeb44f"/>
    <vCard:N rdf:resource="rdf:#07900db7-3f8d-4739-a2b2-4674b17e44e5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7b4fc7ca-3e9a-4ec7-bdb6-974a1e13edf7">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#12ab8b6a-ef76-4930-91e7-778ffe282e2e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ad5a2403-5d45-40cf-8024-be7a273af5bf">
    <vCard:N rdf:resource="rdf:#701b2a27-c4c4-41ba-a5bb-fdc45481add7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3d9dd3d4-5ed0-4b8c-9503-c5558ae2311b">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#b4344afa-e226-4709-9a11-e0fa1cbf07d7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#24d6156d-e8d8-4f00-bcb3-74efedfebdde">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Rice</vCard:Family>
    <vCard:Other>Jeremy</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#12ab8b6a-ef76-4930-91e7-778ffe282e2e">
    <vCard:Given>William</vCard:Given>
    <vCard:Family>Hunter</vCard:Family>
    <vCard:Other>C</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d8a3e810-3dda-4be5-add6-c88a78c4808e">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6a109e1a-b2da-484e-8a62-7215288ae0f5">
    <dc:creator rdf:resource="rdf:#0d6eabe8-dea3-4f69-a621-c296f03732d8"/>
    <dc:title>Comparison of putative cooperative mechanisms in cardiac muscle: length dependence and dynamic responses (model 5)</dc:title>
    <bqs:volume>276</bqs:volume>
    <bqs:first_page>H1734</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#abbb3a6c-2c17-4070-975f-bf8d264ba7a7"/>
    <dcterms:issued rdf:resource="rdf:#ff7af07c-0b24-41f5-9436-c89098d6fdfa"/>
    <bqs:last_page>H1754</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ff7af07c-0b24-41f5-9436-c89098d6fdfa">
    <dcterms:W3CDTF>1999-05</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#07900db7-3f8d-4739-a2b2-4674b17e44e5">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#161612c8-4837-4c90-947d-5e4da72d9dbf">
    <dcterms:W3CDTF>2007-06-05T09:58:57+12:00</dcterms:W3CDTF>
  </rdf:Description>
</rdf:RDF>
</model>