Location: Vangoor, Lebeau, Krsmanovic, Sherman, Catt, Stojilkovic, 2000 @ 57e6746b5725 / vangoor_lebeau_krsmanovic_sherman_catt_stojilkovic_2000.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2006-07-09 07:59:10+12:00
Desc:
committing version01 of vangoor_lebeau_krsmanovic_sherman_catt_stojilkovic_2000
Permanent Source URI:
http://models.cellml.org/workspace/vangoor_lebeau_krsmanovic_sherman_catt_stojilkovic_2000/rawfile/57e6746b57251c751c1b3921fdc6389142ec24d4/vangoor_lebeau_krsmanovic_sherman_catt_stojilkovic_2000.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  vangoor_model_2000.xml

CREATED :  8th 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 Van Goor et al's
2000 mathematical model of amplitude-dependent spike broadening and enhanced Ca2+ signaling in GnRH-secreting neurons.

CHANGES:  
  19/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#" cmeta:id="vangoor_lebeau_krsmanovic_sherman_catt_stojilkovic_2000_version01" name="vangoor_lebeau_krsmanovic_sherman_catt_stojilkovic_2000_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Amplitude-Dependent Spike-Broadening In GnRH-Secreting 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>
During an action potential, depolarisation of the plasma membrane causes voltage-gated Ca<superscript>2+</superscript> channels to open, allowing Ca<superscript>2+</superscript> to diffuse into the cell.  The subsequent increase in cytosolic calcium concentration ([Ca<superscript>2+</superscript>]<subscript>i</subscript>) triggers a number of biochemical events involved in the control of cell function.  The frequency, amplitude and duration of the Ca<superscript>2+</superscript> signals determine the specificity of the cellular response.  However, the exact relationship between the action potential waveform, the gating properties of the voltage-gated Ca<superscript>2+</superscript> channels, the pattern of Ca<superscript>2+</superscript> influx and the ionic mechanisms underlying the the agonist-induced Ca<superscript>2+</superscript> spike broadening has not been characterised. </para>

<para>  
In their 2000 paper, Fredrick Van Goor, Andrew P. LeBeau, Lazar Z. Krsmanovic, Arthur Sherman, Kevin J. Catt and Stanko S. Stojikovic examine the ionic mechanism mediating depolarisation-induced Ca<superscript>2+</superscript> spike broadening in GnRH-secreting (GT1) neurons.  The validity of their experimental results were tested by comparing them with a computational model.  They also used this mathematical model to examine the influence of voltage-gated Ca<superscript>2+</superscript> channel gating properties on the ability of action potential broadening to enhance the I<subscript>Ca</subscript>.
</para>

<para>
Their mathematical model is based exclusively on the on the properties of the individual ionic currents characterised in GT1 neurons.  The ionic currents used were a TTX-sensitive I<subscript>Na</subscript>, L- and T-type I<subscript>Ca</subscript>, a voltage sensitive (delayed-rectifying) I<subscript>K(dr)</subscript>, an M-like I<subscript>K</subscript> (I<subscript>M</subscript>), an inward rectifier I<subscript>K</subscript> (I<subscript>ir</subscript>), and a Ca<superscript>2+</superscript> permeable, voltage-insensitive inward leak, I<subscript>d</subscript> (see <xref linkend="fig_cell_diagram"/> below).
</para>


<para>
In their original model, Van Goor <emphasis>et al</emphasis>. used a Hodgkin-Huxley type equation to represent the sodium current.  However, in this description, I<subscript>Na</subscript> exhibited a large surge of activation during the spike repolarisation, and it also did not spike during broadening, suggesting that the model gating kinetics didn't accurately represent those in GTI neurons.  Consequently, they replaced the Hodgkin-Huxley-like description with Kuo and Bean's 1994 model of I<subscript>Na</subscript> (see <xref linkend="fig_reaction_diagram"/> below).  In this description the Na+ channel is composed of subunits with four possible states: deactivated (D), activated (A), deactivated-inhibited (D*), and activated-inhibited (A*).  The open, conducting state (O) is given by A3 in the reaction scheme.
</para>

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

<para>
<ulink url="http://www.biophysj.org/cgi/content/abstract/79/3/1310">Amplitude-Dependent Spike-Broadening and Enhanced Ca<superscript>2+</superscript> Signaling in GnRH-Secreting Neurons</ulink>, Fredrick Van Goor, Andrew P. LeBeau, Lazar Z. Krsmanovic, Arthur Sherman, Kevin J. Catt and Stanko S. Stojikovic, 2000, <ulink url="http://www.biophysj.org/">
            <emphasis>Biophysical Journal</emphasis>
          </ulink>, 79, 1310-1323.  (<ulink url="http://www.biophysj.org/cgi/content/full/79/3/1310">Full text</ulink> and <ulink url="http://www.biophysj.org/cgi/reprint/79/3/1310.pdf">PDF</ulink> versions of the article are available for Journal Members on the BJ website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=10968994&amp;dopt=Abstract">PubMed ID: 10968994</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_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>diagram of the cell model</title>
    </objectinfo>
    <imagedata fileref="vangoor_2000-celldiagram.png"/>
  </imageobject>
</mediaobject>
</informalfigure>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>diagram of the cell model</title>
    </objectinfo>
    <imagedata fileref="vangoor_2000.png"/>
  </imageobject>
</mediaobject>
</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="per_millisecond">
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>
  
  <units name="picoF">
    <unit units="farad" prefix="pico"/>
  </units>
  
  <units name="picoA">
    <unit units="ampere" prefix="pico"/>
  </units>
  
  <units name="nanoS">
    <unit units="siemens" prefix="nano"/>
  </units>
  
  <units name="millimolar">
    <unit units="mole" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="millimolar_per_ms">
    <unit units="millimolar"/>
    <unit units="millisecond" exponent="-1"/>
  </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="Cm" initial_value="14.0"/>
    <variable units="picoA" name="i_app"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="picoA" public_interface="in" name="i_Na"/>
    <variable units="picoA" public_interface="in" name="i_Ca_L"/>
    <variable units="picoA" public_interface="in" name="i_Ca_T"/>
    <variable units="picoA" public_interface="in" name="i_K_dr"/>
    <variable units="picoA" public_interface="in" name="i_M"/>
    <variable units="picoA" public_interface="in" name="i_ir"/>
    <variable units="picoA" public_interface="in" name="i_d"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="membrane_voltage_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> V </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> i_app </ci>
            <apply>
              <plus/>
              <ci> i_Na </ci>
              <ci> i_Ca_L </ci>
              <ci> i_Ca_T </ci>
              <ci> i_K_dr </ci>
              <ci> i_M </ci>
              <ci> i_ir </ci>
              <ci> i_d </ci>
            </apply>
          </apply>
          <ci> Cm </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <rdf:RDF>
    <rdf:Description rdf:about="#sodium_current">
      <cmeta:comment rdf:parseType="Resource">
        <rdf:value>
          In their original model, Van Goor et al. used a Hodgkin-Huxley type 
          equation to represent the sodium current (see CellML description 
          below).  However, in this description, INa exhibited a large surge of 
          activation during the spike repolarisation, and also did not spike 
          during broadening, suggesting that the model gating kinetics don't 
          accurately represent those in GTI neurons.  
          
          Consequently, they replaced the Hodgkin-Huxley-like description with 
          Kuo and Bean's 1994 model of INa.  In this description the Na+ channel 
          is composed of subunits with four possible states: deactivated (D), 
          activated (A), deactivated-inhibited (D*), and activated-inhibited 
          (A*).  The open, conducting state (O) is given by A3 in the reaction 
          scheme.
      
          Following is the original reaction scheme. 
        </rdf:value>
        
        <rdf:value rdf:parseType="Literal">
          <component name="sodium_current">
            <variable units="picoA" public_interface="out" name="i_Na"/>

            <variable units="nanoS" name="g_Na" initial_value="60.0"/>
            <variable units="millivolt" name="V_Na" initial_value="60.0"/>

            <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
            <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>

            <variable units="dimensionless" private_interface="in" name="m"/>
            <variable units="dimensionless" private_interface="in" name="h"/>

            <math xmlns="http://www.w3.org/1998/Math/MathML">
              <apply id="i_Na_calculation">
                <eq/>
                <ci> i_Na </ci>
                <apply>
                  <times/>
                  <ci> g_Na </ci>
                  <apply>
                    <power/>
                    <ci> m </ci>
                    <cn cellml:units="dimensionless"> 3.0 </cn>
                  </apply>
                  <ci> h </ci>
                  <apply>
                    <minus/>
                    <ci> V </ci>
                    <ci> V_Na </ci>
                  </apply>
                </apply>
              </apply>
            </math>
          </component>

          <component name="sodium_current_m_gate">
            <variable units="dimensionless" public_interface="out" name="m"/>

            <variable units="dimensionless" name="m_infinity"/>
            <variable units="millisecond" name="tau_m"/>

            <variable units="millivolt" public_interface="in" name="V"/>
            <variable units="millisecond" public_interface="in" name="time"/>

            <math xmlns="http://www.w3.org/1998/Math/MathML">
              <apply id="sodium_current_m_gate_m_diff_eq">
                <eq/>
                <apply>
                  <diff/>
                  <bvar>
                    <ci> time </ci>
                  </bvar>
                  <ci> m </ci>
                </apply>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> m_infinity </ci>
                    <ci> m </ci>
                  </apply>
                  <ci> tau_m </ci>
                </apply>
              </apply>

              <apply id="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/>
                          <apply>
                            <plus/>
                            <ci> V </ci>
                            <cn cellml:units="millivolt"> 42.1 </cn>
                          </apply>
                        </apply>
                        <cn cellml:units="millivolt"> 4.3 </cn>
                      </apply>
                    </apply>    
                  </apply>
                </apply>
              </apply>

              <apply id="sodium_current_m_gate_tau_m_calculation">
                <eq/>
                <ci> tau_m </ci>
                <apply>
                  <plus/>
                  <apply>
                    <divide/>
                    <cn cellml:units="millisecond"> 4.3 </cn>
                    <apply>
                      <plus/>
                      <apply>
                        <exp/>
                        <apply>
                          <divide/>
                          <apply>
                            <plus/>
                            <ci> V </ci>
                            <cn cellml:units="millivolt"> 47.0 </cn>
                          </apply>
                          <cn cellml:units="millivolt"> 11.0 </cn>
                        </apply>
                      </apply>
                      <apply>
                        <times/>
                        <cn cellml:units="dimensionless"> 2.0 </cn>
                        <apply>
                          <exp/>
                          <apply>
                            <divide/>
                            <apply>
                              <times/>
                              <cn cellml:units="dimensionless"> -2.0 </cn>
                              <apply>
                                <plus/>
                                <ci> V </ci>
                                <cn cellml:units="millivolt"> 47.0 </cn>
                              </apply>
                            </apply>
                            <cn cellml:units="millivolt"> 11.0 </cn>
                          </apply>
                        </apply>
                      </apply>
                    </apply>
                  </apply>
                  <cn cellml:units="millisecond"> 0.1 </cn>
                </apply>
              </apply>
            </math>
          </component>

          <component name="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="millivolt" public_interface="in" name="V"/>
            <variable units="millisecond" public_interface="in" name="time"/>

            <math xmlns="http://www.w3.org/1998/Math/MathML">
              <apply id="sodium_current_h_gate_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="sodium_current_h_gate_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>
                          <plus/>
                          <ci> V </ci>
                          <cn cellml:units="millivolt"> 68.2 </cn>
                        </apply>
                        <cn cellml:units="millivolt"> 10.8 </cn>
                      </apply>
                    </apply>    
                  </apply>
                </apply>
              </apply>

              <apply id="tau_h_calculation">
                <eq/>
                <ci> tau_h </ci>
                <apply>
                  <divide/>
                  <cn cellml:units="millisecond"> 150.0 </cn>
                  <apply>
                    <plus/>
                    <apply>
                      <exp/>
                      <apply>
                        <divide/>
                        <apply>
                          <plus/>
                          <ci> V </ci>
                          <cn cellml:units="millivolt"> 80.0 </cn>
                        </apply>
                        <cn cellml:units="millivolt"> 19.0 </cn>
                      </apply>
                    </apply>
                    <apply>
                      <times/>
                      <cn cellml:units="dimensionless"> 2.0 </cn>
                      <apply>
                        <exp/>
                        <apply>
                          <divide/>
                          <apply>
                            <times/>
                            <cn cellml:units="dimensionless"> -2.0 </cn>
                            <apply>
                              <plus/>
                              <ci> V </ci>
                              <cn cellml:units="millivolt"> 80.0 </cn>
                            </apply>
                          </apply>
                          <cn cellml:units="millivolt"> 19.0 </cn>
                        </apply>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </math>
          </component>
   
        </rdf:value>
        <!-- The creator of the comment. -->
        <dc:creator rdf:parseType="Resource">
          <vCard:FN>Catherine Lloyd</vCard:FN>
        </dc:creator>
      </cmeta:comment>
    </rdf:Description>
  </rdf:RDF>
  
  <component name="sodium_current">
    <variable units="picoA" public_interface="out" name="i_Na"/>
    
    <variable units="nanoS" name="g_Na" initial_value="60.0"/>
    <variable units="millivolt" name="V_Na" initial_value="60.0"/>
    <variable units="dimensionless" name="O"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millimolar" public_interface="in" name="A"/>
    
     <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>
          <ci> O </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_Na </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="O_calculation">
        <eq/>
        <ci> O </ci>
        <apply>
          <power/>
          <ci> A </ci>
          <cn cellml:units="dimensionless"> 3.0 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="A">
    <variable units="millimolar" public_interface="out" name="A"/>
    
    <variable units="millimolar_per_ms" public_interface="in" name="delta_A_rxn_1"/>
    <variable units="millimolar_per_ms" public_interface="in" name="delta_A_rxn_4"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply id="A_diff_eq">
        <eq/>
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>A</ci> 
        </apply> 
        <apply>
          <plus/>
          <ci>delta_A_rxn_1</ci>
          <ci>delta_A_rxn_4</ci>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="inhibited_A">
    <variable units="millimolar" public_interface="out" name="inhibited_A"/>
    
    <variable units="millimolar_per_ms" public_interface="in" name="delta_inhibited_A_rxn_1"/>
    <variable units="millimolar_per_ms" public_interface="in" name="delta_inhibited_A_rxn_2"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply id="inhibited_A_diff_eq">
        <eq/>
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>inhibited_A</ci> 
        </apply> 
        <apply>
          <plus/>
          <ci>delta_inhibited_A_rxn_1</ci>
          <ci>delta_inhibited_A_rxn_2</ci>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="D">
    <variable units="millimolar" public_interface="out" name="D"/>
    
    <variable units="millimolar_per_ms" public_interface="in" name="delta_D_rxn_3"/>
    <variable units="millimolar_per_ms" public_interface="in" name="delta_D_rxn_4"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply id="D_diff_eq">
        <eq/>
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>D</ci> 
        </apply> 
        <apply>
          <plus/>
          <ci>delta_D_rxn_3</ci>
          <ci>delta_D_rxn_4</ci>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="inhibited_D">
    <variable units="millimolar" public_interface="out" name="inhibited_D"/>
    
    <variable units="millimolar_per_ms" public_interface="in" name="delta_inhibited_D_rxn_2"/>
    <variable units="millimolar_per_ms" public_interface="in" name="delta_inhibited_D_rxn_3"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply id="inhibited_D_diff_eq">
        <eq/>
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>inhibited_D</ci> 
        </apply> 
        <apply>
          <plus/>
          <ci>delta_inhibited_D_rxn_2</ci>
          <ci>delta_inhibited_D_rxn_3</ci>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="reaction1">
    <variable units="millimolar" public_interface="in" name="A"/>
    <variable units="millimolar" public_interface="in" name="inhibited_A"/>
    <variable units="per_millisecond" public_interface="out" name="k1" initial_value="0.3"/>
    <variable units="per_millisecond" public_interface="out" name="k1_" initial_value="0.03"/>
    <variable units="millimolar_per_ms" public_interface="out" name="delta_A_rxn_1"/>
    <variable units="millimolar_per_ms" public_interface="out" name="delta_inhibited_A_rxn_1"/>
    <variable units="millimolar_per_ms" name="rate"/>
    <reaction reversible="yes">
      <variable_ref variable="A">
        <role stoichiometry="1" direction="forward" delta_variable="delta_A_rxn_1" role="reactant"/>
      </variable_ref>
      <variable_ref variable="inhibited_A">
        <role stoichiometry="1" direction="forward" delta_variable="delta_inhibited_A_rxn_1" role="product"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci> rate </ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci>k1</ci>
                  <ci>A</ci>
                </apply>
                <apply>
                  <minus/>
                  <apply>
                    <times/>
                    <ci>k1_</ci>
                    <ci>inhibited_A</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
  </component>
  
  <component name="reaction2">
    <variable units="millimolar" public_interface="in" name="inhibited_A"/>
    <variable units="millimolar" public_interface="in" name="inhibited_D"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="per_millisecond" public_interface="in" name="k1"/>
    <variable units="per_millisecond" public_interface="in" name="k1_"/>
    <variable units="per_millisecond" public_interface="in" name="k3"/>
    <variable units="per_millisecond" public_interface="in" name="k3_"/>
    <variable units="millisecond" public_interface="out" name="alpha"/>
    <variable units="millisecond" public_interface="out" name="beta"/>
    <variable units="dimensionless" name="a"/>
    <variable units="millimolar_per_ms" public_interface="out" name="delta_inhibited_A_rxn_2"/>
    <variable units="millimolar_per_ms" public_interface="out" name="delta_inhibited_D_rxn_2"/>
    <variable units="millimolar_per_ms" name="rate"/>
    <reaction reversible="yes">
      <variable_ref variable="inhibited_A">
        <role stoichiometry="1" direction="forward" delta_variable="delta_inhibited_A_rxn_2" role="reactant"/>
      </variable_ref>
      <variable_ref variable="inhibited_D">
        <role stoichiometry="1" direction="forward" delta_variable="delta_inhibited_D_rxn_2" role="product"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci> rate </ci>
              <apply>
                <plus/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <ci>beta</ci>
                    <ci>inhibited_A</ci>
                  </apply>
                  <ci>a</ci>
                </apply>
                <apply>
                  <minus/>
                  <apply>
                    <times/>
                    <ci>alpha</ci>
                    <ci>a</ci>
                    <ci>inhibited_D</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_calculation">
        <eq/>
        <ci> alpha </ci>
        <apply>
          <divide/>
          <cn cellml:units="millisecond"> 10.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 6.0 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 10.0 </cn>
              </apply>
            </apply>    
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_calculation">
        <eq/>
        <ci> beta </ci>
        <apply>
          <divide/>
          <cn cellml:units="millisecond"> 10.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 54.4 </cn>
                </apply>
                <cn cellml:units="millivolt"> 4.6 </cn>
              </apply>
            </apply>    
          </apply>
        </apply>
      </apply>
      
      <apply id="a_calculation">
        <eq/>
        <ci> a </ci>
        <apply>
          <power/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> k1 </ci>
              <ci> k3_ </ci>
            </apply>
            <apply>
              <times/>
              <ci> k1_ </ci>
              <ci> k3 </ci>
            </apply>
          </apply>
          <cn cellml:units="dimensionless"> 0.5 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="reaction3">
    <variable units="millimolar" public_interface="in" name="D"/>
    <variable units="millimolar" public_interface="in" name="inhibited_D"/>
    <variable units="per_millisecond" public_interface="out" name="k3" initial_value="0.001"/>
    <variable units="per_millisecond" public_interface="out" name="k3_" initial_value="0.01"/>
    <variable units="millimolar_per_ms" public_interface="out" name="delta_D_rxn_3"/>
    <variable units="millimolar_per_ms" public_interface="out" name="delta_inhibited_D_rxn_3"/>
    <variable units="millimolar_per_ms" name="rate"/>
    <reaction reversible="yes">
      <variable_ref variable="D">
        <role stoichiometry="1" direction="forward" delta_variable="delta_D_rxn_3" role="reactant"/>
      </variable_ref>
      <variable_ref variable="inhibited_D">
        <role stoichiometry="1" direction="forward" delta_variable="delta_inhibited_D_rxn_3" role="product"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci> rate </ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci>k3</ci>
                  <ci>D</ci>
                </apply>
                <apply>
                  <minus/>
                  <apply>
                    <times/>
                    <ci>k3_</ci>
                    <ci>inhibited_D</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
  </component>
  
  <component name="reaction4">
    <variable units="millimolar" public_interface="in" name="A"/>
    <variable units="millimolar" public_interface="in" name="D"/>
    <variable units="millisecond" public_interface="in" name="alpha"/>
    <variable units="millisecond" public_interface="in" name="beta"/>
    <variable units="millimolar_per_ms" public_interface="out" name="delta_A_rxn_4"/>
    <variable units="millimolar_per_ms" public_interface="out" name="delta_D_rxn_4"/>
    <variable units="millimolar_per_ms" name="rate"/>
    <reaction reversible="yes">
      <variable_ref variable="A">
        <role stoichiometry="1" direction="forward" delta_variable="delta_A_rxn_4" role="product"/>
      </variable_ref>
      <variable_ref variable="D">
        <role stoichiometry="1" direction="forward" delta_variable="delta_D_rxn_4" role="reactant"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci> rate </ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci>alpha</ci>
                  <ci>D</ci>
                </apply>
                <apply>
                  <minus/>
                  <apply>
                    <times/>
                    <ci>beta</ci>
                    <ci>A</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
  </component>
  
  <component name="L_type_calcium_current">
    <variable units="picoA" public_interface="out" name="i_Ca_L"/>
    <variable units="millivolt" public_interface="out" name="V_Ca" initial_value="100.0"/>
    
    <variable units="nanoS" name="g_Ca_L" initial_value="1.3"/>
    
    <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"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Ca_L_calculation">
        <eq/>
        <ci> i_Ca_L </ci>
        <apply>
          <times/>
          <ci> g_Ca_L </ci>
          <apply>
            <power/>
            <ci> m </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_calcium_current_m_gate">
    <variable units="dimensionless" public_interface="out" name="m"/>
    
    <variable units="dimensionless" name="m_infinity"/>
    <variable units="millisecond" name="tau_m"/>
    <variable units="millivolt" name="Vh" initial_value="40.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="L_type_calcium_current_m_gate_m_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> m </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> m_infinity </ci>
            <ci> m </ci>
          </apply>
          <ci> tau_m </ci>
        </apply>
      </apply>
      
      <apply id="L_type_calcium_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/>
                  <apply>
                    <minus/>
                    <ci> V </ci>
                    <ci> Vh </ci>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 12.0 </cn>
              </apply>
            </apply>    
          </apply>
        </apply>
      </apply>
      
      <apply id="L_type_calcium_current_m_gate_tau_m_calculation">
        <eq/>
        <ci> tau_m </ci>
        <apply>
          <divide/>
          <cn cellml:units="millisecond"> 5.0 </cn>
          <apply>
            <plus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 15.0 </cn>
                </apply>
                <cn cellml:units="millivolt"> 25.0 </cn>
              </apply>
            </apply>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 15.0 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 25.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="T_type_calcium_current">
    <variable units="picoA" public_interface="out" name="i_Ca_T"/>
     
    <variable units="nanoS" name="g_Ca_T" initial_value="0.94"/>
    
    <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="V_Ca"/>
    
    <variable units="dimensionless" private_interface="in" name="m"/>
    <variable units="dimensionless" private_interface="in" name="h"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Ca_T_calculation">
        <eq/>
        <ci> i_Ca_T </ci>
        <apply>
          <times/>
          <ci> g_Ca_T </ci>
          <apply>
            <power/>
            <ci> m </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <ci> h </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="T_type_calcium_current_m_gate">
    <variable units="dimensionless" public_interface="out" name="m"/>
    
    <variable units="dimensionless" name="m_infinity"/>
    <variable units="millisecond" name="tau_m"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="T_type_calcium_current_m_gate_m_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> m </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> m_infinity </ci>
            <ci> m </ci>
          </apply>
          <ci> tau_m </ci>
        </apply>
      </apply>
      
      <apply id="T_type_calcium_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/>
                  <apply>
                    <minus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 56.1 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 10.0 </cn>
              </apply>
            </apply>    
          </apply>
        </apply>
      </apply>
      
      <apply id="T_type_calcium_current_m_gate_tau_m_calculation">
        <eq/>
        <ci> tau_m </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="millisecond"> 7.0 </cn>
            <apply>
              <plus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 50.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 9.0 </cn>
                </apply>
              </apply>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 50.0 </cn>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt"> 9.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="millisecond"> 0.8 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="T_type_calcium_current_h_gate">
    <variable units="dimensionless" public_interface="out" name="h"/>
    
    <variable units="dimensionless" name="h_infinity"/>
    <variable units="millisecond" name="tau_h" initial_value="22.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="T_type_calcium_current_h_gate_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="T_type_calcium_current_h_gate_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>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 86.4 </cn>
                </apply>
                <cn cellml:units="millivolt"> 4.7 </cn>
              </apply>
            </apply>    
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="delayed_rectifier_K_channel_current">
    <variable units="picoA" public_interface="out" name="i_K_dr"/>
    <variable units="millivolt" public_interface="out" name="V_K" initial_value="-80.0"/>
    
    <variable units="nanoS" name="g_K_dr" initial_value="20.0"/>   
   
    <variable units="second" 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"/>
    <variable units="dimensionless" private_interface="in" name="h"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_K_dr_calculation">
        <eq/>
        <ci> i_K_dr </ci>
        <apply>
          <times/>
          <ci> g_K_dr </ci>
          <apply>
            <power/>
            <ci> n </ci>
            <cn cellml:units="dimensionless"> 4.0 </cn>
          </apply>
          <ci> h </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="delayed_rectifier_K_channel_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="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="delayed_rectifier_K_channel_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="delayed_rectifier_K_channel_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/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 25.0 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 15.0 </cn>
              </apply>
            </apply>    
          </apply>
        </apply>
      </apply>
      
      <apply id="delayed_rectifier_K_channel_current_n_gate_tau_n_calculation">
        <eq/>
        <ci> tau_n </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="millisecond"> 15.0 </cn>
            <apply>
              <plus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 30.0 </cn>
                  </apply>
                  <cn cellml:units="millivolt"> 15.0 </cn>
                </apply>
              </apply>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 30.0 </cn>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt"> 15.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="millisecond"> 1.0 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="delayed_rectifier_K_channel_current_h_gate">
    <variable units="dimensionless" public_interface="out" name="h"/>
    
    <variable units="dimensionless" name="h_infinity"/>
    <variable units="millisecond" name="tau_h" initial_value="1000.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="delayed_rectifier_K_channel_current_h_gate_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="delayed_rectifier_K_channel_current_h_gate_infinity_calculation">
        <eq/>
        <ci> h_infinity </ci>
        <apply>
          <plus/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 0.7 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 35.0 </cn>
                    </apply>
                  </apply>
                  <cn cellml:units="millivolt"> 10.0 </cn>
                </apply>
              </apply>
            </apply>    
          </apply>
          <cn cellml:units="dimensionless"> 0.3 </cn>
        </apply>
      </apply>
    </math>
  </component>  
  
  <component name="M_like_K_current">
    <variable units="picoA" public_interface="out" name="i_M"/>
    
    <variable units="nanoS" name="g_M" initial_value="4.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="millivolt" public_interface="in" name="V_K"/>
    
    <variable units="dimensionless" private_interface="in" name="n"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_M_calculation">
        <eq/>
        <ci> i_M </ci>
        <apply>
          <times/>
          <ci> g_M </ci>
          <ci> n </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="M_like_K_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="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="M_like_K_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="M_like_K_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/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 37.0 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 4.0 </cn>
              </apply>
            </apply>    
          </apply>
        </apply>
      </apply>
      
      <apply id="M_like_K_current_n_gate_tau_n_calculation">
        <eq/>
        <ci> tau_n </ci>
        <apply>
          <divide/>
          <cn cellml:units="millisecond"> 80.0 </cn>
          <apply>
            <plus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 30.0 </cn>
                </apply>
                <cn cellml:units="millivolt"> 15.0 </cn>
              </apply>
            </apply>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <apply>
                    <plus/>
                    <ci> V </ci>
                    <cn cellml:units="millivolt"> 30.0 </cn>
                  </apply>
                </apply>
                <cn cellml:units="millivolt"> 15.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="inward_rectifier_K_current">
    <variable units="picoA" public_interface="out" name="i_ir"/>
    
    <variable units="nanoS" name="g_ir" initial_value="1.71"/>   
    
    <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="V_K"/>
    
    <variable units="dimensionless" private_interface="in" name="n_infinity"/>
    
    <math xmlns="http://www.w3.org/1998/irath/irathirL">
      <apply id="i_ir_calculation">
        <eq/>
        <ci> i_ir </ci>
        <apply>
          <times/>
          <ci> g_ir </ci>
          <ci> n_infinity </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="inward_rectifier_K_current_n_gate">
    <variable units="dimensionless" public_interface="out" name="n_infinity"/>   
    <variable units="millivolt" public_interface="in" name="V"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="inward_rectifier_K_current_n_gate_n_infinity_calculation">
        <eq/>
        <ci> n_infinity </ci>
        <apply>
          <plus/>
          <apply>
            <times/>
            <cn cellml:units="dimensionless"> 0.8 </cn>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <plus/>
                      <ci> V </ci>
                      <cn cellml:units="millivolt"> 80.0 </cn>
                    </apply>
                    <cn cellml:units="millivolt"> 12.0 </cn>
                  </apply>
                </apply>    
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless"> 0.2 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="inward_leak_current">
    <variable units="picoA" public_interface="out" name="i_d"/>
    
    <variable units="nanoS" name="g_d" initial_value="0.044"/>   
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millivolt" public_interface="in" name="V_Ca"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_d_calculation">
        <eq/>
        <ci> i_d </ci>
        <apply>
          <times/>
          <ci> g_d </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="sodium_current"/>
      <component_ref component="L_type_calcium_current">
        <component_ref component="L_type_calcium_current_m_gate"/>
      </component_ref>
      <component_ref component="T_type_calcium_current">
        <component_ref component="T_type_calcium_current_m_gate"/>
        <component_ref component="T_type_calcium_current_h_gate"/>
      </component_ref>
      <component_ref component="delayed_rectifier_K_channel_current">
        <component_ref component="delayed_rectifier_K_channel_current_n_gate"/>
        <component_ref component="delayed_rectifier_K_channel_current_h_gate"/>
      </component_ref>
      <component_ref component="M_like_K_current">
        <component_ref component="M_like_K_current_n_gate"/>
      </component_ref>
      <component_ref component="inward_rectifier_K_current">
        <component_ref component="inward_rectifier_K_current_n_gate"/>
      </component_ref>
      <component_ref component="inward_leak_current"/>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="encapsulation"/>
    <component_ref component="L_type_calcium_current">
      <component_ref component="L_type_calcium_current_m_gate"/>
    </component_ref>
    <component_ref component="T_type_calcium_current">
      <component_ref component="T_type_calcium_current_m_gate"/>
      <component_ref component="T_type_calcium_current_h_gate"/>
    </component_ref>
    <component_ref component="delayed_rectifier_K_channel_current">
      <component_ref component="delayed_rectifier_K_channel_current_n_gate"/>
      <component_ref component="delayed_rectifier_K_channel_current_h_gate"/>
    </component_ref>
    <component_ref component="M_like_K_current">
      <component_ref component="M_like_K_current_n_gate"/>
    </component_ref>
    <component_ref component="inward_rectifier_K_current">
      <component_ref component="inward_rectifier_K_current_n_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="sodium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="L_type_calcium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="T_type_calcium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="delayed_rectifier_K_channel_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="M_like_K_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="inward_rectifier_K_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="inward_leak_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="A"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="inhibited_A"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="D"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="inhibited_D"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="membrane" component_1="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="L_type_calcium_current"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Ca_L" variable_1="i_Ca_L"/>
  </connection>
  
  <connection>
    <map_components component_2="membrane" component_1="T_type_calcium_current"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Ca_T" variable_1="i_Ca_T"/>
  </connection>
  
  <connection>
    <map_components component_2="membrane" component_1="delayed_rectifier_K_channel_current"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_K_dr" variable_1="i_K_dr"/>
  </connection>
  
  <connection>
    <map_components component_2="membrane" component_1="M_like_K_current"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_M" variable_1="i_M"/>
  </connection>
  
  <connection>
    <map_components component_2="membrane" component_1="inward_rectifier_K_current"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_ir" variable_1="i_ir"/>
  </connection>
  
  <connection>
    <map_components component_2="membrane" component_1="inward_leak_current"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_d" variable_1="i_d"/>
  </connection>
  
  <connection>
    <map_components component_2="membrane" component_1="reaction2"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="A" component_1="sodium_current"/>
    <map_variables variable_2="A" variable_1="A"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction4" component_1="A"/>
    <map_variables variable_2="A" variable_1="A"/>
    <map_variables variable_2="delta_A_rxn_4" variable_1="delta_A_rxn_4"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction1" component_1="A"/>
    <map_variables variable_2="A" variable_1="A"/>
    <map_variables variable_2="delta_A_rxn_1" variable_1="delta_A_rxn_1"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction1" component_1="inhibited_A"/>
    <map_variables variable_2="inhibited_A" variable_1="inhibited_A"/>
    <map_variables variable_2="delta_inhibited_A_rxn_1" variable_1="delta_inhibited_A_rxn_1"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction2" component_1="inhibited_A"/>
    <map_variables variable_2="inhibited_A" variable_1="inhibited_A"/>
    <map_variables variable_2="delta_inhibited_A_rxn_2" variable_1="delta_inhibited_A_rxn_2"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction4" component_1="D"/>
    <map_variables variable_2="D" variable_1="D"/>
    <map_variables variable_2="delta_D_rxn_4" variable_1="delta_D_rxn_4"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction3" component_1="D"/>
    <map_variables variable_2="D" variable_1="D"/>
    <map_variables variable_2="delta_D_rxn_3" variable_1="delta_D_rxn_3"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction3" component_1="inhibited_D"/>
    <map_variables variable_2="inhibited_D" variable_1="inhibited_D"/>
    <map_variables variable_2="delta_inhibited_D_rxn_3" variable_1="delta_inhibited_D_rxn_3"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction2" component_1="inhibited_D"/>
    <map_variables variable_2="inhibited_D" variable_1="inhibited_D"/>
    <map_variables variable_2="delta_inhibited_D_rxn_2" variable_1="delta_inhibited_D_rxn_2"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction2" component_1="reaction1"/>
    <map_variables variable_2="k1" variable_1="k1"/>
    <map_variables variable_2="k1_" variable_1="k1_"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction2" component_1="reaction3"/>
    <map_variables variable_2="k3" variable_1="k3"/>
    <map_variables variable_2="k3_" variable_1="k3_"/>
  </connection>
  
  <connection>
    <map_components component_2="reaction4" component_1="reaction2"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
    <map_variables variable_2="beta" variable_1="beta"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_calcium_current" component_1="L_type_calcium_current"/>
    <map_variables variable_2="V_Ca" variable_1="V_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="inward_leak_current" component_1="L_type_calcium_current"/>
    <map_variables variable_2="V_Ca" variable_1="V_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="M_like_K_current" component_1="delayed_rectifier_K_channel_current"/>
    <map_variables variable_2="V_K" variable_1="V_K"/>
  </connection>
  
  <connection>
    <map_components component_2="inward_rectifier_K_current" component_1="delayed_rectifier_K_channel_current"/>
    <map_variables variable_2="V_K" variable_1="V_K"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_calcium_current_m_gate" component_1="L_type_calcium_current"/>
    <map_variables variable_2="m" variable_1="m"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_calcium_current_m_gate" component_1="T_type_calcium_current"/>
    <map_variables variable_2="m" variable_1="m"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_calcium_current_h_gate" component_1="T_type_calcium_current"/>
    <map_variables variable_2="h" variable_1="h"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="delayed_rectifier_K_channel_current_n_gate" component_1="delayed_rectifier_K_channel_current"/>
    <map_variables variable_2="n" variable_1="n"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="delayed_rectifier_K_channel_current_h_gate" component_1="delayed_rectifier_K_channel_current"/>
    <map_variables variable_2="h" variable_1="h"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="M_like_K_current_n_gate" component_1="M_like_K_current"/>
    <map_variables variable_2="n" variable_1="n"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="inward_rectifier_K_current_n_gate" component_1="inward_rectifier_K_current"/>
    <map_variables variable_2="n_infinity" variable_1="n_infinity"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>

<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <rdf:Bag rdf:about="rdf:#536cdbb2-9001-459a-8c24-44fbe3ce12ad">
    <rdf:li>Neuron</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#a3f120b9-0581-4c7f-8fc5-6b87f49d562f">
    <rdf:li rdf:resource="rdf:#1807f78f-5475-44b7-86f0-35549b1f47a5"/>
    <rdf:li rdf:resource="rdf:#6a20fa5a-e884-4601-a361-0a7a23f7f95d"/>
    <rdf:li rdf:resource="rdf:#ea92b91b-6384-423d-9b63-5a6043184efc"/>
    <rdf:li rdf:resource="rdf:#fffb4327-c414-4043-a5b8-446db5c0679a"/>
    <rdf:li rdf:resource="rdf:#fefb8af0-16c6-4661-b770-36cc893efdb9"/>
    <rdf:li rdf:resource="rdf:#25c2b5f1-382b-4553-97d7-2e443c48bf9b"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#f51a8011-2574-443c-b86d-ad3e947d4efb">
    <dc:creator rdf:resource="rdf:#ccd3d8b5-5361-4ba9-a423-a50a31804ebe"/>
    <rdf:value>
          This is the CellML description of Van Goor et al's 2000 mathematical 
          model of amplitude-dependent spike broadening and enhanced Ca2+ 
          signaling in GnRH-secreting neurons.
        </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fefb8af0-16c6-4661-b770-36cc893efdb9">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#9234d85c-9d33-475e-85c6-49a0d8138b7c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6a20fa5a-e884-4601-a361-0a7a23f7f95d">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#5c58826d-d388-4580-b385-f6310ee41e55"/>
  </rdf:Description>
  <rdf:Description rdf:about="#vangoor_lebeau_krsmanovic_sherman_catt_stojilkovic_2000_version01">
    <dc:title>
        Van Goor et al's 2000 mathematical model of amplitude-dependent spike 
        broadening and enhanced Ca2+ signaling in GnRH-secreting neurons.
     </dc:title>
    <cmeta:bio_entity>Neuron</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#f51a8011-2574-443c-b86d-ad3e947d4efb"/>
    <bqs:reference rdf:resource="rdf:#87c27679-a2ea-47c0-a25c-e34bfecb9008"/>
    <bqs:reference rdf:resource="rdf:#97c7a49f-79bf-41d0-bc73-a8b97f2d3354"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ccd3d8b5-5361-4ba9-a423-a50a31804ebe">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#171e259a-fd75-4a57-9e86-3bfcf0418b4d">
    <vCard:Given>Arthur</vCard:Given>
    <vCard:Family>Sherman</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0876e16f-0dee-46c9-b639-ccc4dfc21d89">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#75785157-121b-4dfc-9b01-4b2f891b6b70">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#536cdbb2-9001-459a-8c24-44fbe3ce12ad"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#79e7a6d7-2b4a-42ab-a2ca-c9bd79822d9c">
    <vCard:Given>Lazar</vCard:Given>
    <vCard:Family>Krsmanovic</vCard:Family>
    <vCard:Other>Z</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6dac19d0-7dc5-45c9-b125-0026ff1c8ee8">
    <dc:title>Biophysical Journal</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#417a1a58-8f1e-4533-a127-3af6f9f078e8">
    <dcterms:modified rdf:resource="rdf:#845245b5-00d9-4d2b-97ed-9c88b5c2443b"/>
    <rdf:value>
          Made mathml id's unique
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#ff903244-afb3-4f85-9c94-e44a1235ff87"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#25c2b5f1-382b-4553-97d7-2e443c48bf9b">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#91be835b-3104-47d0-b3fb-ca4ef4d606cb"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7fd14a2d-f514-4b4b-8220-a15d9cd9f4a8">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b0f3d4dd-7112-4614-8e35-befb53e3d163">
    <vCard:ORG rdf:resource="rdf:#2e1a808e-eec0-48e9-bba0-586b22f6eb92"/>
    <vCard:EMAIL rdf:resource="rdf:#b0f045c7-3f75-4103-9fe4-60d4c4fc2727"/>
    <vCard:N rdf:resource="rdf:#7fd14a2d-f514-4b4b-8220-a15d9cd9f4a8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ff903244-afb3-4f85-9c94-e44a1235ff87">
    <vCard:N rdf:resource="rdf:#fb815993-03a1-4eaa-8995-d520635db10b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9234d85c-9d33-475e-85c6-49a0d8138b7c">
    <vCard:Given>Kevin</vCard:Given>
    <vCard:Family>Catt</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fb815993-03a1-4eaa-8995-d520635db10b">
    <vCard:Given>Peter</vCard:Given>
    <vCard:Family>Villiger</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b0f045c7-3f75-4103-9fe4-60d4c4fc2727">
    <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:#d9a4a466-d03f-43ff-a8a9-156552b6dc31">
    <dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d9936a78-4c65-4583-a70f-29c9d553a14c">
    <dc:creator rdf:resource="rdf:#a3f120b9-0581-4c7f-8fc5-6b87f49d562f"/>
    <dc:title>
            Amplitude-Dependent Spike-Broadening and Enhanced Ca2+ Signaling in              GnRH-Secreting Neurons
          </dc:title>
    <bqs:volume>79</bqs:volume>
    <bqs:first_page>1310</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#6dac19d0-7dc5-45c9-b125-0026ff1c8ee8"/>
    <dcterms:issued rdf:resource="rdf:#11cdca98-79dc-4577-95e8-de7def01390d"/>
    <bqs:last_page>1323</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#94470040-bffd-4f19-a22d-6b426468800f">
    <vCard:N rdf:resource="rdf:#ba343f9d-2270-46a1-9386-8054ab4eb7af"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d66e41a0-d139-4823-9110-09726f264ddd">
    <dcterms:modified rdf:resource="rdf:#d9a4a466-d03f-43ff-a8a9-156552b6dc31"/>
    <rdf:value>
          Added publication date information.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#94470040-bffd-4f19-a22d-6b426468800f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8dd53c63-8274-459e-9c29-e971600bdc8b">
    <dcterms:W3CDTF>2002-05-08</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e6ed0148-8119-49f8-b447-c539c0db4f27">
    <dcterms:modified rdf:resource="rdf:#58709ae9-167b-4a8f-be27-a213691b51f8"/>
    <rdf:value>
          Added more metadata.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#37157642-9ea6-4178-8fd6-69312c54d0d8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#97c7a49f-79bf-41d0-bc73-a8b97f2d3354">
    <bqs:Pubmed_id>10968994</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#d9936a78-4c65-4583-a70f-29c9d553a14c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ea92b91b-6384-423d-9b63-5a6043184efc">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#79e7a6d7-2b4a-42ab-a2ca-c9bd79822d9c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#663e6803-4a84-4450-b085-c284d019d158">
    <vCard:Given>Fredrick</vCard:Given>
    <vCard:Family>Van Goor</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2e1a808e-eec0-48e9-bba0-586b22f6eb92">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fffb4327-c414-4043-a5b8-446db5c0679a">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#171e259a-fd75-4a57-9e86-3bfcf0418b4d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ba343f9d-2270-46a1-9386-8054ab4eb7af">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#91be835b-3104-47d0-b3fb-ca4ef4d606cb">
    <vCard:Given>Stanko</vCard:Given>
    <vCard:Family>Stojilkovic</vCard:Family>
    <vCard:Other>S</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Bioengineering Institute
      </dc:publisher>
    <cmeta:modification rdf:resource="rdf:#417a1a58-8f1e-4533-a127-3af6f9f078e8"/>
    <cmeta:modification rdf:resource="rdf:#d66e41a0-d139-4823-9110-09726f264ddd"/>
    <cmeta:modification rdf:resource="rdf:#e6ed0148-8119-49f8-b447-c539c0db4f27"/>
    <dcterms:created rdf:resource="rdf:#8dd53c63-8274-459e-9c29-e971600bdc8b"/>
    <dc:creator rdf:resource="rdf:#b0f3d4dd-7112-4614-8e35-befb53e3d163"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#87c27679-a2ea-47c0-a25c-e34bfecb9008">
    <dc:subject rdf:resource="rdf:#75785157-121b-4dfc-9b01-4b2f891b6b70"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#845245b5-00d9-4d2b-97ed-9c88b5c2443b">
    <dcterms:W3CDTF>2005-04-20</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#11cdca98-79dc-4577-95e8-de7def01390d">
    <dcterms:W3CDTF>2000-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1807f78f-5475-44b7-86f0-35549b1f47a5">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#663e6803-4a84-4450-b085-c284d019d158"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#58709ae9-167b-4a8f-be27-a213691b51f8">
    <dcterms:W3CDTF>2002-07-19</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#37157642-9ea6-4178-8fd6-69312c54d0d8">
    <vCard:N rdf:resource="rdf:#0876e16f-0dee-46c9-b639-ccc4dfc21d89"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c58826d-d388-4580-b385-f6310ee41e55">
    <vCard:Given>Andrew</vCard:Given>
    <vCard:Family>LeBeau</vCard:Family>
    <vCard:Other>P</vCard:Other>
  </rdf:Description>
</rdf:RDF>
</model>