Location: Hodgkin & Huxley (1952) model @ ab2962b87966 / BG / INa_HH_open.cellml

Author:
WeiweiAi <wai484@aucklanduni.ac.nz>
Date:
2022-03-31 10:48:40+13:00
Desc:
Add open channel IV of HH models; Add simulation and plot python scripts
Permanent Source URI:
https://models.fieldml.org/workspace/64f/rawfile/ab2962b8796666ad9938a2577611c954f006b5fd/BG/INa_HH_open.cellml

<?xml version='1.0' encoding='UTF-8'?>
<!-- The original open Na channel model. -->
<model name="INa_HH_open" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
    <import xlink:href="../cellLib/Components/units.cellml">
        <units name="mS" units_ref="mS"/>
        <units name="nS" units_ref="nS"/>
        <units name="mV" units_ref="mV"/>
        <units name="mM" units_ref="mM"/>
        <units name="mS_per_cmsq" units_ref="mS_per_cmsq"/>
        <units name="pA" units_ref="pA"/>
        <units name="cmsq" units_ref="cmsq"/>
        <units name="mV_per_ms" units_ref="mV_per_ms"/>
        <units name="J_per_K_mol" units_ref="J_per_K_mol"/>
        <units name="C_per_mmol" units_ref="C_per_mmol"/>
    </import>
    <import xlink:href="../cellLib/BG/units_BG.cellml">
        <units name="fA" units_ref="fA"/>
    </import>
    <import xlink:href="../cellLib/BG/constants_BG.cellml">
        <component component_ref="constants_BG" name="constants_BG"/>
    </import>
    <import xlink:href="default_values_BG.cellml">
        <component component_ref="model_para" name="model_para"/>
    </import>
    <import xlink:href="../cellLib/Components/time.cellml">
        <component component_ref="time_ms" name="time"/>
    </import>
    <import xlink:href="../cellLib/Components/Ionic_currents.cellml">
        <component component_ref="Ionic_currents" name="Ionic_currents"/>
    </import>
    <import xlink:href="../cellLib/Components/E_Nernst.cellml">
        <component component_ref="E_Na" name="E_Na"/>
    </import>
    <import xlink:href="../cellLib/Protocols/Voltage_ramp_protocol.cellml">
        <component component_ref="vramp_protocol_ms" name="vramp_protocol"/>
    </import>
    <!-- Connect the parameters to the protocol definition -->
    <component name="Vstim_para">
        <variable initial_value="-100" name="V_init" public_interface="out" units="mV"/>
        <variable initial_value="1" name="V_rate" public_interface="out" units="mV_per_ms"/>
    </component>
    <connection>
        <map_components component_1="Vstim_para" component_2="vramp_protocol"/>
        <map_variables variable_1="V_init" variable_2="V_init"/>
        <map_variables variable_1="V_rate" variable_2="V_rate"/>
    </connection>
    <connection>
        <map_components component_1="vramp_protocol" component_2="INa_HH_open"/>
        <map_variables variable_1="V" variable_2="V"/>
    </connection>
    <connection>
        <map_components component_1="time" component_2="vramp_protocol"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <component name="INa_HH_open">
        <!--Parameters-->
        <variable name="T" private_interface="out" public_interface="in" units="kelvin"/>
        <variable name="R" private_interface="out" public_interface="in" units="J_per_K_mol"/>
        <variable name="F" private_interface="out" public_interface="in" units="C_per_mmol"/>
        <variable initial_value="484" name="Nao" private_interface="out" units="mM"/>
        <variable initial_value="60.7" name="Nai" private_interface="out" units="mM"/>
        <variable initial_value="120" name="gNa_max" private_interface="out" units="mS_per_cmsq"/>
        <variable initial_value="0.3" name="A_cap" units="cmsq"/>
        <variable name="g_Na" private_interface="out" units="nS"/>
        <variable initial_value="1" name="PO" private_interface="out" units="dimensionless"/>
        <!--Inputs-->
        <variable name="V" private_interface="out" public_interface="in" units="mV"/>
        <!--Outputs-->
        <variable name="I_Na" private_interface="in" units="pA"/>
        <variable name="I_Na_fA" units="fA"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>g_Na</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>gNa_max</ci>
                        <ci>A_cap</ci>
                        <cn cellml:units="nS" type="e-notation">1<sep/>6</cn>
                    </apply>
                    <cn cellml:units="mS">1</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_Na_fA</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>I_Na</ci>
                        <cn cellml:units="fA">1000</cn>
                    </apply>
                    <cn cellml:units="pA">1</cn>
                </apply>
            </apply>
        </math>
    </component>
    <group>
        <relationship_ref relationship="encapsulation"/>
        <component_ref component="INa_HH_open">
            <component_ref component="Ionic_currents"/>
            <component_ref component="E_Na"/>
        </component_ref>
    </group>
    <connection>
        <map_components component_1="INa_HH_open" component_2="Ionic_currents"/>
        <map_variables variable_1="I_Na" variable_2="I"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="g_Na" variable_2="g_max"/>
        <map_variables variable_1="PO" variable_2="PO"/>
    </connection>
    <connection>
        <map_components component_1="model_para" component_2="INa_HH_open"/>
        <map_variables variable_1="T" variable_2="T"/>
    </connection>
    <connection>
        <map_components component_1="INa_HH_open" component_2="E_Na"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="Nai" variable_2="Nai"/>
        <map_variables variable_1="Nao" variable_2="Nao"/>
    </connection>
    <connection>
        <map_components component_1="E_Na" component_2="Ionic_currents"/>
        <map_variables variable_1="E_Na" variable_2="E"/>
    </connection>
    <connection>
        <map_components component_1="INa_HH_open" component_2="constants_BG"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
</model>