<?xml version='1.0' encoding='utf-8'?>
<!-- FILE : butera_model_1999.xml
CREATED : 9th May 2002
LAST MODIFIED : 20th April 2005
AUTHOR : Catherine Lloyd
Bioengineering Institute
The University of Auckland
MODEL STATUS : This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/01/2002 CellML Metadata 1.0 Specification.
DESCRIPTION : This file contains a CellML description of Butera et al's first 1999 mathematical model of respiratory rhythm generation in the pre-Botzinger complex in bursting pacemaker neurons.
CHANGES:
18/07/2002 - CML - Added more metadata.
09/04/2003 - AAC - Added publication date information.
20/04/2005 - PJV - Made MathML id's unique
--><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#" name="butera_rinzel_smith_1999_version01" cmeta:id="butera_rinzel_smith_1999_version01">
<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 is the original unchecked version of the model imported from the previous
CellML model repository, 24-Jan-2006.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
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 (see <xref linkend="fig_cell_diagram1"/> and <xref linkend="fig_cell_diagram2"/> below). 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&db=PubMed&list_uids=10400966&dopt=Abstract">PubMed ID: 10400966</ulink>
</para>
<para>
The raw CellML description of the model can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>
</para>
<informalfigure float="0" id="fig_cell_diagram1">
<mediaobject>
<imageobject>
<objectinfo>
<title>diagram of the first model</title>
</objectinfo>
<imagedata fileref="butera_1999a.png"/>
</imageobject>
</mediaobject>
<caption>The first mathematical model is based on a single-compartment Hodgkin-Huxley type formalism. It is composed of five ionic currents across the plasma membrane: a fast sodium current, I<subscript>Na</subscript>; a delayed rectifier potassium current, I<subscript>K</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).</caption>
</informalfigure>
<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 model appears identical to the first except with the addition of a slow K<superscript>+</superscript> current, I<subscript>KS</subscript>. (The removal of the inactivation term "h" from I<subscript>NaP</subscript> is not visible in the model diagram.)</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<!--
Below, we define some additional units for association with variables and
constants within the model. The identifiers are fairly self-explanatory.
-->
<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>
<!--
The "environment" component is used to declare variables that are used by
all or most of the other components, in this case just "time".
-->
<component name="environment">
<variable units="millisecond" public_interface="out" name="time"/>
</component>
<component name="membrane">
<variable units="millivolt" public_interface="out" name="V"/>
<variable units="picoF" name="C" initial_value="21.0"/>
<variable units="picoA" name="i_app" initial_value="0.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_L"/>
<variable units="picoA" public_interface="in" name="i_tonic_e"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="membrane_voltage_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> V </ci>
</apply>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<ci> i_NaP </ci>
<ci> i_Na </ci>
<ci> i_K </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.0"/>
<variable units="nanoS" name="g_Na" initial_value="28.0"/>
<variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="dimensionless" private_interface="in" name="m_infinity"/>
<variable units="dimensionless" private_interface="in" name="n"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_Na_calculation">
<eq/>
<ci> i_Na </ci>
<apply>
<times/>
<ci> g_Na </ci>
<apply>
<power/>
<ci> m_infinity </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </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.0"/>
<variable units="millivolt" name="omega_m" initial_value="-5.0"/>
<variable units="millivolt" public_interface="in" name="V"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="fast_sodium_current_m_gate_m_infinity_calculation">
<eq/>
<ci> m_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </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"/>
<variable units="dimensionless" name="n_infinity"/>
<variable units="millisecond" name="tau_n"/>
<variable units="millisecond" name="tau_n_max" initial_value="10.0"/>
<variable units="millivolt" name="theta_n" initial_value="-29.0"/>
<variable units="millivolt" name="omega_n" initial_value="-4.0"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="fast_sodium_current_n_gate_n_diff_eq">
<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 id="fast_sodium_current_n_gate_n_infinity_calculation">
<eq/>
<ci> n_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </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 id="fast_sodium_current_n_gate_tau_n_calculation">
<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.0 </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="nanoS" name="g_K" initial_value="11.2"/>
<variable units="millivolt" name="E_K" initial_value="-85.0"/>
<variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="dimensionless" private_interface="in" name="n"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_K_calculation">
<eq/>
<ci> i_K </ci>
<apply>
<times/>
<ci> g_K </ci>
<apply>
<power/>
<ci> n </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="potassium_current_n_gate">
<variable units="dimensionless" public_interface="out" name="n"/>
<variable units="dimensionless" name="n_infinity"/>
<variable units="millisecond" name="tau_n"/>
<variable units="millisecond" name="tau_n_max" initial_value="10.0"/>
<variable units="millivolt" name="theta_n" initial_value="-29.0"/>
<variable units="millivolt" name="omega_n" initial_value="-4.0"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="potassium_current_n_gate_n_diff_eq">
<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 id="potassium_current_n_gate_n_infinity_calculation">
<eq/>
<ci> n_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </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 id="potassium_current_n_gate_tau_n_calculation">
<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.0 </cn>
<ci> omega_n </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"/>
<variable units="dimensionless" private_interface="in" name="h"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_NaP_calculation">
<eq/>
<ci> i_NaP </ci>
<apply>
<times/>
<ci> g_NaP </ci>
<ci> m_infinity </ci>
<ci> h </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.0"/>
<variable units="millivolt" name="omega_m" initial_value="-6.0"/>
<variable units="millivolt" public_interface="in" name="V"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="persistent_sodium_current_m_gate_m_infinity_calculation">
<eq/>
<ci> m_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </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="persistent_sodium_current_h_gate">
<variable units="dimensionless" public_interface="out" name="h"/>
<variable units="dimensionless" name="h_infinity"/>
<variable units="millisecond" name="tau_h"/>
<variable units="millisecond" name="tau_h_max" initial_value="10000.0"/>
<variable units="millivolt" name="theta_h" initial_value="-48.0"/>
<variable units="millivolt" name="omega_h" initial_value="6.0"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="h_diff_eq">
<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>
<apply id="h_infinity_calculation">
<eq/>
<ci> h_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
<ci> theta_h </ci>
</apply>
<ci> omega_h </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_h_calculation">
<eq/>
<ci> tau_h </ci>
<apply>
<divide/>
<ci> tau_h_max </ci>
<apply>
<cosh/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
<ci> theta_h </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> omega_h </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.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 id="i_L_calculation">
<eq/>
<ci> i_L </ci>
<apply>
<times/>
<ci> g_L </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_L </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="tonic_current">
<variable units="picoA" public_interface="out" name="i_tonic_e"/>
<variable units="nanoS" name="g_tonic_e" initial_value="0.0"/>
<variable units="millivolt" name="E_syn_e" initial_value="0.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/MathMtonic_e">
<apply id="i_tonic_e_calculation">
<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="persistent_sodium_current">
<component_ref component="persistent_sodium_current_m_gate"/>
<component_ref component="persistent_sodium_current_h_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="persistent_sodium_current">
<component_ref component="persistent_sodium_current_m_gate"/>
<component_ref component="persistent_sodium_current_h_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="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="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="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="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>
<connection>
<map_components component_2="persistent_sodium_current_h_gate" component_1="persistent_sodium_current"/>
<map_variables variable_2="h" variable_1="h"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<rdf:RDF>
<rdf:Seq rdf:about="rdf:#f7a413ca-e561-4d50-b5b8-2df2a585c2b0">
<rdf:li resource="rdf:#394263b5-c652-4f53-8636-eb84cecc2458"/>
<rdf:li resource="rdf:#d5f46795-de47-4352-a351-1f5750bb0003"/>
<rdf:li resource="rdf:#90029e84-c922-40c5-bcf2-90fd82f66f37"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#15bbd9c7-8c70-4dbc-b9f9-a53b9cdf5233">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#95fe2c41-d201-4f16-90a8-135918df1593">
<dcterms:W3CDTF>2002-05-09</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b48ca180-f74b-402d-bd52-694814dc2279">
<dcterms:W3CDTF>2005-04-20</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#578cd774-5f5d-4fc3-a95d-7e159a096db9">
<dcterms:modified rdf:resource="rdf:#4568d027-e411-4065-822b-21b56cc51897"/>
<rdf:value>
Added publication date information.
</rdf:value>
<cmeta:modifier rdf:resource="rdf:#3cce81f8-8e6f-4dcf-99b6-fd0c725ceab7"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3cce81f8-8e6f-4dcf-99b6-fd0c725ceab7">
<vCard:N rdf:resource="rdf:#a357e230-4a4a-4a3d-95ab-a7bb995e51a2"/>
</rdf:Description>
<rdf:Description rdf:about="#butera_rinzel_smith_1999_version01">
<dc:title>
Butera et al's first 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:#1e70c7fb-dfc7-4b5e-a971-6f06d6a1e7be"/>
<bqs:reference rdf:resource="rdf:#ea99c96a-d446-41b4-be01-4f5c06ad6b6a"/>
<bqs:reference rdf:resource="rdf:#f5f5cdf7-d366-4040-a6a3-4a5bf5fe3a99"/>
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>respiratory</rdf:li>
<rdf:li>electrophysiology</rdf:li>
<rdf:li>central pattern generator</rdf:li>
<rdf:li>neuron</rdf:li>
<rdf:li>electrophysiological</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ea99c96a-d446-41b4-be01-4f5c06ad6b6a">
<bqs:Pubmed_id>10400966</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#cdef9fd9-3619-4861-8a70-1e2a2e1ce0d9"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#da00117b-2ff5-4766-aaed-53dbe7a242ef">
<vCard:ORG rdf:resource="rdf:#15bbd9c7-8c70-4dbc-b9f9-a53b9cdf5233"/>
<vCard:EMAIL rdf:resource="rdf:#9958a102-904d-49f7-ae3b-7a5132bfd16c"/>
<vCard:N rdf:resource="rdf:#12a5f78f-dd88-42ec-a8af-02e09587b438"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2ba081e7-4b2c-440a-a3a3-16e3d61ce5d1">
<vCard:N rdf:resource="rdf:#882a86c6-f385-49f9-baac-dd381486a301"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7fe065c8-f59e-4195-a3e8-561e5bedb167">
<vCard:N rdf:resource="rdf:#73063f81-eaf0-4e22-bb41-28823a3946f9"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9958a102-904d-49f7-ae3b-7a5132bfd16c">
<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:#c4f3c31c-4878-472a-b554-a68aa1a439dc">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1e70c7fb-dfc7-4b5e-a971-6f06d6a1e7be">
<dc:creator rdf:resource="rdf:#c4f3c31c-4878-472a-b554-a68aa1a439dc"/>
<rdf:value>
This is the CellML description of Butera et al's first 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:#394263b5-c652-4f53-8636-eb84cecc2458">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#90045915-0fef-4911-8140-35b19f3c6111"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#461e60b1-7713-42a6-806e-971b0088a97a">
<dcterms:W3CDTF>1999-07</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0e7ec12a-a6c0-4b09-9d75-77e858889dea">
<dcterms:modified rdf:resource="rdf:#68dcbed4-b31e-4069-b2f5-eea01129d31d"/>
<rdf:value>
Added more metadata.
</rdf:value>
<cmeta:modifier rdf:resource="rdf:#7fe065c8-f59e-4195-a3e8-561e5bedb167"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#12a5f78f-dd88-42ec-a8af-02e09587b438">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>
The University of Auckland, Bioengineering Institute
</dc:publisher>
<cmeta:modification rdf:resource="rdf:#0e7ec12a-a6c0-4b09-9d75-77e858889dea"/>
<cmeta:modification rdf:resource="rdf:#578cd774-5f5d-4fc3-a95d-7e159a096db9"/>
<cmeta:modification rdf:resource="rdf:#e73d9d13-3c1c-43a2-a142-2a9a9629802f"/>
<dcterms:created rdf:resource="rdf:#95fe2c41-d201-4f16-90a8-135918df1593"/>
<dc:creator rdf:resource="rdf:#da00117b-2ff5-4766-aaed-53dbe7a242ef"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#882a86c6-f385-49f9-baac-dd381486a301">
<vCard:Given>Peter</vCard:Given>
<vCard:Family>Villiger</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#90029e84-c922-40c5-bcf2-90fd82f66f37">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#9c6f6749-28ca-4a09-96da-b5b7d09671f5"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4568d027-e411-4065-822b-21b56cc51897">
<dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f5f5cdf7-d366-4040-a6a3-4a5bf5fe3a99">
<dc:subject rdf:resource="rdf:#b9ec4f7a-6176-4a1d-a932-9b5885b48afe"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a357e230-4a4a-4a3d-95ab-a7bb995e51a2">
<vCard:Given>Autumn</vCard:Given>
<vCard:Family>Cuellar</vCard:Family>
<vCard:Other>A</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9c6f6749-28ca-4a09-96da-b5b7d09671f5">
<vCard:Given>Jeffrey</vCard:Given>
<vCard:Family>Smith</vCard:Family>
<vCard:Other>C</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#44843408-f019-4caf-b365-63fa03b4d733">
<dc:title>Journal of Neurophysiology</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#73063f81-eaf0-4e22-bb41-28823a3946f9">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#90045915-0fef-4911-8140-35b19f3c6111">
<vCard:Given>Robert</vCard:Given>
<vCard:Family>Butera</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#cdef9fd9-3619-4861-8a70-1e2a2e1ce0d9">
<dc:creator rdf:resource="rdf:#f7a413ca-e561-4d50-b5b8-2df2a585c2b0"/>
<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:#44843408-f019-4caf-b365-63fa03b4d733"/>
<dcterms:issued rdf:resource="rdf:#461e60b1-7713-42a6-806e-971b0088a97a"/>
<bqs:last_page>397</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#68dcbed4-b31e-4069-b2f5-eea01129d31d">
<dcterms:W3CDTF>2002-07-18</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b9ec4f7a-6176-4a1d-a932-9b5885b48afe">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>calcium dymanics</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d9f742c4-5b72-4496-bbea-e0791a58c001">
<vCard:Given>John</vCard:Given>
<vCard:Family>Rinzel</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d5f46795-de47-4352-a351-1f5750bb0003">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#d9f742c4-5b72-4496-bbea-e0791a58c001"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e73d9d13-3c1c-43a2-a142-2a9a9629802f">
<dcterms:modified rdf:resource="rdf:#b48ca180-f74b-402d-bd52-694814dc2279"/>
<rdf:value>
Made MathML id's unique
</rdf:value>
<cmeta:modifier rdf:resource="rdf:#2ba081e7-4b2c-440a-a3a3-16e3d61ce5d1"/>
</rdf:Description>
</rdf:RDF>
</model>