Location: Salem, Saidel, Stanley, Cabrera, 2002 @ e0281b3b5e99 / salem_saidel_stanley_cabrera_2002.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2010-08-05 01:20:43+12:00
Desc:
Updated documentation and fixed broken annotations.
Permanent Source URI:
http://models.cellml.org/workspace/salem_saidel_stanley_cabrera_2002/rawfile/e0281b3b5e99b20e072081571349dcd2af085f6e/salem_saidel_stanley_cabrera_2002.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="salem_model_2002" name="salem_model_2002">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Mechanistic model of myocardial energy metabolism under normal and ischemic conditions</title>
  <author>
    <firstname>Autumn</firstname>
    <othername>A.</othername>
    <surname>Cuellar</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
          This CellML model runs in both COR and OpenCell and we believe it is close to recreating the original published results.  The units have been checked and they are consistent.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: A moderate reduction in coronary blood flow results in decreased myocardial oxygen consumption, accelerated glycolysis, decreased pyruvate oxidation, and lactate accumulation. To quantitatively understand cardiac metabolism during ischemia, we have developed a mechanistic, mathematical model based on biochemical mass balances and reaction kinetics in cardiac cells. By numerical solution of model equations, computer simulations showed the dynamic responses in glucose, fatty acid, glucose-6-phosphate, glycogen, triglyceride, pyruvate, lactate, acetyl-CoA, and free-CoA as well as CO2, O2, phosphocreatine/creatine, nicotinamide adenine dinucleotide (reduced form)/nicotinamide adenine dinucleotide (oxidized form) (NADH/NAD+), and adenosine diphosphate/adenosine triphosphate (ADP/ATP). When myocardial ischemia was simulated by a 60% reduction in coronary blood flow, the model generated myocardial concentrations, uptakes, and fluxes that were consistent with experimental data from in vivo pig studies. After 60 min of ischemia the concentrations of glycogen, phosphocreatine, and ATP were decreased by 60%, 75%, and 50%, respectively. With the onset of ischemia, myocardial lactate concentration increased and the myocardium switched from net consumer to net producer of lactate. Our model predicted a rapid 13-fold increase in NADH/NAD+, but only a twofold increase in the ratio of acetyl-CoA to free-CoA. These findings are consistent with the concept that pyruvate oxidation is inhibited during ischemia partially by the rise in NADH/NAD+.
</para>

<para>
The original paper reference is cited below:
</para>

<para>
Mechanistic model of myocardial energy metabolism under normal and ischemic conditions, Jennifer E. Salem, Gerald M. Saidel, William C. Stanley, and Marco E. Cabrera, 2002, <emphasis>Annals of Biomedical Engineering</emphasis>, 30, 202-216. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=11962772&amp;dopt=Abstract">PubMed ID: 11962772</ulink> 
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>Schematic diagrams of the RyR models</title>
    </objectinfo>
    <imagedata fileref="salem_2002.png"/>
  </imageobject>
</mediaobject>
<caption>The reaction pathway of the Salem <emphasis>et al</emphasis>. model.  Glucose (through the GLUT transporters), lactate, CO<subscript>2</subscript>, O<subscript>2</subscript>, and palmitate (fatty acids) are the only species in this model considered to transfer between the blood and myocytes.</caption>
</informalfigure>

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

  <units name="umol_per_ml">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" prefix="milli" exponent="-1"/>
  </units>
  
  <units name="g_per_ml">
    <unit units="gram"/>
    <unit units="litre" prefix="milli" exponent="-1"/>
  </units>

  <units name="ml_per_g_min">
    <unit units="litre" prefix="milli"/>
    <unit units="gram" exponent="-1"/>
    <unit units="minute" exponent="-1"/>
  </units>
  
  <units name="minute">
    <unit units="second" multiplier="60"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="minute" exponent="-1"/>
  </units>
  
  <units name="flux">
    <unit units="umol_per_g"/>
    <unit units="minute" exponent="-1"/>
  </units>

  
  <component name="environment">
    <variable units="minute" public_interface="out" name="time"/>
  </component>

  <!-- reactions -->

  <component name="v1">
    <variable units="flux" public_interface="out" name="v1"/>
    
    <variable units="flux" name="Vmax" initial_value="5.90"/> 
    <variable units="umol_per_g" name="Km" initial_value="0.05"/> 
     
    <variable units="umol_per_g" public_interface="in" name="GL"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="1"> 
      <apply>
        <eq/> 
        <ci> v1 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <divide/>
            <ci> GL </ci>
            <apply>
              <plus/>
              <ci> Km </ci>
              <ci> GL </ci>
            </apply>
          </apply>
        </apply> 
      </apply> 
    </math> 
  </component>

  <component name="v2">
    <variable units="flux" public_interface="out" name="v2"/>
    
    <variable units="flux" name="Vmax" initial_value="67.6"/> 
    <variable units="umol_per_g" name="Km" initial_value="0.01"/> 
     
    <variable units="umol_per_g" public_interface="in" name="O2"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="2"> 
      <apply>
        <eq/> 
        <ci> v2 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <divide/>
            <ci> O2 </ci>
            <apply>
              <plus/>
              <ci> Km </ci>
              <ci> O2 </ci>
            </apply>
          </apply>
        </apply> 
      </apply> 
    </math> 
  </component>

  <component name="v3">
    <variable units="flux" public_interface="out" name="v3"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="4.90"/> 
    <variable units="first_order_rate_constant" name="k3"/>
    
    <variable units="umol_per_g" public_interface="in" name="FA"/>
    <variable units="dimensionless" public_interface="in" name="PS"/>
    <variable units="dimensionless" public_interface="in" name="PS0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="3">
      <apply>
        <eq/>
        <ci> v3 </ci>
        <apply>
          <times/>
          <ci> k3 </ci>
          <ci> FA </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="4"> 
      <apply>
        <eq/> 
        <ci> k3 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> PS </ci>
              <cn cellml:units="dimensionless"> -1 </cn>
            </apply>
            <apply>
              <plus/>
              <apply>
                <power/>
                <ci> PS0 </ci>
                <cn cellml:units="dimensionless"> -1 </cn>
              </apply>
              <apply>
                <power/>
                <ci> PS </ci>
                <cn cellml:units="dimensionless"> -1 </cn>
              </apply>
            </apply>
          </apply>
        </apply> 
      </apply> 
    </math> 
  </component>

  <component name="v4">
    <variable units="flux" public_interface="out" name="v4"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="21.3"/> 
    <variable units="first_order_rate_constant" name="k4"/>
    <variable units="dimensionless" name="epsilon" initial_value="0.6"/>
    
    <variable units="umol_per_g" public_interface="in" name="FA"/>
    <variable units="dimensionless" public_interface="in" name="RS"/>
    <variable units="dimensionless" public_interface="in" name="RS0"/>
    <variable units="dimensionless" public_interface="in" name="AF"/>
    <variable units="dimensionless" public_interface="in" name="AF0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="5">
      <apply>
        <eq/>
        <ci> v4 </ci>
        <apply>
          <times/>
          <ci> k4 </ci>
          <ci> FA </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="6"> 
      <apply>
        <eq/> 
        <ci> k4 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> AF </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> AF0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> AF </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1 </cn>
                <ci> epsilon </ci>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> RS </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> RS0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> RS </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component>  

  <component name="v5">
    <variable units="flux" public_interface="out" name="v5"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="2.82"/> 
    <variable units="first_order_rate_constant" name="k5"/>
    <variable units="dimensionless" name="epsilon" initial_value="0.254"/>
    
    <variable units="umol_per_g" public_interface="in" name="GP"/>
    <variable units="dimensionless" public_interface="in" name="RS"/>
    <variable units="dimensionless" public_interface="in" name="RS0"/>
    <variable units="dimensionless" public_interface="in" name="PS"/>
    <variable units="dimensionless" public_interface="in" name="PS0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="7">
      <apply>
        <eq/>
        <ci> v5 </ci>
        <apply>
          <times/>
          <ci> k5 </ci>
          <ci> GP </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="8"> 
      <apply>
        <eq/> 
        <ci> k5 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <apply>
                <divide/>
                <ci> PS </ci>
                <apply>
                  <plus/>
                  <ci> PS0 </ci>
                  <ci> PS </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1 </cn>
                <ci> epsilon </ci>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> RS </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> RS0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> RS </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component>  
  
  <component name="v6">
    <variable units="flux" public_interface="out" name="v6"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="3.14"/> 
    <variable units="first_order_rate_constant" name="k6"/>
    <variable units="dimensionless" name="epsilon" initial_value="0.5"/>
    
    <variable units="umol_per_g" public_interface="in" name="GP"/>
    <variable units="dimensionless" public_interface="in" name="CS"/>
    <variable units="dimensionless" public_interface="in" name="CS0"/>
    <variable units="dimensionless" public_interface="in" name="PS"/>
    <variable units="dimensionless" public_interface="in" name="PS0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="9">
      <apply>
        <eq/>
        <ci> v6 </ci>
        <apply>
          <times/>
          <ci> k6 </ci>
          <ci> GP </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="10"> 
      <apply>
        <eq/> 
        <ci> k6 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> PS </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> PS0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> PS </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1 </cn>
                <ci> epsilon </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <apply>
                    <power/>
                    <ci> CS </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <plus/>
                    <apply>
                      <power/>
                      <ci> CS0 </ci>
                      <cn cellml:units="dimensionless"> -1 </cn>
                    </apply>
                    <apply>
                      <power/>
                      <ci> CS </ci>
                      <cn cellml:units="dimensionless"> -1 </cn>
                    </apply>
                  </apply>
                </apply>
                <cn cellml:units="dimensionless"> 2 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component> 

  <component name="v7">
    <variable units="flux" public_interface="out" name="v7"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="0.0162"/> 
    <variable units="first_order_rate_constant" name="k7"/>
    <variable units="dimensionless" name="epsilon" initial_value="0.5"/>
    
    <variable units="umol_per_g" public_interface="in" name="GY"/>
    <variable units="dimensionless" public_interface="in" name="CS"/>
    <variable units="dimensionless" public_interface="in" name="CS0"/>
    <variable units="dimensionless" public_interface="in" name="PS"/>
    <variable units="dimensionless" public_interface="in" name="PS0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="11">
      <apply>
        <eq/>
        <ci> v7 </ci>
        <apply>
          <times/>
          <ci> k7 </ci>
          <ci> GY </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="12"> 
      <apply>
        <eq/> 
        <ci> k7 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <apply>
                <divide/>
                <ci> PS </ci>
                <apply>
                  <plus/>
                  <ci> PS0 </ci>
                  <ci> PS </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1 </cn>
                <ci> epsilon </ci>
              </apply>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> CS </ci>
                  <apply>
                    <plus/>
                    <ci> CS0 </ci>
                    <ci> CS </ci>
                  </apply>
                </apply>
                <cn cellml:units="dimensionless"> 2 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component> 

  <component name="v8">
    <variable units="flux" public_interface="out" name="v8"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="0.005"/> 
    <variable units="first_order_rate_constant" name="k8"/>
    
    <variable units="umol_per_g" public_interface="in" name="TG"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="13">
      <apply>
        <eq/>
        <ci> v8 </ci>
        <apply>
          <times/>
          <ci> k8 </ci>
          <ci> TG </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="14"> 
      <apply>
        <eq/> 
        <ci> k8 </ci> 
        <ci> Vmax </ci>
      </apply>
    </math> 
  </component> 

  <component name="v9">
    <variable units="flux" public_interface="out" name="v9"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="1.8"/> 
    <variable units="first_order_rate_constant" name="k9"/>
    
    <variable units="umol_per_g" public_interface="in" name="PY"/>
    <variable units="dimensionless" public_interface="in" name="RS"/>
    <variable units="dimensionless" public_interface="in" name="RS0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="15">
      <apply>
        <eq/>
        <ci> v9 </ci>
        <apply>
          <times/>
          <ci> k9 </ci>
          <ci> PY </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="16"> 
      <apply>
        <eq/> 
        <ci> k9 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <divide/>
            <ci> RS </ci>
            <apply>
              <plus/>
              <ci> RS0 </ci>
              <ci> RS </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math> 
  </component>

  <component name="v10">
    <variable units="flux" public_interface="out" name="v10"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="12.6"/> 
    <variable units="first_order_rate_constant" name="k10"/>
    <variable units="dimensionless" name="epsilon" initial_value="0.98"/>
    
    <variable units="umol_per_g" public_interface="in" name="PY"/>
    <variable units="dimensionless" public_interface="in" name="RS"/>
    <variable units="dimensionless" public_interface="in" name="RS0"/>
    <variable units="dimensionless" public_interface="in" name="AF"/>
    <variable units="dimensionless" public_interface="in" name="AF0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="17">
      <apply>
        <eq/>
        <ci> v10 </ci>
        <apply>
          <times/>
          <ci> k10 </ci>
          <ci> PY </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="18"> 
      <apply>
        <eq/> 
        <ci> k10 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> RS </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> RS0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> RS </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1 </cn>
                <ci> epsilon </ci>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> AF </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> AF0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> AF </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component> 

  <component name="v11">
    <variable units="flux" public_interface="out" name="v11"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="0.96"/> 
    <variable units="first_order_rate_constant" name="k11"/>
    
    <variable units="umol_per_g" public_interface="in" name="LA"/>
    <variable units="dimensionless" public_interface="in" name="RS"/>
    <variable units="dimensionless" public_interface="in" name="RS0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="19">
      <apply>
        <eq/>
        <ci> v11 </ci>
        <apply>
          <times/>
          <ci> k11 </ci>
          <ci> LA </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="20"> 
      <apply>
        <eq/> 
        <ci> k11 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> RS </ci>
              <cn cellml:units="dimensionless"> -1 </cn>
            </apply>
            <apply>
              <plus/>
              <apply>
                <power/>
                <ci> RS0 </ci>
                <cn cellml:units="dimensionless"> -1 </cn>
              </apply>
              <apply>
                <power/>
                <ci> RS </ci>
                <cn cellml:units="dimensionless"> -1 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component>

  <component name="v12">
    <variable units="flux" public_interface="out" name="v12"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="695.7"/> 
    <variable units="first_order_rate_constant" name="k12"/>
    <variable units="dimensionless" name="epsilon" initial_value="0.75"/>
    
    <variable units="umol_per_g" public_interface="in" name="AC"/>
    <variable units="dimensionless" public_interface="in" name="RS"/>
    <variable units="dimensionless" public_interface="in" name="RS0"/>
    <variable units="dimensionless" public_interface="in" name="PS"/>
    <variable units="dimensionless" public_interface="in" name="PS0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="21">
      <apply>
        <eq/>
        <ci> v12 </ci>
        <apply>
          <times/>
          <ci> k12 </ci>
          <ci> AC </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="22"> 
      <apply>
        <eq/> 
        <ci> k12 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> RS </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> RS0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> RS </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1 </cn>
                <ci> epsilon </ci>
              </apply>
              <apply>
                <divide/>
                <ci> PS </ci>
                <apply>
                  <plus/>
                  <ci> PS0 </ci>
                  <ci> PS </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component>  

  <component name="v13">
    <variable units="flux" public_interface="out" name="v13"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="0.455"/> 
    <variable units="first_order_rate_constant" name="k13"/>
    
    <variable units="umol_per_g" public_interface="in" name="CR"/>
    <variable units="dimensionless" public_interface="in" name="PS"/>
    <variable units="dimensionless" public_interface="in" name="PS0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="23">
      <apply>
        <eq/>
        <ci> v13 </ci>
        <apply>
          <times/>
          <ci> k13 </ci>
          <ci> CR </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="24"> 
      <apply>
        <eq/> 
        <ci> k13 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> PS </ci>
              <cn cellml:units="dimensionless"> -1 </cn>
            </apply>
            <apply>
              <plus/>
              <apply>
                <power/>
                <ci> PS0 </ci>
                <cn cellml:units="dimensionless"> -1 </cn>
              </apply>
              <apply>
                <power/>
                <ci> PS </ci>
                <cn cellml:units="dimensionless"> -1 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component>

  <component name="v14">
    <variable units="flux" public_interface="out" name="v14"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="0.455"/> 
    <variable units="first_order_rate_constant" name="k14"/>
    
    <variable units="umol_per_g" public_interface="in" name="PC"/>
    <variable units="dimensionless" public_interface="in" name="PS"/>
    <variable units="dimensionless" public_interface="in" name="PS0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="25">
      <apply>
        <eq/>
        <ci> v14 </ci>
        <apply>
          <times/>
          <ci> k14 </ci>
          <ci> PC </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="26"> 
      <apply>
        <eq/> 
        <ci> k14 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <divide/>
            <ci> PS </ci>
            <apply>
              <plus/>
              <ci> PS0 </ci>
              <ci> PS </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component> 

  <component name="v15">
    <variable units="flux" public_interface="out" name="v15"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="626.1"/> 
    <variable units="first_order_rate_constant" name="k15"/>
    <variable units="dimensionless" name="epsilon" initial_value="0.669"/>
    
    <variable units="umol_per_g" public_interface="in" name="CoA_pool"/>
    <variable units="dimensionless" public_interface="in" name="RS"/>
    <variable units="dimensionless" public_interface="in" name="RS0"/>
    <variable units="umol_per_g" public_interface="in" name="FC"/>
    <variable units="umol_per_g" public_interface="in" name="FC0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="27">
      <apply>
        <eq/>
        <ci> v15 </ci>
        <apply>
          <times/>
          <ci> k15 </ci>
          <ci> CoA_pool </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="28"> 
      <apply>
        <eq/> 
        <ci> k15 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> RS </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> RS0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> RS </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1 </cn>
                <ci> epsilon </ci>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> FC </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> FC0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> FC </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component> 

  <component name="v16">
    <variable units="flux" public_interface="out" name="v16"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="67.0"/> 
    <variable units="first_order_rate_constant" name="k16"/>
    <variable units="dimensionless" name="epsilon" initial_value="0.775"/>
    
    <variable units="umol_per_g" public_interface="in" name="AC"/>
    <variable units="dimensionless" public_interface="in" name="RS"/>
    <variable units="dimensionless" public_interface="in" name="RS0"/>
    <variable units="umol_per_g" public_interface="in" name="FC"/>
    <variable units="umol_per_g" public_interface="in" name="FC0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="29">
      <apply>
        <eq/>
        <ci> v16 </ci>
        <apply>
          <times/>
          <ci> k16 </ci>
          <ci> AC </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="30"> 
      <apply>
        <eq/> 
        <ci> k16 </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> epsilon </ci>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> RS </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> RS0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> RS </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1 </cn>
                <ci> epsilon </ci>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <power/>
                  <ci> FC </ci>
                  <cn cellml:units="dimensionless"> -1 </cn>
                </apply>
                <apply>
                  <plus/>
                  <apply>
                    <power/>
                    <ci> FC0 </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                  <apply>
                    <power/>
                    <ci> FC </ci>
                    <cn cellml:units="dimensionless"> -1 </cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component>                                                 
  
  <component name="v17">
    <variable units="flux" public_interface="out" name="v17"/>
    
    <variable units="first_order_rate_constant" name="Vmax" initial_value="7.76"/> 
    <variable units="first_order_rate_constant" name="k17"/>
    
    <variable units="umol_per_g" public_interface="in" name="ATP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="31">
      <apply>
        <eq/>
        <ci> v17 </ci>
        <apply>
          <times/>
          <ci> k17 </ci>
          <ci> ATP </ci>
        </apply>
      </apply>    
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="32"> 
      <apply>
        <eq/> 
        <ci> k17 </ci> 
        <ci> Vmax </ci>
      </apply>
    </math> 
  </component>

  <!-- substrate concentrations -->

  <component name="GL" cmeta:id="GL">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="GL">
        <dc:title>GL</dc:title>
        <dcterms:alternative>glucose</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="GL" initial_value="0.998"/>

    <variable units="umol_per_ml" name="aGL" initial_value="4.0"/>
    <variable units="g_per_ml" name="sigmaGL" initial_value="3.76"/>

    <variable units="ml_per_g_min" public_interface="in" name="F"/> 
    <variable units="flux" public_interface="in" name="v1"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="33"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> GL </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> F </ci>
            <apply>
              <minus/>
              <ci> aGL </ci>
              <apply>
                <times/>
                <ci> sigmaGL </ci>
                <ci> GL </ci>
              </apply>
            </apply> 
          </apply> 
          <ci> v1 </ci>
        </apply>
      </apply>  
    </math>  
  </component>

  <component name="FA" cmeta:id="FA">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#FA">
        <dc:title>FA</dc:title>
        <dcterms:alternative>fatty acids</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="FA" initial_value="0.021"/>

    <variable units="umol_per_ml" name="aFA" initial_value="0.5"/>
    <variable units="g_per_ml" name="sigmaFA" initial_value="13.2"/>

    <variable units="ml_per_g_min" public_interface="in" name="F"/> 
    <variable units="flux" public_interface="in" name="v8"/>
    <variable units="flux" public_interface="in" name="v3"/>
    <variable units="flux" public_interface="in" name="v4"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="34"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> FA </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3 </cn>
              <ci> v8 </ci>
            </apply>
            <apply>
              <times/>
              <ci> F </ci>
              <apply>
                <minus/>
                <ci> aFA </ci>
                <apply>
                  <times/>
                  <ci> sigmaFA </ci>
                  <ci> FA </ci>
                </apply>
              </apply> 
            </apply>
          </apply> 
          <apply>
            <plus/>
            <ci> v3 </ci>
            <ci> v4 </ci>
          </apply>
        </apply>
      </apply>  
    </math>  
  </component>

  <component name="GP" cmeta:id="GP">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#GP">
        <dc:title>GP</dc:title>
        <dcterms:alternative>glucose-6-phosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="GP" initial_value="0.171"/>
    
    <variable units="flux" public_interface="in" name="v1"/>
    <variable units="flux" public_interface="in" name="v5"/>
    <variable units="flux" public_interface="in" name="v6"/> 
    <variable units="flux" public_interface="in" name="v7"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="35"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> GP </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> v1 </ci>
            <ci> v7 </ci>
          </apply>
          <apply>
            <plus/>
            <ci> v5 </ci>
            <ci> v6 </ci>
          </apply>
        </apply> 
      </apply> 
    </math>  
  </component>

  <component name="GY" cmeta:id="GY">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#GY">
        <dc:title>GY</dc:title>
        <dcterms:alternative>glycogen</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="GY" initial_value="33.0"/>
    
    <variable units="flux" public_interface="in" name="v6"/> 
    <variable units="flux" public_interface="in" name="v7"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="36"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> GY </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> v6 </ci>
          <ci> v7 </ci>
        </apply> 
      </apply> 
    </math>  
  </component>

  <component name="TG" cmeta:id="TG">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#TG">
        <dc:title>TG</dc:title>
        <dcterms:alternative>triglyceride</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="TG" initial_value="3.96"/>
    
    <variable units="flux" public_interface="in" name="v3"/> 
    <variable units="flux" public_interface="in" name="v8"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="37"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> TG </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless"> 1 </cn>
              <cn cellml:units="dimensionless"> 3 </cn>
            </apply>
            <ci> v3 </ci>
          </apply>
          <ci> v8 </ci>
        </apply> 
      </apply> 
    </math>  
  </component>

  <component name="PY" cmeta:id="PY">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#PY">
        <dc:title>PY</dc:title>
        <dcterms:alternative>pyruvate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="PY" initial_value="0.20"/>
    
    <variable units="flux" public_interface="in" name="v5"/>
    <variable units="flux" public_interface="in" name="v11"/>
    <variable units="flux" public_interface="in" name="v9"/> 
    <variable units="flux" public_interface="in" name="v10"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="38"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> PY </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2 </cn>
              <ci> v5 </ci>
            </apply>
            <ci> v11 </ci>
          </apply>
          <apply>
            <plus/>
            <ci> v9 </ci>
            <ci> v10 </ci>
          </apply>
        </apply> 
      </apply> 
    </math>  
  </component>

  <component name="LA" cmeta:id="LA">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#LA">
        <dc:title>LA</dc:title>
        <dcterms:alternative>lactate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="LA" initial_value="1.98"/>

    <variable units="umol_per_ml" name="aLA" initial_value="1.8"/>
    <variable units="g_per_ml" name="sigmaLA" initial_value="0.51"/>

    <variable units="ml_per_g_min" public_interface="in" name="F"/> 
    <variable units="flux" public_interface="in" name="v9"/>
    <variable units="flux" public_interface="in" name="v11"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="39"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> LA </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> v9 </ci>
            <apply>
              <times/>
              <ci> F </ci>
              <apply>
                <minus/>
                <ci> aLA </ci>
                <apply>
                  <times/>
                  <ci> sigmaLA </ci>
                  <ci> LA </ci>
                </apply>
              </apply> 
            </apply>
          </apply> 
          <ci> v11 </ci>
        </apply>
      </apply>  
    </math>  
  </component>

  <component name="AC" cmeta:id="AC">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#AC">
        <dc:title>AC</dc:title>
        <dcterms:alternative>acetyl-CoA</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="AC" initial_value="0.0046"/>
    
    <variable units="flux" public_interface="in" name="v10"/>
    <variable units="flux" public_interface="in" name="v4"/>
    <variable units="flux" public_interface="in" name="v12"/> 
    <variable units="flux" public_interface="in" name="v16"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="40"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> AC </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> v10 </ci>
            <ci> v4 </ci>
          </apply>
          <apply>
            <plus/>
            <ci> v12 </ci>
            <ci> v16 </ci>
          </apply>
        </apply> 
      </apply> 
    </math>  
  </component>

  <component name="FC" cmeta:id="FC">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#FC">
        <dc:title>FC</dc:title>
        <dcterms:alternative>free-CoA</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="FC" initial_value="0.0088"/>
    
    <variable units="flux" public_interface="in" name="v10"/>
    <variable units="flux" public_interface="in" name="v4"/>
    <variable units="flux" public_interface="in" name="v12"/> 
    <variable units="flux" public_interface="in" name="v15"/>
    <variable units="flux" public_interface="in" name="v16"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="40a"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> FC </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> v12 </ci>
            <ci> v15 </ci>
            <ci> v16 </ci>
          </apply>
          <apply>
            <plus/>
            <ci> v10 </ci>
            <ci> v4 </ci>
          </apply>
        </apply> 
      </apply> 
    </math>  
  </component>

  <component name="CoA_pool" cmeta:id="CoA_pool">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#CoA_pool">
        <dc:title>CoA_pool</dc:title>
        <dcterms:alternative>remainder of the CoA species</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="CoA_pool" initial_value="0.043"/>
    
    <variable units="flux" public_interface="in" name="v15"/>
    <variable units="flux" public_interface="in" name="v16"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="40b"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> CoA_pool </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> v16 </ci>
          <ci> v15 </ci>
        </apply> 
      </apply> 
    </math>  
  </component>

  <component name="O2" cmeta:id="O2">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#O2">
        <dc:title>O2</dc:title>
        <dcterms:alternative>oxygen</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="O2" initial_value="0.963"/>

    <variable units="umol_per_ml" name="aO2" initial_value="6.53"/>
    <variable units="g_per_ml" name="sigmaO2" initial_value="1.0"/>

    <variable units="ml_per_g_min" public_interface="in" name="F"/> 
    <variable units="flux" public_interface="in" name="v2"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="41"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> O2 </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> F </ci>
            <apply>
              <minus/>
              <ci> aO2 </ci>
              <apply>
                <times/>
                <ci> sigmaO2 </ci>
                <ci> O2 </ci>
              </apply>
            </apply> 
          </apply>
          <ci> v2 </ci>
        </apply>
      </apply>  
    </math>  
  </component>

  <component name="CO2" cmeta:id="CO2">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#CO2">
        <dc:title>CO2</dc:title>
        <dcterms:alternative>carbon dioxide</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="CO2" initial_value="20.0"/>

    <variable units="umol_per_ml" name="aCO2" initial_value="15.5"/>
    <variable units="g_per_ml" name="sigmaCO2" initial_value="1.0"/>

    <variable units="ml_per_g_min" public_interface="in" name="F"/> 
    <variable units="flux" public_interface="in" name="v10"/>
    <variable units="flux" public_interface="in" name="v12"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="42"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> CO2 </ci> 
        </apply>
        <apply>
          <plus/>
          <ci> v10 </ci>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2 </cn>
            <ci> v12 </ci>
          </apply>
          <apply>
            <times/>
            <ci> F </ci>
            <apply>
              <minus/>
              <ci> aCO2 </ci>
              <apply>
                <times/>
                <ci> sigmaCO2 </ci>
                <ci> CO2 </ci>
              </apply>
            </apply> 
          </apply>
        </apply>
      </apply>  
    </math>  
  </component>

  <component name="NAD"> 
    <variable units="umol_per_g" public_interface="out" name="NAD" initial_value="1.81"/>
    
    <variable units="flux" public_interface="in" name="v5"/>
    <variable units="flux" public_interface="in" name="v11"/>
    <variable units="flux" public_interface="in" name="v9"/> 
    <variable units="flux" public_interface="in" name="v10"/>
    <variable units="flux" public_interface="in" name="v12"/>
    <variable units="flux" public_interface="in" name="v4"/>
    <variable units="flux" public_interface="in" name="v2"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="43"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> NAD </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> v9 </ci>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2 </cn>
              <ci> v2 </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2 </cn>
              <ci> v5 </ci>
            </apply>
            <ci> v10 </ci>
            <ci> v11 </ci>
            <apply>
              <times/>
              <apply>
                <divide/>
                <cn cellml:units="dimensionless"> 11 </cn>
                <cn cellml:units="dimensionless"> 3 </cn>
              </apply>
              <ci> v12 </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <divide/>
                <cn cellml:units="dimensionless"> 35 </cn>
                <cn cellml:units="dimensionless"> 3 </cn>
              </apply>
              <ci> v4 </ci>
            </apply>
          </apply>
        </apply> 
      </apply> 
    </math>  
  </component>

  <component name="ADP"> 
    <variable units="umol_per_g" public_interface="out" name="ADP" initial_value="0.90"/>
    
    <variable units="flux" public_interface="in" name="v1"/>
    <variable units="flux" public_interface="in" name="v5"/>
    <variable units="flux" public_interface="in" name="v6"/> 
    <variable units="flux" public_interface="in" name="v12"/>
    <variable units="flux" public_interface="in" name="v4"/>
    <variable units="flux" public_interface="in" name="v3"/>
    <variable units="flux" public_interface="in" name="v2"/>
    <variable units="flux" public_interface="in" name="v14"/>
    <variable units="flux" public_interface="in" name="v13"/>
    <variable units="flux" public_interface="in" name="v17"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="44"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> ADP </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> v1 </ci>
            <ci> v6 </ci>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2 </cn>
              <ci> v4 </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2 </cn>
              <ci> v3 </ci>
            </apply>
            <ci> v13 </ci>
            <ci> v17 </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 3 </cn>
              <ci> v5 </ci>
            </apply>
            <ci> v12 </ci>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 6 </cn>
              <ci> v2 </ci>
            </apply>
            <ci> v14 </ci>
          </apply>
        </apply> 
      </apply> 
    </math>  
  </component>

  <component name="PC" cmeta:id="PC">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#PC">
        <dc:title>PC</dc:title>
        <dcterms:alternative>phosphocreatine</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="PC" initial_value="8.80"/>

    <variable units="flux" public_interface="in" name="v13"/>
    <variable units="flux" public_interface="in" name="v14"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="45"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> PC </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> v13 </ci>
          <ci> v14 </ci>
        </apply>
      </apply>  
    </math>  
  </component>

  <component name="CR" cmeta:id="CR">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#CR">
        <dc:title>CR</dc:title>
        <dcterms:alternative>creatine</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="umol_per_g" public_interface="out" name="CR" initial_value="3.5"/>

    <variable units="flux" public_interface="in" name="v13"/>
    <variable units="flux" public_interface="in" name="v14"/>
    <variable units="minute" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="45a"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> CR </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> v14 </ci>
          <ci> v13 </ci>
        </apply>
      </apply>  
    </math>  
  </component>

  <component name="ATP">
    <variable units="umol_per_g" public_interface="out" name="ATP" initial_value="4.5"/>
  </component>

  <component name="NADH">
    <variable units="umol_per_g" public_interface="out" name="NADH" initial_value="0.19"/>
  </component>

  <component name="CS0">
    <variable units="dimensionless" public_interface="out" name="CS0" initial_value="1.0"/>
  </component>

  <component name="PS0">
    <variable units="dimensionless" public_interface="out" name="PS0" initial_value="0.2"/>
  </component>

  <component name="RS0">
    <variable units="dimensionless" public_interface="out" name="RS0" initial_value="0.111"/>
  </component>

  <component name="AF0">
    <variable units="dimensionless" public_interface="out" name="AF0" initial_value="0.523"/>
  </component>

  <component name="FC0">
    <variable units="umol_per_g" public_interface="out" name="FC0" initial_value="0.0088"/>
  </component>

  <component name="CS"> 
    <variable units="dimensionless" public_interface="out" name="CS"/>

    <variable units="umol_per_g" public_interface="in" name="CR"/>
    <variable units="umol_per_g" public_interface="in" name="PC"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="46"> 
      <apply>
        <eq/> 
        <ci> CS </ci> 
        <apply>
          <divide/>
          <ci> CR </ci>
          <ci> PC </ci>
        </apply>
      </apply>  
    </math>  
  </component> 

  <component name="RS"> 
    <variable units="dimensionless" public_interface="out" name="RS"/>

    <variable units="umol_per_g" public_interface="in" name="NADH"/>
    <variable units="umol_per_g" public_interface="in" name="NAD"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="47"> 
      <apply>
        <eq/> 
        <ci> RS </ci> 
        <apply>
          <divide/>
          <ci> NADH </ci>
          <ci> NAD </ci>
        </apply>
      </apply>  
    </math>  
  </component> 

  <component name="PS"> 
    <variable units="dimensionless" public_interface="out" name="PS"/>

    <variable units="umol_per_g" public_interface="in" name="ADP"/>
    <variable units="umol_per_g" public_interface="in" name="ATP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="48"> 
      <apply>
        <eq/> 
        <ci> PS </ci> 
        <apply>
          <divide/>
          <ci> ADP </ci>
          <ci> ATP </ci>
        </apply>
      </apply>  
    </math>  
  </component> 

  <component name="AF"> 
    <variable units="dimensionless" public_interface="out" name="AF"/>

    <variable units="umol_per_g" public_interface="in" name="AC"/>
    <variable units="umol_per_g" public_interface="in" name="FC"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="49"> 
      <apply>
        <eq/> 
        <ci> AF </ci> 
        <apply>
          <divide/>
          <ci> AC </ci>
          <ci> FC </ci>
        </apply>
      </apply>  
    </math>  
  </component> 

  <component name="model_parameters">
    <variable units="ml_per_g_min" public_interface="out" name="F"/>
    <variable units="minute" public_interface="in" name="time"/>

    <math xmlns="http://www.w3.org/1998/Math/MathML" id="50"> 
      <apply>
        <eq/>
        <ci> F </ci>
        <piecewise>
          <piece>
            <cn cellml:units="ml_per_g_min"> 1.0 </cn>
            <apply>
              <and/>
              <apply>
                <gt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 0.0 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="minute"> 5.0 </cn>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="ml_per_g_min"> 0.4 </cn>
          </otherwise>
        </piecewise>
      </apply>
    </math>
  </component>


  <connection>
    <map_components component_2="GL" component_1="v1"/>
    <map_variables variable_2="v1" variable_1="v1"/>
    <map_variables variable_2="GL" variable_1="GL"/>
  </connection>

  <connection>
    <map_components component_2="O2" component_1="v2"/>
    <map_variables variable_2="v2" variable_1="v2"/>
    <map_variables variable_2="O2" variable_1="O2"/>
  </connection>

  <connection>
    <map_components component_2="FA" component_1="v3"/>
    <map_variables variable_2="v3" variable_1="v3"/>
    <map_variables variable_2="FA" variable_1="FA"/>
  </connection>

  <connection>
    <map_components component_2="PS" component_1="v3"/>
    <map_variables variable_2="PS" variable_1="PS"/>
  </connection>

  <connection>
    <map_components component_2="PS0" component_1="v3"/>
    <map_variables variable_2="PS0" variable_1="PS0"/>
  </connection>

  <connection>
    <map_components component_2="FA" component_1="v4"/>
    <map_variables variable_2="v4" variable_1="v4"/>
    <map_variables variable_2="FA" variable_1="FA"/>
  </connection>

  <connection>
    <map_components component_2="RS" component_1="v4"/>
    <map_variables variable_2="RS" variable_1="RS"/>
  </connection>

  <connection>
    <map_components component_2="RS0" component_1="v4"/>
    <map_variables variable_2="RS0" variable_1="RS0"/>
  </connection>

  <connection>
    <map_components component_2="AF" component_1="v4"/>
    <map_variables variable_2="AF" variable_1="AF"/>
  </connection>

  <connection>
    <map_components component_2="AF0" component_1="v4"/>
    <map_variables variable_2="AF0" variable_1="AF0"/>
  </connection>

  <connection>
    <map_components component_2="GP" component_1="v5"/>
    <map_variables variable_2="v5" variable_1="v5"/>
    <map_variables variable_2="GP" variable_1="GP"/>
  </connection>

  <connection>
    <map_components component_2="RS" component_1="v5"/>
    <map_variables variable_2="RS" variable_1="RS"/>
  </connection>

  <connection>
    <map_components component_2="RS0" component_1="v5"/>
    <map_variables variable_2="RS0" variable_1="RS0"/>
  </connection>

  <connection>
    <map_components component_2="PS" component_1="v5"/>
    <map_variables variable_2="PS" variable_1="PS"/>
  </connection>

  <connection>
    <map_components component_2="PS0" component_1="v5"/>
    <map_variables variable_2="PS0" variable_1="PS0"/>
  </connection>

  <connection>
    <map_components component_2="GP" component_1="v6"/>
    <map_variables variable_2="v6" variable_1="v6"/>
    <map_variables variable_2="GP" variable_1="GP"/>
  </connection>

  <connection>
    <map_components component_2="CS" component_1="v6"/>
    <map_variables variable_2="CS" variable_1="CS"/>
  </connection>

  <connection>
    <map_components component_2="CS0" component_1="v6"/>
    <map_variables variable_2="CS0" variable_1="CS0"/>
  </connection>

  <connection>
    <map_components component_2="PS" component_1="v6"/>
    <map_variables variable_2="PS" variable_1="PS"/>
  </connection>

  <connection>
    <map_components component_2="PS0" component_1="v6"/>
    <map_variables variable_2="PS0" variable_1="PS0"/>
  </connection>

  <connection>
    <map_components component_2="GY" component_1="v7"/>
    <map_variables variable_2="v7" variable_1="v7"/>
    <map_variables variable_2="GY" variable_1="GY"/>
  </connection>

  <connection>
    <map_components component_2="CS" component_1="v7"/>
    <map_variables variable_2="CS" variable_1="CS"/>
  </connection>

  <connection>
    <map_components component_2="CS0" component_1="v7"/>
    <map_variables variable_2="CS0" variable_1="CS0"/>
  </connection>

  <connection>
    <map_components component_2="PS" component_1="v7"/>
    <map_variables variable_2="PS" variable_1="PS"/>
  </connection>

  <connection>
    <map_components component_2="PS0" component_1="v7"/>
    <map_variables variable_2="PS0" variable_1="PS0"/>
  </connection>

  <connection>
    <map_components component_2="TG" component_1="v8"/>
    <map_variables variable_2="v8" variable_1="v8"/>
    <map_variables variable_2="TG" variable_1="TG"/>
  </connection>

  <connection>
    <map_components component_2="PY" component_1="v9"/>
    <map_variables variable_2="v9" variable_1="v9"/>
    <map_variables variable_2="PY" variable_1="PY"/>
  </connection>

  <connection>
    <map_components component_2="RS" component_1="v9"/>
    <map_variables variable_2="RS" variable_1="RS"/>
  </connection>

  <connection>
    <map_components component_2="RS0" component_1="v9"/>
    <map_variables variable_2="RS0" variable_1="RS0"/>
  </connection>

  <connection>
    <map_components component_2="PY" component_1="v10"/>
    <map_variables variable_2="v10" variable_1="v10"/>
    <map_variables variable_2="PY" variable_1="PY"/>
  </connection>

  <connection>
    <map_components component_2="RS" component_1="v10"/>
    <map_variables variable_2="RS" variable_1="RS"/>
  </connection>

  <connection>
    <map_components component_2="RS0" component_1="v10"/>
    <map_variables variable_2="RS0" variable_1="RS0"/>
  </connection>

  <connection>
    <map_components component_2="AF" component_1="v10"/>
    <map_variables variable_2="AF" variable_1="AF"/>
  </connection>

  <connection>
    <map_components component_2="AF0" component_1="v10"/>
    <map_variables variable_2="AF0" variable_1="AF0"/>
  </connection>

  <connection>
    <map_components component_2="LA" component_1="v11"/>
    <map_variables variable_2="v11" variable_1="v11"/>
    <map_variables variable_2="LA" variable_1="LA"/>
  </connection>

  <connection>
    <map_components component_2="RS" component_1="v11"/>
    <map_variables variable_2="RS" variable_1="RS"/>
  </connection>

  <connection>
    <map_components component_2="RS0" component_1="v11"/>
    <map_variables variable_2="RS0" variable_1="RS0"/>
  </connection>

  <connection>
    <map_components component_2="AC" component_1="v12"/>
    <map_variables variable_2="v12" variable_1="v12"/>
    <map_variables variable_2="AC" variable_1="AC"/>
  </connection>

  <connection>
    <map_components component_2="RS" component_1="v12"/>
    <map_variables variable_2="RS" variable_1="RS"/>
  </connection>

  <connection>
    <map_components component_2="RS0" component_1="v12"/>
    <map_variables variable_2="RS0" variable_1="RS0"/>
  </connection>

  <connection>
    <map_components component_2="PS" component_1="v12"/>
    <map_variables variable_2="PS" variable_1="PS"/>
  </connection>

  <connection>
    <map_components component_2="PS0" component_1="v12"/>
    <map_variables variable_2="PS0" variable_1="PS0"/>
  </connection>

  <connection>
    <map_components component_2="CR" component_1="v13"/>
    <map_variables variable_2="v13" variable_1="v13"/>
    <map_variables variable_2="CR" variable_1="CR"/>
  </connection>

  <connection>
    <map_components component_2="PS" component_1="v13"/>
    <map_variables variable_2="PS" variable_1="PS"/>
  </connection>

  <connection>
    <map_components component_2="PS0" component_1="v13"/>
    <map_variables variable_2="PS0" variable_1="PS0"/>
  </connection>

  <connection>
    <map_components component_2="PC" component_1="v14"/>
    <map_variables variable_2="v14" variable_1="v14"/>
    <map_variables variable_2="PC" variable_1="PC"/>
  </connection>

  <connection>
    <map_components component_2="PS" component_1="v14"/>
    <map_variables variable_2="PS" variable_1="PS"/>
  </connection>

  <connection>
    <map_components component_2="PS0" component_1="v14"/>
    <map_variables variable_2="PS0" variable_1="PS0"/>
  </connection>

  <connection>
    <map_components component_2="CoA_pool" component_1="v15"/>
    <map_variables variable_2="v15" variable_1="v15"/>
    <map_variables variable_2="CoA_pool" variable_1="CoA_pool"/>
  </connection>

  <connection>
    <map_components component_2="RS" component_1="v15"/>
    <map_variables variable_2="RS" variable_1="RS"/>
  </connection>

  <connection>
    <map_components component_2="RS0" component_1="v15"/>
    <map_variables variable_2="RS0" variable_1="RS0"/>
  </connection>

  <connection>
    <map_components component_2="FC" component_1="v15"/>
    <map_variables variable_2="v15" variable_1="v15"/>
    <map_variables variable_2="FC" variable_1="FC"/>
  </connection>

  <connection>
    <map_components component_2="FC0" component_1="v15"/>
    <map_variables variable_2="FC0" variable_1="FC0"/>
  </connection>

  <connection>
    <map_components component_2="AC" component_1="v16"/>
    <map_variables variable_2="v16" variable_1="v16"/>
    <map_variables variable_2="AC" variable_1="AC"/>
  </connection>

  <connection>
    <map_components component_2="RS" component_1="v16"/>
    <map_variables variable_2="RS" variable_1="RS"/>
  </connection>

  <connection>
    <map_components component_2="RS0" component_1="v16"/>
    <map_variables variable_2="RS0" variable_1="RS0"/>
  </connection>

  <connection>
    <map_components component_2="FC" component_1="v16"/>
    <map_variables variable_2="v16" variable_1="v16"/>
    <map_variables variable_2="FC" variable_1="FC"/>
  </connection>

  <connection>
    <map_components component_2="FC0" component_1="v16"/>
    <map_variables variable_2="FC0" variable_1="FC0"/>
  </connection>

  <connection>
    <map_components component_2="ATP" component_1="v17"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="GL"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="FA"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="O2"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="CO2"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="LA"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="v1" component_1="GP"/>
    <map_variables variable_2="v1" variable_1="v1"/>
  </connection>

  <connection>
    <map_components component_2="v7" component_1="GP"/>
    <map_variables variable_2="v7" variable_1="v7"/>
  </connection>

  <connection>
    <map_components component_2="v6" component_1="GY"/>
    <map_variables variable_2="v6" variable_1="v6"/>
  </connection>

  <connection>
    <map_components component_2="v3" component_1="TG"/>
    <map_variables variable_2="v3" variable_1="v3"/>
  </connection>

  <connection>
    <map_components component_2="v5" component_1="PY"/>
    <map_variables variable_2="v5" variable_1="v5"/>
  </connection>

  <connection>
    <map_components component_2="v11" component_1="PY"/>
    <map_variables variable_2="v11" variable_1="v11"/>
  </connection>

  <connection>
    <map_components component_2="v9" component_1="LA"/>
    <map_variables variable_2="v9" variable_1="v9"/>
  </connection>

  <connection>
    <map_components component_2="v10" component_1="AC"/>
    <map_variables variable_2="v10" variable_1="v10"/>
  </connection>

  <connection>
    <map_components component_2="v4" component_1="AC"/>
    <map_variables variable_2="v4" variable_1="v4"/>
  </connection>

  <connection>
    <map_components component_2="v10" component_1="FC"/>
    <map_variables variable_2="v10" variable_1="v10"/>
  </connection>

  <connection>
    <map_components component_2="v4" component_1="FC"/>
    <map_variables variable_2="v4" variable_1="v4"/>
  </connection>

  <connection>
    <map_components component_2="v12" component_1="FC"/>
    <map_variables variable_2="v12" variable_1="v12"/>
  </connection>
  
  <connection>
    <map_components component_2="v8" component_1="FA"/>
    <map_variables variable_2="v8" variable_1="v8"/>
  </connection>

  <connection>
    <map_components component_2="v16" component_1="CoA_pool"/>
    <map_variables variable_2="v16" variable_1="v16"/>
  </connection>

  <connection>
    <map_components component_2="v10" component_1="CO2"/>
    <map_variables variable_2="v10" variable_1="v10"/>
  </connection>

  <connection>
    <map_components component_2="v12" component_1="CO2"/>
    <map_variables variable_2="v12" variable_1="v12"/>
  </connection>

  <connection>
    <map_components component_2="v10" component_1="NAD"/>
    <map_variables variable_2="v10" variable_1="v10"/>
  </connection>

  <connection>
    <map_components component_2="v11" component_1="NAD"/>
    <map_variables variable_2="v11" variable_1="v11"/>
  </connection>

  <connection>
    <map_components component_2="v5" component_1="NAD"/>
    <map_variables variable_2="v5" variable_1="v5"/>
  </connection>

  <connection>
    <map_components component_2="v9" component_1="NAD"/>
    <map_variables variable_2="v9" variable_1="v9"/>
  </connection>

  <connection>
    <map_components component_2="v12" component_1="NAD"/>
    <map_variables variable_2="v12" variable_1="v12"/>
  </connection>

  <connection>
    <map_components component_2="v4" component_1="NAD"/>
    <map_variables variable_2="v4" variable_1="v4"/>
  </connection>

  <connection>
    <map_components component_2="v2" component_1="NAD"/>
    <map_variables variable_2="v2" variable_1="v2"/>
  </connection>

  <connection>
    <map_components component_2="v1" component_1="ADP"/>
    <map_variables variable_2="v1" variable_1="v1"/>
  </connection>

  <connection>
    <map_components component_2="v5" component_1="ADP"/>
    <map_variables variable_2="v5" variable_1="v5"/>
  </connection>

  <connection>
    <map_components component_2="v6" component_1="ADP"/>
    <map_variables variable_2="v6" variable_1="v6"/>
  </connection>

  <connection>
    <map_components component_2="v12" component_1="ADP"/>
    <map_variables variable_2="v12" variable_1="v12"/>
  </connection>

  <connection>
    <map_components component_2="v4" component_1="ADP"/>
    <map_variables variable_2="v4" variable_1="v4"/>
  </connection>

  <connection>
    <map_components component_2="v3" component_1="ADP"/>
    <map_variables variable_2="v3" variable_1="v3"/>
  </connection>

  <connection>
    <map_components component_2="v2" component_1="ADP"/>
    <map_variables variable_2="v2" variable_1="v2"/>
  </connection>

  <connection>
    <map_components component_2="v14" component_1="ADP"/>
    <map_variables variable_2="v14" variable_1="v14"/>
  </connection>

  <connection>
    <map_components component_2="v13" component_1="ADP"/>
    <map_variables variable_2="v13" variable_1="v13"/>
  </connection>

  <connection>
    <map_components component_2="v17" component_1="ADP"/>
    <map_variables variable_2="v17" variable_1="v17"/>
  </connection>

  <connection>
    <map_components component_2="v13" component_1="PC"/>
    <map_variables variable_2="v13" variable_1="v13"/>
  </connection>

  <connection>
    <map_components component_2="v14" component_1="CR"/>
    <map_variables variable_2="v14" variable_1="v14"/>
  </connection>

  <connection>
    <map_components component_2="CR" component_1="CS"/>
    <map_variables variable_2="CR" variable_1="CR"/>
  </connection>

  <connection>
    <map_components component_2="PC" component_1="CS"/>
    <map_variables variable_2="PC" variable_1="PC"/>
  </connection>

  <connection>
    <map_components component_2="NAD" component_1="RS"/>
    <map_variables variable_2="NAD" variable_1="NAD"/>
  </connection>

  <connection>
    <map_components component_2="NADH" component_1="RS"/>
    <map_variables variable_2="NADH" variable_1="NADH"/>
  </connection>

  <connection>
    <map_components component_2="ADP" component_1="PS"/>
    <map_variables variable_2="ADP" variable_1="ADP"/>
  </connection>

  <connection>
    <map_components component_2="ATP" component_1="PS"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
  </connection>

  <connection>
    <map_components component_2="AC" component_1="AF"/>
    <map_variables variable_2="AC" variable_1="AC"/>
  </connection>

  <connection>
    <map_components component_2="FC" component_1="AF"/>
    <map_variables variable_2="FC" variable_1="FC"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="GL"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="FA"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="GP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="GY"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="TG"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="PY"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="LA"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="AC"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="FC"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="CoA_pool"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="O2"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="CO2"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="NAD"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="ADP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="PC"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="CR"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

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




<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#ff3cdf4a-868a-4db1-a3e1-c96d25414785">
    <rdf:li>cardiac</rdf:li>
    <rdf:li>myocardium</rdf:li>
    <rdf:li>energy</rdf:li>
    <rdf:li>metabolism</rdf:li>
    <rdf:li>ischemic</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#e7987719-3527-4179-b3db-9c5ca707c2c3">
    <rdf:li rdf:resource="rdf:#d07b288f-2a02-42a1-a4da-b5b75c096636"/>
    <rdf:li rdf:resource="rdf:#b1557055-2cbf-47a1-9f67-c98b12ed9e2c"/>
    <rdf:li rdf:resource="rdf:#d7b92f77-9aca-4eb4-b56f-86f05d2f5ee7"/>
    <rdf:li rdf:resource="rdf:#753efb47-36b1-4a76-8ce7-2a9f436a8e04"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#24e76f24-cd3c-46de-b6e7-2abe0d311b7f">
    <bqs:Pubmed_id>11962772</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#7192cbf1-f3e7-4fc2-8c8c-ff5c210f1563"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Auckland Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#f68ec13e-c811-4e97-bfd7-1d7714b683fe"/>
    <dcterms:created rdf:resource="rdf:#ec4348b2-c82a-40ca-8400-bf844a5a1aa9"/>
    <dc:creator rdf:resource="rdf:#396467af-6cc1-4809-b494-8e189bb7ae58"/>
    <cmeta:modification rdf:resource="rdf:#9f6bd8c2-48d6-4bbd-9ff7-4211e35cd1a7"/>
    <cmeta:modification rdf:resource="rdf:#ca7c087a-2c9a-4652-a222-e48ba7d3b471"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2b34792b-703b-43db-b4d0-b01fc61af286">
    <vCard:N rdf:resource="rdf:#686cbbab-3783-4607-9686-d9e6815ab594"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#753efb47-36b1-4a76-8ce7-2a9f436a8e04">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#d4b05f00-1e2c-42f1-857f-64ed3631339a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ec4348b2-c82a-40ca-8400-bf844a5a1aa9">
    <dcterms:W3CDTF>2002-08-20T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d07b288f-2a02-42a1-a4da-b5b75c096636">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#9691c86f-b7ee-490f-841c-e81729d5f0ea"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2dafe33a-bf85-4de2-bdbd-bec4dc39c361">
    <dcterms:W3CDTF>2002-02</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b9bb65e9-3484-4597-a476-28bc18f1dae3">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2af03eef-c410-43f5-bb0b-a94cd99facc9">
    <vCard:Given>Gerald</vCard:Given>
    <vCard:Family>Saidel</vCard:Family>
    <vCard:Other>M</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7fd515d7-3830-435d-867a-a96f86a7070c">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#ff3cdf4a-868a-4db1-a3e1-c96d25414785"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9691c86f-b7ee-490f-841c-e81729d5f0ea">
    <vCard:Given>Jennifer</vCard:Given>
    <vCard:Family>Salem</vCard:Family>
    <vCard:Other>E</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#64d84aeb-fd5a-48b7-ab02-a8aaefd818c8">
    <dc:title>Annals of Biomedical Engineering</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#45da7636-e506-4438-b9c7-820ba3286f5c">
    <vCard:FN>Autumn Cuellar</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7192cbf1-f3e7-4fc2-8c8c-ff5c210f1563">
    <dc:creator rdf:resource="rdf:#e7987719-3527-4179-b3db-9c5ca707c2c3"/>
    <dc:title>Mechanistic Model of Myocardial Energy Metabolism Under Normal and &#13;
            Ischemic Conditions</dc:title>
    <bqs:volume>30</bqs:volume>
    <bqs:first_page>202</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#64d84aeb-fd5a-48b7-ab02-a8aaefd818c8"/>
    <dcterms:issued rdf:resource="rdf:#2dafe33a-bf85-4de2-bdbd-bec4dc39c361"/>
    <bqs:last_page>216</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c166dceb-bfbd-4f08-8e24-8909f1be51dd">
    <dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6760a3b3-896c-4722-95e2-6571489d00b6">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#51940d62-3f3e-421b-b8b7-a0045eb80fa4">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>a.cuellar@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#salem_model_2002">
    <dc:title>
        The Salem et al Mechanistic Model of Myocardial Energy Metabolism Under
        Normal and Ischemic Conditions, 2002
      </dc:title>
    <cmeta:comment rdf:resource="rdf:#5c02a91b-1538-4e13-bbf6-d2b2b507b5ce"/>
    <bqs:reference rdf:resource="rdf:#24e76f24-cd3c-46de-b6e7-2abe0d311b7f"/>
    <bqs:reference rdf:resource="rdf:#e354e0a0-6f4e-4e36-bb27-a45005af7d60"/>
    <cmeta:species>Mammalia</cmeta:species>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#686cbbab-3783-4607-9686-d9e6815ab594">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d7b92f77-9aca-4eb4-b56f-86f05d2f5ee7">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#e53ef42e-760b-45ca-89b5-8f1aee4f6670"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e354e0a0-6f4e-4e36-bb27-a45005af7d60">
    <dc:subject rdf:resource="rdf:#7fd515d7-3830-435d-867a-a96f86a7070c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b1557055-2cbf-47a1-9f67-c98b12ed9e2c">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#2af03eef-c410-43f5-bb0b-a94cd99facc9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d4b05f00-1e2c-42f1-857f-64ed3631339a">
    <vCard:Given>Marco</vCard:Given>
    <vCard:Family>Cabrera</vCard:Family>
    <vCard:Other>E</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f68ec13e-c811-4e97-bfd7-1d7714b683fe">
    <dc:creator rdf:resource="rdf:#5f776c58-0e40-4cdd-9dbc-6e4922423a18"/>
    <rdf:value>This CellML model has been recoded to remove the reaction element.  The model runs in both COR and PCEnv and we believe it is close to recreating the original published results.  The units have been checked and they are consistent.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#230c4dd0-bb67-4a8c-a3a2-db2da33b57ac">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#04575492-59b9-4440-97f0-9468b4b6ab3b">
    <vCard:N rdf:resource="rdf:#b9bb65e9-3484-4597-a476-28bc18f1dae3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9f6bd8c2-48d6-4bbd-9ff7-4211e35cd1a7">
    <dcterms:modified rdf:resource="rdf:#c166dceb-bfbd-4f08-8e24-8909f1be51dd"/>
    <rdf:value>
          Added publication date information.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#2b34792b-703b-43db-b4d0-b01fc61af286"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a389274e-93e8-4732-ac69-4f7de2afd7e9">
    <dcterms:W3CDTF>2008-01-16T09:08:29+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c02a91b-1538-4e13-bbf6-d2b2b507b5ce">
    <dc:creator rdf:resource="rdf:#45da7636-e506-4438-b9c7-820ba3286f5c"/>
    <rdf:value>This file contains a CellML description of Salem et al's 2002 mathematical model of the mechanical properties of myocardial energy metabolism under normal and ischemic conditions.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5f776c58-0e40-4cdd-9dbc-6e4922423a18">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#396467af-6cc1-4809-b494-8e189bb7ae58">
    <vCard:ORG rdf:resource="rdf:#6760a3b3-896c-4722-95e2-6571489d00b6"/>
    <vCard:EMAIL rdf:resource="rdf:#51940d62-3f3e-421b-b8b7-a0045eb80fa4"/>
    <vCard:N rdf:resource="rdf:#230c4dd0-bb67-4a8c-a3a2-db2da33b57ac"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e53ef42e-760b-45ca-89b5-8f1aee4f6670">
    <vCard:Given>William</vCard:Given>
    <vCard:Family>Stanley</vCard:Family>
    <vCard:Other>C</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ca7c087a-2c9a-4652-a222-e48ba7d3b471">
    <dcterms:modified rdf:resource="rdf:#a389274e-93e8-4732-ac69-4f7de2afd7e9"/>
    <rdf:value>Completely recoded the model to remove the reaction element.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#04575492-59b9-4440-97f0-9468b4b6ab3b"/>
  </rdf:Description>
</rdf:RDF>
</model>