Location: Mosekilde, Lading, Yanchuk, Maistrenko, 2001 @ aa0e676c403f / mosekilde_lading_yanchuk_maistrenko_2001.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2009-06-17 15:12:26+12:00
Desc:
committing version01 of mosekilde_lading_yanchuk_maistrenko_2001
Permanent Source URI:
http://models.cellml.org/workspace/mosekilde_lading_yanchuk_maistrenko_2001/rawfile/aa0e676c403f6edc747b9a481fcfbb039baab391/mosekilde_lading_yanchuk_maistrenko_2001.cellml

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

CREATED :  14th January 2003

LAST MODIFIED : 9th April 2003

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/01/2002 CellML Metadata 1.0 Specification.

DESCRIPTION :  This file contains a CellML description of Mosekilde et al's
2001 mathematical model of bursting pancreatic beta cells.

CHANGES:  
  09/04/2003 - AAC - Added publication date information.  
    
--><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="mosekilde_lading_yanchuk_maistrenko_2001_version01" name="mosekilde_lading_yanchuk_maistrenko_2001_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Bifurcation Structure of a Model of Bursting Pancreatic Cells</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This model contains partial differentials and as such can not currently be solved by existing CellML tools.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
Through the continuous activity of ion pumps in the plasma membrane, relatively steep transmembrane electrochemical gradients are maintained.  Because the transmembrane distributions of many ions such as Na<superscript>+</superscript>, K<superscript>+</superscript> and Ca<superscript>2+</superscript> are far from equilibrium, cells are inclined towards excitable electrical activity, or the membrane potential, determined by relative charge distribution, exhibits complicated patterns of slow and fast oscillations associated with variations in the transmembrane ionic currents.  These dynamics play an essential role in intracellular function and intercellular communication.  
</para>

<para>
When exposed to a threshold concentration of glucose, pancreatic beta-cells exhibit a complicated pattern of electrical activity.  Bursts of action potential spikes (the "active" phase) are observed, separated by a "silent" phase of membrane repolarisation.  Insulin secretion is dependent on the relative amount of time the cells spend in the active state, and this time increases with the concentration of extracellular glucose.  Bursting dynamics control the influx of Ca<superscript>2+</superscript> into the cell, and in turn, Ca<superscript>2+</superscript> triggers the release of insulin.  So, via this mechanism, the bursting dynamics of the pancreatic beta-cells regulate their response to varying glucose concentrations.  Below a threshold glucose concentration (about 5mM), the cells do not burst.  Above a certain concentration (&gt; 22 mM), the cells display continuous spiking, and insulin secretion becomes saturated. 
</para>

<para>
Mathematical models of pancreatic beta-cells are usually based on the Hodgkin-Huxley formalism (for the original model see <ulink url="${HTML_EXMPL_HHSA_INTRO}">The Hodgkin-Huxley Squid Axon Model, 1952</ulink>), with elaborations to account for Ca<superscript>2+</superscript> storage, glucose metabolism, and regulation by ATP and various hormones.  Several models have been published (see <ulink url="${HTML_REPOSITORY_INTRODUCTION}">the CellML repository</ulink> for more pancreatic beta-cell and calcium oscillation models), but although some have been established for some years, their bifurcation structure is so complicated it is not yet fully understood.  The goal of the 2001 Mosekilde <emphasis>et al.</emphasis> paper analysed here was to give an accurate, but simplified account of the bifurcation structure of the individual beta-cell (see <xref linkend="fig_cell_diagram"/> below).
</para>

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

<para>
<ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&amp;_udi=B6T2K-4447MJD-J&amp;_user=140507&amp;_coverDate=12%2F31%2F2001&amp;_alid=70440630&amp;_rdoc=1&amp;_fmt=summary&amp;_orig=search&amp;_qd=1&amp;_cdi=4921&amp;_sort=d&amp;_docanchor=&amp;wchp=dGLbVlz-lSzBk&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=dc19bf053e438dae9918fa3c279f2f7e">Bifurcation structure of a model of bursting pancreatic cells</ulink>, Erik Mosekilde, Brian Lading, Sergiy Yanchuk and Yuri Maistrenko, 2001, <ulink url="http://www.sciencedirect.com/science?_ob=JournalURL&amp;_cdi=4921&amp;_auth=y&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=511bbadc6ebfd16e10115a252026e228">
            <emphasis>Biosystems</emphasis>
          </ulink>, 63, 3-13.  (PDF and <ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&amp;_udi=B6T2K-4447MJD-J&amp;_coverDate=12%2F31%2F2001&amp;_alid=70440630&amp;_rdoc=1&amp;_fmt=&amp;_orig=search&amp;_qd=1&amp;_cdi=4921&amp;_sort=d&amp;wchp=dGLbVlz-lSzBk&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=63f03834b470d0be17a65ca137c26f84">full text</ulink> versions of the article are available to subscribers on the journal website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=11595325&amp;dopt=Abstract">PubMed ID: 11595325</ulink>
</para>

<para>
The raw CellML description of the model can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>diagram of the model</title>
    </objectinfo>
    <imagedata fileref="mosekilde_2001.png"/>
  </imageobject>
</mediaobject>
<caption>A schematic representation of the ionic currents described by the Mosekilde <emphasis>et al.</emphasis> 2001 pancreatic beta-cell model.  The model includes a calcium current, I <subscript>Ca<superscript>2+</superscript>
            </subscript>; a potassium current, I<subscript>K</subscript>; and a slow current, I<subscript>S</subscript>.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
  
  
  <!--
    Below, we define some additional units for association with variables and
    constants within the model. 
  -->
  
  <units name="per_second">
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>
  
  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="microS_per_cm2">
    <unit units="siemens" prefix="micro"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </units>
  
  <units name="nanoA_per_cm2">
    <unit units="ampere" prefix="nano"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </units>
  
  <component name="environment">
    <variable units="second" public_interface="out" name="time"/>
  </component>
  
  <component name="membrane">
    <variable units="millivolt" public_interface="out" name="V"/>
    <variable units="second" public_interface="out" name="tau" initial_value="0.02"/>
             
    <variable units="second" public_interface="in" name="time"/>
    <variable units="nanoA_per_cm2" public_interface="in" name="i_K"/>
    <variable units="nanoA_per_cm2" public_interface="in" name="i_Ca"/>
    <variable units="nanoA_per_cm2" public_interface="in" name="i_S"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="membrane_voltage_diff_eq">
        <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_S </ci>
            </apply>
          </apply>
          <ci> tau </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_current">
    <variable units="nanoA_per_cm2" public_interface="out" name="i_Ca"/>
    <variable units="dimensionless" public_interface="out" name="J11"/>
   
    <variable units="millivolt" name="V_Ca" initial_value="25.0"/>
    <variable units="microS_per_cm2" name="g_Ca" initial_value="3.6"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="millivolt" public_interface="in" name="V_K"/>
    <variable units="dimensionless" private_interface="in" name="m_infinity"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Ca_calculation">
        <eq/>
        <ci> i_Ca </ci>
        <apply>
          <times/>
          <ci> g_Ca </ci>
          <ci> m_infinity </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_Ca </ci>
          </apply>
        </apply>
      </apply>

      <apply id="J11_calculation">
        <eq/>
        <ci> J11 </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <ci> g_Ca </ci>
            </apply>
            <apply>
              <partialdiff/>
              <bvar>
                <ci> V </ci>
              </bvar>
              <ci> m_infinity </ci>
            </apply>
            <apply>
              <minus/>
              <ci> V </ci>
              <ci> V_Ca </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> g_Ca </ci>
            <ci> m_infinity </ci>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci> V_Ca </ci>
                <ci> V_K </ci>
              </apply>
              <apply>
                <minus/>
                <ci> V </ci>
                <ci> V_K </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_current_m_gate">
    <variable units="dimensionless" public_interface="out" name="m_infinity"/>
    
    <variable units="millivolt" name="theta_m" initial_value="12.0"/>
    <variable units="millivolt" name="V_m" initial_value="-20.0"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="m_infinity_calculation">
        <eq/>
        <ci> m_infinity </ci>
        <apply>
          <power/>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> V_m </ci>
                  <ci> V </ci>
                </apply>
                <ci> theta_m </ci>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless"> -1.0 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="potassium_current">
    <variable units="nanoA_per_cm2" public_interface="out" name="i_K"/>
    <variable units="millivolt" public_interface="out" name="V_K" initial_value="-75.0"/>
    <variable units="dimensionless" public_interface="out" name="J12"/>
    <variable units="dimensionless" public_interface="out" private_interface="in" name="J21"/>
    <variable units="dimensionless" public_interface="out" private_interface="in" name="J22"/>
    
    <variable units="microS_per_cm2" name="g_K" initial_value="10.0"/>
    
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="second" public_interface="in" private_interface="out" name="tau"/>
    
    <variable units="dimensionless" private_interface="in" name="n"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_K_calculation">
        <eq/>
        <ci> i_K </ci>
        <apply>
          <times/>
          <ci> g_K </ci>
          <ci> n </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_K </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="J12_calculation">
        <eq/>
        <ci> J12 </ci>
        <apply>
          <times/>
          <apply>
            <minus/>
            <ci> g_K </ci>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="potassium_current_n_gate">
    <variable units="dimensionless" public_interface="out" name="n"/>
    <variable units="dimensionless" public_interface="out" name="J21"/>
    <variable units="dimensionless" public_interface="out" name="J22"/>
    
    <variable units="dimensionless" name="n_infinity"/>
    <variable units="millivolt" name="theta_n" initial_value="5.6"/>
    <variable units="dimensionless" name="delta" initial_value="0.85"/>
    <variable units="millivolt" name="V_n" initial_value="-16.0"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="tau"/>
    <variable units="second" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="n_infinity_calculation">
        <eq/>
        <ci> n_infinity </ci>
        <apply>
          <power/>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> V_n </ci>
                  <ci> V </ci>
                </apply>
                <ci> theta_n </ci>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless"> -1.0 </cn>
        </apply>
      </apply>
      
      <apply id="n_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> n </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> delta </ci>
            <apply>
              <minus/>
              <ci> n_infinity </ci>
              <ci> n </ci>
            </apply>
          </apply>
          <ci> tau </ci>
        </apply>
      </apply>
      
      <apply id="J21_calculation">
        <eq/>
        <ci> J21 </ci>
        <apply>
          <times/>
          <ci> delta </ci>
          <apply>
            <partialdiff/>
            <bvar>
              <ci> V </ci>
            </bvar>
            <ci> n_infinity </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="J22_calculation">
        <eq/>
        <ci> J22 </ci>
        <apply>
          <minus/>
          <ci> delta </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_current">
    <variable units="nanoA_per_cm2" public_interface="out" name="i_S"/> 
    <variable units="second" private_interface="out" name="tau_S" initial_value="35.0"/>
    <variable units="dimensionless" public_interface="out" name="J13"/>
    <variable units="dimensionless" public_interface="out" name="J33"/>
    <variable units="dimensionless" public_interface="out" private_interface="in" name="J31"/>
    
    <variable units="microS_per_cm2" name="g_S" initial_value="4.0"/>
    <variable units="dimensionless" name="k_S"/>
    
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="millivolt" public_interface="in" name="V_K"/>
    <variable units="dimensionless" private_interface="in" name="S"/>
    <variable units="second" public_interface="in" name="tau"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_S_calculation">
        <eq/>
        <ci> i_S </ci>
        <apply>
          <times/>
          <ci> g_S </ci>
          <ci> S </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_K </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="J13_calculation">
        <eq/>
        <ci> J13 </ci>
        <apply>
          <times/>
          <apply>
            <minus/>
            <ci> g_S </ci>
          </apply>
          <ci> k_S </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_K </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="k_S_calculation">
        <eq/>
        <ci> k_S </ci>
        <apply>
          <divide/>
          <ci> tau </ci>
          <ci> tau_S </ci>
        </apply>
      </apply>
      
      <apply id="J33_calculation">
        <eq/>
        <ci> J33 </ci>
        <apply>
          <minus/>
          <ci> k_S </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="slow_current_S_gate">
    <variable units="dimensionless" public_interface="out" name="S"/>
    <variable units="dimensionless" public_interface="out" name="J31"/>
    
    <variable units="dimensionless" name="S_infinity"/>
    <variable units="millivolt" name="theta_S" initial_value="10.0"/>
    <variable units="millivolt" name="V_S" initial_value="-38.34"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="second" public_interface="in" name="tau_S"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="S_infinity_calculation">
        <eq/>
        <ci> S_infinity </ci>
        <apply>
          <power/>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> V_S </ci>
                  <ci> V </ci>
                </apply>
                <ci> theta_S </ci>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless"> -1.0 </cn>
        </apply>
      </apply>
      
      <apply id="S_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> S </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> S_infinity </ci>
            <ci> S </ci>
          </apply>
          <ci> tau_S </ci>
        </apply>
      </apply>
      
      <apply id="J31_calculation">
        <eq/>
        <ci> J31 </ci>
        <apply>
          <partialdiff/>
          <bvar>
            <ci> V </ci>
          </bvar>
          <ci> S_infinity </ci>
        </apply>
      </apply>
    </math>
  </component> 
  
  <component name="jacobian_matrix">
    <variable units="dimensionless" name="J"/>
    
    <variable units="dimensionless" public_interface="in" name="J11"/>
    <variable units="dimensionless" public_interface="in" name="J12"/>
    <variable units="dimensionless" public_interface="in" name="J13"/>
    <variable units="dimensionless" public_interface="in" name="J21"/>
    <variable units="dimensionless" public_interface="in" name="J22"/>
    <variable units="dimensionless" public_interface="in" name="J31"/>
    <variable units="dimensionless" public_interface="in" name="J33"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="J_calculation">
        <eq/>
        <ci> J </ci>
        <matrix>
          <matrixrow>
            <ci> J11 </ci>
            <ci> J12 </ci>
            <ci> J13 </ci>
          </matrixrow>
          <matrixrow>
            <ci> J21 </ci>
            <ci> J22 </ci>
            <cn cellml:units="dimensionless"> 0 </cn>
          </matrixrow>
          <matrixrow>
            <ci> J31 </ci>
            <cn cellml:units="dimensionless"> 0 </cn>
            <ci> J33 </ci> 
          </matrixrow>
        </matrix>
      </apply>
    </math>
  </component>
  
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="calcium_current">
        <component_ref component="calcium_current_m_gate"/>
      </component_ref>
      <component_ref component="potassium_current">
        <component_ref component="potassium_current_n_gate"/>
      </component_ref>
      <component_ref component="slow_current">
        <component_ref component="slow_current_S_gate"/>
      </component_ref>
      <component_ref component="jacobian_matrix"/>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="encapsulation"/> 
    <component_ref component="calcium_current">
      <component_ref component="calcium_current_m_gate"/>
    </component_ref>
    <component_ref component="potassium_current">
      <component_ref component="potassium_current_n_gate"/>
    </component_ref>
    <component_ref component="slow_current">
      <component_ref component="slow_current_S_gate"/>
    </component_ref>
  </group>
      
  <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="calcium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="slow_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="membrane" component_1="calcium_current"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Ca" variable_1="i_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="membrane" component_1="potassium_current"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_K" variable_1="i_K"/>
    <map_variables variable_2="tau" variable_1="tau"/>
  </connection>  
  
  <connection>
    <map_components component_2="membrane" component_1="slow_current"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_S" variable_1="i_S"/>
    <map_variables variable_2="tau" variable_1="tau"/>
  </connection>
  
  <connection>
    <map_components component_2="jacobian_matrix" component_1="calcium_current"/>
    <map_variables variable_2="J11" variable_1="J11"/>
  </connection>  
  
  <connection>
    <map_components component_2="slow_current" component_1="potassium_current"/>
    <map_variables variable_2="V_K" variable_1="V_K"/>
  </connection>  
  
  <connection>
    <map_components component_2="calcium_current" component_1="potassium_current"/>
    <map_variables variable_2="V_K" variable_1="V_K"/>
  </connection>  
  
  <connection>
    <map_components component_2="jacobian_matrix" component_1="potassium_current"/>
    <map_variables variable_2="J12" variable_1="J12"/>
    <map_variables variable_2="J21" variable_1="J21"/>
    <map_variables variable_2="J22" variable_1="J22"/>
  </connection>
  
  <connection>
    <map_components component_2="jacobian_matrix" component_1="slow_current"/>
    <map_variables variable_2="J13" variable_1="J13"/>
    <map_variables variable_2="J31" variable_1="J31"/>
    <map_variables variable_2="J33" variable_1="J33"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_current_m_gate" component_1="calcium_current"/>
    <map_variables variable_2="m_infinity" variable_1="m_infinity"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection> 
  
  <connection>
    <map_components component_2="potassium_current_n_gate" component_1="potassium_current"/>
    <map_variables variable_2="n" variable_1="n"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="tau" variable_1="tau"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="J21" variable_1="J21"/>
    <map_variables variable_2="J22" variable_1="J22"/>
  </connection>
  
  <connection>
    <map_components component_2="slow_current_S_gate" component_1="slow_current"/>
    <map_variables variable_2="S" variable_1="S"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="tau_S" variable_1="tau_S"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="J31" variable_1="J31"/>
  </connection>   


<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#fd7b97ae-ffdf-48dd-b306-c473c86858c9">
    <rdf:li>calcium dynamics</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>beta cell</rdf:li>
    <rdf:li>Pancreatic Beta-Cell</rdf:li>
    <rdf:li>pancreas</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#c3faec9f-3018-4305-950e-be3dcb8edc12">
    <rdf:li rdf:resource="rdf:#c24ff00e-585d-4acc-b5b1-436b23803f23"/>
    <rdf:li rdf:resource="rdf:#15d524f8-6e4a-4540-ace6-c52c5b0e1663"/>
    <rdf:li rdf:resource="rdf:#df2d4b7a-54a8-42e2-b6f9-aed0d292b225"/>
    <rdf:li rdf:resource="rdf:#db1de282-66d7-4bdc-ad64-43733d419371"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#db1de282-66d7-4bdc-ad64-43733d419371">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#cb11a1e5-702a-43be-8e55-71b4b43c8690"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d4a04805-9fcb-4333-a170-c8754d4c24e6">
    <dc:title>Biosystems</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#15d524f8-6e4a-4540-ace6-c52c5b0e1663">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#b21408b4-0335-4ec4-8336-e595898f69aa"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#347a770d-4cd6-46cc-8f2e-d18d94ad32c8">
    <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:#83e73c1b-7e75-46a6-a90c-569963c22562">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4513d6fd-0d79-4af8-83e0-259fb6bce781">
    <dc:subject rdf:resource="rdf:#fa520c71-cbfb-4a93-911a-7efb99c7b230"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Bioengineering Institute
      </dc:publisher>
    <cmeta:modification rdf:resource="rdf:#491370b3-8caf-4258-b087-f9b0bde9fa28"/>
    <dcterms:created rdf:resource="rdf:#a3894cc9-a303-4ea8-94e7-e5aeddf79d0d"/>
    <dc:creator rdf:resource="rdf:#6f732587-9bb4-4951-851a-c5b5421ae3ae"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#99865786-f410-416f-a7be-92f11a9b911e">
    <dcterms:W3CDTF>2001-11</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#622b2450-e962-4071-8a87-af8dd1301253">
    <dc:creator rdf:resource="rdf:#c3faec9f-3018-4305-950e-be3dcb8edc12"/>
    <dc:title>
            Bifurcation structure of a model of bursting pancreatic cells
          </dc:title>
    <bqs:volume>63</bqs:volume>
    <bqs:first_page>3</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#d4a04805-9fcb-4333-a170-c8754d4c24e6"/>
    <dcterms:issued rdf:resource="rdf:#99865786-f410-416f-a7be-92f11a9b911e"/>
    <bqs:last_page>13</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b21408b4-0335-4ec4-8336-e595898f69aa">
    <vCard:Given>Brian</vCard:Given>
    <vCard:Family>Lading</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1751e80d-f892-40a3-b14c-000ed56eb1f0">
    <vCard:Given>Erik</vCard:Given>
    <vCard:Family>Mosekilde</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#491370b3-8caf-4258-b087-f9b0bde9fa28">
    <dcterms:modified rdf:resource="rdf:#6ad764dd-9c13-40c8-bf22-1dfe486db4ca"/>
    <rdf:value>
          Added publication date information.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#a6d19d6a-a594-482e-a2c8-62e9bd788cb6"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a6d19d6a-a594-482e-a2c8-62e9bd788cb6">
    <vCard:N rdf:resource="rdf:#da1567b6-fe87-4497-84af-10bbcf40c549"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1d41ac58-bff1-4483-a173-7186dd9c9de9">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6ad764dd-9c13-40c8-bf22-1dfe486db4ca">
    <dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#df2d4b7a-54a8-42e2-b6f9-aed0d292b225">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#35e92aa4-660a-4d13-a1b0-19f3259f7904"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6f732587-9bb4-4951-851a-c5b5421ae3ae">
    <vCard:ORG rdf:resource="rdf:#22d928f3-6249-4205-bc3a-2e38ebae5eed"/>
    <vCard:EMAIL rdf:resource="rdf:#347a770d-4cd6-46cc-8f2e-d18d94ad32c8"/>
    <vCard:N rdf:resource="rdf:#83e73c1b-7e75-46a6-a90c-569963c22562"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a3894cc9-a303-4ea8-94e7-e5aeddf79d0d">
    <dcterms:W3CDTF>2003-01-14</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#00d39e01-5f94-4446-b6df-cd340d732195">
    <dc:creator rdf:resource="rdf:#1d41ac58-bff1-4483-a173-7186dd9c9de9"/>
    <rdf:value>
          This is the CellML description of Mosekilde et al's 2001 mathematical 
          model of bursting pancreatic beta cells.
        </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#da1567b6-fe87-4497-84af-10bbcf40c549">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b0ed9a82-b241-4eb0-a7af-25c7ca5e54fc">
    <bqs:Pubmed_id>11595325</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#622b2450-e962-4071-8a87-af8dd1301253"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#35e92aa4-660a-4d13-a1b0-19f3259f7904">
    <vCard:Given>Sergiy</vCard:Given>
    <vCard:Family>Yanchuk</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#22d928f3-6249-4205-bc3a-2e38ebae5eed">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c24ff00e-585d-4acc-b5b1-436b23803f23">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#1751e80d-f892-40a3-b14c-000ed56eb1f0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fa520c71-cbfb-4a93-911a-7efb99c7b230">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#fd7b97ae-ffdf-48dd-b306-c473c86858c9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cb11a1e5-702a-43be-8e55-71b4b43c8690">
    <vCard:Given>Yuri</vCard:Given>
    <vCard:Family>Maistrenko</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="#mosekilde_lading_yanchuk_maistrenko_2001_version01">
    <dc:title>
        Mosekilde et al's 2001 mathematical model of bursting pancreatic beta 
        cells.
      </dc:title>
    <cmeta:bio_entity>Pancreatic Beta-cell</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#00d39e01-5f94-4446-b6df-cd340d732195"/>
    <bqs:reference rdf:resource="rdf:#4513d6fd-0d79-4af8-83e0-259fb6bce781"/>
    <bqs:reference rdf:resource="rdf:#b0ed9a82-b241-4eb0-a7af-25c7ca5e54fc"/>
  </rdf:Description>
</rdf:RDF>
</model>