Location: Butera, Rinzel, Smith, 1999 @ 7d39300d356a / butera_rinzel_smith_1999.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2008-05-30 02:51:16+12:00
Desc:
committing version04 of butera_rinzel_smith_1999
Permanent Source URI:
https://staging.physiomeproject.org/workspace/butera_rinzel_smith_1999/rawfile/7d39300d356ac290c1447217b33db3f591770bab/butera_rinzel_smith_1999.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 28/05/2008 at 1:46:31 at p.m. using:

COR (0.9.31.955)
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_rinzel_smith_1999_version02" name="butera_rinzel_smith_1999_version02">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Models Of Respiratory Rhythm Generation In The Pre-Botzinger Complex. I. Bursting Pacemaker Neurons</title>
				<author>
					<firstname>Catherine</firstname>
					<surname>Lloyd</surname>
					<affiliation>
						<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
            This published paper contains two minimal models to describe the behaviour of pacemaker neurons. This CellML model (version 4) represents model 2 from the paper: where bursting behaviour bursting arises via a fast-activating persistent Na+ current (INaP) and the slow activation of a K+ current (IKS). In both models, action potentials are generated via fast Na+ and K+ currents. The two models also differ in a few parameters to facilitate a rigorous comparison of the two different burst-generating mechanisms.
          </para>
          <para>
            This model is known to run in both PCEnv and COR to recreate the published results.  Note that some initial conditions were added which were required but were not in the published paper.  To the best of our knowledge these appear to be correct.
          </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
Inhalation and exhalation movements associated with respiration in mammals are generated by networks of neurons in the lower brain stem that produce a rhythmic pattern of electrical activity.  The principal neuronal kernel for rhythm generation has been located in the pre-Botzinger complex, a subregion of the ventro-lateral medulla.  In order to understand respiratory rhythm generation, it is necessary to consider mechanisms incorporating intrinsic cellular pacemaker properties.  These mechanisms are captured by a hybrid pacemaker model (Smith 1997; Smith <emphasis>et al</emphasis>. 1995), in which a rhythm arises from the dynamic interactions of both intrinsic and synaptic properties within a bilaterally distributed population of coupled bursting pacemaker neurons.  ("Bursting" refers to a complicated pattern of electrical activity.  Bursts of action potential spikes (the "active" phase) are observed, separated by a "silent" phase of membrane repolarisation).
</para>
				<para>  
In their 1999 paper, Robert J. Butera, J<subscript>R</subscript>., John Rinzel and Jeffrey C. Smith present two computational versions of this hybrid pacemaker-network model.  In the first model, bursting arises via fast activation and slow inactivation of a persistent Na<superscript>+</superscript> current, I<subscript>NaP</subscript>.  In the second model, bursting arises via a fast-activating persistent Na<superscript>+</superscript> current, I<subscript>NaP</subscript>,(the inactivation term "h" has been removed) and slow activation of a K<superscript>+</superscript> current, I<subscript>KS</subscript>.  In both models, action potentials are generated via fast Na<superscript>+</superscript> and K<superscript>+</superscript> currents.  Both models are consistent with experimental data, and the authors suggest several experimental tests to demonstrate the validity of either model and to differentiate between the two mechanisms.
</para>
				<para>
The complete original paper reference is cited below:
</para>
				<para>
					<ulink url="http://jn.physiology.org/cgi/content/abstract/82/1/382">Models of Respiratory Rhythm Generation in the Pre-Bötzinger Complex. I. Bursting Pacemaker Neurons</ulink>, Robert J. Butera, Jr., John Rinzel and Jeffrey C. Smith, 1999, <ulink url="http://jn.physiology.org/">
						<emphasis>Journal of Neurophysiology</emphasis>
					</ulink>, 81, 382-397.  (<ulink url="http://jn.physiology.org/cgi/content/full/82/1/382">Full text</ulink> and <ulink url="http://jn.physiology.org/cgi/reprint/82/1/382.pdf">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=10400966&amp;dopt=Abstract">PubMed ID: 10400966</ulink>
				</para>
				<informalfigure float="0" id="fig_cell_diagram2">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>diagram of the first model</title>
							</objectinfo>
							<imagedata fileref="butera_1999b.png"/>
						</imageobject>
					</mediaobject>
					<caption>The second mathematical model is based on a single-compartment Hodgkin-Huxley type formalism.  It is composed of six ionic currents across the plasma membrane: a fast sodium current, I<subscript>Na</subscript>; a delayed rectifier potassium current, I<subscript>K</subscript>; a slow K<superscript>+</superscript> current, I<subscript>KS</subscript>; a persistent sodium current,  I<subscript>NaP</subscript>; a passive leakage current, I<subscript>L</subscript>; and a tonic current, I<subscript>tonic_e</subscript> (although this last current is considered to be inactive in these models).  Note that the removal of the inactivation term "h" from I<subscript>NaP</subscript> is not visible in the model diagram.</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="picoA">
      <unit units="ampere" prefix="nano"/>
   </units>
   <units name="nanoS">
      <unit units="siemens" prefix="nano"/>
   </units>
   <units name="picoF">
      <unit units="farad" prefix="pico"/>
   </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="1"/>
      <variable units="picoF" name="C" initial_value="21"/>
      <variable units="picoA" name="i_app" initial_value="0"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="picoA" public_interface="in" name="i_NaP"/>
      <variable units="picoA" public_interface="in" name="i_Na"/>
      <variable units="picoA" public_interface="in" name="i_K"/>
      <variable units="picoA" public_interface="in" name="i_KS"/>
      <variable units="picoA" public_interface="in" name="i_L"/>
      <variable units="picoA" public_interface="in" name="i_tonic_e"/>
      <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_NaP</ci>
                        <ci>i_Na</ci>
                        <ci>i_K</ci>
                        <ci>i_KS</ci>
                        <ci>i_L</ci>
                        <ci>i_tonic_e</ci>
                     </apply>
                  </apply>
                  <ci>i_app</ci>
               </apply>
               <ci>C</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current">
      <variable units="picoA" public_interface="out" name="i_Na"/>
      <variable units="millivolt" public_interface="out" name="E_Na" initial_value="50"/>
      <variable units="nanoS" name="g_Na" initial_value="28"/>
      <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_infinity"/>
      <variable units="dimensionless" private_interface="in" name="n"/>
      <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_infinity</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>n</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_m_gate">
      <variable units="dimensionless" public_interface="out" name="m_infinity"/>
      <variable units="millivolt" name="theta_m" initial_value="-34"/>
      <variable units="millivolt" name="omega_m" initial_value="-5"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <ci>theta_m</ci>
                        </apply>
                        <ci>omega_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_sodium_current_n_gate">
      <variable units="dimensionless" public_interface="out" name="n" initial_value="0.003"/>
      <variable units="dimensionless" name="n_infinity"/>
      <variable units="millisecond" name="tau_n"/>
      <variable units="millisecond" name="tau_n_max" initial_value="10"/>
      <variable units="millivolt" name="theta_n" initial_value="-29"/>
      <variable units="millivolt" name="omega_n" initial_value="-4"/>
      <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>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>n</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>n_infinity</ci>
                  <ci>n</ci>
               </apply>
               <ci>tau_n</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>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/>
                           <ci>V</ci>
                           <ci>theta_n</ci>
                        </apply>
                        <ci>omega_n</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_n</ci>
            <apply>
               <divide/>
               <ci>tau_n_max</ci>
               <apply>
                  <cosh/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <ci>theta_n</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>omega_n</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="potassium_current">
      <variable units="picoA" public_interface="out" name="i_K"/>
      <variable units="millivolt" public_interface="out" name="E_K" initial_value="-85"/>
      <variable units="nanoS" name="g_K" initial_value="11.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="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="potassium_current_n_gate">
      <variable units="dimensionless" public_interface="out" name="n" initial_value="0.003"/>
      <variable units="dimensionless" name="n_infinity"/>
      <variable units="millisecond" name="tau_n"/>
      <variable units="millisecond" name="tau_n_max" initial_value="10"/>
      <variable units="millivolt" name="theta_n" initial_value="-29"/>
      <variable units="millivolt" name="omega_n" initial_value="-4"/>
      <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>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>n</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>n_infinity</ci>
                  <ci>n</ci>
               </apply>
               <ci>tau_n</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>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/>
                           <ci>V</ci>
                           <ci>theta_n</ci>
                        </apply>
                        <ci>omega_n</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_n</ci>
            <apply>
               <divide/>
               <ci>tau_n_max</ci>
               <apply>
                  <cosh/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <ci>theta_n</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>omega_n</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_potassium_current">
      <variable units="picoA" public_interface="out" name="i_KS"/>
      <variable units="nanoS" name="g_KS" initial_value="5.6"/>
      <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="k"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_KS</ci>
            <apply>
               <times/>
               <ci>g_KS</ci>
               <ci>k</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_potassium_current_k_gate">
      <variable units="dimensionless" public_interface="out" name="k" initial_value="0.01"/>
      <variable units="dimensionless" name="k_infinity"/>
      <variable units="millisecond" name="tau_k"/>
      <variable units="millisecond" name="tau_k_max" initial_value="10000"/>
      <variable units="millivolt" name="theta_k" initial_value="-38"/>
      <variable units="millivolt" name="omega_k" initial_value="-6"/>
      <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>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>k</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>k_infinity</ci>
                  <ci>k</ci>
               </apply>
               <ci>tau_k</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <ci>theta_k</ci>
                        </apply>
                        <ci>omega_k</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_k</ci>
            <apply>
               <divide/>
               <ci>tau_k_max</ci>
               <apply>
                  <cosh/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <ci>theta_k</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>omega_k</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="persistent_sodium_current">
      <variable units="picoA" public_interface="out" name="i_NaP"/>
      <variable units="nanoS" name="g_NaP" initial_value="2.8"/>
      <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_Na"/>
      <variable units="dimensionless" private_interface="in" name="m_infinity"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_NaP</ci>
            <apply>
               <times/>
               <ci>g_NaP</ci>
               <ci>m_infinity</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Na</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="persistent_sodium_current_m_gate">
      <variable units="dimensionless" public_interface="out" name="m_infinity"/>
      <variable units="millivolt" name="theta_m" initial_value="-40"/>
      <variable units="millivolt" name="omega_m" initial_value="-6"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <ci>theta_m</ci>
                        </apply>
                        <ci>omega_m</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="leakage_current">
      <variable units="picoA" public_interface="out" name="i_L"/>
      <variable units="nanoS" name="g_L" initial_value="2.8"/>
      <variable units="millivolt" name="E_L" initial_value="-65"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <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="tonic_current">
      <variable units="picoA" public_interface="out" name="i_tonic_e"/>
      <variable units="nanoS" name="g_tonic_e" initial_value="0"/>
      <variable units="millivolt" name="E_syn_e" initial_value="0"/>
      <variable units="millisecond" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_tonic_e</ci>
            <apply>
               <times/>
               <ci>g_tonic_e</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_syn_e</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="fast_sodium_current">
            <component_ref component="fast_sodium_current_m_gate"/>
            <component_ref component="fast_sodium_current_n_gate"/>
         </component_ref>
         <component_ref component="potassium_current">
            <component_ref component="potassium_current_n_gate"/>
         </component_ref>
         <component_ref component="slow_potassium_current">
            <component_ref component="slow_potassium_current_k_gate"/>
         </component_ref>
         <component_ref component="persistent_sodium_current">
            <component_ref component="persistent_sodium_current_m_gate"/>
         </component_ref>
         <component_ref component="leakage_current"/>
         <component_ref component="tonic_current"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="fast_sodium_current">
         <component_ref component="fast_sodium_current_m_gate"/>
         <component_ref component="fast_sodium_current_n_gate"/>
      </component_ref>
      <component_ref component="potassium_current">
         <component_ref component="potassium_current_n_gate"/>
      </component_ref>
      <component_ref component="slow_potassium_current">
         <component_ref component="slow_potassium_current_k_gate"/>
      </component_ref>
      <component_ref component="persistent_sodium_current">
         <component_ref component="persistent_sodium_current_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="fast_sodium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="slow_potassium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="persistent_sodium_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="leakage_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="tonic_current"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="fast_sodium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_Na" variable_1="i_Na"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="potassium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K" variable_1="i_K"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="slow_potassium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_KS" variable_1="i_KS"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="persistent_sodium_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_NaP" variable_1="i_NaP"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="leakage_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_L" variable_1="i_L"/>
   </connection>
   <connection>
      <map_components component_2="membrane" component_1="tonic_current"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_tonic_e" variable_1="i_tonic_e"/>
   </connection>
   <connection>
      <map_components component_2="persistent_sodium_current" component_1="fast_sodium_current"/>
      <map_variables variable_2="E_Na" variable_1="E_Na"/>
   </connection>
   <connection>
      <map_components component_2="slow_potassium_current" component_1="potassium_current"/>
      <map_variables variable_2="E_K" variable_1="E_K"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_m_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="m_infinity" variable_1="m_infinity"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="fast_sodium_current_n_gate" component_1="fast_sodium_current"/>
      <map_variables variable_2="n" variable_1="n"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="potassium_current_n_gate" component_1="potassium_current"/>
      <map_variables variable_2="n" variable_1="n"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="slow_potassium_current_k_gate" component_1="slow_potassium_current"/>
      <map_variables variable_2="k" variable_1="k"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="persistent_sodium_current_m_gate" component_1="persistent_sodium_current"/>
      <map_variables variable_2="m_infinity" variable_1="m_infinity"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#11d48b49-9c76-4579-8aad-15d94dade852">
    <rdf:li>electrophysiological</rdf:li>
    <rdf:li>neuron</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#citationAuthorsSeq">
    <rdf:li rdf:resource="rdf:#author1Vcard"/>
    <rdf:li rdf:resource="rdf:#author2Vcard"/>
    <rdf:li rdf:resource="rdf:#author3Vcard"/>
  </rdf:Seq>
  <rdf:Seq rdf:about="rdf:#20352c5f-770c-4ebc-942a-c3937011566a">
    <rdf:li rdf:resource="rdf:#9306fd13-cf2f-4afc-807f-e487279f4796"/>
    <rdf:li rdf:resource="rdf:#8a9e7e18-32e5-490e-b562-8b78943c9466"/>
    <rdf:li rdf:resource="rdf:#fad4a924-670c-4ffc-9c17-b8bf51a46a9f"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#b013bf39-3e95-4d9b-89b7-4dae8879b257">
    <dcterms:W3CDTF>2003-03-28</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#68c80246-e345-4f02-beb9-baafe99ab53c"/>
    <dcterms:created rdf:resource="rdf:#a802ea57-edc3-467b-97df-0ff66207e3a2"/>
    <dc:creator rdf:resource="rdf:#c904c27c-6950-4cd9-b797-a734408cbce2"/>
    <cmeta:modification rdf:resource="rdf:#0ecd6368-9e19-4064-92c3-973f604c17ea"/>
    <cmeta:modification rdf:resource="rdf:#63cccde8-5cc0-4549-9f4c-53264d73f8d1"/>
    <cmeta:modification rdf:resource="rdf:#82121e35-eddf-4d90-a784-6ccb0ec22428"/>
    <cmeta:modification rdf:resource="rdf:#b3e01c8f-980c-4b5d-871f-c9cd84e0e317"/>
    <cmeta:modification rdf:resource="rdf:#bb71921f-1a5f-458c-97bd-e960d45db3c5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#33b10fd4-7264-40ee-b07d-3a988d0e8b2a">
    <vCard:Given>Geoffrey</vCard:Given>
    <vCard:Family>Nunns</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#68fcf6a9-7520-4e92-98b9-ff2c8dd1a63e">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#227354f2-8bb6-49ba-8e70-a87262858157">
    <dcterms:W3CDTF>1999-07-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3fe3f7da-036d-447a-9af4-daefa3dc78c0">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9b1caf23-61f5-4e29-8e08-b38b506638d1">
    <dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author3VcardN">
    <vCard:Given>Jeffrey</vCard:Given>
    <vCard:Family>Smith</vCard:Family>
    <vCard:Other>C</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#72931518-64d1-47ca-b192-e44b05773d77">
    <dc:subject rdf:resource="rdf:#10276f05-a56d-456a-a816-412d7802dbc7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0a82f688-cd55-4e89-bb03-69e2f0175210">
    <vCard:N rdf:resource="rdf:#33b10fd4-7264-40ee-b07d-3a988d0e8b2a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#591cba3a-aba4-48a5-89bd-15ea34cf4cf1">
    <dcterms:W3CDTF>2002-07-18</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8a9e7e18-32e5-490e-b562-8b78943c9466">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#e34646b9-5e65-49fc-988a-9b762b1827d7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fad4a924-670c-4ffc-9c17-b8bf51a46a9f">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#62003e04-ccd5-44d5-8ccb-042b2b83f95d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#63cccde8-5cc0-4549-9f4c-53264d73f8d1">
    <dcterms:modified rdf:resource="rdf:#9946545a-fa42-4000-95eb-e3849f316cb2"/>
    <rdf:value>
          Made MathML id's unique
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#9292a948-e772-45fa-a751-3043045c2485"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#10276f05-a56d-456a-a816-412d7802dbc7">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#11d48b49-9c76-4579-8aad-15d94dade852"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ab1ea60d-8a86-445d-b8b5-563d3af3c96a">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9292a948-e772-45fa-a751-3043045c2485">
    <vCard:N rdf:resource="rdf:#9992748e-2193-4841-bfb4-aa9b62c599bc"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#56565406-7709-40e5-a055-b762c58ee618">
    <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:#0ecd6368-9e19-4064-92c3-973f604c17ea">
    <dcterms:modified rdf:resource="rdf:#e1494f2a-bd91-4131-9b36-350cdd2f8fea"/>
    <rdf:value>Defined initial values for variables defined by ODEs. While these did not appear in the paper, they accurately replicate published results. </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#0a82f688-cd55-4e89-bb03-69e2f0175210"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ba7565d7-495f-49d1-832d-06d63d54fe22">
    <bqs:Pubmed_id>10400966</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#cfcf968d-a332-431e-8654-6052180bc0d3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9946545a-fa42-4000-95eb-e3849f316cb2">
    <dcterms:W3CDTF>2005-04-20</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author2VcardN">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Rinzel</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c904c27c-6950-4cd9-b797-a734408cbce2">
    <vCard:ORG rdf:resource="rdf:#3fe3f7da-036d-447a-9af4-daefa3dc78c0"/>
    <vCard:EMAIL rdf:resource="rdf:#56565406-7709-40e5-a055-b762c58ee618"/>
    <vCard:N rdf:resource="rdf:#68fcf6a9-7520-4e92-98b9-ff2c8dd1a63e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b73d79cf-2907-4dca-b7b6-6b6e38c6fbd0">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3e63aaac-04dc-481a-9c91-05d52589e2cb">
    <dc:creator rdf:resource="rdf:#ab1ea60d-8a86-445d-b8b5-563d3af3c96a"/>
    <rdf:value>This is the CellML description of Butera et al's second 1999 mathematical model of respiratory rhythm generation in the pre-Botzinger complex in bursting pacemaker neurons.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#574cb709-4de2-4ede-994a-a5f8774f6f15">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e34646b9-5e65-49fc-988a-9b762b1827d7">
    <vCard:Given>John</vCard:Given>
    <vCard:Family>Rinzel</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#68c80246-e345-4f02-beb9-baafe99ab53c">
    <dc:creator rdf:resource="rdf:#b73d79cf-2907-4dca-b7b6-6b6e38c6fbd0"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author2Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author2VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e1494f2a-bd91-4131-9b36-350cdd2f8fea">
    <dcterms:W3CDTF>2008-05-30T14:46:41+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a802ea57-edc3-467b-97df-0ff66207e3a2">
    <dcterms:W3CDTF>2002-05-09T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#82121e35-eddf-4d90-a784-6ccb0ec22428">
    <dcterms:modified rdf:resource="rdf:#591cba3a-aba4-48a5-89bd-15ea34cf4cf1"/>
    <rdf:value>
          Added more metadata.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#0576f4d7-093f-44cb-b43b-4e2b2afe8658"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d29cfd60-3fcc-4268-bb51-9c0d889a2c99">
    <dc:title>Journal of Neurophysiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="#butera_rinzel_smith_1999_version02">
    <dc:title>
         Butera et al's second 1999 mathematical model of respiratory rhythm 
         generation in the pre-Botzinger complex in bursting pacemaker neurons.
      </dc:title>
    <cmeta:bio_entity>pacemaker neuron</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#3e63aaac-04dc-481a-9c91-05d52589e2cb"/>
    <bqs:reference rdf:resource="rdf:#72931518-64d1-47ca-b192-e44b05773d77"/>
    <bqs:reference rdf:resource="rdf:#ba7565d7-495f-49d1-832d-06d63d54fe22"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#62003e04-ccd5-44d5-8ccb-042b2b83f95d">
    <vCard:Given>Jeffrey</vCard:Given>
    <vCard:Family>Smith</vCard:Family>
    <vCard:Other>C</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b3e01c8f-980c-4b5d-871f-c9cd84e0e317">
    <dcterms:modified rdf:resource="rdf:#b013bf39-3e95-4d9b-89b7-4dae8879b257"/>
    <rdf:value>
          Changed the name and cmeta:id of the model so that they were not the 
          same as those of the other Butera 1999 model.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#9d0dbbac-44a2-4f4b-8f04-0e876e255cc9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9ad04996-44b1-4497-addf-b986765ad561">
    <vCard:Given>Robert</vCard:Given>
    <vCard:Family>Butera</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cfcf968d-a332-431e-8654-6052180bc0d3">
    <dc:creator rdf:resource="rdf:#20352c5f-770c-4ebc-942a-c3937011566a"/>
    <dc:title>Models of Respiratory Rhythm Generation in the Pre-Botzinger Complex. I. Bursting Pacemaker Neurons</dc:title>
    <bqs:volume>81</bqs:volume>
    <bqs:first_page>382</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#d29cfd60-3fcc-4268-bb51-9c0d889a2c99"/>
    <dcterms:issued rdf:resource="rdf:#227354f2-8bb6-49ba-8e70-a87262858157"/>
    <bqs:last_page>397</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author1Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author1VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author1VcardN">
    <vCard:Given>Robert</vCard:Given>
    <vCard:Family>Butera</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9992748e-2193-4841-bfb4-aa9b62c599bc">
    <vCard:Given>Peter</vCard:Given>
    <vCard:Family>Villiger</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#db6ca05b-a4f4-416d-96a5-ec5820f72ff4">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#96dd05c7-cdcd-44aa-8983-2bf596025cd9">
    <vCard:N rdf:resource="rdf:#36d41e4f-5232-4c7b-80cb-5949a004cef7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author3Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author3VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0576f4d7-093f-44cb-b43b-4e2b2afe8658">
    <vCard:N rdf:resource="rdf:#db6ca05b-a4f4-416d-96a5-ec5820f72ff4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9d0dbbac-44a2-4f4b-8f04-0e876e255cc9">
    <vCard:N rdf:resource="rdf:#574cb709-4de2-4ede-994a-a5f8774f6f15"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#36d41e4f-5232-4c7b-80cb-5949a004cef7">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9306fd13-cf2f-4afc-807f-e487279f4796">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#9ad04996-44b1-4497-addf-b986765ad561"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bb71921f-1a5f-458c-97bd-e960d45db3c5">
    <dcterms:modified rdf:resource="rdf:#9b1caf23-61f5-4e29-8e08-b38b506638d1"/>
    <rdf:value>
          Added publication date information.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#96dd05c7-cdcd-44aa-8983-2bf596025cd9"/>
  </rdf:Description>
</rdf:RDF>
</model>