Location: Vendelin, Kongas, Saks, 2000 @ 54152205044e / vendelin_kongas_saks_2000.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2010-08-07 04:34:59+12:00
Desc:
Updated documentation and fixed broken annotations.
Permanent Source URI:
https://staging.physiomeproject.org/workspace/vendelin_kongas_saks_2000/rawfile/54152205044e2616b9016b3f9f606fa9fb7104dd/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>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
        This CellML model has been checked in both OpenCell and COR.  However the model contains partial differential equations which are currently beyond the scope of the CellML API.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: The purpose of this study is to investigate theoretically which intracellular factors may be important for regulation of mitochondrial respiration in working heart cells in vivo. We have developed a model that describes quantitatively the published experimental data on dependence of the rate of oxygen consumption and metabolic state of working isolated perfused rat heart on workload over its physiological range (Williamson JR, Ford G, Illingworth J, Safer B. Circ Res 38, Suppl I, I39-I51, 1976). Analysis of this model shows that for phosphocreatine, creatine, and ATP the equilibrium assumption is an acceptable approximation with respect to their diffusion in the intracellular bulk water phase. However, the ADP concentration changes in the contraction cycle in a nonequilibrium workload-dependent manner, showing the existence of the intracellular concentration gradients. The model shows that workload-dependent alteration of ADP concentration in the compartmentalized creatine kinase system may be taken, together with the changes in P(i) concentration, to be among the major components of the metabolic feedback signal for regulation of respiration in muscle cells.
</para>

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

<para>
Regulation of mitochondrial respiration in heart cells analyzed by reaction-diffusion model of energy transfer, Marko Vendelin, Olav Kongas, and Valdur Saks, 2000,             <emphasis>American Journal of Physiology - Cell Physiology</emphasis>, 278, C747-C764. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/10751324">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.</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"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>electrophysiology</rdf:li>
        <rdf:li>cardiac</rdf:li>
        <rdf:li>mitochondria</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <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.</rdf:value>
  </rdf:Description>
</rdf:RDF>
</model>