Rendering of the source text

<?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_2004" name="warashina_2004">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modeling of stimulation-secretion coupling in a chromaffin cell</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
         This CellML model description is currently incomplete. In order to fully describe the original model it would be necessary to create 6 different compartments and then decribe the diffusion of calcium between these.  Further there are other parameter values which are missing.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: We constructed a chromaffin cell model for analysis of stimulation-secretion coupling in computer simulation studies. The model includes mechanisms involved in the excitatory synapse, voltage-dependent Na(+), K(+) and Ca(2+) channels, Ca(2+)-activated K(+) channels (SK type), buffered Ca(2+) diffusion, Ca(2+) extrusion, fluorescent Ca(2+) indicators and Ca(2+)-triggered exocytosis. Calculations of the modeled mechanisms were carried out using the NEURON simulation environment (Hines and Carnevale, Neural Computation 9:1179-1209, 1997). A set of parameter values was determined so as to fit basic experimental results reported in the literature. The model was also applied to simulate our experimental results obtained from chromaffin cells in the perfused rat adrenal medulla. Observed profiles of Ca(2+)responses induced by electrically stimulating the splanchnic nerve with various frequencies (1-50 Hz) were adequately simulated with minor readjustments of parameter values for Ca(2+)influx and extrusion. Secretory responses measured at the same time as the Ca(2+)responses were also simulated with consideration of a time constant to detect catecholamines in the experiment. Similarly, model simulations reproduced both Ca(2+)responses and secretory responses evoked by elevations of the extracellular K(+) concentration for different periods. The results suggest that the presented model provides a useful tool for analyzing and predicting quantitative relations in various events occurring in stimulation-secretion coupling in chromaffin cells.
</para>

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

<para>
Modeling of stimulation-secretion coupling in a chromaffin cell, A. Warashina and T. Ogura, 2004, <emphasis>European Journal of Physiology</emphasis>, 448, 369-391.  <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>

</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" initial_value="-75.0"/>
              
    <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="microS" name="G" initial_value="0.001"/>
    <variable units="microS" name="L" initial_value="0.016"/>
    <variable units="microS" name="D" initial_value="0.000"/>
    <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"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <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_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> G </ci>
        </apply>
		<apply>
          <minus/>
          <apply>
            <times/>
            <ci> L </ci>
            <apply>
              <divide/>
			  <cn cellml:units="dimensionless"> 1 </cn>
			  <ci> tau1 </ci>
			</apply>
          </apply>
          <apply>
            <times/>
            <ci> D </ci>
            <apply>
              <divide/>
			  <cn cellml:units="dimensionless"> 1 </cn>
			  <ci> tau2 </ci>
			</apply>
          </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" initial_value="0.001"/>
    
    <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" initial_value="0.001"/>
    
    <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" initial_value="0.001"/>
    
    <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" initial_value="0.001"/>
    
    <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/>
              <cn cellml:units="millivolt"> 2.0 </cn>
			  <ci> V </ci>
            </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" initial_value="0.000"/>
    <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" initial_value="0.000"/>
    <variable units="millisecond" name="tau_s" initial_value="0.000"/>
    
    <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" initial_value="0.001"/>
    
    <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="micrometre2_per_millisecond" name="DCa" initial_value="0.22"/>
    <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="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"  initial_value="0.12"/>
    
    <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" initial_value="0.1"/>
              
    <variable units="micrometre2" name="A0" initial_value="706.86"/>
    <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" initial_value="1767.15"/>
    <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="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="excitatory_synapse"/>
    <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_2004">
    <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>