Location: Calcium Dynamics Annotations @ 4d6408caaf3e / winograd_destexhe_sanchezvives_2008 / model / winograd_destexhe_sanchezvives_2008.cellml

Author:
neilstephen2001 <npar276@aucklanduni.ac.nz>
Date:
2021-01-26 14:35:03+13:00
Desc:
modifying rdf files
Permanent Source URI:
http://models.cellml.org/workspace/7f1/rawfile/4d6408caaf3e7b1b86d33d53d26105360ef95347/winograd_destexhe_sanchezvives_2008/model/winograd_destexhe_sanchezvives_2008.cellml

<?xml version="1.0" encoding="UTF-8"?>
<model xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:semsim="http://www.bhi.washington.edu/SemSim#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="winograd_destexhe_sanchezvives_2008" cmeta:id="winograd_destexhe_sanchez_vives_2008">
  <documentation xmlns="http://cellml.org/tmp-documentation">
    <article>
      <articleinfo>
        <title>Hyperpolarization-Activated Graded Persistent Activity in the Prefrontal Cortex</title>
        <author>
          <firstname>Geoffrey</firstname>
          <surname>Nunns</surname>
          <affiliation>
            <shortaffil>Bioengineering Institute, The University of Auckland</shortaffil>
          </affiliation>
        </author>
      </articleinfo>
      <section id="sec_status">
        <title>Model Status</title>
        <para>This model runs in PCEnv and is based off the NEURON code for this paper. The units are continuous throughout but the model does not recreate the published results. It utilizes DAE solving techniques and advances in these could result in better output.</para>
      </section>
      <sect1 id="sec_structure">
        <title>Model Structure</title>
        <para>Abstract: We describe a phenomenon of hyperpolarization-activated graded persistent activity (HAGPA) in prefrontal cortex neurons. Successive hyperpolarizing pulses induced increasingly higher rates of tonic firing that remained stable for tens of seconds, allowing the neuron to retain a memory of the previous history of stimulation. This phenomenon occurred at the cellular level and in the absence of neuromodulators. Neurons with HAGPA had a sag during hyperpolarization, and blocking h-current eliminated the sag and prevented HAGPA, suggesting that the activation of this hyperpolarization-activated cationic current was necessary for the occurrence of the phenomenon. A single-neuron biophysical model including h-current modulation by intracellular calcium was able to display HAGPA. This form of neuronal memory not only allows the transformation of inhibition into an increase of firing rate, but also endows neurons with a mechanism to compute the properties of successive inputs into persistent activity, thus solving a difficult computational problem.</para>
        <para>The original paper reference is cited below:</para>
        <para>
          Hyperpolarization-Activated Graded Persistent Activity in the Prefrontal Cortex, Milena Winograd, Alain Destexhe, Maria V. Sanchez-Vives, 2008,
          <emphasis>Proceedings of the National Academy of the Sciences of the United States of America</emphasis>
          , 105, 7298-7303.
          <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=18474856&amp;dopt=Abstract">PubMed ID: 18474856</ulink>
        </para>
        <informalfigure float="0" id="fig_reaction_diagram">
          <mediaobject>
            <imageobject>
              <objectinfo>
                <title>model diagram</title>
              </objectinfo>
              <imagedata fileref="winograd_2008.png" />
            </imageobject>
          </mediaobject>
          <caption>Schematic diagram of the neuron model.</caption>
        </informalfigure>
      </sect1>
    </article>
  </documentation>
  <units name="coulomb_per_mole">
    <unit units="coulomb" />
    <unit exponent="-1.0" units="mole" />
  </units>
  <units name="mole_liter_per_millimole_per_cm3">
    <unit exponent="-1.0" prefix="milli" units="mole" />
    <unit exponent="-3.0" units="centimeter" />
    <unit units="liter" />
    <unit units="mole" />
  </units>
  <units name="coulomb_per_cm_cubed">
    <unit exponent="-3.0" units="centimeter" />
    <unit units="coulomb" />
  </units>
  <units name="volt_per_millivolt">
    <unit exponent="-1.0" units="millivolt" />
    <unit units="volt" />
  </units>
  <units name="mF_per_cm_squared">
    <unit exponent="-2.0" units="centimeter" />
    <unit prefix="milli" units="farad" />
  </units>
  <units name="per_second">
    <unit exponent="-1.0" units="second" />
  </units>
  <units name="per_second_per_millivolt">
    <unit units="per_second" />
    <unit units="per_millivolt" />
  </units>
  <units name="per_millivolt_per_second">
    <unit units="per_millivolt" />
    <unit units="per_second" />
  </units>
  <units name="fixer">
    <unit exponent="3.0" units="centimeter" />
    <unit exponent="-1.0" prefix="milli" units="ampere" />
    <unit exponent="-1.0" units="mole" />
    <unit exponent="-1.0" units="liter" />
    <unit prefix="milli" units="mole" />
    <unit units="ampere" />
  </units>
  <units name="milliampere_per_millisiemens_per_millivolt">
    <unit prefix="milli" units="ampere" />
    <unit exponent="-1.0" prefix="milli" units="siemens" />
    <unit exponent="-1.0" prefix="milli" units="volt" />
  </units>
  <units name="cm_per_second">
    <unit prefix="centi" units="meter" />
    <unit exponent="-1.0" units="second" />
  </units>
  <units name="millivolts_per_volts">
    <unit exponent="-1.0" units="volt" />
    <unit units="millivolt" />
  </units>
  <units name="millisiemens_per_cm_squared">
    <unit prefix="milli" units="siemens" />
    <unit exponent="-2.0" units="centimeter" />
  </units>
  <units name="millivolt">
    <unit prefix="milli" units="volt" />
  </units>
  <units name="mM">
    <unit exponent="-1.0" units="liter" />
    <unit prefix="milli" units="mole" />
  </units>
  <units name="millicoulomb_per_coulomb">
    <unit prefix="milli" units="coulomb" />
    <unit exponent="-1.0" units="coulomb" />
  </units>
  <units name="per_millivolt">
    <unit exponent="-1.0" units="millivolt" />
  </units>
  <units name="mA_per_cm_squared">
    <unit prefix="milli" units="ampere" />
    <unit exponent="-2.0" units="centimeter" />
  </units>
  <units name="centimeter">
    <unit prefix="centi" units="meter" />
  </units>
  <units name="milliampere_per_cm_squared">
    <unit prefix="milli" units="ampere" />
    <unit exponent="-2.0" units="centimeter" />
  </units>
  <units name="mM_per_second">
    <unit exponent="-1.0" units="second" />
    <unit units="mM" />
  </units>
  <units name="joule_per_mole_per_kelvin">
    <unit exponent="-1.0" units="mole" />
    <unit units="joule" />
    <unit exponent="-1.0" units="kelvin" />
  </units>
  <component name="environment">
    <variable cmeta:id="environment.time" name="time" public_interface="out" units="second" />
  </component>
  <component name="membrane">
    <variable cmeta:id="membrane.V" initial_value="-70" name="V" public_interface="out" units="millivolt" />
    <variable cmeta:id="membrane.V_T" initial_value="-55" name="V_T" public_interface="out" units="millivolt" />
    <variable cmeta:id="membrane.V_S" initial_value="0" name="V_S" public_interface="out" units="millivolt" />
    <variable cmeta:id="membrane.C_m" initial_value="1e-3" name="C_m" units="mF_per_cm_squared" />
    <variable cmeta:id="membrane.F" initial_value="96489" name="F" public_interface="out" units="coulomb_per_mole" />
    <variable cmeta:id="membrane.R" initial_value="8.314" name="R" public_interface="out" units="joule_per_mole_per_kelvin" />
    <variable cmeta:id="membrane.T" initial_value="296.65" name="T" public_interface="out" units="kelvin" />
    <variable cmeta:id="membrane.I_leak" name="I_leak" public_interface="in" units="milliampere_per_cm_squared" />
    <variable cmeta:id="membrane.I_Na" name="I_Na" public_interface="in" units="milliampere_per_cm_squared" />
    <variable cmeta:id="membrane.I_KD" name="I_KD" public_interface="in" units="milliampere_per_cm_squared" />
    <variable cmeta:id="membrane.I_KM" name="I_KM" public_interface="in" units="milliampere_per_cm_squared" />
    <variable cmeta:id="membrane.I_CaL" name="I_CaL" public_interface="in" units="milliampere_per_cm_squared" />
    <variable cmeta:id="membrane.I_h" name="I_h" public_interface="in" units="milliampere_per_cm_squared" />
    <variable cmeta:id="membrane.I_app" name="I_app" public_interface="in" units="milliampere_per_cm_squared" />
    <variable cmeta:id="membrane.time" name="time" public_interface="in" units="second" />
    <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>
            <times />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolts_per_volts">1e-3</cn>
            <apply>
              <minus />
              <apply>
                <minus />
                <apply>
                  <minus />
                  <apply>
                    <minus />
                    <apply>
                      <minus />
                      <apply>
                        <plus />
                        <ci>I_app</ci>
                        <apply>
                          <minus />
                          <ci>I_leak</ci>
                        </apply>
                      </apply>
                      <ci>I_Na</ci>
                    </apply>
                    <ci>I_KD</ci>
                  </apply>
                  <ci>I_KM</ci>
                </apply>
                <ci>I_CaL</ci>
              </apply>
              <ci>I_h</ci>
            </apply>
          </apply>
          <ci>C_m</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="stimulus_protocol">
    <variable cmeta:id="stimulus_protocol.I_app" name="I_app" public_interface="out" units="milliampere_per_cm_squared" />
    <variable cmeta:id="stimulus_protocol.i_stimStart" initial_value="5" name="i_stimStart" units="second" />
    <variable cmeta:id="stimulus_protocol.i_stimEnd" initial_value="9" name="i_stimEnd" units="second" />
    <variable cmeta:id="stimulus_protocol.i_stimAmplitude" initial_value="-0.3" name="i_stimAmplitude" units="milliampere_per_cm_squared" />
    <variable cmeta:id="stimulus_protocol.time" name="time" public_interface="in" units="second" />
    <variable cmeta:id="stimulus_protocol.tau" name="tau" units="second" />
    <variable cmeta:id="stimulus_protocol.period" initial_value="9" name="period" units="second" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>tau</ci>
        <apply>
          <minus />
          <ci>time</ci>
          <apply>
            <times />
            <ci>period</ci>
            <apply>
              <floor />
              <apply>
                <divide />
                <ci>time</ci>
                <ci>period</ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>I_app</ci>
        <piecewise>
          <piece>
            <ci>i_stimAmplitude</ci>
            <apply>
              <and />
              <apply>
                <geq />
                <ci>tau</ci>
                <ci>i_stimStart</ci>
              </apply>
              <apply>
                <leq />
                <ci>tau</ci>
                <ci>i_stimEnd</ci>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="milliampere_per_cm_squared">0</cn>
          </otherwise>
        </piecewise>
      </apply>
    </math>
  </component>
  <component name="I_leak">
    <variable cmeta:id="I_leak.I_leak" name="I_leak" public_interface="out" units="milliampere_per_cm_squared" />
    <variable cmeta:id="I_leak.g_leak" initial_value="1" name="g_leak" units="millisiemens_per_cm_squared" />
    <variable cmeta:id="I_leak.E_leak" initial_value="-70" name="E_leak" units="millivolt" />
    <variable cmeta:id="I_leak.V" name="V" public_interface="in" units="millivolt" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>I_leak</ci>
        <apply>
          <times />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="milliampere_per_millisiemens_per_millivolt">1000</cn>
          <ci>g_leak</ci>
          <apply>
            <minus />
            <ci>V</ci>
            <ci>E_leak</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="I_Na">
    <variable cmeta:id="I_Na.I_Na" name="I_Na" public_interface="out" units="milliampere_per_cm_squared" />
    <variable cmeta:id="I_Na.g_Na" initial_value="70" name="g_Na" units="millisiemens_per_cm_squared" />
    <variable cmeta:id="I_Na.E_Na" initial_value="50" name="E_Na" units="millivolt" />
    <variable cmeta:id="I_Na.V_S" name="V_S" public_interface="in" private_interface="out" units="millivolt" />
    <variable cmeta:id="I_Na.V_T" name="V_T" public_interface="in" private_interface="out" units="millivolt" />
    <variable cmeta:id="I_Na.V" name="V" public_interface="in" private_interface="out" units="millivolt" />
    <variable cmeta:id="I_Na.time" name="time" public_interface="in" private_interface="out" units="second" />
    <variable cmeta:id="I_Na.m" name="m" private_interface="in" units="dimensionless" />
    <variable cmeta:id="I_Na.h" name="h" private_interface="in" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>I_Na</ci>
        <apply>
          <times />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="milliampere_per_millisiemens_per_millivolt">1000</cn>
          <ci>g_Na</ci>
          <apply>
            <power />
            <ci>m</ci>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
          </apply>
          <ci>h</ci>
          <apply>
            <minus />
            <ci>V</ci>
            <ci>E_Na</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="Na_m_gate">
    <variable cmeta:id="Na_m_gate.m" initial_value="0" name="m" public_interface="out" units="dimensionless" />
    <variable cmeta:id="Na_m_gate.V_T" name="V_T" public_interface="in" units="millivolt" />
    <variable cmeta:id="Na_m_gate.V" name="V" public_interface="in" units="millivolt" />
    <variable cmeta:id="Na_m_gate.time" name="time" public_interface="in" units="second" />
    <variable cmeta:id="Na_m_gate.alpha" name="alpha" units="per_second" />
    <variable cmeta:id="Na_m_gate.beta" name="beta" units="per_second" />
    <variable cmeta:id="Na_m_gate.tau_m" name="tau_m" units="second" />
    <variable cmeta:id="Na_m_gate.m_inf" name="m_inf" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>alpha</ci>
        <piecewise>
          <piece>
            <apply>
              <times />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second_per_millivolt">0.32</cn>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <plus />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">13</cn>
                      <ci>V_T</ci>
                    </apply>
                    <ci>V</ci>
                  </apply>
                  <apply>
                    <times />
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt />
              <apply>
                <abs />
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <plus />
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">13</cn>
                      <ci>V_T</ci>
                    </apply>
                    <ci>V</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
                </apply>
              </apply>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1e-6</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide />
              <apply>
                <times />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second_per_millivolt">0.32</cn>
                <apply>
                  <minus />
                  <apply>
                    <plus />
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">13</cn>
                    <ci>V_T</ci>
                  </apply>
                  <ci>V</ci>
                </apply>
              </apply>
              <apply>
                <minus />
                <apply>
                  <exp />
                  <apply>
                    <divide />
                    <apply>
                      <minus />
                      <apply>
                        <plus />
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">13</cn>
                        <ci>V_T</ci>
                      </apply>
                      <ci>V</ci>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
                  </apply>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <ci>beta</ci>
        <piecewise>
          <piece>
            <apply>
              <times />
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second_per_millivolt">0.28</cn>
              </apply>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <minus />
                      <apply>
                        <minus />
                        <ci>V</ci>
                        <ci>V_T</ci>
                      </apply>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                    </apply>
                  </apply>
                  <apply>
                    <times />
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt />
              <apply>
                <abs />
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <minus />
                      <apply>
                        <minus />
                        <ci>V</ci>
                        <ci>V_T</ci>
                      </apply>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                    </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                </apply>
              </apply>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1e-6</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide />
              <apply>
                <times />
                <apply>
                  <minus />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second_per_millivolt">0.28</cn>
                </apply>
                <apply>
                  <minus />
                  <apply>
                    <minus />
                    <ci>V</ci>
                    <ci>V_T</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                </apply>
              </apply>
              <apply>
                <minus />
                <apply>
                  <exp />
                  <apply>
                    <divide />
                    <apply>
                      <minus />
                      <apply>
                        <minus />
                        <apply>
                          <minus />
                          <ci>V</ci>
                          <ci>V_T</ci>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                      </apply>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                  </apply>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <ci>tau_m</ci>
        <apply>
          <divide />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
          <apply>
            <plus />
            <ci>alpha</ci>
            <ci>beta</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>m_inf</ci>
        <apply>
          <divide />
          <ci>alpha</ci>
          <apply>
            <plus />
            <ci>alpha</ci>
            <ci>beta</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>m</ci>
        </apply>
        <apply>
          <divide />
          <apply>
            <minus />
            <ci>m_inf</ci>
            <ci>m</ci>
          </apply>
          <ci>tau_m</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="Na_h_gate">
    <variable cmeta:id="Na_h_gate.h" initial_value="0" name="h" public_interface="out" units="dimensionless" />
    <variable cmeta:id="Na_h_gate.V_S" name="V_S" public_interface="in" units="millivolt" />
    <variable cmeta:id="Na_h_gate.V_T" name="V_T" public_interface="in" units="millivolt" />
    <variable cmeta:id="Na_h_gate.V" name="V" public_interface="in" units="millivolt" />
    <variable cmeta:id="Na_h_gate.time" name="time" public_interface="in" units="second" />
    <variable cmeta:id="Na_h_gate.alpha_h" name="alpha_h" units="per_second" />
    <variable cmeta:id="Na_h_gate.beta_h" name="beta_h" units="per_second" />
    <variable cmeta:id="Na_h_gate.h_inf" name="h_inf" units="dimensionless" />
    <variable cmeta:id="Na_h_gate.tau_h" name="tau_h" units="second" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>alpha_h</ci>
        <apply>
          <times />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">0.128</cn>
          <apply>
            <exp />
            <apply>
              <divide />
              <apply>
                <minus />
                <apply>
                  <plus />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">17</cn>
                  <ci>V_T</ci>
                  <ci>V_S</ci>
                </apply>
                <ci>V</ci>
              </apply>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">18</cn>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>beta_h</ci>
        <apply>
          <divide />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">4</cn>
          <apply>
            <plus />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            <apply>
              <exp />
              <apply>
                <divide />
                <apply>
                  <minus />
                  <apply>
                    <plus />
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                    <ci>V_S</ci>
                    <ci>V_T</ci>
                  </apply>
                  <ci>V</ci>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>tau_h</ci>
        <apply>
          <divide />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
          <apply>
            <plus />
            <ci>alpha_h</ci>
            <ci>beta_h</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>h_inf</ci>
        <apply>
          <divide />
          <ci>alpha_h</ci>
          <apply>
            <plus />
            <ci>alpha_h</ci>
            <ci>beta_h</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>h</ci>
        </apply>
        <apply>
          <divide />
          <apply>
            <minus />
            <ci>h_inf</ci>
            <ci>h</ci>
          </apply>
          <ci>tau_h</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="I_KD">
    <variable cmeta:id="I_KD.I_KD" name="I_KD" public_interface="out" units="milliampere_per_cm_squared" />
    <variable cmeta:id="I_KD.g_KD" initial_value="7" name="g_KD" units="millisiemens_per_cm_squared" />
    <variable cmeta:id="I_KD.E_K" initial_value="-95" name="E_K" public_interface="out" units="millivolt" />
    <variable cmeta:id="I_KD.V_T" name="V_T" public_interface="in" private_interface="out" units="millivolt" />
    <variable cmeta:id="I_KD.V" name="V" public_interface="in" private_interface="out" units="millivolt" />
    <variable cmeta:id="I_KD.time" name="time" public_interface="in" private_interface="out" units="second" />
    <variable cmeta:id="I_KD.n" name="n" private_interface="in" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>I_KD</ci>
        <apply>
          <times />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="milliampere_per_millisiemens_per_millivolt">1000</cn>
          <ci>g_KD</ci>
          <apply>
            <power />
            <ci>n</ci>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">4</cn>
          </apply>
          <apply>
            <minus />
            <ci>V</ci>
            <ci>E_K</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="KD_n_gate">
    <variable cmeta:id="KD_n_gate.n" initial_value="0" name="n" public_interface="out" units="dimensionless" />
    <variable cmeta:id="KD_n_gate.V_T" name="V_T" public_interface="in" units="millivolt" />
    <variable cmeta:id="KD_n_gate.V" name="V" public_interface="in" units="millivolt" />
    <variable cmeta:id="KD_n_gate.time" name="time" public_interface="in" units="second" />
    <variable cmeta:id="KD_n_gate.alpha_n" name="alpha_n" units="per_second" />
    <variable cmeta:id="KD_n_gate.beta_n" name="beta_n" units="per_second" />
    <variable cmeta:id="KD_n_gate.tau_n" name="tau_n" units="second" />
    <variable cmeta:id="KD_n_gate.n_inf" name="n_inf" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>alpha_n</ci>
        <piecewise>
          <piece>
            <apply>
              <times />
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second_per_millivolt">0.032</cn>
              </apply>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <minus />
                      <ci>V</ci>
                      <ci>V_T</ci>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">15</cn>
                  </apply>
                  <apply>
                    <times />
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt />
              <apply>
                <abs />
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <minus />
                      <ci>V</ci>
                      <ci>V_T</ci>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">15</cn>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                </apply>
              </apply>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1e-6</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide />
              <apply>
                <times />
                <apply>
                  <minus />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second_per_millivolt">0.032</cn>
                </apply>
                <apply>
                  <minus />
                  <apply>
                    <minus />
                    <ci>V</ci>
                    <ci>V_T</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">15</cn>
                </apply>
              </apply>
              <apply>
                <minus />
                <apply>
                  <exp />
                  <apply>
                    <divide />
                    <apply>
                      <minus />
                      <apply>
                        <minus />
                        <ci>V</ci>
                        <ci>V_T</ci>
                      </apply>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">15</cn>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                  </apply>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <ci>beta_n</ci>
        <piecewise>
          <piece>
            <apply>
              <times />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second_per_millivolt">0.5</cn>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
              <apply>
                <plus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <minus />
                      <ci>V</ci>
                      <ci>V_T</ci>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                  </apply>
                  <apply>
                    <times />
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt />
              <apply>
                <abs />
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <minus />
                      <apply>
                        <minus />
                        <ci>V</ci>
                        <ci>V_T</ci>
                      </apply>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                    </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                </apply>
              </apply>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1e-6</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide />
              <apply>
                <times />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second_per_millivolt">0.5</cn>
                <apply>
                  <minus />
                  <apply>
                    <minus />
                    <apply>
                      <minus />
                      <ci>V</ci>
                      <ci>V_T</ci>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                  </apply>
                </apply>
              </apply>
              <apply>
                <minus />
                <apply>
                  <exp />
                  <apply>
                    <divide />
                    <apply>
                      <minus />
                      <apply>
                        <minus />
                        <apply>
                          <minus />
                          <ci>V</ci>
                          <ci>V_T</ci>
                        </apply>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                      </apply>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">40</cn>
                  </apply>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <ci>tau_n</ci>
        <apply>
          <divide />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
          <apply>
            <plus />
            <ci>alpha_n</ci>
            <ci>beta_n</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>n_inf</ci>
        <apply>
          <divide />
          <ci>alpha_n</ci>
          <apply>
            <plus />
            <ci>alpha_n</ci>
            <ci>beta_n</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>n</ci>
        </apply>
        <apply>
          <divide />
          <apply>
            <minus />
            <ci>n_inf</ci>
            <ci>n</ci>
          </apply>
          <ci>tau_n</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="I_KM">
    <variable cmeta:id="I_KM.I_KM" name="I_KM" public_interface="out" units="milliampere_per_cm_squared" />
    <variable cmeta:id="I_KM.g_KM" initial_value="0.004" name="g_KM" units="millisiemens_per_cm_squared" />
    <variable cmeta:id="I_KM.E_K" name="E_K" public_interface="in" units="millivolt" />
    <variable cmeta:id="I_KM.V" name="V" public_interface="in" private_interface="out" units="millivolt" />
    <variable cmeta:id="I_KM.time" name="time" public_interface="in" private_interface="out" units="second" />
    <variable cmeta:id="I_KM.p" name="p" private_interface="in" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>I_KM</ci>
        <apply>
          <times />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="milliampere_per_millisiemens_per_millivolt">1000</cn>
          <ci>g_KM</ci>
          <ci>p</ci>
          <apply>
            <minus />
            <ci>V</ci>
            <ci>E_K</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="KM_p_gate">
    <variable cmeta:id="KM_p_gate.p" initial_value="0" name="p" public_interface="out" units="dimensionless" />
    <variable cmeta:id="KM_p_gate.V" name="V" public_interface="in" units="millivolt" />
    <variable cmeta:id="KM_p_gate.time" name="time" public_interface="in" units="second" />
    <variable cmeta:id="KM_p_gate.p_inf" name="p_inf" units="dimensionless" />
    <variable cmeta:id="KM_p_gate.tau_p" name="tau_p" units="second" />
    <variable cmeta:id="KM_p_gate.tau_max" initial_value="4" name="tau_max" units="second" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>p_inf</ci>
        <piecewise>
          <piece>
            <apply>
              <divide />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
              <apply>
                <plus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <apply>
                  <exp />
                  <apply>
                    <divide />
                    <apply>
                      <minus />
                      <apply>
                        <plus />
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">35</cn>
                      </apply>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and />
              <apply>
                <lt />
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <plus />
                      <ci>V</ci>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">35</cn>
                    </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">25</cn>
              </apply>
              <apply>
                <gt />
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <plus />
                      <ci>V</ci>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">35</cn>
                    </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                </apply>
                <apply>
                  <minus />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">25</cn>
                </apply>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <ci>tau_p</ci>
        <piecewise>
          <piece>
            <apply>
              <divide />
              <ci>tau_max</ci>
              <apply>
                <plus />
                <apply>
                  <times />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3.3</cn>
                  <apply>
                    <exp />
                    <apply>
                      <divide />
                      <apply>
                        <plus />
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">35</cn>
                      </apply>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">20</cn>
                    </apply>
                  </apply>
                </apply>
                <apply>
                  <exp />
                  <apply>
                    <divide />
                    <apply>
                      <minus />
                      <apply>
                        <plus />
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">35</cn>
                      </apply>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">20</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <and />
              <apply>
                <lt />
                <apply>
                  <divide />
                  <apply>
                    <plus />
                    <ci>V</ci>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">35</cn>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">20</cn>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">25</cn>
              </apply>
              <apply>
                <gt />
                <apply>
                  <divide />
                  <apply>
                    <plus />
                    <ci>V</ci>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">35</cn>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">20</cn>
                </apply>
                <apply>
                  <minus />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">25</cn>
                </apply>
              </apply>
            </apply>
          </piece>
          <otherwise>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>p</ci>
        </apply>
        <apply>
          <divide />
          <apply>
            <minus />
            <ci>p_inf</ci>
            <ci>p</ci>
          </apply>
          <ci>tau_p</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="I_CaL">
    <variable cmeta:id="I_CaL.I_CaL" name="I_CaL" public_interface="out" units="milliampere_per_cm_squared" />
    <variable cmeta:id="I_CaL.P_Ca" initial_value="2.76e-4" name="P_Ca" units="cm_per_second" />
    <variable cmeta:id="I_CaL.V" name="V" public_interface="in" private_interface="out" units="millivolt" />
    <variable cmeta:id="I_CaL.time" name="time" public_interface="in" private_interface="out" units="second" />
    <variable cmeta:id="I_CaL.G" name="G" public_interface="in" units="coulomb_per_cm_cubed" />
    <variable cmeta:id="I_CaL.q" name="q" private_interface="in" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>I_CaL</ci>
        <apply>
          <times />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millicoulomb_per_coulomb">1000</cn>
          <ci>P_Ca</ci>
          <apply>
            <power />
            <ci>q</ci>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
          </apply>
          <ci>G</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="CaL_q_gate">
    <variable cmeta:id="CaL_q_gate.q" initial_value="0.00247262" name="q" public_interface="out" units="dimensionless" />
    <variable cmeta:id="CaL_q_gate.V" name="V" public_interface="in" units="millivolt" />
    <variable cmeta:id="CaL_q_gate.time" name="time" public_interface="in" units="second" />
    <variable cmeta:id="CaL_q_gate.alpha_q" name="alpha_q" units="per_second" />
    <variable cmeta:id="CaL_q_gate.beta_q" name="beta_q" units="per_second" />
    <variable cmeta:id="CaL_q_gate.tau_q" name="tau_q" units="second" />
    <variable cmeta:id="CaL_q_gate.q_inf" name="q_inf" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>alpha_q</ci>
        <apply>
          <divide />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">6.32</cn>
          <apply>
            <plus />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            <apply>
              <exp />
              <apply>
                <divide />
                <apply>
                  <minus />
                  <apply>
                    <minus />
                    <ci>V</ci>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5</cn>
                  </apply>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">13.89</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>beta_q</ci>
        <piecewise>
          <piece>
            <apply>
              <times />
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second">0.02</cn>
              <apply>
                <plus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">5.36</cn>
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">1.31</cn>
                    <ci>V</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">2</cn>
                </apply>
              </apply>
            </apply>
            <apply>
              <lt />
              <apply>
                <abs />
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">1.31</cn>
                    <ci>V</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5.36</cn>
                </apply>
              </apply>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1e-6</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <divide />
              <apply>
                <times />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_second_per_millivolt">0.02</cn>
                <apply>
                  <minus />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">1.31</cn>
                  <ci>V</ci>
                </apply>
              </apply>
              <apply>
                <minus />
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                <apply>
                  <exp />
                  <apply>
                    <divide />
                    <apply>
                      <minus />
                      <ci>V</ci>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">1.31</cn>
                    </apply>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5.36</cn>
                  </apply>
                </apply>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
      <apply>
        <eq />
        <ci>tau_q</ci>
        <apply>
          <divide />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
          <apply>
            <plus />
            <ci>alpha_q</ci>
            <ci>beta_q</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>q_inf</ci>
        <apply>
          <divide />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
          <apply>
            <plus />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            <apply>
              <exp />
              <apply>
                <divide />
                <apply>
                  <plus />
                  <ci>V</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                </apply>
                <apply>
                  <minus />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">10</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>q</ci>
        </apply>
        <apply>
          <divide />
          <apply>
            <minus />
            <ci>q_inf</ci>
            <ci>q</ci>
          </apply>
          <ci>tau_q</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="G_nonlin">
    <variable cmeta:id="G_nonlin.G" name="G" public_interface="out" units="coulomb_per_cm_cubed" />
    <variable cmeta:id="G_nonlin.Z" initial_value="2" name="Z" units="dimensionless" />
    <variable cmeta:id="G_nonlin.Ca_o" initial_value="2" name="Ca_o" units="mM" />
    <variable cmeta:id="G_nonlin.V" name="V" public_interface="in" units="millivolt" />
    <variable cmeta:id="G_nonlin.time" name="time" public_interface="in" private_interface="out" units="second" />
    <variable cmeta:id="G_nonlin.F" name="F" public_interface="in" private_interface="out" units="coulomb_per_mole" />
    <variable cmeta:id="G_nonlin.R" name="R" public_interface="in" private_interface="out" units="joule_per_mole_per_kelvin" />
    <variable cmeta:id="G_nonlin.T" name="T" public_interface="in" private_interface="out" units="kelvin" />
    <variable cmeta:id="G_nonlin.I_CaL" name="I_CaL" public_interface="in" private_interface="out" units="milliampere_per_cm_squared" />
    <variable cmeta:id="G_nonlin.Ca_i" name="Ca_i" public_interface="out" private_interface="in" units="mM" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>G</ci>
        <apply>
          <divide />
          <apply>
            <times />
            <apply>
              <divide />
              <apply>
                <times />
                <apply>
                  <power />
                  <ci>Z</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                </apply>
                <apply>
                  <power />
                  <ci>F</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="volt_per_millivolt">1e-3</cn>
                <ci>V</ci>
              </apply>
              <apply>
                <times />
                <ci>R</ci>
                <ci>T</ci>
              </apply>
            </apply>
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mole_liter_per_millimole_per_cm3">1e-6</cn>
            <apply>
              <minus />
              <ci>Ca_i</ci>
              <apply>
                <times />
                <ci>Ca_o</ci>
                <apply>
                  <exp />
                  <apply>
                    <divide />
                    <apply>
                      <times />
                      <ci>Z</ci>
                      <ci>F</ci>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="volt_per_millivolt">1e-3</cn>
                      <ci>V</ci>
                    </apply>
                    <apply>
                      <times />
                      <ci>R</ci>
                      <ci>T</ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            <apply>
              <exp />
              <apply>
                <divide />
                <apply>
                  <times />
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="volt_per_millivolt">1e-3</cn>
                  <ci>Z</ci>
                  <ci>F</ci>
                  <ci>V</ci>
                </apply>
                <apply>
                  <times />
                  <ci>R</ci>
                  <ci>T</ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="dCa_i_dt">
    <variable cmeta:id="dCa_i_dt.Ca_i" initial_value="100e-6" name="Ca_i" public_interface="out" units="mM" />
    <variable cmeta:id="dCa_i_dt.Ca_inf" initial_value="100e-6" name="Ca_inf" units="mM" />
    <variable cmeta:id="dCa_i_dt.tau_r" initial_value="17e-3" name="tau_r" units="second" />
    <variable cmeta:id="dCa_i_dt.d" initial_value="1e-4" name="d" units="centimeter" />
    <variable cmeta:id="dCa_i_dt.time" name="time" public_interface="in" units="second" />
    <variable cmeta:id="dCa_i_dt.drive_channel" name="drive_channel" units="mM_per_second" />
    <variable cmeta:id="dCa_i_dt.k" initial_value="0.1" name="k" units="fixer" />
    <variable cmeta:id="dCa_i_dt.I_CaL" name="I_CaL" public_interface="in" units="milliampere_per_cm_squared" />
    <variable cmeta:id="dCa_i_dt.F" name="F" public_interface="in" units="coulomb_per_mole" />
    <variable cmeta:id="dCa_i_dt.R" name="R" public_interface="in" units="joule_per_mole_per_kelvin" />
    <variable cmeta:id="dCa_i_dt.T" name="T" public_interface="in" units="kelvin" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>drive_channel</ci>
        <apply>
          <divide />
          <apply>
            <times />
            <ci>k</ci>
            <ci>I_CaL</ci>
          </apply>
          <apply>
            <times />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
            <ci>F</ci>
            <ci>d</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <apply>
          <diff />
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Ca_i</ci>
        </apply>
        <piecewise>
          <piece>
            <apply>
              <divide />
              <apply>
                <minus />
                <ci>Ca_inf</ci>
                <ci>Ca_i</ci>
              </apply>
              <ci>tau_r</ci>
            </apply>
            <apply>
              <leq />
              <ci>drive_channel</ci>
              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="mM_per_second">0</cn>
            </apply>
          </piece>
          <otherwise>
            <apply>
              <plus />
              <ci>drive_channel</ci>
              <apply>
                <divide />
                <apply>
                  <minus />
                  <ci>Ca_inf</ci>
                  <ci>Ca_i</ci>
                </apply>
                <ci>tau_r</ci>
              </apply>
            </apply>
          </otherwise>
        </piecewise>
      </apply>
    </math>
  </component>
  <component name="I_h">
    <variable cmeta:id="I_h.I_h" name="I_h" public_interface="out" units="milliampere_per_cm_squared" />
    <variable cmeta:id="I_h.m" name="m" units="dimensionless" />
    <variable cmeta:id="I_h.E_h" initial_value="-20" name="E_h" units="millivolt" />
    <variable cmeta:id="I_h.g_hbar" initial_value="0.02" name="g_hbar" units="millisiemens_per_cm_squared" />
    <variable cmeta:id="I_h.cac" initial_value="0.006" name="cac" units="mM" />
    <variable cmeta:id="I_h.V" name="V" public_interface="in" private_interface="out" units="millivolt" />
    <variable cmeta:id="I_h.V_S" initial_value="0" name="V_S" private_interface="out" units="millivolt" />
    <variable cmeta:id="I_h.Ca_i" name="Ca_i" public_interface="in" private_interface="out" units="mM" />
    <variable cmeta:id="I_h.o_1" name="o_1" private_interface="in" units="dimensionless" />
    <variable cmeta:id="I_h.o_2" name="o_2" private_interface="in" units="dimensionless" />
    <variable cmeta:id="I_h.g_inc" initial_value="2" name="g_inc" units="dimensionless" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>m</ci>
        <apply>
          <plus />
          <ci>o_1</ci>
          <ci>g_inc</ci>
          <ci>o_2</ci>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>I_h</ci>
        <apply>
          <times />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="milliampere_per_millisiemens_per_millivolt">1000</cn>
          <ci>g_hbar</ci>
          <ci>m</ci>
          <apply>
            <minus />
            <ci>V</ci>
            <ci>E_h</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <component name="kinetic">
    <variable cmeta:id="kinetic.p_0" name="p_0" units="dimensionless" />
    <variable cmeta:id="kinetic.p_1" name="p_1" public_interface="out" units="dimensionless" />
    <variable cmeta:id="kinetic.c_1" name="c_1" units="dimensionless" />
    <variable cmeta:id="kinetic.o_1" name="o_1" public_interface="out" units="dimensionless" />
    <variable cmeta:id="kinetic.o_2" name="o_2" public_interface="out" units="dimensionless" />
    <variable cmeta:id="kinetic.alpha" name="alpha" public_interface="in" units="dimensionless" />
    <variable cmeta:id="kinetic.beta" name="beta" public_interface="in" units="dimensionless" />
    <variable cmeta:id="kinetic.k_1Ca" name="k_1Ca" public_interface="in" units="per_second" />
    <variable cmeta:id="kinetic.k_2" name="k_2" public_interface="in" units="per_second" />
    <variable cmeta:id="kinetic.k_3p" name="k_3p" public_interface="in" units="per_second" />
    <variable cmeta:id="kinetic.k_4" name="k_4" public_interface="in" units="per_second" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>p_0</ci>
        <apply>
          <divide />
          <apply>
            <times />
            <ci>p_1</ci>
            <ci>k_2</ci>
          </apply>
          <ci>k_1Ca</ci>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>p_1</ci>
        <apply>
          <minus />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
          <ci>p_0</ci>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>c_1</ci>
        <apply>
          <times />
          <apply>
            <divide />
            <ci>beta</ci>
            <ci>alpha</ci>
          </apply>
          <ci>o_1</ci>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>o_1</ci>
        <apply>
          <times />
          <apply>
            <divide />
            <ci>k_4</ci>
            <ci>k_3p</ci>
          </apply>
          <ci>o_2</ci>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>o_2</ci>
        <apply>
          <minus />
          <apply>
            <minus />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            <ci>c_1</ci>
          </apply>
          <ci>o_1</ci>
        </apply>
      </apply>
    </math>
  </component>
  <component name="rate_constants">
    <variable cmeta:id="rate_constants.h_inf" name="h_inf" units="second" />
    <variable cmeta:id="rate_constants.tau_s" name="tau_s" units="second" />
    <variable cmeta:id="rate_constants.alpha" name="alpha" public_interface="out" units="dimensionless" />
    <variable cmeta:id="rate_constants.beta" name="beta" public_interface="out" units="dimensionless" />
    <variable cmeta:id="rate_constants.k_1Ca" name="k_1Ca" public_interface="out" units="per_second" />
    <variable cmeta:id="rate_constants.k_3p" name="k_3p" public_interface="out" units="per_second" />
    <variable cmeta:id="rate_constants.V" name="V" public_interface="in" units="millivolt" />
    <variable cmeta:id="rate_constants.V_S" name="V_S" public_interface="in" units="millivolt" />
    <variable cmeta:id="rate_constants.p_1" name="p_1" public_interface="in" units="dimensionless" />
    <variable cmeta:id="rate_constants.k_2" initial_value="0.1" name="k_2" public_interface="out" units="per_second" />
    <variable cmeta:id="rate_constants.P_c" initial_value="0.01" name="P_c" units="dimensionless" />
    <variable cmeta:id="rate_constants.k_4" initial_value="1" name="k_4" public_interface="out" units="per_second" />
    <variable cmeta:id="rate_constants.n_Ca" initial_value="4" name="n_Ca" units="dimensionless" />
    <variable cmeta:id="rate_constants.n_exp" initial_value="1" name="n_exp" units="dimensionless" />
    <variable cmeta:id="rate_constants.p_C" initial_value="0.01" name="p_C" units="dimensionless" />
    <variable cmeta:id="rate_constants.Ca_c" initial_value="0.006" name="Ca_c" units="mM" />
    <variable cmeta:id="rate_constants.Ca_i" name="Ca_i" public_interface="in" units="mM" />
    <variable cmeta:id="rate_constants.tau_m" initial_value="20e-3" name="tau_m" units="second" />
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq />
        <ci>h_inf</ci>
        <apply>
          <divide />
          <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
          <apply>
            <plus />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
            <apply>
              <exp />
              <apply>
                <divide />
                <apply>
                  <minus />
                  <apply>
                    <plus />
                    <ci>V</ci>
                    <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">75</cn>
                  </apply>
                  <ci>V_S</ci>
                </apply>
                <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">5.5</cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>tau_s</ci>
        <apply>
          <plus />
          <ci>tau_m</ci>
          <apply>
            <divide />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1000</cn>
            <apply>
              <plus />
              <apply>
                <exp />
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <plus />
                      <ci>V</ci>
                      <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">71.55</cn>
                    </apply>
                    <ci>V_S</ci>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">14.2</cn>
                </apply>
              </apply>
              <apply>
                <exp />
                <apply>
                  <divide />
                  <apply>
                    <minus />
                    <apply>
                      <minus />
                      <apply>
                        <plus />
                        <ci>V</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">89</cn>
                      </apply>
                      <ci>V_S</ci>
                    </apply>
                  </apply>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">11.6</cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>alpha</ci>
        <apply>
          <divide />
          <ci>h_inf</ci>
          <ci>tau_s</ci>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>beta</ci>
        <apply>
          <divide />
          <apply>
            <minus />
            <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
            <ci>h_inf</ci>
          </apply>
          <ci>tau_s</ci>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>k_1Ca</ci>
        <apply>
          <times />
          <ci>k_2</ci>
          <apply>
            <power />
            <apply>
              <divide />
              <ci>Ca_i</ci>
              <ci>Ca_c</ci>
            </apply>
            <ci>n_Ca</ci>
          </apply>
        </apply>
      </apply>
      <apply>
        <eq />
        <ci>k_3p</ci>
        <apply>
          <times />
          <ci>k_4</ci>
          <apply>
            <power />
            <apply>
              <divide />
              <ci>p_1</ci>
              <ci>p_C</ci>
            </apply>
            <ci>n_exp</ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  <group>
    <relationship_ref relationship="encapsulation" />
    <component_ref component="I_h">
      <component_ref component="rate_constants" />
      <component_ref component="kinetic" />
    </component_ref>
    <component_ref component="I_KM">
      <component_ref component="KM_p_gate" />
    </component_ref>
    <component_ref component="G_nonlin">
      <component_ref component="dCa_i_dt" />
    </component_ref>
    <component_ref component="I_Na">
      <component_ref component="Na_h_gate" />
      <component_ref component="Na_m_gate" />
    </component_ref>
    <component_ref component="I_CaL">
      <component_ref component="CaL_q_gate" />
    </component_ref>
    <component_ref component="I_KD">
      <component_ref component="KD_n_gate" />
    </component_ref>
  </group>
  <connection>
    <map_components component_1="environment" component_2="membrane" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="I_CaL" component_2="G_nonlin" />
    <map_variables variable_1="I_CaL" variable_2="I_CaL" />
  </connection>
  <connection>
    <map_components component_1="Na_m_gate" component_2="I_Na" />
    <map_variables variable_1="m" variable_2="m" />
  </connection>
  <connection>
    <map_components component_1="G_nonlin" component_2="I_CaL" />
    <map_variables variable_1="G" variable_2="G" />
  </connection>
  <connection>
    <map_components component_1="stimulus_protocol" component_2="membrane" />
    <map_variables variable_1="I_app" variable_2="I_app" />
  </connection>
  <connection>
    <map_components component_1="membrane" component_2="I_KM" />
    <map_variables variable_1="V" variable_2="V" />
  </connection>
  <connection>
    <map_components component_1="I_Na" component_2="Na_h_gate" />
    <map_variables variable_1="V_S" variable_2="V_S" />
    <map_variables variable_1="V" variable_2="V" />
    <map_variables variable_1="time" variable_2="time" />
    <map_variables variable_1="V_T" variable_2="V_T" />
  </connection>
  <connection>
    <map_components component_1="I_h" component_2="membrane" />
    <map_variables variable_1="I_h" variable_2="I_h" />
  </connection>
  <connection>
    <map_components component_1="I_CaL" component_2="membrane" />
    <map_variables variable_1="I_CaL" variable_2="I_CaL" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="I_Na" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="I_CaL" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="dCa_i_dt" component_2="G_nonlin" />
    <map_variables variable_1="Ca_i" variable_2="Ca_i" />
  </connection>
  <connection>
    <map_components component_1="membrane" component_2="I_h" />
    <map_variables variable_1="V" variable_2="V" />
  </connection>
  <connection>
    <map_components component_1="I_KD" component_2="KD_n_gate" />
    <map_variables variable_1="V_T" variable_2="V_T" />
    <map_variables variable_1="V" variable_2="V" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="I_h" component_2="rate_constants" />
    <map_variables variable_1="V_S" variable_2="V_S" />
    <map_variables variable_1="V" variable_2="V" />
    <map_variables variable_1="Ca_i" variable_2="Ca_i" />
  </connection>
  <connection>
    <map_components component_1="kinetic" component_2="I_h" />
    <map_variables variable_1="o_1" variable_2="o_1" />
    <map_variables variable_1="o_2" variable_2="o_2" />
  </connection>
  <connection>
    <map_components component_1="CaL_q_gate" component_2="I_CaL" />
    <map_variables variable_1="q" variable_2="q" />
  </connection>
  <connection>
    <map_components component_1="I_KD" component_2="membrane" />
    <map_variables variable_1="I_KD" variable_2="I_KD" />
  </connection>
  <connection>
    <map_components component_1="I_KM" component_2="KM_p_gate" />
    <map_variables variable_1="time" variable_2="time" />
    <map_variables variable_1="V" variable_2="V" />
  </connection>
  <connection>
    <map_components component_1="G_nonlin" component_2="dCa_i_dt" />
    <map_variables variable_1="R" variable_2="R" />
    <map_variables variable_1="T" variable_2="T" />
    <map_variables variable_1="time" variable_2="time" />
    <map_variables variable_1="I_CaL" variable_2="I_CaL" />
    <map_variables variable_1="F" variable_2="F" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="I_KM" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="G_nonlin" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="KM_p_gate" component_2="I_KM" />
    <map_variables variable_1="p" variable_2="p" />
  </connection>
  <connection>
    <map_components component_1="G_nonlin" component_2="I_h" />
    <map_variables variable_1="Ca_i" variable_2="Ca_i" />
  </connection>
  <connection>
    <map_components component_1="I_KM" component_2="membrane" />
    <map_variables variable_1="I_KM" variable_2="I_KM" />
  </connection>
  <connection>
    <map_components component_1="I_leak" component_2="membrane" />
    <map_variables variable_1="I_leak" variable_2="I_leak" />
  </connection>
  <connection>
    <map_components component_1="Na_h_gate" component_2="I_Na" />
    <map_variables variable_1="h" variable_2="h" />
  </connection>
  <connection>
    <map_components component_1="rate_constants" component_2="kinetic" />
    <map_variables variable_1="k_1Ca" variable_2="k_1Ca" />
    <map_variables variable_1="alpha" variable_2="alpha" />
    <map_variables variable_1="beta" variable_2="beta" />
    <map_variables variable_1="k_4" variable_2="k_4" />
    <map_variables variable_1="k_3p" variable_2="k_3p" />
    <map_variables variable_1="k_2" variable_2="k_2" />
  </connection>
  <connection>
    <map_components component_1="membrane" component_2="I_leak" />
    <map_variables variable_1="V" variable_2="V" />
  </connection>
  <connection>
    <map_components component_1="KD_n_gate" component_2="I_KD" />
    <map_variables variable_1="n" variable_2="n" />
  </connection>
  <connection>
    <map_components component_1="membrane" component_2="I_Na" />
    <map_variables variable_1="V" variable_2="V" />
    <map_variables variable_1="V_T" variable_2="V_T" />
    <map_variables variable_1="V_S" variable_2="V_S" />
  </connection>
  <connection>
    <map_components component_1="membrane" component_2="I_KD" />
    <map_variables variable_1="V" variable_2="V" />
    <map_variables variable_1="V_T" variable_2="V_T" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="I_KD" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="I_CaL" component_2="CaL_q_gate" />
    <map_variables variable_1="V" variable_2="V" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="membrane" component_2="I_CaL" />
    <map_variables variable_1="V" variable_2="V" />
  </connection>
  <connection>
    <map_components component_1="environment" component_2="stimulus_protocol" />
    <map_variables variable_1="time" variable_2="time" />
  </connection>
  <connection>
    <map_components component_1="I_Na" component_2="Na_m_gate" />
    <map_variables variable_1="V" variable_2="V" />
    <map_variables variable_1="time" variable_2="time" />
    <map_variables variable_1="V_T" variable_2="V_T" />
  </connection>
  <connection>
    <map_components component_1="I_Na" component_2="membrane" />
    <map_variables variable_1="I_Na" variable_2="I_Na" />
  </connection>
  <connection>
    <map_components component_1="membrane" component_2="G_nonlin" />
    <map_variables variable_1="V" variable_2="V" />
    <map_variables variable_1="T" variable_2="T" />
    <map_variables variable_1="F" variable_2="F" />
    <map_variables variable_1="R" variable_2="R" />
  </connection>
  <connection>
    <map_components component_1="I_KD" component_2="I_KM" />
    <map_variables variable_1="E_K" variable_2="E_K" />
  </connection>
  <connection>
    <map_components component_1="kinetic" component_2="rate_constants" />
    <map_variables variable_1="p_1" variable_2="p_1" />
  </connection>
  <rdf:RDF xmlns:bqmodel="http://biomodels.net/model-qualifiers/" xmlns:bqbiol="http://biomodels.net/biology-qualifiers/" xmlns:ro="http://www.obofoundry.org/ro/ro.owl#" xmlns:opb="http://bhi.washington.edu/OPB#">
    <rdf:Description rdf:about="#rate_constants.k_4">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.tau_r">
      <dcterms:description>calcium removal time constant</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KM.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#stimulus_protocol.period">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.I_Na">
      <dcterms:description>sodium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_Na.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.F">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00089" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KD.n">
      <dcterms:description>delayed rectifier potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_h.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.I_leak">
      <dcterms:description>leakage current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.p_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in second messenger-bound P1 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_Na.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_0">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#stimulus_protocol.I_app">
      <dcterms:description>externally applied stimulus current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_CaL.q">
      <dcterms:description>L-type calcium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.I_h">
      <dcterms:description>hyperpolarisation-activated cationic current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.h_inf">
      <dcterms:description>steady state function of sodium current inactivation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_leak.E_leak">
      <dcterms:description>leakage reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.Ca_i">
      <bqbiol:isPropertyOf>
        <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#entity_2">
          <bqbiol:is rdf:resource="https://identifiers.org/chebi/CHEBI:29108" />
        </rdf:Description>
      </bqbiol:isPropertyOf>
      <dcterms:description>concentration of intracelullar calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.F">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00089" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_Na.m">
      <dcterms:description>sodium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_Na.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KM.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#CaL_q_gate.beta_q">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_h.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KM.E_K">
      <dcterms:description>leakage reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_Na.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.p_0">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in second messenger-bound P0 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.alpha">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#stimulus_protocol.i_stimAmplitude">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00318" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KM_p_gate.tau_p">
      <dcterms:description>activation time constant of slow potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.C_m">
      <dcterms:description>membrane capacitance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_Na.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_KD.E_K">
      <dcterms:description>leakage reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.R">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00410" />
    </rdf:Description>
    <rdf:Description rdf:about="#KD_n_gate.beta_n">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.Ca_o">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_22">
          <semsim:physicalPropertyOf>
            <rdf:Description rdf:about="#entity_0">
              <ro:part_of>
                <rdf:Description rdf:about="#entity_1">
                  <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/fma/FMA:70022" />
                </rdf:Description>
              </ro:part_of>
              <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/chebi/CHEBI:29108" />
            </rdf:Description>
          </semsim:physicalPropertyOf>
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>concentration of extracelullar calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_h.E_h">
      <dcterms:description>hyperpolarisation-activated cationic current reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KM_p_gate.p">
      <dcterms:description>slow potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.o_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in open O1 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_h.cac">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_30">
          <semsim:physicalPropertyOf>
            <rdf:Description rdf:about="#entity_2">
              <ro:part_of>
                <rdf:Description rdf:about="#entity_3">
                  <ro:part_of>
                    <rdf:Description rdf:about="#entity_4">
                      <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/fma/FMA:224850" />
                    </rdf:Description>
                  </ro:part_of>
                  <semsim:hasPhysicalDefinition rdf:resource="http://identifiers.org/fma/FMA:226054" />
                </rdf:Description>
              </ro:part_of>
              <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/chebi/CHEBI:29108" />
            </rdf:Description>
          </semsim:physicalPropertyOf>
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>calcium half-activation constant</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.beta">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_h.Ca_i">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_32">
          <semsim:physicalPropertyOf rdf:resource="#entity_2" />
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>concentration of intracelullar calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.Ca_o">
      <bqbiol:isPropertyOf>
        <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#entity_0">
          <bqbiol:is rdf:resource="https://identifiers.org/chebi/CHEBI:29108" />
        </rdf:Description>
      </bqbiol:isPropertyOf>
      <dcterms:description>concentration of extracelullar calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.k_2">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#CaL_q_gate.q_inf">
      <dcterms:description>steady state function of L-type calcium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.T">
      <dcterms:description>absolute temperature</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.T">
      <dcterms:description>absolute temperature</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_leak.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.Ca_i">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_35">
          <semsim:physicalPropertyOf rdf:resource="#entity_2" />
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>concentration of intracelullar calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KD_n_gate.tau_n">
      <dcterms:description>activation time constant of delayed rectifier potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_m_gate.m">
      <dcterms:description>sodium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.I_CaL">
      <dcterms:description>L-type calcium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.alpha">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KD_n_gate.alpha_n">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_leak.g_leak">
      <dcterms:description>leakage current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_h.Ca_i">
      <bqbiol:isPropertyOf rdf:resource="winograd_destexhe_sanchezvives_2008.cellml#entity_2" />
      <dcterms:description>concentration of intracelullar calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_m_gate.m_inf">
      <dcterms:description>steady state function of sodium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_h.g_hbar">
      <dcterms:description>hyperpolarisation-activated cationic current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.I_CaL">
      <dcterms:description>L-type calcium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.R">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00410" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_h_gate.h_inf">
      <dcterms:description>steady state function of sodium current inactivation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.d">
      <dcterms:description>shell depth beneath membrane</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#CaL_q_gate.beta_q">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_leak.g_leak">
      <dcterms:description>leakage current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KD_n_gate.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_15">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.alpha">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#entity_1">
      <bqbiol:is rdf:resource="http://identifiers.org/fma/FMA:70022" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_leak.E_leak">
      <dcterms:description>leakage reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_h_gate.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#winograd_destexhe_sanchez_vives_2008" bqmodel:isDescribedBy="18474856" semsim:modelName="winograd_destexhe_sanchez_vives_2008">
      <semsim:modelName>winograd_destexhe_sanchezvives_2008</semsim:modelName>
    </rdf:Description>
    <rdf:Description rdf:about="#stimulus_protocol.i_stimAmplitude">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00318" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#CaL_q_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KM_p_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#Na_m_gate.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_11">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.T">
      <dcterms:description>absolute temperature</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_m_gate.m_inf">
      <dcterms:description>steady state function of sodium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KM.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_h_gate.h">
      <dcterms:description>sodium current inactivation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_CaL.P_Ca">
      <dcterms:description>L-type calcium channel permeability to calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_CaL.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.k">
      <dcterms:description>proportionality constant</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.p_0">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in second messenger-bound P0 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.Ca_inf">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_27">
          <semsim:physicalPropertyOf rdf:resource="#entity_2" />
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>equilibrium concentration of intracellular calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_h.o_2">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in open O2 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.Ca_i">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_25">
          <semsim:physicalPropertyOf rdf:resource="#entity_2" />
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>concentration of intracelullar calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KM_p_gate.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_18">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_CaL.q">
      <dcterms:description>L-type calcium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KD_n_gate.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_h_gate.tau_h">
      <dcterms:description>time constant of sodium current inactivation gate</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_h_gate.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KD.g_KD">
      <dcterms:description>delayed rectifier potassium current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_Na.E_Na">
      <dcterms:description>sodium reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.h_inf">
      <dcterms:description>steady state function of sodium current inactivation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#CaL_q_gate.q">
      <dcterms:description>L-type calcium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_h.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.I_Na">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_2">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>sodium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KM_p_gate.p_inf">
      <dcterms:description>steady state function of slow potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KM.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_17">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#CaL_q_gate.q_inf">
      <dcterms:description>steady state function of L-type calcium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.I_CaL">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_28">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>L-type calcium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KD_n_gate.tau_n">
      <dcterms:description>activation time constant of delayed rectifier potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.T">
      <dcterms:description>absolute temperature</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.o_2">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in open O2 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.k_3p">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#environment.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_Na.E_Na">
      <dcterms:description>sodium reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.beta">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.I_CaL">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_5">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>L-type calcium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_h_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_h.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_31">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KD.n">
      <dcterms:description>delayed rectifier potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_h_gate.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_12">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.I_app">
      <dcterms:description>externally applied stimulus current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_leak.I_leak">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_7">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>leakage current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.I_KM">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_4">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>slow potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KM_p_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.k">
      <dcterms:description>proportionality constant</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_23">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KD_n_gate.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.k_1Ca">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.Ca_c">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_34">
          <semsim:physicalPropertyOf rdf:resource="#entity_2" />
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>calcium half-activation constant</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_m_gate.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KM_p_gate.tau_max">
      <dcterms:description>maximum activation time constant of slow potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KM_p_gate.tau_max">
      <dcterms:description>maximum activation time constant of slow potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.k_1Ca">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.I_KD">
      <dcterms:description>delayed rectifier potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.k_3p">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.k_1Ca">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.o_2">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in open O2 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_Na.h">
      <dcterms:description>sodium current inactivation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.F">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00089" />
    </rdf:Description>
    <rdf:Description rdf:about="#stimulus_protocol.period">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#Na_h_gate.h">
      <dcterms:description>sodium current inactivation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_m_gate.tau_m">
      <dcterms:description>time constant of sodium current activation gate</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.Ca_i">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_26">
          <semsim:physicalPropertyOf rdf:resource="#entity_2" />
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00340" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>concentration of intracelullar calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_Na.I_Na">
      <dcterms:description>sodium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.p_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in second messenger-bound P1 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KD_n_gate.n_inf">
      <dcterms:description>steady state function of delayed rectifier potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_h.I_h">
      <dcterms:description>hyperpolarisation-activated cationic current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#stimulus_protocol.tau">
      <dcterms:description>time constant of externally applied stimulus current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_m_gate.alpha">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.o_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in open O1 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#stimulus_protocol.tau">
      <dcterms:description>time constant of externally applied stimulus current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.T">
      <dcterms:description>absolute temperature</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KD.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_CaL.I_CaL">
      <dcterms:description>L-type calcium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#CaL_q_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KM_p_gate.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_m_gate.beta">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#Na_m_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_h_gate.beta_h">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.R">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00410" />
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.d">
      <dcterms:description>shell depth beneath membrane</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KD_n_gate.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#stimulus_protocol.i_stimStart">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.I_app">
      <dcterms:description>externally applied stimulus current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.G">
      <dcterms:description>function of voltage and ionic concentrations</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_CaL.G">
      <dcterms:description>function of voltage and ionic concentrations</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KD.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KM.g_KM">
      <dcterms:description>slow potassium current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_Na.g_Na">
      <dcterms:description>sodium current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_CaL.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#Na_h_gate.alpha_h">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.I_h">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_6">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>hyperpolarisation-activated cationic current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_m_gate.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.k_2">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#stimulus_protocol.I_app">
      <dcterms:description>externally applied stimulus current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KM.g_KM">
      <dcterms:description>slow potassium current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KD.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KD.g_KD">
      <dcterms:description>delayed rectifier potassium current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_CaL.P_Ca">
      <dcterms:description>L-type calcium channel permeability to calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_leak.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_8">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_CaL.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.F">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00089" />
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.Z">
      <dcterms:description>valence of calcium ion</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.c_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in closed state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_m_gate.tau_m">
      <dcterms:description>time constant of sodium current activation gate</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_m_gate.alpha">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.I_CaL">
      <dcterms:description>L-type calcium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#entity_3">
      <bqbiol:is rdf:resource="http://identifiers.org/fma/FMA:226054" />
    </rdf:Description>
    <rdf:Description rdf:about="#KD_n_gate.alpha_n">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_h.o_2">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in open O2 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.k_3p">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_Na.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.T">
      <dcterms:description>absolute temperature</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.k_4">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_KD.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KM.E_K">
      <dcterms:description>leakage reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_Na.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_h.E_h">
      <dcterms:description>hyperpolarisation-activated cationic current reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#environment.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#KD_n_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#KD_n_gate.n">
      <dcterms:description>delayed rectifier potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_Na.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_10">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_h.cac">
      <bqbiol:isPropertyOf rdf:resource="winograd_destexhe_sanchezvives_2008.cellml#entity_2" />
      <dcterms:description>calcium half-activation constant</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_h_gate.alpha_h">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_h.o_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in open O1 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_CaL.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_20">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_h_gate.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.Ca_c">
      <bqbiol:isPropertyOf rdf:resource="winograd_destexhe_sanchezvives_2008.cellml#entity_2" />
      <dcterms:description>calcium half-activation constant</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_CaL.G">
      <dcterms:description>function of voltage and ionic concentrations</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KM.p">
      <dcterms:description>slow potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.F">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00089" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_Na.I_Na">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_9">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>sodium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KM.I_KM">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_16">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>slow potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.k_4">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#Na_h_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_Na.m">
      <dcterms:description>sodium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#CaL_q_gate.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_21">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_h.o_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in open O1 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KD.E_K">
      <dcterms:description>leakage reversal potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KD.I_KD">
      <dcterms:description>delayed rectifier potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KD.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#stimulus_protocol.i_stimEnd">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_33">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#stimulus_protocol.i_stimStart">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.k_2">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#CaL_q_gate.tau_q">
      <dcterms:description>time constant of L-type calcium current activation gate</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#CaL_q_gate.alpha_q">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_CaL.I_CaL">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_19">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>L-type calcium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_Na.h">
      <dcterms:description>sodium current inactivation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.k_4">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.F">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00089" />
    </rdf:Description>
    <rdf:Description rdf:about="#CaL_q_gate.q">
      <dcterms:description>L-type calcium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#membrane.I_KM">
      <dcterms:description>slow potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#stimulus_protocol.i_stimEnd">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.R">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00410" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KM.p">
      <dcterms:description>slow potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.I_leak">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_1">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>leakage current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KD_n_gate.n_inf">
      <dcterms:description>steady state function of delayed rectifier potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.R">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00410" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KD_n_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#Na_h_gate.h_inf">
      <dcterms:description>steady state function of sodium current inactivation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.k_3p">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.p_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in second messenger-bound P1 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.alpha">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_Na.g_Na">
      <dcterms:description>sodium current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#G_nonlin.I_CaL">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_24">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>L-type calcium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.G">
      <dcterms:description>function of voltage and ionic concentrations</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_h_gate.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.R">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_00410" />
    </rdf:Description>
    <rdf:Description rdf:about="#rate_constants.k_2">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#Na_m_gate.m">
      <dcterms:description>sodium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#CaL_q_gate.tau_q">
      <dcterms:description>time constant of L-type calcium current activation gate</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_m_gate.beta">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.C_m">
      <dcterms:description>membrane capacitance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_h_gate.beta_h">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="#I_h.I_h">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_29">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>hyperpolarisation-activated cationic current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KM_p_gate.tau_p">
      <dcterms:description>activation time constant of slow potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_h_gate.tau_h">
      <dcterms:description>time constant of sodium current inactivation gate</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#Na_m_gate.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.beta">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#dCa_i_dt.tau_r">
      <dcterms:description>calcium removal time constant</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_h_gate.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#kinetic.c_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in closed state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#Na_m_gate.V_T">
      <dcterms:description>threshold voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#CaL_q_gate.V">
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.I_KD">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_3">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>delayed rectifier potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#winograd_destexhe_sanchez_vives_2008" bqmodel:isDescribedBy="18474856" />
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.Ca_i">
      <bqbiol:isPropertyOf rdf:resource="winograd_destexhe_sanchezvives_2008.cellml#entity_2" />
      <dcterms:description>concentration of intracelullar calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_h.g_hbar">
      <dcterms:description>hyperpolarisation-activated cationic current conductance</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#kinetic.k_1Ca">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#G_nonlin.Z">
      <dcterms:description>valence of calcium ion</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.p_1">
      <dcterms:description>fraction of hyperpolarisation-activated cationic channels in second messenger-bound P1 state</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KD.V">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_14">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00506" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>membrane potential</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#CaL_q_gate.alpha_q">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KD_n_gate.beta_n">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_leak.I_leak">
      <dcterms:description>leakage current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KM_p_gate.p_inf">
      <dcterms:description>steady state function of slow potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#KM_p_gate.p">
      <dcterms:description>slow potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#stimulus_protocol.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#I_KM.I_KM">
      <dcterms:description>slow potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#stimulus_protocol.time">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01023" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#rate_constants.beta">
      <bqbiol:is rdf:resource="https://identifiers.org/opb/OPB_01296" />
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#KD_n_gate.n">
      <dcterms:description>delayed rectifier potassium current activation variable</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#membrane.V_S">
      <dcterms:description>spike voltage</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.Ca_inf">
      <bqbiol:isPropertyOf rdf:resource="winograd_destexhe_sanchezvives_2008.cellml#entity_2" />
      <dcterms:description>equilibrium concentration of intracellular calcium</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="#I_KD.I_KD">
      <semsim:isComputationalComponentFor>
        <rdf:Description rdf:about="#property_13">
          <semsim:hasPhysicalDefinition rdf:resource="https://identifiers.org/opb/OPB_00318" />
        </rdf:Description>
      </semsim:isComputationalComponentFor>
      <dcterms:description>delayed rectifier potassium current</dcterms:description>
    </rdf:Description>
    <rdf:Description rdf:about="winograd_destexhe_sanchezvives_2008.cellml#dCa_i_dt.Ca_i">
      <bqbiol:isPropertyOf rdf:resource="winograd_destexhe_sanchezvives_2008.cellml#entity_2" />
      <dcterms:description>concentration of intracelullar calcium</dcterms:description>
    </rdf:Description>
  </rdf:RDF>
</model>