Location: Rice, Wang, Bers, Detombe, 2008 @ d90a594a7a24 / rice_wang_bers_detombe_2008.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2008-07-07 22:31:04+12:00
Desc:
committing version01 of rice_wang_bers_detombe_2008
Permanent Source URI:
https://models.physiomeproject.org/workspace/rice_wang_bers_detombe_2008/rawfile/d90a594a7a249398295acbbc7d38bfc5e27069d5/rice_wang_bers_detombe_2008.cellml

<?xml version='1.0' encoding='utf-8'?>
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="rice_model_2008" name="rice_model_2008">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations</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>
        This cellML model is known to run in both PCEnv and COR.  The units have been checked and they are consistent.
          </para>
				<para>
Abstract: We develop a point model of the cardiac myofilament (MF) to simulate wide variety of experimental muscle characterizations including Force-Ca relations and twitches under isometric, isosarcometric, isotonic and auxotonic conditions. Complex MF behaviors are difficult to model because spatial interactions cannot be directly implemented as ordinary differential equations (ODEs). We therefore allow phenomenological approximations with careful consideration to the relationships with the underlying biophysical mechanisms. We describe new formulations that avoid mean-field approximations found in most existing MF models. To increase the scope and applicability of the model, we include length- and temperature-dependent effects that play important roles in MF response. We have also included a representation of passive restoring forces to simulate isolated cell shortening protocols. Possessing both computational efficiency and the ability to simulate a wide variety of muscle responses, the MF representation is well-suited for coupling to existing cardiac cell models of electrophysiology and Ca-handling mechanisms. To illustrate this suitability, the MF model is coupled to the Chicago rabbit cardiomyocyte model. The combined model generates realistic appearing action potentials, intracellular Ca transients, and cell shortening signals. The combined model also demonstrates that the feedback effects of force on Ca binding to troponin can modify the cystolic Ca transient.
				</para>
				<para>
Note that the current CellMl model describes the rat model of myofilamentcontraction. One of the goals of the published work by Rice et al. was to develop a model of the myofilaments which could be coupled to existing models of cardiac electrophysiology and Ca-handling mechanisms.  Similarly, it is our intention to be able to couple the CellML model of myofilament contration with cardiac models of electrophysiology and calcium change. To achieve this we will first have to modify the rat myofilament model to convert it to a rabbit myofilament model. This can be achieved by decreasing the transition rates in the crossbridge cycle by a factor of 5 to simulate the changes in myosin isoforms (rat is predominately V1 while rabbit is V3) and species-based changes in rate.  We can then couple the model to the <ulink url="http://www.cellml.org/models/shannon_wang_puglisi_weber_bers_2004_version04">Chicago rabbit ventricular myocyte model</ulink> by using the cytosolic Ca concentration ([Ca]c) from the Chicago model as the input to the myofilament model.  Other minor modifications are also required, as described in the paper below:
				</para>
				<para>
					<ulink url="http://www.biophysj.org/cgi/content/abstract/biophysj.107.119487v1">Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations</ulink>, John Jeremy Rice, Fei Wang, Donald M. Bers and Pieter P. de Tombe, 2008, <ulink url="http://www.biophysj.org/">
						<emphasis>Biophysical Journal</emphasis>
					</ulink>.  (A <ulink url="http://www.biophysj.org/cgi/rapidpdf/biophysj.107.119487v1">PDF</ulink> version of the article is available to journal subscribers on the <emphasis>Biophysical Journal</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=18234826&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 18234826</ulink>
				</para>
		<informalfigure float="0" id="fig_reaction_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>model diagram</title>
							</objectinfo>
							<imagedata fileref="rice_2008.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the model.</caption>
				</informalfigure> 
							</sect1>
		</article>
	</documentation>
	   <units base_units="yes" name="unit_normalised_force"/>
   <units name="unit_normalised_force_millisecond2_per_micrometre">
      <unit units="unit_normalised_force"/>
      <unit units="millisecond" exponent="2"/>
      <unit units="micrometre" exponent="-1"/>
   </units>
   <units name="unit_normalised_force_millisecond">
      <unit units="unit_normalised_force"/>
      <unit units="millisecond"/>
   </units>
   <units name="unit_normalised_force_millisecond_per_micrometre">
      <unit units="unit_normalised_force"/>
      <unit units="millisecond"/>
      <unit units="micrometre" exponent="-1"/>
   </units>
   <units name="unit_normalised_force_per_micrometre">
      <unit units="unit_normalised_force"/>
      <unit units="micrometre" exponent="-1"/>
   </units>
   <units name="micrometre">
      <unit units="metre" prefix="micro"/>
   </units>
   <units name="per_micrometre">
      <unit units="micrometre" exponent="-1"/>
   </units>
   <units name="micrometre_per_millisecond">
      <unit units="micrometre"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="millisecond_per_micrometre">
      <unit units="millisecond"/>
      <unit units="micrometre" exponent="-1"/>
   </units>
   <units name="millisecond">
      <unit units="second" prefix="milli"/>
   </units>
   <units name="millisecond2_per_micrometre">
      <unit units="millisecond" exponent="2"/>
      <unit units="micrometre" exponent="-1"/>
   </units>
   <units name="micromolar">
      <unit units="mole" prefix="micro"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="micromolar_per_millisecond">
      <unit units="micromolar"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="first_order_rate_constant">
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="second_order_rate_constant">
      <unit units="micromolar" exponent="-1"/>
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="millinewton_per_millimetre2">
      <unit units="newton" prefix="milli"/>
      <unit units="metre" prefix="milli" exponent="-2"/>
   </units>
   <units name="millinewton_micrometre_per_millimetre2">
      <unit units="newton" prefix="milli"/>
      <unit units="metre" prefix="micro"/>
      <unit units="metre" prefix="milli" exponent="-2"/>
   </units>
   <units name="millinewton_micrometre_per_millimetre2_per_millisecond">
      <unit units="newton" prefix="milli"/>
      <unit units="metre" prefix="micro"/>
      <unit units="metre" prefix="milli" exponent="-2"/>
      <unit units="second" prefix="milli" exponent="-1"/>
   </units>
   <component name="environment">
      <variable units="millisecond" public_interface="out" name="time"/>
   </component>
   <component name="sarcomere_geometry">
      <variable units="dimensionless" public_interface="out" name="SOVFThick"/>
      <variable units="dimensionless" public_interface="out" name="SOVFThin"/>
      <variable units="micrometre" name="sovr_ze"/>
      <variable units="micrometre" name="sovr_cle"/>
      <variable units="micrometre" name="len_sovr"/>
      <variable units="micrometre" public_interface="in" name="SLmax"/>
      <variable units="micrometre" public_interface="in" name="SLmin"/>
      <variable units="micrometre" public_interface="in" name="len_thin"/>
      <variable units="micrometre" public_interface="in" name="len_thick"/>
      <variable units="micrometre" public_interface="in" name="len_hbare"/>
      <variable units="micrometre" public_interface="in" name="SL"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>sovr_ze</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <ci>len_thick</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <lt/>
                     <apply>
                        <divide/>
                        <ci>len_thick</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <apply>
                        <divide/>
                        <ci>SL</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <ci>SL</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>sovr_cle</ci>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <apply>
                        <divide/>
                        <ci>SL</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <apply>
                        <minus/>
                        <ci>SL</ci>
                        <ci>len_thin</ci>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <apply>
                        <minus/>
                        <apply>
                           <divide/>
                           <ci>SL</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <minus/>
                           <ci>SL</ci>
                           <ci>len_thin</ci>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <ci>len_hbare</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <ci>len_hbare</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>len_sovr</ci>
            <apply>
               <minus/>
               <ci>sovr_ze</ci>
               <ci>sovr_cle</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>SOVFThick</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>len_sovr</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>len_thick</ci>
                  <ci>len_hbare</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>SOVFThin</ci>
            <apply>
               <divide/>
               <ci>len_sovr</ci>
               <ci>len_thin</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_binding_to_troponin_to_thin_filament_regulation">
      <variable units="dimensionless" public_interface="out" name="TRPNCaL" initial_value="0.5"/>
      <variable units="dimensionless" public_interface="out" name="TRPNCaH" initial_value="0.5"/>
      <variable units="first_order_rate_constant" public_interface="out" name="dTRPNCaL"/>
      <variable units="first_order_rate_constant" public_interface="out" name="dTRPNCaH"/>
      <variable units="first_order_rate_constant" public_interface="out" name="kn_pT"/>
      <variable units="first_order_rate_constant" public_interface="out" name="kp_nT"/>
      <variable units="second_order_rate_constant" name="konT"/>
      <variable units="first_order_rate_constant" name="koffLT"/>
      <variable units="first_order_rate_constant" name="koffHT"/>
      <variable units="dimensionless" name="Qkon" initial_value="1.5"/>
      <variable units="dimensionless" name="Qkoff" initial_value="1.3"/>
      <variable units="dimensionless" name="Qkn_p" initial_value="1.6"/>
      <variable units="dimensionless" name="Qkp_n" initial_value="1.6"/>
      <variable units="second_order_rate_constant" name="kon" initial_value="50e-3"/>
      <variable units="first_order_rate_constant" name="koffL" initial_value="250e-3"/>
      <variable units="first_order_rate_constant" name="koffH" initial_value="25e-3"/>
      <variable units="dimensionless" name="perm50" initial_value="0.5"/>
      <variable units="dimensionless" name="nperm" initial_value="15"/>
      <variable units="first_order_rate_constant" name="kn_p" initial_value="500e-3"/>
      <variable units="first_order_rate_constant" name="kp_n" initial_value="50e-3"/>
      <variable units="dimensionless" name="koffmod" initial_value="1"/>
      <variable units="dimensionless" name="Tropreg"/>
      <variable units="dimensionless" name="permtot"/>
      <variable units="dimensionless" name="inprmt"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="celsius" public_interface="in" name="TmpC"/>
      <variable units="micromolar" public_interface="in" name="Cai"/>
      <variable units="dimensionless" public_interface="in" name="SOVFThin"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>konT</ci>
            <apply>
               <times/>
               <ci>kon</ci>
               <apply>
                  <power/>
                  <ci>Qkon</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>koffLT</ci>
            <apply>
               <times/>
               <ci>koffL</ci>
               <ci>koffmod</ci>
               <apply>
                  <power/>
                  <ci>Qkoff</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>koffHT</ci>
            <apply>
               <times/>
               <ci>koffH</ci>
               <ci>koffmod</ci>
               <apply>
                  <power/>
                  <ci>Qkoff</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>kn_pT</ci>
            <apply>
               <times/>
               <ci>kn_p</ci>
               <ci>permtot</ci>
               <apply>
                  <power/>
                  <ci>Qkn_p</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>kp_nT</ci>
            <apply>
               <times/>
               <ci>kp_n</ci>
               <ci>inprmt</ci>
               <apply>
                  <power/>
                  <ci>Qkp_n</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>permtot</ci>
            <apply>
               <power/>
               <apply>
                  <abs/>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <ci>perm50</ci>
                              <ci>Tropreg</ci>
                           </apply>
                           <ci>nperm</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="dimensionless">0.5</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>inprmt</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>permtot</ci>
                  </apply>
                  <apply>
                     <lt/>
                     <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>permtot</ci>
                     </apply>
                     <cn cellml:units="dimensionless">100</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">100</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>TRPNCaL</ci>
            </apply>
            <ci>dTRPNCaL</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dTRPNCaL</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>konT</ci>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>TRPNCaL</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>koffLT</ci>
                  <ci>TRPNCaL</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>TRPNCaH</ci>
            </apply>
            <ci>dTRPNCaH</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dTRPNCaH</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>konT</ci>
                  <ci>Cai</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>TRPNCaH</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>koffHT</ci>
                  <ci>TRPNCaH</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Tropreg</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>SOVFThin</ci>
                  </apply>
                  <ci>TRPNCaL</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>SOVFThin</ci>
                  <ci>TRPNCaH</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="thin_filament_regulation_and_crossbridge_cycling_rates">
      <variable units="first_order_rate_constant" public_interface="out" name="fappT"/>
      <variable units="first_order_rate_constant" public_interface="out" name="gappT"/>
      <variable units="first_order_rate_constant" public_interface="out" name="hfT"/>
      <variable units="first_order_rate_constant" public_interface="out" name="hbT"/>
      <variable units="first_order_rate_constant" public_interface="out" name="gxbT"/>
      <variable units="first_order_rate_constant" public_interface="out" name="fapp" initial_value="500e-3"/>
      <variable units="first_order_rate_constant" public_interface="out" name="gapp" initial_value="70e-3"/>
      <variable units="first_order_rate_constant" public_interface="out" name="hf" initial_value="2000e-3"/>
      <variable units="first_order_rate_constant" public_interface="out" name="hb" initial_value="400e-3"/>
      <variable units="first_order_rate_constant" public_interface="out" name="gxb" initial_value="70e-3"/>
      <variable units="dimensionless" name="gslmod" initial_value="6"/>
      <variable units="dimensionless" name="hfmd"/>
      <variable units="dimensionless" name="hbmd"/>
      <variable units="dimensionless" name="hfmdc" initial_value="5"/>
      <variable units="dimensionless" name="hbmdc" initial_value="0"/>
      <variable units="dimensionless" name="sigmap" initial_value="8"/>
      <variable units="dimensionless" name="sigman" initial_value="1"/>
      <variable units="dimensionless" name="xbmodsp" initial_value="1"/>
      <variable units="dimensionless" name="Qfapp" initial_value="6.25"/>
      <variable units="dimensionless" name="Qgapp" initial_value="2.5"/>
      <variable units="dimensionless" name="Qhf" initial_value="6.25"/>
      <variable units="dimensionless" name="Qhb" initial_value="6.25"/>
      <variable units="dimensionless" name="Qgxb" initial_value="6.25"/>
      <variable units="dimensionless" name="gxbmd"/>
      <variable units="dimensionless" name="gapslmd"/>
      <variable units="micrometre" public_interface="in" name="x_0"/>
      <variable units="celsius" public_interface="in" name="TmpC"/>
      <variable units="dimensionless" public_interface="in" name="SOVFThick"/>
      <variable units="micrometre" public_interface="in" name="xXBpostr"/>
      <variable units="micrometre" public_interface="in" name="xXBprer"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>fappT</ci>
            <apply>
               <times/>
               <ci>fapp</ci>
               <ci>xbmodsp</ci>
               <apply>
                  <power/>
                  <ci>Qfapp</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>gapslmd</ci>
            <apply>
               <plus/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>SOVFThick</ci>
                  </apply>
                  <ci>gslmod</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>gappT</ci>
            <apply>
               <times/>
               <ci>gapp</ci>
               <ci>gapslmd</ci>
               <ci>xbmodsp</ci>
               <apply>
                  <power/>
                  <ci>Qgapp</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>hfmd</ci>
            <apply>
               <exp/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>xXBprer</ci>
                     </apply>
                     <apply>
                        <abs/>
                        <ci>xXBprer</ci>
                     </apply>
                  </apply>
                  <ci>hfmdc</ci>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <ci>xXBprer</ci>
                        <ci>x_0</ci>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>hbmd</ci>
            <apply>
               <exp/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>xXBpostr</ci>
                        <ci>x_0</ci>
                     </apply>
                     <apply>
                        <abs/>
                        <apply>
                           <minus/>
                           <ci>xXBpostr</ci>
                           <ci>x_0</ci>
                        </apply>
                     </apply>
                  </apply>
                  <ci>hbmdc</ci>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>xXBpostr</ci>
                           <ci>x_0</ci>
                        </apply>
                        <ci>x_0</ci>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>hfT</ci>
            <apply>
               <times/>
               <ci>hf</ci>
               <ci>hfmd</ci>
               <ci>xbmodsp</ci>
               <apply>
                  <power/>
                  <ci>Qhf</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>hbT</ci>
            <apply>
               <times/>
               <ci>hb</ci>
               <ci>hbmd</ci>
               <ci>xbmodsp</ci>
               <apply>
                  <power/>
                  <ci>Qhb</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>gxbmd</ci>
            <piecewise>
               <piece>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>sigmap</ci>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>x_0</ci>
                                 <ci>xXBpostr</ci>
                              </apply>
                              <ci>x_0</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>xXBpostr</ci>
                     <ci>x_0</ci>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>sigman</ci>
                        <apply>
                           <power/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>xXBpostr</ci>
                                 <ci>x_0</ci>
                              </apply>
                              <ci>x_0</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>gxbT</ci>
            <apply>
               <times/>
               <ci>gxb</ci>
               <ci>gxbmd</ci>
               <ci>xbmodsp</ci>
               <apply>
                  <power/>
                  <ci>Qgxb</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>TmpC</ci>
                        <cn cellml:units="celsius">37</cn>
                     </apply>
                     <cn cellml:units="celsius">10</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="regulation_and_crossbridge_cycling_state_equations">
      <variable units="dimensionless" public_interface="out" name="XBpostr" initial_value="0.01"/>
      <variable units="dimensionless" public_interface="out" name="XBprer" initial_value="0.01"/>
      <variable units="first_order_rate_constant" public_interface="out" name="dXBpostr"/>
      <variable units="first_order_rate_constant" public_interface="out" name="dXBprer"/>
      <variable units="dimensionless" name="N_NoXB" initial_value="0.99"/>
      <variable units="dimensionless" name="P_NoXB" initial_value="0.01"/>
      <variable units="dimensionless" name="P"/>
      <variable units="dimensionless" name="N" initial_value="0.97"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="first_order_rate_constant" public_interface="in" name="fappT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gappT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="hfT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="hbT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gxbT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="kn_pT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="kp_nT"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>N_NoXB</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kp_nT</ci>
                  <ci>P_NoXB</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>kn_pT</ci>
                  <ci>N_NoXB</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_NoXB</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kn_pT</ci>
                  <ci>N_NoXB</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>kp_nT</ci>
                  <ci>P_NoXB</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>N</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>kp_nT</ci>
                  <ci>P</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>kn_pT</ci>
                  <ci>N</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>XBprer</ci>
            </apply>
            <ci>dXBprer</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dXBprer</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>fappT</ci>
                     <ci>P</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hbT</ci>
                     <ci>XBpostr</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>gappT</ci>
                     <ci>XBprer</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hfT</ci>
                     <ci>XBprer</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>XBpostr</ci>
            </apply>
            <ci>dXBpostr</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dXBpostr</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>hfT</ci>
                  <ci>XBprer</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>hbT</ci>
                     <ci>XBpostr</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxbT</ci>
                     <ci>XBpostr</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>P</ci>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>N</ci>
                  </apply>
                  <ci>XBprer</ci>
               </apply>
               <ci>XBpostr</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="mean_strain_of_strongly_bound_states">
      <variable units="micrometre" public_interface="out" name="xXBpostr" initial_value="0"/>
      <variable units="micrometre" public_interface="out" name="xXBprer" initial_value="0.007"/>
      <variable units="micrometre_per_millisecond" public_interface="out" name="dxXBpostr"/>
      <variable units="micrometre_per_millisecond" public_interface="out" name="dxXBprer"/>
      <variable units="dimensionless" name="xPsi" initial_value="2"/>
      <variable units="dimensionless" name="dutyprer"/>
      <variable units="dimensionless" name="dutypostr"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="micrometre_per_millisecond" public_interface="in" name="dSL"/>
      <variable units="micrometre" public_interface="in" name="x_0"/>
      <variable units="first_order_rate_constant" public_interface="in" name="fappT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gappT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="hfT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="hbT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gxbT"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xXBprer</ci>
            </apply>
            <ci>dxXBprer</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dxXBprer</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <ci>dSL</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>xPsi</ci>
                     <ci>dutyprer</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>fappT</ci>
                        <apply>
                           <minus/>
                           <ci>xXBprer</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>hbT</ci>
                        <apply>
                           <minus/>
                           <ci>xXBpostr</ci>
                           <apply>
                              <plus/>
                              <ci>x_0</ci>
                              <ci>xXBprer</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>xXBpostr</ci>
            </apply>
            <ci>dxXBpostr</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dxXBpostr</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <ci>dSL</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>xPsi</ci>
                     <ci>dutypostr</ci>
                  </apply>
                  <ci>hfT</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>xXBprer</ci>
                        <ci>x_0</ci>
                     </apply>
                     <ci>xXBpostr</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dutyprer</ci>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>hbT</ci>
                     <ci>fappT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxbT</ci>
                     <ci>fappT</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>fappT</ci>
                     <ci>hfT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxbT</ci>
                     <ci>hfT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxbT</ci>
                     <ci>gappT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hbT</ci>
                     <ci>fappT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hbT</ci>
                     <ci>gappT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxbT</ci>
                     <ci>fappT</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dutypostr</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>fappT</ci>
                  <ci>hfT</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>fappT</ci>
                     <ci>hfT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxbT</ci>
                     <ci>hfT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxbT</ci>
                     <ci>gappT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hbT</ci>
                     <ci>fappT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hbT</ci>
                     <ci>gappT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxbT</ci>
                     <ci>fappT</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="normalised_active_and_passive_force">
      <variable units="dimensionless" public_interface="out" name="SSXBpostr"/>
      <variable units="dimensionless" public_interface="out" name="SSXBprer"/>
      <variable units="micrometre_per_millisecond" public_interface="out" name="dSL"/>
      <variable units="micrometre" public_interface="out" name="SL" initial_value="2"/>
      <variable units="micrometre" public_interface="out" name="SLmin" initial_value="1.4"/>
      <variable units="micrometre" public_interface="out" name="SLmax" initial_value="2.4"/>
      <variable units="millinewton_per_millimetre2" name="kxb" initial_value="120"/>
      <variable units="millinewton_micrometre_per_millimetre2" name="Fnordv"/>
      <variable units="millinewton_micrometre_per_millimetre2" name="force"/>
      <variable units="unit_normalised_force" name="active"/>
      <variable units="unit_normalised_force" name="ppforce"/>
      <variable units="unit_normalised_force" name="ppforce_t"/>
      <variable units="unit_normalised_force" name="ppforce_c"/>
      <variable units="unit_normalised_force" name="preload"/>
      <variable units="unit_normalised_force" name="afterload"/>
      <variable units="unit_normalised_force_millisecond" name="intf" initial_value="0"/>
      <variable units="micrometre" name="SL_c" initial_value="2.25"/>
      <variable units="micrometre" name="SLrest" initial_value="1.85"/>
      <variable units="micrometre" name="SLset" initial_value="1.9"/>
      <variable units="unit_normalised_force" name="constant_afterload" initial_value="0.001"/>
      <variable units="unit_normalised_force" name="PCon_t" initial_value="0.002"/>
      <variable units="per_micrometre" name="PExp_t" initial_value="10"/>
      <variable units="unit_normalised_force" name="PCon_c" initial_value="0.02"/>
      <variable units="per_micrometre" name="PExp_c" initial_value="70"/>
      <variable units="unit_normalised_force_millisecond2_per_micrometre" name="massf" initial_value="0.00005e6"/>
      <variable units="unit_normalised_force_millisecond_per_micrometre" name="visc" initial_value="0.003e3"/>
      <variable units="unit_normalised_force_per_micrometre" name="KSE" initial_value="1"/>
      <variable units="micrometre" public_interface="in" name="x_0"/>
      <variable units="first_order_rate_constant" public_interface="in" name="fapp"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gapp"/>
      <variable units="first_order_rate_constant" public_interface="in" name="hf"/>
      <variable units="first_order_rate_constant" public_interface="in" name="hb"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gxb"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="first_order_rate_constant" public_interface="in" name="fappT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gappT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="hfT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="hbT"/>
      <variable units="first_order_rate_constant" public_interface="in" name="gxbT"/>
      <variable units="dimensionless" public_interface="in" name="XBpostr"/>
      <variable units="dimensionless" public_interface="in" name="XBprer"/>
      <variable units="micrometre" public_interface="in" name="xXBpostr"/>
      <variable units="micrometre" public_interface="in" name="xXBprer"/>
      <variable units="dimensionless" public_interface="in" name="SOVFThick"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>SSXBprer</ci>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>hb</ci>
                     <ci>fapp</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxb</ci>
                     <ci>fapp</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>fapp</ci>
                     <ci>hf</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxb</ci>
                     <ci>hf</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxb</ci>
                     <ci>gapp</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hb</ci>
                     <ci>fapp</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hb</ci>
                     <ci>gapp</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxb</ci>
                     <ci>fapp</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>SSXBpostr</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>fapp</ci>
                  <ci>hf</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>fapp</ci>
                     <ci>hf</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxb</ci>
                     <ci>hf</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxb</ci>
                     <ci>gapp</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hb</ci>
                     <ci>fapp</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>hb</ci>
                     <ci>gapp</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>gxb</ci>
                     <ci>fapp</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Fnordv</ci>
            <apply>
               <times/>
               <ci>kxb</ci>
               <ci>x_0</ci>
               <ci>SSXBpostr</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>force</ci>
            <apply>
               <times/>
               <ci>kxb</ci>
               <ci>SOVFThick</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>xXBpostr</ci>
                     <ci>XBpostr</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>xXBprer</ci>
                     <ci>XBprer</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>active</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="unit_normalised_force">1</cn>
                  <ci>force</ci>
               </apply>
               <ci>Fnordv</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ppforce_t</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>SL</ci>
                     <ci>SLrest</ci>
                  </apply>
                  <apply>
                     <abs/>
                     <apply>
                        <minus/>
                        <ci>SL</ci>
                        <ci>SLrest</ci>
                     </apply>
                  </apply>
               </apply>
               <ci>PCon_t</ci>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>PExp_t</ci>
                        <apply>
                           <abs/>
                           <apply>
                              <minus/>
                              <ci>SL</ci>
                              <ci>SLrest</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ppforce_c</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <ci>PCon_c</ci>
                     <apply>
                        <minus/>
                        <apply>
                           <exp/>
                           <apply>
                              <times/>
                              <ci>PExp_c</ci>
                              <apply>
                                 <abs/>
                                 <apply>
                                    <minus/>
                                    <ci>SL</ci>
                                    <ci>SL_c</ci>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>SL</ci>
                     <ci>SL_c</ci>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="unit_normalised_force">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>ppforce</ci>
            <apply>
               <plus/>
               <ci>ppforce_t</ci>
               <ci>ppforce_c</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>preload</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <abs/>
                     <apply>
                        <minus/>
                        <ci>SLset</ci>
                        <ci>SLrest</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>SLset</ci>
                     <ci>SLrest</ci>
                  </apply>
               </apply>
               <ci>PCon_t</ci>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <ci>PExp_t</ci>
                        <apply>
                           <abs/>
                           <apply>
                              <minus/>
                              <ci>SLset</ci>
                              <ci>SLrest</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>afterload</ci>
            <piecewise>
               <piece>
                  <ci>constant_afterload</ci>
                  <apply>
                     <gt/>
                     <ci>SL</ci>
                     <ci>SL_c</ci>
                  </apply>
               </piece>
               <piece>
                  <apply>
                     <times/>
                     <ci>KSE</ci>
                     <apply>
                        <minus/>
                        <ci>SLset</ci>
                        <ci>SL</ci>
                     </apply>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>SL</ci>
                     <ci>SL_c</ci>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="unit_normalised_force">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>SL</ci>
            </apply>
            <ci>dSL</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dSL</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <plus/>
                        <ci>intf</ci>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>SLset</ci>
                              <ci>SL</ci>
                           </apply>
                           <ci>visc</ci>
                        </apply>
                     </apply>
                     <ci>massf</ci>
                  </apply>
                  <apply>
                     <and/>
                     <apply>
                        <leq/>
                        <ci>SL</ci>
                        <ci>SLmax</ci>
                     </apply>
                     <apply>
                        <gt/>
                        <ci>SL</ci>
                        <ci>SLmin</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="micrometre_per_millisecond">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>intf</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>preload</ci>
                  <ci>afterload</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>ppforce</ci>
                  <ci>active</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding">
      <variable units="dimensionless" name="FrSBXB"/>
      <variable units="first_order_rate_constant" name="dFrSBXB"/>
      <variable units="micrometre_per_millisecond" name="dsovr_ze"/>
      <variable units="micrometre_per_millisecond" name="dsovr_cle"/>
      <variable units="micrometre_per_millisecond" name="dlen_sovr"/>
      <variable units="first_order_rate_constant" name="dSOVFThick"/>
      <variable units="first_order_rate_constant" name="dSOVFThin"/>
      <variable units="millinewton_per_millimetre2" name="kxb" initial_value="120"/>
      <variable units="millinewton_micrometre_per_millimetre2_per_millisecond" name="dforce"/>
      <variable units="micromolar" name="Trop_conc" initial_value="70"/>
      <variable units="micromolar" name="TropTot"/>
      <variable units="micromolar_per_millisecond" name="dTropTot"/>
      <variable units="micrometre_per_millisecond" public_interface="in" name="dSL"/>
      <variable units="micrometre" public_interface="in" name="SL"/>
      <variable units="first_order_rate_constant" public_interface="in" name="dXBpostr"/>
      <variable units="first_order_rate_constant" public_interface="in" name="dXBprer"/>
      <variable units="micrometre_per_millisecond" public_interface="in" name="dxXBpostr"/>
      <variable units="micrometre_per_millisecond" public_interface="in" name="dxXBprer"/>
      <variable units="dimensionless" public_interface="in" name="XBpostr"/>
      <variable units="dimensionless" public_interface="in" name="XBprer"/>
      <variable units="micrometre" public_interface="in" name="xXBpostr"/>
      <variable units="micrometre" public_interface="in" name="xXBprer"/>
      <variable units="dimensionless" public_interface="in" name="SSXBpostr"/>
      <variable units="dimensionless" public_interface="in" name="SSXBprer"/>
      <variable units="dimensionless" public_interface="in" name="SOVFThick"/>
      <variable units="dimensionless" public_interface="in" name="SOVFThin"/>
      <variable units="micrometre" public_interface="in" name="len_thin"/>
      <variable units="micrometre" public_interface="in" name="len_thick"/>
      <variable units="micrometre" public_interface="in" name="len_hbare"/>
      <variable units="dimensionless" public_interface="in" name="TRPNCaL"/>
      <variable units="dimensionless" public_interface="in" name="TRPNCaH"/>
      <variable units="first_order_rate_constant" public_interface="in" name="dTRPNCaL"/>
      <variable units="first_order_rate_constant" public_interface="in" name="dTRPNCaH"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>FrSBXB</ci>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <ci>XBpostr</ci>
                  <ci>XBprer</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>SSXBpostr</ci>
                  <ci>SSXBprer</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dFrSBXB</ci>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <ci>dXBpostr</ci>
                  <ci>dXBprer</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>SSXBpostr</ci>
                  <ci>SSXBprer</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dsovr_ze</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">0.5</cn>
                     </apply>
                     <ci>dSL</ci>
                  </apply>
                  <apply>
                     <lt/>
                     <ci>SL</ci>
                     <ci>len_thick</ci>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="micrometre_per_millisecond">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>dsovr_cle</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">0.5</cn>
                     </apply>
                     <ci>dSL</ci>
                  </apply>
                  <apply>
                     <gt/>
                     <apply>
                        <minus/>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>len_thin</ci>
                        </apply>
                        <ci>SL</ci>
                     </apply>
                     <ci>len_hbare</ci>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="micrometre_per_millisecond">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>dlen_sovr</ci>
            <apply>
               <minus/>
               <ci>dsovr_ze</ci>
               <ci>dsovr_cle</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dSOVFThin</ci>
            <apply>
               <divide/>
               <ci>dlen_sovr</ci>
               <ci>len_thin</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dSOVFThick</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">2</cn>
                  <ci>dlen_sovr</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>len_thick</ci>
                  <ci>len_hbare</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>TropTot</ci>
            <apply>
               <times/>
               <ci>Trop_conc</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>SOVFThin</ci>
                     </apply>
                     <ci>TRPNCaL</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>SOVFThin</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>FrSBXB</ci>
                           <ci>TRPNCaH</ci>
                        </apply>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                              <ci>FrSBXB</ci>
                           </apply>
                           <ci>TRPNCaL</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dTropTot</ci>
            <apply>
               <times/>
               <ci>Trop_conc</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <ci>dSOVFThin</ci>
                     </apply>
                     <ci>TRPNCaL</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>SOVFThin</ci>
                     </apply>
                     <ci>dTRPNCaL</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>dSOVFThin</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>FrSBXB</ci>
                           <ci>TRPNCaH</ci>
                        </apply>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <cn cellml:units="dimensionless">1</cn>
                              <ci>FrSBXB</ci>
                           </apply>
                           <ci>TRPNCaL</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>SOVFThin</ci>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>dFrSBXB</ci>
                              <ci>TRPNCaH</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>FrSBXB</ci>
                              <ci>dTRPNCaH</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <minus/>
                                 <cn cellml:units="dimensionless">1</cn>
                                 <ci>FrSBXB</ci>
                              </apply>
                              <ci>dTRPNCaL</ci>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>dFrSBXB</ci>
                           <ci>TRPNCaL</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dforce</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>kxb</ci>
                  <ci>dSOVFThick</ci>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>xXBpostr</ci>
                        <ci>XBpostr</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>xXBprer</ci>
                        <ci>XBprer</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>kxb</ci>
                  <ci>SOVFThick</ci>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>dxXBpostr</ci>
                        <ci>XBpostr</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>xXBpostr</ci>
                        <ci>dXBpostr</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>dxXBprer</ci>
                        <ci>XBprer</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>xXBprer</ci>
                        <ci>dXBprer</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="equation_for_simulated_calcium_transient">
      <variable units="micromolar" public_interface="out" name="Cai"/>
      <variable units="dimensionless" name="beta"/>
      <variable units="millisecond" name="tau1" initial_value="20"/>
      <variable units="millisecond" name="tau2" initial_value="110"/>
      <variable units="millisecond" name="start_time" initial_value="5"/>
      <variable units="micromolar" name="Ca_amplitude" initial_value="1.45"/>
      <variable units="micromolar" name="Ca_diastolic" initial_value="0.09"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>beta</ci>
            <apply>
               <minus/>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>tau1</ci>
                     <ci>tau2</ci>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <divide/>
                           <ci>tau1</ci>
                           <ci>tau2</ci>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <power/>
                  <apply>
                     <divide/>
                     <ci>tau1</ci>
                     <ci>tau2</ci>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <divide/>
                           <ci>tau2</ci>
                           <ci>tau1</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Cai</ci>
            <piecewise>
               <piece>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>Ca_amplitude</ci>
                              <ci>Ca_diastolic</ci>
                           </apply>
                           <ci>beta</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <minus/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>start_time</ci>
                                    </apply>
                                 </apply>
                                 <ci>tau1</ci>
                              </apply>
                           </apply>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <minus/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>start_time</ci>
                                    </apply>
                                 </apply>
                                 <ci>tau2</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <ci>Ca_diastolic</ci>
                  </apply>
                  <apply>
                     <gt/>
                     <ci>time</ci>
                     <ci>start_time</ci>
                  </apply>
               </piece>
               <otherwise>
                  <ci>Ca_diastolic</ci>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="model_parameters">
      <variable units="micrometre" public_interface="out" name="len_thin" initial_value="1.2"/>
      <variable units="micrometre" public_interface="out" name="len_thick" initial_value="1.65"/>
      <variable units="micrometre" public_interface="out" name="len_hbare" initial_value="0.1"/>
      <variable units="celsius" public_interface="out" name="TmpC" initial_value="24"/>
      <variable units="micrometre" public_interface="out" name="x_0" initial_value="0.007"/>
   </component>
   <connection>
      <map_components component_2="environment" component_1="Ca_binding_to_troponin_to_thin_filament_regulation"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="regulation_and_crossbridge_cycling_state_equations"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="mean_strain_of_strongly_bound_states"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="normalised_active_and_passive_force"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="equation_for_simulated_calcium_transient"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="normalised_active_and_passive_force" component_1="sarcomere_geometry"/>
      <map_variables variable_2="SLmax" variable_1="SLmax"/>
      <map_variables variable_2="SLmin" variable_1="SLmin"/>
      <map_variables variable_2="SL" variable_1="SL"/>
      <map_variables variable_2="SOVFThick" variable_1="SOVFThick"/>
   </connection>
   <connection>
      <map_components component_2="model_parameters" component_1="sarcomere_geometry"/>
      <map_variables variable_2="len_thin" variable_1="len_thin"/>
      <map_variables variable_2="len_thick" variable_1="len_thick"/>
      <map_variables variable_2="len_hbare" variable_1="len_hbare"/>
   </connection>
   <connection>
      <map_components component_2="model_parameters" component_1="Ca_binding_to_troponin_to_thin_filament_regulation"/>
      <map_variables variable_2="TmpC" variable_1="TmpC"/>
   </connection>
   <connection>
      <map_components component_2="equation_for_simulated_calcium_transient" component_1="Ca_binding_to_troponin_to_thin_filament_regulation"/>
      <map_variables variable_2="Cai" variable_1="Cai"/>
   </connection>
   <connection>
      <map_components component_2="sarcomere_geometry" component_1="Ca_binding_to_troponin_to_thin_filament_regulation"/>
      <map_variables variable_2="SOVFThin" variable_1="SOVFThin"/>
   </connection>
   <connection>
      <map_components component_2="model_parameters" component_1="thin_filament_regulation_and_crossbridge_cycling_rates"/>
      <map_variables variable_2="x_0" variable_1="x_0"/>
      <map_variables variable_2="TmpC" variable_1="TmpC"/>
   </connection>
   <connection>
      <map_components component_2="sarcomere_geometry" component_1="thin_filament_regulation_and_crossbridge_cycling_rates"/>
      <map_variables variable_2="SOVFThick" variable_1="SOVFThick"/>
   </connection>
   <connection>
      <map_components component_2="mean_strain_of_strongly_bound_states" component_1="thin_filament_regulation_and_crossbridge_cycling_rates"/>
      <map_variables variable_2="xXBpostr" variable_1="xXBpostr"/>
      <map_variables variable_2="xXBprer" variable_1="xXBprer"/>
      <map_variables variable_2="fappT" variable_1="fappT"/>
      <map_variables variable_2="gappT" variable_1="gappT"/>
      <map_variables variable_2="hfT" variable_1="hfT"/>
      <map_variables variable_2="hbT" variable_1="hbT"/>
      <map_variables variable_2="gxbT" variable_1="gxbT"/>
   </connection>
   <connection>
      <map_components component_2="thin_filament_regulation_and_crossbridge_cycling_rates" component_1="regulation_and_crossbridge_cycling_state_equations"/>
      <map_variables variable_2="fappT" variable_1="fappT"/>
      <map_variables variable_2="gappT" variable_1="gappT"/>
      <map_variables variable_2="hfT" variable_1="hfT"/>
      <map_variables variable_2="hbT" variable_1="hbT"/>
      <map_variables variable_2="gxbT" variable_1="gxbT"/>
   </connection>
   <connection>
      <map_components component_2="Ca_binding_to_troponin_to_thin_filament_regulation" component_1="regulation_and_crossbridge_cycling_state_equations"/>
      <map_variables variable_2="kn_pT" variable_1="kn_pT"/>
      <map_variables variable_2="kp_nT" variable_1="kp_nT"/>
   </connection>
   <connection>
      <map_components component_2="model_parameters" component_1="mean_strain_of_strongly_bound_states"/>
      <map_variables variable_2="x_0" variable_1="x_0"/>
   </connection>
   <connection>
      <map_components component_2="thin_filament_regulation_and_crossbridge_cycling_rates" component_1="normalised_active_and_passive_force"/>
      <map_variables variable_2="fapp" variable_1="fapp"/>
      <map_variables variable_2="gapp" variable_1="gapp"/>
      <map_variables variable_2="hf" variable_1="hf"/>
      <map_variables variable_2="hb" variable_1="hb"/>
      <map_variables variable_2="gxb" variable_1="gxb"/>
      <map_variables variable_2="fappT" variable_1="fappT"/>
      <map_variables variable_2="gappT" variable_1="gappT"/>
      <map_variables variable_2="hfT" variable_1="hfT"/>
      <map_variables variable_2="hbT" variable_1="hbT"/>
      <map_variables variable_2="gxbT" variable_1="gxbT"/>
   </connection>
   <connection>
      <map_components component_2="regulation_and_crossbridge_cycling_state_equations" component_1="normalised_active_and_passive_force"/>
      <map_variables variable_2="XBpostr" variable_1="XBpostr"/>
      <map_variables variable_2="XBprer" variable_1="XBprer"/>
   </connection>
   <connection>
      <map_components component_2="mean_strain_of_strongly_bound_states" component_1="normalised_active_and_passive_force"/>
      <map_variables variable_2="xXBpostr" variable_1="xXBpostr"/>
      <map_variables variable_2="xXBprer" variable_1="xXBprer"/>
      <map_variables variable_2="dSL" variable_1="dSL"/>
   </connection>
   <connection>
      <map_components component_2="model_parameters" component_1="normalised_active_and_passive_force"/>
      <map_variables variable_2="x_0" variable_1="x_0"/>
   </connection>
   <connection>
      <map_components component_2="model_parameters" component_1="calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding"/>
      <map_variables variable_2="len_thin" variable_1="len_thin"/>
      <map_variables variable_2="len_thick" variable_1="len_thick"/>
      <map_variables variable_2="len_hbare" variable_1="len_hbare"/>
   </connection>
   <connection>
      <map_components component_2="Ca_binding_to_troponin_to_thin_filament_regulation" component_1="calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding"/>
      <map_variables variable_2="TRPNCaL" variable_1="TRPNCaL"/>
      <map_variables variable_2="TRPNCaH" variable_1="TRPNCaH"/>
      <map_variables variable_2="dTRPNCaL" variable_1="dTRPNCaL"/>
      <map_variables variable_2="dTRPNCaH" variable_1="dTRPNCaH"/>
   </connection>
   <connection>
      <map_components component_2="regulation_and_crossbridge_cycling_state_equations" component_1="calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding"/>
      <map_variables variable_2="XBpostr" variable_1="XBpostr"/>
      <map_variables variable_2="XBprer" variable_1="XBprer"/>
      <map_variables variable_2="dXBpostr" variable_1="dXBpostr"/>
      <map_variables variable_2="dXBprer" variable_1="dXBprer"/>
   </connection>
   <connection>
      <map_components component_2="mean_strain_of_strongly_bound_states" component_1="calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding"/>
      <map_variables variable_2="xXBpostr" variable_1="xXBpostr"/>
      <map_variables variable_2="xXBprer" variable_1="xXBprer"/>
      <map_variables variable_2="dxXBpostr" variable_1="dxXBpostr"/>
      <map_variables variable_2="dxXBprer" variable_1="dxXBprer"/>
   </connection>
   <connection>
      <map_components component_2="normalised_active_and_passive_force" component_1="calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding"/>
      <map_variables variable_2="SSXBpostr" variable_1="SSXBpostr"/>
      <map_variables variable_2="SSXBprer" variable_1="SSXBprer"/>
      <map_variables variable_2="SL" variable_1="SL"/>
      <map_variables variable_2="dSL" variable_1="dSL"/>
   </connection>
   <connection>
      <map_components component_2="sarcomere_geometry" component_1="calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding"/>
      <map_variables variable_2="SOVFThin" variable_1="SOVFThin"/>
      <map_variables variable_2="SOVFThick" variable_1="SOVFThick"/>
   </connection>
<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#566d983b-56b3-42f6-a590-f7d68d809dc0">
    <rdf:li>cardiac</rdf:li>
    <rdf:li>crossbridge</rdf:li>
    <rdf:li>heart</rdf:li>
    <rdf:li>muscle</rdf:li>
    <rdf:li>myofilament</rdf:li>
    <rdf:li>actin</rdf:li>
    <rdf:li>myosin</rdf:li>
    <rdf:li>calcium</rdf:li>
    <rdf:li>troponin</rdf:li>
    <rdf:li>tropomyosin</rdf:li>
    <rdf:li>cooperativity</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#36951891-93b4-40cf-9a25-160b1c6946a2">
    <rdf:li rdf:resource="rdf:#f3a04525-2b12-4772-9e23-cd93c67e9e78"/>
    <rdf:li rdf:resource="rdf:#051cc6a3-f44a-4518-b07b-2ab5d89e597a"/>
    <rdf:li rdf:resource="rdf:#0de0bd2d-04d9-4890-a28b-cadd70e5aa32"/>
    <rdf:li rdf:resource="rdf:#c2cc0b51-3e8d-4b12-8b69-962ab631d207"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#62b099d5-de9a-4737-9c69-835d3e9ca3ca">
    <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="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#7b266c4c-de53-46d8-941a-2774821064dd"/>
    <dcterms:created rdf:resource="rdf:#e0e183ed-ea1e-4d7c-973d-78cf014528df"/>
    <dc:creator rdf:resource="rdf:#3a56bba3-ed4f-493d-8f97-a40ca56450cb"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d0bebae0-d725-466b-87fe-94e7e6d0a0e9">
    <bqs:Pubmed_id>18234826</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#2bff6762-8c0a-4aff-ae92-c6364a9b8e7c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#af0b201a-d2c8-48ed-8e0e-f399e5fb12c7">
    <vCard:Given>Pieter</vCard:Given>
    <vCard:Family>de Tombe</vCard:Family>
    <vCard:Other>P</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9877a75d-4984-41e7-9944-e52d8bcc8eb5">
    <dc:creator rdf:resource="rdf:#70a62a5c-d04e-4d49-ab97-3412f2fbc7a4"/>
    <rdf:value>This is the CellML description of Rice et al.'s 2008 model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#af9b6f52-32c9-45bf-a9ad-8bf0941e40d7">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3a56bba3-ed4f-493d-8f97-a40ca56450cb">
    <vCard:ORG rdf:resource="rdf:#7165023a-3752-4722-89a5-7df81bb54bed"/>
    <vCard:EMAIL rdf:resource="rdf:#62b099d5-de9a-4737-9c69-835d3e9ca3ca"/>
    <vCard:N rdf:resource="rdf:#0178c64c-9cee-4eea-b339-154504b893b8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0178c64c-9cee-4eea-b339-154504b893b8">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="#rice_model_2008">
    <dc:title>
        Rice et al.'s 2008 model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations.
      </dc:title>
    <cmeta:comment rdf:resource="rdf:#9877a75d-4984-41e7-9944-e52d8bcc8eb5"/>
    <bqs:reference rdf:resource="rdf:#6a666545-2193-4623-b2aa-d35ce0c9f4d5"/>
    <bqs:reference rdf:resource="rdf:#d0bebae0-d725-466b-87fe-94e7e6d0a0e9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0de0bd2d-04d9-4890-a28b-cadd70e5aa32">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#45dcdcbc-faab-4d0e-a666-c3e209e2177d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7165023a-3752-4722-89a5-7df81bb54bed">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e0e183ed-ea1e-4d7c-973d-78cf014528df">
    <dcterms:W3CDTF>2008-07-08T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c217df42-16db-4e0d-a592-7863faeb51f6">
    <dc:title>Biophysical Journal</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0d84e6f0-0e6a-4c1d-bbbe-2359bdc5370b">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Rice</vCard:Family>
    <vCard:Other>Jeremy</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f3a04525-2b12-4772-9e23-cd93c67e9e78">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#0d84e6f0-0e6a-4c1d-bbbe-2359bdc5370b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6a666545-2193-4623-b2aa-d35ce0c9f4d5">
    <dc:subject rdf:resource="rdf:#a8a75408-610d-48a5-a438-9e813e7069e1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7b266c4c-de53-46d8-941a-2774821064dd">
    <dc:creator rdf:resource="rdf:#af9b6f52-32c9-45bf-a9ad-8bf0941e40d7"/>
    <rdf:value>This CellML model is known to run in both PCEnv and COR.  The units have been checked and they are consistent.  The CellML model was created using a combination of both the published paper and also the author's Matlab code.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#45dcdcbc-faab-4d0e-a666-c3e209e2177d">
    <vCard:Given>Donald</vCard:Given>
    <vCard:Family>Bers</vCard:Family>
    <vCard:Other>M</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2bff6762-8c0a-4aff-ae92-c6364a9b8e7c">
    <dc:creator rdf:resource="rdf:#36951891-93b4-40cf-9a25-160b1c6946a2"/>
    <dc:title>Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations</dc:title>
    <bqs:volume/>
    <bqs:first_page/>
    <bqs:Journal rdf:resource="rdf:#c217df42-16db-4e0d-a592-7863faeb51f6"/>
    <dcterms:issued rdf:resource="rdf:#966bee92-5d08-4009-8f7f-2e7b6f1a9d2d"/>
    <bqs:last_page/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a8a75408-610d-48a5-a438-9e813e7069e1">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#566d983b-56b3-42f6-a590-f7d68d809dc0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#051cc6a3-f44a-4518-b07b-2ab5d89e597a">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#3185c990-f72c-4fb1-a94e-aeb89023a28c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3185c990-f72c-4fb1-a94e-aeb89023a28c">
    <vCard:Given>Fei</vCard:Given>
    <vCard:Family>Wang</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#966bee92-5d08-4009-8f7f-2e7b6f1a9d2d">
    <dcterms:W3CDTF>2008</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#70a62a5c-d04e-4d49-ab97-3412f2fbc7a4">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c2cc0b51-3e8d-4b12-8b69-962ab631d207">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#af0b201a-d2c8-48ed-8e0e-f399e5fb12c7"/>
  </rdf:Description>
</rdf:RDF>
</model>