- Author:
- Hanne <Hanne@hanne-nielsens-macbook.local>
- Date:
- 2010-07-13 13:49:48+12:00
- Desc:
- Added xul and session file
- Permanent Source URI:
- http://models.cellml.org/workspace/purvis_smith_koizumi_butera_2007/rawfile/d75208ca718c8b58c3628de76570d65b18a312e0/purvis_smith_koizumi_butera_2007_npc.cellml
<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 17/07/2008 at 1:23:14 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="purvis_2007" name="purvis_2007">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Intrinsic Bursters Increase the Robustness of Rythm Generation in an Excitatory Network</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This CellML model runs in OpenCell and COR to reproduce the published results (Figure 2A where g_tonic = 0.2 nS and g_NaP = 1.5 nS, note that E_L = -57.5 mv in accordance with the original Butera 1999 model). The CellML model should be run for a duration of 20000 ms with a step size of 0.1 ms and a high point density of 100000 points/graph. This model represents a single non-pacemaker cell.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
Abstract: The pre-Botzinger complex (pBC) is a vital subcircuit of the respiratory central pattern generator. Although the existence of neurons with pacemaker-like bursting properties in this network is not questioned, their role in network rhythmogenesis is unresolved. Modeling is ideally suited to address this debate because of the ease with which biophysical parameters of individual cells and network architecture can be manipulated. We modeled the parameter variability of experimental data from pBC bursting pacemaker and nonpacemaker neurons using a modified version of our previously developed pBC neuron and network models. To investigate the role of pacemakers in networkwide rhythmogenesis, we simulated networks of these neurons and varied the fraction of the population made up of pacemakers. For each number of pacemaker neurons, we varied the amount of tonic drive to the network and measured the frequency of synchronous networkwide bursting produced. Both excitatory networks with all-to-all coupling and sparsely connected networks were explored for several levels of synaptic coupling strength. Networks containing only nonpacemakers were able to produce networkwide bursting, but with a low probability of bursting and low input and output ranges. The results indicate that inclusion of pacemakers in an excitatory network increases robustness of the network by more than tripling the input and output ranges compared with networks containing no pacemakers. The largest increase in dynamic range occurs when the number of pacemakers in the network is greater than 20% of the population. Experimental tests of the model predictions are proposed.</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="purvis_2007.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram of the cell model.</caption>
</informalfigure>
<para>
The original paper reference is cited below:
</para>
<para>
Intrinsic Bursters Increase the Robustness of Rythm Generation in an Excitatory Network, L.K. Purvis, J.C. Smith, H. Koizumi, R.J. Butera 2007, <emphasis>Journal of Neurophysiology</emphasis>, 97, 1515-1526. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=17167061&dopt=Abstract">PubMed ID: 17167061 </ulink>
</para>
</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="nanoS">
<unit units="siemens" prefix="nano"/>
</units>
<units name="nanoF">
<unit units="farad" prefix="nano"/>
</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="-50.0"/>
<variable units="nanoF" name="C" initial_value="0.021"/>
<variable units="nanoA" name="i_app" initial_value="0"/>
<variable units="nanoA" public_interface="in" name="i_tonic_e"/>
<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="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>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>
</apply>
</apply>
<ci>i_tonic_e</ci>
<ci>i_app</ci>
</apply>
<ci>C</ci>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current">
<variable units="nanoA" 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" public_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>
<divide/>
<cn cellml:units="volt">1</cn>
<cn cellml:units="millivolt">1000</cn>
</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="potassium_current">
<variable units="nanoA" 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"/>
<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" public_interface="out" 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>
<divide/>
<cn cellml:units="volt">1</cn>
<cn cellml:units="millivolt">1000</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.01"/>
<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="persistent_sodium_current">
<variable units="nanoA" public_interface="out" name="i_NaP"/>
<variable units="nanoS" name="g_NaP" initial_value="1.5"/>
<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>
<eq/>
<ci>i_NaP</ci>
<apply>
<times/>
<ci>g_NaP</ci>
<ci>m_infinity</ci>
<ci>h</ci>
<apply>
<divide/>
<cn cellml:units="volt">1</cn>
<cn cellml:units="millivolt">1000</cn>
</apply>
<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="-45.1"/>
<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="persistent_sodium_current_h_gate">
<variable units="dimensionless" public_interface="out" name="h" initial_value="0.92"/>
<variable units="dimensionless" name="h_infinity"/>
<variable units="millisecond" name="tau_h"/>
<variable units="millisecond" name="tau_h_max" initial_value="10000"/>
<variable units="millivolt" name="theta_h" initial_value="-53"/>
<variable units="millivolt" name="omega_h" 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>h</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>h_infinity</ci>
<ci>h</ci>
</apply>
<ci>tau_h</ci>
</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>
<minus/>
<ci>V</ci>
<ci>theta_h</ci>
</apply>
<ci>omega_h</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<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</cn>
<ci>omega_h</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="leakage_current">
<variable units="nanoA" public_interface="out" name="i_leak"/>
<variable units="nanoS" name="g_leak" initial_value="2.2"/>
<variable units="millivolt" name="E_leak" initial_value="-57.5"/>
<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_leak</ci>
<apply>
<times/>
<ci>g_leak</ci>
<apply>
<divide/>
<cn cellml:units="volt">1</cn>
<cn cellml:units="millivolt">1000</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_leak</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="tonic_current">
<variable units="nanoA" public_interface="out" name="i_tonic_e"/>
<variable units="nanoS" name="g_tonic_e" initial_value="0.2"/>
<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>
<divide/>
<cn cellml:units="volt">1</cn>
<cn cellml:units="millivolt">1000</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_syn_e</ci>
</apply>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="fast_sodium_current">
<component_ref component="fast_sodium_current_m_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_leak" variable_1="i_leak"/>
</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="potassium_current" component_1="fast_sodium_current"/>
<map_variables variable_2="n" variable_1="n"/>
</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="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 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
<rdf:Description rdf:about="">
<dc:title>Intrinsic Bursters Increase the Robustness of Rythm Generation in an Excitatory Network: Single Non-Pacemaker Cell</dc:title>
<dc:creator rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Lloyd</vCard:Family>
<vCard:Given>Catherine</vCard:Given>
<vCard:Other>May</vCard:Other>
</vCard:N>
<vCard:EMAIL rdf:parseType="Resource">
<rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
</vCard:EMAIL>
<vCard:ORG rdf:parseType="Resource">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
</vCard:ORG>
</dc:creator>
</rdf:Description>
<rdf:Description rdf:about="#purvis_2007">
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>electrophysiology</rdf:li>
<rdf:li>neuron</rdf:li>
<rdf:li>neurobiology</rdf:li>
<rdf:li>pacemaker</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<bqs:reference rdf:parseType="Resource">
<bqs:Pubmed_id>17167061</bqs:Pubmed_id>
<bqs:JournalArticle rdf:parseType="Resource">
<dc:creator>
<rdf:Seq>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Purvis</vCard:Family>
<vCard:Given>L</vCard:Given>
<vCard:Other>K</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Smith</vCard:Family>
<vCard:Given>J</vCard:Given>
<vCard:Other>C</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Koizumi</vCard:Family>
<vCard:Given>H</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Butera</vCard:Family>
<vCard:Given>R</vCard:Given>
<vCard:Other>J</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
</rdf:Seq>
</dc:creator>
<dc:title>
Intrinsic Bursters Increase the Robustness of Rythm Generation in an Excitatory Network
</dc:title>
<dcterms:issued rdf:parseType="Resource">
<dcterms:W3CDTF>2007-02</dcterms:W3CDTF>
</dcterms:issued>
<bqs:Journal rdf:parseType="Resource">
<dc:title>Journal of Neurophysiology</dc:title>
</bqs:Journal>
<bqs:volume>97</bqs:volume>
<bqs:first_page>1515</bqs:first_page>
<bqs:last_page>1526</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
</model>
