Location: Beard, 2005 @ 8528ecaa1770 / beard_2005.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2008-06-12 03:09:32+12:00
Desc:
committing version02 of beard_2005
Permanent Source URI:
http://models.cellml.org/workspace/beard_2005/rawfile/8528ecaa17708dc9078eb4e038fd631a461633cf/beard_2005.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 12/06/2008 at 2:51:59 at p.m. using:

COR (0.9.31.955)
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 PCENV but whether it replicates experimental results is yet to be determined. 
          </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>
					<ulink url="http://compbiol.plosjournals.org/perlserv/?request=get-document&amp;doi=10.1371/journal.pcbi.0010036">A Biophysical Model of the Mitochondrial Respiratory System and Oxidative Phosphorylation</ulink>, Daniel A. Beard, 2005, <ulink url="http://compbiol.plosjournals.org/">
						<emphasis>PLoS Computational Biology</emphasis>
					</ulink>, 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>
			</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="millisecond">
      <unit units="second" prefix="milli"/>
   </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="millisecond" 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">1e-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">1e-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">1e-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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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="millisecond" 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:#a6ce8bbc-5d87-471d-b4c7-d71bcae7f5df">
    <rdf:li rdf:resource="rdf:#6f31e477-db53-4ca1-9afa-665a219b9652"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#d8ba8df5-3533-47e4-af37-04ae3280c428">
    <vCard:Given>Geoffrey</vCard:Given>
    <vCard:Family>Nunns</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#961280a6-ded8-466d-b636-4193f50a32bc">
    <vCard:Given>Peter</vCard:Given>
    <vCard:Family>Villager</vCard:Family>
    <vCard:Other>John</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#80f81b84-6b74-49df-9e85-dceb34fce228">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Bioengineering Research Group</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#70a020dd-3a56-4ff5-972d-bf465e9e10c1">
    <bqs:Pubmed_id>16163394</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#48021d00-dfb5-4d42-8633-66f228ae6e97"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Research Group</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#dbad833e-5dc0-4cb6-97c8-f2bae1ad1a84"/>
    <dcterms:created rdf:resource="rdf:#d7c740a0-a3c0-4e33-ac66-859d288f9124"/>
    <dc:creator rdf:resource="rdf:#aae0bded-f661-44c9-849a-ae600c8ca640"/>
    <cmeta:modification rdf:resource="rdf:#7e3b3219-d28c-45be-9505-6eec1d6b93b2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2cfcc773-6ee1-4a7f-8812-8a0013d5d25d">
    <dc:creator rdf:resource="rdf:#1af8e08b-332e-45c9-80f2-3efc4f43f5b0"/>
    <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:#6f31e477-db53-4ca1-9afa-665a219b9652">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#61126b0d-1928-486c-8a68-5c18ddc69a1d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#949ff7f4-6a1e-47cb-8f08-6ed94411614b">
    <vCard:N rdf:resource="rdf:#d8ba8df5-3533-47e4-af37-04ae3280c428"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bc406da4-c974-4de7-86ce-6cd2167c2112">
    <dcterms:W3CDTF>2008-06-12T15:06:58+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1290811b-c879-4ead-a04e-954f71f5c877">
    <vCard:FN>Geoff Nunns</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2ee38178-e31e-4f30-96fd-1fa77485e20a">
    <dc:title>PLoS Compututational Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dbad833e-5dc0-4cb6-97c8-f2bae1ad1a84">
    <dc:creator rdf:resource="rdf:#1290811b-c879-4ead-a04e-954f71f5c877"/>
    <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:#7e3b3219-d28c-45be-9505-6eec1d6b93b2">
    <dcterms:modified rdf:resource="rdf:#bc406da4-c974-4de7-86ce-6cd2167c2112"/>
    <rdf:value>This version has updated units and will now run in COR and PCENV.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#949ff7f4-6a1e-47cb-8f08-6ed94411614b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#61126b0d-1928-486c-8a68-5c18ddc69a1d">
    <vCard:Given>Daniel</vCard:Given>
    <vCard:Family>Beard</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d7c740a0-a3c0-4e33-ac66-859d288f9124">
    <dcterms:W3CDTF>2005-09-09T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#28be1893-04a5-457e-abfb-e8e87535506b">
    <dcterms:W3CDTF>2005-09-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#48021d00-dfb5-4d42-8633-66f228ae6e97">
    <dc:creator rdf:resource="rdf:#a6ce8bbc-5d87-471d-b4c7-d71bcae7f5df"/>
    <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:#2ee38178-e31e-4f30-96fd-1fa77485e20a"/>
    <dcterms:issued rdf:resource="rdf:#28be1893-04a5-457e-abfb-e8e87535506b"/>
    <bqs:last_page/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aae0bded-f661-44c9-849a-ae600c8ca640">
    <vCard:ORG rdf:resource="rdf:#80f81b84-6b74-49df-9e85-dceb34fce228"/>
    <vCard:EMAIL rdf:resource="rdf:#8f818abe-355f-4ed1-b5ea-7e9d5dac8237"/>
    <vCard:N rdf:resource="rdf:#961280a6-ded8-466d-b636-4193f50a32bc"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8f818abe-355f-4ed1-b5ea-7e9d5dac8237">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1af8e08b-332e-45c9-80f2-3efc4f43f5b0">
    <vCard:FN>Geoff Nunns</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="#beard_model_2005">
    <bqs:reference rdf:resource="rdf:#70a020dd-3a56-4ff5-972d-bf465e9e10c1"/>
    <cmeta:comment rdf:resource="rdf:#2cfcc773-6ee1-4a7f-8812-8a0013d5d25d"/>
  </rdf:Description>
</rdf:RDF>
</model>