Rendering of the source text

<?xml version="1.0"?>
<!--
This CellML file was generated on 6/01/2010 at 10:28:40 at a.m. using:

COR (0.9.31.1333)
Copyright 2002-2010 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#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" name="bertram_2006" cmeta:id="bertram_2006">
	
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>A simplified model for mitochondrial ATP production</title>
				<author>
					<firstname>Tessa</firstname>
					<surname>Paris</surname>
					<affiliation>
						<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para> This CellML model is based on the author's original code (published on Richard Bertram's website). The CellML model runs in both OpenCell and COR and can reproduce some of the results from some parts of figures 3 to 9.</para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>ABSTRACT: Most of the adenosine triphosphate (ATP) synthesized during glucose metabolism is produced in the mitochondria through oxidativephosphorylation. This is a complex reaction powered by the proton gradient across the mitochondrial inner membrane, which isgenerated by mitochondrial respiration. A detailed model of this reaction, which includes dynamic equations for the key mitochondrialvariables, was developed earlier by Magnus and Keizer. However, this model is extraordinarily complicated. We develop a simpler modelthat captures the behavior of the original model but is easier to use and to understand. We then use it to investigate the mitochondrialresponses to glycolytic and calcium input. WWe use the model to explain experimental observations of the opposite effects of raising cytosolic Ca(2+)in low and high glucose, and to predict the effects of a mutation in the mitochondrial enzyme nicotinamide nucleotide transhydrogenase (Nnt) in pancreatic beta-cells.</para>
				
				<para>The original paper reference is cited below:</para>
				
				<para>A simplified model for mitochondrial ATP production, Richard Bertram, Morten Gram Pedersenb, Dan S. Lucianic and Arthur Shermand, 2006 <emphasis>Journal of Theoretical Biology</emphasis>
                              , 243 (2006), 575<ulink url="http://www.ncbi.nlm.nih.gov/pubmed/16945388?">PubMed ID: 16945388</ulink>
				</para>
				
				<informalfigure float="0" id="fig_cell_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>cell schematic for the model</title>
							</objectinfo>
							<imagedata fileref="bertram_2006.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the fluxes and reactions used in the model.</caption>
				</informalfigure>
			</sect1>
		</article>
	</documentation>   
	
	<units name="mM">
      <unit prefix="milli" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="uM">
      <unit prefix="micro" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="mV">
      <unit prefix="milli" units="volt"/>
   </units>
   <units name="ms">
      <unit prefix="milli" units="second"/>
   </units>
   <units name="per_ms">
      <unit exponent="-1" units="ms"/>
   </units>
   <units name="uM_per_ms">
      <unit units="uM"/>
      <unit exponent="-1" units="ms"/>
   </units>
   <units name="per_uM_per_ms_per_mV">
      <unit exponent="-1" units="uM"/>
      <unit exponent="-1" units="ms"/>
      <unit exponent="-1" units="mV"/>
   </units>
   <units name="uM_per_mV">
      <unit units="uM"/>
      <unit exponent="-1" units="mV"/>
   </units>
   <units name="per_uM_per_ms">
      <unit exponent="-1" units="uM"/>
      <unit exponent="-1" units="ms"/>
   </units>
   <units name="per_mV">
      <unit exponent="-1" units="mV"/>
   </units>
   <units name="uM_per_ms_per_mV">
      <unit units="uM"/>
      <unit exponent="-1" units="ms"/>
      <unit exponent="-1" units="mV"/>
   </units>
   <units name="C_per_mol">
      <unit units="coulomb"/>
      <unit exponent="-1" units="mole"/>
   </units>
   <units name="mJ_per_K_per_mol">
      <unit prefix="milli" units="joule"/>
      <unit exponent="-1" units="kelvin"/>
      <unit exponent="-1" units="mole"/>
   </units>
   <component name="Environment">
      <variable cmeta:id="Environment_time" name="time" public_interface="out" units="ms"/>
   </component>
   <component name="Mitochondrial_variables">
      <variable initial_value="0.6" name="NADHm" public_interface="out" units="mM"/>
      <variable name="Jpdh" public_interface="in" units="uM_per_ms"/>
      <variable name="Jo" public_interface="in" units="uM_per_ms"/>
      <variable initial_value="0.001" name="gamma" units="dimensionless"/>
      <variable name="NADm" public_interface="out" units="mM"/>
      <variable initial_value="7.4" name="ADPm" public_interface="out" units="mM"/>
      <variable initial_value="93" name="PSIm" public_interface="out" units="mV"/>
      <variable initial_value="0.1" name="Cam" public_interface="out" units="uM"/>
      <variable initial_value="10" name="NADtot" units="mM"/>
      <variable initial_value="0.01" name="fm" units="dimensionless"/>
      <variable initial_value="1.8" name="Cmito" units="uM_per_mV"/>
      <variable initial_value="15" name="Amtot" units="mM"/>
      <variable name="ATPm" public_interface="out" units="mM"/>
      <variable name="RATm" public_interface="out" units="dimensionless"/>
      <variable name="Jant" public_interface="in" units="uM_per_ms"/>
      <variable name="Jf1f0" public_interface="in" units="uM_per_ms"/>
      <variable name="Jh_res" public_interface="in" units="uM_per_ms"/>
      <variable name="Jh_atp" public_interface="in" units="uM_per_ms"/>
      <variable name="Jh_leak" public_interface="in" units="uM_per_ms"/>
      <variable name="Jnaca" public_interface="in" units="uM_per_ms"/>
      <variable name="Juni" public_interface="in" units="uM_per_ms"/>
      <variable name="time" public_interface="in" units="ms"/>
      <variable name="Jmito" public_interface="in" units="uM_per_ms"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>NADHm</ci>
            </apply>
            <apply>
               <times/>
               <ci>gamma</ci>
               <apply>
                  <minus/>
                  <ci>Jpdh</ci>
                  <ci>Jo</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>NADm</ci>
            <apply>
               <minus/>
               <ci>NADtot</ci>
               <ci>NADHm</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ATPm</ci>
            <apply>
               <minus/>
               <ci>Amtot</ci>
               <ci>ADPm</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ADPm</ci>
            </apply>
            <apply>
               <times/>
               <ci>gamma</ci>
               <apply>
                  <minus/>
                  <ci>Jant</ci>
                  <ci>Jf1f0</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>RATm</ci>
            <apply>
               <divide/>
               <ci>ATPm</ci>
               <ci>ADPm</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>PSIm</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>Jh_res</ci>
                              <ci>Jh_atp</ci>
                           </apply>
                           <ci>Jant</ci>
                        </apply>
                        <ci>Jh_leak</ci>
                     </apply>
                     <ci>Jnaca</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>Juni</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <ci>Cmito</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Cam</ci>
            </apply>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>fm</ci>
               </apply>
               <ci>Jmito</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Cytosol">
      <variable name="time" public_interface="in" units="ms"/>
      <variable name="c" public_interface="in" units="uM"/>
      <variable name="ATPc" public_interface="out" units="mM"/>
      <variable initial_value="1850" name="ADPc" units="uM"/>
      <variable initial_value="2500" name="Ac_tot" units="uM"/>
      <variable name="Jant" public_interface="in" units="uM_per_ms"/>
      <variable initial_value="0.00005" name="khyd" units="dimensionless"/>
      <variable initial_value="0.00005" name="Jhydbas" units="mM"/>
      <variable name="Jh_yd" units="mM"/>
      <variable name="delta" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>delta</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3.9</cn>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">53.2</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>ATPc</ci>
            <apply>
               <minus/>
               <ci>Ac_tot</ci>
               <ci>ADPc</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jh_yd</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>khyd</ci>
                        <ci>c</ci>
                     </apply>
                     <ci>Jhydbas</ci>
                  </apply>
                  <ci>ATPc</ci>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="uM">1</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>ADPc</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>delta</ci>
                  </apply>
                  <ci>Jant</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>Jh_yd</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_ms">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Clamp_protocol">
      <variable name="time" private_interface="out" public_interface="in" units="ms"/>
      <variable initial_value="0.1" name="c" public_interface="out" units="uM"/>
      <variable name="Fproto" private_interface="in" units="uM"/>
      <variable name="Cproto" private_interface="in" units="uM"/>
      <variable initial_value="0.5" name="FBP" public_interface="out" units="uM"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>FBP</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Fproto</ci>
                  <ci>FBP</ci>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0.1</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>c</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Cproto</ci>
                  <ci>c</ci>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0.1</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Clamp_parameters">
      <variable name="time" public_interface="in" units="ms"/>
      <variable initial_value="1" name="Fhold" units="uM"/>
      <variable initial_value="5" name="Ftest" units="uM"/>
      <variable initial_value="90000" name="Fton" units="ms"/>
      <variable initial_value="330000" name="Ftoff" units="ms"/>
      <variable name="Fproto" public_interface="out" units="uM"/>
      <variable name="Cproto" public_interface="out" units="uM"/>
      <variable initial_value="0.1" name="Chold" units="uM"/>
      <variable initial_value="0.1" name="Ctest" units="uM"/>
      <variable initial_value="120000" name="Cton1" units="ms"/>
      <variable initial_value="180000" name="Cton2" units="ms"/>
      <variable initial_value="240000" name="Cton3" units="ms"/>
      <variable initial_value="150000" name="Ctoff1" units="ms"/>
      <variable initial_value="210000" name="Ctoff2" units="ms"/>
      <variable initial_value="270000" name="Ctoff3" units="ms"/>
      <variable name="pulse1" units="uM"/>
      <variable name="pulse2" units="uM"/>
      <variable name="pulse3" units="uM"/>
      <variable name="heav_on" units="dimensionless"/>
      <variable name="heav_off" units="dimensionless"/>
      <variable name="heav_Cton1" units="dimensionless"/>
      <variable name="heav_Cton2" units="dimensionless"/>
      <variable name="heav_Cton3" units="dimensionless"/>
      <variable name="heav_Ctoff1" units="dimensionless"/>
      <variable name="heav_Ctoff2" units="dimensionless"/>
      <variable name="heav_Ctoff3" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>pulse1</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>Ctest</ci>
                  <ci>Chold</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>heav_Cton1</ci>
                  <ci>heav_Ctoff1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>pulse2</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>Ctest</ci>
                  <ci>Chold</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>heav_Cton2</ci>
                  <ci>heav_Ctoff2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>pulse3</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>Ctest</ci>
                  <ci>Chold</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>heav_Cton3</ci>
                  <ci>heav_Ctoff3</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Cproto</ci>
            <apply>
               <plus/>
               <ci>Chold</ci>
               <ci>pulse1</ci>
               <ci>pulse2</ci>
               <ci>pulse3</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Fproto</ci>
            <apply>
               <plus/>
               <ci>Fhold</ci>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>Ftest</ci>
                     <ci>Fhold</ci>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>heav_on</ci>
                     <ci>heav_off</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>heav_on</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Fton</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_off</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Ftoff</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Cton1</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Cton1</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Cton2</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Cton2</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Cton3</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Cton3</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Ctoff1</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Ctoff1</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Ctoff2</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Ctoff2</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Ctoff3</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Ctoff3</ci>
                     </apply>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="J_variables">
      <variable cmeta:id="J_variables_Juni" name="Juni" public_interface="out" units="uM_per_ms"/>
      <variable cmeta:id="J_variables_Jpdh" name="Jpdh" public_interface="out" units="uM_per_ms"/>
      <variable cmeta:id="J_variables_Jnaca" name="Jnaca" public_interface="out" units="uM_per_ms"/>
      <variable name="Jmito" public_interface="out" units="uM_per_ms"/>
      <variable cmeta:id="J_variables_Jh_leak" name="Jh_leak" public_interface="out" units="uM_per_ms"/>
      <variable cmeta:id="J_variables_Jh_atp" name="Jh_atp" public_interface="out" units="uM_per_ms"/>
      <variable cmeta:id="J_variables_Jo" name="Jo" public_interface="out" units="uM_per_ms"/>
      <variable cmeta:id="J_variables_Jh_res" name="Jh_res" public_interface="out" units="uM_per_ms"/>
      <variable cmeta:id="J_variables_Jf1f0" name="Jf1f0" public_interface="out" units="uM_per_ms"/>
      <variable cmeta:id="J_variables_Jant" name="Jant" public_interface="out" units="uM_per_ms"/>
      <variable cmeta:id="J_variables_Jgpdh" name="Jgpdh" units="uM_per_ms"/>
      <variable name="c" public_interface="in" units="uM"/>
      <variable name="PSIm" public_interface="in" units="mV"/>
      <variable name="Cam" public_interface="in" units="uM"/>
      <variable name="NADHm" public_interface="in" units="mM"/>
      <variable name="p1" private_interface="in" units="dimensionless"/>
      <variable name="p2" private_interface="in" units="dimensionless"/>
      <variable name="p3" private_interface="in" units="uM"/>
      <variable name="p21" private_interface="in" units="per_uM_per_ms_per_mV"/>
      <variable name="p22" private_interface="in" units="per_uM_per_ms"/>
      <variable name="p23" private_interface="in" units="uM_per_ms"/>
      <variable name="p24" private_interface="in" units="per_mV"/>
      <variable name="p4" private_interface="in" units="uM_per_ms"/>
      <variable name="p5" private_interface="in" units="mM"/>
      <variable name="p6" private_interface="in" units="mV"/>
      <variable name="p7" private_interface="in" units="mV"/>
      <variable name="p8" private_interface="in" units="uM_per_ms"/>
      <variable name="p9" private_interface="in" units="mM"/>
      <variable name="p10" private_interface="in" units="mV"/>
      <variable name="p11" private_interface="in" units="mV"/>
      <variable name="p12" private_interface="in" units="uM_per_ms"/>
      <variable name="p13" private_interface="in" units="mM"/>
      <variable name="p14" private_interface="in" units="mV"/>
      <variable name="p15" private_interface="in" units="mV"/>
      <variable name="p16" private_interface="in" units="uM_per_ms"/>
      <variable name="p17" private_interface="in" units="uM_per_ms_per_mV"/>
      <variable name="p18" private_interface="in" units="uM_per_ms"/>
      <variable name="p19" private_interface="in" units="uM_per_ms"/>
      <variable name="p20" private_interface="in" units="dimensionless"/>
      <variable name="MM1" units="uM_per_ms"/>
      <variable name="MM2" units="uM_per_ms"/>
      <variable name="b13" units="uM_per_ms"/>
      <variable name="b2" units="uM_per_ms"/>
      <variable name="NADm" public_interface="in" units="mM"/>
      <variable name="ATPm" public_interface="in" units="mM"/>
      <variable name="ATPc" public_interface="in" units="mM"/>
      <variable name="RATm" public_interface="in" units="dimensionless"/>
      <variable initial_value="0.037" name="FRT" units="per_mV"/>
      <variable initial_value="0.0005" name="kgpdh" units="uM_per_ms"/>
      <variable name="FBP" public_interface="in" units="uM"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Juni</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>p21</ci>
                     <ci>PSIm</ci>
                  </apply>
                  <ci>p22</ci>
               </apply>
               <apply>
                  <power/>
                  <ci>c</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jpdh</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <ci>p1</ci>
                  <apply>
                     <plus/>
                     <ci>p2</ci>
                     <apply>
                        <divide/>
                        <ci>NADHm</ci>
                        <ci>NADm</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <ci>Cam</ci>
                  <apply>
                     <plus/>
                     <ci>p3</ci>
                     <ci>Cam</ci>
                  </apply>
               </apply>
               <ci>Jgpdh</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jnaca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>p23</ci>
                     <ci>Cam</ci>
                  </apply>
                  <ci>c</ci>
               </apply>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <ci>p24</ci>
                     <ci>PSIm</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jmito</ci>
            <apply>
               <minus/>
               <ci>Jnaca</ci>
               <ci>Juni</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jh_leak</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>p17</ci>
                  <ci>PSIm</ci>
               </apply>
               <ci>p18</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>MM1</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>p4</ci>
                  <ci>NADHm</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>p5</ci>
                  <ci>NADHm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>MM2</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>p8</ci>
                  <ci>NADHm</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>p9</ci>
                  <ci>NADHm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jo</ci>
            <apply>
               <divide/>
               <ci>MM1</ci>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>PSIm</ci>
                           <ci>p6</ci>
                        </apply>
                        <ci>p7</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jh_res</ci>
            <apply>
               <divide/>
               <ci>MM2</ci>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>PSIm</ci>
                           <ci>p10</ci>
                        </apply>
                        <ci>p11</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>b13</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>p12</ci>
                  <ci>p13</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>p13</ci>
                  <ci>ATPm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jh_atp</ci>
            <apply>
               <divide/>
               <ci>b13</ci>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>p14</ci>
                           <ci>PSIm</ci>
                        </apply>
                        <ci>p15</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>b2</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>p16</ci>
                  <ci>p13</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>p13</ci>
                  <ci>ATPm</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jf1f0</ci>
            <apply>
               <divide/>
               <ci>b2</ci>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>p14</ci>
                           <ci>PSIm</ci>
                        </apply>
                        <ci>p15</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jant</ci>
            <apply>
               <divide/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>p19</ci>
                     <ci>RATm</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>RATm</ci>
                     <ci>p20</ci>
                  </apply>
               </apply>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.5</cn>
                     </apply>
                     <ci>FRT</ci>
                     <ci>PSIm</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Jgpdh</ci>
            <apply>
               <times/>
               <ci>kgpdh</ci>
               <apply>
                  <root/>
                  <apply>
                     <divide/>
                     <ci>FBP</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="uM">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Parameters">
      <variable initial_value="400" name="p1" public_interface="out" units="dimensionless"/>
      <variable initial_value="1" name="p2" public_interface="out" units="dimensionless"/>
      <variable initial_value="0.01" name="p3" public_interface="out" units="uM"/>
      <variable initial_value="0.6" name="p4" public_interface="out" units="uM_per_ms"/>
      <variable initial_value="0.1" name="p5" public_interface="out" units="mM"/>
      <variable initial_value="177" name="p6" public_interface="out" units="mV"/>
      <variable initial_value="5" name="p7" public_interface="out" units="mV"/>
      <variable initial_value="7" name="p8" public_interface="out" units="uM_per_ms"/>
      <variable initial_value="0.1" name="p9" public_interface="out" units="mM"/>
      <variable initial_value="177" name="p10" public_interface="out" units="mV"/>
      <variable initial_value="5" name="p11" public_interface="out" units="mV"/>
      <variable initial_value="120" name="p12" public_interface="out" units="uM_per_ms"/>
      <variable initial_value="10" name="p13" public_interface="out" units="mM"/>
      <variable initial_value="190" name="p14" public_interface="out" units="mV"/>
      <variable initial_value="8.5" name="p15" public_interface="out" units="mV"/>
      <variable initial_value="35" name="p16" public_interface="out" units="uM_per_ms"/>
      <variable initial_value="0.002" name="p17" public_interface="out" units="uM_per_ms_per_mV"/>
      <variable initial_value="-0.03" name="p18" public_interface="out" units="uM_per_ms"/>
      <variable initial_value="0.35" name="p19" public_interface="out" units="uM_per_ms"/>
      <variable initial_value="2" name="p20" public_interface="out" units="dimensionless"/>
      <variable initial_value="0.01" name="p21" public_interface="out" units="per_uM_per_ms_per_mV"/>
      <variable initial_value="1.1" name="p22" public_interface="out" units="per_uM_per_ms"/>
      <variable initial_value="0.001" name="p23" public_interface="out" units="uM_per_ms"/>
      <variable initial_value="0.016" name="p24" public_interface="out" units="per_mV"/>
   </component>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="J_variables">
         <component_ref component="Parameters"/>
      </component_ref>
      <component_ref component="Clamp_protocol">
         <component_ref component="Clamp_parameters"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="Mitochondrial_variables" component_2="J_variables"/>
      <map_variables variable_1="PSIm" variable_2="PSIm"/>
      <map_variables variable_1="Cam" variable_2="Cam"/>
      <map_variables variable_1="NADHm" variable_2="NADHm"/>
      <map_variables variable_1="Jo" variable_2="Jo"/>
      <map_variables variable_1="Jh_res" variable_2="Jh_res"/>
      <map_variables variable_1="ATPm" variable_2="ATPm"/>
      <map_variables variable_1="Jh_atp" variable_2="Jh_atp"/>
      <map_variables variable_1="Jant" variable_2="Jant"/>
      <map_variables variable_1="Jpdh" variable_2="Jpdh"/>
      <map_variables variable_1="Jf1f0" variable_2="Jf1f0"/>
      <map_variables variable_1="Jnaca" variable_2="Jnaca"/>
      <map_variables variable_1="Jh_leak" variable_2="Jh_leak"/>
      <map_variables variable_1="Juni" variable_2="Juni"/>
      <map_variables variable_1="Jmito" variable_2="Jmito"/>
      <map_variables variable_1="NADm" variable_2="NADm"/>
      <map_variables variable_1="RATm" variable_2="RATm"/>
   </connection>
   <connection>
      <map_components component_1="Parameters" component_2="J_variables"/>
      <map_variables variable_1="p1" variable_2="p1"/>
      <map_variables variable_1="p2" variable_2="p2"/>
      <map_variables variable_1="p3" variable_2="p3"/>
      <map_variables variable_1="p4" variable_2="p4"/>
      <map_variables variable_1="p5" variable_2="p5"/>
      <map_variables variable_1="p6" variable_2="p6"/>
      <map_variables variable_1="p7" variable_2="p7"/>
      <map_variables variable_1="p8" variable_2="p8"/>
      <map_variables variable_1="p9" variable_2="p9"/>
      <map_variables variable_1="p10" variable_2="p10"/>
      <map_variables variable_1="p11" variable_2="p11"/>
      <map_variables variable_1="p12" variable_2="p12"/>
      <map_variables variable_1="p13" variable_2="p13"/>
      <map_variables variable_1="p14" variable_2="p14"/>
      <map_variables variable_1="p15" variable_2="p15"/>
      <map_variables variable_1="p16" variable_2="p16"/>
      <map_variables variable_1="p17" variable_2="p17"/>
      <map_variables variable_1="p18" variable_2="p18"/>
      <map_variables variable_1="p19" variable_2="p19"/>
      <map_variables variable_1="p20" variable_2="p20"/>
      <map_variables variable_1="p21" variable_2="p21"/>
      <map_variables variable_1="p22" variable_2="p22"/>
      <map_variables variable_1="p23" variable_2="p23"/>
      <map_variables variable_1="p24" variable_2="p24"/>
   </connection>
   <connection>
      <map_components component_1="Environment" component_2="Mitochondrial_variables"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Environment" component_2="Clamp_protocol"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Environment" component_2="Cytosol"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Clamp_protocol" component_2="J_variables"/>
      <map_variables variable_1="FBP" variable_2="FBP"/>
      <map_variables variable_1="c" variable_2="c"/>
   </connection>
   <connection>
      <map_components component_1="Clamp_protocol" component_2="Cytosol"/>
      <map_variables variable_1="c" variable_2="c"/>
   </connection>
   <connection>
      <map_components component_1="Cytosol" component_2="J_variables"/>
      <map_variables variable_1="ATPc" variable_2="ATPc"/>
      <map_variables variable_1="Jant" variable_2="Jant"/>
   </connection>
   <connection>
      <map_components component_1="Clamp_protocol" component_2="Clamp_parameters"/>
      <map_variables variable_1="Fproto" variable_2="Fproto"/>
      <map_variables variable_1="Cproto" variable_2="Cproto"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   
   <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
    
	<rdf:Description rdf:about="">
      <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
        <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
          <vCard:Family>Paris</vCard:Family>
          <vCard:Given>Tessa</vCard:Given>
		</vCard:N>
        <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
          <rdf:value>tpar054@aucklanduni.ac.nz</rdf:value>
          <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
        </vCard:EMAIL>
        <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
          <vCard:Orgname>The University of Auckland</vCard:Orgname>
          <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
        </vCard:ORG>
      </dc:creator>
  
      <dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
        <dcterms:W3CDTF>2010-03-22</dcterms:W3CDTF>
      </dcterms:created>
    </rdf:Description>
    
    <rdf:Description rdf:about="#bertram_2006">
      
      <bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
        <dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>electrophysiology</rdf:li>
			  <rdf:li>metabolism</rdf:li>
			  <rdf:li>calcium</rdf:li>
			  <rdf:li>mitochondria</rdf:li>
			  <rdf:li>ATP</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>

      <bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
        <bqs:Pubmed_id>16945388</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/">
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Bertram</vCard:Family>
                    <vCard:Given>R</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Pedersen</vCard:Family>
                    <vCard:Given>M</vCard:Given>
					<vCard:Other>G</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Luciani</vCard:Family>
                    <vCard:Given>D</vCard:Given>
					<vCard:Other>S</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Sherman</vCard:Family>
                    <vCard:Given>A</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
            A simplified model for mitochondrial ATP production
          </dc:title>
          <dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
            <dcterms:W3CDTF>2006-12-26</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Journal of Theoretical Biology</dc:title>
          </bqs:Journal>
		  <bqs:volume>243</bqs:volume>
          <bqs:first_page>575</bqs:first_page>
          <bqs:last_page>586</bqs:last_page>
        </bqs:JournalArticle>
      </bqs:reference>
    </rdf:Description>
  </rdf:RDF>
</model>