Location: Hodgkin, Huxley, 1952 @ 326854e3bbcb / hodgkin_huxley_1952.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2006-08-21 04:59:40+12:00
Desc:
committing version02 of hodgkin_huxley_1952
Permanent Source URI:
https://models.cellml.org/workspace/hodgkin_huxley_1952/rawfile/326854e3bbcb10fa69521f4d4df1c9cad764671a/hodgkin_huxley_1952.cellml

<?xml version='1.0' encoding='utf-8'?>
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="1952_hodgkin_huxley" name="hodgkin_huxley_1952_version02">
    
  

  <!--
    Generally, we want to move away from initial/default values being
    stored in the model directly. But until we are using CellML 1.1
    it is probably quite useful to include the values in the model.
  -->

  <!-- Global units -->
  <units name="mV">
    <unit units="volt" prefix="milli"/>
  </units>
  <units name="uApmmsq">
    <unit units="ampere" prefix="micro"/>
    <unit units="metre" prefix="milli" exponent="-2"/>
  </units>
  <units name="uApmmcu">
    <unit units="ampere" prefix="micro"/>
    <unit units="metre" prefix="milli" exponent="-3"/>
  </units>
  <units name="uFpmmsq">
    <unit units="farad" prefix="micro"/>
    <unit units="metre" prefix="milli" exponent="-2"/>
  </units>
  <units name="mSpmmsq">
    <unit units="siemens" prefix="milli"/>
    <unit units="metre" prefix="milli" exponent="-2"/>
  </units>
  <units name="ms">
    <unit units="second" prefix="milli"/>
  </units>

  <component cmeta:id="interface" name="interface">
  
    
  
    <!-- Variables we expect to be set/controlled externally -->
    <variable units="ms" private_interface="out" name="t"/>
    <variable units="uFpmmsq" private_interface="out" name="Cm" initial_value="0.01"/>
    <variable units="mV" private_interface="out" name="Er" initial_value="-65.0"/>
    <variable units="mSpmmsq" private_interface="out" name="gNa_max" initial_value="120.0e-2"/>
    <variable units="mSpmmsq" private_interface="out" name="gK_max" initial_value="36.0e-2"/>
    <variable units="mSpmmsq" private_interface="out" name="gleak_max" initial_value="0.3e-2"/>
    <variable units="uApmmcu" private_interface="out" name="Istim"/>
    
    <!-- Variables we want to make available externally -->
    <variable units="mV" public_interface="out" private_interface="in" name="Vm"/>
    <variable units="dimensionless" public_interface="out" private_interface="in" name="m"/>
    <variable units="dimensionless" public_interface="out" private_interface="in" name="h"/>
    <variable units="dimensionless" public_interface="out" private_interface="in" name="n"/>
    <variable units="uApmmsq" public_interface="out" private_interface="in" name="INa"/>
    <variable units="uApmmsq" public_interface="out" private_interface="in" name="IK"/>
    <variable units="uApmmsq" public_interface="out" private_interface="in" name="Ileak"/>
    
    

    <variable units="uApmmcu" public_interface="out" name="IStimC"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="IStim_for_cmiss_eq">
      <apply id="IStim_for_cmiss">
        <eq/>
        <ci>IStimC</ci>
        <ci>Istim</ci>
      </apply>
    </math>
    
  </component> <!--interface-->

  <component cmeta:id="membrane" name="membrane">
  
    

    <!-- Inputs -->
    <variable units="ms" public_interface="in" private_interface="out" name="t"/>
    <variable units="uFpmmsq" public_interface="in" private_interface="out" name="Cm"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Er"/>
    <variable units="mSpmmsq" public_interface="in" private_interface="out" name="gNa_max"/>
    <variable units="mSpmmsq" public_interface="in" private_interface="out" name="gK_max"/>
    <variable units="mSpmmsq" public_interface="in" private_interface="out" name="gleak_max"/>
    <variable units="uApmmcu" public_interface="in" private_interface="out" name="Istim"/>
    
    <!-- Outputs computed here -->
    <variable units="mV" public_interface="out" private_interface="out" name="Vm" initial_value="-65.0"/>
    
    <!-- Outputs coming from encapsulated components -->
    <variable units="dimensionless" public_interface="out" private_interface="in" name="m"/>
    <variable units="dimensionless" public_interface="out" private_interface="in" name="h"/>
    <variable units="dimensionless" public_interface="out" private_interface="in" name="n"/>
    <variable units="uApmmsq" public_interface="out" private_interface="in" name="INa"/>
    <variable units="uApmmsq" public_interface="out" private_interface="in" name="IK"/>
    <variable units="uApmmsq" public_interface="out" private_interface="in" name="Ileak"/> 
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Vm_diff_calculation_eq">
      <apply id="Vm_diff_calculation">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>Vm</ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci>Istim</ci>
            <apply>
              <plus/>
              <ci>INa</ci>
              <ci>IK</ci>
              <ci>Ileak</ci>
            </apply>
          </apply>
          <ci>Cm</ci>
        </apply>
      </apply>
    </math>
  </component>
  <!--membrane-->
  
  <component cmeta:id="sodium_channel" name="sodium_channel">

    

    <!-- Inputs -->
    <variable units="ms" public_interface="in" private_interface="out" name="t"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Er"/>
    <variable units="mSpmmsq" public_interface="in" private_interface="out" name="gNa_max"/>
    
    <!-- Outputs computed here -->
    <variable units="uApmmsq" public_interface="out" name="INa"/>
    
    <!-- Outputs from encapsulated components -->
    <variable units="dimensionless" public_interface="out" private_interface="in" name="m_gate"/>
    <variable units="dimensionless" public_interface="out" private_interface="in" name="h"/>
    
    <!-- Local variables -->
    <variable units="mSpmmsq" name="gNa"/>
    <variable units="mV" name="ENa"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="INa_calculation_eq">
      <apply id="INa_calculation">
        <eq/>
        <ci>INa</ci>
        <apply>
          <times/>
          <ci>gNa</ci>
          <apply>
            <minus/>
            <ci>Vm</ci>
            <ci>ENa</ci>
          </apply>
        </apply>
      </apply>
    </math>
      
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="gNa_calculation_eq">
      <apply id="gNa_calculation">
        <eq/>
        <ci>gNa</ci>
        <apply>
          <times/>
          <ci>gNa_max</ci>
          <apply>
            <power/>
            <ci>m_gate</ci>
            <cn cellml:units="dimensionless">3.0</cn>
          </apply>
          <ci>h</ci>
        </apply>
      </apply>
      
    </math>
      
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="ENa_calculation_eq">
      <apply id="ENa_calculation">
        <eq/>
        <ci>ENa</ci>
        <apply>
          <plus/>
          <ci>Er</ci>
          <cn cellml:units="mV">115.0</cn>
        </apply>
      </apply>
    </math>
  </component>
  <!--sodium_channel-->
  
  <component cmeta:id="sodium_channel_m_gate" name="sodium_channel_m_gate">

    

    <!-- Local units -->
    <units name="pms">
      <unit units="ms" exponent="-1"/>
    </units>
    
    <!-- Inputs -->
    <variable units="ms" public_interface="in" private_interface="out" name="t"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Er"/>
    
    <!-- Outputs computed here -->
    <variable units="dimensionless" public_interface="out" name="m" initial_value="0.0"/>
    
    <!-- Local variables -->
    <variable units="pms" name="alpha"/>
    <variable units="pms" name="beta"/>
    <variable units="mV" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="alpha_m_calculation_eq">
      <apply id="alpha_m_calculation">
        <eq/>
        <ci>alpha</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="pms">0.1</cn>
            <apply>
              <minus/>
              <cn cellml:units="mV">25.0</cn>
              <ci>V</ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.1 </cn>
                <apply>
                  <minus/>
                  <cn cellml:units="mV">25.0</cn>
                  <ci>V</ci>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">1.0</cn>
          </apply>
        </apply>
      </apply>
    </math>

    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="beta_m_calculation_eq">
      <apply id="beta_m_calculation">
        <eq/>
        <ci>beta</ci>
        <apply>
          <times/>
          <cn cellml:units="pms">4.0</cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>V</ci>
              </apply>
              <cn cellml:units="mV">18.0</cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>

    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="dm_dt_eq">
      <apply id="dm_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>m</ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci>alpha</ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1.0</cn>
              <ci>m</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>beta</ci>
            <ci>m</ci>
          </apply>
        </apply>
      </apply>
    </math>

    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="m_V_eq">
      <apply id="m_V">
        <eq/>
        <ci>V</ci>
        <apply>
          <minus/>
          <ci>Vm</ci>
          <ci>Er</ci>
        </apply>
      </apply>
    </math>
  </component>
  <!--sodium_channel_m_gate-->
  
  <component cmeta:id="sodium_channel_h_gate" name="sodium_channel_h_gate">
  
    

    <!-- Local units -->
    <units name="pms">
      <unit units="ms" exponent="-1"/>
    </units>
    
    <!-- Inputs -->
    <variable units="ms" public_interface="in" private_interface="out" name="t"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Er"/>
    
    <!-- Outputs computed here -->
    <variable units="dimensionless" public_interface="out" name="h" initial_value="0.6"/>
    
    <!-- Local variables -->
    <variable units="pms" name="alpha"/>
    <variable units="pms" name="beta"/>
    <variable units="mV" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="beta_h_calculation_eq">
      <apply id="beta_h_calculation">
        <eq/>
        <ci>beta</ci>
        <apply>
          <divide/>
          <cn cellml:units="pms">1.0</cn>
          <apply>
            <plus/>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless">0.1</cn>
                <apply>
                  <minus/>
                  <cn cellml:units="mV">30.0</cn>
                  <ci>V</ci>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">1.0</cn>
          </apply>
        </apply>
      </apply>
    </math>

    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="alpha_h_calculation_eq">
      <apply id="alpha_h_calculation">
        <eq/>
        <ci>alpha</ci>
        <apply>
          <times/>
          <cn cellml:units="pms">0.07</cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>V</ci>
              </apply>
              <cn cellml:units="mV">20.0</cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>

    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="dh_dt_eq">
      <apply id="dh_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>h</ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci>alpha</ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1.0</cn>
              <ci>h</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>beta</ci>
            <ci>h</ci>
          </apply>
        </apply>
      </apply>
    </math>
      
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="h_V_eq">
      <apply id="h_V">
        <eq/>
        <ci>V</ci>
        <apply>
          <minus/>
          <ci>Vm</ci>
          <ci>Er</ci>
        </apply>
      </apply>
    </math>
  </component>
  <!--sodium_channel_h_gate-->
    
  <component cmeta:id="potassium_channel" name="potassium_channel">
  
    

    <!-- Inputs -->
    <variable units="ms" public_interface="in" private_interface="out" name="t"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Er"/>
    <variable units="mSpmmsq" public_interface="in" private_interface="out" name="gK_max"/>
    
    <!-- Outputs computed here -->
    <variable units="uApmmsq" public_interface="out" name="IK"/>
    
    <!-- Outputs from encapsulated components -->
    <variable units="dimensionless" public_interface="out" private_interface="in" name="n"/>
    
    <!-- Local variables -->
    <variable units="mSpmmsq" name="gK"/>
    <variable units="mV" name="EK"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="IK_calculation_eq">
      <apply id="IK_calculation">
        <eq/>
        <ci>IK</ci>
        <apply>
          <times/>
          <ci>gK</ci>
          <apply>
            <minus/>
            <ci>Vm</ci>
            <ci>EK</ci>
          </apply>
        </apply>
      </apply>
    </math>
      
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="gK_calculation_eq">
      <apply id="gK_calculation">
        <eq/>
        <ci>gK</ci>
        <apply>
          <times/>
          <ci>gK_max</ci>
          <apply>
            <power/>
            <ci>n</ci>
            <cn cellml:units="dimensionless">4.0</cn>
          </apply>
        </apply>
      </apply>
    </math>
      
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="EK_calculation_eq">      
      <apply id="EK_calculation">
        <eq/>
        <ci>EK</ci>
        <apply>
          <minus/>
          <ci>Er</ci>
          <cn cellml:units="mV">12.0</cn>
        </apply>
      </apply>
    </math>
  </component>
  <!--potassium_channel-->
  
  <component cmeta:id="potassium_channel_n_gate" name="potassium_channel_n_gate">
  
    

    <!-- Local units -->
    <units name="pms">
      <unit units="ms" exponent="-1"/>
    </units>
    
    <!-- Inputs -->
    <variable units="ms" public_interface="in" private_interface="out" name="t"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Er"/>
    
    <!-- Outputs computed here -->
    <variable units="dimensionless" public_interface="out" name="n" initial_value="0.3"/>
    
    <!-- Local variables -->
    <variable units="pms" name="alpha"/>
    <variable units="pms" name="beta"/>
    <variable units="mV" name="V"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="alpha_n_calculation_eq">
      <apply id="alpha_n_calculation">
        <eq/>
        <ci>alpha</ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="pms">0.01</cn>
            <apply>
              <minus/>
              <cn cellml:units="mV">10.0</cn>
              <ci>V</ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.1 </cn>
                <apply>
                  <minus/>
                  <cn cellml:units="mV">10.0</cn>
                  <ci>V</ci>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless">1.0</cn>
          </apply>
        </apply>
      </apply>
    </math>

    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="beta_n_calculation_eq">
      <apply id="beta_n_calculation">
        <eq/>
        <ci>beta</ci>
        <apply>
          <times/>
          <cn cellml:units="pms">0.125</cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <ci>V</ci>
              </apply>
              <cn cellml:units="mV">80.0</cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>

    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="dn_dt_eq">
      <apply id="dn_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>t</ci>
          </bvar>
          <ci>n</ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci>alpha</ci>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless">1.0</cn>
              <ci>n</ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci>beta</ci>
            <ci>n</ci>
          </apply>
        </apply>
      </apply>
    </math>
      
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="n_V_eq">
      <apply id="n_V">
        <eq/>
        <ci>V</ci>
        <apply>
          <minus/>
          <ci>Vm</ci>
          <ci>Er</ci>
        </apply>
      </apply>
    </math>
  </component>
  <!--potassium_channel_n_gate-->
  
  <component cmeta:id="leak_channel" name="leak_channel">

    
    
    <!-- Inputs -->
    <variable units="mV" public_interface="in" private_interface="out" name="Vm"/>
    <variable units="mV" public_interface="in" private_interface="out" name="Er"/>
    <variable units="mSpmmsq" public_interface="in" private_interface="out" name="gleak_max"/>
    
    <!-- Outputs computed here -->
    <variable units="uApmmsq" public_interface="out" name="Ileak"/>
    
    <!-- Local variables -->
    <variable units="mV" name="Eleak"/>
    
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Ileak_calculation_eq">
      <apply id="Ileak_calculation">
        <eq/>
        <ci>Ileak</ci>
        <apply>
          <times/>
          <ci>gleak_max</ci>
          <apply>
            <minus/>
            <ci>Vm</ci>
            <ci>Eleak</ci>
          </apply>
        </apply>
      </apply>
    </math>
      
    
    <math xmlns="http://www.w3.org/1998/Math/MathML" cmeta:id="Eleak_calculation_eq">
      <apply id="Eleak_calculation">
        <eq/>
        <ci>Eleak</ci>
        <apply>
          <plus/>
          <ci>Er</ci>
          <cn cellml:units="mV">10.613</cn>
        </apply>
      </apply>
    </math>
  </component>
  <!--leak_channel-->
  
  <connection>
    <map_components component_2="membrane" component_1="interface"/>
    <map_variables variable_2="t" variable_1="t"/>
    <map_variables variable_2="Cm" variable_1="Cm"/>
    <map_variables variable_2="Er" variable_1="Er"/>
    <map_variables variable_2="gNa_max" variable_1="gNa_max"/>
    <map_variables variable_2="gK_max" variable_1="gK_max"/>
    <map_variables variable_2="gleak_max" variable_1="gleak_max"/>
    <map_variables variable_2="Istim" variable_1="Istim"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="m" variable_1="m"/>
    <map_variables variable_2="h" variable_1="h"/>
    <map_variables variable_2="n" variable_1="n"/>
    <map_variables variable_2="INa" variable_1="INa"/>
    <map_variables variable_2="IK" variable_1="IK"/>
    <map_variables variable_2="Ileak" variable_1="Ileak"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_channel" component_1="membrane"/>
    <map_variables variable_2="t" variable_1="t"/>
    <map_variables variable_2="Er" variable_1="Er"/>
    <map_variables variable_2="gNa_max" variable_1="gNa_max"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="m_gate" variable_1="m"/>
    <map_variables variable_2="h" variable_1="h"/>
    <map_variables variable_2="INa" variable_1="INa"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_channel_m_gate" component_1="sodium_channel"/>
    <map_variables variable_2="t" variable_1="t"/>
    <map_variables variable_2="Er" variable_1="Er"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="m" variable_1="m_gate"/>
  </connection>
  
  <connection>
    <map_components component_2="sodium_channel_h_gate" component_1="sodium_channel"/>
    <map_variables variable_2="t" variable_1="t"/>
    <map_variables variable_2="Er" variable_1="Er"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="h" variable_1="h"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_channel" component_1="membrane"/>
    <map_variables variable_2="t" variable_1="t"/>
    <map_variables variable_2="Er" variable_1="Er"/>
    <map_variables variable_2="gK_max" variable_1="gK_max"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="n" variable_1="n"/>
    <map_variables variable_2="IK" variable_1="IK"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_channel_n_gate" component_1="potassium_channel"/>
    <map_variables variable_2="t" variable_1="t"/>
    <map_variables variable_2="Er" variable_1="Er"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="n" variable_1="n"/>
  </connection>
  
  <connection>
    <map_components component_2="leak_channel" component_1="membrane"/>
    <map_variables variable_2="Er" variable_1="Er"/>
    <map_variables variable_2="gleak_max" variable_1="gleak_max"/>
    <map_variables variable_2="Vm" variable_1="Vm"/>
    <map_variables variable_2="Ileak" variable_1="Ileak"/>
  </connection>
  
  <group>
    <relationship_ref relationship="encapsulation"/>
    <component_ref component="interface">
      <component_ref component="membrane">
        <component_ref component="sodium_channel">
          <component_ref component="sodium_channel_m_gate"/>
          <component_ref component="sodium_channel_h_gate"/>
        </component_ref>
        <component_ref component="potassium_channel">
          <component_ref component="potassium_channel_n_gate"/>
        </component_ref>
        <component_ref component="leak_channel"/>
      </component_ref>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="sodium_channel">
        <component_ref component="sodium_channel_m_gate"/>
        <component_ref component="sodium_channel_h_gate"/>
      </component_ref>
      <component_ref component="potassium_channel">
        <component_ref component="potassium_channel_n_gate"/>
      </component_ref>
      <component_ref component="leak_channel"/>
    </component_ref>    
  </group>







<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#citationAuthorsSeq">
    <rdf:li rdf:resource="rdf:#author1Vcard"/>
    <rdf:li rdf:resource="rdf:#author2Vcard"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#97221e0f-38d6-44d3-9139-0915b24e0805">
    <rdf:value>
            Calculation of the closing rate of the activation gate for the
            potassium channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2069ade8-05f4-4ac3-8b61-3825fb83471e">
    <rdf:value>
            Calculation of the potassium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#360e8bd9-354f-4e18-9602-2d7c172fb788">
    <rdf:value>
            Calculation of the sodium channel's current equilibrium potential.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#beta_h_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#c617a128-348e-4d67-9e4a-7c00d98a59a6"/>
  </rdf:Description>
  <rdf:Description rdf:about="#interface">
    <cmeta:comment rdf:resource="rdf:#1aa41f67-af3e-4665-abc2-3016c525c2d7"/>
  </rdf:Description>
  <rdf:Description rdf:about="#ENa_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#360e8bd9-354f-4e18-9602-2d7c172fb788"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6a310dc0-299d-443c-9754-a3d2b3751ed2">
    <dc:creator rdf:resource="rdf:#citationAuthorsSeq"/>
    <dc:title>A quantitative description of membrane current and its application to conductance and excitation in nerve</dc:title>
    <bqs:volume>117</bqs:volume>
    <bqs:first_page>500</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#88a878fe-067e-409a-82f0-336b4751d2ce"/>
    <dcterms:issued rdf:resource="rdf:#c64902f3-3c52-4fe3-b963-8de7042482c0"/>
    <bqs:last_page>544</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author2Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author2VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author1Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author1VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="#gNa_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#206c34f0-3609-4107-83ce-553a1d4d0033"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7ed94e38-add3-411d-b22a-609d22061947">
    <rdf:value>
            This is the voltage dependent activation gate for the sodium 
            channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#Ileak_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#a245fb62-f072-41a7-acbf-301a38260e93"/>
  </rdf:Description>
  <rdf:Description rdf:about="#dn_dt_eq">
    <cmeta:comment rdf:resource="rdf:#a3256a0c-40a9-4251-af11-f5b0515b6891"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dcterms:created rdf:resource="rdf:#c15e62dc-7222-4dbc-b076-95eacf59b717"/>
    <dc:creator rdf:resource="rdf:#ec976bfd-1119-441c-857d-c1d590d82ddd"/>
  </rdf:Description>
  <rdf:Description rdf:about="#alpha_n_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#a984b10a-2780-4e9c-acff-edb8e3bb30fb"/>
  </rdf:Description>
  <rdf:Description rdf:about="#INa_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#c0db5eec-4571-47ab-965e-93839a909728"/>
  </rdf:Description>
  <rdf:Description rdf:about="#Vm_diff_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#5f5b7f05-490e-43ae-9d0f-e8da76c23425"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#725c8afe-f091-4d91-b04e-5d69d803aed4">
    <bqs:Pubmed_id>11773311</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#6a310dc0-299d-443c-9754-a3d2b3751ed2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author2VcardN">
    <vCard:Given>A</vCard:Given>
    <vCard:Family>Huxley</vCard:Family>
    <vCard:Other>F</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#206c34f0-3609-4107-83ce-553a1d4d0033">
    <rdf:value>
            Calculation of the sodium channel's current conductance based on 
            the activation and inactivation gates.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5f5b7f05-490e-43ae-9d0f-e8da76c23425">
    <rdf:value>
            This is the main equation for the model, describing the dependency
            of the time course of tranmembrane potential on the three ionic
            currents plus any applied stimulus current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ec976bfd-1119-441c-857d-c1d590d82ddd">
    <vCard:ORG rdf:resource="rdf:#52b5d41f-6334-4688-8e30-178823e5229d"/>
    <vCard:EMAIL rdf:resource="rdf:#8486516d-83f5-4429-8ba8-b69e46596be7"/>
    <vCard:N rdf:resource="rdf:#7c5bac4f-43a7-47d9-b701-98697532ca8d"/>
  </rdf:Description>
  <rdf:Description rdf:about="#beta_n_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#97221e0f-38d6-44d3-9139-0915b24e0805"/>
  </rdf:Description>
  <rdf:Description rdf:about="#membrane">
    <cmeta:comment rdf:resource="rdf:#2e6fbe8c-b57d-48dd-8b83-64c090bce016"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d042d5c8-76c9-48b8-b99b-7a6fc72a4f1a">
    <rdf:value>This is a CellML version of the model published in 1952 by Hodgkin and Huxley on the electrical activation of squid axons.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5bd52e00-dbad-4ee2-ab2e-68584cf803cc">
    <rdf:value>
            The sodium channel is the first of three membrane bound ionic 
            channel proteins in this model. It is primarily responsible
            for the upstroke of the action potential.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8486516d-83f5-4429-8ba8-b69e46596be7">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>d.nickerson@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#16953998-ed7f-4d07-8fd0-396be1864840">
    <rdf:value>
            A small leakage current, made up by chloride and other ions, with a
            constant conductance.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#sodium_channel_m_gate">
    <cmeta:comment rdf:resource="rdf:#7ed94e38-add3-411d-b22a-609d22061947"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c64902f3-3c52-4fe3-b963-8de7042482c0">
    <dcterms:W3CDTF>1952-01-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3f059b29-7155-4c9d-9b6d-bc70d21958ae">
    <rdf:value>
            Calculation of the opening rate of the activation gate of the sodium
            channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#98fbc281-a0ae-4df1-a10d-689e023d64d9">
    <rdf:value>
            This is a dummy equation that we simply use to make grabbing the
            value in CMISS much easier.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d2e52a53-7b61-4fd5-8eaa-f13dd9a17db8">
    <rdf:value>
            Calculation of the equilibrium potential of the leak channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#dh_dt_eq">
    <cmeta:comment rdf:resource="rdf:#cc96406c-d1f8-4d7f-bf97-264cd3d562bd"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a3256a0c-40a9-4251-af11-f5b0515b6891">
    <rdf:value>
            The rate of change in the activation gate of the potassium channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a0ac5b78-5653-430b-93ac-53a8cfc5f4e8">
    <rdf:value>
            Calculation of the potassium channel conductance.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a245fb62-f072-41a7-acbf-301a38260e93">
    <rdf:value>
            Calculation of the leakage current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#alpha_m_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#3f059b29-7155-4c9d-9b6d-bc70d21958ae"/>
  </rdf:Description>
  <rdf:Description rdf:about="#potassium_channel_n_gate">
    <cmeta:comment rdf:resource="rdf:#99d95ddb-2c92-484f-8314-e3ced2a42834"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c4954a42-0d5d-41c7-abc0-b24e9be53794">
    <rdf:value>
            Calculation of the potassium channel equilibrium potential.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c617a128-348e-4d67-9e4a-7c00d98a59a6">
    <rdf:value>
            The closing rate of the inactivation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#sodium_channel_h_gate">
    <cmeta:comment rdf:resource="rdf:#42aefda5-8a2a-4c16-84dd-174ab762429f"/>
  </rdf:Description>
  <rdf:Description rdf:about="#gK_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#a0ac5b78-5653-430b-93ac-53a8cfc5f4e8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c0db5eec-4571-47ab-965e-93839a909728">
    <rdf:value>
            Calculation of the sodium current.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1aa41f67-af3e-4665-abc2-3016c525c2d7">
    <rdf:value>
            We'll use this component as the `interface' to the model, all 
            other components are hidden via encapsulation in this component.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#sodium_channel">
    <cmeta:comment rdf:resource="rdf:#5bd52e00-dbad-4ee2-ab2e-68584cf803cc"/>
  </rdf:Description>
  <rdf:Description rdf:about="#beta_m_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#803d5066-3967-40dd-844a-93a342265dfb"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author1VcardN">
    <vCard:Given>A</vCard:Given>
    <vCard:Family>Hodgkin</vCard:Family>
    <vCard:Other>L</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#42aefda5-8a2a-4c16-84dd-174ab762429f">
    <rdf:value>
            This is the voltage dependent inactivation gate for the sodium
            channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#IK_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#2069ade8-05f4-4ac3-8b61-3825fb83471e"/>
  </rdf:Description>
  <rdf:Description rdf:about="#1952_hodgkin_huxley">
    <bqs:reference rdf:resource="rdf:#725c8afe-f091-4d91-b04e-5d69d803aed4"/>
    <cmeta:comment rdf:resource="rdf:#d042d5c8-76c9-48b8-b99b-7a6fc72a4f1a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a984b10a-2780-4e9c-acff-edb8e3bb30fb">
    <rdf:value>
            Calculation of the opening rate of the activation gate for the
            potassium channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#IStim_for_cmiss_eq">
    <cmeta:comment rdf:resource="rdf:#98fbc281-a0ae-4df1-a10d-689e023d64d9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#52b5d41f-6334-4688-8e30-178823e5229d">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="#EK_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#c4954a42-0d5d-41c7-abc0-b24e9be53794"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#88a878fe-067e-409a-82f0-336b4751d2ce">
    <dc:title>Journal of Physiology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cc96406c-d1f8-4d7f-bf97-264cd3d562bd">
    <rdf:value>
            The rate of change of the inactivation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#803d5066-3967-40dd-844a-93a342265dfb">
    <rdf:value>
            Calculation of the closing rate of the activation gate of the sodium
            channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#alpha_h_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#4e5a9a69-6ec0-490a-ac7a-386d14747f25"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7c5bac4f-43a7-47d9-b701-98697532ca8d">
    <vCard:Given>David</vCard:Given>
    <vCard:Family>Nickerson</vCard:Family>
    <vCard:Other>P</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7493db27-66b9-496d-98e7-b7bc8fca632f">
    <rdf:value>
            The potassium channel is the second of three membrane bound ionic
            channels defined by this model. It is primarily responsible for the
            repolarisation phase of the action potential.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#potassium_channel">
    <cmeta:comment rdf:resource="rdf:#7493db27-66b9-496d-98e7-b7bc8fca632f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2e6fbe8c-b57d-48dd-8b83-64c090bce016">
    <rdf:value>
            The membrane component is the "real" parent component which 
            physically contains its child components. This component defines
            the transmembrane potential.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#Eleak_calculation_eq">
    <cmeta:comment rdf:resource="rdf:#d2e52a53-7b61-4fd5-8eaa-f13dd9a17db8"/>
  </rdf:Description>
  <rdf:Description rdf:about="#dm_dt_eq">
    <cmeta:comment rdf:resource="rdf:#498240df-fe17-4929-93ea-c39b182e167a"/>
  </rdf:Description>
  <rdf:Description rdf:about="#leak_channel">
    <cmeta:comment rdf:resource="rdf:#16953998-ed7f-4d07-8fd0-396be1864840"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#498240df-fe17-4929-93ea-c39b182e167a">
    <rdf:value>
            The rate of change of the activation gate of the sodium channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#99d95ddb-2c92-484f-8314-e3ced2a42834">
    <rdf:value>
            The voltage dependent activation gate for the potassium channel.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4e5a9a69-6ec0-490a-ac7a-386d14747f25">
    <rdf:value>
            The opening rate of the inactivation gate.
          </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c15e62dc-7222-4dbc-b076-95eacf59b717">
    <dcterms:W3CDTF>2003-01-31</dcterms:W3CDTF>
  </rdf:Description>
</rdf:RDF>
</model>