Location: Purvis, Butera, 2005 @ f33430047c6f / purvis_butera_2005.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2008-07-17 02:18:21+12:00
Desc:
committing version02 of purvis_butera_2005
Permanent Source URI:
http://models.cellml.org/workspace/purvis_butera_2005/rawfile/f33430047c6f6c2f195381b678b87499bd2f3475/purvis_butera_2005.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 17/07/2008 at 2:14:10 at p.m. using:

COR (0.9.31.979)
Copyright 2002-2008 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk

CellML 1.0 was used to generate this model
http://www.CellML.org/
--><model xmlns="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="butera_purvis_2004" name="butera_purvis_2004">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Ionic Current Model of Hypoglossal Motoneuron</title>
				<author>
					<firstname>Geoffrey</firstname>
					<surname>Nunns</surname>
					<affiliation>
						<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
           This CellML model is known to run in both PCEnv and COR and the units have been checked and are consistent. However the CellML model is curently unable to perfectly recreate the published results (we are several orders of magnitude out with certain variables). Further curation is required and we are endevouring to seek the help of the original model authors.
          </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
Abstract: The hypoglossal motoneuron (HM) innervates the tongue and along with its roles in mastication and swallowing, the tongue plays an important role in breathing (Bartlett <emphasis>et al.</emphasis> 1990; Lowe 1980; reviewed in Sawczuk and Mosier 2001). The muscle of the tongue contributes to effective breathing by maintaining a patent upper airway and problems with the control of HMs can cause the breathing disorder obstructive sleep apnea (Fogel <emphasis>et al.</emphasis> 2000; Remmers <emphasis>et al.</emphasis> 1978). For this reason, extensive experimental work has been done on these motoneurons. In this communication, Butera and Smith present a computational model of a HM that provides a framework for analyzing the role of ionic currents in these neurons.
</para>
				<para>  
In this paper, Liston Purvis and Robert Butera have developed a single-compartment, electrophysiological, hypoglossal motoneuron (HM) model based primarily on experimental data from neonatal rat HMs. The model is able to reproduce the fine features of the HM action potential: the fast afterhyperpolarization, the afterdepolarization, and the medium-duration afterhyperpolarization (mAHP). The model also reproduces the repetitive firing properties seen in neonatal HMs and replicates the neuron's response to pharmacological experiments. The model was used to study the role of specific ionic currents in HM firing and how variations in the densities of these currents may account for age-dependent changes in excitability seen in HMs. By varying the density of a fast inactivating calcium current, the model alternates between accelerating and adapting firing patterns. Modeling the age-dependent increase in H current density accounts for the decrease in mAHP duration observed experimentally, but does not fully account for the decrease in input resistance. An increase in the density of the voltage-dependent potassium currents and the H current is required to account for the decrease in input resistance. These changes also account for the age-dependent decrease in action potential duration.</para>
				<para>
The complete original paper reference is cited below:
</para>
				<para>
					<ulink url="http://jn.physiology.org/cgi/content/abstract/93/2/723">Ionic Current Model of a Hypoglossal Motoneuron</ulink>, Robert J. Butera, Jr., Liston K. Purvis, 2005, <ulink url="http://jn.physiology.org/">
						<emphasis>Journal of Neurophysiology</emphasis>
					</ulink>, 93, 723-733.  (<ulink url="http://intl-jn.physiology.org/cgi/content/full/93/2/723">Full text</ulink> and <ulink url="http://jn.physiology.org/cgi/reprint/93/2/723">PDF</ulink> versions of the article are available for Journal Members on the JN website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=15653786&amp;dopt=Abstract">PubMed ID: 15653786</ulink>
				</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="purvis_2005.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the CellML model.</caption>
</informalfigure>

			</sect1>
		</article>
	</documentation>   <units name="millisecond">
      <unit units="second" prefix="milli"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="nanoA">
      <unit units="ampere" prefix="nano"/>
   </units>
   <units name="uS">
      <unit units="siemens" prefix="micro"/>
   </units>
   <units name="uF">
      <unit units="farad" prefix="nano"/>
   </units>
   <units name="uM">
      <unit units="mole" prefix="micro"/>
      <unit units="liter" exponent="-1"/>
   </units>
   <units name="per_ms">
      <unit units="millisecond" exponent="-1"/>
   </units>
   <units name="uM_per_nanocoulomb">
      <unit units="uM"/>
      <unit units="coulomb" prefix="nano" 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="-71.847"/>
      <variable units="uF" name="C" initial_value="0.04"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="nanoA" public_interface="in" name="i_Na"/>
      <variable units="nanoA" public_interface="in" name="i_NaP"/>
      <variable units="nanoA" public_interface="in" name="i_K"/>
      <variable units="nanoA" public_interface="in" name="i_leak"/>
      <variable units="nanoA" public_interface="in" name="i_T"/>
      <variable units="nanoA" public_interface="in" name="i_N"/>
      <variable units="nanoA" public_interface="in" name="i_P"/>
      <variable units="nanoA" public_interface="in" name="i_SK"/>
      <variable units="nanoA" public_interface="in" name="i_A"/>
      <variable units="nanoA" public_interface="in" name="i_H"/>
      <variable units="nanoA" public_interface="in" name="i_app"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>i_Na</ci>
                        <ci>i_NaP</ci>
                        <ci>i_K</ci>
                        <ci>i_leak</ci>
                        <ci>i_T</ci>
                        <ci>i_N</ci>
                        <ci>i_P</ci>
                        <ci>i_SK</ci>
                        <ci>i_A</ci>
                        <ci>i_H</ci>
                     </apply>
                  </apply>
                  <ci>i_app</ci>
               </apply>
               <ci>C</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current">
      <variable units="nanoA" public_interface="out" name="i_Na"/>
      <variable units="uS" name="g_Na" initial_value="0.7"/>
      <variable units="millivolt" public_interface="out" name="E_Na" initial_value="60"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <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>
            <eq/>
            <ci>i_Na</ci>
            <apply>
               <times/>
               <ci>g_Na</ci>
               <apply>
                  <power/>
                  <ci>m</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current_h_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="h" initial_value="0.981"/>
      <variable units="millivolt" name="theta_h" initial_value="44.1"/>
      <variable units="millivolt" name="sigma_h" initial_value="7"/>
      <variable units="millivolt" name="theta_1" initial_value="35"/>
      <variable units="millivolt" name="sigma_1" initial_value="4"/>
      <variable units="millivolt" name="sigma_2" initial_value="25"/>
      <variable units="millisecond" name="tau_h"/>
      <variable units="dimensionless" name="h_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau_h</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="millisecond">3.5</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_1</ci>
                           </apply>
                           <ci>sigma_1</ci>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <ci>theta_1</ci>
                              </apply>
                           </apply>
                           <ci>sigma_2</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millisecond">1</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>h_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <ci>theta_h</ci>
                        </apply>
                        <ci>sigma_h</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h_infinity</ci>
                  <ci>h</ci>
               </apply>
               <ci>tau_h</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="sodium_current_m_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.015"/>
      <variable units="millivolt" name="theta_m" initial_value="36"/>
      <variable units="millivolt" name="sigma_m" initial_value="8.5"/>
      <variable units="millisecond" name="tau_m" initial_value="0.1"/>
      <variable units="dimensionless" name="m_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <ci>theta_m</ci>
                        </apply>
                        <apply>
                           <minus/>
                           <ci>sigma_m</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>m_infinity</ci>
                  <ci>m</ci>
               </apply>
               <ci>tau_m</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="persistent_sodium_current">
      <variable units="nanoA" public_interface="out" name="i_NaP"/>
      <variable units="uS" name="g_NaP" initial_value="0.05"/>
      <variable units="millivolt" public_interface="in" name="E_Na"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <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>
            <eq/>
            <ci>i_NaP</ci>
            <apply>
               <times/>
               <ci>g_NaP</ci>
               <ci>m</ci>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="persistent_sodium_current_h_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="h" initial_value="0.797"/>
      <variable units="millivolt" name="theta_h" initial_value="65"/>
      <variable units="millivolt" name="sigma_h" initial_value="5"/>
      <variable units="millisecond" name="tau_h" initial_value="150"/>
      <variable units="dimensionless" name="h_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>h_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <ci>theta_h</ci>
                        </apply>
                        <ci>sigma_h</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h_infinity</ci>
                  <ci>h</ci>
               </apply>
               <ci>tau_h</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="persistent_sodium_current_m_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.002"/>
      <variable units="millivolt" name="theta_m" initial_value="47.1"/>
      <variable units="millivolt" name="sigma_m" initial_value="4.1"/>
      <variable units="millisecond" name="tau_m" initial_value="0.1"/>
      <variable units="dimensionless" name="m_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_m</ci>
                           </apply>
                        </apply>
                        <ci>sigma_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>m_infinity</ci>
                  <ci>m</ci>
               </apply>
               <ci>tau_m</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="delayed_rectifier_current">
      <variable units="nanoA" public_interface="out" name="i_K"/>
      <variable units="uS" name="g_K" initial_value="1.3"/>
      <variable units="millivolt" public_interface="out" name="E_K" initial_value="-80"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="dimensionless" private_interface="in" name="n"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K</ci>
            <apply>
               <times/>
               <ci>g_K</ci>
               <apply>
                  <power/>
                  <ci>n</ci>
                  <cn cellml:units="dimensionless">4</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="delayed_rectifier_current_n_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="n" initial_value="0.158"/>
      <variable units="millivolt" name="theta_n" initial_value="30"/>
      <variable units="millivolt" name="sigma_n" initial_value="25"/>
      <variable units="millivolt" name="theta_1" initial_value="30"/>
      <variable units="millivolt" name="sigma_1" initial_value="40"/>
      <variable units="millivolt" name="sigma_2" initial_value="50"/>
      <variable units="millisecond" name="tau_n"/>
      <variable units="dimensionless" name="n_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau_n</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="millisecond">2.5</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_1</ci>
                           </apply>
                           <ci>sigma_1</ci>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <ci>theta_1</ci>
                              </apply>
                           </apply>
                           <ci>sigma_2</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millisecond">0.01</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>n_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_n</ci>
                           </apply>
                        </apply>
                        <ci>sigma_n</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>n</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>n_infinity</ci>
                  <ci>n</ci>
               </apply>
               <ci>tau_n</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="leak_current">
      <variable units="nanoA" public_interface="out" name="i_leak"/>
      <variable units="uS" name="g_leak" initial_value="0.005"/>
      <variable units="millivolt" name="E_leak" initial_value="-50"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_leak</ci>
            <apply>
               <times/>
               <ci>g_leak</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_leak</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="LVA_calcium_current">
      <variable units="nanoA" public_interface="out" name="i_T"/>
      <variable units="uS" name="g_T" initial_value="0.1"/>
      <variable units="millivolt" public_interface="out" name="E_Ca" initial_value="40"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <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>
            <eq/>
            <ci>i_T</ci>
            <apply>
               <times/>
               <ci>g_T</ci>
               <ci>m</ci>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="LVA_calcium_current_m_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.001"/>
      <variable units="millivolt" name="theta_m" initial_value="38"/>
      <variable units="millivolt" name="sigma_m" initial_value="5"/>
      <variable units="millivolt" name="theta_1" initial_value="28"/>
      <variable units="millivolt" name="sigma_1" initial_value="25"/>
      <variable units="millivolt" name="sigma_2" initial_value="70"/>
      <variable units="millisecond" name="tau_m"/>
      <variable units="dimensionless" name="m_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau_m</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="millisecond">5</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_1</ci>
                           </apply>
                           <ci>sigma_1</ci>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <ci>theta_1</ci>
                              </apply>
                           </apply>
                           <ci>sigma_2</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millisecond">2</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_m</ci>
                           </apply>
                        </apply>
                        <ci>sigma_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>m_infinity</ci>
                  <ci>m</ci>
               </apply>
               <ci>tau_m</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="LVA_calcium_current_h_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="h" initial_value="0.562"/>
      <variable units="millivolt" name="theta_h" initial_value="70.1"/>
      <variable units="millivolt" name="sigma_h" initial_value="7"/>
      <variable units="millivolt" name="theta_1" initial_value="70"/>
      <variable units="millivolt" name="sigma_1" initial_value="65"/>
      <variable units="millisecond" name="tau_h"/>
      <variable units="dimensionless" name="h_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau_h</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="millisecond">20</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_1</ci>
                           </apply>
                           <ci>sigma_1</ci>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <ci>theta_1</ci>
                              </apply>
                           </apply>
                           <ci>sigma_1</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millisecond">1</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>h_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <ci>theta_h</ci>
                        </apply>
                        <ci>sigma_h</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h_infinity</ci>
                  <ci>h</ci>
               </apply>
               <ci>tau_h</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="N_HVA_calcium_current">
      <variable units="nanoA" public_interface="out" name="i_N"/>
      <variable units="uS" name="g_N" initial_value="0.05"/>
      <variable units="millivolt" public_interface="in" name="E_Ca"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <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>
            <eq/>
            <ci>i_N</ci>
            <apply>
               <times/>
               <ci>g_N</ci>
               <ci>m</ci>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="N_HVA_calcium_current_m_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.001"/>
      <variable units="millivolt" name="theta_m" initial_value="30"/>
      <variable units="millivolt" name="sigma_m" initial_value="6"/>
      <variable units="millisecond" name="tau_m" initial_value="5"/>
      <variable units="dimensionless" name="m_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_m</ci>
                           </apply>
                        </apply>
                        <ci>sigma_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>m_infinity</ci>
                  <ci>m</ci>
               </apply>
               <ci>tau_m</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="N_HVA_calcium_current_h_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="h" initial_value="0.649"/>
      <variable units="millivolt" name="theta_h" initial_value="70"/>
      <variable units="millivolt" name="sigma_h" initial_value="3"/>
      <variable units="millisecond" name="tau_h" initial_value="25"/>
      <variable units="dimensionless" name="h_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>h_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <ci>theta_h</ci>
                        </apply>
                        <ci>sigma_h</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h_infinity</ci>
                  <ci>h</ci>
               </apply>
               <ci>tau_h</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="P_HVA_calcium_current">
      <variable units="nanoA" public_interface="out" name="i_P"/>
      <variable units="uS" name="g_P" initial_value="0.05"/>
      <variable units="millivolt" public_interface="in" name="E_Ca"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="dimensionless" private_interface="in" name="m"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_P</ci>
            <apply>
               <times/>
               <ci>g_P</ci>
               <ci>m</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="P_HVA_calcium_current_m_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0"/>
      <variable units="millivolt" name="theta_m" initial_value="17"/>
      <variable units="millivolt" name="sigma_m" initial_value="3"/>
      <variable units="millisecond" name="tau_m" initial_value="10"/>
      <variable units="dimensionless" name="m_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_m</ci>
                           </apply>
                        </apply>
                        <ci>sigma_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>m_infinity</ci>
                  <ci>m</ci>
               </apply>
               <ci>tau_m</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_dependent_potassium_current">
      <variable units="nanoA" public_interface="out" name="i_SK"/>
      <variable units="uS" name="g_SK" initial_value="0.3"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="nanoA" public_interface="in" private_interface="out" name="i_T"/>
      <variable units="nanoA" public_interface="in" private_interface="out" name="i_N"/>
      <variable units="nanoA" public_interface="in" private_interface="out" name="i_P"/>
      <variable units="dimensionless" private_interface="in" name="z"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_SK</ci>
            <apply>
               <times/>
               <ci>g_SK</ci>
               <apply>
                  <power/>
                  <ci>z</ci>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_dependent_potassium_current_z_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="nanoA" public_interface="in" name="i_T"/>
      <variable units="nanoA" public_interface="in" name="i_N"/>
      <variable units="nanoA" public_interface="in" name="i_P"/>
      <variable units="dimensionless" public_interface="out" name="z" initial_value="0"/>
      <variable units="uM_per_nanocoulomb" name="K1" initial_value="-500"/>
      <variable units="per_ms" name="K2" initial_value="0.04"/>
      <variable units="uM" name="Ca_conc" initial_value="0.0604"/>
      <variable units="millisecond" name="tau_z" initial_value="1"/>
      <variable units="dimensionless" name="z_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_conc</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <cn cellml:units="second">1</cn>
                     <cn cellml:units="millisecond">1000</cn>
                  </apply>
                  <ci>K1</ci>
                  <apply>
                     <plus/>
                     <ci>i_T</ci>
                     <ci>i_N</ci>
                     <ci>i_P</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>K2</ci>
                  <ci>Ca_conc</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>z_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <power/>
                     <apply>
                        <divide/>
                        <cn cellml:units="uM">0.003</cn>
                        <ci>Ca_conc</ci>
                     </apply>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>z</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>z_infinity</ci>
                  <ci>z</ci>
               </apply>
               <ci>tau_z</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_transient_potassium_current">
      <variable units="nanoA" public_interface="out" name="i_A"/>
      <variable units="uS" name="g_A" initial_value="1"/>
      <variable units="millivolt" public_interface="in" name="E_K"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <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>
            <eq/>
            <ci>i_A</ci>
            <apply>
               <times/>
               <ci>g_A</ci>
               <ci>m</ci>
               <ci>h</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_transient_potassium_current_m_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.057"/>
      <variable units="millivolt" name="theta_m" initial_value="27"/>
      <variable units="millivolt" name="sigma_m" initial_value="16"/>
      <variable units="millivolt" name="theta_1" initial_value="40"/>
      <variable units="millivolt" name="theta_2" initial_value="74"/>
      <variable units="millivolt" name="sigma_1" initial_value="5"/>
      <variable units="millivolt" name="sigma_2" initial_value="7.5"/>
      <variable units="millisecond" name="tau_m"/>
      <variable units="dimensionless" name="m_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau_m</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="millisecond">1</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_1</ci>
                           </apply>
                           <ci>sigma_1</ci>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <ci>theta_2</ci>
                              </apply>
                           </apply>
                           <ci>sigma_2</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millisecond">0.37</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_m</ci>
                           </apply>
                        </apply>
                        <ci>sigma_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>m_infinity</ci>
                  <ci>m</ci>
               </apply>
               <ci>tau_m</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_transient_potassium_current_h_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="h" initial_value="0.287"/>
      <variable units="millivolt" name="theta_h" initial_value="80"/>
      <variable units="millivolt" name="sigma_h" initial_value="11"/>
      <variable units="millisecond" name="tau_h" initial_value="20"/>
      <variable units="dimensionless" name="h_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>h_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <ci>theta_h</ci>
                        </apply>
                        <ci>sigma_h</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h_infinity</ci>
                  <ci>h</ci>
               </apply>
               <ci>tau_h</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarization_activated_current">
      <variable units="nanoA" public_interface="out" name="i_H"/>
      <variable units="uS" name="g_H" initial_value="0.005"/>
      <variable units="millivolt" name="E_H" initial_value="-38.8"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
      <variable units="dimensionless" private_interface="in" name="m"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_H</ci>
            <apply>
               <times/>
               <ci>g_H</ci>
               <ci>m</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_H</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="hyperpolarization_activated_current_m_gate">
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="dimensionless" public_interface="out" name="m" initial_value="0.182"/>
      <variable units="millivolt" name="theta_m" initial_value="79.8"/>
      <variable units="millivolt" name="sigma_m" initial_value="5.3"/>
      <variable units="millivolt" name="theta_1" initial_value="70"/>
      <variable units="millivolt" name="sigma_1" initial_value="11"/>
      <variable units="millisecond" name="tau_m"/>
      <variable units="dimensionless" name="m_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau_m</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn cellml:units="millisecond">1</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <plus/>
                              <ci>V</ci>
                              <ci>theta_1</ci>
                           </apply>
                           <ci>sigma_1</ci>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <apply>
                                 <plus/>
                                 <ci>V</ci>
                                 <ci>theta_1</ci>
                              </apply>
                           </apply>
                           <ci>sigma_1</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn cellml:units="millisecond">50</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <plus/>
                           <ci>V</ci>
                           <ci>theta_m</ci>
                        </apply>
                        <ci>sigma_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>m</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>m_infinity</ci>
                  <ci>m</ci>
               </apply>
               <ci>tau_m</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="stimulus_protocol">
      <variable units="nanoA" public_interface="out" name="i_app"/>
      <variable units="millisecond" name="i_stimStart" initial_value="0.5"/>
      <variable units="millisecond" name="i_stimEnd" initial_value="0.9"/>
      <variable units="nanoA" name="i_stimAmplitude" initial_value="1"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_app</ci>
            <piecewise>
               <piece>
                  <ci>i_stimAmplitude</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>time</ci>
                        <ci>i_stimStart</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <ci>time</ci>
                        <ci>i_stimEnd</ci>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="nanoA">0</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="sodium_current">
         <component_ref component="sodium_current_h_gate"/>
         <component_ref component="sodium_current_m_gate"/>
      </component_ref>
      <component_ref component="persistent_sodium_current">
         <component_ref component="persistent_sodium_current_h_gate"/>
         <component_ref component="persistent_sodium_current_m_gate"/>
      </component_ref>
      <component_ref component="delayed_rectifier_current">
         <component_ref component="delayed_rectifier_current_n_gate"/>
      </component_ref>
      <component_ref component="LVA_calcium_current">
         <component_ref component="LVA_calcium_current_h_gate"/>
         <component_ref component="LVA_calcium_current_m_gate"/>
      </component_ref>
      <component_ref component="N_HVA_calcium_current">
         <component_ref component="N_HVA_calcium_current_h_gate"/>
         <component_ref component="N_HVA_calcium_current_m_gate"/>
      </component_ref>
      <component_ref component="P_HVA_calcium_current">
         <component_ref component="P_HVA_calcium_current_m_gate"/>
      </component_ref>
      <component_ref component="calcium_dependent_potassium_current">
         <component_ref component="calcium_dependent_potassium_current_z_gate"/>
      </component_ref>
      <component_ref component="fast_transient_potassium_current">
         <component_ref component="fast_transient_potassium_current_h_gate"/>
         <component_ref component="fast_transient_potassium_current_m_gate"/>
      </component_ref>
      <component_ref component="hyperpolarization_activated_current">
         <component_ref component="hyperpolarization_activated_current_m_gate"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_2="membrane" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="persistent_sodium_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="delayed_rectifier_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="leak_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="LVA_calcium_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="N_HVA_calcium_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="P_HVA_calcium_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dependent_potassium_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="fast_transient_potassium_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="hyperpolarization_activated_current" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="stimulus_protocol" component_1="environment"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="stimulus_protocol" component_1="membrane"/>
      <map_variables variable_2="i_app" variable_1="i_app"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current_h_gate" component_1="sodium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="h" variable_1="h"/>
   </connection>
   <connection>
      <map_components component_2="sodium_current_m_gate" component_1="sodium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="m" variable_1="m"/>
   </connection>
   <connection>
      <map_components component_2="persistent_sodium_current" component_1="sodium_current"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="persistent_sodium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_NaP" variable_1="i_NaP"/>
   </connection>
   <connection>
      <map_components component_2="persistent_sodium_current_h_gate" component_1="persistent_sodium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="h" variable_1="h"/>
   </connection>
   <connection>
      <map_components component_2="persistent_sodium_current_m_gate" component_1="persistent_sodium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="m" variable_1="m"/>
   </connection>
   <connection>
      <map_components component_2="delayed_rectifier_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K" variable_1="i_K"/>
   </connection>
   <connection>
      <map_components component_2="delayed_rectifier_current_n_gate" component_1="delayed_rectifier_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="n" variable_1="n"/>
   </connection>
   <connection>
      <map_components component_2="leak_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_leak" variable_1="i_leak"/>
   </connection>
   <connection>
      <map_components component_2="LVA_calcium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_T" variable_1="i_T"/>
   </connection>
   <connection>
      <map_components component_2="LVA_calcium_current_m_gate" component_1="LVA_calcium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="m" variable_1="m"/>
   </connection>
   <connection>
      <map_components component_2="LVA_calcium_current_h_gate" component_1="LVA_calcium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="h" variable_1="h"/>
   </connection>
   <connection>
      <map_components component_2="N_HVA_calcium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_N" variable_1="i_N"/>
   </connection>
   <connection>
      <map_components component_2="N_HVA_calcium_current_m_gate" component_1="N_HVA_calcium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="m" variable_1="m"/>
   </connection>
   <connection>
      <map_components component_2="N_HVA_calcium_current_h_gate" component_1="N_HVA_calcium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="h" variable_1="h"/>
   </connection>
   <connection>
      <map_components component_2="P_HVA_calcium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_P" variable_1="i_P"/>
   </connection>
   <connection>
      <map_components component_2="P_HVA_calcium_current_m_gate" component_1="P_HVA_calcium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="m" variable_1="m"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dependent_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_SK" variable_1="i_SK"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dependent_potassium_current_z_gate" component_1="calcium_dependent_potassium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="z" variable_1="z"/>
      <map_variables variable_2="i_T" variable_1="i_T"/>
      <map_variables variable_2="i_N" variable_1="i_N"/>
      <map_variables variable_2="i_P" variable_1="i_P"/>
   </connection>
   <connection>
      <map_components component_2="LVA_calcium_current" component_1="calcium_dependent_potassium_current"/>
      <map_variables variable_2="i_T" variable_1="i_T"/>
   </connection>
   <connection>
      <map_components component_2="N_HVA_calcium_current" component_1="calcium_dependent_potassium_current"/>
      <map_variables variable_2="i_N" variable_1="i_N"/>
   </connection>
   <connection>
      <map_components component_2="P_HVA_calcium_current" component_1="calcium_dependent_potassium_current"/>
      <map_variables variable_2="i_P" variable_1="i_P"/>
   </connection>
   <connection>
      <map_components component_2="fast_transient_potassium_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_A" variable_1="i_A"/>
   </connection>
   <connection>
      <map_components component_2="fast_transient_potassium_current_m_gate" component_1="fast_transient_potassium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="m" variable_1="m"/>
   </connection>
   <connection>
      <map_components component_2="fast_transient_potassium_current_h_gate" component_1="fast_transient_potassium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="h" variable_1="h"/>
   </connection>
   <connection>
      <map_components component_2="hyperpolarization_activated_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_H" variable_1="i_H"/>
   </connection>
   <connection>
      <map_components component_2="hyperpolarization_activated_current_m_gate" component_1="hyperpolarization_activated_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="m" variable_1="m"/>
   </connection>
   <connection>
      <map_components component_2="N_HVA_calcium_current" component_1="LVA_calcium_current"/>
      <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
   </connection>
   <connection>
      <map_components component_2="P_HVA_calcium_current" component_1="LVA_calcium_current"/>
      <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
   </connection>
   <connection>
      <map_components component_2="calcium_dependent_potassium_current" component_1="delayed_rectifier_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="fast_transient_potassium_current" component_1="delayed_rectifier_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#8b1978f5-4405-4f4d-b904-af40ceb0fd4a">
    <rdf:li rdf:resource="rdf:#9e6882d7-6d79-4ea2-ace3-e84080d5ae48"/>
    <rdf:li rdf:resource="rdf:#a6e11a9d-5e0f-476c-a35e-99b777b1dfbd"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#206126f2-1c75-4de3-8dcf-ae7ed109c156">
    <dcterms:W3CDTF>2008-07-17T14:14:40+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#48d473ab-bf4e-4724-b3f8-b31a5aa3d5d7">
    <vCard:FN>Geoff Nunns</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#30a31d97-f602-421e-aa27-18fed465f455">
    <vCard:Orgname>Auckland Bioengineering Institute</vCard:Orgname>
    <vCard:Orgunit>CellML Team</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bc043ac6-6685-4510-aa65-3d0853b0ab92">
    <vCard:ORG rdf:resource="rdf:#30a31d97-f602-421e-aa27-18fed465f455"/>
    <vCard:EMAIL rdf:resource="rdf:#29216841-4983-4b2c-b509-ad23be935a99"/>
    <vCard:N rdf:resource="rdf:#efb8aefa-5853-4e5a-ada0-4acca2b658ef"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#87b021c8-0019-4b85-990d-0a7f8316944d">
    <dcterms:modified rdf:resource="rdf:#206126f2-1c75-4de3-8dcf-ae7ed109c156"/>
    <rdf:value>Fixed "not exactly equivalent but dimensionally equivalent" errors.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#33e024e1-01fe-4630-834d-c4d672c9da78"/>
  </rdf:Description>
  <rdf:Description rdf:about="#butera_purvis_2004">
    <bqs:reference rdf:resource="rdf:#0374b293-bd96-4711-8837-863da4cb0825"/>
    <cmeta:comment rdf:resource="rdf:#dc951c89-3e4b-4fc8-82bd-564c5da002c4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dc951c89-3e4b-4fc8-82bd-564c5da002c4">
    <dc:creator rdf:resource="rdf:#5efe7532-c46c-4924-94ad-0755ae564c0c"/>
    <rdf:value>This CellML model is known to run in both PCEnv and COR and the units have been checked and are consistent. However the CellML model is curently unable to perfectly recreate the published results (we are several orders of magnitude out with certain variables). Further curation is required and we are endevouring to seek the help of the original model authors.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a6e11a9d-5e0f-476c-a35e-99b777b1dfbd">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#278539e8-acab-484c-8141-07b6b4c651cd"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8ae42df4-ef58-465e-9e64-502c1573bce9">
    <dc:creator rdf:resource="rdf:#48d473ab-bf4e-4724-b3f8-b31a5aa3d5d7"/>
    <rdf:value>This CellML model is known to run in both PCEnv and COR and the units have been checked and are consistent. However the CellML model is curently unable to perfectly recreate the published results (we are several orders of magnitude out with certain variables). Further curation is required and we are endevouring to seek the help of the original model authors.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0374b293-bd96-4711-8837-863da4cb0825">
    <bqs:Pubmed_id>15653786</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#9f5236bb-358f-46e3-978d-07af149fe984"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher/>
    <cmeta:comment rdf:resource="rdf:#8ae42df4-ef58-465e-9e64-502c1573bce9"/>
    <dcterms:created rdf:resource="rdf:#02040c9c-1efc-4ffc-ae82-5e2ae8bfb463"/>
    <dc:creator rdf:resource="rdf:#bc043ac6-6685-4510-aa65-3d0853b0ab92"/>
    <cmeta:modification rdf:resource="rdf:#87b021c8-0019-4b85-990d-0a7f8316944d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#efb8aefa-5853-4e5a-ada0-4acca2b658ef">
    <vCard:Given>Geoffrey</vCard:Given>
    <vCard:Family>Nunns</vCard:Family>
    <vCard:Other>Rogan</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#33e024e1-01fe-4630-834d-c4d672c9da78">
    <vCard:N rdf:resource="rdf:#a6ce75b8-d758-49c6-9a0c-8829d3907745"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2eaddb83-3aa5-47d9-9a66-1af32773010e">
    <dcterms:W3CDTF>2005-00-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#209e456f-249e-4c27-9414-a0509255fd87">
    <dc:title>Journal of Neurophysiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#02040c9c-1efc-4ffc-ae82-5e2ae8bfb463">
    <dcterms:W3CDTF>2008-07-17T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9e6882d7-6d79-4ea2-ace3-e84080d5ae48">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#3f3159e0-93cc-464a-99cf-eea1d3f5b1d8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#29216841-4983-4b2c-b509-ad23be935a99">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>gnunns1@jhem.jhu.edu</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#278539e8-acab-484c-8141-07b6b4c651cd">
    <vCard:Given>Robert</vCard:Given>
    <vCard:Family>Butera</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9f5236bb-358f-46e3-978d-07af149fe984">
    <dc:creator rdf:resource="rdf:#8b1978f5-4405-4f4d-b904-af40ceb0fd4a"/>
    <dc:title>Ionic Current Model of a Hypoglossal Motoneuron</dc:title>
    <bqs:volume>93</bqs:volume>
    <bqs:first_page>723</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#209e456f-249e-4c27-9414-a0509255fd87"/>
    <dcterms:issued rdf:resource="rdf:#2eaddb83-3aa5-47d9-9a66-1af32773010e"/>
    <bqs:last_page>733</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a6ce75b8-d758-49c6-9a0c-8829d3907745">
    <vCard:Given>Geoffrey</vCard:Given>
    <vCard:Family>Nunns</vCard:Family>
    <vCard:Other>Rogan</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3f3159e0-93cc-464a-99cf-eea1d3f5b1d8">
    <vCard:Given>Liston</vCard:Given>
    <vCard:Family>Purvis</vCard:Family>
    <vCard:Other>K</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5efe7532-c46c-4924-94ad-0755ae564c0c">
    <vCard:FN>Geoff Nunns</vCard:FN>
  </rdf:Description>
</rdf:RDF>
</model>