- 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/IK_HH_open.cellml
<?xml version='1.0' encoding='UTF-8'?>
<!-- The original open K channel model. -->
<model name="IK_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_K" name="E_K"/>
</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="IK_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="IK_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="20" name="Ko" private_interface="out" units="mM"/>
<variable initial_value="489.5" name="Ki" private_interface="out" units="mM"/>
<variable initial_value="36" name="gK_max" private_interface="out" units="mS_per_cmsq"/>
<variable initial_value="0.3" name="A_cap" units="cmsq"/>
<variable name="g_K" 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_K" private_interface="in" units="pA"/>
<variable name="I_K_fA" units="fA"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>g_K</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>gK_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_K_fA</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>I_K</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="IK_HH_open">
<component_ref component="Ionic_currents"/>
<component_ref component="E_K"/>
</component_ref>
</group>
<connection>
<map_components component_1="IK_HH_open" component_2="Ionic_currents"/>
<map_variables variable_1="I_K" variable_2="I"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="g_K" variable_2="g_max"/>
<map_variables variable_1="PO" variable_2="PO"/>
</connection>
<connection>
<map_components component_1="model_para" component_2="IK_HH_open"/>
<map_variables variable_1="T" variable_2="T"/>
</connection>
<connection>
<map_components component_1="IK_HH_open" component_2="E_K"/>
<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="Ki" variable_2="Ki"/>
<map_variables variable_1="Ko" variable_2="Ko"/>
</connection>
<connection>
<map_components component_1="E_K" component_2="Ionic_currents"/>
<map_variables variable_1="E_K" variable_2="E"/>
</connection>
<connection>
<map_components component_1="IK_HH_open" component_2="constants_BG"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
</model>
