Location: Beard, 2005 @ 43b415148771 / beard_2005.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2009-11-27 14:34:18+13:00
Desc:
Added images in ai and svg format, removed non pub med references
Permanent Source URI:
https://models.cellml.org/workspace/beard_2005/rawfile/43b415148771ba435cc838ef1332acb9123f9af7/beard_2005.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 17/07/2008 at 1:28:07 at p.m. using:

COR (0.9.31.979)
Copyright 2002-2008 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk

CellML 1.0 was used to generate this model
http://www.CellML.org/
-->
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="beard_model_2005" name="beard_model_2005_version_2">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>A Biophysical Model of the Mitochondrial Respiratory System and Oxidative Phosphorylation, Beard, 2005</title>
				<author>
					<firstname>David</firstname>
					<surname>Nickerson</surname>
					<affiliation>
						<shortaffil>Bioengineering Institute, The University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
            This model has been edited to correct dimensional errors and checked for other problems. It runs in COR and OpenCell but whether it replicates experimental results is yet to be determined. ValidateCellML confirms this model as valid CellML.
          </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>A computational model for the mitochondrial respiratory chain that appropriately balances mass, charge, and free energy transduction is introduced and analyzed based on a previously published set of data measured on isolated cardiac mitochondria. The basic components included in the model are the reactions at complexes I, III, and IV of the electron transport system, ATP synthesis at F1F0 ATPase, substrate transporters including adenine nucleotide translocase and the phosphate\u2013hydrogen co-transporter, and cation fluxes across the inner membrane including fluxes through the K+/H+ antiporter and passive H+ and K+ permeation. Estimation of 16 adjustable parameter values is based on fitting model simulations to nine independent data curves. The identified model is further validated by comparison to additional datasets measured from mitochondria isolated from rat heart and liver and observed at low oxygen concentration. To obtain reasonable fits to the available data, it is necessary to incorporate inorganic-phosphate-dependent activation of the dehydrogenase activity and the electron transport system. Specifically, it is shown that a model incorporating phosphate-dependent activation of complex III is able to reasonably reproduce the observed data. The resulting validated and verified model provides a foundation for building larger and more complex systems models and investigating complex physiological and pathophysiological interactions in cardiac energetics.    
</para>
				<para>
The complete original paper reference is cited below:
</para>
				<para>
A Biophysical Model of the Mitochondrial Respiratory System and Oxidative Phosphorylation, Daniel A. Beard, 2005, <<emphasis>PLoS Computational Biology</emphasis>, September 9, 1(4). <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=16163394&amp;dopt=Citation">PubMed ID: 16163394</ulink>
				</para>

<para>
Illustration of the Components Included in the Model of Mitochondrial Oxidative Phosphorylation:
</para>

<informalfigure float="0" id="fig_reaction_diagram1">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="beard_2005a.png"/>
  </imageobject>
</mediaobject>
<caption>A) Schematic diagram of the major components of the electron transport system, which transfers reducing potential from NADH to oxygen, and the F1F0 ATPase, which transduces energy from proton motive force to ATP, are illustrated. Complexes I, III, and IV are labeled C1, C3, and C4, respectively.
</caption>
</informalfigure>

<informalfigure float="0" id="fig_reaction_diagram2">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="beard_2005b.png"/>
  </imageobject>
</mediaobject>
<caption>B) Schematic diagram of the substrate transport process included in the model, including the ANT and PiHt on the inner membrane, and passive permeation of ATP, ADP, AMP, and phosphate across the outer membrane. The AK reaction in the IM space is shown.
</caption>
</informalfigure>

<informalfigure float="0" id="fig_reaction_diagram3">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="beard_2005c.png"/>
  </imageobject>
</mediaobject>
<caption>C) Schematic diagram of the transporters for hydrogen and potassium ions on the inner membrane, including K+/H+ antiporter and passive proton and potassium fluxes. It is assumed that these cations rapidly equilibrate across the outer membrane.
</caption>
</informalfigure>
			</sect1>
		</article>
	</documentation>
	   <units base_units="yes" name="cell"/>
   <units base_units="yes" name="cytoplasm"/>
   <units base_units="yes" name="mitochondria"/>
   <units name="cyto_per_cell">
      <unit units="cytoplasm"/>
      <unit units="cell" exponent="-1"/>
   </units>
   <units name="mito_per_cell">
      <unit units="mitochondria"/>
      <unit units="cell" exponent="-1"/>
   </units>
   <units name="l_cyto_per_l_mito">
      <unit units="cytoplasm"/>
      <unit units="mitochondria" exponent="-1"/>
   </units>
   <units name="l_water_per_l_mito">
      <unit units="mitochondria" exponent="-1"/>
   </units>
   <units name="mg_per_mitochondria">
      <unit units="gram" prefix="milli"/>
      <unit units="mitochondria" exponent="-1"/>
   </units>
   <units name="molar">
      <unit units="mole"/>
      <unit units="liter" exponent="-1"/>
   </units>
   <units name="molar_per_second">
      <unit units="molar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="molar_squared">
      <unit units="molar" exponent="2"/>
   </units>
   <units name="molar_cubed">
      <unit units="molar" exponent="3"/>
   </units>
   <units name="micromolar">
      <unit units="molar" prefix="micro"/>
   </units>
   <units name="per_molar">
      <unit units="molar" exponent="-1"/>
   </units>
   <units name="micron_per_second">
      <unit units="metre" prefix="-6"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_micron">
      <unit units="metre" prefix="-6" exponent="-1"/>
      <unit units="mitochondria" exponent="-1"/>
   </units>
   <units name="kilojoule_per_mole">
      <unit units="joule" prefix="kilo"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="kilojoule_per_mole_per_millivolt">
      <unit units="kilojoule_per_mole"/>
      <unit units="millivolt" exponent="-1"/>
   </units>
   <units name="molar_per_cyto">
      <unit units="mole"/>
      <unit units="liter" exponent="-1"/>
      <unit units="mitochondria" exponent="-1"/>
   </units>
   <units name="molar_per_mito">
      <unit units="mole"/>
      <unit units="liter" exponent="-1"/>
      <unit units="mitochondria" exponent="-1"/>
   </units>
   <units name="mole_per_second_per_l_cyto">
      <unit units="molar_per_cyto"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="mole_per_second_per_l_mito">
      <unit units="molar_per_mito"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="mole_per_second_per_l_mito_per_molar">
      <unit units="mole_per_second_per_l_mito"/>
      <unit units="molar" exponent="-1"/>
   </units>
   <units name="mole_per_second_per_l_mito_per_molar_per_half_molar">
      <unit units="mole_per_second_per_l_mito_per_molar"/>
      <unit units="molar" exponent="-0.5"/>
   </units>
   <units name="mole_per_second_per_l_mito_per_molar_per_molar">
      <unit units="mole_per_second_per_l_mito_per_molar"/>
      <unit units="molar" exponent="-1"/>
   </units>
   <units name="mole_per_l_mito_per_millivolt">
      <unit units="molar_per_cyto"/>
      <unit units="millivolt" exponent="-1"/>
   </units>
   <units name="mole_per_second_per_l_mito_per_molar_per_millivolt">
      <unit units="mole_per_second_per_l_mito_per_molar"/>
      <unit units="millivolt" exponent="-1"/>
   </units>
   <component name="Environment">
      <variable units="second" public_interface="out" name="t"/>
      <variable units="molar" public_interface="out" name="Mg_tot" initial_value="0.005"/>
      <variable units="molar" public_interface="out" name="Pi_e" initial_value="0.001"/>
      <variable units="molar" public_interface="out" name="ADP_e" initial_value="0.0013"/>
   </component>
   <component name="Fixed_parameters">
      <variable units="kilojoule_per_mole" public_interface="out" name="RT" initial_value="2.4734"/>
      <variable units="kilojoule_per_mole_per_millivolt" public_interface="out" name="F" initial_value="0.096485"/>
      <variable units="dimensionless" public_interface="out" name="n_A" initial_value="3"/>
      <variable units="kilojoule_per_mole" public_interface="out" name="dG_C1o" initial_value="-69.37"/>
      <variable units="kilojoule_per_mole" public_interface="out" name="dG_C3o" initial_value="-32.53"/>
      <variable units="kilojoule_per_mole" public_interface="out" name="dG_C4o" initial_value="-122.94"/>
      <variable units="kilojoule_per_mole" public_interface="out" name="dG_F1o" initial_value="36.03"/>
      <variable units="dimensionless" public_interface="out" name="pH_e" initial_value="7.1"/>
      <variable units="molar" public_interface="out" name="H_e"/>
      <variable units="molar" public_interface="out" name="K_e" initial_value="0.15"/>
      <variable units="molar" public_interface="out" name="ATP_e" initial_value="0"/>
      <variable units="molar" public_interface="out" name="AMP_e" initial_value="0"/>
      <variable units="molar" public_interface="out" name="k_dHPi" initial_value="1.7782754100389227e-7"/>
      <variable units="molar" public_interface="out" name="k_dHatp" initial_value="3.311311214825911e-7"/>
      <variable units="molar" public_interface="out" name="k_dHadp" initial_value="5.128613839913649e-7"/>
      <variable units="molar" public_interface="out" name="K_DT" initial_value="2.4e-5"/>
      <variable units="molar" public_interface="out" name="K_DD" initial_value="3.47e-4"/>
      <variable units="dimensionless" public_interface="out" name="K_AK" initial_value="0.4331"/>
      <variable units="l_water_per_l_mito" public_interface="out" name="W_m" initial_value="0.72376"/>
      <variable units="l_water_per_l_mito" public_interface="out" name="W_x"/>
      <variable units="l_water_per_l_mito" public_interface="out" name="W_i"/>
      <variable units="per_micron" public_interface="out" name="gamma" initial_value="5.99"/>
      <variable units="molar" public_interface="out" name="Ctot" initial_value="0.0027"/>
      <variable units="molar" public_interface="out" name="Qtot" initial_value="0.00135"/>
      <variable units="molar" public_interface="out" name="NADtot" initial_value="0.00297"/>
      <variable units="molar" public_interface="out" name="H_i"/>
      <variable units="molar" public_interface="out" name="K_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>H_e</ci>
            <apply>
               <times/>
               <cn cellml:units="molar">1</cn>
               <apply>
                  <power/>
                  <cn cellml:units="dimensionless">10</cn>
                  <apply>
                     <minus/>
                     <ci>pH_e</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>W_x</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.9</cn>
               <ci>W_m</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>W_i</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.1</cn>
               <ci>W_m</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>H_i</ci>
            <ci>H_e</ci>
         </apply>
         <apply>
            <eq/>
            <ci>K_i</ci>
            <ci>K_e</ci>
         </apply>
      </math>
   </component>
   <component name="Adjustable_parameters">
      <variable units="molar" public_interface="out" name="k_Pi1" initial_value="1.3413e-4"/>
      <variable units="molar" public_interface="out" name="k_Pi2" initial_value="6.7668e-4"/>
      <variable units="molar" public_interface="out" name="k_Pi3" initial_value="1.9172e-4"/>
      <variable units="molar" public_interface="out" name="k_Pi4" initial_value="0.02531"/>
      <variable units="molar" public_interface="out" name="k_PiH" initial_value="4.5082e-4"/>
      <variable units="dimensionless" public_interface="out" name="r" initial_value="4.5807"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="out" name="x_DH" initial_value="0.09183"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="out" name="x_C1" initial_value="0.36923"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="out" name="x_C3" initial_value="0.091737"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="out" name="x_C4" initial_value="3.2562e-5"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="out" name="x_F1" initial_value="150.93"/>
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="x_ANT" initial_value="0.0079204"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="out" name="x_Pi1" initial_value="339430"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="out" name="x_KH" initial_value="2.9802e7"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_millivolt" public_interface="out" name="x_Hle" initial_value="250"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_millivolt" public_interface="out" name="x_K" initial_value="0"/>
      <variable units="molar" public_interface="out" name="k_mADP" initial_value="3.5e-6"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="out" name="x_AK" initial_value="0"/>
      <variable units="micron_per_second" public_interface="out" name="p_A" initial_value="85"/>
      <variable units="molar" public_interface="out" name="k_O2" initial_value="1.2e-4"/>
      <variable units="dimensionless" public_interface="out" name="r_buff" initial_value="4.253"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="out" name="x_MgA" initial_value="1000000"/>
      <variable units="micron_per_second" public_interface="out" name="x_Pi2" initial_value="327"/>
   </component>
   <component name="Proton_motive_force">
      <variable units="kilojoule_per_mole" public_interface="out" name="dG_H"/>
      <variable units="kilojoule_per_mole_per_millivolt" public_interface="in" name="F"/>
      <variable units="millivolt" public_interface="in" name="dPsi"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="RT"/>
      <variable units="molar" public_interface="in" name="H_i"/>
      <variable units="molar" public_interface="in" name="H_x"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>dG_H</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>F</ci>
                  <ci>dPsi</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>RT</ci>
                  <apply>
                     <ln/>
                     <apply>
                        <divide/>
                        <ci>H_i</ci>
                        <ci>H_x</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Dehydrogenase_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_DH"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="in" name="x_DH"/>
      <variable units="molar" public_interface="in" name="k_Pi1"/>
      <variable units="molar" public_interface="in" name="k_Pi2"/>
      <variable units="dimensionless" public_interface="in" name="r"/>
      <variable units="molar" public_interface="in" name="NAD_x"/>
      <variable units="molar" public_interface="in" name="NADH_x"/>
      <variable units="molar" public_interface="in" name="Pi_x"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_DH</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>x_DH</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>r</ci>
                        <ci>NAD_x</ci>
                     </apply>
                     <ci>NADH_x</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Pi_x</ci>
                        <ci>k_Pi1</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>Pi_x</ci>
                     <ci>k_Pi2</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Electron_flux_complex_I">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_C1"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="in" name="x_C1"/>
      <variable units="kilojoule_per_mole" name="dG_C1op"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="dG_C1o"/>
      <variable units="molar" public_interface="in" name="NADH_x"/>
      <variable units="molar" public_interface="in" name="NAD_x"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="dG_H"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="RT"/>
      <variable units="molar" public_interface="in" name="Q"/>
      <variable units="molar" public_interface="in" name="QH2"/>
      <variable units="molar" public_interface="in" name="H_x"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>dG_C1op</ci>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <ci>dG_C1o</ci>
                  <apply>
                     <times/>
                     <ci>RT</ci>
                     <apply>
                        <ln/>
                        <apply>
                           <divide/>
                           <ci>H_x</ci>
                           <cn cellml:units="molar" type="e-notation">1<sep/>-7</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>RT</ci>
                  <apply>
                     <ln/>
                     <apply>
                        <divide/>
                        <ci>Q</ci>
                        <ci>QH2</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_C1</ci>
            <apply>
               <times/>
               <ci>x_C1</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>dG_C1op</ci>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">4</cn>
                                    <ci>dG_H</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <ci>RT</ci>
                        </apply>
                     </apply>
                     <ci>NADH_x</ci>
                  </apply>
                  <ci>NAD_x</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Electron_flux_complex_III">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_C3"/>
      <variable units="molar" public_interface="in" name="k_Pi3"/>
      <variable units="molar" public_interface="in" name="k_Pi4"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="in" name="x_C3"/>
      <variable units="kilojoule_per_mole" name="dG_C3op"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="dG_C3o"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="dG_H"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="RT"/>
      <variable units="kilojoule_per_mole_per_millivolt" public_interface="in" name="F"/>
      <variable units="millivolt" public_interface="in" name="dPsi"/>
      <variable units="molar" public_interface="in" name="Cox"/>
      <variable units="molar" public_interface="in" name="Cred"/>
      <variable units="molar" public_interface="in" name="Q"/>
      <variable units="molar" public_interface="in" name="QH2"/>
      <variable units="molar" public_interface="in" name="H_x"/>
      <variable units="molar" public_interface="in" name="Pi_x"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>dG_C3op</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>dG_C3o</ci>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>RT</ci>
                     <apply>
                        <ln/>
                        <apply>
                           <divide/>
                           <ci>H_x</ci>
                           <cn cellml:units="molar" type="e-notation">1<sep/>-7</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>RT</ci>
                  <apply>
                     <ln/>
                     <apply>
                        <divide/>
                        <ci>QH2</ci>
                        <ci>Q</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_C3</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>x_C3</ci>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <divide/>
                           <ci>Pi_x</ci>
                           <ci>k_Pi3</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <apply>
                        <divide/>
                        <ci>Pi_x</ci>
                        <ci>k_Pi4</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <plus/>
                                    <ci>dG_C3op</ci>
                                    <apply>
                                       <times/>
                                       <cn cellml:units="dimensionless">4</cn>
                                       <ci>dG_H</ci>
                                    </apply>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">2</cn>
                                    <ci>F</ci>
                                    <ci>dPsi</ci>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">2</cn>
                              <ci>RT</ci>
                           </apply>
                        </apply>
                     </apply>
                     <ci>Cox</ci>
                  </apply>
                  <ci>Cred</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Electron_flux_complex_IV">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_C4"/>
      <variable units="molar" public_interface="in" name="k_O2"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="in" name="x_C4"/>
      <variable units="kilojoule_per_mole" name="dG_C4op"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="dG_C4o"/>
      <variable units="molar" public_interface="in" name="Cred"/>
      <variable units="molar" public_interface="in" name="Ctot"/>
      <variable units="molar" public_interface="in" name="O2"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="dG_H"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="RT"/>
      <variable units="molar" public_interface="in" name="Cox"/>
      <variable units="kilojoule_per_mole_per_millivolt" public_interface="in" name="F"/>
      <variable units="millivolt" public_interface="in" name="dPsi"/>
      <variable units="molar" public_interface="in" name="H_x"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>dG_C4op</ci>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <ci>dG_C4o</ci>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>RT</ci>
                     <apply>
                        <ln/>
                        <apply>
                           <divide/>
                           <ci>H_x</ci>
                           <cn cellml:units="molar" type="e-notation">1<sep/>-7</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>RT</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <ln/>
                     <apply>
                        <divide/>
                        <ci>O2</ci>
                        <cn cellml:units="molar">1</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_C4</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <ci>x_C4</ci>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <ci>k_O2</ci>
                              <ci>O2</ci>
                           </apply>
                        </apply>
                     </apply>
                     <ci>Cred</ci>
                  </apply>
                  <ci>Ctot</ci>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <ci>dG_C4op</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>dG_H</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">2</cn>
                              <ci>RT</ci>
                           </apply>
                        </apply>
                     </apply>
                     <ci>Cred</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Cox</ci>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>F</ci>
                              <ci>dPsi</ci>
                           </apply>
                           <ci>RT</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ATP_synthesis_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_F1"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="in" name="x_F1"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="dG_F1o"/>
      <variable units="molar" public_interface="in" name="K_DD"/>
      <variable units="molar" public_interface="in" name="K_DT"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="RT"/>
      <variable units="dimensionless" public_interface="in" name="n_A"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="dG_H"/>
      <variable units="molar" public_interface="in" name="ADP_mx"/>
      <variable units="molar" public_interface="in" name="Pi_x"/>
      <variable units="molar" public_interface="in" name="ATP_mx"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_F1</ci>
            <apply>
               <times/>
               <ci>x_F1</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <exp/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <minus/>
                                    <apply>
                                       <minus/>
                                       <ci>dG_F1o</ci>
                                       <apply>
                                          <times/>
                                          <ci>n_A</ci>
                                          <ci>dG_H</ci>
                                       </apply>
                                    </apply>
                                 </apply>
                                 <ci>RT</ci>
                              </apply>
                           </apply>
                           <ci>K_DD</ci>
                        </apply>
                        <ci>K_DT</ci>
                     </apply>
                     <ci>ADP_mx</ci>
                     <ci>Pi_x</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>ATP_mx</ci>
                     <cn cellml:units="molar">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ANT_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_ANT"/>
      <variable units="millivolt" public_interface="out" name="Psi_x"/>
      <variable units="millivolt" public_interface="out" name="Psi_i"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="x_ANT"/>
      <variable units="molar" public_interface="in" name="ADP_fi"/>
      <variable units="molar" public_interface="in" name="ATP_fi"/>
      <variable units="kilojoule_per_mole_per_millivolt" public_interface="in" name="F"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="RT"/>
      <variable units="millivolt" public_interface="in" name="dPsi"/>
      <variable units="molar" public_interface="in" name="ADP_fx"/>
      <variable units="molar" public_interface="in" name="ATP_fx"/>
      <variable units="molar" public_interface="in" name="k_mADP"/>
      <variable units="molar" name="mincond" initial_value="1e-12"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Psi_x</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">0.65</cn>
               </apply>
               <ci>dPsi</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Psi_i</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">0.35</cn>
               <ci>dPsi</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_ANT</ci>
            <piecewise>
               <piece>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>x_ANT</ci>
                        <apply>
                           <minus/>
                           <apply>
                              <divide/>
                              <ci>ADP_fi</ci>
                              <apply>
                                 <plus/>
                                 <ci>ADP_fi</ci>
                                 <apply>
                                    <times/>
                                    <ci>ATP_fi</ci>
                                    <apply>
                                       <exp/>
                                       <apply>
                                          <divide/>
                                          <apply>
                                             <times/>
                                             <apply>
                                                <minus/>
                                                <ci>F</ci>
                                             </apply>
                                             <ci>Psi_i</ci>
                                          </apply>
                                          <ci>RT</ci>
                                       </apply>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                           <apply>
                              <divide/>
                              <ci>ADP_fx</ci>
                              <apply>
                                 <plus/>
                                 <ci>ADP_fx</ci>
                                 <apply>
                                    <times/>
                                    <ci>ATP_fx</ci>
                                    <apply>
                                       <exp/>
                                       <apply>
                                          <divide/>
                                          <apply>
                                             <times/>
                                             <apply>
                                                <minus/>
                                                <ci>F</ci>
                                             </apply>
                                             <ci>Psi_x</ci>
                                          </apply>
                                          <ci>RT</ci>
                                       </apply>
                                    </apply>
                                 </apply>
                              </apply>
                           </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <divide/>
                           <ci>k_mADP</ci>
                           <ci>ADP_fi</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <or/>
                     <apply>
                        <gt/>
                        <ci>ADP_fi</ci>
                        <ci>mincond</ci>
                     </apply>
                     <apply>
                        <gt/>
                        <ci>ATP_fi</ci>
                        <ci>mincond</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="mole_per_second_per_l_mito">0</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="MgATPx_binding_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_MgATPx"/>
      <variable units="molar" public_interface="out" name="ATP_fx"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="in" name="x_MgA"/>
      <variable units="molar" public_interface="in" name="ATP_x"/>
      <variable units="molar" public_interface="in" name="ATP_mx"/>
      <variable units="molar" public_interface="in" name="Mg_x"/>
      <variable units="molar" public_interface="in" name="K_DT"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>ATP_fx</ci>
            <apply>
               <minus/>
               <ci>ATP_x</ci>
               <ci>ATP_mx</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_MgATPx</ci>
            <apply>
               <times/>
               <ci>x_MgA</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>ATP_fx</ci>
                     <ci>Mg_x</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_DT</ci>
                     <ci>ATP_mx</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="MgADPx_binding_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_MgADPx"/>
      <variable units="molar" public_interface="out" name="ADP_fx"/>
      <variable units="molar" public_interface="in" name="ADP_x"/>
      <variable units="molar" public_interface="in" name="ADP_mx"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="in" name="x_MgA"/>
      <variable units="molar" public_interface="in" name="K_DD"/>
      <variable units="molar" public_interface="in" name="Mg_x"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>ADP_fx</ci>
            <apply>
               <minus/>
               <ci>ADP_x</ci>
               <ci>ADP_mx</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_MgADPx</ci>
            <apply>
               <times/>
               <ci>x_MgA</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>ADP_fx</ci>
                     <ci>Mg_x</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_DD</ci>
                     <ci>ADP_mx</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="MgATPi_binding_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_MgATPi"/>
      <variable units="molar" public_interface="out" name="ATP_fi"/>
      <variable units="molar" public_interface="in" name="ATP_i"/>
      <variable units="molar" public_interface="in" name="ATP_mi"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="in" name="x_MgA"/>
      <variable units="molar" public_interface="in" name="K_DT"/>
      <variable units="molar" public_interface="in" name="Mg_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>ATP_fi</ci>
            <apply>
               <minus/>
               <ci>ATP_i</ci>
               <ci>ATP_mi</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_MgATPi</ci>
            <apply>
               <times/>
               <ci>x_MgA</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>ATP_fi</ci>
                     <ci>Mg_i</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_DT</ci>
                     <ci>ATP_mi</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="MgADPi_binding_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_MgADPi"/>
      <variable units="molar" public_interface="out" name="ADP_fi"/>
      <variable units="molar" public_interface="in" name="ADP_i"/>
      <variable units="molar" public_interface="in" name="ADP_mi"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="in" name="x_MgA"/>
      <variable units="molar" public_interface="in" name="K_DD"/>
      <variable units="molar" public_interface="in" name="Mg_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>ADP_fi</ci>
            <apply>
               <minus/>
               <ci>ADP_i</ci>
               <ci>ADP_mi</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_MgADPi</ci>
            <apply>
               <times/>
               <ci>x_MgA</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>ADP_fi</ci>
                     <ci>Mg_i</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_DD</ci>
                     <ci>ADP_mi</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ATP_substrate_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_ATP"/>
      <variable units="molar" public_interface="in" name="ATP_e"/>
      <variable units="per_micron" public_interface="in" name="gamma"/>
      <variable units="micron_per_second" public_interface="in" name="p_A"/>
      <variable units="molar" public_interface="in" name="ATP_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_ATP</ci>
            <apply>
               <times/>
               <ci>gamma</ci>
               <ci>p_A</ci>
               <apply>
                  <minus/>
                  <ci>ATP_e</ci>
                  <ci>ATP_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ADP_substrate_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_ADP"/>
      <variable units="per_micron" public_interface="in" name="gamma"/>
      <variable units="micron_per_second" public_interface="in" name="p_A"/>
      <variable units="molar" public_interface="in" name="ADP_e"/>
      <variable units="molar" public_interface="in" name="ADP_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_ADP</ci>
            <apply>
               <times/>
               <ci>gamma</ci>
               <ci>p_A</ci>
               <apply>
                  <minus/>
                  <ci>ADP_e</ci>
                  <ci>ADP_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="AMP_substrate_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_AMP"/>
      <variable units="molar" public_interface="in" name="AMP_e"/>
      <variable units="per_micron" public_interface="in" name="gamma"/>
      <variable units="micron_per_second" public_interface="in" name="p_A"/>
      <variable units="molar" public_interface="in" name="AMP_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_AMP</ci>
            <apply>
               <times/>
               <ci>gamma</ci>
               <ci>p_A</ci>
               <apply>
                  <minus/>
                  <ci>AMP_e</ci>
                  <ci>AMP_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Pi_substrate_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_Pi2"/>
      <variable units="micron_per_second" public_interface="in" name="x_Pi2"/>
      <variable units="per_micron" public_interface="in" name="gamma"/>
      <variable units="molar" public_interface="in" name="Pi_e"/>
      <variable units="molar" public_interface="in" name="Pi_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_Pi2</ci>
            <apply>
               <times/>
               <ci>gamma</ci>
               <ci>x_Pi2</ci>
               <apply>
                  <minus/>
                  <ci>Pi_e</ci>
                  <ci>Pi_i</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Phosphate_hydrogen_cotransporter_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_Pi1"/>
      <variable units="molar" public_interface="out" name="H2PIi"/>
      <variable units="molar" public_interface="out" name="H2PIx"/>
      <variable units="molar" public_interface="in" name="Pi_i"/>
      <variable units="molar" public_interface="in" name="H_i"/>
      <variable units="molar" public_interface="in" name="k_dHPi"/>
      <variable units="molar" public_interface="in" name="Pi_x"/>
      <variable units="molar" public_interface="in" name="H_x"/>
      <variable units="mole_per_second_per_l_mito_per_molar" public_interface="in" name="x_Pi1"/>
      <variable units="molar" public_interface="in" name="k_PiH"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>H2PIi</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Pi_i</ci>
                  <ci>H_i</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>H_i</ci>
                  <ci>k_dHPi</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>H2PIx</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Pi_x</ci>
                  <ci>H_x</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>H_x</ci>
                  <ci>k_dHPi</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_Pi1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>x_Pi1</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>H_x</ci>
                        <ci>H2PIi</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>H_i</ci>
                        <ci>H2PIx</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <ci>H2PIi</ci>
                  <ci>k_PiH</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Adenylate_kinase_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_AKi"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="in" name="x_AK"/>
      <variable units="dimensionless" public_interface="in" name="K_AK"/>
      <variable units="molar" public_interface="in" name="ADP_i"/>
      <variable units="molar" public_interface="in" name="AMP_i"/>
      <variable units="molar" public_interface="in" name="ATP_i"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_AKi</ci>
            <apply>
               <times/>
               <ci>x_AK</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>K_AK</ci>
                     <ci>ADP_i</ci>
                     <ci>ADP_i</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>AMP_i</ci>
                     <ci>ATP_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Hydrogen_leak_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_Hle"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_millivolt" public_interface="in" name="x_Hle"/>
      <variable units="millivolt" public_interface="in" name="dPsi"/>
      <variable units="molar" public_interface="in" name="H_i"/>
      <variable units="kilojoule_per_mole_per_millivolt" public_interface="in" name="F"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="RT"/>
      <variable units="molar" public_interface="in" name="H_x"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_Hle</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>x_Hle</ci>
                  <ci>dPsi</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>H_i</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>F</ci>
                                 <ci>dPsi</ci>
                              </apply>
                              <ci>RT</ci>
                           </apply>
                        </apply>
                     </apply>
                     <ci>H_x</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>F</ci>
                           <ci>dPsi</ci>
                        </apply>
                        <ci>RT</ci>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Passive_potassium_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_K"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_millivolt" public_interface="in" name="x_K"/>
      <variable units="millivolt" public_interface="in" name="dPsi"/>
      <variable units="molar" public_interface="in" name="K_i"/>
      <variable units="kilojoule_per_mole_per_millivolt" public_interface="in" name="F"/>
      <variable units="kilojoule_per_mole" public_interface="in" name="RT"/>
      <variable units="molar" public_interface="in" name="K_x"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_K</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>x_K</ci>
                  <ci>dPsi</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>K_i</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <ci>F</ci>
                                 <ci>dPsi</ci>
                              </apply>
                              <ci>RT</ci>
                           </apply>
                        </apply>
                     </apply>
                     <ci>K_x</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>F</ci>
                           <ci>dPsi</ci>
                        </apply>
                        <ci>RT</ci>
                     </apply>
                  </apply>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Potassium_hydrogen_flux">
      <variable units="mole_per_second_per_l_mito" public_interface="out" name="J_KH"/>
      <variable units="mole_per_second_per_l_mito_per_molar_per_molar" public_interface="in" name="x_KH"/>
      <variable units="molar" public_interface="in" name="K_i"/>
      <variable units="molar" public_interface="in" name="K_x"/>
      <variable units="molar" public_interface="in" name="H_i"/>
      <variable units="molar" public_interface="in" name="H_x"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_KH</ci>
            <apply>
               <times/>
               <ci>x_KH</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>K_i</ci>
                     <ci>H_x</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>K_x</ci>
                     <ci>H_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="NAD_x_concentration">
      <variable units="molar" public_interface="out" name="NAD_x"/>
      <variable units="molar" public_interface="in" name="NADH_x"/>
      <variable units="molar" public_interface="in" name="NADtot"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>NAD_x</ci>
            <apply>
               <minus/>
               <ci>NADtot</ci>
               <ci>NADH_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Q_concentration">
      <variable units="molar" public_interface="out" name="Q"/>
      <variable units="molar" public_interface="in" name="Qtot"/>
      <variable units="molar" public_interface="in" name="QH2"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Q</ci>
            <apply>
               <minus/>
               <ci>Qtot</ci>
               <ci>QH2</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Cox_concentration">
      <variable units="molar" public_interface="out" name="Cox"/>
      <variable units="molar" public_interface="in" name="Ctot"/>
      <variable units="molar" public_interface="in" name="Cred"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Cox</ci>
            <apply>
               <minus/>
               <ci>Ctot</ci>
               <ci>Cred</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ADP_binding">
      <variable units="molar" public_interface="out" name="ADP_me"/>
      <variable units="molar" public_interface="out" name="ADP_fe"/>
      <variable units="molar" public_interface="in" name="ADP_e"/>
      <variable units="molar" public_interface="in" name="K_DD"/>
      <variable units="molar" public_interface="in" name="Mg_tot"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>ADP_fe</ci>
            <apply>
               <minus/>
               <ci>ADP_e</ci>
               <ci>ADP_me</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ADP_me</ci>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>K_DD</ci>
                     <ci>ADP_e</ci>
                     <ci>Mg_tot</ci>
                  </apply>
                  <apply>
                     <power/>
                     <apply>
                        <minus/>
                        <apply>
                           <power/>
                           <apply>
                              <plus/>
                              <ci>K_DD</ci>
                              <ci>ADP_e</ci>
                              <ci>Mg_tot</ci>
                           </apply>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">4</cn>
                           <ci>Mg_tot</ci>
                           <ci>ADP_e</ci>
                        </apply>
                     </apply>
                     <cn cellml:units="dimensionless">0.5</cn>
                  </apply>
               </apply>
               <cn cellml:units="dimensionless">2</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Mg_binding">
      <variable units="molar" public_interface="out" name="Mg_e"/>
      <variable units="molar" public_interface="in" name="Mg_tot"/>
      <variable units="molar" public_interface="out" name="Mg_i"/>
      <variable units="molar" public_interface="in" name="ADP_me"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Mg_e</ci>
            <apply>
               <minus/>
               <ci>Mg_tot</ci>
               <ci>ADP_me</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Mg_i</ci>
            <ci>Mg_e</ci>
         </apply>
      </math>
   </component>
   <component name="dH_x_dt">
      <variable units="molar" public_interface="out" name="H_x" initial_value="6.309573444801932e-8"/>
      <variable units="dimensionless" public_interface="in" name="r_buff"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_DH"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C3"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C4"/>
      <variable units="dimensionless" public_interface="in" name="n_A"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_F1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_Pi1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_Hle"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_KH"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>H_x</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>r_buff</ci>
                  </apply>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <apply>
                                 <minus/>
                                 <apply>
                                    <minus/>
                                    <ci>J_DH</ci>
                                 </apply>
                                 <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">3</cn>
                                    <ci>J_C1</ci>
                                 </apply>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">4</cn>
                                 <ci>J_C3</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <cn cellml:units="dimensionless">2</cn>
                              <ci>J_C4</ci>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>n_A</ci>
                              <cn cellml:units="dimensionless">1</cn>
                           </apply>
                           <ci>J_F1</ci>
                        </apply>
                        <apply>
                           <times/>
                           <cn cellml:units="dimensionless">2</cn>
                           <ci>J_Pi1</ci>
                        </apply>
                        <ci>J_Hle</ci>
                     </apply>
                     <ci>J_KH</ci>
                  </apply>
               </apply>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dK_x_dt">
      <variable units="molar" public_interface="out" name="K_x" initial_value="0.14"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_KH"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_K"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>K_x</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <ci>J_KH</ci>
                  <ci>J_K</ci>
               </apply>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dMg_x_dt">
      <variable units="molar" public_interface="out" name="Mg_x" initial_value="0.005"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_MgATPx"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_MgADPx"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>Mg_x</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <ci>J_MgATPx</ci>
                  </apply>
                  <ci>J_MgADPx</ci>
               </apply>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dNADH_x_dt">
      <variable units="molar" public_interface="out" name="NADH_x" initial_value="0.0015"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_DH"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C1"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>NADH_x</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>J_DH</ci>
                  <ci>J_C1</ci>
               </apply>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dQH2_dt">
      <variable units="molar" public_interface="out" name="QH2" initial_value="8e-4"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C3"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>QH2</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>J_C1</ci>
                  <ci>J_C3</ci>
               </apply>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dCred_dt">
      <variable units="molar" public_interface="out" name="Cred" initial_value="0.001"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C3"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C4"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_i"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>Cred</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>J_C3</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>J_C4</ci>
                  </apply>
               </apply>
               <ci>W_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dATP_x_dt">
      <variable units="molar" public_interface="out" name="ATP_x" initial_value="0"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_F1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_ANT"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>ATP_x</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>J_F1</ci>
                  <ci>J_ANT</ci>
               </apply>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dADP_x_dt">
      <variable units="molar" public_interface="out" name="ADP_x" initial_value="0.01"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_F1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_ANT"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>ADP_x</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <ci>J_F1</ci>
                  </apply>
                  <ci>J_ANT</ci>
               </apply>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dATP_mx_dt">
      <variable units="molar" public_interface="out" name="ATP_mx" initial_value="0"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_MgATPx"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>ATP_mx</ci>
            </apply>
            <apply>
               <divide/>
               <ci>J_MgATPx</ci>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dADP_mx_dt">
      <variable units="molar" public_interface="out" name="ADP_mx" initial_value="0"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_MgADPx"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>ADP_mx</ci>
            </apply>
            <apply>
               <divide/>
               <ci>J_MgADPx</ci>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dPi_x_dt">
      <variable units="molar" public_interface="out" name="Pi_x" initial_value="0.001"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_F1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_Pi1"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_x"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>Pi_x</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <ci>J_F1</ci>
                  </apply>
                  <ci>J_Pi1</ci>
               </apply>
               <ci>W_x</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dATP_i_dt">
      <variable units="molar" public_interface="out" name="ATP_i" initial_value="0"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_ATP"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_ANT"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_AKi"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_i"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>ATP_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <ci>J_ATP</ci>
                  <ci>J_ANT</ci>
                  <ci>J_AKi</ci>
               </apply>
               <ci>W_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dADP_i_dt">
      <variable units="molar" public_interface="out" name="ADP_i" initial_value="0"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_ADP"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_ANT"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_AKi"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_i"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>ADP_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <ci>J_ADP</ci>
                     <ci>J_ANT</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>J_AKi</ci>
                  </apply>
               </apply>
               <ci>W_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dAMP_i_dt">
      <variable units="molar" public_interface="out" name="AMP_i" initial_value="0"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_AMP"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_AKi"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_i"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>AMP_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <ci>J_AMP</ci>
                  <ci>J_AKi</ci>
               </apply>
               <ci>W_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dATP_mi_dt">
      <variable units="molar" public_interface="out" name="ATP_mi" initial_value="0"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_MgATPi"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_i"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>ATP_mi</ci>
            </apply>
            <apply>
               <divide/>
               <ci>J_MgATPi</ci>
               <ci>W_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dADP_mi_dt">
      <variable units="molar" public_interface="out" name="ADP_mi" initial_value="0"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_MgADPi"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_i"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>ADP_mi</ci>
            </apply>
            <apply>
               <divide/>
               <ci>J_MgADPi</ci>
               <ci>W_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dPi_i_dt">
      <variable units="molar" public_interface="out" name="Pi_i" initial_value="0.001"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_Pi1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_Pi2"/>
      <variable units="l_water_per_l_mito" public_interface="in" name="W_i"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>Pi_i</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <ci>J_Pi1</ci>
                  </apply>
                  <ci>J_Pi2</ci>
               </apply>
               <ci>W_i</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dPsi_dt">
      <variable units="millivolt" public_interface="out" name="dPsi" initial_value="160"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C3"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_C4"/>
      <variable units="dimensionless" public_interface="in" name="n_A"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_F1"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_ANT"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_Hle"/>
      <variable units="mole_per_second_per_l_mito" public_interface="in" name="J_K"/>
      <variable units="second" public_interface="in" name="t"/>
      <variable units="mole_per_l_mito_per_millivolt" name="C_im" initial_value="1e-6"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>dPsi</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">4</cn>
                                 <ci>J_C1</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">2</cn>
                                 <ci>J_C3</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <cn cellml:units="dimensionless">4</cn>
                                 <ci>J_C4</ci>
                              </apply>
                           </apply>
                           <apply>
                              <times/>
                              <ci>n_A</ci>
                              <ci>J_F1</ci>
                           </apply>
                        </apply>
                        <ci>J_ANT</ci>
                     </apply>
                     <ci>J_Hle</ci>
                  </apply>
                  <ci>J_K</ci>
               </apply>
               <ci>C_im</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="dO2_dt">
      <variable units="molar" public_interface="out" name="O2" initial_value="2.6e-5"/>
      <variable units="second" public_interface="in" name="t"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>t</ci>
               </bvar>
               <ci>O2</ci>
            </apply>
            <cn cellml:units="molar_per_second">0</cn>
         </apply>
      </math>
   </component>
   <connection>
      <map_components component_2="MgATPx_binding_flux" component_1="dATP_x_dt"/>
      <map_variables variable_2="ATP_x" variable_1="ATP_x"/>
   </connection>
   <connection>
      <map_components component_2="MgATPx_binding_flux" component_1="dATP_mx_dt"/>
      <map_variables variable_2="ATP_mx" variable_1="ATP_mx"/>
      <map_variables variable_2="J_MgATPx" variable_1="J_MgATPx"/>
   </connection>
   <connection>
      <map_components component_2="MgADPx_binding_flux" component_1="dADP_x_dt"/>
      <map_variables variable_2="ADP_x" variable_1="ADP_x"/>
   </connection>
   <connection>
      <map_components component_2="MgADPx_binding_flux" component_1="dADP_mx_dt"/>
      <map_variables variable_2="ADP_mx" variable_1="ADP_mx"/>
      <map_variables variable_2="J_MgADPx" variable_1="J_MgADPx"/>
   </connection>
   <connection>
      <map_components component_2="MgATPi_binding_flux" component_1="dATP_i_dt"/>
      <map_variables variable_2="ATP_i" variable_1="ATP_i"/>
   </connection>
   <connection>
      <map_components component_2="MgATPi_binding_flux" component_1="dATP_mi_dt"/>
      <map_variables variable_2="ATP_mi" variable_1="ATP_mi"/>
      <map_variables variable_2="J_MgATPi" variable_1="J_MgATPi"/>
   </connection>
   <connection>
      <map_components component_2="MgADPi_binding_flux" component_1="dADP_i_dt"/>
      <map_variables variable_2="ADP_i" variable_1="ADP_i"/>
   </connection>
   <connection>
      <map_components component_2="MgADPi_binding_flux" component_1="dADP_mi_dt"/>
      <map_variables variable_2="ADP_mi" variable_1="ADP_mi"/>
      <map_variables variable_2="J_MgADPi" variable_1="J_MgADPi"/>
   </connection>
   <connection>
      <map_components component_2="ADP_binding" component_1="Fixed_parameters"/>
      <map_variables variable_2="K_DD" variable_1="K_DD"/>
   </connection>
   <connection>
      <map_components component_2="ADP_binding" component_1="Environment"/>
      <map_variables variable_2="Mg_tot" variable_1="Mg_tot"/>
      <map_variables variable_2="ADP_e" variable_1="ADP_e"/>
   </connection>
   <connection>
      <map_components component_2="Mg_binding" component_1="Environment"/>
      <map_variables variable_2="Mg_tot" variable_1="Mg_tot"/>
   </connection>
   <connection>
      <map_components component_2="Mg_binding" component_1="ADP_binding"/>
      <map_variables variable_2="ADP_me" variable_1="ADP_me"/>
   </connection>
   <connection>
      <map_components component_2="Proton_motive_force" component_1="Fixed_parameters"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="RT" variable_1="RT"/>
      <map_variables variable_2="H_i" variable_1="H_i"/>
   </connection>
   <connection>
      <map_components component_2="Proton_motive_force" component_1="dPsi_dt"/>
      <map_variables variable_2="dPsi" variable_1="dPsi"/>
   </connection>
   <connection>
      <map_components component_2="Proton_motive_force" component_1="dH_x_dt"/>
      <map_variables variable_2="H_x" variable_1="H_x"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_I" component_1="Fixed_parameters"/>
      <map_variables variable_2="dG_C1o" variable_1="dG_C1o"/>
      <map_variables variable_2="RT" variable_1="RT"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_I" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_C1" variable_1="x_C1"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_I" component_1="Q_concentration"/>
      <map_variables variable_2="Q" variable_1="Q"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_I" component_1="dQH2_dt"/>
      <map_variables variable_2="QH2" variable_1="QH2"/>
      <map_variables variable_2="J_C1" variable_1="J_C1"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_I" component_1="dH_x_dt"/>
      <map_variables variable_2="H_x" variable_1="H_x"/>
      <map_variables variable_2="J_C1" variable_1="J_C1"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_I" component_1="dNADH_x_dt"/>
      <map_variables variable_2="NADH_x" variable_1="NADH_x"/>
      <map_variables variable_2="J_C1" variable_1="J_C1"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_I" component_1="NAD_x_concentration"/>
      <map_variables variable_2="NAD_x" variable_1="NAD_x"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_I" component_1="Proton_motive_force"/>
      <map_variables variable_2="dG_H" variable_1="dG_H"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_III" component_1="dH_x_dt"/>
      <map_variables variable_2="H_x" variable_1="H_x"/>
      <map_variables variable_2="J_C3" variable_1="J_C3"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_III" component_1="Q_concentration"/>
      <map_variables variable_2="Q" variable_1="Q"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_III" component_1="dQH2_dt"/>
      <map_variables variable_2="QH2" variable_1="QH2"/>
      <map_variables variable_2="J_C3" variable_1="J_C3"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_III" component_1="dPi_x_dt"/>
      <map_variables variable_2="Pi_x" variable_1="Pi_x"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_III" component_1="Cox_concentration"/>
      <map_variables variable_2="Cox" variable_1="Cox"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_III" component_1="dCred_dt"/>
      <map_variables variable_2="Cred" variable_1="Cred"/>
      <map_variables variable_2="J_C3" variable_1="J_C3"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_III" component_1="Adjustable_parameters"/>
      <map_variables variable_2="k_Pi3" variable_1="k_Pi3"/>
      <map_variables variable_2="k_Pi4" variable_1="k_Pi4"/>
      <map_variables variable_2="x_C3" variable_1="x_C3"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_III" component_1="Fixed_parameters"/>
      <map_variables variable_2="dG_C3o" variable_1="dG_C3o"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="RT" variable_1="RT"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_IV" component_1="Fixed_parameters"/>
      <map_variables variable_2="dG_C4o" variable_1="dG_C4o"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="RT" variable_1="RT"/>
      <map_variables variable_2="Ctot" variable_1="Ctot"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_IV" component_1="Adjustable_parameters"/>
      <map_variables variable_2="k_O2" variable_1="k_O2"/>
      <map_variables variable_2="x_C4" variable_1="x_C4"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_IV" component_1="dH_x_dt"/>
      <map_variables variable_2="H_x" variable_1="H_x"/>
      <map_variables variable_2="J_C4" variable_1="J_C4"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_IV" component_1="dO2_dt"/>
      <map_variables variable_2="O2" variable_1="O2"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_IV" component_1="dCred_dt"/>
      <map_variables variable_2="Cred" variable_1="Cred"/>
      <map_variables variable_2="J_C4" variable_1="J_C4"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_IV" component_1="Cox_concentration"/>
      <map_variables variable_2="Cox" variable_1="Cox"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_IV" component_1="dPsi_dt"/>
      <map_variables variable_2="dPsi" variable_1="dPsi"/>
      <map_variables variable_2="J_C4" variable_1="J_C4"/>
   </connection>
   <connection>
      <map_components component_2="ATP_synthesis_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_F1" variable_1="x_F1"/>
   </connection>
   <connection>
      <map_components component_2="ATP_synthesis_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="K_DD" variable_1="K_DD"/>
      <map_variables variable_2="K_DT" variable_1="K_DT"/>
      <map_variables variable_2="RT" variable_1="RT"/>
      <map_variables variable_2="n_A" variable_1="n_A"/>
      <map_variables variable_2="dG_F1o" variable_1="dG_F1o"/>
   </connection>
   <connection>
      <map_components component_2="ATP_synthesis_flux" component_1="Proton_motive_force"/>
      <map_variables variable_2="dG_H" variable_1="dG_H"/>
   </connection>
   <connection>
      <map_components component_2="ATP_synthesis_flux" component_1="dADP_mx_dt"/>
      <map_variables variable_2="ADP_mx" variable_1="ADP_mx"/>
   </connection>
   <connection>
      <map_components component_2="ATP_synthesis_flux" component_1="dATP_mx_dt"/>
      <map_variables variable_2="ATP_mx" variable_1="ATP_mx"/>
   </connection>
   <connection>
      <map_components component_2="ATP_synthesis_flux" component_1="dPi_x_dt"/>
      <map_variables variable_2="Pi_x" variable_1="Pi_x"/>
      <map_variables variable_2="J_F1" variable_1="J_F1"/>
   </connection>
   <connection>
      <map_components component_2="Dehydrogenase_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_DH" variable_1="x_DH"/>
      <map_variables variable_2="k_Pi1" variable_1="k_Pi1"/>
      <map_variables variable_2="k_Pi2" variable_1="k_Pi2"/>
      <map_variables variable_2="r" variable_1="r"/>
   </connection>
   <connection>
      <map_components component_2="NAD_x_concentration" component_1="Dehydrogenase_flux"/>
      <map_variables variable_2="NAD_x" variable_1="NAD_x"/>
   </connection>
   <connection>
      <map_components component_2="Dehydrogenase_flux" component_1="dNADH_x_dt"/>
      <map_variables variable_2="NADH_x" variable_1="NADH_x"/>
      <map_variables variable_2="J_DH" variable_1="J_DH"/>
   </connection>
   <connection>
      <map_components component_2="Dehydrogenase_flux" component_1="dPi_x_dt"/>
      <map_variables variable_2="Pi_x" variable_1="Pi_x"/>
   </connection>
   <connection>
      <map_components component_2="ANT_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_ANT" variable_1="x_ANT"/>
      <map_variables variable_2="k_mADP" variable_1="k_mADP"/>
   </connection>
   <connection>
      <map_components component_2="ANT_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="RT" variable_1="RT"/>
   </connection>
   <connection>
      <map_components component_2="ANT_flux" component_1="MgATPx_binding_flux"/>
      <map_variables variable_2="ATP_fx" variable_1="ATP_fx"/>
   </connection>
   <connection>
      <map_components component_2="ANT_flux" component_1="MgADPx_binding_flux"/>
      <map_variables variable_2="ADP_fx" variable_1="ADP_fx"/>
   </connection>
   <connection>
      <map_components component_2="ANT_flux" component_1="MgATPi_binding_flux"/>
      <map_variables variable_2="ATP_fi" variable_1="ATP_fi"/>
   </connection>
   <connection>
      <map_components component_2="ANT_flux" component_1="MgADPi_binding_flux"/>
      <map_variables variable_2="ADP_fi" variable_1="ADP_fi"/>
   </connection>
   <connection>
      <map_components component_2="ANT_flux" component_1="dPsi_dt"/>
      <map_variables variable_2="dPsi" variable_1="dPsi"/>
      <map_variables variable_2="J_ANT" variable_1="J_ANT"/>
   </connection>
   <connection>
      <map_components component_2="Phosphate_hydrogen_cotransporter_flux" component_1="dH_x_dt"/>
      <map_variables variable_2="H_x" variable_1="H_x"/>
      <map_variables variable_2="J_Pi1" variable_1="J_Pi1"/>
   </connection>
   <connection>
      <map_components component_2="Phosphate_hydrogen_cotransporter_flux" component_1="dPi_i_dt"/>
      <map_variables variable_2="Pi_i" variable_1="Pi_i"/>
      <map_variables variable_2="J_Pi1" variable_1="J_Pi1"/>
   </connection>
   <connection>
      <map_components component_2="Phosphate_hydrogen_cotransporter_flux" component_1="dPi_x_dt"/>
      <map_variables variable_2="Pi_x" variable_1="Pi_x"/>
      <map_variables variable_2="J_Pi1" variable_1="J_Pi1"/>
   </connection>
   <connection>
      <map_components component_2="Phosphate_hydrogen_cotransporter_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="H_i" variable_1="H_i"/>
      <map_variables variable_2="k_dHPi" variable_1="k_dHPi"/>
   </connection>
   <connection>
      <map_components component_2="Phosphate_hydrogen_cotransporter_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_Pi1" variable_1="x_Pi1"/>
      <map_variables variable_2="k_PiH" variable_1="k_PiH"/>
   </connection>
   <connection>
      <map_components component_2="Hydrogen_leak_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_Hle" variable_1="x_Hle"/>
   </connection>
   <connection>
      <map_components component_2="Hydrogen_leak_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="H_i" variable_1="H_i"/>
      <map_variables variable_2="RT" variable_1="RT"/>
   </connection>
   <connection>
      <map_components component_2="Hydrogen_leak_flux" component_1="dPsi_dt"/>
      <map_variables variable_2="dPsi" variable_1="dPsi"/>
      <map_variables variable_2="J_Hle" variable_1="J_Hle"/>
   </connection>
   <connection>
      <map_components component_2="Hydrogen_leak_flux" component_1="dH_x_dt"/>
      <map_variables variable_2="H_x" variable_1="H_x"/>
      <map_variables variable_2="J_Hle" variable_1="J_Hle"/>
   </connection>
   <connection>
      <map_components component_2="Potassium_hydrogen_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="K_i" variable_1="K_i"/>
      <map_variables variable_2="H_i" variable_1="H_i"/>
   </connection>
   <connection>
      <map_components component_2="Potassium_hydrogen_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_KH" variable_1="x_KH"/>
   </connection>
   <connection>
      <map_components component_2="Potassium_hydrogen_flux" component_1="dK_x_dt"/>
      <map_variables variable_2="K_x" variable_1="K_x"/>
      <map_variables variable_2="J_KH" variable_1="J_KH"/>
   </connection>
   <connection>
      <map_components component_2="Potassium_hydrogen_flux" component_1="dH_x_dt"/>
      <map_variables variable_2="H_x" variable_1="H_x"/>
      <map_variables variable_2="J_KH" variable_1="J_KH"/>
   </connection>
   <connection>
      <map_components component_2="Passive_potassium_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="K_i" variable_1="K_i"/>
      <map_variables variable_2="F" variable_1="F"/>
      <map_variables variable_2="RT" variable_1="RT"/>
   </connection>
   <connection>
      <map_components component_2="Passive_potassium_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_K" variable_1="x_K"/>
   </connection>
   <connection>
      <map_components component_2="Passive_potassium_flux" component_1="dPsi_dt"/>
      <map_variables variable_2="dPsi" variable_1="dPsi"/>
      <map_variables variable_2="J_K" variable_1="J_K"/>
   </connection>
   <connection>
      <map_components component_2="Passive_potassium_flux" component_1="dK_x_dt"/>
      <map_variables variable_2="K_x" variable_1="K_x"/>
      <map_variables variable_2="J_K" variable_1="J_K"/>
   </connection>
   <connection>
      <map_components component_2="Adenylate_kinase_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="K_AK" variable_1="K_AK"/>
   </connection>
   <connection>
      <map_components component_2="Adjustable_parameters" component_1="Adenylate_kinase_flux"/>
      <map_variables variable_2="x_AK" variable_1="x_AK"/>
   </connection>
   <connection>
      <map_components component_2="Adenylate_kinase_flux" component_1="dADP_i_dt"/>
      <map_variables variable_2="ADP_i" variable_1="ADP_i"/>
      <map_variables variable_2="J_AKi" variable_1="J_AKi"/>
   </connection>
   <connection>
      <map_components component_2="Adenylate_kinase_flux" component_1="dAMP_i_dt"/>
      <map_variables variable_2="AMP_i" variable_1="AMP_i"/>
      <map_variables variable_2="J_AKi" variable_1="J_AKi"/>
   </connection>
   <connection>
      <map_components component_2="Adenylate_kinase_flux" component_1="dATP_i_dt"/>
      <map_variables variable_2="ATP_i" variable_1="ATP_i"/>
      <map_variables variable_2="J_AKi" variable_1="J_AKi"/>
   </connection>
   <connection>
      <map_components component_2="AMP_substrate_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="AMP_e" variable_1="AMP_e"/>
      <map_variables variable_2="gamma" variable_1="gamma"/>
   </connection>
   <connection>
      <map_components component_2="AMP_substrate_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="p_A" variable_1="p_A"/>
   </connection>
   <connection>
      <map_components component_2="AMP_substrate_flux" component_1="dAMP_i_dt"/>
      <map_variables variable_2="AMP_i" variable_1="AMP_i"/>
      <map_variables variable_2="J_AMP" variable_1="J_AMP"/>
   </connection>
   <connection>
      <map_components component_2="ADP_substrate_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="gamma" variable_1="gamma"/>
   </connection>
   <connection>
      <map_components component_2="ADP_substrate_flux" component_1="Environment"/>
      <map_variables variable_2="ADP_e" variable_1="ADP_e"/>
   </connection>
   <connection>
      <map_components component_2="ADP_substrate_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="p_A" variable_1="p_A"/>
   </connection>
   <connection>
      <map_components component_2="ADP_substrate_flux" component_1="dADP_i_dt"/>
      <map_variables variable_2="ADP_i" variable_1="ADP_i"/>
      <map_variables variable_2="J_ADP" variable_1="J_ADP"/>
   </connection>
   <connection>
      <map_components component_2="ATP_substrate_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="ATP_e" variable_1="ATP_e"/>
      <map_variables variable_2="gamma" variable_1="gamma"/>
   </connection>
   <connection>
      <map_components component_2="ATP_substrate_flux" component_1="dATP_i_dt"/>
      <map_variables variable_2="ATP_i" variable_1="ATP_i"/>
      <map_variables variable_2="J_ATP" variable_1="J_ATP"/>
   </connection>
   <connection>
      <map_components component_2="Pi_substrate_flux" component_1="Environment"/>
      <map_variables variable_2="Pi_e" variable_1="Pi_e"/>
   </connection>
   <connection>
      <map_components component_2="Pi_substrate_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="gamma" variable_1="gamma"/>
   </connection>
   <connection>
      <map_components component_2="Pi_substrate_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_Pi2" variable_1="x_Pi2"/>
   </connection>
   <connection>
      <map_components component_2="Pi_substrate_flux" component_1="dPi_i_dt"/>
      <map_variables variable_2="Pi_i" variable_1="Pi_i"/>
      <map_variables variable_2="J_Pi2" variable_1="J_Pi2"/>
   </connection>
   <connection>
      <map_components component_2="MgATPx_binding_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="K_DT" variable_1="K_DT"/>
   </connection>
   <connection>
      <map_components component_2="MgATPx_binding_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_MgA" variable_1="x_MgA"/>
   </connection>
   <connection>
      <map_components component_2="MgATPx_binding_flux" component_1="dMg_x_dt"/>
      <map_variables variable_2="Mg_x" variable_1="Mg_x"/>
      <map_variables variable_2="J_MgATPx" variable_1="J_MgATPx"/>
   </connection>
   <connection>
      <map_components component_2="MgADPx_binding_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="K_DD" variable_1="K_DD"/>
   </connection>
   <connection>
      <map_components component_2="MgADPx_binding_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_MgA" variable_1="x_MgA"/>
   </connection>
   <connection>
      <map_components component_2="MgADPx_binding_flux" component_1="dMg_x_dt"/>
      <map_variables variable_2="Mg_x" variable_1="Mg_x"/>
      <map_variables variable_2="J_MgADPx" variable_1="J_MgADPx"/>
   </connection>
   <connection>
      <map_components component_2="MgATPi_binding_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="K_DT" variable_1="K_DT"/>
   </connection>
   <connection>
      <map_components component_2="MgATPi_binding_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_MgA" variable_1="x_MgA"/>
   </connection>
   <connection>
      <map_components component_2="MgATPi_binding_flux" component_1="Mg_binding"/>
      <map_variables variable_2="Mg_i" variable_1="Mg_i"/>
   </connection>
   <connection>
      <map_components component_2="MgADPi_binding_flux" component_1="Fixed_parameters"/>
      <map_variables variable_2="K_DD" variable_1="K_DD"/>
   </connection>
   <connection>
      <map_components component_2="MgADPi_binding_flux" component_1="Adjustable_parameters"/>
      <map_variables variable_2="x_MgA" variable_1="x_MgA"/>
   </connection>
   <connection>
      <map_components component_2="MgADPi_binding_flux" component_1="Mg_binding"/>
      <map_variables variable_2="Mg_i" variable_1="Mg_i"/>
   </connection>
   <connection>
      <map_components component_2="dH_x_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dH_x_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="n_A" variable_1="n_A"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="dH_x_dt" component_1="Adjustable_parameters"/>
      <map_variables variable_2="r_buff" variable_1="r_buff"/>
   </connection>
   <connection>
      <map_components component_2="dH_x_dt" component_1="Dehydrogenase_flux"/>
      <map_variables variable_2="J_DH" variable_1="J_DH"/>
   </connection>
   <connection>
      <map_components component_2="dH_x_dt" component_1="ATP_synthesis_flux"/>
      <map_variables variable_2="J_F1" variable_1="J_F1"/>
   </connection>
   <connection>
      <map_components component_2="Environment" component_1="dK_x_dt"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dK_x_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="dMg_x_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dMg_x_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="dNADH_x_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dNADH_x_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="dQH2_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dQH2_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="dCred_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dCred_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_i" variable_1="W_i"/>
   </connection>
   <connection>
      <map_components component_2="dO2_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dATP_x_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dATP_x_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="dATP_x_dt" component_1="ATP_synthesis_flux"/>
      <map_variables variable_2="J_F1" variable_1="J_F1"/>
   </connection>
   <connection>
      <map_components component_2="ANT_flux" component_1="dATP_x_dt"/>
      <map_variables variable_2="J_ANT" variable_1="J_ANT"/>
   </connection>
   <connection>
      <map_components component_2="dADP_x_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dADP_x_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="dADP_x_dt" component_1="ATP_synthesis_flux"/>
      <map_variables variable_2="J_F1" variable_1="J_F1"/>
   </connection>
   <connection>
      <map_components component_2="ANT_flux" component_1="dADP_x_dt"/>
      <map_variables variable_2="J_ANT" variable_1="J_ANT"/>
   </connection>
   <connection>
      <map_components component_2="Environment" component_1="dATP_mx_dt"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dATP_mx_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="Environment" component_1="dADP_mx_dt"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dADP_mx_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="Environment" component_1="dPi_x_dt"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dPi_x_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_x" variable_1="W_x"/>
   </connection>
   <connection>
      <map_components component_2="Environment" component_1="dATP_i_dt"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dATP_i_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_i" variable_1="W_i"/>
   </connection>
   <connection>
      <map_components component_2="dATP_i_dt" component_1="ANT_flux"/>
      <map_variables variable_2="J_ANT" variable_1="J_ANT"/>
   </connection>
   <connection>
      <map_components component_2="dADP_i_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dADP_i_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_i" variable_1="W_i"/>
   </connection>
   <connection>
      <map_components component_2="dADP_i_dt" component_1="ANT_flux"/>
      <map_variables variable_2="J_ANT" variable_1="J_ANT"/>
   </connection>
   <connection>
      <map_components component_2="dAMP_i_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dAMP_i_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_i" variable_1="W_i"/>
   </connection>
   <connection>
      <map_components component_2="dATP_mi_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dATP_mi_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_i" variable_1="W_i"/>
   </connection>
   <connection>
      <map_components component_2="dADP_mi_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dADP_mi_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_i" variable_1="W_i"/>
   </connection>
   <connection>
      <map_components component_2="dPi_i_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dPi_i_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="W_i" variable_1="W_i"/>
   </connection>
   <connection>
      <map_components component_2="dPsi_dt" component_1="Environment"/>
      <map_variables variable_2="t" variable_1="t"/>
   </connection>
   <connection>
      <map_components component_2="dPsi_dt" component_1="Fixed_parameters"/>
      <map_variables variable_2="n_A" variable_1="n_A"/>
   </connection>
   <connection>
      <map_components component_2="dPsi_dt" component_1="Electron_flux_complex_I"/>
      <map_variables variable_2="J_C1" variable_1="J_C1"/>
   </connection>
   <connection>
      <map_components component_2="dPsi_dt" component_1="Electron_flux_complex_III"/>
      <map_variables variable_2="J_C3" variable_1="J_C3"/>
      <map_variables variable_2="dPsi" variable_1="dPsi"/>
   </connection>
   <connection>
      <map_components component_2="ATP_synthesis_flux" component_1="dPsi_dt"/>
      <map_variables variable_2="J_F1" variable_1="J_F1"/>
   </connection>
   <connection>
      <map_components component_2="NAD_x_concentration" component_1="Fixed_parameters"/>
      <map_variables variable_2="NADtot" variable_1="NADtot"/>
   </connection>
   <connection>
      <map_components component_2="NAD_x_concentration" component_1="dNADH_x_dt"/>
      <map_variables variable_2="NADH_x" variable_1="NADH_x"/>
   </connection>
   <connection>
      <map_components component_2="Q_concentration" component_1="Fixed_parameters"/>
      <map_variables variable_2="Qtot" variable_1="Qtot"/>
   </connection>
   <connection>
      <map_components component_2="Q_concentration" component_1="dQH2_dt"/>
      <map_variables variable_2="QH2" variable_1="QH2"/>
   </connection>
   <connection>
      <map_components component_2="Cox_concentration" component_1="Fixed_parameters"/>
      <map_variables variable_2="Ctot" variable_1="Ctot"/>
   </connection>
   <connection>
      <map_components component_2="Cox_concentration" component_1="dCred_dt"/>
      <map_variables variable_2="Cred" variable_1="Cred"/>
   </connection>
   <connection>
      <map_components component_2="Proton_motive_force" component_1="Electron_flux_complex_III"/>
      <map_variables variable_2="dG_H" variable_1="dG_H"/>
   </connection>
   <connection>
      <map_components component_2="Electron_flux_complex_IV" component_1="Proton_motive_force"/>
      <map_variables variable_2="dG_H" variable_1="dG_H"/>
   </connection>
   <connection>
      <map_components component_2="Adjustable_parameters" component_1="ATP_substrate_flux"/>
      <map_variables variable_2="p_A" variable_1="p_A"/>
   </connection>
<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#a74e020f-db63-46fb-a32d-17aae23fd10a">
    <rdf:li rdf:resource="rdf:#56dbf4b0-d78c-4cbb-93e4-7603c65b3a06"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#15829d8d-0717-4c07-a249-a20f31b92658">
    <dc:creator rdf:resource="rdf:#4ad2dfe9-6eab-4941-bdc8-83d5c4d60a90"/>
    <rdf:value>This model has been edited to correct dimensional errors and checked for other problems. It runs in COR and PCENV but whether it replicates experimental results is yet to be determined.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4ad2dfe9-6eab-4941-bdc8-83d5c4d60a90">
    <vCard:FN>Geoff Nunns</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Research Group</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#15829d8d-0717-4c07-a249-a20f31b92658"/>
    <dcterms:created rdf:resource="rdf:#3116256b-abf8-4da2-952c-86e693c76550"/>
    <dc:creator rdf:resource="rdf:#9998c1fe-b38a-4a3a-9d54-91a05bc4400c"/>
    <cmeta:modification rdf:resource="rdf:#d57e9f1d-a302-4061-8456-048a239e4870"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aac21dd2-1f68-445a-9fea-c0d4c2c3fc5c">
    <vCard:N rdf:resource="rdf:#0f617f85-5f1e-4365-af17-67e01aed621d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#56dbf4b0-d78c-4cbb-93e4-7603c65b3a06">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#37d25a29-0ad5-479d-85a7-5730fb29b55c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9998c1fe-b38a-4a3a-9d54-91a05bc4400c">
    <vCard:ORG rdf:resource="rdf:#2cf3c72c-54b5-4d94-b0d8-29651e70da7d"/>
    <vCard:EMAIL rdf:resource="rdf:#46a48b1a-83ad-43fb-958e-5bd80c5a409e"/>
    <vCard:N rdf:resource="rdf:#f3f10030-cdb2-446f-b95d-ad17b9c7576c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0f617f85-5f1e-4365-af17-67e01aed621d">
    <vCard:Given>Geoffrey</vCard:Given>
    <vCard:Family>Nunns</vCard:Family>
    <vCard:Other>Rogan</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7634e3ea-7ab7-41e1-8426-7597614c9eeb">
    <dc:creator rdf:resource="rdf:#ced5e5e5-2a4c-4fae-b6e3-f905b7f7bb0c"/>
    <rdf:value>This model has been edited to correct dimensional errors and checked for other problems. It runs in COR and PCENV but whether it replicates experimental results is yet to be determined.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ced5e5e5-2a4c-4fae-b6e3-f905b7f7bb0c">
    <vCard:FN>Geoff Nunns</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c2e5dfc2-27b8-4572-80f7-a32acded401d">
    <bqs:Pubmed_id>16163394</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#1a689e5f-0fe8-497a-aca6-2b0922850804"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d57e9f1d-a302-4061-8456-048a239e4870">
    <dcterms:modified rdf:resource="rdf:#a5dcd4e6-33a3-45b4-b862-943e3aec224e"/>
    <rdf:value>Fixed &quot;not exactly equivalent but dimensionally equivalent&quot; errors.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#aac21dd2-1f68-445a-9fea-c0d4c2c3fc5c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a5dcd4e6-33a3-45b4-b862-943e3aec224e">
    <dcterms:W3CDTF>2008-07-17T14:26:02+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4c7055ad-2b89-4233-9d6b-ce8e25323f53">
    <dc:title>PLoS Compututational Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="#beard_model_2005">
    <bqs:reference rdf:resource="rdf:#c2e5dfc2-27b8-4572-80f7-a32acded401d"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>mitochondria</rdf:li>
        <rdf:li>respiration</rdf:li>
        <rdf:li>metabolism</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <cmeta:comment rdf:resource="rdf:#7634e3ea-7ab7-41e1-8426-7597614c9eeb"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1a689e5f-0fe8-497a-aca6-2b0922850804">
    <dc:creator rdf:resource="rdf:#a74e020f-db63-46fb-a32d-17aae23fd10a"/>
    <dc:title>A Biophysical Model of the Mitochondrial Respiratory System and Oxidative Phosphorylation</dc:title>
    <bqs:volume/>
    <bqs:first_page/>
    <bqs:Journal rdf:resource="rdf:#4c7055ad-2b89-4233-9d6b-ce8e25323f53"/>
    <dcterms:issued rdf:resource="rdf:#3085d114-3642-400c-b8b4-d3055a056a6f"/>
    <bqs:last_page/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#37d25a29-0ad5-479d-85a7-5730fb29b55c">
    <vCard:Given>Daniel</vCard:Given>
    <vCard:Family>Beard</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3085d114-3642-400c-b8b4-d3055a056a6f">
    <dcterms:W3CDTF>2005-09-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3116256b-abf8-4da2-952c-86e693c76550">
    <dcterms:W3CDTF>2005-09-09T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f3f10030-cdb2-446f-b95d-ad17b9c7576c">
    <vCard:Given>Peter</vCard:Given>
    <vCard:Family>Villager</vCard:Family>
    <vCard:Other>John</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2cf3c72c-54b5-4d94-b0d8-29651e70da7d">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Bioengineering Research Group</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#46a48b1a-83ad-43fb-958e-5bd80c5a409e">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value/>
  </rdf:Description>
</rdf:RDF>
</model>