Location: Bertram, Pedersen, Luciani and Sherman 2006 @ ebb767eb7c72 / bertram_2006.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2010-03-22 15:19:43+13:00
Desc:
Added metadata to Tessa's model.
Permanent Source URI:
https://models.cellml.org/workspace/bertram_pedersen_luciani_sherman/rawfile/ebb767eb7c7229714b511146cdba6ebd5b498ea9/bertram_2006.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
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: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:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="bertram_2006" cmeta:id="bertram_2006" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
	
	<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 xmlns="http://www.cellml.org/cellml/1.0#" name="mM">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="uM">
      <unit units="mole" prefix="micro"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="mV">
      <unit units="volt" prefix="milli"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="ms">
      <unit units="second" prefix="milli"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_ms">
      <unit units="ms" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="uM_per_ms">
      <unit units="uM"/>
      <unit units="ms" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_uM_per_ms_per_mV">
      <unit units="uM" exponent="-1"/>
      <unit units="ms" exponent="-1"/>
      <unit units="mV" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="uM_per_mV">
      <unit units="uM"/>
      <unit units="mV" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_uM_per_ms">
      <unit units="uM" exponent="-1"/>
      <unit units="ms" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_mV">
      <unit units="mV" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="uM_per_ms_per_mV">
      <unit units="uM"/>
      <unit units="ms" exponent="-1"/>
      <unit units="mV" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="C_per_mol">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="mJ_per_K_per_mol">
      <unit units="joule" prefix="milli"/>
      <unit units="kelvin" exponent="-1"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Environment">
      <variable name="time" units="ms" public_interface="out"/>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Mitochondrial_variables">
      <variable name="NADHm" units="mM" initial_value="0.6" public_interface="out"/>
      <variable name="Jpdh" units="uM_per_ms" public_interface="in"/>
      <variable name="Jo" units="uM_per_ms" public_interface="in"/>
      <variable name="gamma" units="dimensionless" initial_value="0.001"/>
      <variable name="NADm" units="mM" public_interface="out"/>
      <variable name="ADPm" units="mM" initial_value="7.4" public_interface="out"/>
      <variable name="PSIm" units="mV" initial_value="93" public_interface="out"/>
      <variable name="Cam" units="uM" initial_value="0.1" public_interface="out"/>
      <variable name="NADtot" units="mM" initial_value="10"/>
      <variable name="fm" units="dimensionless" initial_value="0.01"/>
      <variable name="Cmito" units="uM_per_mV" initial_value="1.8"/>
      <variable name="Amtot" units="mM" initial_value="15"/>
      <variable name="ATPm" units="mM" public_interface="out"/>
      <variable name="RATm" units="dimensionless" public_interface="out"/>
      <variable name="Jant" units="uM_per_ms" public_interface="in"/>
      <variable name="Jf1f0" units="uM_per_ms" public_interface="in"/>
      <variable name="Jh_res" units="uM_per_ms" public_interface="in"/>
      <variable name="Jh_atp" units="uM_per_ms" public_interface="in"/>
      <variable name="Jh_leak" units="uM_per_ms" public_interface="in"/>
      <variable name="Jnaca" units="uM_per_ms" public_interface="in"/>
      <variable name="Juni" units="uM_per_ms" public_interface="in"/>
      <variable name="time" units="ms" public_interface="in"/>
      <variable name="Jmito" units="uM_per_ms" public_interface="in"/>
      <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 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 xmlns="http://www.cellml.org/cellml/1.0#" name="Cytosol">
      <variable name="time" units="ms" public_interface="in"/>
      <variable name="c" units="uM" public_interface="in"/>
      <variable name="ATPc" units="mM" public_interface="out"/>
      <variable name="ADPc" units="uM" initial_value="1850"/>
      <variable name="Ac_tot" units="uM" initial_value="2500"/>
      <variable name="Jant" units="uM_per_ms" public_interface="in"/>
      <variable name="khyd" units="dimensionless" initial_value="0.00005"/>
      <variable name="Jhydbas" units="mM" initial_value="0.00005"/>
      <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 cellml:units="dimensionless">3.9</cn>
               <cn 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 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 cellml:units="per_ms">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Clamp_protocol">
      <variable name="time" units="ms" public_interface="in" private_interface="out"/>
      <variable name="c" units="uM" initial_value="0.1" public_interface="out"/>
      <variable name="Fproto" units="uM" private_interface="in"/>
      <variable name="Cproto" units="uM" private_interface="in"/>
      <variable name="FBP" units="uM" initial_value="0.5" public_interface="out"/>
      <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 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 cellml:units="ms">0.1</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Clamp_parameters">
      <variable name="time" units="ms" public_interface="in"/>
      <variable name="Fhold" units="uM" initial_value="1"/>
      <variable name="Ftest" units="uM" initial_value="5"/>
      <variable name="Fton" units="ms" initial_value="90000"/>
      <variable name="Ftoff" units="ms" initial_value="330000"/>
      <variable name="Fproto" units="uM" public_interface="out"/>
      <variable name="Cproto" units="uM" public_interface="out"/>
      <variable name="Chold" units="uM" initial_value="0.1"/>
      <variable name="Ctest" units="uM" initial_value="0.1"/>
      <variable name="Cton1" units="ms" initial_value="120000"/>
      <variable name="Cton2" units="ms" initial_value="180000"/>
      <variable name="Cton3" units="ms" initial_value="240000"/>
      <variable name="Ctoff1" units="ms" initial_value="150000"/>
      <variable name="Ctoff2" units="ms" initial_value="210000"/>
      <variable name="Ctoff3" units="ms" initial_value="270000"/>
      <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 cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Fton</ci>
                     </apply>
                     <cn cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_off</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Ftoff</ci>
                     </apply>
                     <cn cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Cton1</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Cton1</ci>
                     </apply>
                     <cn cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Cton2</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Cton2</ci>
                     </apply>
                     <cn cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Cton3</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Cton3</ci>
                     </apply>
                     <cn cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Ctoff1</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Ctoff1</ci>
                     </apply>
                     <cn cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Ctoff2</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Ctoff2</ci>
                     </apply>
                     <cn cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>heav_Ctoff3</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <geq/>
                     <apply>
                        <minus/>
                        <ci>time</ci>
                        <ci>Ctoff3</ci>
                     </apply>
                     <cn cellml:units="ms">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="J_variables">
      <variable name="Juni" units="uM_per_ms" public_interface="out"/>
      <variable name="Jpdh" units="uM_per_ms" public_interface="out"/>
      <variable name="Jnaca" units="uM_per_ms" public_interface="out"/>
      <variable name="Jmito" units="uM_per_ms" public_interface="out"/>
      <variable name="Jh_leak" units="uM_per_ms" public_interface="out"/>
      <variable name="Jh_atp" units="uM_per_ms" public_interface="out"/>
      <variable name="Jo" units="uM_per_ms" public_interface="out"/>
      <variable name="Jh_res" units="uM_per_ms" public_interface="out"/>
      <variable name="Jf1f0" units="uM_per_ms" public_interface="out"/>
      <variable name="Jant" units="uM_per_ms" public_interface="out"/>
      <variable name="Jgpdh" units="uM_per_ms"/>
      <variable name="c" units="uM" public_interface="in"/>
      <variable name="PSIm" units="mV" public_interface="in"/>
      <variable name="Cam" units="uM" public_interface="in"/>
      <variable name="NADHm" units="mM" public_interface="in"/>
      <variable name="p1" units="dimensionless" private_interface="in"/>
      <variable name="p2" units="dimensionless" private_interface="in"/>
      <variable name="p3" units="uM" private_interface="in"/>
      <variable name="p21" units="per_uM_per_ms_per_mV" private_interface="in"/>
      <variable name="p22" units="per_uM_per_ms" private_interface="in"/>
      <variable name="p23" units="uM_per_ms" private_interface="in"/>
      <variable name="p24" units="per_mV" private_interface="in"/>
      <variable name="p4" units="uM_per_ms" private_interface="in"/>
      <variable name="p5" units="mM" private_interface="in"/>
      <variable name="p6" units="mV" private_interface="in"/>
      <variable name="p7" units="mV" private_interface="in"/>
      <variable name="p8" units="uM_per_ms" private_interface="in"/>
      <variable name="p9" units="mM" private_interface="in"/>
      <variable name="p10" units="mV" private_interface="in"/>
      <variable name="p11" units="mV" private_interface="in"/>
      <variable name="p12" units="uM_per_ms" private_interface="in"/>
      <variable name="p13" units="mM" private_interface="in"/>
      <variable name="p14" units="mV" private_interface="in"/>
      <variable name="p15" units="mV" private_interface="in"/>
      <variable name="p16" units="uM_per_ms" private_interface="in"/>
      <variable name="p17" units="uM_per_ms_per_mV" private_interface="in"/>
      <variable name="p18" units="uM_per_ms" private_interface="in"/>
      <variable name="p19" units="uM_per_ms" private_interface="in"/>
      <variable name="p20" units="dimensionless" private_interface="in"/>
      <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" units="mM" public_interface="in"/>
      <variable name="ATPm" units="mM" public_interface="in"/>
      <variable name="ATPc" units="mM" public_interface="in"/>
      <variable name="RATm" units="dimensionless" public_interface="in"/>
      <variable name="FRT" units="per_mV" initial_value="0.037"/>
      <variable name="kgpdh" units="uM_per_ms" initial_value="0.0005"/>
      <variable name="FBP" units="uM" public_interface="in"/>
      <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 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 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 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 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 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 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 cellml:units="uM">1</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Parameters">
      <variable name="p1" units="dimensionless" initial_value="400" public_interface="out"/>
      <variable name="p2" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="p3" units="uM" initial_value="0.01" public_interface="out"/>
      <variable name="p4" units="uM_per_ms" initial_value="0.6" public_interface="out"/>
      <variable name="p5" units="mM" initial_value="0.1" public_interface="out"/>
      <variable name="p6" units="mV" initial_value="177" public_interface="out"/>
      <variable name="p7" units="mV" initial_value="5" public_interface="out"/>
      <variable name="p8" units="uM_per_ms" initial_value="7" public_interface="out"/>
      <variable name="p9" units="mM" initial_value="0.1" public_interface="out"/>
      <variable name="p10" units="mV" initial_value="177" public_interface="out"/>
      <variable name="p11" units="mV" initial_value="5" public_interface="out"/>
      <variable name="p12" units="uM_per_ms" initial_value="120" public_interface="out"/>
      <variable name="p13" units="mM" initial_value="10" public_interface="out"/>
      <variable name="p14" units="mV" initial_value="190" public_interface="out"/>
      <variable name="p15" units="mV" initial_value="8.5" public_interface="out"/>
      <variable name="p16" units="uM_per_ms" initial_value="35" public_interface="out"/>
      <variable name="p17" units="uM_per_ms_per_mV" initial_value="0.002" public_interface="out"/>
      <variable name="p18" units="uM_per_ms" initial_value="-0.03" public_interface="out"/>
      <variable name="p19" units="uM_per_ms" initial_value="0.35" public_interface="out"/>
      <variable name="p20" units="dimensionless" initial_value="2" public_interface="out"/>
      <variable name="p21" units="per_uM_per_ms_per_mV" initial_value="0.01" public_interface="out"/>
      <variable name="p22" units="per_uM_per_ms" initial_value="1.1" public_interface="out"/>
      <variable name="p23" units="uM_per_ms" initial_value="0.001" public_interface="out"/>
      <variable name="p24" units="per_mV" initial_value="0.016" public_interface="out"/>
   </component>
   <group xmlns="http://www.cellml.org/cellml/1.0#">
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Environment" component_2="Mitochondrial_variables"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Environment" component_2="Clamp_protocol"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Environment" component_2="Cytosol"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Clamp_protocol" component_2="Cytosol"/>
      <map_variables variable_1="c" variable_2="c"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <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#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
    
	<rdf:Description rdf:about="">
      <dc:creator rdf:parseType="Resource">
        <vCard:N rdf:parseType="Resource">
          <vCard:Family>Paris</vCard:Family>
          <vCard:Given>Tessa</vCard:Given>
		</vCard:N>
        <vCard:EMAIL 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 rdf:parseType="Resource">
          <vCard:Orgname>The University of Auckland</vCard:Orgname>
          <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
        </vCard:ORG>
      </dc:creator>
  
      <dcterms:created rdf:parseType="Resource">
        <dcterms:W3CDTF>2010-03-22</dcterms:W3CDTF>
      </dcterms:created>
    </rdf:Description>
    
    <rdf:Description rdf:about="#bertram_2006">
      
      <bqs:reference rdf:parseType="Resource">
        <dc:subject rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>electrophysiology</rdf:li>
			  <rdf:li>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 rdf:parseType="Resource">
        <bqs:Pubmed_id>16945388</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator>
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N 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 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 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 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>
            A simplified model for mitochondrial ATP production
          </dc:title>
          <dcterms:issued rdf:parseType="Resource">
            <dcterms:W3CDTF>2006-12-26</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title>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>