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

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2009-06-17 15:51:03+12:00
Desc:
committing version02 of salem_saidel_stanley_cabrera_2002
Permanent Source URI:
https://models.physiomeproject.org/workspace/salem_saidel_stanley_cabrera_2002/rawfile/246f86dcb49c603087c31b00f98d0db7c5c76ab3/salem_saidel_stanley_cabrera_2002.cellml

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

CREATED :  20th August 2002

LAST MODIFIED :  9th April 2003

AUTHOR :  Autumn A. Cuellar
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the CellML Metadata 1.0 Specification released on 16th
January, 2002.

DESCRIPTION :  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.

CHANGES: 
  09/04/2003 - AAC - Added publication date information.  
--><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>Mathematical Mechanistic Model of Myocardial Metabolism during Ischemia</title>
  <author>
    <firstname>Autumn</firstname>
    <othername>A.</othername>
    <surname>Cuellar</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
          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.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
In myocardial ischemia, or the reduction of blood flow to the heart, myocardial oxygen consumption and pyruvate oxidation decreases, glycolysis becomes accelerated, and lactate accumulates.  To better understand myocardial metabolism under ischemic conditions, Jennifer Salem <emphasis>et al.</emphasis> developed a quantitative mechanistic model in their 2002 paper (see the figure below).  They modelled a 60% reduction of coronary blood flow and found that at the onset of ischemia, the concentration of lactate in myocytes increased, and the myocardium switched from a net consumer of lactate to a net producer.  After sixty minutes of ischemia, the glycogen concentration was reduced by 60%, phosphocreatine by 75%, and ATP by 50%.  Furthermore, their model predicts a rapid 13-fold increase in the ratio of NADH to NAD<superscript>+</superscript>, but only a twofold increase in the ratio of acetyl-CoA to free-CoA.  These predictions are consistent with the idea that pyruvate oxidation is inhibited during ischemia partially by the rise in NADH/NAD<superscript>+</superscript>.
</para>

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

<para>
<ulink url="http://ojps.aip.org/getabs/servlet/GetabsServlet?prog=normal&amp;id=ABMECF000030000002000202000001&amp;idtype=cvips&amp;gifs=Yes">Mechanistic Model of Myocardial Energy Metabolism Under Normal and Ischemic Conditions</ulink>, Jennifer E. Salem, Gerald M. Saidel, William C. Stanley, and Marco E. Cabrera, 2002, <ulink url="http://www.ruf.rice.edu/~abme/annals.html">
            <emphasis>Annals of Biomedical Engineering</emphasis>
          </ulink>, 30, 202-216. (A PDF of the article is available to subscribers on the <emphasis>Annals of Biomedical Engineering</emphasis> website.)  <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, 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>The 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>