<?xml version="1.0" encoding="UTF-8"?>
<channelml xmlns="http://morphml.org/channelml/schema" 
           xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" 
           xmlns:meta="http://morphml.org/metadata/schema" 
           xsi:schemaLocation="http://morphml.org/channelml/schema  http://www.neuroml.org/NeuroMLValidator/NeuroMLFiles/Schemata/v1.8.1/Level2/ChannelML_v1.8.1.xsd"
           units="Physiological Units">

    <meta:notes>ChannelML file based on Traub et al. 2003</meta:notes>

    <channel_type name="kc_fast" density="yes">
        
        <status value="stable">
            <meta:comment>Verified equivalence of NEURON and GENESIS mapping to orig NEURON mod impl at 0.02ms dt with current pulse</meta:comment>
            <meta:issue>In original mod, m is initialised to 0, as opposed to minf at t=0. ChannelML impl corrects this</meta:issue>
            <meta:contributor>
                <meta:name>Padraig Gleeson</meta:name>
            </meta:contributor>
            <meta:contributor>
                <meta:name>Yoana Dimitrova</meta:name>
            </meta:contributor>
        </status>
        
        <meta:notes>Fast voltage and [Ca2+] dependent K conductance (BK channel). Based on NEURON port of FRB L2/3 model from Traub et al 2003. Same channel used in Traub et al 2005</meta:notes>
        
        <meta:authorList>
            <meta:modelAuthor>
                <meta:name>Roger D Traub</meta:name>
                <meta:comment>Conversion of FORTRAN model to NEURON. See Traub et al 2003 for detailed origin of channels</meta:comment>
            </meta:modelAuthor>
            <meta:modelAuthor>
                <meta:name>Maciej Lazarewicz</meta:name>
                <meta:comment>Conversion of FORTRAN model to NEURON. See Traub et al 2003 for detailed origin of channels</meta:comment>
            </meta:modelAuthor>
            <meta:modelTranslator>
                <meta:name>Padraig Gleeson</meta:name>
                <meta:institution>UCL</meta:institution>
                <meta:email>p.gleeson - at - ucl.ac.uk</meta:email>
            </meta:modelTranslator>
            <meta:modelTranslator>
                <meta:name>Yoana Dimitrova</meta:name>
                <meta:institution>UCL</meta:institution>
            </meta:modelTranslator>
        </meta:authorList>

        <meta:publication>
            <meta:fullTitle>Roger D. Traub, Eberhard H. Buhl, Tengis Gloveli, and Miles A. Whittington                
Fast Rhythmic Bursting Can Be Induced in Layer 2/3 Cortical Neurons by Enhancing Persistent Na+ Conductance or by Blocking BK Channels
J Neurophysiol 89: 909-921, 2003</meta:fullTitle>
            <meta:pubmedRef>http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=pubmed&amp;dopt=Abstract&amp;list_uids=12574468</meta:pubmedRef>
        </meta:publication>

        <meta:publication>
            <meta:fullTitle>Roger D. Traub, Diego Contreras, Mark O. Cunningham, Hilary Murray, Fiona E. N. LeBeau, Anita Roopun, Andrea Bibbig, W. Bryan Wilent, Michael J. Higley, and Miles A. Whittington
Single-column thalamocortical network model exhibiting gamma oscillations, sleep spindles, and epileptogenic bursts.
J. Neurophysiol. 93, 2194-2232, 2005</meta:fullTitle>
            <meta:pubmedRef>http://www.ncbi.nlm.nih.gov/pubmed/15525801?dopt=Abstract</meta:pubmedRef>
        </meta:publication>
        
        <meta:neuronDBref>
            <meta:modelName>K channels</meta:modelName>
            <meta:uri>http://senselab.med.yale.edu/senselab/NeuronDB/channelGene2.htm#table3</meta:uri>
        </meta:neuronDBref>
        
        <meta:modelDBref>
            <meta:modelName>Mechanisms of fast rhythmic bursting in a layer 2/3 cortical neuron (Traub et al 2003)</meta:modelName>
            <meta:uri>http://senselab.med.yale.edu/ModelDB/ShowModel.asp?model=20756</meta:uri>
        </meta:modelDBref>

        <meta:modelDBref>
            <meta:modelName>A single column thalamocortical network model (Traub et al 2005)</meta:modelName>
            <meta:uri>http://senselab.med.yale.edu/ModelDb/ShowModel.asp?model=45539</meta:uri>
        </meta:modelDBref>

        <parameters>
            <parameter name="speedup_factor" value ="2"/>  
        </parameters>
    
        <current_voltage_relation cond_law="ohmic" ion="k" default_gmax="12" default_erev="-95">
            
             <!-- New element in v1.7.1, used for time independent concentration factor -->
            <conc_factor  ion="ca" charge="2" variable_name="ca_conc" expr="ca_conc &lt; 0.00025 ? (ca_conc / 0.00025) : 1" min_conc="0" max_conc="0.001"/>

            <gate name="m" instances="1">
                
                <closed_state id="m0"/>
                <open_state id="m"/>
                     
                <transition name="alpha" from="m0" to="m" expr_form="generic" expr="v &lt; -10 ? speedup_factor * (2 / 37.95) * ( exp ( ((v + 50 ) / 11) - (( v + 53.5) / 27) ) ) : speedup_factor * 2 * (exp ( ( (-1 * v) - 53.5 ) / 27 ))" />
                    
                <transition name="beta" from="m" to="m0" expr_form="generic" expr="v &lt; -10 ? speedup_factor * 2 * (exp ( ( (-1 * v) - 53.5 ) / 27 )) - alpha : speedup_factor * 0.0" />
                    
            </gate>
            
        </current_voltage_relation> 
             
        <impl_prefs>
            <comment>This will ensure more points in calcium dependence table</comment>
            <table_settings max_v="140" min_v="-120" table_divisions="1041"/>
        </impl_prefs>
        
    </channel_type>
</channelml>

