Location: Warashina, Ogura, 2004 @ 0084a08956b5 / warashina_ogura_2004.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2009-06-17 16:21:29+12:00
Desc:
committing version01 of warashina_ogura_2004
Permanent Source URI:
https://models.cellml.org/workspace/warashina_ogura_2004/rawfile/0084a08956b59d8e17535af20ef5aad3bfff4b67/warashina_ogura_2004.cellml

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

CREATED : 3rd February 2004 

LAST MODIFIED : 27th April 2005

AUTHOR : Catherine Lloyd
         Bioengineering Instute
         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 Warashina and Ogura's 2004 mathematical model of stimulation-secretion coupling in a chromaffin cell.

CHANGES: 
  27/04/2005 - PJV - Made MathML id's unique. Changed unit dimension
                     to correct inconsistancy. Removed undefined
                     initial values.
--><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="warashina_ogura_2004_version01" name="warashina_ogura_2004_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modelling Stimulation-Sectretion Coupling in a Chromaffin Cell</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 is the original unchecked version of the model imported from the previous
            CellML model repository, 24-Jan-2006.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
Adrenal chromaffin cells have been used extensively for studying the mechanisms involved in stimulation-secretion coupling, namely electrophysiological progesses, signal transduction pathways, calcium homeostasis, and exocytotic processes.  However, despite the accululation of large quantities of experimental data, a clear understanding of how these processes interact is lacking.  In the Warashina and Ogura 2004 publication described here, the authors present a mathematical model which describes a chromaffin cell (see the figure below).  This model is able to simulate the cell's secretory responses, including the process of stimulation-secretion coupling.
</para>

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

<para>
<ulink url="http://springerlink.metapress.com/app/home/contribution.asp?wasp=acecbdaqwg6yph41qjft&amp;referrer=parent&amp;backto=issue,5,16;journal,5,110;linkingpublicationresults,1:100448,1">Modeling of stimulation-secretion coupling in a chromaffin cell</ulink>, A. Warashina and T. Ogura, 2004, <ulink url="http://www.springerlink.com/app/home/journal.asp?wasp=g2e72bwqwj3txlfgxldv&amp;referrer=backto&amp;backto=searchpublicationsresults,3,3;&amp;absoluteposition=1#A1">
            <emphasis>European Journal of Physiology</emphasis>
          </ulink>, 448, 369-391.  (<ulink url="http://springerlink.metapress.com/media/988TNJWHWL5HJFH4UC5V/Contributions/F/B/A/4/FBA4F5TX6P4PEC61_html/fulltext.html">A Full text (HTML)</ulink> and PDF versions of the article are available to subscribers on the <emphasis>European Journal of Physiology</emphasis> website.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=14730449&amp;dopt=Abstract">PubMed ID: 14730449</ulink> 
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram</title>
    </objectinfo>
    <imagedata fileref="warashina_2004.png"/>
  </imageobject>
</mediaobject>
<caption>The model chromaffin cell includes mechanisms involved in the excitatory synapse, Hodgkin-Huxley type Na<superscript>+</superscript> and K<superscript>+</superscript> channels, volatage-dependent Ca<superscript>2+</superscript> channels, Ca<superscript>2+</superscript>-dependent K<superscript>+</superscript> channels, Ca<superscript>2+</superscript> extrusion, buffered diffusion of intracellular Ca<superscript>2+</superscript>, Ca<superscript>2+</superscript> indication with a fluorescent dye, and Ca<superscript>2+</superscript>-triggered secretion.</caption>
</informalfigure>

<para>
Model simulations were carried out using the <ulink url="http://www.neuron.yale.edu/">NEURON</ulink> simulation environment.  The simulations were able to reproduce several of the experimental results reported in the literature.
</para>

</sect1>
</article>
</documentation>
  
  
  
  <units name="per_millisecond">
    <unit units="second" prefix="milli" exponent="-1"/>
  </units>
  
  <units name="per_millimolar_per_millisecond">
    <unit units="millimolar" exponent="-1"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>
  
  <units name="per_millivolt_millisecond">
    <unit units="millivolt" exponent="-1"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="millisecond">
    <unit units="second" prefix="milli"/>
  </units>
  
  <units name="micrometre_per_second">
    <unit units="micrometre"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="metre_per_second">
    <unit units="metre"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="micrometre2_per_millisecond">
    <unit units="micrometre" exponent="2"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="micromolar_per_millisecond">
    <unit units="micromolar"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="millimolar_per_millisecond">
    <unit units="millimolar"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="per_micromolar3_per_millisecond">
    <unit units="micromolar" exponent="-3"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="microS">
    <unit units="siemens" prefix="micro"/>
  </units>
 
  <units name="milliS_per_cm2">
    <unit units="siemens" prefix="milli"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </units>
  
  <units name="microF_per_cm2">
    <unit units="farad" prefix="micro"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </units>
  
  <units name="microA_per_cm2">
    <unit units="ampere" prefix="micro"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </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="nanomolar">
    <unit units="mole" prefix="nano"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="micrometre">
    <unit units="metre" prefix="micro"/>
  </units>
  
  <units name="micrometre2">
    <unit units="metre" prefix="micro" exponent="2"/>
  </units>
  
  <units name="micrometre3">
    <unit units="metre" prefix="micro" exponent="3"/>
  </units>
  
  <units name="joule_per_mole_kelvin">
    <unit units="joule"/>
    <unit units="mole" exponent="-1"/>
    <unit units="kelvin" exponent="-1"/>
  </units>
  
  <units name="coulomb_per_mole">
    <unit units="coulomb"/>
    <unit units="mole" exponent="-1"/>
  </units>
    
  
  <component name="environment">
    <variable units="millisecond" public_interface="out" name="time"/>
  </component>
  
  <component name="membrane">
    <variable units="millivolt" public_interface="out" name="V"/>
              
    <variable units="microF_per_cm2" name="Cm" initial_value="1.0"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_Na"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_sp"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_K"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_l"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_Ca"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_SK"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_ex"/>
     
    <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_sp  </ci>
              <ci> i_Na </ci>
              <ci> i_K </ci>
              <ci> i_l </ci>
              <ci> i_Ca </ci>
              <ci> i_SK </ci>
              <ci> i_ex </ci>
            </apply>
          </apply>
          <ci> Cm </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="excitatory_synapse">
    <variable units="microA_per_cm2" public_interface="out" name="i_sp"/>

    <variable units="milliS_per_cm2" name="G"/>
    <variable units="microS" name="L" initial_value="0.016"/>
    <variable units="microS" name="D"/>
    <variable units="millivolt" name="E_rev" initial_value="-7.0"/>
    <variable units="millisecond" name="tau1" initial_value="1.0"/>
    <variable units="millisecond" name="tau2" initial_value="6.3"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_sp_calculation">
        <eq/>
        <ci> i_sp </ci>
        <apply>
          <times/>
          <ci> G </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_rev </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="G_calculation">
        <eq/>
        <ci> G </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> L </ci>
            <ci> tau1 </ci>
          </apply>
          <apply>
            <times/>
            <ci> D </ci>
            <ci> tau2 </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sodium_channel">
    <variable units="microA_per_cm2" public_interface="out" name="i_Na"/>

    <variable units="milliS_per_cm2" name="g_Na" initial_value="40.0"/>
    <variable units="millivolt" name="E_Na" initial_value="50.0"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    
    <variable units="dimensionless" private_interface="in" name="m"/>
    <variable units="dimensionless" private_interface="in" name="h"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Na_calculation">
        <eq/>
        <ci> i_Na </ci>
        <apply>
          <times/>
          <ci> g_Na </ci>
          <apply>
            <power/>
            <ci> m </ci>
            <cn cellml:units="dimensionless"> 3.0 </cn>
          </apply>
          <ci> h </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Na </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sodium_channel_m_gate">
    <variable units="dimensionless" public_interface="out" name="m"/>
    
    <variable units="per_millisecond" name="alpha_m"/>
    <variable units="per_millisecond" name="beta_m"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="sodium_channel_m_gate_alpha_m_calculation">
        <eq/>
        <ci> alpha_m </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millisecond"> 0.109 </cn>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 38.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci> V </ci> 
                    <cn cellml:units="millivolt"> 38.0 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 5.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_channel_m_gate_beta_m_calculation">
        <eq/>
        <ci> beta_m </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millisecond"> 0.0744 </cn>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 38.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> V </ci> 
                  <cn cellml:units="millivolt"> 38.0 </cn>
                </apply>
                <cn cellml:units="millivolt"> 5.0 </cn>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="sodium_channel_m_gate_dm_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> m </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_m </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> m </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_m </ci>
            <ci> m </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="sodium_channel_h_gate">
    <variable units="dimensionless" public_interface="out" name="h"/>
    
    <variable units="per_millisecond" name="alpha_h"/>
    <variable units="per_millisecond" name="beta_h"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_h_calculation">
        <eq/>
        <ci> alpha_h </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millisecond"> 0.0192 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <apply>
                  <plus/>
                  <cn cellml:units="millivolt"> 55.0 </cn>
                  <ci> V </ci> 
                </apply>
              </apply>
              <cn cellml:units="millivolt"> 15.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_h_calculation">
        <eq/>
        <ci> beta_h </ci>
        <apply>
          <divide/>
          <cn cellml:units="per_millisecond"> 2.48 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt"> 17.0 </cn>
                  <ci> V </ci> 
                </apply>
                <cn cellml:units="millivolt"> 21.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="dh_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> h </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_h </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> h </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_h </ci>
            <ci> h </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="potassium_channel">
    <variable units="microA_per_cm2" public_interface="out" name="i_K"/>

    <variable units="milliS_per_cm2" name="g_K" initial_value="1.2"/>
    
    <variable units="millisecond" 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="E_K"/>
    
    <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>
          <apply>
            <power/>
            <ci> n </ci>
            <cn cellml:units="dimensionless"> 4.0 </cn>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="potassium_channel_n_gate">
    <variable units="dimensionless" public_interface="out" name="n"/>
    
    <variable units="per_millisecond" name="alpha_n"/>
    <variable units="per_millisecond" name="beta_n"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_n_calculation">
        <eq/>
        <ci> alpha_n </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millisecond"> 0.012 </cn>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 45.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <cn cellml:units="millivolt"> 45.0 </cn>
                    <ci> V </ci> 
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 5.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_n_calculation">
        <eq/>
        <ci> beta_n </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millisecond"> 0.204 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <apply>
                  <plus/>
                  <ci> V </ci> 
                  <cn cellml:units="millivolt"> 50.0 </cn>
                </apply>
              </apply>
              <cn cellml:units="millivolt"> 40.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="dn_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> n </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_n </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> n </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_n </ci>
            <ci> n </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="leak_channel">
    <variable units="microA_per_cm2" public_interface="out" name="i_l"/>

    <variable units="milliS_per_cm2" name="g_l" initial_value="1.6"/>
    <variable units="millivolt" name="E_l" initial_value="-65.0"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_l_calculation">
        <eq/>
        <ci> i_l </ci>
        <apply>
          <times/>
          <ci> g_l </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_l </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_channel">
    <variable units="microA_per_cm2" public_interface="out" name="i_Ca"/>

    <variable units="micrometre_per_second" name="PCa" initial_value="2.0"/>
    <variable units="dimensionless" name="ghk"/>
    <variable units="millimolar" name="Cao" initial_value="2.0"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="millimolar" public_interface="in" name="Ca0"/>
    <variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    
    <variable units="dimensionless" private_interface="in" name="m"/>
    <variable units="dimensionless" private_interface="in" name="h"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Ca_calculation">
        <eq/>
        <ci> i_Ca </ci>
        <apply>
          <times/>
          <ci> PCa </ci>
          <apply>
            <power/>
            <ci> m </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <ci> h </ci>
          <ci> ghk </ci>
        </apply>
      </apply>
      
      <apply id="ghk_calculation">
        <eq/>
        <ci> ghk </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> V </ci>
              <apply>
                <power/>
                <apply>
                  <times/>
                  <cn cellml:units="dimensionless"> 2.0 </cn>
                  <ci> F </ci>
                </apply>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Cao </ci>
              <apply>
                <times/>
                <ci> Ca0 </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 2.0 </cn>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_channel_m_gate">
    <variable units="dimensionless" public_interface="out" name="m"/>
    
    <variable units="per_millisecond" name="alpha_m"/>
    <variable units="per_millisecond" name="beta_m"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="calcium_channel_m_gate_alpha_m_calculation">
        <eq/>
        <ci> alpha_m </ci>
        <apply>
          <divide/>
          <cn cellml:units="per_millisecond"> 5.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> -0.075 </cn>
                <apply>
                  <minus/>
                  <ci> V </ci> 
                  <cn cellml:units="millivolt"> 5.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="calcium_channel_m_gate_beta_m_calculation">
        <eq/>
        <ci> beta_m </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millisecond"> 0.05 </cn>
            <apply>
              <minus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> V </ci> 
                  <cn cellml:units="millivolt"> 2.0 </cn>
                </apply>
                <cn cellml:units="millivolt"> 8.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="calcium_channel_m_gate_dm_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> m </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_m </ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> m </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> beta_m </ci>
            <ci> m </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_channel_h_gate">
    <variable units="dimensionless" public_interface="out" name="h"/>
    
    <variable units="dimensionless" name="Fr"/>
    <variable units="dimensionless" name="Ff" initial_value="0.79"/>
    <variable units="dimensionless" name="Fm" initial_value="0.0"/>
    <variable units="dimensionless" name="Fs" initial_value="0.0"/>
    <variable units="millisecond" name="tau_f" initial_value="500.0"/>
    <variable units="millisecond" name="tau_m"/>
    <variable units="millisecond" name="tau_s"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="h_calculation">
        <eq/>
        <ci> h </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> Ff </ci>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/> 
                  <ci> time </ci>
                </apply>
                <ci> tau_f </ci> 
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> Fm </ci>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/> 
                  <ci> time </ci>
                </apply>
                <ci> tau_m </ci> 
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> Fs </ci>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/> 
                  <ci> time </ci>
                </apply>
                <ci> tau_s </ci> 
              </apply>
            </apply>
          </apply>
          <ci> Fr </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_dependent_potassium_channel">
    <variable units="microA_per_cm2" public_interface="out" name="i_SK"/>

    <variable units="milliS_per_cm2" name="g_SK" initial_value="3.0"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="millimolar" public_interface="in" private_interface="out" name="Ca0"/>
    <variable units="millivolt" public_interface="in" name="E_K"/>
    
    <variable units="dimensionless" private_interface="in" name="m"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_SK_calculation">
        <eq/>
        <ci> i_SK </ci>
        <apply>
          <times/>
          <ci> g_SK </ci>
          <apply>
            <power/>
            <ci> m </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_dependent_potassium_channel_m_gate">
    <variable units="dimensionless" public_interface="out" name="m"/>
    
    <variable units="dimensionless" name="m_infinity"/>
    <variable units="metre_per_second" name="beta_k" initial_value="0.1"/>
    <variable units="dimensionless" name="n_SK" initial_value="1.7"/>
    <variable units="millisecond" name="tau_m"/>
    <variable units="millimolar" name="Car"/>
    <variable units="nanomolar" name="Cac" initial_value="700.0"/>
    
    <variable units="millimolar" public_interface="in" name="Ca0"/>
    <variable units="millisecond" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="m_infinity_calculation">
        <eq/>
        <ci> m_infinity </ci>
        <apply>
          <divide/>
          <ci> Car </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <ci> Car </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="tau_m_calculation">
        <eq/>
        <ci> tau_m </ci>
        <apply>
          <divide/>
          <cn cellml:units="millisecond"> 1.0 </cn>
          <apply>
            <times/>
            <ci> beta_k </ci>
            <apply>
              <plus/>
              <cn cellml:units="millimolar"> 1.0 </cn>
              <ci> Car </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="Car_calculation">
        <eq/>
        <ci> Car </ci>
        <apply>
          <power/>
          <apply>
            <divide/>
            <ci> Ca0 </ci>
            <ci> Cac </ci> 
          </apply>
          <ci> n_SK </ci> 
        </apply>
      </apply>
      
      <apply id="calcium_dependent_potassium_channel_m_gate_dm_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> m </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <minus/>
              <ci> m </ci>
              <ci> m_infinity </ci>
            </apply>
          </apply>
          <ci> tau_m </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="buffered_calcium_diffusion">
    <variable units="millimolar" public_interface="out" name="Bu"/>
    <variable units="millimolar" public_interface="out" name="CaBu"/>
    <variable units="per_millimolar_per_millisecond" public_interface="out" name="k1b" initial_value="100.0"/>
    <variable units="per_millisecond" public_interface="out" name="k2b" initial_value="10.0"/>
    <variable units="per_millisecond" public_interface="out" name="kfc"/>
    <variable units="per_millisecond" public_interface="out" name="kbc"/>
    
    <variable units="millimolar" name="Cai"/>
    <variable units="millimolar" name="Cai2"/>
    <variable units="micrometre2_per_millisecond" name="DCa" initial_value="0.22"/>
    <variable units="micromolar_per_millisecond" name="FCai"/>
    <variable units="millimolar" name="TBu" initial_value="4.0"/>
    
    <variable units="micrometre" public_interface="in" name="delta_r"/>
    <variable units="micrometre2" public_interface="in" name="Ai2"/>
    <variable units="micrometre3" public_interface="in" name="vol_i"/>
     
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Cai_calculation">
        <eq/>
        <ci> Cai </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> CaBu </ci>
            <ci> k2b </ci>
          </apply>
          <apply>
            <times/>
            <ci> Cai </ci>
            <ci> Bu </ci>
            <ci> k1b </ci>
          </apply> 
        </apply>
      </apply>
      
      <apply id="Bu_calculation">
        <eq/>
        <ci> Bu </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> CaBu </ci>
            <ci> k2b </ci>
          </apply>
          <apply>
            <times/>
            <ci> Cai </ci>
            <ci> Bu </ci>
            <ci> k1b </ci>
          </apply> 
        </apply>
      </apply>
      
      <apply id="CaBu_calculation">
        <eq/>
        <ci> CaBu </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> Cai </ci>
            <ci> Bu </ci>
            <ci> k1b </ci>
          </apply>
          <apply>
            <times/>
            <ci> CaBu </ci>
            <ci> k2b </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="FCai_calculation">
        <eq/>
        <ci> FCai </ci>
        <apply>
          <times/>
          <apply>
            <minus/>
            <ci> Cai2 </ci>
            <ci> Cai </ci>
          </apply>
          <ci> DCa </ci>
          <apply>
            <divide/>
            <ci> Ai2 </ci>
            <ci> delta_r </ci>
          </apply> 
        </apply>
      </apply>
      
      <apply id="kfc_calculation">
        <eq/>
        <ci> kfc </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> DCa </ci>
            <ci> Ai2 </ci>
          </apply>
          <apply>
            <times/>
            <ci> delta_r </ci>
            <ci> vol_i </ci>
          </apply> 
        </apply>
      </apply>
      
      <apply id="kbc_calculation">
        <eq/>
        <ci> kbc </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> DCa </ci>
            <ci> Ai2 </ci>
          </apply>
          <apply>
            <times/>
            <ci> delta_r </ci>
            <ci> vol_i </ci>
          </apply> 
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_extrusion">
    <variable units="microA_per_cm2" public_interface="out" name="i_ex"/>
    
    <variable units="microA_per_cm2" name="Vmax" initial_value="3.0"/>
    <variable units="nanomolar" name="Kex" initial_value="830.0"/>
    <variable units="nanomolar" name="Ca_rest" initial_value="100.0"/>
    
    <variable units="millimolar" public_interface="in" name="Ca0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_ex_calculation">
        <eq/>
        <ci> i_ex </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> Vmax </ci>
            <apply>
              <divide/>
              <ci> Ca0 </ci>
              <apply>
                <plus/>
                <ci> Kex </ci>
                <ci> Ca0 </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> Vmax </ci>
            <apply>
              <divide/>
              <ci> Ca_rest </ci>
              <apply>
                <plus/>
                <ci> Kex </ci>
                <ci> Ca_rest </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_triggered_exocytosis">
    <variable units="millimolar_per_millisecond" name="SR"/>
    <variable units="per_millisecond" name="k1f"/>
    <variable units="per_millisecond" name="k1b"/>
    <variable units="per_millisecond" name="k2"/>
    <variable units="per_millisecond" name="a1" initial_value="9.0E-6"/>
    <variable units="per_millisecond" name="a2" initial_value="5.5E-6"/>
    <variable units="per_micromolar3_per_millisecond" name="a3" initial_value="2.3E-5"/>
    <variable units="micromolar" name="b1" initial_value="1.2"/>
    <variable units="dimensionless" name="B"/>
    <variable units="dimensionless" name="A" initial_value="2000.0"/>
    <variable units="dimensionless" name="C"/>
    
    <variable units="millimolar" public_interface="in" name="Ca0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="k1f_calculation">
        <eq/>
        <ci> k1f </ci>
        <apply>
          <times/>
          <ci> a1 </ci>
          <apply>
            <divide/>
            <ci> Ca0 </ci>
            <apply>
              <plus/>
              <ci> b1 </ci>
              <ci> Ca0 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="k1b_calculation">
        <eq/>
        <ci> k1b </ci>
        <ci> a2 </ci>
      </apply>
      
      <apply id="k2_calculation">
        <eq/>
        <ci> k2 </ci>
        <apply>
          <times/>
          <ci> a3 </ci>
          <apply>
            <power/>
            <ci> Ca0 </ci>
            <cn cellml:units="dimensionless"> 3.0 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="B_calculation">
        <eq/>
        <ci> B </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> A </ci>
            <ci> k1f </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> B </ci>
              <ci> k1b </ci>
            </apply>
            <apply>
              <times/>
              <ci> B </ci>
              <ci> k2 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="C_calculation">
        <eq/>
        <ci> C </ci>
        <apply>
          <times/>
          <ci> B </ci>
          <ci> k2 </ci>
        </apply>
      </apply>
      
      <apply id="SR_calculation">
        <eq/>
        <ci> SR </ci>
        <apply>
          <times/>
          <ci> k2 </ci>
          <ci> B </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_indicator">
    <variable units="millimolar" public_interface="out" name="CaDy"/>
    <variable units="millimolar" public_interface="out" name="Dy"/>
    <variable units="per_millimolar_per_millisecond" public_interface="out" name="k1d" initial_value="500.0"/>
    <variable units="per_millisecond" public_interface="out" name="k2d"/>
    
    <variable units="millimolar" name="FDy"/>
    <variable units="millimolar" name="Dyi2"/>
    <variable units="millimolar" name="Cai"/>
    <variable units="millimolar" name="FCaDy"/>
    <variable units="micromolar" name="TDy" initial_value="5.0"/>
    <variable units="micrometre2_per_millisecond" name="Ddy" initial_value="0.05"/>
    <variable units="nanomolar" name="Dkd" initial_value="240.0"/>
    <variable units="dimensionless" name="F0" initial_value="0.36"/>
    <variable units="per_millisecond" name="kfd"/>
    <variable units="per_millisecond" name="kbd"/>
    
    <variable units="micrometre" public_interface="in" name="delta_r"/>
    <variable units="micrometre2" public_interface="in" name="Ai2"/>
    <variable units="micrometre3" public_interface="in" name="vol_i"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Dy_calculation">
        <eq/>
        <ci> Dy </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> Dyi2 </ci>
            <ci> kbd </ci>
          </apply>
          <apply>
            <times/>
            <ci> Dy </ci>
            <ci> kfd </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="kfd_calculation">
        <eq/>
        <ci> kfd </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> Ddy </ci>
            <ci> Ai2 </ci>
          </apply>
          <apply>
            <times/>
            <ci> delta_r </ci>
            <ci> vol_i </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="kbd_calculation">
        <eq/>
        <ci> kbd </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> Ddy </ci>
            <ci> Ai2 </ci>
          </apply>
          <apply>
            <times/>
            <ci> delta_r </ci>
            <ci> vol_i </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="CaDy_calculation">
        <eq/>
        <ci> CaDy </ci>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> Dy </ci>
            <ci> Cai </ci>
            <ci> k1d </ci>
          </apply>
          <apply>
            <times/>
            <ci> CaDy </ci>
            <ci> k2d </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="FDy_calculation">
        <eq/>
        <ci> FDy </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <ci> F0 </ci>
            </apply>
            <ci> FCaDy </ci>
          </apply>
          <ci> F0 </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_concentration">
    <variable units="millimolar" public_interface="out" name="Ca0"/>
              
    <variable units="micrometre2" name="A0"/>
    <variable units="micrometre3" name="vol0"/>
    <variable units="millimolar" name="Ca1"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_Ca"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_ex"/>
    <variable units="millimolar" public_interface="in" name="CaBu"/>
    <variable units="millimolar" public_interface="in" name="Bu"/>
    <variable units="millimolar" public_interface="in" name="CaDy"/>
    <variable units="millimolar" public_interface="in" name="Dy"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    <variable units="per_millimolar_per_millisecond" public_interface="in" name="k1d"/>
    <variable units="per_millisecond" public_interface="in" name="k2d"/>
    <variable units="per_millimolar_per_millisecond" public_interface="in" name="k1b"/>
    <variable units="per_millisecond" public_interface="in" name="k2b"/>
    <variable units="per_millisecond" public_interface="in" name="kfc"/>
    <variable units="per_millisecond" public_interface="in" name="kbc"/>
     
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Ca0_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca0 </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> kbc  </ci>
              <ci> Ca1 </ci>
            </apply>
            <apply>
              <times/>
              <ci> k2b  </ci>
              <ci> CaBu </ci>
            </apply>
            <apply>
              <times/>
              <ci> k2d  </ci>
              <ci> CaDy </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <divide/>
              <apply>
                <times/>
                <apply>
                  <plus/>
                  <ci> i_Ca </ci>
                  <ci> i_ex </ci>
                </apply>
                <ci> A0 </ci>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 2.0 </cn>
                <ci> F </ci>
                <ci> vol0 </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> kfc  </ci>
              <ci> Ca0 </ci>
            </apply>
            <apply>
              <times/>
              <ci> k1b  </ci>
              <ci> Bu </ci>
              <ci> Ca0 </ci>
            </apply>
            <apply>
              <times/>
              <ci> k1d  </ci>
              <ci> Dy </ci>
              <ci> Ca0 </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="model_parameters">
    <variable units="joule_per_mole_kelvin" public_interface="out" name="R" initial_value="8.3144"/>
    <variable units="kelvin" public_interface="out" name="T" initial_value="310.15"/>
    <variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96485"/>
    <variable units="micrometre" public_interface="out" name="delta_r"/>
    <variable units="micrometre2" public_interface="out" name="Ai2"/>
    <variable units="micrometre3" public_interface="out" name="vol_i"/>
    <variable units="millivolt" public_interface="out" name="E_K" initial_value="-80.0"/>
    
    <variable units="micrometre" name="diameter" initial_value="15.0"/>
    <variable units="dimensionless" name="Nc" initial_value="6.0"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="delta_r_calculation">
        <eq/>
        <ci> delta_r </ci>
        <apply>
          <divide/>
          <ci> diameter </ci>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 2.0 </cn>
            <apply>
              <minus/>
              <ci> Nc </ci>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply> 
        </apply>
      </apply>
    </math>
  </component>
  
  
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="excitatory_synapse"/>
      <component_ref component="calcium_channel">
        <component_ref component="calcium_channel_m_gate"/>
        <component_ref component="calcium_channel_h_gate"/>
      </component_ref>
      <component_ref component="sodium_channel">
        <component_ref component="sodium_channel_m_gate"/>
        <component_ref component="sodium_channel_h_gate"/>
      </component_ref>
      <component_ref component="potassium_channel">
        <component_ref component="potassium_channel_n_gate"/>
      </component_ref>
      <component_ref component="calcium_dependent_potassium_channel">
        <component_ref component="calcium_dependent_potassium_channel_m_gate"/>
      </component_ref>
      <component_ref component="buffered_calcium_diffusion"/>
      <component_ref component="calcium_extrusion"/>
      <component_ref component="calcium_triggered_exocytosis"/>
      <component_ref component="calcium_indicator"/>
      <component_ref component="calcium_concentration"/>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="encapsulation"/>
    <component_ref component="calcium_channel">
      <component_ref component="calcium_channel_m_gate"/>
      <component_ref component="calcium_channel_h_gate"/>
    </component_ref>
    <component_ref component="sodium_channel">
      <component_ref component="sodium_channel_m_gate"/>
      <component_ref component="sodium_channel_h_gate"/>
    </component_ref>
    <component_ref component="potassium_channel">
      <component_ref component="potassium_channel_n_gate"/>
    </component_ref>
    <component_ref component="calcium_dependent_potassium_channel">
      <component_ref component="calcium_dependent_potassium_channel_m_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="sodium_channel"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="potassium_channel"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="calcium_channel"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="calcium_dependent_potassium_channel"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="calcium_concentration"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="excitatory_synapse" component_1="membrane"/>
    <map_variables variable_2="i_sp" variable_1="i_sp"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_channel" 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="sodium_channel" component_1="membrane"/>
    <map_variables variable_2="i_Na" variable_1="i_Na"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_channel" 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="calcium_dependent_potassium_channel" component_1="membrane"/>
    <map_variables variable_2="i_SK" variable_1="i_SK"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="leak_channel" component_1="membrane"/>
    <map_variables variable_2="i_l" variable_1="i_l"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_extrusion" component_1="membrane"/>
    <map_variables variable_2="i_ex" variable_1="i_ex"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_channel" component_1="model_parameters"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
  </connection>
  
  <connection>
    <map_components component_2="buffered_calcium_diffusion" component_1="model_parameters"/>
    <map_variables variable_2="delta_r" variable_1="delta_r"/>
    <map_variables variable_2="Ai2" variable_1="Ai2"/>
    <map_variables variable_2="vol_i" variable_1="vol_i"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_indicator" component_1="model_parameters"/>
    <map_variables variable_2="delta_r" variable_1="delta_r"/>
    <map_variables variable_2="Ai2" variable_1="Ai2"/>
    <map_variables variable_2="vol_i" variable_1="vol_i"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_concentration" component_1="model_parameters"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_channel" component_1="model_parameters"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_dependent_potassium_channel" component_1="model_parameters"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_channel" component_1="calcium_concentration"/>
    <map_variables variable_2="Ca0" variable_1="Ca0"/>
    <map_variables variable_2="i_Ca" variable_1="i_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_dependent_potassium_channel" component_1="calcium_concentration"/>
    <map_variables variable_2="Ca0" variable_1="Ca0"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_extrusion" component_1="calcium_concentration"/>
    <map_variables variable_2="Ca0" variable_1="Ca0"/>
    <map_variables variable_2="i_ex" variable_1="i_ex"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_triggered_exocytosis" component_1="calcium_concentration"/>
    <map_variables variable_2="Ca0" variable_1="Ca0"/>
  </connection>
  
  <connection>
    <map_components component_2="buffered_calcium_diffusion" component_1="calcium_concentration"/>
    <map_variables variable_2="Bu" variable_1="Bu"/>
    <map_variables variable_2="CaBu" variable_1="CaBu"/>
    <map_variables variable_2="k1b" variable_1="k1b"/>
    <map_variables variable_2="k2b" variable_1="k2b"/>
    <map_variables variable_2="kfc" variable_1="kfc"/>
    <map_variables variable_2="kbc" variable_1="kbc"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_indicator" component_1="calcium_concentration"/>
    <map_variables variable_2="Dy" variable_1="Dy"/>
    <map_variables variable_2="CaDy" variable_1="CaDy"/>
    <map_variables variable_2="k1d" variable_1="k1d"/>
    <map_variables variable_2="k2d" variable_1="k2d"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_channel_m_gate" component_1="sodium_channel"/>
    <map_variables variable_2="m" variable_1="m"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_channel_h_gate" component_1="sodium_channel"/>
    <map_variables variable_2="h" variable_1="h"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_channel_n_gate" component_1="potassium_channel"/>
    <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="calcium_channel_m_gate" component_1="calcium_channel"/>
    <map_variables variable_2="m" variable_1="m"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_channel_h_gate" component_1="calcium_channel"/>
    <map_variables variable_2="h" variable_1="h"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_dependent_potassium_channel_m_gate" component_1="calcium_dependent_potassium_channel"/>
    <map_variables variable_2="m" variable_1="m"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="Ca0" variable_1="Ca0"/>
  </connection>
  

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#0224cf43-6d94-472a-89a9-2e966d999857">
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>signal transduction</rdf:li>
    <rdf:li>stimulation-secretion coupling</rdf:li>
    <rdf:li>chromaffin cell</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#081b814b-71a4-438d-b720-41eabfff7b16">
    <rdf:li rdf:resource="rdf:#d74d9c39-ab5b-49d5-b856-f1a28ca59cc1"/>
    <rdf:li rdf:resource="rdf:#e35a67c4-aac4-4018-8b04-aed939d1eced"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#0966ebeb-54fb-4746-b98b-6f9fc506dd27">
    <dc:creator rdf:resource="rdf:#a0968c21-c7c3-401a-bd14-4741694ca78c"/>
    <rdf:value>
          This is the CellML description of Warashina and Ogura's 2004 mathematical model of stimulation-secretion coupling in a chromaffin cell.
        </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d74d9c39-ab5b-49d5-b856-f1a28ca59cc1">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#1c8c2a39-870a-4d9f-b57a-0560dad53f24"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#24d62153-a9d6-4420-8f13-63bcbbd4e66b">
    <dc:subject rdf:resource="rdf:#023d6573-2506-4f3e-b2f6-8bec53bbe45f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6335ff5d-092d-47c1-ae58-2471c971d717">
    <vCard:ORG rdf:resource="rdf:#cbd53560-cc49-421c-bb8d-db797070f95a"/>
    <vCard:EMAIL rdf:resource="rdf:#85eb1f66-4657-40f7-8573-fba6cf30f443"/>
    <vCard:N rdf:resource="rdf:#91df4d46-3d07-42b0-94fc-320129f7ab58"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#92f97b4a-dc37-4059-b7e5-f90413ff9437">
    <dcterms:W3CDTF>2004-02-03</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5cf3148b-cdba-4601-95f2-6930779a2b63">
    <dc:title>Pflugers Archive European Journal of Physiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#85eb1f66-4657-40f7-8573-fba6cf30f443">
    <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:#a0968c21-c7c3-401a-bd14-4741694ca78c">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#023d6573-2506-4f3e-b2f6-8bec53bbe45f">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#0224cf43-6d94-472a-89a9-2e966d999857"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1c8c2a39-870a-4d9f-b57a-0560dad53f24">
    <vCard:Given>A</vCard:Given>
    <vCard:Family>Warashina</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ecf648fb-51d9-40ab-b3b3-4825547dc7ff">
    <bqs:Pubmed_id>14730449</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#7cf1a76d-def5-4d34-a72c-a5adc979764d"/>
  </rdf:Description>
  <rdf:Description rdf:about="#warashina_ogura_2004_version01">
    <dc:title>
        Warashina and Ogura's 2004 mathematical model of stimulation-secretion coupling in a chromaffin cell.
      </dc:title>
    <cmeta:bio_entity>chromaffin cell</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#0966ebeb-54fb-4746-b98b-6f9fc506dd27"/>
    <bqs:reference rdf:resource="rdf:#24d62153-a9d6-4420-8f13-63bcbbd4e66b"/>
    <bqs:reference rdf:resource="rdf:#ecf648fb-51d9-40ab-b3b3-4825547dc7ff"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#83fea202-fde3-4bed-a9bd-475dffc73a17">
    <dcterms:W3CDTF>2004-01-17</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cbd53560-cc49-421c-bb8d-db797070f95a">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#101374dd-7674-4d77-85a3-fcaf5acc0644">
    <vCard:Given>T</vCard:Given>
    <vCard:Family>Ogura</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#91df4d46-3d07-42b0-94fc-320129f7ab58">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7cf1a76d-def5-4d34-a72c-a5adc979764d">
    <dc:creator rdf:resource="rdf:#081b814b-71a4-438d-b720-41eabfff7b16"/>
    <dc:title>
            Modeling of simulation-secretion coupling in a chromaffin cell
          </dc:title>
    <bqs:volume>Online First</bqs:volume>
    <bqs:first_page>1</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#5cf3148b-cdba-4601-95f2-6930779a2b63"/>
    <dcterms:issued rdf:resource="rdf:#83fea202-fde3-4bed-a9bd-475dffc73a17"/>
    <bqs:last_page>8</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e35a67c4-aac4-4018-8b04-aed939d1eced">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#101374dd-7674-4d77-85a3-fcaf5acc0644"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Bioengineering Institute
      </dc:publisher>
    <dcterms:created rdf:resource="rdf:#92f97b4a-dc37-4059-b7e5-f90413ff9437"/>
    <dc:creator rdf:resource="rdf:#6335ff5d-092d-47c1-ae58-2471c971d717"/>
  </rdf:Description>
</rdf:RDF>
</model>