Location: Vendelin, Kongas, Saks, 2000 @ 717bc42691a3 / vendelin_kongas_saks_2000.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2009-03-01 22:51:07+13:00
Desc:
committing version01 of vendelin_kongas_saks_2000
Permanent Source URI:
https://staging.physiomeproject.org/workspace/vendelin_kongas_saks_2000/rawfile/717bc42691a3c7ff2095aa4b60f0d59e811e1025/vendelin_kongas_saks_2000.cellml

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

CREATED :  6th November 2003

LAST MODIFIED : 21st April 2005

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/1/02 CellML Metadata 1.0 Specification.

DESCRIPTION :  This file contains a CellML description of Vendelin et al.'s mathematical model of mitochondrial respiration in heart cells.

CHANGES:  
    21/04/2005 - PJV - Changed unit dimensions to correct inconsistancy
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" 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:NS1="http://purl.org/dc/terms/" xmlns:NS2="http://purl.org/dc/elements/1.1/" xmlns:NS3="http://imc.org/vCard/3.0#" xmlns:NS4="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:NS5="http://www.cellml.org/metadata/simulation/1.0#" xmlns:NS6="http://www.cellml.org/metadata/1.0#" xmlns:NS7="http://www.cellml.org/bqs/1.0#" cmeta:id="vendelin_model_2000" name="vendelin_model_2000">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Creatine Kinase in Energy Metabolic Signalling in Muscle</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</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 checked in both PCEnv and COR.  However the model contains partial differential equations which are currently beyond the scope of the tools.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>
<para>
Traditionally the conceptual framework of cellular energy metabolism has been based on the assumption that the reactions involved are at equilibrium. However, this assumption is based on little experimental evidence, and the experiments that have been carried out have all been in resting muscle - where the equilibrium assumption is more likely to be met. More recently this assumption of equilibrium has been analysed using mathematical models of compartmentalised energy transfer, and these studies have shown that the assumption is not totally valid. In a pumping heart, the creatine kinase (CK) reaction is out of equilibrium for most of the contraction cycle.
</para>
<para>
This observation led to two questions:
1) For which reactions, if any, is the assumption of equilibrium reasonable? and
2) Do concentration gradients exist for metabolites, and how steep are these gradients?
Understanding these two points is important to the understanding of the intracellular factors which regulate the rate of mitochondrial oxidative phosphorylation in the heart cells in vivo.</para>

<para>
In the 2000 Vendelin et al. publication described here, the authors develop a mathematical model that describes quantitatively the experimental data on how the rate of oxygen consumption and the metabolic state of the working isolated perfused rat heart depend on workload over its physiological range (see Figure 1 below). They use the model of Aliev and Saks (1997) as a foundation, and then supplement this with the model of oxidative phosphorylation by Korzeniewski and Froncisz (1991, 1998), in order to account for the changes in mitochondrial membrane potential and to calculate the rates of oxygen consumption for different workloads.
</para>

<para>
Model simulation results suggest that the assumption of reaction equilibrium for phosphocreatine, creatine, and ATP diffusion in the intracellular bulk water phase is acceptable. However, due to the nonequilibrium state of the CK reaction, intracellular ADP concentration gradients exist. Together with changes in inorganic phosphate (Pi) concentration, localised changes in ADP concentration (oscillations) are thought to be a major component of the metabolic feedback signal for regulation of respiration in muscle cells.
</para>

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

<para>
<ulink url="http://ajpcell.physiology.org/cgi/content/abstract/278/4/C747">Regulation of mitochondrial respiration in heart cells analyzed by reaction-diffusion model of energy transfer</ulink>, Marko Vendelin, Olav Kongas, and Valdur Saks, 2000, <ulink url="http://ajpcell.physiology.org/">
            <emphasis>American Journal of Physiology - Cell Physiology</emphasis>, 278, C747-C764.
          </ulink>  (<ulink url="http://ajpcell.physiology.org/cgi/content/full/278/4/C747">Full text (HTML)</ulink> and <ulink url="http://ajpcell.physiology.org/cgi/reprint/278/4/C747">PDF</ulink> versions of the article are available on the <emphasis>American Journal of Physiology </emphasis>website.) <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/10751324?dopt=Abstract">PubMed ID: 10751324</ulink> 
</para>

<informalfigure float="0">
<mediaobject>
  <imageobject>
    <imagedata fileref="vendelin_2000.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the two-dimensional compartmentalised energy transfer model of cardiac cells: the components of the model. The total diffusion pathway is 1.2 micrometres, including the outer mitochondrial membrane, the layer of myoplasm and the myofibril. Creatine kinase (CK) and myosin magnesium ATPase activities are distributed non-uniformly along the myofibril.</caption>
</informalfigure>

</sect1>
</article>
</documentation>          
  
  <units name="millisecond">
    <unit units="second" prefix="milli"/>
  </units>
  
  <units name="molar">
    <unit units="mole"/>
    <unit units="litre" exponent="-1"/>
  </units>
   
  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="micromolar2_per_s">
    <unit units="micromolar" exponent="2"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="micrometre">
    <unit units="metre" prefix="micro"/>
  </units>
  
  <units name="micrometre_per_s">
    <unit units="micrometre"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>
  
  <units name="per_millivolt">
    <unit units="millivolt" exponent="-1"/>
  </units>
  
  <units name="flux">
    <unit units="micromolar"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="second_order_rate_constant">
    <unit units="micromolar" exponent="-1"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="kilojoule_per_mole">
    <unit units="joule" prefix="kilo"/>
    <unit units="mole" exponent="-1"/>
  </units>
  
  <units name="kilojoule_per_mole_kelvin">
    <unit units="joule" prefix="kilo"/>
    <unit units="mole" exponent="-1"/>
    <unit units="kelvin" exponent="-1"/>
  </units>
  
  <units name="kilojoule_per_mole_millivolt">
    <unit units="joule" prefix="kilo"/>
    <unit units="mole" exponent="-1"/>
    <unit units="millivolt" exponent="-1"/>
  </units>
  
  
  
  <component name="environment">
    <variable units="millisecond" public_interface="out" name="time"/>
  </component>
  
  <component cmeta:id="ATP" name="ATP">
    <rdf:RDF>
      <rdf:Description rdf:about="ATP">
        <dc:title>ATP</dc:title>
        <dcterms:alternative>adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="ATP"/>
    
    <variable units="dimensionless" public_interface="in" name="G_CK"/>
    <variable units="flux" public_interface="in" name="v_CK"/>
    <variable units="dimensionless" public_interface="in" name="G_AK"/>
    <variable units="flux" public_interface="in" name="v_AK"/>
    <variable units="dimensionless" public_interface="in" name="G_H"/>
    <variable units="flux" public_interface="in" name="H_ATP"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_ATP"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>ATP</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> D_ATP </ci>
                <ci> ATP </ci>
              </apply>
              <apply>
                <times/>
                <ci> G_CK </ci>
                <ci> v_CK </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> G_AK </ci>
              <ci> v_AK </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> G_H </ci>
            <ci> H_ATP </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="ADP" name="ADP">
    <rdf:RDF>
      <rdf:Description rdf:about="ADP">
        <dc:title>ADP</dc:title>
        <dcterms:alternative>adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="ADP"/>
    
    <variable units="dimensionless" public_interface="in" name="G_CK"/>
    <variable units="flux" public_interface="in" name="v_CK"/>
    <variable units="dimensionless" public_interface="in" name="G_AK"/>
    <variable units="flux" public_interface="in" name="v_AK"/>
    <variable units="dimensionless" public_interface="in" name="G_H"/>
    <variable units="flux" public_interface="in" name="H_ATP"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_ADP"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>ADP</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> D_ADP </ci>
                <ci> ADP </ci>
              </apply>
              <apply>
                <times/>
                <ci> G_CK </ci>
                <ci> v_CK </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> G_AK </ci>
              <ci> v_AK </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> G_H </ci>
            <ci> H_ATP </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="AMP" name="AMP">
    <rdf:RDF>
      <rdf:Description rdf:about="AMP">
        <dc:title>AMP</dc:title>
        <dcterms:alternative>adenosine monophosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="AMP"/>
    
    <variable units="dimensionless" public_interface="in" name="G_AK"/>
    <variable units="flux" public_interface="in" name="v_AK"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_AMP"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>AMP</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> D_AMP </ci>
            <ci> AMP </ci>
          </apply>
          <apply>
            <times/>
            <ci> G_AK </ci>
            <ci> v_AK </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="PCr" name="PCr">
    <rdf:RDF>
      <rdf:Description rdf:about="PCr">
        <dc:title>PCr</dc:title>
        <dcterms:alternative>phosphocreatine</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="PCr"/>
    
    <variable units="dimensionless" public_interface="in" name="G_CK"/>
    <variable units="flux" public_interface="in" name="v_CK"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_PCr"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>PCr</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> D_PCr </ci>
            <ci> PCr </ci>
          </apply>
          <apply>
            <times/>
            <ci> G_CK </ci>
            <ci> v_CK </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="Cr" name="Cr">
    <rdf:RDF>
      <rdf:Description rdf:about="Cr">
        <dc:title>Cr</dc:title>
        <dcterms:alternative>creatine</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="Cr"/>
    
    <variable units="dimensionless" public_interface="in" name="G_CK"/>
    <variable units="flux" public_interface="in" name="v_CK"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_Cr"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Cr</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> D_Cr </ci>
            <ci> Cr </ci>
          </apply>
          <apply>
            <times/>
            <ci> G_CK </ci>
            <ci> v_CK </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="Pi" name="Pi">
    <rdf:RDF>
      <rdf:Description rdf:about="Pi">
        <dc:title>Pi</dc:title>
        <dcterms:alternative>inorganic phosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="Pi"/>
    
    <variable units="dimensionless" public_interface="in" name="G_H"/>
    <variable units="flux" public_interface="in" name="H_ATP"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_Pi"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Pi</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> D_Pi </ci>
            <ci> Pi </ci>
          </apply>
          <apply>
            <times/>
            <ci> G_H </ci>
            <ci> H_ATP </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="v_CK"> 
    <variable units="flux" public_interface="out" name="v_CK"/>
    
    <variable units="micromolar" name="Kia" initial_value="9.0E2"/>
    <variable units="micromolar" name="Kib" initial_value="3.49E4"/>
    <variable units="micromolar" name="Kc"/>
    <variable units="micromolar" name="Kb" initial_value="1.55E4"/>
    <variable units="micromolar" name="Kic" initial_value="2.224E2"/>
    <variable units="micromolar" name="KIb" initial_value="3.49E4"/>
    <variable units="micromolar" name="Kid" initial_value="4.73E3"/>
    <variable units="micromolar" name="Kd" initial_value="1.67E3"/>
    <variable units="flux" name="V1" initial_value="1.30E4"/>
    <variable units="flux" name="V1_" initial_value="5.48E4"/>
    <variable units="dimensionless" name="DenCK"/>
   
    <variable units="micromolar" public_interface="in" name="mATP"/>
    <variable units="micromolar" public_interface="in" name="Cr"/>
    <variable units="micromolar" public_interface="in" name="PCr"/>
    <variable units="micromolar" public_interface="in" name="mADP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/>
        <ci> v_CK </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> V1 </ci>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> mATP </ci>
                  <ci> Cr </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> Kia </ci>
                  <ci> Kb </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> V1_ </ci>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> mADP </ci>
                  <ci> PCr </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> Kic </ci>
                  <ci> Kd </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <ci> DenCK </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <ci> DenCK </ci>
        <apply>
          <plus/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <divide/>
            <ci> Cr </ci>
            <ci> Kib </ci>
          </apply>
          <apply>
            <divide/>
            <ci> PCr </ci>
            <ci> Kid </ci>
          </apply>
          <apply>
            <times/>
            <ci> mATP </ci>
            <apply>
              <plus/>
              <apply>
                <divide/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> Kia </ci>
              </apply>
              <apply>
                <divide/>
                <ci> Cr </ci>
                <apply>
                  <times/>
                  <ci> Kia </ci>
                  <ci> Kb </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> mADP </ci>
            <apply>
              <plus/>
              <apply>
                <divide/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> Kic </ci>
              </apply>
              <apply>
                <divide/>
                <ci> PCr </ci>
                <apply>
                  <times/>
                  <ci> Kid </ci>
                  <ci> Kc </ci>
                </apply>
              </apply>
              <apply>
                <divide/>
                <ci> Cr </ci>
                <apply>
                  <times/>
                  <ci> Kic </ci>
                  <ci> KIb </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <ci> Kc </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> Kic </ci>
            <ci> Kd </ci>
          </apply>
          <ci> Kid </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component cmeta:id="fATP" name="fATP">
    <rdf:RDF>
      <rdf:Description rdf:about="fATP">
        <dc:title>fATP</dc:title>
        <dcterms:alternative>total Mg-free adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fATP"/>
    
    <variable units="micromolar" name="KDT"/>
    
    <variable units="micromolar" public_interface="in" name="ATP"/>
    <variable units="micromolar" public_interface="in" name="Mg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fATP</ci> 
        <apply>
          <divide/>
          <ci> ATP </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mg </ci>
              <ci> KDT </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="fATPx" name="fATPx">
    <rdf:RDF>
      <rdf:Description rdf:about="fATPx">
        <dc:title>fATPx</dc:title>
        <dcterms:alternative>mitochondrial matrix Mg-free adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fATPx"/>
    
    <variable units="micromolar" name="KDTx" initial_value="17.0"/>
    
    <variable units="micromolar" public_interface="in" name="ATPx"/>
    <variable units="micromolar" public_interface="in" name="Mgx"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fATPx</ci> 
        <apply>
          <divide/>
          <ci> ATPx </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mgx </ci>
              <ci> KDTx </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="fATPext" name="fATPext">
    <rdf:RDF>
      <rdf:Description rdf:about="fATPext">
        <dc:title>fATPext</dc:title>
        <dcterms:alternative>Cytoplasmic Mg-free adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fATPext"/>
    
    <variable units="micromolar" name="KDText" initial_value="24.0"/>
    
    <variable units="micromolar" public_interface="in" name="ATPext"/>
    <variable units="micromolar" public_interface="in" name="Mgext"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fATPext</ci> 
        <apply>
          <divide/>
          <ci> ATPext </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mgext </ci>
              <ci> KDText </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="fATPg" name="fATPg">
    <rdf:RDF>
      <rdf:Description rdf:about="fATPg">
        <dc:title>fATPg</dc:title>
        <dcterms:alternative>microcompartment Mg-free adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fATPg"/>
    
    <variable units="micromolar" name="KDTg"/>
    
    <variable units="micromolar" public_interface="in" name="ATPg"/>
    <variable units="micromolar" public_interface="in" name="Mgg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fATPg</ci> 
        <apply>
          <divide/>
          <ci> ATPg </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mgg </ci>
              <ci> KDTg </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="fATPi" name="fATPi">
    <rdf:RDF>
      <rdf:Description rdf:about="fATPi">
        <dc:title>fATPi</dc:title>
        <dcterms:alternative>mitochondrial IM space Mg-free adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fATPi"/>
    
    <variable units="micromolar" name="KDTi"/>
    
    <variable units="micromolar" public_interface="in" name="ATPi"/>
    <variable units="micromolar" public_interface="in" name="Mgi"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fATPi</ci> 
        <apply>
          <divide/>
          <ci> ATPi </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mgi </ci>
              <ci> KDTi </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mATP" name="mATP">
    <rdf:RDF>
      <rdf:Description rdf:about="mATP">
        <dc:title>mATP</dc:title>
        <dcterms:alternative>total Mg-bound adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="mATP"/>
    
    <variable units="micromolar" public_interface="in" name="ATP"/>
    <variable units="micromolar" public_interface="in" name="fATP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mATP</ci> 
        <apply>
          <minus/>
          <ci> ATP </ci>
          <ci> fATP </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mATPx" name="mATPx">
    <rdf:RDF>
      <rdf:Description rdf:about="mATPx">
        <dc:title>mATPx</dc:title>
        <dcterms:alternative>mitochondrial matrix Mg-bound adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" name="mATPx"/>
    
    <variable units="micromolar" public_interface="in" name="ATPx"/>
    <variable units="micromolar" public_interface="in" name="fATPx"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mATPx</ci> 
        <apply>
          <minus/>
          <ci> ATPx </ci>
          <ci> fATPx </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mATPext" name="mATPext">
    <rdf:RDF>
      <rdf:Description rdf:about="mATPext">
        <dc:title>mATPext</dc:title>
        <dcterms:alternative>myoplasmic Mg-bound adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" name="mATPext"/>
    
    <variable units="micromolar" public_interface="in" name="ATPext"/>
    <variable units="micromolar" public_interface="in" name="fATPext"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mATPext</ci> 
        <apply>
          <minus/>
          <ci> ATPext </ci>
          <ci> fATPext </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mATPg" name="mATPg">
    <rdf:RDF>
      <rdf:Description rdf:about="mATPg">
        <dc:title>mATPg</dc:title>
        <dcterms:alternative>microcompartment Mg-bound adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="mATPg"/>
    
    <variable units="micromolar" public_interface="in" name="ATPg"/>
    <variable units="micromolar" public_interface="in" name="fATPg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mATPg</ci> 
        <apply>
          <minus/>
          <ci> ATPg </ci>
          <ci> fATPg </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mATPi" name="mATPi">
    <rdf:RDF>
      <rdf:Description rdf:about="mATPi">
        <dc:title>mATPi</dc:title>
        <dcterms:alternative>mitochondrial IM space Mg-bound adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="mATPi"/>
    
    <variable units="micromolar" public_interface="in" name="ATPi"/>
    <variable units="micromolar" public_interface="in" name="fATPi"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mATPi</ci> 
        <apply>
          <minus/>
          <ci> ATPi </ci>
          <ci> fATPi </ci>
        </apply>
      </apply>    
    </math>
  </component>
 
  <component cmeta:id="fADP" name="fADP">
    <rdf:RDF>
      <rdf:Description rdf:about="fADP">
        <dc:title>fADP</dc:title>
        <dcterms:alternative>total Mg-free adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fADP"/>
    
    <variable units="micromolar" name="KDD"/>
    
    <variable units="micromolar" public_interface="in" name="ADP"/>
    <variable units="micromolar" public_interface="in" name="Mg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fADP</ci> 
        <apply>
          <divide/>
          <ci> ADP </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mg </ci>
              <ci> KDD </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="fADPx" name="fADPx">
    <rdf:RDF>
      <rdf:Description rdf:about="fADPx">
        <dc:title>fADPx</dc:title>
        <dcterms:alternative>mitochondrial matrix Mg-free adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fADPx"/>
    
    <variable units="micromolar" name="KDDx" initial_value="282.0"/>
    
    <variable units="micromolar" public_interface="in" name="ADPx"/>
    <variable units="micromolar" public_interface="in" name="Mgx"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fADPx</ci> 
        <apply>
          <divide/>
          <ci> ADPx </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mgx </ci>
              <ci> KDDx </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="fADPext" name="fADPext">
    <rdf:RDF>
      <rdf:Description rdf:about="fADPext">
        <dc:title>fADPext</dc:title>
        <dcterms:alternative>myoplasmic Mg-free adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fADPext"/>
    
    <variable units="micromolar" name="KDDext" initial_value="347.0"/>
    
    <variable units="micromolar" public_interface="in" name="ADPext"/>
    <variable units="micromolar" public_interface="in" name="Mgext"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fADPext</ci> 
        <apply>
          <divide/>
          <ci> ADPext </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mgext </ci>
              <ci> KDDext </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="fADPg" name="fADPg">
    <rdf:RDF>
      <rdf:Description rdf:about="fADPg">
        <dc:title>fADPg</dc:title>
        <dcterms:alternative>microspace Mg-free adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fADPg"/>
    
    <variable units="micromolar" name="KDDg"/>
    
    <variable units="micromolar" public_interface="in" name="ADPg"/>
    <variable units="micromolar" public_interface="in" name="Mgg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fADPg</ci> 
        <apply>
          <divide/>
          <ci> ADPg </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mgg </ci>
              <ci> KDDg </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="fADPi" name="fADPi">
    <rdf:RDF>
      <rdf:Description rdf:about="fADPi">
        <dc:title>fADPi</dc:title>
        <dcterms:alternative>mitochondrial IM space Mg-free adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="fADPi"/>
    
    <variable units="micromolar" name="KDDi"/>
    
    <variable units="micromolar" public_interface="in" name="ADPi"/>
    <variable units="micromolar" public_interface="in" name="Mgi"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>fADPi</ci> 
        <apply>
          <divide/>
          <ci> ADPi </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <divide/>
              <ci> Mgi </ci>
              <ci> KDDi </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mADP" name="mADP">
    <rdf:RDF>
      <rdf:Description rdf:about="mADP">
        <dc:title>mADP</dc:title>
        <dcterms:alternative>total Mg-bound adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="mADP"/>
    
    <variable units="micromolar" public_interface="in" name="ADP"/>
    <variable units="micromolar" public_interface="in" name="fADP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mADP</ci> 
        <apply>
          <minus/>
          <ci> ADP </ci>
          <ci> fADP </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mADPx" name="mADPx">
    <rdf:RDF>
      <rdf:Description rdf:about="mADPx">
        <dc:title>mADPx</dc:title>
        <dcterms:alternative>mitochondrial matrix Mg-bound adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" name="mADPx"/>
    
    <variable units="micromolar" public_interface="in" name="ADPx"/>
    <variable units="micromolar" public_interface="in" name="fADPx"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mADPx</ci> 
        <apply>
          <minus/>
          <ci> ADPx </ci>
          <ci> fADPx </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mADPext" name="mADPext">
    <rdf:RDF>
      <rdf:Description rdf:about="mADPext">
        <dc:title>mADPext</dc:title>
        <dcterms:alternative>myoplasmic Mg-bound adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" name="mADPext"/>
    
    <variable units="micromolar" public_interface="in" name="ADPext"/>
    <variable units="micromolar" public_interface="in" name="fADPext"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mADPext</ci> 
        <apply>
          <minus/>
          <ci> ADPext </ci>
          <ci> fADPext </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mADPg" name="mADPg">
    <rdf:RDF>
      <rdf:Description rdf:about="mADPg">
        <dc:title>mADPg</dc:title>
        <dcterms:alternative>microspace Mg-bound adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="mADPg"/>
    
    <variable units="micromolar" public_interface="in" name="ADPg"/>
    <variable units="micromolar" public_interface="in" name="fADPg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mADPg</ci> 
        <apply>
          <minus/>
          <ci> ADPg </ci>
          <ci> fADPg </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="mADPi" name="mADPi">
    <rdf:RDF>
      <rdf:Description rdf:about="mADPi">
        <dc:title>mADPi</dc:title>
        <dcterms:alternative>mitochondrial IM space Mg-bound adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="mADPi"/>
    
    <variable units="micromolar" public_interface="in" name="ADPi"/>
    <variable units="micromolar" public_interface="in" name="fADPi"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>mADPi</ci> 
        <apply>
          <minus/>
          <ci> ADPi </ci>
          <ci> fADPi </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="v_AK"> 
    <variable units="flux" public_interface="out" name="v_AK"/>
    
    <variable units="second_order_rate_constant" name="kfa" initial_value="63.8"/>
    <variable units="second_order_rate_constant" name="kba" initial_value="1.68"/>
    
    <variable units="micromolar" public_interface="in" name="fADP"/>
    <variable units="micromolar" public_interface="in" name="mATP"/>
    <variable units="micromolar" public_interface="in" name="AMP"/>
    <variable units="micromolar" public_interface="in" name="mADP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/>
        <ci> v_AK </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kfa </ci>
            <ci> fADP </ci>
            <ci> mADP </ci>
          </apply>
          <apply>
            <times/>
            <ci> kba </ci>
            <ci> mATP </ci>
            <ci> AMP </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="v_AKmit"> 
    <variable units="flux" public_interface="out" name="v_AKmit"/>
    
    <variable units="second_order_rate_constant" name="kfa" initial_value="63.8"/>
    <variable units="second_order_rate_constant" name="kba" initial_value="1.68"/>
    
    <variable units="micromolar" public_interface="in" name="fADP"/>
    <variable units="micromolar" public_interface="in" name="mATP"/>
    <variable units="micromolar" public_interface="in" name="AMP"/>
    <variable units="micromolar" public_interface="in" name="mADP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/>
        <ci> v_AKmit </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> kfa </ci>
            <ci> fADP </ci>
            <ci> mADP </ci>
          </apply>
          <apply>
            <times/>
            <ci> kba </ci>
            <ci> mATP </ci>
            <ci> AMP </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="H_ATP" name="H_ATP">
    <rdf:RDF>
      <rdf:Description rdf:about="H_ATP">
        <dc:title>H_ATP</dc:title>
        <dcterms:alternative>ATP hydrolysis function</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="flux" public_interface="out" name="H_ATP"/>
    
    <variable units="flux" name="H_ATPmax" initial_value="3.626E4"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>H_ATP</ci>
        <piecewise>
          <piece>
            <apply>
              <times/> 
              <apply>
                <divide/>
                <ci> time </ci>
                <cn cellml:units="millisecond"> 30.0 </cn>
              </apply>
              <ci> H_ATPmax </ci>
            </apply>
            <apply>
              <and/>
              <apply>
                <geq/>
                <ci> time </ci> 
                <cn cellml:units="millisecond"> 0.0 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> time </ci> 
                <cn cellml:units="millisecond"> 30.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <times/> 
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="millisecond"> 60.0 </cn>
                  <ci> time </ci>
                </apply>
                <cn cellml:units="millisecond"> 30.0 </cn>
              </apply>
              <ci> H_ATPmax </ci>
            </apply>
            <apply>
              <and/>
              <apply>
                <geq/>
                <ci> time </ci> 
                <cn cellml:units="millisecond"> 30.0 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> time </ci> 
                <cn cellml:units="millisecond"> 60.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="flux"> 0.0 </cn>
            <apply>
              <and/>
              <apply>
                <geq/>
                <ci> time </ci> 
                <cn cellml:units="millisecond"> 60.0 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> time </ci> 
                <cn cellml:units="millisecond"> 180.0 </cn>
              </apply>
            </apply>
          </piece>
        </piecewise>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="G_CK" name="G_CK">
    <rdf:RDF>
      <rdf:Description rdf:about="G_CK">
        <dc:title>G_CK</dc:title>
        <dcterms:alternative>spatial distribution of CK</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="dimensionless" public_interface="out" name="G_CK"/>
    
    <variable units="micrometre" name="x"/>
    <variable units="micrometre" name="y"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>G_CK</ci>
        <piecewise>
          <piece>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.94 </cn> 
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <divide/> 
                  <ci> x </ci>
                  <cn cellml:units="micrometre"> 0.375 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.375 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> y </ci> 
                <cn cellml:units="micrometre"> 1.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 1.61 </cn>
              <apply>
                <divide/> 
                <apply>
                  <minus/>
                  <ci> x </ci>
                  <cn cellml:units="micrometre"> 0.375 </cn>
                </apply>
                <cn cellml:units="micrometre"> 0.375 </cn>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <geq/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.375 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.75 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> y </ci> 
                <cn cellml:units="micrometre"> 1.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="dimensionless"> 1.61 </cn>
            <apply>
              <or/>
              <apply>
                <geq/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.375 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> y </ci> 
                <cn cellml:units="micrometre"> 1.0 </cn>
              </apply>
            </apply>
          </piece>
        </piecewise>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="G_AK" name="G_AK">
    <rdf:RDF>
      <rdf:Description rdf:about="G_AK">
        <dc:title>G_AK</dc:title>
        <dcterms:alternative>spatial distribution of AK</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="dimensionless" public_interface="out" name="G_AK"/>
    
    <variable units="dimensionless" public_interface="in" name="G_CK"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>G_AK</ci>
        <ci>G_CK</ci>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="G_H" name="G_H">
    <rdf:RDF>
      <rdf:Description rdf:about="G_H">
        <dc:title>G_H</dc:title>
        <dcterms:alternative>spatial distribution of ATPase</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="dimensionless" public_interface="out" name="G_H"/>
    
    <variable units="micrometre" name="x"/>
    <variable units="micrometre" name="y"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>G_H</ci>
        <piecewise>
          <piece>
            <apply>
              <times/>
              <apply>
                <divide/> 
                <cn cellml:units="dimensionless"> 4.0 </cn> 
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <apply>
                <divide/> 
                <ci> x </ci>
                <cn cellml:units="micrometre"> 0.125 </cn>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <lt/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.125 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> y </ci> 
                <cn cellml:units="micrometre"> 1.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <divide/> 
              <cn cellml:units="dimensionless"> 4.0 </cn> 
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
            <apply>
              <and/>
              <apply>
                <geq/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.125 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.75 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> y </ci> 
                <cn cellml:units="micrometre"> 1.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <apply>
              <times/>
              <apply>
                <divide/> 
                <cn cellml:units="dimensionless"> 4.0 </cn> 
                <cn cellml:units="dimensionless"> 3.0 </cn>
              </apply>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <divide/>
                  <apply>
                    <minus/> 
                    <ci> x </ci>
                    <cn cellml:units="micrometre"> 0.75 </cn>
                  </apply>
                  <cn cellml:units="micrometre"> 0.125 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <and/>
              <apply>
                <geq/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.75 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.875 </cn>
              </apply>
              <apply>
                <lt/>
                <ci> y </ci> 
                <cn cellml:units="micrometre"> 1.0 </cn>
              </apply>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="dimensionless"> 0.0 </cn>
            <apply>
              <or/>
              <apply>
                <geq/>
                <ci> x </ci> 
                <cn cellml:units="micrometre"> 0.875 </cn>
              </apply>
              <apply>
                <geq/>
                <ci> y </ci> 
                <cn cellml:units="micrometre"> 1.0 </cn>
              </apply>
            </apply>
          </piece>
        </piecewise>
      </apply>    
    </math>
  </component>
  
  <component name="diffusion_through_mitochondrial_membrane">
    <variable units="flux" name="ATP_diff"/>
    <variable units="flux" name="ADP_diff"/>
    <variable units="flux" name="AMP_diff"/>
    <variable units="flux" name="PCr_diff"/>
    <variable units="flux" name="Cr_diff"/>
    <variable units="flux" name="Pi_diff"/> 
    
    <variable units="micromolar" public_interface="in" name="ATP"/>
    <variable units="micromolar" public_interface="in" name="ADP"/>
    <variable units="micromolar" public_interface="in" name="AMP"/>
    <variable units="micromolar" public_interface="in" name="PCr"/>
    <variable units="micromolar" public_interface="in" name="Cr"/>
    <variable units="micromolar" public_interface="in" name="Pi"/>
    <variable units="micromolar" public_interface="in" name="ATPi"/>
    <variable units="micromolar" public_interface="in" name="ADPi"/>
    <variable units="micromolar" public_interface="in" name="AMPi"/>
    <variable units="micromolar" public_interface="in" name="PCri"/>
    <variable units="micromolar" public_interface="in" name="Cri"/>
    <variable units="micromolar" public_interface="in" name="Pii"/>
    <variable units="micrometre_per_s" public_interface="in" name="R_ATP"/>
    <variable units="micrometre_per_s" public_interface="in" name="R_ADP"/>
    <variable units="micrometre_per_s" public_interface="in" name="R_AMP"/>
    <variable units="micrometre_per_s" public_interface="in" name="R_PCr"/>
    <variable units="micrometre_per_s" public_interface="in" name="R_Cr"/>
    <variable units="micrometre_per_s" public_interface="in" name="R_Pi"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_ATP"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_ADP"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_AMP"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_PCr"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_Cr"/>
    <variable units="micromolar2_per_s" public_interface="in" name="D_Pi"/>
    
    <variable units="dimensionless" name="n"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>n</ci>
          </bvar> 
          <ci>ATP_diff</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> R_ATP </ci>
            <apply>
              <minus/>
              <ci> ATPi </ci>
              <ci> ATP </ci>
            </apply>
          </apply>
          <ci> D_ATP </ci>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>n</ci>
          </bvar> 
          <ci>ADP_diff</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> R_ADP </ci>
            <apply>
              <minus/>
              <ci> ADPi </ci>
              <ci> ADP </ci>
            </apply>
          </apply>
          <ci> D_ADP </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>n</ci>
          </bvar> 
          <ci>AMP_diff</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> R_AMP </ci>
            <apply>
              <minus/>
              <ci> AMPi </ci>
              <ci> AMP </ci>
            </apply>
          </apply>
          <ci> D_AMP </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>n</ci>
          </bvar> 
          <ci>PCr_diff</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> R_PCr </ci>
            <apply>
              <minus/>
              <ci> PCri </ci>
              <ci> PCr </ci>
            </apply>
          </apply>
          <ci> D_PCr </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>n</ci>
          </bvar> 
          <ci>Cr_diff</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> R_Cr </ci>
            <apply>
              <minus/>
              <ci> Cri </ci>
              <ci> Cr </ci>
            </apply>
          </apply>
          <ci> D_Cr </ci>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <apply>
          <partialdiff/> 
          <bvar>
            <ci>n</ci>
          </bvar> 
          <ci>Pi_diff</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> R_Pi </ci>
            <apply>
              <minus/>
              <ci> Pii </ci>
              <ci> Pi </ci>
            </apply>
          </apply>
          <ci> D_Pi </ci>
        </apply>
      </apply>      
    </math>
  </component>
  
  <component cmeta:id="ATPi" name="ATPi">
    <rdf:RDF>
      <rdf:Description rdf:about="ATPi">
        <dc:title>ATPi</dc:title>
        <dcterms:alternative>IM space adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="ATPi"/>
    
    <variable units="micrometre_per_s" public_interface="in" name="R_ATP"/>
    <variable units="dimensionless" public_interface="in" name="FVi"/>
    <variable units="flux" public_interface="in" name="v_ANT"/>
    <variable units="flux" public_interface="in" name="v_MiCK"/>
    <variable units="flux" public_interface="in" name="v_AKmit"/>
    <variable units="micromolar" public_interface="in" name="ATP"/>
    <variable units="micrometre" public_interface="in" name="Li"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>ATPi</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <apply>
                  <divide/>
                  <ci> R_ATP </ci>
                  <ci> Li </ci>
                </apply>
                <apply>
                  <minus/>
                  <ci> ATP </ci>
                  <ci> ATPi </ci>
                </apply>
              </apply>
              <apply>
                <divide/>
                <ci> v_ANT </ci>
                <ci> FVi </ci>
              </apply>
            </apply>
            <apply>
              <divide/>
              <ci> v_MiCK </ci>
              <ci> FVi </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> v_AKmit </ci>
            <ci> FVi </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="ADPi" name="ADPi">
    <rdf:RDF>
      <rdf:Description rdf:about="ADPi">
        <dc:title>ADPi</dc:title>
        <dcterms:alternative>IM space adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="ADPi"/>
    
    <variable units="micrometre_per_s" public_interface="in" name="R_ADP"/>
    <variable units="dimensionless" public_interface="in" name="FVi"/>
    <variable units="flux" public_interface="in" name="v_ANT"/>
    <variable units="flux" public_interface="in" name="v_MiCK"/>
    <variable units="flux" public_interface="in" name="v_AKmit"/>
    <variable units="micromolar" public_interface="in" name="ADP"/>
    <variable units="micrometre" public_interface="in" name="Li"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>ADPi</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <apply>
                <divide/>
                <ci> R_ADP </ci>
                <ci> Li </ci>
              </apply>
              <apply>
                <minus/>
                <ci> ADP </ci>
                <ci> ADPi </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> v_AKmit </ci>
                </apply>
                <ci> FVi </ci>
              </apply>
              <apply>
                <divide/>
                <ci> v_ANT </ci>
                <ci> FVi </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> v_MiCK </ci>
            <ci> FVi </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="AMPi" name="AMPi">
    <rdf:RDF>
      <rdf:Description rdf:about="AMPi">
        <dc:title>AMPi</dc:title>
        <dcterms:alternative>IM space adenosine monophosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="AMPi"/>
    
    <variable units="micrometre_per_s" public_interface="in" name="R_AMP"/>
    <variable units="dimensionless" public_interface="in" name="FVi"/>
    <variable units="flux" public_interface="in" name="v_AKmit"/>
    <variable units="micromolar" public_interface="in" name="AMP"/>
    <variable units="micrometre" public_interface="in" name="Li"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>AMPi</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> R_AMP </ci>
              <ci> Li </ci>
            </apply>
            <apply>
              <minus/>
              <ci> AMP </ci>
              <ci> AMPi </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> v_AKmit </ci>
            <ci> FVi </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="PCri" name="PCri">
    <rdf:RDF>
      <rdf:Description rdf:about="PCri">
        <dc:title>PCri</dc:title>
        <dcterms:alternative>IM space phosphocreatine</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="PCri"/>
    
    <variable units="micrometre_per_s" public_interface="in" name="R_PCr"/>
    <variable units="dimensionless" public_interface="in" name="FVi"/>
    <variable units="flux" public_interface="in" name="v_AKmit"/>
    <variable units="micromolar" public_interface="in" name="PCr"/>
    <variable units="micrometre" public_interface="in" name="Li"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>PCri</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> R_PCr </ci>
              <ci> Li </ci>
            </apply>
            <apply>
              <minus/>
              <ci> PCr </ci>
              <ci> PCri </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> v_AKmit </ci>
            <ci> FVi </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="Cri" name="Cri">
    <rdf:RDF>
      <rdf:Description rdf:about="Cri">
        <dc:title>Cri</dc:title>
        <dcterms:alternative>IM space creatine</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="Cri"/>
    
    <variable units="micrometre_per_s" public_interface="in" name="R_Cr"/>
    <variable units="dimensionless" public_interface="in" name="FVi"/>
    <variable units="flux" public_interface="in" name="v_MiCK"/>
    <variable units="micromolar" public_interface="in" name="Cr"/>
    <variable units="micrometre" public_interface="in" name="Li"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Cri</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> R_Cr </ci>
              <ci> Li </ci>
            </apply>
            <apply>
              <minus/>
              <ci> Cr </ci>
              <ci> Cri </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> v_MiCK </ci>
            <ci> FVi </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="Pii" name="Pii">
    <rdf:RDF>
      <rdf:Description rdf:about="Pii">
        <dc:title>Pii</dc:title>
        <dcterms:alternative>IM space creatine</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="Pii"/>
    
    <variable units="micrometre_per_s" public_interface="in" name="R_Pi"/>
    <variable units="dimensionless" public_interface="in" name="FVi"/>
    <variable units="flux" public_interface="in" name="v_PI"/>
    <variable units="micromolar" public_interface="in" name="Pi"/>
    <variable units="micrometre" public_interface="in" name="Li"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Pii</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> R_Pi </ci>
              <ci> Li </ci>
            </apply>
            <apply>
              <minus/>
              <ci> Pi </ci>
              <ci> Pii </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> v_PI </ci>
            <ci> FVi </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="v_MiCK">
    <variable units="flux" public_interface="out" name="v_MiCK"/> 
    
    <variable units="flux" name="v_MiCK_G"/>
    <variable units="flux" name="v_MiCK_I"/>
    <variable units="flux" name="V1" initial_value="7.7E3"/>
    <variable units="flux" name="V1_" initial_value="2.97E4"/>
    <variable units="micromolar" name="Kia_G" initial_value="7.12E2"/>
    <variable units="micromolar" name="Kia" initial_value="1.61E2"/>
    <variable units="micromolar" name="Kic_G" initial_value="1.42E2"/>
    <variable units="micromolar" name="Kic" initial_value="1.10E4"/>
    <variable units="micromolar" name="Kd" initial_value="5.0E2"/>
    <variable units="micromolar" name="Kb" initial_value="5.2E3"/>
    <variable units="dimensionless" name="DenMiCK"/>
    <variable units="micromolar" name="Kib" initial_value="2.6E4"/>
    <variable units="micromolar" name="Kid" initial_value="1.6E3"/>
    <variable units="micromolar" name="KIb" initial_value="2.6E4"/>
     
    <variable units="micromolar" public_interface="in" name="PCri"/>
    <variable units="micromolar" public_interface="in" name="Cri"/>
    <variable units="micromolar" public_interface="in" name="mATPg"/>
    <variable units="micromolar" public_interface="in" name="mATPi"/>
    <variable units="micromolar" public_interface="in" name="mADPg"/>
    <variable units="micromolar" public_interface="in" name="mADPi"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci> v_MiCK </ci>
        <apply>
          <plus/>
          <ci> v_MiCK_G </ci>
          <ci> v_MiCK_I </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <ci> v_MiCK_G </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> V1 </ci>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> mATPg </ci>
                  <ci> Cri </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> Kia_G </ci>
                  <ci> Kb </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> V1_ </ci>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> mADPg </ci>
                  <ci> PCri </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> Kic_G </ci>
                  <ci> Kd </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <ci> DenMiCK </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <ci> v_MiCK_I </ci>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> V1 </ci>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> mATPi </ci>
                  <ci> Cri </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> Kia </ci>
                  <ci> Kb </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> V1_ </ci>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> mADPi </ci>
                  <ci> PCri </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> Kic </ci>
                  <ci> Kd </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <ci> DenMiCK </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <ci> DenMiCK </ci>
        <apply>
          <plus/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <divide/>
            <ci> Cri </ci>
            <ci> Kib </ci>
          </apply>
          <apply>
            <divide/>
            <ci> PCri </ci>
            <ci> Kid </ci>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <divide/>
                <ci> mATPg </ci>
                <ci> Kia_G </ci>
              </apply>
              <apply>
                <divide/>
                <ci> mATPi </ci>
                <ci> Kia </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> Cri </ci>
                <ci> Kb </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <divide/>
                <ci> mADPg </ci>
                <ci> Kic_G </ci>
              </apply>
              <apply>
                <divide/>
                <ci> mADPi </ci>
                <ci> Kic </ci>
              </apply>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> PCri </ci>
                <ci> Kd </ci>
              </apply>
              <apply>
                <divide/>
                <ci> Cri </ci>
                <ci> KIb </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>
  
  <component name="v_exch"> 
    <variable units="flux" public_interface="out" name="v_exch_ATP"/>
    <variable units="flux" public_interface="out" name="v_exch_ADP"/>
    
    <variable units="first_order_rate_constant" name="R_exch_ATP" initial_value="1.88"/>
    <variable units="first_order_rate_constant" name="R_exch_ADP" initial_value="1.88"/>
     
    <variable units="micromolar" public_interface="in" name="ATPi"/>
    <variable units="micromolar" public_interface="in" name="ADPi"/>
    <variable units="micromolar" public_interface="in" name="ATPg"/>
    <variable units="micromolar" public_interface="in" name="ADPg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci> v_exch_ATP </ci>
        <apply>
          <times/>
          <ci> R_exch_ATP </ci>
          <apply>
            <minus/>
            <ci> ATPi </ci>
            <ci> ATPg </ci>
          </apply>      
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <ci> v_exch_ADP </ci>
        <apply>
          <times/>
          <ci> R_exch_ADP </ci>
          <apply>
            <minus/>
            <ci> ADPi </ci>
            <ci> ADPg </ci>
          </apply>      
        </apply>
      </apply>
    </math>
  </component>
  
  <component cmeta:id="UQ" name="UQ">
    <rdf:RDF>
      <rdf:Description rdf:about="UQ">
        <dc:title>UQ</dc:title>
        <dcterms:alternative>oxidised coenzyme Q</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="UQ"/>
    
    <variable units="dimensionless" public_interface="in" name="FVx"/>
    <variable units="flux" public_interface="in" name="v_C1"/>
    <variable units="flux" public_interface="in" name="v_C3"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>UQ</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> v_C3 </ci>
            <ci> v_C1 </ci>
          </apply>
          <ci> FVx </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="c3plus" name="c3plus">
    <rdf:RDF>
      <rdf:Description rdf:about="c3plus">
        <dc:title>c3plus</dc:title>
        <dcterms:alternative>oxidised cytochrome c</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="c3plus"/>
    
    <variable units="dimensionless" public_interface="in" name="FVx"/>
    <variable units="flux" public_interface="in" name="v_O2"/>
    <variable units="flux" public_interface="in" name="v_C3"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>c3plus</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <apply>
              <minus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> v_O2 </ci>
              </apply>
              <ci> v_C3 </ci>
            </apply>
          </apply>
          <ci> FVx </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="NAD" name="NAD">
    <rdf:RDF>
      <rdf:Description rdf:about="NAD">
        <dc:title>NAD</dc:title>
        <dcterms:alternative>matrix NAD</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="NAD"/>
    
    <variable units="flux" public_interface="in" name="v_dh"/>
    <variable units="dimensionless" public_interface="in" name="FVx"/>
    <variable units="flux" public_interface="in" name="v_C1"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>NAD</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> v_C1 </ci>
            <ci> v_dh </ci>
          </apply>
          <ci> FVx </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component cmeta:id="Hx" name="Hx">
    <rdf:RDF>
      <rdf:Description rdf:about="Hx">
        <dc:title>Hx</dc:title>
        <dcterms:alternative>matrix H</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="Hx"/>
    <variable units="flux" public_interface="out" name="delta_Hx"/>
    
    <variable units="dimensionless" public_interface="in" name="na"/>
    <variable units="molar" public_interface="in" name="r_buff"/>
    <variable units="dimensionless" public_interface="in" name="FVx"/>
    <variable units="flux" public_interface="in" name="v_O2"/>
    <variable units="flux" public_interface="in" name="v_C1"/>
    <variable units="flux" public_interface="in" name="v_sn"/>
    <variable units="flux" public_interface="in" name="v_C3"/>
    <variable units="flux" public_interface="in" name="v_ANT"/>
    <variable units="flux" public_interface="in" name="v_PI"/>
    <variable units="flux" public_interface="in" name="v_leak"/>
    <variable units="dimensionless" public_interface="in" name="u"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Hx</ci> 
        </apply>
        <ci>delta_Hx</ci>
      </apply> 
        
      <apply>
        <eq/>
        <ci>delta_Hx</ci>  
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <minus/>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <apply>
                    <plus/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> u </ci>
                  </apply>
                  <ci> v_O2 </ci>
                </apply>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 4.0 </cn>
                  <ci> v_C1 </ci>
                </apply>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <cn cellml:units="dimensionless"> 4.0 </cn>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 2.0 </cn>
                      <ci> u </ci>
                    </apply>
                  </apply>
                  <ci> v_C3 </ci>
                </apply>
              </apply>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci> na </ci>
                  <ci> v_sn </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> v_ANT </ci>
                  <ci> u </ci>
                </apply>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <cn cellml:units="dimensionless"> 1.0 </cn>
                    <ci> u </ci>
                  </apply>
                  <ci> v_PI </ci>
                </apply>
                <ci> v_leak </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> FVx </ci>
            <ci> r_buff </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="ATPx" name="ATPx">
    <rdf:RDF>
      <rdf:Description rdf:about="ATPx">
        <dc:title>ATPx</dc:title>
        <dcterms:alternative>matrix ATP</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="ATPx"/>
    
    <variable units="dimensionless" public_interface="in" name="FVx"/>
    <variable units="flux" public_interface="in" name="v_sn"/>
    <variable units="flux" public_interface="in" name="v_ANT"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>ATPx</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> v_sn </ci>
            <ci> v_ANT </ci>
          </apply>
          <ci> FVx </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="Pix" name="Pix">
    <rdf:RDF>
      <rdf:Description rdf:about="Pix">
        <dc:title>Pix</dc:title>
        <dcterms:alternative>matrix Pi</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF> 
    <variable units="micromolar" public_interface="out" name="Pix"/>
    
    <variable units="dimensionless" public_interface="in" name="FVx"/>
    <variable units="flux" public_interface="in" name="v_sn"/>
    <variable units="flux" public_interface="in" name="v_PI"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Pix</ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> v_PI </ci>
            <ci> v_sn </ci>
          </apply>
          <ci> FVx </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="oxidative_phosphorylation"> 
    <variable units="flux" public_interface="out" name="v_dh"/>
    <variable units="flux" public_interface="out" name="v_C1"/>
    <variable units="flux" public_interface="out" name="v_C3"/>
    <variable units="flux" public_interface="out" name="v_O2"/>
    <variable units="flux" public_interface="out" name="v_sn"/>
    <variable units="flux" public_interface="out" name="v_leak"/>
    <variable units="flux" public_interface="out" name="v_PI"/>
    
    <variable units="flux" name="k_L1" initial_value="2.1895E3"/>
    <variable units="per_millivolt" name="k_L2" initial_value="0.003"/>
    <variable units="flux" name="k_dh" initial_value="5.75E4"/>
    <variable units="flux" name="k_sn" initial_value="2.96E4"/>
    <variable units="flux" name="k_C1" initial_value="2.0E2"/>
    <variable units="flux" name="k_C3" initial_value="1.142E2"/>
    <variable units="second_order_rate_constant" name="k_C4" initial_value="2.351"/>
    <variable units="micromolar" name="k_mO" initial_value="150.0"/>
    <variable units="dimensionless" name="k_mN" initial_value="100.0"/>
    <variable units="dimensionless" name="pD" initial_value="0.8"/>
    <variable units="micromolar" name="O2" initial_value="2.4E2"/>
    <variable units="second_order_rate_constant" name="vfPI" initial_value="138.5"/>
    <variable units="second_order_rate_constant" name="vbPI" initial_value="138.5"/>
    <variable units="dimensionless" name="pKa" initial_value="6.8"/>
    
    <variable units="dimensionless" public_interface="in" name="delta_p"/>
    <variable units="dimensionless" public_interface="in" name="delta_Gsn"/>
    <variable units="millivolt" public_interface="in" name="Z"/>
    <variable units="millivolt" public_interface="in" name="En"/>
    <variable units="millivolt" public_interface="in" name="Eu"/>
    <variable units="millivolt" public_interface="in" name="Ec"/>
    <variable units="millivolt" public_interface="in" name="Ea"/>
    <variable units="micromolar" public_interface="in" name="NADH"/>
    <variable units="micromolar" public_interface="in" name="NAD"/>
    <variable units="micromolar" public_interface="in" name="c2plus"/>
    <variable units="micromolar" public_interface="in" name="a2plus"/>
    <variable units="micromolar" public_interface="in" name="Pii"/>
    <variable units="micromolar" public_interface="in" name="Pix"/>
    <variable units="micromolar" public_interface="in" name="Hext"/>
    <variable units="micromolar" public_interface="in" name="Hx"/>
    <variable units="dimensionless" public_interface="in" name="pHx"/>
    <variable units="dimensionless" public_interface="in" name="pHext"/>
    <variable units="dimensionless" public_interface="in" name="u"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/>
        <ci> v_dh </ci>
        <apply>
          <divide/>
          <ci> k_dh </ci>
          <apply>
            <power/>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> k_mN </ci>
                <apply>
                  <divide/>
                  <ci>NAD</ci>
                  <ci>NADH</ci>
                </apply>
              </apply> 
            </apply>
            <ci> pD </ci>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <ci> v_C1 </ci>
        <apply>
          <times/>
          <ci> k_C1 </ci>
          <apply>
            <minus/>
            <ci> Eu </ci>
            <apply>
              <plus/>
              <ci> En </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> delta_p </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>   
      
      <apply>
        <eq/>
        <ci> v_C3 </ci>
        <apply>
          <times/>
          <ci> k_C3 </ci>
          <apply>
            <minus/>
            <ci> Ec </ci>
            <apply>
              <plus/>
              <ci> Eu </ci>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> u </ci>
                </apply>
                <ci> delta_p </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>  
      
      <apply>
        <eq/>
        <ci> v_O2 </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.5 </cn>
          <ci> k_C4 </ci>
          <ci> a2plus </ci>
          <ci> c2plus </ci>
          <apply>
            <divide/>
            <ci> O2 </ci>
            <apply>
              <plus/>
              <ci> O2 </ci>
              <ci> k_mO </ci>
            </apply>
          </apply>
        </apply>
      </apply>  
      
      <apply>
        <eq/>
        <ci> v_sn </ci>
        <apply>
          <times/>
          <ci> k_sn </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <power/>
                <cn cellml:units="dimensionless"> 10.0 </cn>
                <apply>
                  <divide/>
                  <ci> delta_Gsn </ci>
                  <ci> Z </ci>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
            <apply>
              <plus/>
              <apply>
                <power/>
                <cn cellml:units="dimensionless"> 10.0 </cn>
                <apply>
                  <divide/>
                  <ci> delta_Gsn </ci>
                  <ci> Z </ci>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply>
        <eq/>
        <ci> v_PI </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> vfPI </ci>
            <ci> Hext </ci>
            <apply>
              <divide/>
              <ci> Pii </ci>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 10.0 </cn>
                  <apply>
                    <minus/>
                    <ci> pHext </ci>
                    <ci> pKa </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> vbPI </ci>
            <ci> Hx </ci>
            <apply>
              <divide/>
              <ci> Pix </ci>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <power/>
                  <cn cellml:units="dimensionless"> 10.0 </cn>
                  <apply>
                    <minus/>
                    <ci> pHx </ci>
                    <ci> pKa </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>     
    
      <apply>
        <eq/>
        <ci> v_leak </ci>
        <apply>
          <times/>
          <ci> k_L1 </ci>
          <apply>
            <minus/>
            <apply>
              <power/>
              <cn cellml:units="dimensionless"> 10.0 </cn>
              <apply>
                <times/>
                <ci> k_L2 </ci>
                <ci> delta_p </ci>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>
    </math>     
  </component>
  
  <component name="NADH"> 
    <variable units="micromolar" public_interface="out" name="NADH"/>
    <variable units="micromolar" name="NAD_tot" initial_value="2.97E3"/>
    <variable units="micromolar" public_interface="in" name="NAD"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>NADH</ci> 
        <apply>
          <minus/>
          <ci> NAD_tot </ci>
          <ci> NAD </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="UQH2">
    <variable units="micromolar" public_interface="out" name="UQH2"/>
    <variable units="micromolar" name="UQ_tot" initial_value="1.35E3"/>
    <variable units="micromolar" public_interface="in" name="UQ"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>UQH2</ci> 
        <apply>
          <minus/>
          <ci> UQ_tot </ci>
          <ci> UQ </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="c2plus"> 
    <variable units="micromolar" name="ctot" initial_value="2.7E2"/>
    <variable units="micromolar" public_interface="out" name="c2plus"/>
    <variable units="micromolar" public_interface="in" name="c3plus"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>c2plus</ci> 
        <apply>
          <minus/>
          <ci> ctot </ci>
          <ci> c3plus </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="a2plus"> 
    <variable units="micromolar" name="atot" initial_value="1.35E2"/>
    <variable units="micromolar" public_interface="out" name="a2plus"/>
    <variable units="micromolar" public_interface="in" name="a3_2"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>a2plus</ci> 
        <apply>
          <minus/>
          <ci> atot </ci>
          <ci> a3_2 </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="ADPx"> 
    <variable units="micromolar" name="Ax" initial_value="1.0E4"/>
    <variable units="micromolar" public_interface="in" name="ATPx"/>
    <variable units="micromolar" public_interface="out" name="ADPx"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> ADPx </ci> 
        <apply>
          <minus/>
          <ci> Ax </ci>
          <ci> ATPx </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="pHx"> 
    <variable units="dimensionless" public_interface="out" name="pHx"/>
    <variable units="dimensionless" public_interface="in" name="eta"/>
    <variable units="micromolar" public_interface="in" name="Hx"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> pHx </ci> 
        <apply>
          <minus/>
          <apply>
            <log/>
            <apply>
              <times/>
              <ci> eta </ci>
              <ci> Hx </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="delta_pH"> 
    <variable units="dimensionless" public_interface="out" name="delta_pH"/>
    <variable units="dimensionless" public_interface="in" name="pHx"/>
    <variable units="dimensionless" public_interface="in" name="pHext"/>
    <variable units="millivolt" public_interface="in" name="Z"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> delta_pH </ci>
        <apply>
          <times/>
          <ci> Z </ci> 
          <apply>
            <minus/>
            <ci> pHext </ci>
            <ci> pHx </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="delta_p"> 
    <variable units="dimensionless" public_interface="out" name="delta_p"/>
    <variable units="dimensionless" public_interface="in" name="u"/>
    <variable units="dimensionless" public_interface="in" name="delta_pH"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> delta_p </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> u </ci>
            </apply>
          </apply>
          <ci> delta_pH </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="delta_psi"> 
    <variable units="dimensionless" public_interface="out" name="delta_psi"/>
    <variable units="dimensionless" public_interface="in" name="delta_p"/>
    <variable units="dimensionless" public_interface="in" name="delta_pH"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> delta_psi </ci>
        <apply>
          <minus/>
          <apply>
            <minus/>
            <ci> delta_p </ci>
            <ci> delta_pH </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Hext"> 
    <variable units="micromolar" public_interface="out" name="Hext"/>
    <variable units="dimensionless" public_interface="in" name="eta"/>
    <variable units="dimensionless" public_interface="in" name="pHext"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> Hext </ci>
        <apply>
          <times/>
          <ci> eta </ci>
          <apply>
            <power/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <minus/>
              <ci> pHext </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="r_buff">
    <variable units="molar" public_interface="out" name="r_buff"/> 
    
    <variable units="molar" name="r_buff0" initial_value="0.022"/>
    
    <variable units="dimensionless" public_interface="in" name="delta_pH"/> 
    <variable units="flux" public_interface="in" name="delta_Hx"/>
    <variable units="dimensionless" public_interface="in" name="eta"/>
    <variable units="dimensionless" public_interface="in" name="pHx"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> r_buff </ci>
        <apply>
          <times/>
          <ci> r_buff0 </ci>
          <ci> eta </ci>
          <apply>
            <divide/>
            <ci> delta_pH </ci>
            <ci> delta_Hx </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="Z"> 
    <variable units="millivolt" public_interface="out" name="Z"/>
    
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="kilojoule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kilojoule_per_mole_millivolt" public_interface="in" name="F"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> Z </ci>
        <apply>
          <times/>
          <apply>
            <ln/>
            <cn cellml:units="dimensionless"> 10.0 </cn>
          </apply>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="v_ANT"> 
    <variable units="flux" public_interface="out" name="v_ANT" initial_value="5.70E3"/>
    
    <variable units="flux" name="v_ANT_XG"/>
    <variable units="flux" name="v_ANT_XI"/>
    <variable units="micromolar" name="Kg" initial_value="4.25"/>
    <variable units="micromolar" name="Ki" initial_value="4.25"/>
   
    <variable units="millivolt" public_interface="in" name="Z"/>
    <variable units="dimensionless" public_interface="in" name="delta_psi"/>
    <variable units="micromolar" public_interface="in" name="fADPi"/>
    <variable units="micromolar" public_interface="in" name="fATPi"/>
    <variable units="micromolar" public_interface="in" name="fADPg"/>
    <variable units="micromolar" public_interface="in" name="fATPg"/>
    <variable units="micromolar" public_interface="in" name="fADPx"/>
    <variable units="micromolar" public_interface="in" name="fATPx"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/>
        <ci> v_ANT </ci>
        <apply>
          <plus/>
          <ci> v_ANT_XG </ci>
          <ci> v_ANT_XI </ci>
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <ci> v_ANT_XG </ci>
        <apply>
          <times/>
          <ci> v_ANT </ci>
          <apply>
            <divide/>
            <ci> fADPg </ci>
            <apply>
              <times/>
              <ci> Kg </ci>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 10.0 </cn>
                  <ci> fADPg </ci>
                </apply>
                <apply>
                  <plus/>
                  <ci> Kg </ci>
                  <apply>
                    <divide/>
                    <ci> fADPi </ci>
                    <ci> Ki </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <divide/>
              <ci> fADPg </ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci> fATPg </ci>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 10.0 </cn>
                    <apply>
                      <divide/>
                      <apply>
                        <times/>
                        <cn cellml:units="dimensionless"> -0.35 </cn>
                        <ci> delta_psi </ci>
                      </apply>
                      <ci> Z </ci>
                    </apply>
                  </apply>
                </apply>
                <ci> fADPg </ci>
              </apply>
            </apply>
            <apply>
              <divide/>
              <ci> fADPx </ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci> fATPx </ci>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 10.0 </cn>
                    <apply>
                      <divide/>
                      <apply>
                        <times/>
                        <cn cellml:units="dimensionless"> -0.65 </cn>
                        <ci> delta_psi </ci>
                      </apply>
                      <ci> Z </ci>
                    </apply>
                  </apply>
                </apply>
                <ci> fADPx </ci>
              </apply>
            </apply>
          </apply>        
        </apply>
      </apply>
      
      <apply>
        <eq/>
        <ci> v_ANT_XI </ci>
        <apply>
          <times/>
          <ci> v_ANT </ci>
          <apply>
            <divide/>
            <ci> fADPi </ci>
            <apply>
              <times/>
              <ci> Ki </ci>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <cn cellml:units="dimensionless"> 10.0 </cn>
                  <ci> fADPg </ci>
                </apply>
                <apply>
                  <plus/>
                  <ci> Kg </ci>
                  <apply>
                    <divide/>
                    <ci> fADPi </ci>
                    <ci> Ki </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <divide/>
              <ci> fADPi </ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci> fATPi </ci>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 10.0 </cn>
                    <apply>
                      <divide/>
                      <apply>
                        <times/>
                        <cn cellml:units="dimensionless"> -0.35 </cn>
                        <ci> delta_psi </ci>
                      </apply>
                      <ci> Z </ci>
                    </apply>
                  </apply>
                </apply>
                <ci> fADPi </ci>
              </apply>
            </apply>
            <apply>
              <divide/>
              <ci> fADPx </ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci> fATPx </ci>
                  <apply>
                    <power/>
                    <cn cellml:units="dimensionless"> 10.0 </cn>
                    <apply>
                      <divide/>
                      <apply>
                        <times/>
                        <cn cellml:units="dimensionless"> -0.65 </cn>
                        <ci> delta_psi </ci>
                      </apply>
                      <ci> Z </ci>
                    </apply>
                  </apply>
                </apply>
                <ci> fADPx </ci>
              </apply>
            </apply>
          </apply>        
        </apply>
      </apply>
    </math>
  </component>
 
  <component name="delta_Gsn"> 
    <variable units="dimensionless" public_interface="out" name="delta_Gsn"/> 
    
    <variable units="kilojoule_per_mole" name="delta_Gsn0" initial_value="31.9"/>
    
    <variable units="dimensionless" public_interface="in" name="delta_p"/>
    <variable units="dimensionless" public_interface="in" name="na"/>
    <variable units="millivolt" public_interface="in" name="Z"/>
    <variable units="dimensionless" public_interface="in" name="eta"/>
    <variable units="kilojoule_per_mole_millivolt" public_interface="in" name="F"/>
    <variable units="micromolar" public_interface="in" name="ADPx"/>
    <variable units="micromolar" public_interface="in" name="ATPx"/>
    <variable units="micromolar" public_interface="in" name="Pix"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> delta_Gsn </ci>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> na </ci>
              <ci> delta_p </ci>
            </apply>
            <apply>
              <divide/>
              <ci> delta_Gsn0 </ci>
              <ci> F </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> Z </ci>
            <apply>
              <log/>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> eta </ci>
                  <ci> ATPx </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> ADPx </ci>
                  <ci> Pix </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="redox_potentials"> 
    <variable units="millivolt" public_interface="out" name="En"/>
    <variable units="millivolt" public_interface="out" name="Eu"/>
    <variable units="millivolt" public_interface="out" name="Ec"/>
    <variable units="millivolt" public_interface="out" name="Ea"/>
    <variable units="micromolar" public_interface="out" name="a3_2"/>
    
    <variable units="millivolt" name="En0" initial_value="-320.0"/>
    <variable units="millivolt" name="Eu0" initial_value="85.0"/>
    <variable units="millivolt" name="Ec0" initial_value="250.0"/>
    <variable units="millivolt" name="Ea0" initial_value="540.0"/> 
    
    <variable units="dimensionless" public_interface="in" name="delta_p"/>
    <variable units="dimensionless" public_interface="in" name="u"/>
    <variable units="millivolt" public_interface="in" name="Z"/>
    <variable units="micromolar" public_interface="in" name="NAD"/>
    <variable units="micromolar" public_interface="in" name="NADH"/>
    <variable units="micromolar" public_interface="in" name="UQ"/>
    <variable units="micromolar" public_interface="in" name="UQH2"/>
    <variable units="micromolar" public_interface="in" name="c3plus"/>
    <variable units="micromolar" public_interface="in" name="c2plus"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> En </ci>
        <apply>
          <plus/>
          <ci> En0 </ci>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> Z </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <log/>
              <apply>
                <divide/>
                <ci> NAD </ci>
                <ci> NADH </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>  
      
      <apply>
        <eq/> 
        <ci> Eu </ci>
        <apply>
          <plus/>
          <ci> Eu0 </ci>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> Z </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <log/>
              <apply>
                <divide/>
                <ci> UQ </ci>
                <ci> UQH2 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply>
        <eq/> 
        <ci> Ec </ci>
        <apply>
          <plus/>
          <ci> Ec0 </ci>
          <apply>
            <times/>
            <apply>
              <divide/>
              <ci> Z </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <log/>
              <apply>
                <divide/>
                <ci> c3plus </ci>
                <ci> c2plus </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>  
      
      <apply>
        <eq/> 
        <ci> Ea </ci>
        <apply>
          <plus/>
          <ci> Ec </ci>
          <apply>
            <times/>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> u </ci>
            </apply>
            <ci> delta_p </ci>
          </apply>
        </apply>
      </apply>    
      
      <apply>
        <eq/> 
        <ci> a3_2 </ci>
        <apply>
          <power/>
          <cn cellml:units="micromolar"> 10.0 </cn>
          <apply>
            <times/>
            <apply>
              <minus/>
              <ci> Ea </ci>
              <ci> Ea0 </ci>
            </apply>
            <ci> Z </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="diffusion_constants"> 
    <variable units="micromolar2_per_s" public_interface="out" name="D_ATP" initial_value="1.45E2"/>
    <variable units="micromolar2_per_s" public_interface="out" name="D_ADP" initial_value="1.45E2"/>
    <variable units="micromolar2_per_s" public_interface="out" name="D_AMP" initial_value="1.45E2"/>
    <variable units="micromolar2_per_s" public_interface="out" name="D_PCr" initial_value="2.6E2"/>
    <variable units="micromolar2_per_s" public_interface="out" name="D_Cr" initial_value="2.6E2"/>
    <variable units="micromolar2_per_s" public_interface="out" name="D_Pi" initial_value="3.27E2"/>
  </component>
  
  <component name="Mg_buffering"> 
    <variable units="micromolar" public_interface="out" name="Mg"/>
    <variable units="micromolar" public_interface="out" name="Mgg"/>
    <variable units="micromolar" public_interface="out" name="Mgi"/>
    <variable units="micromolar" public_interface="out" name="Mgext" initial_value="1.0E3"/>
    <variable units="micromolar" public_interface="out" name="Mgx" initial_value="380.0"/>
    <variable units="micromolar" public_interface="out" name="ATPg"/>
    <variable units="micromolar" public_interface="out" name="ADPg"/>
    <variable units="micromolar" public_interface="out" name="ATPext"/>
    <variable units="micromolar" public_interface="out" name="ADPext"/>
  </component>
  
  <component name="H_buffering"> 
    <variable units="dimensionless" public_interface="out" name="pHext" initial_value="6.9"/>
  </component>
  
  <component name="diffusion_through_mitochondrial_membrane_parameters"> 
    <variable units="micrometre" public_interface="out" name="Li" initial_value="0.1"/>
    <variable units="micrometre_per_s" public_interface="out" name="R_ATP" initial_value="0.16"/>
    <variable units="micrometre_per_s" public_interface="out" name="R_ADP" initial_value="0.16"/>
    <variable units="micrometre_per_s" public_interface="out" name="R_AMP" initial_value="0.16"/>
    <variable units="micrometre_per_s" public_interface="out" name="R_PCr" initial_value="260.0"/>
    <variable units="micrometre_per_s" public_interface="out" name="R_Cr" initial_value="260.0"/>
    <variable units="micrometre_per_s" public_interface="out" name="R_Pi" initial_value="327.0"/>
  </component>
  
  <component name="general_constants"> 
    <variable units="kilojoule_per_mole_kelvin" public_interface="out" name="R" initial_value="0.0083"/>
    <variable units="kelvin" public_interface="out" name="T" initial_value="298.0"/>
    <variable units="kilojoule_per_mole_millivolt" public_interface="out" name="F" initial_value="0.0965"/>
    <variable units="dimensionless" public_interface="out" name="u" initial_value="0.861"/>
    <variable units="dimensionless" public_interface="out" name="eta" initial_value="10E6"/>
    <variable units="dimensionless" public_interface="out" name="na" initial_value="2.5"/>
  </component>
  
  <component name="fractional_volumes"> 
    <variable units="dimensionless" public_interface="out" name="FVx" initial_value="0.1875"/>
    <variable units="dimensionless" public_interface="out" name="FVi" initial_value="0.0625"/>
  </component>
  
  
  
  <connection>
    <map_components component_2="G_CK" component_1="ATP"/>
    <map_variables variable_2="G_CK" variable_1="G_CK"/>
  </connection>
  
  <connection>
    <map_components component_2="v_CK" component_1="ATP"/>
    <map_variables variable_2="v_CK" variable_1="v_CK"/>
  </connection>
  
  <connection>
    <map_components component_2="G_AK" component_1="ATP"/>
    <map_variables variable_2="G_AK" variable_1="G_AK"/>
  </connection>
  
  <connection>
    <map_components component_2="v_AK" component_1="ATP"/>
    <map_variables variable_2="v_AK" variable_1="v_AK"/>
  </connection>
  
  <connection>
    <map_components component_2="G_H" component_1="ATP"/>
    <map_variables variable_2="G_H" variable_1="G_H"/>
  </connection>
  
  <connection>
    <map_components component_2="H_ATP" component_1="ATP"/>
    <map_variables variable_2="H_ATP" variable_1="H_ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="diffusion_constants" component_1="ATP"/>
    <map_variables variable_2="D_ATP" variable_1="D_ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="G_CK" component_1="ADP"/>
    <map_variables variable_2="G_CK" variable_1="G_CK"/>
  </connection>
  
  <connection>
    <map_components component_2="v_CK" component_1="ADP"/>
    <map_variables variable_2="v_CK" variable_1="v_CK"/>
  </connection>
  
  <connection>
    <map_components component_2="G_AK" component_1="ADP"/>
    <map_variables variable_2="G_AK" variable_1="G_AK"/>
  </connection>
  
  <connection>
    <map_components component_2="v_AK" component_1="ADP"/>
    <map_variables variable_2="v_AK" variable_1="v_AK"/>
  </connection>
  
  <connection>
    <map_components component_2="G_H" component_1="ADP"/>
    <map_variables variable_2="G_H" variable_1="G_H"/>
  </connection>
  
  <connection>
    <map_components component_2="H_ATP" component_1="ADP"/>
    <map_variables variable_2="H_ATP" variable_1="H_ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="diffusion_constants" component_1="ADP"/>
    <map_variables variable_2="D_ADP" variable_1="D_ADP"/>
  </connection>
  
  <connection>
    <map_components component_2="G_AK" component_1="AMP"/>
    <map_variables variable_2="G_AK" variable_1="G_AK"/>
  </connection>
  
  <connection>
    <map_components component_2="v_AK" component_1="AMP"/>
    <map_variables variable_2="v_AK" variable_1="v_AK"/>
    <map_variables variable_2="AMP" variable_1="AMP"/>
  </connection>
  
  <connection>
    <map_components component_2="diffusion_constants" component_1="AMP"/>
    <map_variables variable_2="D_AMP" variable_1="D_AMP"/>
  </connection>
  
  <connection>
    <map_components component_2="G_CK" component_1="G_AK"/>
    <map_variables variable_2="G_CK" variable_1="G_CK"/>
  </connection>
  
  <connection>
    <map_components component_2="G_CK" component_1="PCr"/>
    <map_variables variable_2="G_CK" variable_1="G_CK"/>
  </connection>
  
  <connection>
    <map_components component_2="v_CK" component_1="PCr"/>
    <map_variables variable_2="v_CK" variable_1="v_CK"/>
    <map_variables variable_2="PCr" variable_1="PCr"/>
  </connection>
  
  <connection>
    <map_components component_2="diffusion_constants" component_1="PCr"/>
    <map_variables variable_2="D_PCr" variable_1="D_PCr"/>
  </connection>
  
  <connection>
    <map_components component_2="G_CK" component_1="Cr"/>
    <map_variables variable_2="G_CK" variable_1="G_CK"/>
  </connection>
  
  <connection>
    <map_components component_2="v_CK" component_1="Cr"/>
    <map_variables variable_2="v_CK" variable_1="v_CK"/>
    <map_variables variable_2="Cr" variable_1="Cr"/>
  </connection>
  
  <connection>
    <map_components component_2="diffusion_constants" component_1="Cr"/>
    <map_variables variable_2="D_Cr" variable_1="D_Cr"/>
  </connection>
  
  <connection>
    <map_components component_2="G_H" component_1="Pi"/>
    <map_variables variable_2="G_H" variable_1="G_H"/>
  </connection>
  
  <connection>
    <map_components component_2="H_ATP" component_1="Pi"/>
    <map_variables variable_2="H_ATP" variable_1="H_ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="diffusion_constants" component_1="Pi"/>
    <map_variables variable_2="D_Pi" variable_1="D_Pi"/>
  </connection>
  
  <connection>
    <map_components component_2="mATP" component_1="v_CK"/>
    <map_variables variable_2="mATP" variable_1="mATP"/>
  </connection>
  
  <connection>
    <map_components component_2="mADP" component_1="v_CK"/>
    <map_variables variable_2="mADP" variable_1="mADP"/>
  </connection>
  
  <connection>
    <map_components component_2="ATP" component_1="fATP"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="fATP"/>
    <map_variables variable_2="Mg" variable_1="Mg"/>
  </connection>
  
  <connection>
    <map_components component_2="fATP" component_1="mATP"/>
    <map_variables variable_2="fATP" variable_1="fATP"/>
  </connection>
  
  <connection>
    <map_components component_2="ATP" component_1="mATP"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="ADP" component_1="fADP"/>
    <map_variables variable_2="ADP" variable_1="ADP"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="fADP"/>
    <map_variables variable_2="Mg" variable_1="Mg"/>
  </connection>
  
  <connection>
    <map_components component_2="fADP" component_1="mADP"/>
    <map_variables variable_2="fADP" variable_1="fADP"/>
  </connection>
  
  <connection>
    <map_components component_2="ADP" component_1="mADP"/>
    <map_variables variable_2="ADP" variable_1="ADP"/>
  </connection>
  
  <connection>
    <map_components component_2="fADP" component_1="v_AK"/>
    <map_variables variable_2="fADP" variable_1="fADP"/>
  </connection>
  
  <connection>
    <map_components component_2="mATP" component_1="v_AK"/>
    <map_variables variable_2="mATP" variable_1="mATP"/>
  </connection>
  
  <connection>
    <map_components component_2="mADP" component_1="v_AK"/>
    <map_variables variable_2="mADP" variable_1="mADP"/>
  </connection>
  
  <connection>
    <map_components component_2="fADP" component_1="v_AKmit"/>
    <map_variables variable_2="fADP" variable_1="fADP"/>
  </connection>
  
  <connection>
    <map_components component_2="mATP" component_1="v_AKmit"/>
    <map_variables variable_2="mATP" variable_1="mATP"/>
  </connection>
  
  <connection>
    <map_components component_2="mADP" component_1="v_AKmit"/>
    <map_variables variable_2="mADP" variable_1="mADP"/>
  </connection>
  
  <connection>
    <map_components component_2="AMP" component_1="v_AKmit"/>
    <map_variables variable_2="AMP" variable_1="AMP"/>
  </connection>
  
  <connection>
    <map_components component_2="ATP" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPi" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="ATPi" variable_1="ATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="ADP" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="ADP" variable_1="ADP"/>
  </connection>
  
  <connection>
    <map_components component_2="ADPi" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="ADPi" variable_1="ADPi"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPi" component_1="diffusion_through_mitochondrial_membrane_parameters"/>
    <map_variables variable_2="R_ATP" variable_1="R_ATP"/>
    <map_variables variable_2="Li" variable_1="Li"/>
  </connection>
  
  <connection>
    <map_components component_2="ADPi" component_1="diffusion_through_mitochondrial_membrane_parameters"/>
    <map_variables variable_2="R_ADP" variable_1="R_ADP"/>
    <map_variables variable_2="Li" variable_1="Li"/>
  </connection>
  
  <connection>
    <map_components component_2="AMPi" component_1="diffusion_through_mitochondrial_membrane_parameters"/>
    <map_variables variable_2="R_AMP" variable_1="R_AMP"/>
    <map_variables variable_2="Li" variable_1="Li"/>
  </connection>
  
  <connection>
    <map_components component_2="PCri" component_1="diffusion_through_mitochondrial_membrane_parameters"/>
    <map_variables variable_2="R_PCr" variable_1="R_PCr"/>
    <map_variables variable_2="Li" variable_1="Li"/>
  </connection>
  
  <connection>
    <map_components component_2="Cri" component_1="diffusion_through_mitochondrial_membrane_parameters"/>
    <map_variables variable_2="R_Cr" variable_1="R_Cr"/>
    <map_variables variable_2="Li" variable_1="Li"/>
  </connection>
  
  <connection>
    <map_components component_2="Pii" component_1="diffusion_through_mitochondrial_membrane_parameters"/>
    <map_variables variable_2="R_Pi" variable_1="R_Pi"/>
    <map_variables variable_2="Li" variable_1="Li"/>
  </connection>
  
  <connection>
    <map_components component_2="diffusion_through_mitochondrial_membrane" component_1="diffusion_through_mitochondrial_membrane_parameters"/>
    <map_variables variable_2="R_Pi" variable_1="R_Pi"/>
    <map_variables variable_2="R_Cr" variable_1="R_Cr"/>
    <map_variables variable_2="R_PCr" variable_1="R_PCr"/>
    <map_variables variable_2="R_AMP" variable_1="R_AMP"/>
    <map_variables variable_2="R_ADP" variable_1="R_ADP"/>
    <map_variables variable_2="R_ATP" variable_1="R_ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="AMP" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="AMP" variable_1="AMP"/>
  </connection>
  
  <connection>
    <map_components component_2="AMPi" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="AMPi" variable_1="AMPi"/>
  </connection>
  
  <connection>
    <map_components component_2="PCr" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="PCr" variable_1="PCr"/>
  </connection>
  
  <connection>
    <map_components component_2="PCri" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="PCri" variable_1="PCri"/>
  </connection>
  
  <connection>
    <map_components component_2="Cr" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="Cr" variable_1="Cr"/>
  </connection>
  
  <connection>
    <map_components component_2="Cri" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="Cri" variable_1="Cri"/>
  </connection>
  
  <connection>
    <map_components component_2="Pi" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="Pi" variable_1="Pi"/>
  </connection>
  
  <connection>
    <map_components component_2="Pii" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="Pii" variable_1="Pii"/>
  </connection>
  
  <connection>
    <map_components component_2="diffusion_constants" component_1="diffusion_through_mitochondrial_membrane"/>
    <map_variables variable_2="D_ATP" variable_1="D_ATP"/>
    <map_variables variable_2="D_ADP" variable_1="D_ADP"/>
    <map_variables variable_2="D_AMP" variable_1="D_AMP"/>
    <map_variables variable_2="D_PCr" variable_1="D_PCr"/>
    <map_variables variable_2="D_Cr" variable_1="D_Cr"/>
    <map_variables variable_2="D_Pi" variable_1="D_Pi"/>
  </connection>
  
  <connection>
    <map_components component_2="ATP" component_1="ATPi"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
  </connection>
  
  <connection>
    <map_components component_2="v_ANT" component_1="ATPi"/>
    <map_variables variable_2="v_ANT" variable_1="v_ANT"/>
  </connection>
  
  <connection>
    <map_components component_2="v_MiCK" component_1="ATPi"/>
    <map_variables variable_2="v_MiCK" variable_1="v_MiCK"/>
  </connection>
  
  <connection>
    <map_components component_2="v_AKmit" component_1="ATPi"/>
    <map_variables variable_2="v_AKmit" variable_1="v_AKmit"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="ATPi"/>
    <map_variables variable_2="FVi" variable_1="FVi"/>
  </connection>
  
  <connection>
    <map_components component_2="ADP" component_1="ADPi"/>
    <map_variables variable_2="ADP" variable_1="ADP"/>
  </connection>
  
  <connection>
    <map_components component_2="v_ANT" component_1="ADPi"/>
    <map_variables variable_2="v_ANT" variable_1="v_ANT"/>
  </connection>
  
  <connection>
    <map_components component_2="v_MiCK" component_1="ADPi"/>
    <map_variables variable_2="v_MiCK" variable_1="v_MiCK"/>
  </connection>
  
  <connection>
    <map_components component_2="v_AKmit" component_1="ADPi"/>
    <map_variables variable_2="v_AKmit" variable_1="v_AKmit"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="ADPi"/>
    <map_variables variable_2="FVi" variable_1="FVi"/>
  </connection>
  
  <connection>
    <map_components component_2="AMP" component_1="AMPi"/>
    <map_variables variable_2="AMP" variable_1="AMP"/>
  </connection>
  
  <connection>
    <map_components component_2="v_AKmit" component_1="AMPi"/>
    <map_variables variable_2="v_AKmit" variable_1="v_AKmit"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="AMPi"/>
    <map_variables variable_2="FVi" variable_1="FVi"/>
  </connection>
  
  <connection>
    <map_components component_2="PCr" component_1="PCri"/>
    <map_variables variable_2="PCr" variable_1="PCr"/>
  </connection>
  
  <connection>
    <map_components component_2="v_AKmit" component_1="PCri"/>
    <map_variables variable_2="v_AKmit" variable_1="v_AKmit"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="PCri"/>
    <map_variables variable_2="FVi" variable_1="FVi"/>
  </connection>
  
  <connection>
    <map_components component_2="Cr" component_1="Cri"/>
    <map_variables variable_2="Cr" variable_1="Cr"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="Cri"/>
    <map_variables variable_2="FVi" variable_1="FVi"/>
  </connection>
  
  <connection>
    <map_components component_2="v_MiCK" component_1="Cri"/>
    <map_variables variable_2="v_MiCK" variable_1="v_MiCK"/>
    <map_variables variable_2="Cri" variable_1="Cri"/>
  </connection>
  
  <connection>
    <map_components component_2="Pi" component_1="Pii"/>
    <map_variables variable_2="Pi" variable_1="Pi"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="Pii"/>
    <map_variables variable_2="FVi" variable_1="FVi"/>
  </connection>
  
  <connection>
    <map_components component_2="oxidative_phosphorylation" component_1="Pii"/>
    <map_variables variable_2="v_PI" variable_1="v_PI"/>
    <map_variables variable_2="Pii" variable_1="Pii"/>
  </connection>
  
  <connection>
    <map_components component_2="PCri" component_1="v_MiCK"/>
    <map_variables variable_2="PCri" variable_1="PCri"/>
  </connection>
  
  <connection>
    <map_components component_2="mATPg" component_1="v_MiCK"/>
    <map_variables variable_2="mATPg" variable_1="mATPg"/>
  </connection>
  
  <connection>
    <map_components component_2="mATPi" component_1="v_MiCK"/>
    <map_variables variable_2="mATPi" variable_1="mATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="mADPg" component_1="v_MiCK"/>
    <map_variables variable_2="mADPg" variable_1="mADPg"/>
  </connection>
  
  <connection>
    <map_components component_2="mADPi" component_1="v_MiCK"/>
    <map_variables variable_2="mADPi" variable_1="mADPi"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPi" component_1="v_exch"/>
    <map_variables variable_2="ATPi" variable_1="ATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="v_exch"/>
    <map_variables variable_2="ATPg" variable_1="ATPg"/>
    <map_variables variable_2="ADPg" variable_1="ADPg"/>
  </connection>
  
  <connection>
    <map_components component_2="ADPi" component_1="v_exch"/>
    <map_variables variable_2="ADPi" variable_1="ADPi"/>
  </connection>
  
  <connection>
    <map_components component_2="delta_p" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="delta_p" variable_1="delta_p"/>
  </connection>
  
  <connection>
    <map_components component_2="delta_Gsn" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="delta_Gsn" variable_1="delta_Gsn"/>
  </connection>
  
  <connection>
    <map_components component_2="Z" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="Z" variable_1="Z"/>
  </connection>
  
  <connection>
    <map_components component_2="redox_potentials" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="Ea" variable_1="Ea"/>
    <map_variables variable_2="Eu" variable_1="Eu"/>
    <map_variables variable_2="Ec" variable_1="Ec"/>
    <map_variables variable_2="En" variable_1="En"/>
  </connection>
  
  <connection>
    <map_components component_2="NADH" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="NADH" variable_1="NADH"/>
  </connection>
  
  <connection>
    <map_components component_2="NAD" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="NAD" variable_1="NAD"/>
    <map_variables variable_2="v_C1" variable_1="v_C1"/>
    <map_variables variable_2="v_dh" variable_1="v_dh"/>
  </connection>
  
  <connection>
    <map_components component_2="c2plus" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="c2plus" variable_1="c2plus"/>
  </connection>
  
  <connection>
    <map_components component_2="a2plus" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="a2plus" variable_1="a2plus"/>
  </connection>
  
  <connection>
    <map_components component_2="Hext" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="Hext" variable_1="Hext"/>
  </connection>
  
  <connection>
    <map_components component_2="pHx" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="pHx" variable_1="pHx"/>
  </connection>
  
  <connection>
    <map_components component_2="H_buffering" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="pHext" variable_1="pHext"/>
  </connection>
  
  <connection>
    <map_components component_2="general_constants" component_1="oxidative_phosphorylation"/>
    <map_variables variable_2="u" variable_1="u"/>
  </connection>
 
  <connection>
    <map_components component_2="fractional_volumes" component_1="UQ"/>
    <map_variables variable_2="FVx" variable_1="FVx"/>
  </connection>
  
  <connection>
    <map_components component_2="oxidative_phosphorylation" component_1="UQ"/>
    <map_variables variable_2="v_C1" variable_1="v_C1"/>
    <map_variables variable_2="v_C3" variable_1="v_C3"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="c3plus"/>
    <map_variables variable_2="FVx" variable_1="FVx"/>
  </connection>
  
  <connection>
    <map_components component_2="oxidative_phosphorylation" component_1="c3plus"/>
    <map_variables variable_2="v_O2" variable_1="v_O2"/>
    <map_variables variable_2="v_C3" variable_1="v_C3"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="Hx"/>
    <map_variables variable_2="FVx" variable_1="FVx"/>
  </connection>
  
  <connection>
    <map_components component_2="general_constants" component_1="Hx"/>
    <map_variables variable_2="u" variable_1="u"/>
    <map_variables variable_2="na" variable_1="na"/>
  </connection>
  
  <connection>
    <map_components component_2="oxidative_phosphorylation" component_1="Hx"/>
    <map_variables variable_2="v_O2" variable_1="v_O2"/>
    <map_variables variable_2="v_C1" variable_1="v_C1"/>
    <map_variables variable_2="v_C3" variable_1="v_C3"/>
    <map_variables variable_2="v_PI" variable_1="v_PI"/>
    <map_variables variable_2="v_sn" variable_1="v_sn"/>
    <map_variables variable_2="v_leak" variable_1="v_leak"/>
    <map_variables variable_2="Hx" variable_1="Hx"/>
  </connection>
  
  <connection>
    <map_components component_2="v_ANT" component_1="Hx"/>
    <map_variables variable_2="v_ANT" variable_1="v_ANT"/>
  </connection>
  
  <connection>
    <map_components component_2="r_buff" component_1="Hx"/>
    <map_variables variable_2="r_buff" variable_1="r_buff"/>
    <map_variables variable_2="delta_Hx" variable_1="delta_Hx"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="NAD"/>
    <map_variables variable_2="FVx" variable_1="FVx"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="ATPx"/>
    <map_variables variable_2="FVx" variable_1="FVx"/>
  </connection>
  
  <connection>
    <map_components component_2="oxidative_phosphorylation" component_1="ATPx"/>
    <map_variables variable_2="v_sn" variable_1="v_sn"/>
  </connection>
  
  <connection>
    <map_components component_2="v_ANT" component_1="ATPx"/>
    <map_variables variable_2="v_ANT" variable_1="v_ANT"/>
  </connection>
  
  <connection>
    <map_components component_2="fractional_volumes" component_1="Pix"/>
    <map_variables variable_2="FVx" variable_1="FVx"/>
  </connection>
  
  <connection>
    <map_components component_2="oxidative_phosphorylation" component_1="Pix"/>
    <map_variables variable_2="v_PI" variable_1="v_PI"/>
    <map_variables variable_2="Pix" variable_1="Pix"/>
    <map_variables variable_2="v_sn" variable_1="v_sn"/>
  </connection>
  
  <connection>
    <map_components component_2="NAD" component_1="NADH"/>
    <map_variables variable_2="NAD" variable_1="NAD"/>
  </connection>
  
  <connection>
    <map_components component_2="UQ" component_1="UQH2"/>
    <map_variables variable_2="UQ" variable_1="UQ"/>
  </connection>
  
  <connection>
    <map_components component_2="c3plus" component_1="c2plus"/>
    <map_variables variable_2="c3plus" variable_1="c3plus"/>
  </connection>
  
  <connection>
    <map_components component_2="redox_potentials" component_1="a2plus"/>
    <map_variables variable_2="a3_2" variable_1="a3_2"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPx" component_1="ADPx"/>
    <map_variables variable_2="ATPx" variable_1="ATPx"/>
  </connection>
  
  <connection>
    <map_components component_2="Hx" component_1="pHx"/>
    <map_variables variable_2="Hx" variable_1="Hx"/>
  </connection>
  
  <connection>
    <map_components component_2="general_constants" component_1="pHx"/>
    <map_variables variable_2="eta" variable_1="eta"/>
  </connection>
  
  <connection>
    <map_components component_2="H_buffering" component_1="delta_pH"/>
    <map_variables variable_2="pHext" variable_1="pHext"/>
  </connection>
  
  <connection>
    <map_components component_2="pHx" component_1="delta_pH"/>
    <map_variables variable_2="pHx" variable_1="pHx"/>
  </connection>
  
  <connection>
    <map_components component_2="Z" component_1="delta_pH"/>
    <map_variables variable_2="Z" variable_1="Z"/>
  </connection>
  
  <connection>
    <map_components component_2="delta_pH" component_1="delta_p"/>
    <map_variables variable_2="delta_pH" variable_1="delta_pH"/>
  </connection>
  
  <connection>
    <map_components component_2="general_constants" component_1="delta_p"/>
    <map_variables variable_2="u" variable_1="u"/>
  </connection>
  
  <connection>
    <map_components component_2="delta_p" component_1="delta_psi"/>
    <map_variables variable_2="delta_p" variable_1="delta_p"/>
  </connection>
  
  <connection>
    <map_components component_2="delta_pH" component_1="delta_psi"/>
    <map_variables variable_2="delta_pH" variable_1="delta_pH"/>
  </connection>
  
  <connection>
    <map_components component_2="general_constants" component_1="Hext"/>
    <map_variables variable_2="eta" variable_1="eta"/>
  </connection>
  
  <connection>
    <map_components component_2="H_buffering" component_1="Hext"/>
    <map_variables variable_2="pHext" variable_1="pHext"/>
  </connection>
  
  <connection>
    <map_components component_2="general_constants" component_1="r_buff"/>
    <map_variables variable_2="eta" variable_1="eta"/>
  </connection>
  
  <connection>
    <map_components component_2="pHx" component_1="r_buff"/>
    <map_variables variable_2="pHx" variable_1="pHx"/>
  </connection>
  
  <connection>
    <map_components component_2="delta_pH" component_1="r_buff"/>
    <map_variables variable_2="delta_pH" variable_1="delta_pH"/>
  </connection>
  
  <connection>
    <map_components component_2="general_constants" component_1="Z"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="T" variable_1="T"/>
  </connection>
  
  <connection>
    <map_components component_2="fADPg" component_1="v_ANT"/>
    <map_variables variable_2="fADPg" variable_1="fADPg"/>
  </connection>
  
  <connection>
    <map_components component_2="fATPg" component_1="v_ANT"/>
    <map_variables variable_2="fATPg" variable_1="fATPg"/>
  </connection>
  
  <connection>
    <map_components component_2="fADPi" component_1="v_ANT"/>
    <map_variables variable_2="fADPi" variable_1="fADPi"/>
  </connection>
  
  <connection>
    <map_components component_2="fADPx" component_1="v_ANT"/>
    <map_variables variable_2="fADPx" variable_1="fADPx"/>
  </connection>
  
  <connection>
    <map_components component_2="fATPx" component_1="v_ANT"/>
    <map_variables variable_2="fATPx" variable_1="fATPx"/>
  </connection>
  
  <connection>
    <map_components component_2="fATPi" component_1="v_ANT"/>
    <map_variables variable_2="fATPi" variable_1="fATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="delta_psi" component_1="v_ANT"/>
    <map_variables variable_2="delta_psi" variable_1="delta_psi"/>
  </connection>
  
  <connection>
    <map_components component_2="Z" component_1="v_ANT"/>
    <map_variables variable_2="Z" variable_1="Z"/>
  </connection>
  
  <connection>
    <map_components component_2="delta_p" component_1="delta_Gsn"/>
    <map_variables variable_2="delta_p" variable_1="delta_p"/>
  </connection>
  
  <connection>
    <map_components component_2="Z" component_1="delta_Gsn"/>
    <map_variables variable_2="Z" variable_1="Z"/>
  </connection>
  
  <connection>
    <map_components component_2="general_constants" component_1="delta_Gsn"/>
    <map_variables variable_2="eta" variable_1="eta"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="na" variable_1="na"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPx" component_1="delta_Gsn"/>
    <map_variables variable_2="ATPx" variable_1="ATPx"/>
  </connection>
  
  <connection>
    <map_components component_2="ADPx" component_1="delta_Gsn"/>
    <map_variables variable_2="ADPx" variable_1="ADPx"/>
  </connection>
  
  <connection>
    <map_components component_2="Pix" component_1="delta_Gsn"/>
    <map_variables variable_2="Pix" variable_1="Pix"/>
  </connection>
  
  <connection>
    <map_components component_2="NAD" component_1="redox_potentials"/>
    <map_variables variable_2="NAD" variable_1="NAD"/>
  </connection>
  
  <connection>
    <map_components component_2="NADH" component_1="redox_potentials"/>
    <map_variables variable_2="NADH" variable_1="NADH"/>
  </connection>
  
  <connection>
    <map_components component_2="UQ" component_1="redox_potentials"/>
    <map_variables variable_2="UQ" variable_1="UQ"/>
  </connection>
  
  <connection>
    <map_components component_2="UQH2" component_1="redox_potentials"/>
    <map_variables variable_2="UQH2" variable_1="UQH2"/>
  </connection>
  
  <connection>
    <map_components component_2="c3plus" component_1="redox_potentials"/>
    <map_variables variable_2="c3plus" variable_1="c3plus"/>
  </connection>
  
  <connection>
    <map_components component_2="c2plus" component_1="redox_potentials"/>
    <map_variables variable_2="c2plus" variable_1="c2plus"/>
  </connection>
  
  <connection>
    <map_components component_2="Z" component_1="redox_potentials"/>
    <map_variables variable_2="Z" variable_1="Z"/>
  </connection>
  
  <connection>
    <map_components component_2="delta_p" component_1="redox_potentials"/>
    <map_variables variable_2="delta_p" variable_1="delta_p"/>
  </connection>
  
  <connection>
    <map_components component_2="general_constants" component_1="redox_potentials"/>
    <map_variables variable_2="u" variable_1="u"/>
  </connection>
  
  <connection>
    <map_components component_2="ATP" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="ADP" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="AMP" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="PCr" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="Cr" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="Pi" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="H_ATP" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPi" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="ADPi" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="AMPi" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="PCri" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="Cri" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="Pii" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="UQ" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="c3plus" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="NAD" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="Hx" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPx" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="Pix" component_1="environment"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPi" component_1="fATPi"/>
    <map_variables variable_2="ATPi" variable_1="ATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="fATPi"/>
    <map_variables variable_2="Mgi" variable_1="Mgi"/>
  </connection>
  
  <connection>
    <map_components component_2="fATPi" component_1="mATPi"/>
    <map_variables variable_2="fATPi" variable_1="fATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPi" component_1="mATPi"/>
    <map_variables variable_2="ATPi" variable_1="ATPi"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="fATPg"/>
    <map_variables variable_2="ATPg" variable_1="ATPg"/>
    <map_variables variable_2="Mgg" variable_1="Mgg"/>
  </connection>
  
  <connection>
    <map_components component_2="fATPg" component_1="mATPg"/>
    <map_variables variable_2="fATPg" variable_1="fATPg"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="mATPg"/>
    <map_variables variable_2="ATPg" variable_1="ATPg"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="fATPext"/>
    <map_variables variable_2="Mgext" variable_1="Mgext"/>
    <map_variables variable_2="ATPext" variable_1="ATPext"/>
  </connection>
  
  <connection>
    <map_components component_2="fATPext" component_1="mATPext"/>
    <map_variables variable_2="fATPext" variable_1="fATPext"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="mATPext"/>
    <map_variables variable_2="ATPext" variable_1="ATPext"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPx" component_1="fATPx"/>
    <map_variables variable_2="ATPx" variable_1="ATPx"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="fATPx"/>
    <map_variables variable_2="Mgx" variable_1="Mgx"/>
  </connection>
  
  <connection>
    <map_components component_2="fATPx" component_1="mATPx"/>
    <map_variables variable_2="fATPx" variable_1="fATPx"/>
  </connection>
  
  <connection>
    <map_components component_2="ATPx" component_1="mATPx"/>
    <map_variables variable_2="ATPx" variable_1="ATPx"/>
  </connection>
  
  <connection>
    <map_components component_2="ADPi" component_1="fADPi"/>
    <map_variables variable_2="ADPi" variable_1="ADPi"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="fADPi"/>
    <map_variables variable_2="Mgi" variable_1="Mgi"/>
  </connection>
  
  <connection>
    <map_components component_2="fADPi" component_1="mADPi"/>
    <map_variables variable_2="fADPi" variable_1="fADPi"/>
  </connection>
  
  <connection>
    <map_components component_2="ADPi" component_1="mADPi"/>
    <map_variables variable_2="ADPi" variable_1="ADPi"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="fADPg"/>
    <map_variables variable_2="ADPg" variable_1="ADPg"/>
    <map_variables variable_2="Mgg" variable_1="Mgg"/>
  </connection>
  
  <connection>
    <map_components component_2="fADPg" component_1="mADPg"/>
    <map_variables variable_2="fADPg" variable_1="fADPg"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="mADPg"/>
    <map_variables variable_2="ADPg" variable_1="ADPg"/>
  </connection>
 
  <connection>
    <map_components component_2="Mg_buffering" component_1="fADPext"/>
    <map_variables variable_2="Mgext" variable_1="Mgext"/>
    <map_variables variable_2="ADPext" variable_1="ADPext"/>
  </connection>
  
  <connection>
    <map_components component_2="fADPext" component_1="mADPext"/>
    <map_variables variable_2="fADPext" variable_1="fADPext"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="mADPext"/>
    <map_variables variable_2="ADPext" variable_1="ADPext"/>
  </connection>
  
  <connection>
    <map_components component_2="ADPx" component_1="fADPx"/>
    <map_variables variable_2="ADPx" variable_1="ADPx"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg_buffering" component_1="fADPx"/>
    <map_variables variable_2="Mgx" variable_1="Mgx"/>
  </connection>
  
  <connection>
    <map_components component_2="fADPx" component_1="mADPx"/>
    <map_variables variable_2="fADPx" variable_1="fADPx"/>
  </connection>
  
  <connection>
    <map_components component_2="ADPx" component_1="mADPx"/>
    <map_variables variable_2="ADPx" variable_1="ADPx"/>
  </connection>
 
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <rdf:Seq rdf:about="rdf:#4896477e-09fc-4eea-b86b-b0b7c3bf6f6a">
    <rdf:li rdf:resource="rdf:#00ed4bb7-c6d4-4a9a-910f-99f24952a534"/>
    <rdf:li rdf:resource="rdf:#08f4b652-9a26-43d6-9177-9d3990290e8f"/>
    <rdf:li rdf:resource="rdf:#d9d66ae3-6ee8-48a7-8be8-6f60caf6297a"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#342003d4-6540-4f12-9cba-9529929153d7">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#00ed4bb7-c6d4-4a9a-910f-99f24952a534">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#0a8a3db7-074d-4665-8a12-1317c20e6985"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0a8a3db7-074d-4665-8a12-1317c20e6985">
    <vCard:Given>Marko</vCard:Given>
    <vCard:Family>Vendelin</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher/>
    <cmeta:comment rdf:resource="rdf:#a2fbb275-0f36-4bf3-b8fe-fe3e17442f7a"/>
    <dcterms:created rdf:resource="rdf:#13d51ce0-2810-4166-aaf7-b0379de682c5"/>
    <dc:creator rdf:resource="rdf:#2dd4e21b-bdf6-47c5-84cb-db8c9233e3a9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aa02c3e7-beee-485b-a6b4-55e16bdda695">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e2a1cc7d-6556-4eb9-94da-388db21d5bdf">
    <dc:creator rdf:resource="rdf:#7b732ec9-829c-47da-b5e4-743b45cc9f33"/>
    <rdf:value>This file contains a CellML description of Vendelin et al.'s mathematical model of mitochondrial respiration in heart cells.&#13;
&#13;
</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#08f4b652-9a26-43d6-9177-9d3990290e8f">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#3056cc03-9ff0-4cee-a2d9-e85d5311c40e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3fb7acbf-98de-47e8-a33c-8cb31db23841">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6c7e1c92-4a2b-4138-a1fa-374874ec33fc">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#54181a00-667e-48ec-b2b4-2edcd2f14270">
    <dc:creator rdf:resource="rdf:#4896477e-09fc-4eea-b86b-b0b7c3bf6f6a"/>
    <dc:title>Regulation of mitochondrial respiration in heart cells analyzed by reaction-diffusion model of energy transfer</dc:title>
    <bqs:volume>278</bqs:volume>
    <bqs:first_page>C747</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#5f04142f-265c-4f37-bc64-c0f5c8cef6b0"/>
    <dcterms:issued rdf:resource="rdf:#02d70ae3-8b64-452c-88e0-e87cbf8f80f5"/>
    <bqs:last_page>C764</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3056cc03-9ff0-4cee-a2d9-e85d5311c40e">
    <vCard:Given>Olav</vCard:Given>
    <vCard:Family>Kongas</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2dd4e21b-bdf6-47c5-84cb-db8c9233e3a9">
    <vCard:ORG rdf:resource="rdf:#342003d4-6540-4f12-9cba-9529929153d7"/>
    <vCard:EMAIL rdf:resource="rdf:#3fb7acbf-98de-47e8-a33c-8cb31db23841"/>
    <vCard:N rdf:resource="rdf:#aa02c3e7-beee-485b-a6b4-55e16bdda695"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#13d51ce0-2810-4166-aaf7-b0379de682c5">
    <dcterms:W3CDTF>2003-11-06T00:00:00+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d9d66ae3-6ee8-48a7-8be8-6f60caf6297a">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#cae93553-7cfd-4226-aa1b-167340e7b6f7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#02d70ae3-8b64-452c-88e0-e87cbf8f80f5">
    <dcterms:W3CDTF>2000-00-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#vendelin_model_2000">
    <bqs:reference rdf:resource="rdf:#8906ea57-8fdd-4832-a5f0-59c845ab6426"/>
    <cmeta:comment rdf:resource="rdf:#e2a1cc7d-6556-4eb9-94da-388db21d5bdf"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cae93553-7cfd-4226-aa1b-167340e7b6f7">
    <vCard:Given>Valdur</vCard:Given>
    <vCard:Family>Saks</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7b732ec9-829c-47da-b5e4-743b45cc9f33">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8906ea57-8fdd-4832-a5f0-59c845ab6426">
    <bqs:Pubmed_id>10751324</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#54181a00-667e-48ec-b2b4-2edcd2f14270"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5f04142f-265c-4f37-bc64-c0f5c8cef6b0">
    <dc:title>American Journal of Physiology Cell Physiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a2fbb275-0f36-4bf3-b8fe-fe3e17442f7a">
    <dc:creator rdf:resource="rdf:#6c7e1c92-4a2b-4138-a1fa-374874ec33fc"/>
    <rdf:value>This CellML model has been checked in both PCEnv and COR.  However the model contains partial differential equations which are currently beyond the scope of the tools.&#13;
&#13;
</rdf:value>
  </rdf:Description>
</rdf:RDF>
</model>