Rendering of the source text

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  pedersen_model_2005.xml

CREATED :  27th July 2007

LAST MODIFIED : 27th July 2007

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.1 Specification.

DESCRIPTION :  This file contains a CellML description of Pedersen, Bertram and Sherman's 2005 mathematical model of metabolically driven insulin secretion in pancreatic beta cells.

CHANGES:  
  
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="pedersen_model_2005" name="pedersen_bertram_sherman_model_2005_version01">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Intra- and inter-islet synchronization of metabolically driven insulin secretion</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
Insulin secretion from pancreatic beta-cells is known to be pulsatile, with a period of 5-10 milliseconds.  It has been suggested that this pulsatile behaviour stems from oscillations in the glycolysis pathway, mediated by the enzyme phosphofructokinase (PFK).  In turn, this rhythmical enzymatic behaviour induces periodic activity of ATP-dependent potassium channels (K(ATP)-channels) in the plasma membrane.  Insulin pulsatility is impaired in diabetic patients, which are also relatively insulin insensitive.  Target tissues, such as the liver, are more sensitive to pulsatile insulin than to constant levels, and therefore an improved understanding of the mechanisms underlying pulsatile insulin secretion is important in the search for a potential therapeutic treatment for diabetes.   
</para>

<para>
The link between metabolism and calcium ion influx leading to insulin secretion is provided by the electrical activity of the beta-cells, which has a characteristic behaviour known as <emphasis>bursting</emphasis>.  A burst consists of an active phase of spiking followed by a silent phase of membrane repolarisation.  During the active burst calcium ions diffuse into the cell through voltage-gated calcium channels, and the subsequent rise in the cytosolic calcium concentration results in the secretion of insulin (see diagram below). 
</para>

<para>
In the paper described here, Pedersen <emphasis>et al.</emphasis> present a mathematical model of metabolically driven insulin secretion.  This model is based on one which has been published previously by Bertram <emphasis>et al.</emphasis> (2004), which has also been described in CellML and can be found in the model repository.  Model simulation data indicate that electrical coupling is sufficient to synchronise electrical bursting activity and metabolic oscillations.
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="pedersen_2005.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the pathways described by the mathematical model.  Extracellular glucose (Ge) enters the pancreatic beta-cell through GLUT-2 transporters in the plasma membrane, and is subsequently broken down through a series of enzyme-catalysed reactions in the glycolysis pathway.  The metabolic products of the glycolysis pathway feed into the mitchondria, the site of aerobic respiration, where they are used to produce ATP.  ATP represents the common substrate, linking the metabolic (glycolysis) and electrical (membrane potential) components of the mathematical model by regulating the flow of potassium ions through the K(ATP)-channels in the plasma membrane.  In turn, these control regulate the membrane potential and calcium ion flow leading to insulin secretion.  Insulin acts to lower the plasma glucose concentration though the functions of the liver (represented by a dashed line).</caption>
</informalfigure>

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

<para>
<ulink url="http://www.biophysj.org/cgi/content/abstract/89/1/107">Intra- and inter-islet synchronization of metabolically driven insulin secretion</ulink>, Morten Gram Pedersen, Richard Bertram, and Arthur Sherman, 2005, <ulink url="http://www.biophysj.org/">
            <emphasis>Biophysical Journal</emphasis>
          </ulink>, 89, 107-119.  (<ulink url="http://www.biophysj.org/cgi/content/full/89/1/107A">Full text</ulink> and <ulink url="http://www.biophysj.org/cgi/reprint/89/1/107">PDF</ulink> versions of the article are available to journal subscribers on the <emphasis>Biophysical Journal</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=15834002&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 15834002</ulink>
</para>

<para>
The authors highlight that the original model code can be downloaded <ulink url="http://mrb.niddk.nih.gov/sherman">here</ulink>.  The CellML translation of this model represents the core model and it will run in the PCEnv software.  However the pulsatile behaviour of the model, as published in the original paper can not be recreated.  To achieve these results I suspect several cell models have to be run simultaneously, generating many bursts which are clustered together into what is termed <emphasis>compound bursting</emphasis>.
</para>

</sect1>
</article>
</documentation> 





  <units name="millisecond">
    <unit units="second" prefix="milli"/>
  </units>

  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>

  <units name="femtofarad">
    <unit units="farad" prefix="femto"/>
  </units>

  <units name="picosiemens">
    <unit units="siemens" prefix="pico"/>
  </units>

  <units name="picoampere">
    <unit units="ampere" prefix="pico"/>
  </units>

  <units name="millimolar">
    <unit units="mole" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="per_micromolar">
    <unit units="micromolar" exponent="-1"/>
  </units>

  <units name="millimolar_millisecond">
    <unit units="millimolar"/>
    <unit units="millisecond" exponent="-1"/>
  </units>

  <units name="micromolar_millisecond">
    <unit units="micromolar"/>
    <unit units="millisecond" exponent="-1"/>
  </units>

  <units name="first_order_rate_constant">
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  
  <component name="environment">
    <variable units="millisecond" public_interface="out" name="time"/>
  </component>
  
  <component cmeta:id="Gi" name="Gi">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="Gi">
        <dc:title>Gi</dc:title>
        <dcterms:alternative>intracellular glucose</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="millimolar" public_interface="out" name="Gi"/>
    
    <variable units="millimolar_millisecond" public_interface="in" name="Jglut"/>
    <variable units="millimolar_millisecond" public_interface="in" name="Jgk"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> Gi </ci> 
        </apply>
        <apply>
          <minus/>
          <ci> Jglut </ci>
          <ci> Jgk </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="Ge" name="Ge">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="Ge">
        <dc:title>Ge</dc:title>
        <dcterms:alternative>extracellular glucose</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="millimolar" public_interface="out" name="Ge" initial_value="7.0"/>
    
    <variable units="millimolar" name="G_infinity"/>
    <variable units="millisecond" name="tau_G" initial_value="450000"/>
    <variable units="millimolar" name="Gmax" initial_value="15.0"/>
    <variable units="millimolar" name="Gmin" initial_value="1.0"/>
    <variable units="dimensionless" name="SG" initial_value="1.0"/>
    <variable units="dimensionless" name="I_" initial_value="5.0"/>
    
    <variable units="dimensionless" public_interface="in" name="I"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> Ge </ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> G_infinity </ci>
            <ci> Ge </ci>
          </apply>
          <ci> tau_G </ci>
        </apply>
      </apply> 

      <apply>
        <eq/>
        <ci> G_infinity </ci>
        <apply>
          <plus/>
          <ci> Gmin </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Gmax </ci>
              <ci> Gmin </ci>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> I </ci>
                    <ci> I_ </ci>
                  </apply>
                  <ci> SG </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="G6P" name="G6P">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="G6P">
        <dc:title>G6P</dc:title>
        <dcterms:alternative>glucose 6-phosphate</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="millimolar" public_interface="out" name="G6P" initial_value="200.0"/>
    
    <variable units="dimensionless" public_interface="in" name="kappa"/>
    <variable units="micromolar_millisecond" public_interface="in" name="JPFK"/>
    <variable units="millimolar_millisecond" public_interface="in" name="Jgk"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> G6P </ci> 
        </apply>
        <apply>
          <times/>
          <ci> kappa </ci>
          <apply>
            <minus/>
            <ci> Jgk </ci>
            <ci> JPFK </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="FBP" name="FBP">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="FBP">
        <dc:title>FBP</dc:title>
        <dcterms:alternative>fructose 1-6-bisphosphate</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar" public_interface="out" name="FBP" initial_value="40.0"/>
    
    <variable units="dimensionless" public_interface="in" name="kappa"/>
    <variable units="micromolar_millisecond" public_interface="in" name="JPFK"/>
    <variable units="micromolar_millisecond" public_interface="in" name="JGPDH"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> FBP </ci> 
        </apply>
        <apply>
          <times/>
          <ci> kappa </ci>
          <apply>
            <minus/>
            <ci> JPFK </ci>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.5 </cn>
              <ci> JGPDH </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="F6P" name="F6P">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="F6P">
        <dc:title>F6P</dc:title>
        <dcterms:alternative>fructose 6-phosphate</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="millimolar" public_interface="out" name="F6P"/>
    
    <variable units="millimolar" public_interface="in" name="G6P"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> F6P </ci> 
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.3 </cn>
          <ci> G6P </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="Jglut" name="Jglut">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="Jglut">
        <dc:title>Jglut</dc:title>
        <dcterms:alternative>rate of GLUT-2 facilitated glucose transport</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="millimolar_millisecond" public_interface="out" name="Jglut"/>
    
    <variable units="millimolar" name="Kglut" initial_value="7.0"/>
    <variable units="millimolar_millisecond" name="Vglut" initial_value="8.0"/>
    
    <variable units="millimolar" public_interface="in" name="Gi"/>
    <variable units="millimolar" public_interface="in" name="Ge"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> Jglut </ci> 
        <apply>
          <times/>
          <ci> Vglut </ci>
          <apply>
            <divide/>
            <apply>
              <times/>
              <apply>
                <minus/>
                <ci> Ge </ci>
                <ci> Gi </ci>
              </apply>
              <ci> Kglut </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <plus/>
                <ci> Kglut </ci>
                <ci> Ge </ci>
              </apply>
              <apply>
                <minus/>
                <ci> Kglut </ci>
                <ci> Gi </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="Jgk" name="Jgk">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="Jgk">
        <dc:title>Jgk</dc:title>
        <dcterms:alternative>rate of the glucokinase reaction</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="millimolar_millisecond" public_interface="out" name="Jgk"/>
    
    <variable units="millimolar" name="Kgk" initial_value="7.0"/>
    <variable units="millimolar_millisecond" name="Vgk" initial_value="0.8"/>
    <variable units="dimensionless" name="ngk" initial_value="4.0"/>
    
    <variable units="millimolar" public_interface="in" name="Gi"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> Jgk </ci> 
        <apply>
          <times/>
          <ci> Vgk </ci>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> Gi </ci>
              <ci> ngk </ci>
            </apply>
            <apply>
              <plus/>
              <apply>
                <power/>
                <ci> Kgk </ci>
                <ci> ngk </ci>
              </apply>
              <apply>
                <power/>
                <ci> Gi </ci>
                <ci> ngk </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="JGPDH" name="JGPDH">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="JGPDH">
        <dc:title>JGPDH</dc:title>
        <dcterms:alternative>rate of the glyceride 3-P dehydrogenase reaction</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar_millisecond" public_interface="out" name="JGPDH"/>
    
    <variable units="micromolar" public_interface="in" name="FBP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> JGPDH </ci> 
        <apply>
          <times/>
          <cn cellml:units="micromolar_millisecond"> 0.2 </cn>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> FBP </ci>
              <cn cellml:units="micromolar"> 1.0 </cn>
            </apply>
            <cn cellml:units="dimensionless"> 0.5 </cn>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component cmeta:id="JPFK" name="JPFK">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="JPFK">
        <dc:title>JPFK</dc:title>
        <dcterms:alternative>rate of the phosphofructokinase reaction</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar_millisecond" public_interface="out" name="JPFK"/>
    
    <variable units="dimensionless" name="lambda" initial_value="0.06"/>
    <variable units="micromolar_millisecond" name="Vmax" initial_value="2.0"/>
    
    <variable units="dimensionless" public_interface="in" name="w0000"/>
    <variable units="dimensionless" public_interface="in" name="w1000"/>
    <variable units="dimensionless" public_interface="in" name="w0100"/>
    <variable units="dimensionless" public_interface="in" name="w0010"/>
    <variable units="dimensionless" public_interface="in" name="w0001"/>
    <variable units="dimensionless" public_interface="in" name="w1100"/>
    <variable units="dimensionless" public_interface="in" name="w1010"/>
    <variable units="dimensionless" public_interface="in" name="w1001"/>
    <variable units="dimensionless" public_interface="in" name="w0101"/>
    <variable units="dimensionless" public_interface="in" name="w0011"/>
    <variable units="dimensionless" public_interface="in" name="w0110"/>
    <variable units="dimensionless" public_interface="in" name="w1110"/>
    <variable units="dimensionless" public_interface="in" name="w0111"/>
    <variable units="dimensionless" public_interface="in" name="w1011"/>
    <variable units="dimensionless" public_interface="in" name="w1101"/>
    <variable units="dimensionless" public_interface="in" name="w1111"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> JPFK </ci> 
        <apply>
          <times/>
          <ci> Vmax </ci>
          <apply>
            <divide/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <ci> lambda </ci>
                </apply>
                <ci> w1110 </ci>
              </apply>
              <apply>
                <times/>
                <ci> lambda </ci>
                <apply>
                  <plus/>
                  <ci> w0010 </ci>
                  <ci> w1010 </ci>
                  <ci> w0011 </ci>
                  <ci> w1110 </ci>
                  <ci> w0111 </ci>
                  <ci> w1011 </ci>
                  <ci> w1101 </ci>
                  <ci> w1111 </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <plus/>
              <ci> w0000 </ci>
              <ci> w1000 </ci>
              <ci> w0100 </ci>
              <ci> w0010 </ci>
              <ci> w0001 </ci>
              <ci> w1100 </ci>
              <ci> w1010 </ci>
              <ci> w1001 </ci>
              <ci> w0101 </ci>
              <ci> w0011 </ci>
              <ci> w0110 </ci>
              <ci> w1110 </ci>
              <ci> w0111 </ci>
              <ci> w1011 </ci>
              <ci> w1101 </ci>
              <ci> w1111 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="w">  
    <variable units="dimensionless" public_interface="out" name="w0000"/>
    <variable units="dimensionless" public_interface="out" name="w1000"/>
    <variable units="dimensionless" public_interface="out" name="w0100"/>
    <variable units="dimensionless" public_interface="out" name="w0010"/>
    <variable units="dimensionless" public_interface="out" name="w0001"/>
    <variable units="dimensionless" public_interface="out" name="w1100"/>
    <variable units="dimensionless" public_interface="out" name="w1010"/>
    <variable units="dimensionless" public_interface="out" name="w1001"/>
    <variable units="dimensionless" public_interface="out" name="w0101"/>
    <variable units="dimensionless" public_interface="out" name="w0011"/>
    <variable units="dimensionless" public_interface="out" name="w1110"/>
    <variable units="dimensionless" public_interface="out" name="w0111"/>
    <variable units="dimensionless" public_interface="out" name="w1011"/>
    <variable units="dimensionless" public_interface="out" name="w1101"/>
    <variable units="dimensionless" public_interface="out" name="w1111"/>
    <variable units="dimensionless" public_interface="out" name="w0110"/>
    
    <variable units="dimensionless" name="f13" initial_value="0.02"/>
    <variable units="dimensionless" name="f23" initial_value="0.2"/>
    <variable units="dimensionless" name="f41" initial_value="20.0"/>
    <variable units="dimensionless" name="f42" initial_value="20.0"/>
    <variable units="dimensionless" name="f43" initial_value="20.0"/>
    <variable units="micromolar" name="K1" initial_value="30.0"/>
    <variable units="micromolar" name="K2" initial_value="1.0"/>
    <variable units="micromolar" name="K3" initial_value="50000.0"/>
    <variable units="micromolar" name="K4" initial_value="1000.0"/>
    
    <variable units="micromolar" public_interface="in" name="AMP"/>
    <variable units="micromolar" public_interface="in" name="FBP"/>
    <variable units="millimolar" public_interface="in" name="F6P"/>
    <variable units="micromolar" public_interface="in" name="ATP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> w0000 </ci> 
        <cn cellml:units="dimensionless"> 1.0 </cn>
      </apply>

      <apply>
        <eq/> 
        <ci> w1000 </ci> 
        <apply>
          <power/>
          <apply>
            <divide/>
            <ci> AMP </ci>
            <ci> K1 </ci>
          </apply>
          <cn cellml:units="dimensionless"> 1.0 </cn>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w0100 </ci>
        <apply>
          <power/>
          <apply>
            <divide/>
            <ci> FBP </ci>
            <ci> K2 </ci>
          </apply>
          <cn cellml:units="dimensionless"> 1.0 </cn>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w0010 </ci>
        <apply>
          <power/>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> F6P </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <ci> K3 </ci>
          </apply>
          <cn cellml:units="dimensionless"> 1.0 </cn>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w0001 </ci>
        <apply>
          <power/>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> ATP </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <ci> K4 </ci>
          </apply>
          <cn cellml:units="dimensionless"> 1.0 </cn>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w1100 </ci>
        <apply>
          <times/>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> AMP </ci>
              <ci> K1 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> FBP </ci>
              <ci> K2 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>        

      <apply>
        <eq/> 
        <ci> w1010 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <ci> f13 </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <cn cellml:units="dimensionless"> 1.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> AMP </ci>
              <ci> K1 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> F6P </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K3 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w1001 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <ci> f41 </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <cn cellml:units="dimensionless"> 1.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> AMP </ci>
              <ci> K1 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> ATP </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K4 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w0101 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <ci> f42 </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <cn cellml:units="dimensionless"> 1.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> FBP </ci>
              <ci> K2 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> ATP </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K4 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w0011 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <ci> f43 </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <cn cellml:units="dimensionless"> 1.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> F6P </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K3 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> ATP </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K4 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w0110 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <ci> f23 </ci>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <cn cellml:units="dimensionless"> 1.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> FBP </ci>
              <ci> K2 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> F6P </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K3 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w1110 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <ci> f13 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f23 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> AMP </ci>
              <ci> K1 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> FBP </ci>
              <ci> K2 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> F6P </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K3 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply> 
   
      <apply>
        <eq/> 
        <ci> w0111 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <ci> f23 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f42 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f43 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> FBP </ci>
              <ci> K2 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> F6P </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K3 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> ATP </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K4 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w1011 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <ci> f13 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f41 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f43 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> AMP </ci>
              <ci> K1 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> F6P </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K3 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> ATP </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K4 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w1101 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <ci> f41 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f42 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> AMP </ci>
              <ci> K1 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> FBP </ci>
              <ci> K2 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> ATP </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K4 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> w1111 </ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <times/>
              <apply>
                <power/>
                <ci> f13 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f23 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f41 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f42 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
              <apply>
                <power/>
                <ci> f43 </ci>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> AMP </ci>
              <ci> K1 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <ci> FBP </ci>
              <ci> K2 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> F6P </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K3 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
          <apply>
            <power/>
            <apply>
              <divide/>
              <apply>
                <power/>
                <ci> ATP </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <ci> K4 </ci>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>                                                                                                 
    </math>
  </component>

  <component cmeta:id="ATP" name="ATP">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="ATP">
        <dc:title>ATP</dc:title>
        <dcterms:alternative>adenosine triphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar" public_interface="out" name="ATP"/>
    
    <variable units="micromolar" name="Atot" initial_value="3000.0"/>
    
	<variable units="micromolar" public_interface="in" name="ADP"/>
	<variable units="micromolar" public_interface="in" name="AMP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> ATP </ci>
        <apply>
          <minus/>
          <ci> Atot </ci>
          <apply>
            <plus/>
            <ci> ADP </ci>
			<ci> AMP </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="AMP" name="AMP">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="AMP">
        <dc:title>AMP</dc:title>
        <dcterms:alternative>adenosine monophosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar" public_interface="out" name="AMP"/>
    
    <variable units="micromolar" public_interface="in" name="ATP"/>
    <variable units="micromolar" public_interface="in" name="ADP"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> AMP </ci> 
        <apply>
          <divide/>
          <apply>
            <power/>
            <ci> ADP </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <ci> ATP </ci>
        </apply>
      </apply>
    </math>
  </component>


  <component cmeta:id="ADP" name="ADP">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="ADP">
        <dc:title>ADP</dc:title>
        <dcterms:alternative>adenosine diphosphate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar" public_interface="out" name="ADP" initial_value="780.0"/>
    
    <variable units="millisecond" name="tau_a" initial_value="300000.0"/>
    <variable units="dimensionless" name="r" initial_value="0.5"/>
    <variable units="micromolar" name="r1" initial_value="0.35"/>
    
    <variable units="dimensionless" public_interface="in" name="gamma"/>
    <variable units="micromolar" public_interface="in" name="Ca"/>
    <variable units="micromolar" public_interface="in" name="ATP"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> ADP </ci> 
        </apply>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> tau_a </ci>
          </apply>
          <apply>
            <minus/>
            <ci> ATP </ci>
            <apply>
              <times/>
              <ci> ADP </ci>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <apply>
                    <plus/>
                    <ci> r </ci>
                    <ci> gamma </ci>
                  </apply>
                  <apply>
                    <minus/>
                    <cn cellml:units="dimensionless"> 1.0 </cn>
                    <apply>
                      <divide/>
                      <ci> Ca </ci>
                      <ci> r1 </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="gamma">  
    <variable units="dimensionless" public_interface="out" name="gamma"/>
    
    <variable units="dimensionless" name="v_gamma" initial_value="2.2"/>
    <variable units="micromolar_millisecond" name="k_gamma" initial_value="0.1"/>
    
    <variable units="micromolar_millisecond" public_interface="in" name="JGPDH"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> gamma </ci> 
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> v_gamma </ci>
            <ci> JGPDH </ci>
          </apply>
          <apply>
            <plus/>
            <ci> k_gamma </ci>
            <ci> JGPDH </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="membrane">
    <variable units="millivolt" public_interface="out" name="v" initial_value="-60.0"/>
    
    <variable units="femtofarad" name="cm" initial_value="5300.0"/>
    
    <variable units="picoampere" public_interface="in" name="I_Ca"/>
    <variable units="picoampere" public_interface="in" name="I_K"/>
    <variable units="picoampere" public_interface="in" name="I_K_Ca"/>
    <variable units="picoampere" public_interface="in" name="I_K_ATP"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> v </ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <ci> I_K </ci>
              <ci> I_Ca </ci>
              <ci> I_K_Ca </ci>
              <ci> I_K_ATP </ci>
            </apply>
          </apply>
          <ci> cm </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="I_K" name="I_K">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="I_K">
        <dc:title>I_K</dc:title>
        <dcterms:alternative>voltage-dependent potassium current</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="picoampere" public_interface="out" name="I_K"/>
    
    <variable units="picosiemens" name="gK" initial_value="2700.0"/>   
   
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" name="vK"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="v"/>
    <variable units="dimensionless" private_interface="in" name="n"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci> I_K </ci>
        <apply>
          <times/>
          <ci> gK </ci>
          <ci> n </ci>
          <apply>
            <minus/>
            <ci> v </ci>
            <ci> vK </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="n" name="n">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="n">
        <dc:title>n</dc:title>
        <dcterms:alternative>voltage-dependent potassium current n gate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="dimensionless" public_interface="out" name="n"/>
    
    <variable units="dimensionless" name="n_infinity"/>
    <variable units="millisecond" name="tau_n"/>   
   
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="v"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> n </ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> n_infinity </ci>
            <ci> n </ci>
          </apply>
          <ci> tau_n </ci>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> tau_n </ci> 
        <apply>
          <divide/>
          <cn cellml:units="millisecond"> 1.0 </cn>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.035 </cn>
            <apply>
              <cosh/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> v </ci>
                  <cn cellml:units="millivolt"> 16.0 </cn>
                </apply>
                <cn cellml:units="millivolt"> 22.4 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> n_infinity </ci> 
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci> v </ci>
                    <cn cellml:units="millivolt"> 16.0 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 5.6 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="I_Ca" name="I_Ca">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="I_Ca">
        <dc:title>I_Ca</dc:title>
        <dcterms:alternative>voltage-dependent calcium current</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="picoampere" public_interface="out" name="I_Ca"/>
    
    <variable units="picosiemens" name="gCa" initial_value="1000.0"/>   
   
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="vCa"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="v"/>
    <variable units="dimensionless" private_interface="in" name="m_infinity"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci> I_Ca </ci>
        <apply>
          <times/>
          <ci> gCa </ci>
          <ci> m_infinity </ci>
          <apply>
            <minus/>
            <ci> v </ci>
            <ci> vCa </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="m" name="m">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="m">
        <dc:title>m</dc:title>
        <dcterms:alternative>voltage-dependent calcium current m gate</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="dimensionless" public_interface="out" name="m_infinity"/>
    
    <variable units="millivolt" public_interface="in" name="v"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/> 
        <ci> m_infinity </ci> 
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci> v </ci>
                    <cn cellml:units="millivolt"> 20.0 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 12.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="I_K_Ca" name="I_K_Ca">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="I_K_Ca">
        <dc:title>I_K_Ca</dc:title>
        <dcterms:alternative>calcium-activated potassium current</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="picoampere" public_interface="out" name="I_K_Ca"/>
    
    <variable units="picosiemens" name="gkCa_" initial_value="400.0"/>
    <variable units="picosiemens" name="gkCa"/>
    <variable units="micromolar" name="KD" initial_value="0.5"/>   
   
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="vK"/>
    <variable units="micromolar" public_interface="in" name="Ca"/>
    <variable units="millivolt" public_interface="in" name="v"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci> I_K_Ca </ci>
        <apply>
          <times/>
          <ci> gkCa </ci>
          <apply>
            <minus/>
            <ci> v </ci>
            <ci> vK </ci>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/>
        <ci> gkCa </ci>
        <apply>
          <times/>
          <ci> gkCa_ </ci>
          <apply>
            <divide/>
            <apply>
              <power/>
              <ci> Ca </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <plus/>
              <apply>
                <power/>
                <ci> KD </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
              <apply>
                <power/>
                <ci> Ca </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="I_K_ATP" name="I_K_ATP">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="I_K_ATP">
        <dc:title>I_K_ATP</dc:title>
        <dcterms:alternative>ATP-sensitive potassium current</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="picoampere" public_interface="out" name="I_K_ATP"/>
    
    <variable units="picosiemens" name="gkATP_" initial_value="40000.0"/>
    <variable units="picosiemens" name="gkATP"/>
    <variable units="dimensionless" name="o_infinity"/>
    <variable units="micromolar" name="MgADP"/> 
    <variable units="micromolar" name="ADP3_"/> 
    <variable units="micromolar" name="ATP4_"/>   
   
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="vK"/>
    <variable units="micromolar" public_interface="in" name="ATP"/>
    <variable units="micromolar" public_interface="in" name="ADP"/>
    <variable units="millivolt" public_interface="in" name="v"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <ci> I_K_ATP </ci>
        <apply>
          <times/>
          <ci> gkATP </ci>
          <apply>
            <minus/>
            <ci> v </ci>
            <ci> vK </ci>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/>
        <ci> gkATP </ci>
        <apply>
          <times/>
          <ci> gkATP_ </ci>
          <ci> o_infinity </ci>
        </apply>
      </apply>
          
      <apply>
        <eq/>
        <ci> o_infinity </ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.08 </cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> MgADP </ci>
                  </apply>
                  <cn cellml:units="micromolar"> 17.0 </cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.89 </cn>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> MgADP </ci>
                  <cn cellml:units="micromolar"> 17.0 </cn>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <apply>
              <power/>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> MgADP </ci>
                  </apply>
                  <cn cellml:units="micromolar"> 17.0 </cn>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <divide/>
                <ci> ADP3_ </ci>
                <cn cellml:units="micromolar"> 26.0 </cn>
              </apply>
              <apply>
                <divide/>
                <ci> ATP4_ </ci>
                <cn cellml:units="micromolar"> 1.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/>
        <ci> MgADP </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.165 </cn>
          <ci> ADP </ci>
        </apply>
      </apply>

      <apply>
        <eq/>
        <ci> ADP3_ </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.135 </cn>
          <ci> ADP </ci>
        </apply>
      </apply>

      <apply>
        <eq/>
        <ci> ATP4_ </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 0.005 </cn>
          <ci> ATP </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="Ca" name="Ca">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="Ca">
        <dc:title>Ca</dc:title>
        <dcterms:alternative>calcium</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar" public_interface="out" name="Ca" initial_value="0.1"/>

    <variable units="dimensionless" name="fcyt" initial_value="0.01"/>
    
    <variable units="micromolar_millisecond" public_interface="in" name="Jmem"/>
    <variable units="micromolar_millisecond" public_interface="in" name="Jer"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> Ca </ci> 
        </apply>
        <apply>
          <times/>
          <ci> fcyt </ci>
          <apply>
            <plus/>
            <ci> Jmem </ci>
            <ci> Jer </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="Caer" name="Caer">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="Caer">
        <dc:title>Caer</dc:title>
        <dcterms:alternative>calcium</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar" public_interface="out" name="Caer" initial_value="185.0"/>

    <variable units="dimensionless" name="fer" initial_value="0.01"/>
    <variable units="dimensionless" name="Vcyt_Ver" initial_value="31.0"/>
    
    <variable units="micromolar_millisecond" public_interface="in" name="Jer"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> Caer </ci> 
        </apply>
        <apply>
          <times/>
          <apply>
            <minus/>
            <ci> fer </ci>
          </apply>
          <ci> Vcyt_Ver </ci>
          <ci> Jer </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="Jmem" name="Jmem">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="Jmem">
        <dc:title>Jmem</dc:title>
        <dcterms:alternative>rate of calcium flux across the plasma membrane</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar_millisecond" public_interface="out" name="Jmem"/>
    
    <variable units="first_order_rate_constant" name="kPMCA" initial_value="0.18"/>
    <variable units="micromolar_millisecond" name="alpha" initial_value="4.5E-6"/>
    
    <variable units="micromolar" public_interface="in" name="Ca"/>
    <variable units="picoampere" public_interface="in" name="I_Ca"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> Jmem </ci>
        <apply>
          <minus/>
          <apply>
            <plus/>  
            <apply>
              <times/>
              <ci> alpha </ci>
              <ci> I_Ca </ci>
            </apply>
            <apply>
              <times/>
              <ci> kPMCA </ci>
              <ci> Ca </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="Jer" name="Jer">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="Jer">
        <dc:title>Jer</dc:title>
        <dcterms:alternative>rate of calcium flux out of the endoplasmic reticulum</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar_millisecond" public_interface="out" name="Jer"/>
    
    <variable units="micromolar_millisecond" public_interface="in" name="Jleak"/>
    <variable units="micromolar_millisecond" public_interface="in" name="JSERCA"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> Jer </ci> 
        <apply>
          <minus/>
          <ci> Jleak </ci>
          <ci> JSERCA </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="JSERCA" name="JSERCA">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="JSERCA">
        <dc:title>JSERCA</dc:title>
        <dcterms:alternative>rate of calcium flux into the endoplasmic reticulum through the SERCA pumps</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar_millisecond" public_interface="out" name="JSERCA"/>
    
    <variable units="first_order_rate_constant" name="kSERCA" initial_value="0.4"/>
    
    <variable units="micromolar" public_interface="in" name="Ca"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> JSERCA </ci> 
        <apply>
          <times/>
          <ci> kSERCA </ci>
          <ci> Ca </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="Jleak" name="Jleak">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="Jleak">
        <dc:title>Jleak</dc:title>
        <dcterms:alternative>rate of calcium leakage out of the endoplasmic reticulum</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="micromolar_millisecond" public_interface="out" name="Jleak"/>
    
    <variable units="first_order_rate_constant" name="pleak" initial_value="0.0002"/>
    
    <variable units="micromolar" public_interface="in" name="Ca"/>
    <variable units="micromolar" public_interface="in" name="Caer"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> Jleak </ci> 
        <apply>
          <times/>
          <ci> pleak </ci>
          <apply>
            <minus/>
            <ci> Caer </ci>
            <ci> Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

   <component cmeta:id="I" name="I">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="I">
        <dc:title>I</dc:title>
        <dcterms:alternative>insulin</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="dimensionless" public_interface="out" name="I"/>

    <variable units="dimensionless" name="I_infinity"/>
    <variable units="millisecond" name="tau_I" initial_value="10000.0"/>
    <variable units="per_micromolar" name="I_slope" initial_value="210.0"/>
    <variable units="micromolar" name="Canull" initial_value="0.055"/>
    
    <variable units="micromolar" public_interface="in" name="Ca"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> I </ci> 
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> I_infinity </ci>
            <ci> I </ci>
          </apply>
          <ci> tau_I </ci>
        </apply>
      </apply>

      <apply>
        <eq/>
        <ci> I_infinity </ci>
        <piecewise>
          <piece>
            <apply>
              <times/>
              <ci> I_slope </ci>
              <apply>
                <minus/>
                <ci> Ca </ci>
                <ci> Canull </ci>
              </apply>
            </apply>
            <apply>
              <geq/>
              <ci> Ca </ci>
              <ci> Canull </ci>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="dimensionless"> 0.0 </cn>
            <apply>
              <lt/>
              <ci> Ca </ci>
              <ci> Canull </ci>
            </apply>
          </piece>
        </piecewise>
      </apply>
    </math>
  </component>               

  <component name="model_parameters"> 
    <variable units="dimensionless" public_interface="out" name="kappa" initial_value="0.005"/>
    <variable units="millivolt" public_interface="out" name="vK" initial_value="-75.0"/>
    <variable units="millivolt" public_interface="out" name="vCa" initial_value="25.0"/>
  </component>

  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="I_Ca">
        <component_ref component="m"/>
      </component_ref>
      <component_ref component="I_K">
        <component_ref component="n"/>
      </component_ref>
      <component_ref component="I_K_Ca"/>
      <component_ref component="I_K_ATP"/>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="encapsulation"/>
    <component_ref component="I_Ca">
      <component_ref component="m"/>
    </component_ref>
    <component_ref component="I_K">
      <component_ref component="n"/>
    </component_ref>
  </group>
  
  <connection>
    <map_components component_2="environment" component_1="Gi"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="Ge"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="G6P"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="FBP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="ADP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="membrane"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="I_K"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="I_Ca"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="I_K_Ca"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="I_K_ATP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="Ca"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="Caer"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="I"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="Jgk" component_1="Gi"/>
    <map_variables variable_2="Gi" variable_1="Gi"/>
    <map_variables variable_2="Jgk" variable_1="Jgk"/>
  </connection> 

  <connection>
    <map_components component_2="Jglut" component_1="Gi"/>
    <map_variables variable_2="Gi" variable_1="Gi"/>
    <map_variables variable_2="Jglut" variable_1="Jglut"/>
  </connection>

  <connection>
    <map_components component_2="Jgk" component_1="G6P"/>
    <map_variables variable_2="Jgk" variable_1="Jgk"/>
  </connection> 

  <connection>
    <map_components component_2="JPFK" component_1="G6P"/>
    <map_variables variable_2="JPFK" variable_1="JPFK"/>
  </connection> 

  <connection>
    <map_components component_2="model_parameters" component_1="G6P"/>
    <map_variables variable_2="kappa" variable_1="kappa"/>
  </connection>

  <connection>
    <map_components component_2="F6P" component_1="G6P"/>
    <map_variables variable_2="G6P" variable_1="G6P"/>
  </connection>

  <connection>
    <map_components component_2="JGPDH" component_1="FBP"/>
    <map_variables variable_2="FBP" variable_1="FBP"/>
    <map_variables variable_2="JGPDH" variable_1="JGPDH"/>
  </connection> 

  <connection>
    <map_components component_2="JPFK" component_1="FBP"/>
    <map_variables variable_2="JPFK" variable_1="JPFK"/>
  </connection> 

  <connection>
    <map_components component_2="model_parameters" component_1="FBP"/>
    <map_variables variable_2="kappa" variable_1="kappa"/>
  </connection>

  <connection>
    <map_components component_2="Ge" component_1="Jglut"/>
    <map_variables variable_2="Ge" variable_1="Ge"/>
  </connection>

  <connection>
    <map_components component_2="w" component_1="JPFK"/>
    <map_variables variable_2="w0000" variable_1="w0000"/>
    <map_variables variable_2="w0001" variable_1="w0001"/>
    <map_variables variable_2="w0010" variable_1="w0010"/>
    <map_variables variable_2="w0100" variable_1="w0100"/>
    <map_variables variable_2="w0101" variable_1="w0101"/>
    <map_variables variable_2="w0110" variable_1="w0110"/>
    <map_variables variable_2="w0111" variable_1="w0111"/>
    <map_variables variable_2="w1000" variable_1="w1000"/>
    <map_variables variable_2="w1001" variable_1="w1001"/>
    <map_variables variable_2="w1010" variable_1="w1010"/>
    <map_variables variable_2="w1011" variable_1="w1011"/>
    <map_variables variable_2="w1100" variable_1="w1100"/>
    <map_variables variable_2="w1101" variable_1="w1101"/>
    <map_variables variable_2="w1110" variable_1="w1110"/>
    <map_variables variable_2="w1111" variable_1="w1111"/>
    <map_variables variable_2="w0011" variable_1="w0011"/>
  </connection>

  <connection>
    <map_components component_2="AMP" component_1="w"/>
    <map_variables variable_2="AMP" variable_1="AMP"/>
  </connection>

  <connection>
    <map_components component_2="ATP" component_1="w"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
  </connection>

  <connection>
    <map_components component_2="F6P" component_1="w"/>
    <map_variables variable_2="F6P" variable_1="F6P"/>
  </connection>

  <connection>
    <map_components component_2="FBP" component_1="w"/>
    <map_variables variable_2="FBP" variable_1="FBP"/>
  </connection>

  <connection>
    <map_components component_2="AMP" component_1="ATP"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
	<map_variables variable_2="AMP" variable_1="AMP"/>
  </connection>

  <connection>
    <map_components component_2="ADP" component_1="ATP"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
    <map_variables variable_2="ADP" variable_1="ADP"/>
  </connection>

  <connection>
    <map_components component_2="ADP" component_1="AMP"/>
    <map_variables variable_2="ADP" variable_1="ADP"/>
  </connection>

  <connection>
    <map_components component_2="gamma" component_1="ADP"/>
    <map_variables variable_2="gamma" variable_1="gamma"/>
  </connection>

  <connection>
    <map_components component_2="Ca" component_1="ADP"/>
    <map_variables variable_2="Ca" variable_1="Ca"/>
  </connection>

  <connection>
    <map_components component_2="JGPDH" component_1="gamma"/>
    <map_variables variable_2="JGPDH" variable_1="JGPDH"/>
  </connection>   

  <connection>
    <map_components component_2="I_Ca" component_1="membrane"/>
    <map_variables variable_2="I_Ca" variable_1="I_Ca"/>
    <map_variables variable_2="v" variable_1="v"/>
  </connection>

  <connection>
    <map_components component_2="I_K" component_1="membrane"/>
    <map_variables variable_2="I_K" variable_1="I_K"/>
    <map_variables variable_2="v" variable_1="v"/>
  </connection>

  <connection>
    <map_components component_2="I_K_Ca" component_1="membrane"/>
    <map_variables variable_2="I_K_Ca" variable_1="I_K_Ca"/>
    <map_variables variable_2="v" variable_1="v"/>
  </connection>

  <connection>
    <map_components component_2="I_K_ATP" component_1="membrane"/>
    <map_variables variable_2="I_K_ATP" variable_1="I_K_ATP"/>
    <map_variables variable_2="v" variable_1="v"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="I_K"/>
    <map_variables variable_2="vK" variable_1="vK"/>
  </connection>

  <connection>
    <map_components component_2="I_K" component_1="n"/>
    <map_variables variable_2="n" variable_1="n"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="v" variable_1="v"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="I_Ca"/>
    <map_variables variable_2="vCa" variable_1="vCa"/>
  </connection>

  <connection>
    <map_components component_2="I_Ca" component_1="m"/>
    <map_variables variable_2="m_infinity" variable_1="m_infinity"/>
    <map_variables variable_2="v" variable_1="v"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="I_K_Ca"/>
    <map_variables variable_2="vK" variable_1="vK"/>
  </connection>

  <connection>
    <map_components component_2="Ca" component_1="I_K_Ca"/>
    <map_variables variable_2="Ca" variable_1="Ca"/>
  </connection>  

  <connection>
    <map_components component_2="model_parameters" component_1="I_K_ATP"/>
    <map_variables variable_2="vK" variable_1="vK"/>
  </connection>

  <connection>
    <map_components component_2="ATP" component_1="I_K_ATP"/>
    <map_variables variable_2="ATP" variable_1="ATP"/>
  </connection>

  <connection>
    <map_components component_2="ADP" component_1="I_K_ATP"/>
    <map_variables variable_2="ADP" variable_1="ADP"/>
  </connection>

  <connection>
    <map_components component_2="Jmem" component_1="Ca"/>
    <map_variables variable_2="Ca" variable_1="Ca"/>
    <map_variables variable_2="Jmem" variable_1="Jmem"/>
  </connection>

  <connection>
    <map_components component_2="Jer" component_1="Ca"/>
    <map_variables variable_2="Jer" variable_1="Jer"/>
  </connection>

  <connection>
    <map_components component_2="Jer" component_1="Caer"/>
    <map_variables variable_2="Jer" variable_1="Jer"/>
  </connection>

  <connection>
    <map_components component_2="Jleak" component_1="Jer"/>
    <map_variables variable_2="Jleak" variable_1="Jleak"/>
  </connection>

  <connection>
    <map_components component_2="JSERCA" component_1="Jer"/>
    <map_variables variable_2="JSERCA" variable_1="JSERCA"/>
  </connection>

  <connection>
    <map_components component_2="JSERCA" component_1="Ca"/>
    <map_variables variable_2="Ca" variable_1="Ca"/>
  </connection>

  <connection>
    <map_components component_2="Jleak" component_1="Ca"/>
    <map_variables variable_2="Ca" variable_1="Ca"/>
  </connection>

  <connection>
    <map_components component_2="Jleak" component_1="Caer"/>
    <map_variables variable_2="Caer" variable_1="Caer"/>
  </connection>

  <connection>
    <map_components component_2="I" component_1="Ca"/>
    <map_variables variable_2="Ca" variable_1="Ca"/>
  </connection>

  <connection>
    <map_components component_2="I" component_1="Ge"/>
    <map_variables variable_2="I" variable_1="I"/>
  </connection>

  <connection>
    <map_components component_2="Jmem" component_1="I_Ca"/>
    <map_variables variable_2="I_Ca" variable_1="I_Ca"/>
  </connection>

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#45ceff5c-0914-4676-a569-6b29103af117">
    <rdf:li>endocrine</rdf:li>
    <rdf:li>pancreatic beta cell</rdf:li>
    <rdf:li>insulin</rdf:li>
    <rdf:li>glycolysis</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#67e8a27a-052d-49dc-b5b2-6068a2480b8d">
    <rdf:li rdf:resource="rdf:#bad9a37c-c302-4516-9dc4-e3886a3774b4"/>
    <rdf:li rdf:resource="rdf:#3a313d2c-057a-4029-ad10-72a3d3316205"/>
    <rdf:li rdf:resource="rdf:#561591b1-4918-4773-b7c9-cce0870a442f"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#bad9a37c-c302-4516-9dc4-e3886a3774b4">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#251d50c6-5d8f-4401-b532-c9bc7a52a482"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#f78326c8-0175-4441-8ced-4d8a20ce6154"/>
    <dcterms:created rdf:resource="rdf:#a5cac126-fdd4-4552-b4ca-e7f23bcca187"/>
    <dc:creator rdf:resource="rdf:#a16c727a-ef30-4485-af66-a57dec279703"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3a313d2c-057a-4029-ad10-72a3d3316205">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#40170525-6f91-45fd-9611-ca0454c5c319"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#91155613-aa77-471d-b5af-0cfe5c8f4efc">
    <dc:creator rdf:resource="rdf:#e0a58ba9-933c-4da9-9055-842ae36f646d"/>
    <rdf:value>This is a CellML description of Pedersen, Bertram and Sherman's 2005 mathematical model of metabolically driven insulin secretion in pancreatic beta cells.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e0a58ba9-933c-4da9-9055-842ae36f646d">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#561591b1-4918-4773-b7c9-cce0870a442f">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#22d2c901-28d4-4313-90f0-72a856a59b42"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9e456342-fc16-4a73-adf0-870ed54e0728">
    <dcterms:W3CDTF>2005-07</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#22d2c901-28d4-4313-90f0-72a856a59b42">
    <vCard:Given>Arthur</vCard:Given>
    <vCard:Family>Sherman</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d9a7b236-fbb9-4ebd-8624-54be153d24fa">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a16c727a-ef30-4485-af66-a57dec279703">
    <vCard:ORG rdf:resource="rdf:#5e5d3d82-4a62-4369-9e8a-126b8fa9a5d6"/>
    <vCard:EMAIL rdf:resource="rdf:#d9a7b236-fbb9-4ebd-8624-54be153d24fa"/>
    <vCard:N rdf:resource="rdf:#ca3df698-8191-48cf-b75e-e1ce519f13a5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f78326c8-0175-4441-8ced-4d8a20ce6154">
    <dc:creator rdf:resource="rdf:#1945f0b5-ae27-46c3-a0dc-19b5b6306be9"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e436f2b2-66b8-4d7b-bb5d-8bbfe375e658">
    <bqs:Pubmed_id>15834002</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#e4332344-ead5-4639-b20b-50e02d355aa7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c07164e3-16a2-4035-874d-ee74fa37738b">
    <dc:title>Biophysical Journal</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#251d50c6-5d8f-4401-b532-c9bc7a52a482">
    <vCard:Given>Morten</vCard:Given>
    <vCard:Family>Pedersen</vCard:Family>
    <vCard:Other>Gram</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1945f0b5-ae27-46c3-a0dc-19b5b6306be9">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ca3df698-8191-48cf-b75e-e1ce519f13a5">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5e5d3d82-4a62-4369-9e8a-126b8fa9a5d6">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a5cac126-fdd4-4552-b4ca-e7f23bcca187">
    <dcterms:W3CDTF>2007-07-27T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#40170525-6f91-45fd-9611-ca0454c5c319">
    <vCard:Given>Richard</vCard:Given>
    <vCard:Family>Bertram</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="#pedersen_model_2005">
    <dc:title>
        Pedersen, Bertram and Sherman's 2005 mathematical model of metabolically driven insulin secretion in pancreatic beta cells.
      </dc:title>
    <cmeta:comment rdf:resource="rdf:#91155613-aa77-471d-b5af-0cfe5c8f4efc"/>
    <bqs:reference rdf:resource="rdf:#e436f2b2-66b8-4d7b-bb5d-8bbfe375e658"/>
    <bqs:reference rdf:resource="rdf:#e6907c35-8f0a-42d4-9e4f-37e0ed8f3543"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e6907c35-8f0a-42d4-9e4f-37e0ed8f3543">
    <dc:subject rdf:resource="rdf:#6b5f4407-96eb-4dd6-9e00-f47fb39a243d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e4332344-ead5-4639-b20b-50e02d355aa7">
    <dc:creator rdf:resource="rdf:#67e8a27a-052d-49dc-b5b2-6068a2480b8d"/>
    <dc:title>Intra- and inter-islet synchronization of metabolically driven insulin secretion</dc:title>
    <bqs:volume>89</bqs:volume>
    <bqs:first_page>107</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#c07164e3-16a2-4035-874d-ee74fa37738b"/>
    <dcterms:issued rdf:resource="rdf:#9e456342-fc16-4a73-adf0-870ed54e0728"/>
    <bqs:last_page>119</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6b5f4407-96eb-4dd6-9e00-f47fb39a243d">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#45ceff5c-0914-4676-a569-6b29103af117"/>
  </rdf:Description>
</rdf:RDF>
</model>