<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE NEURON SYSTEM "http://www.neuron.yale.edu/neuron/DTDs/neuron.dtd">
<NEURON version="5.7">
<!-- ./nak.xml -->
<!-- XML file generated by NEURON 5.6 ModelViewer -->
<!-- Authors: Michael Hines and Sushil Kambampati -->
<!-- Yale University -->
<!-- Date: Sat Apr 30 09:48:17 EDT 2005
 -->
<model>
<morphology xmlns="http://morphml.org/namespace">
  <setOfPoints>
  <point id="1" x="0" y="0" z="0" diam="10"/>
  <point id="2" x="1.59155" y="0" z="0" diam="10"/>
  <point id="3" x="3.1831" y="0" z="0" diam="10"/>
  </setOfPoints>
<cell>
  <setOfSegments name="soma">
    <segment id="1" proximal="1" distal="2"/>
    <segment id="2" proximal="2" distal="3" parent="1"/>
  </setOfSegments>
</cell>
</morphology>
<biophysics>
  <ion name="na" charge="1"/>
  <ion name="k" charge="1"/>
  <channel_type name="nahh" density="yes">
    <current_voltage_relation>
      <ohmic ion="na">
        <conductance default_gmax="0.12">
        <gate power="3">
          <state name="m" fraction="1"/>
        </gate>
        <gate power="1">
          <state name="h" fraction="1"/>
        </gate>
        </conductance>
      </ohmic>
    </current_voltage_relation>
    <hh_gate state="m">
      <transition>
        <voltage_gate infinity_tau_format="no">
          <alpha type="linoid" expr="A*x/(1 - exp(-x)) where x = k*v-d)">
            <parameter name="A" value="1"/>
            <parameter name="k" value="0.1"/>
            <parameter name="d" value="-40"/>
          </alpha>
          <beta type="exponential" expr="A*exp(k*(v-d))">
            <parameter name="A" value="4"/>
            <parameter name="k" value="-0.0555556"/>
            <parameter name="d" value="-65"/>
          </beta>
        </voltage_gate>
      </transition>
    </hh_gate>
    <hh_gate state="h">
      <transition>
        <voltage_gate infinity_tau_format="no">
          <alpha type="exponential" expr="A*exp(k*(v-d))">
            <parameter name="A" value="0.07"/>
            <parameter name="k" value="-0.05"/>
            <parameter name="d" value="-65"/>
          </alpha>
          <beta type="sigmoid" expr="A/(1 + exp(-k*(d-v))">
            <parameter name="A" value="1"/>
            <parameter name="k" value="-0.1"/>
            <parameter name="d" value="-35"/>
          </beta>
        </voltage_gate>
      </transition>
    </hh_gate>
  </channel_type>
  <channel_type name="khh" density="yes">
    <current_voltage_relation>
      <ohmic ion="k">
        <conductance default_gmax="0.036">
        <gate power="1">
          <state name="n4" fraction="1"/>
        </gate>
        </conductance>
      </ohmic>
    </current_voltage_relation>
    <ks_gate>
      <state name="n0"/>
      <state name="n1"/>
      <state name="n2"/>
      <state name="n3"/>
      <state name="n4"/>
      <transition src="n0" target="n1">
        <voltage_gate infinity_tau_format="no">
          <alpha type="linoid" expr="A*x/(1 - exp(-x)) where x = k*v-d)">
            <parameter name="A" value="0.4"/>
            <parameter name="k" value="0.1"/>
            <parameter name="d" value="-55"/>
          </alpha>
          <beta type="exponential" expr="A*exp(k*(v-d))">
            <parameter name="A" value="0.125"/>
            <parameter name="k" value="-0.0125"/>
            <parameter name="d" value="-65"/>
          </beta>
        </voltage_gate>
      </transition>
      <transition src="n1" target="n2">
        <voltage_gate infinity_tau_format="no">
          <alpha type="linoid" expr="A*x/(1 - exp(-x)) where x = k*v-d)">
            <parameter name="A" value="0.3"/>
            <parameter name="k" value="0.1"/>
            <parameter name="d" value="-55"/>
          </alpha>
          <beta type="exponential" expr="A*exp(k*(v-d))">
            <parameter name="A" value="0.25"/>
            <parameter name="k" value="-0.0125"/>
            <parameter name="d" value="-65"/>
          </beta>
        </voltage_gate>
      </transition>
      <transition src="n2" target="n3">
        <voltage_gate infinity_tau_format="no">
          <alpha type="linoid" expr="A*x/(1 - exp(-x)) where x = k*v-d)">
            <parameter name="A" value="0.2"/>
            <parameter name="k" value="0.1"/>
            <parameter name="d" value="-55"/>
          </alpha>
          <beta type="exponential" expr="A*exp(k*(v-d))">
            <parameter name="A" value="0.375"/>
            <parameter name="k" value="-0.0125"/>
            <parameter name="d" value="-65"/>
          </beta>
        </voltage_gate>
      </transition>
      <transition src="n3" target="n4">
        <voltage_gate infinity_tau_format="no">
          <alpha type="linoid" expr="A*x/(1 - exp(-x)) where x = k*v-d)">
            <parameter name="A" value="0.1"/>
            <parameter name="k" value="0.1"/>
            <parameter name="d" value="-55"/>
          </alpha>
          <beta type="exponential" expr="A*exp(k*(v-d))">
            <parameter name="A" value="0.5"/>
            <parameter name="k" value="-0.0125"/>
            <parameter name="d" value="-65"/>
          </beta>
        </voltage_gate>
      </transition>
    </ks_gate>
  </channel_type>
  <channel_type name="leak" density="yes">
    <current_voltage_relation>
      <ohmic NonSpecific="yes" default_erev="-54.3">
        <conductance default_gmax="0.0003">
        </conductance>
      </ohmic>
    </current_voltage_relation>
  </channel_type>
<cell name="soma">
  <mechanism name="Ra">
    <parameter name="" value="35.4">
      <section name="soma"/>
    </parameter>
  </mechanism>
  <mechanism name="capacitance">
    <parameter name="cm" value="1">
      <section name="soma"/>
    </parameter>
  </mechanism>
  <mechanism name="na_ion">
    <parameter name="ena" value="50">
      <section name="soma"/>
    </parameter>
  </mechanism>
  <mechanism name="k_ion">
    <parameter name="ek" value="-77">
      <section name="soma"/>
    </parameter>
  </mechanism>
  <mechanism name="nahh">
    <parameter name="gmax_nahh" value="0.12">
      <section name="soma"/>
    </parameter>
  </mechanism>
  <mechanism name="khh">
    <parameter name="gmax_khh" value="0.036">
      <section name="soma"/>
    </parameter>
  </mechanism>
  <mechanism name="leak">
    <parameter name="gmax_leak" value="0.0003">
      <section name="soma"/>
    </parameter>
    <parameter name="e_leak" value="-54.3">
      <section name="soma"/>
    </parameter>
  </mechanism>
</cell>
</biophysics>
</model>
</NEURON>
