<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 3/18/2013 at 3:27:44 at PM using:
COR (0.9.31.1409)
Copyright 2002-2013 Dr Alan Garny
http://cor.physiol.ox.ac.uk/ - cor@physiol.ox.ac.uk
CellML 1.0 was used to generate this model
http://www.cellml.org/
-->
<model name="francis_miriam_middleton_2013_realistic" cmeta:id="francis_miriam_middleton_2013_realistic" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<units name="ms">
<unit units="second" prefix="milli"/>
</units>
<units name="mM">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="mM4">
<unit units="mM" exponent="4"/>
</units>
<units name="cm">
<unit units="metre" prefix="centi"/>
</units>
<units name="per_cm">
<unit units="metre" prefix="centi" exponent="-1"/>
</units>
<units name="mV">
<unit units="volt" prefix="milli"/>
</units>
<units name="pl">
<unit units="litre" prefix="pico"/>
</units>
<units name="pF">
<unit units="farad" prefix="pico"/>
</units>
<units name="pF_per_sqcm">
<unit units="farad" prefix="pico"/>
<unit units="cm" exponent="-2"/>
</units>
<units name="sqcm">
<unit units="cm" exponent="2"/>
</units>
<units name="pA">
<unit units="ampere" prefix="pico"/>
</units>
<units name="coulomb_per_mole">
<unit units="coulomb"/>
<unit units="mole" exponent="-1"/>
</units>
<units name="mJ_per_mole_kelvin">
<unit units="joule" prefix="milli"/>
<unit units="mole" exponent="-1"/>
<unit units="kelvin" exponent="-1"/>
</units>
<units name="mM_per_ms">
<unit units="mM"/>
<unit units="ms" exponent="-1"/>
</units>
<units name="per_mMsq_ms">
<unit units="mM" exponent="-2"/>
<unit units="ms" exponent="-1"/>
</units>
<units name="per_mM_ms">
<unit units="mM" exponent="-1"/>
<unit units="ms" exponent="-1"/>
</units>
<units name="per_ms">
<unit units="ms" exponent="-1"/>
</units>
<units name="pA_per_mM">
<unit units="pA"/>
<unit units="mM" exponent="-1"/>
</units>
<units name="nS">
<unit units="pA"/>
<unit units="mV" exponent="-1"/>
</units>
<component name="environment">
<variable name="t" units="ms" public_interface="out"/>
<variable name="R" units="mJ_per_mole_kelvin" initial_value="8314.472"/>
<variable name="T" units="kelvin" initial_value="310.15"/>
<variable name="F" units="coulomb_per_mole" initial_value="96485.30929" public_interface="out"/>
<variable name="V_tau" units="mV" public_interface="out"/>
<variable name="Ca_o" units="mM" initial_value="1.8" public_interface="out"/>
<variable name="Na_o" units="mM" initial_value="137" public_interface="out"/>
<variable name="K_o" units="mM" initial_value="5.4" public_interface="out"/>
<variable name="vol_cyt" units="pl" public_interface="out"/>
<variable name="vol_pmu" units="pl" initial_value="5" public_interface="out"/>
<variable name="fr_cyt" units="dimensionless" initial_value="0.5"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>vol_cyt</ci>
<apply>
<times/>
<ci>fr_cyt</ci>
<ci>vol_pmu</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_tau</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="Membrane">
<variable name="t" units="ms" public_interface="in"/>
<variable name="V" units="mV" initial_value="-49.42" public_interface="out"/>
<variable name="V_tau" units="mV" public_interface="in"/>
<variable name="VD" units="dimensionless" public_interface="out"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="vol_cyt" units="pl" public_interface="in"/>
<variable name="vol_pmu" units="pl" public_interface="in"/>
<variable name="C_sp" units="pF_per_sqcm" initial_value="0.9e6"/>
<variable name="Ca_o" units="mM" public_interface="in"/>
<variable name="Na_o" units="mM" public_interface="in"/>
<variable name="K_o" units="mM" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="Na_i" units="mM" public_interface="in"/>
<variable name="K_i" units="mM" public_interface="in"/>
<variable name="V_Ca" units="dimensionless" public_interface="out"/>
<variable name="V_Na" units="dimensionless" public_interface="out"/>
<variable name="V_K" units="dimensionless" public_interface="out"/>
<variable name="A_pmu" units="sqcm"/>
<variable name="SVR_pmu" units="per_cm" initial_value="1.6667e4"/>
<variable name="J_Ca" units="mM_per_ms" public_interface="in"/>
<variable name="J_Na" units="mM_per_ms" public_interface="in"/>
<variable name="J_K" units="mM_per_ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>VD</ci>
<apply>
<divide/>
<ci>V</ci>
<ci>V_tau</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>A_pmu</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>SVR_pmu</ci>
<ci>vol_pmu</ci>
<cn type="e-notation" cellml:units="cm">1 <sep/>
-3</cn>
<cn type="e-notation" cellml:units="cm">1 <sep/>
-3</cn>
<cn type="e-notation" cellml:units="cm">1 <sep/>
-3</cn>
</apply>
<cn cellml:units="pl">1</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>V_Ca</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<ln/>
<apply>
<divide/>
<ci>Ca_o</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_Na</ci>
<apply>
<ln/>
<apply>
<divide/>
<ci>Na_o</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_K</ci>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_o</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>V</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>F</ci>
<ci>vol_cyt</ci>
</apply>
<apply>
<times/>
<ci>C_sp</ci>
<ci>A_pmu</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>J_Na</ci>
<ci>J_K</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>J_Ca</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Cytosol">
<variable name="t" units="ms" public_interface="in"/>
<variable name="Ca_i" units="mM" initial_value="0.000188" public_interface="out"/>
<variable name="Na_i" units="mM" initial_value="4.6876" public_interface="out"/>
<variable name="K_i" units="mM" initial_value="126.05893" public_interface="out"/>
<variable name="atp" units="mM" initial_value="2" public_interface="out"/>
<variable name="J_Ca" units="mM_per_ms" public_interface="in"/>
<variable name="J_Na" units="mM_per_ms" public_interface="in"/>
<variable name="J_K" units="mM_per_ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>Ca_i</ci>
</apply>
<ci>J_Ca</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>Na_i</ci>
</apply>
<ci>J_Na</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>K_i</ci>
</apply>
<ci>J_K</ci>
</apply>
</math>
</component>
<component name="calcium_dynamics">
<variable name="J_Ca" units="mM_per_ms" public_interface="out"/>
<variable name="J_ca" units="mM_per_ms"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="I_CaL" units="pA" public_interface="in"/>
<variable name="I_pmca" units="pA" public_interface="in"/>
<variable name="I_xm" units="pA" public_interface="in"/>
<variable name="vol_cyt" units="pl" public_interface="in"/>
<variable name="J_calb" units="mM_per_ms" public_interface="in"/>
<variable name="J_cam" units="mM_per_ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_Ca</ci>
<apply>
<minus/>
<ci>J_ca</ci>
<apply>
<plus/>
<ci>J_calb</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">4</cn>
<ci>J_cam</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_ca</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
<ci>vol_cyt</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_CaL</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>I_pmca</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>I_xm</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_dynamics">
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="J_Na" units="mM_per_ms" public_interface="out"/>
<variable name="I_Na" units="pA" public_interface="in"/>
<variable name="I_Nalk" units="pA" public_interface="in"/>
<variable name="I_NaHCN" units="pA" public_interface="in"/>
<variable name="I_xm" units="pA" public_interface="in"/>
<variable name="I_nk" units="pA" public_interface="in"/>
<variable name="vol_cyt" units="pl" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<times/>
<ci>F</ci>
<ci>vol_cyt</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>I_nk</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>I_xm</ci>
</apply>
<ci>I_Na</ci>
<ci>I_Nalk</ci>
<ci>I_NaHCN</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="potassium_dynamics">
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="J_K" units="mM_per_ms" public_interface="out"/>
<variable name="vol_cyt" units="pl" public_interface="in"/>
<variable name="I_K" units="pA"/>
<variable name="I_Ksk" units="pA" public_interface="in"/>
<variable name="I_Kdr" units="pA" public_interface="in"/>
<variable name="I_Kir" units="pA" public_interface="in"/>
<variable name="I_nk" units="pA" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_K</ci>
<apply>
<plus/>
<ci>I_Ksk</ci>
<ci>I_Kdr</ci>
<ci>I_Kir</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<times/>
<ci>F</ci>
<ci>vol_cyt</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>I_K</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>I_nk</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_buffer_dynamics">
<variable name="t" units="ms" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="Calb" units="mM" initial_value="0.0026"/>
<variable name="Cam" units="mM" initial_value="0.0222"/>
<variable name="Calbtot" units="mM" initial_value="0.005"/>
<variable name="Camtot" units="mM" initial_value="0.0235"/>
<variable name="CaCalb" units="mM"/>
<variable name="CaCam" units="mM" public_interface="out"/>
<variable name="J_calb" units="mM_per_ms" public_interface="out"/>
<variable name="J_cam" units="mM_per_ms" public_interface="out"/>
<variable name="kcal_1" units="per_mM_ms" initial_value="10"/>
<variable name="kcal_2" units="per_ms" initial_value="2e-3"/>
<variable name="kcam_cb" units="per_ms"/>
<variable name="kcam_cd" units="per_ms" initial_value="0.003"/>
<variable name="kcam_nb" units="per_ms"/>
<variable name="kcam_nd" units="per_ms" initial_value="3"/>
<variable name="alpha_cam" units="per_ms"/>
<variable name="beta_cam" units="per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_calb</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>kcal_1</ci>
<ci>Calb</ci>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<ci>kcal_2</ci>
<ci>CaCalb</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_cam</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_cam</ci>
<ci>Cam</ci>
</apply>
<apply>
<times/>
<ci>beta_cam</ci>
<ci>CaCam</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>kcam_cb</ci>
<apply>
<times/>
<cn cellml:units="per_mMsq_ms">12000</cn>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>kcam_nb</ci>
<apply>
<times/>
<cn type="e-notation" cellml:units="per_mMsq_ms">3.7 <sep/>
6</cn>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_cam</ci>
<apply>
<times/>
<ci>kcam_cb</ci>
<ci>kcam_nb</ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>kcam_cb</ci>
<ci>kcam_nd</ci>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>kcam_cd</ci>
<ci>kcam_nd</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_cam</ci>
<apply>
<times/>
<ci>kcam_cd</ci>
<ci>kcam_nd</ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>kcam_cb</ci>
<ci>kcam_nd</ci>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>kcam_cd</ci>
<ci>kcam_nd</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CaCalb</ci>
<apply>
<minus/>
<ci>Calbtot</ci>
<ci>Calb</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>CaCam</ci>
<apply>
<minus/>
<ci>Camtot</ci>
<ci>Cam</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>Calb</ci>
</apply>
<apply>
<minus/>
<ci>J_calb</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>Cam</ci>
</apply>
<apply>
<minus/>
<ci>J_cam</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel">
<variable name="t" units="ms" public_interface="in"/>
<variable name="V" units="mV" public_interface="in"/>
<variable name="VD" units="dimensionless" public_interface="in"/>
<variable name="V_Ca" units="dimensionless" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="Ca_o" units="mM" public_interface="in"/>
<variable name="I_CaL" units="pA" public_interface="out"/>
<variable name="m_cal" units="dimensionless" initial_value="0.006271"/>
<variable name="h_cal" units="dimensionless"/>
<variable name="g_cal" units="pA_per_mM" initial_value="2101.2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_CaL</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>g_cal</ci>
<ci>m_cal</ci>
<ci>h_cal</ci>
<apply>
<root/>
<apply>
<times/>
<ci>Ca_i</ci>
<ci>Ca_o</ci>
</apply>
</apply>
<apply>
<sinh/>
<apply>
<minus/>
<ci>VD</ci>
<ci>V_Ca</ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<sinh/>
<ci>VD</ci>
</apply>
<ci>VD</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>m_cal</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">15</cn>
</apply>
</apply>
<cn cellml:units="mV">7</cn>
</apply>
</apply>
</apply>
</apply>
<ci>m_cal</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="ms">7.68</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">65</cn>
</apply>
<cn cellml:units="mV">17.33</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="ms">0.7231</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>h_cal</ci>
<apply>
<divide/>
<cn cellml:units="mM">0.00045</cn>
<apply>
<plus/>
<cn cellml:units="mM">0.00045</cn>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="transient_Na_channel">
<variable name="t" units="ms" public_interface="in"/>
<variable name="VD" units="dimensionless" public_interface="in"/>
<variable name="V_Na" units="dimensionless" public_interface="in"/>
<variable name="Na_i" units="mM" public_interface="in"/>
<variable name="Na_o" units="mM" public_interface="in"/>
<variable name="I_Na" units="pA" public_interface="out"/>
<variable name="g_na" units="pA_per_mM" initial_value="907.68"/>
<variable name="O_na" units="dimensionless"/>
<variable name="m_na" units="dimensionless" initial_value="0.0952"/>
<variable name="h_na" units="dimensionless" initial_value="0.1848"/>
<variable name="A_mna" units="per_ms" initial_value="1.9651"/>
<variable name="B_mna" units="per_ms" initial_value="0.0424"/>
<variable name="A_hna" units="per_ms" initial_value="9.566e-5"/>
<variable name="B_hna" units="per_ms" initial_value="0.5296"/>
<variable name="za_mna" units="dimensionless" initial_value="1.7127"/>
<variable name="zb_mna" units="dimensionless" initial_value="1.5581"/>
<variable name="za_hna" units="dimensionless" initial_value="-2.4317"/>
<variable name="zb_hna" units="dimensionless" initial_value="-1.1868"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>O_na</ci>
<apply>
<times/>
<apply>
<power/>
<ci>m_na</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>h_na</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>m_na</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>A_mna</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>za_mna</ci>
<ci>VD</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>m_na</ci>
</apply>
</apply>
<apply>
<times/>
<ci>B_mna</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>zb_mna</ci>
</apply>
<ci>VD</ci>
</apply>
</apply>
<ci>m_na</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>h_na</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>A_hna</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>za_hna</ci>
<ci>VD</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>h_na</ci>
</apply>
</apply>
<apply>
<times/>
<ci>B_hna</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>zb_hna</ci>
</apply>
<ci>VD</ci>
</apply>
</apply>
<ci>h_na</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Na</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>g_na</ci>
<ci>O_na</ci>
<apply>
<root/>
<apply>
<times/>
<ci>Na_i</ci>
<ci>Na_o</ci>
</apply>
</apply>
<apply>
<sinh/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<minus/>
<ci>VD</ci>
<ci>V_Na</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<sinh/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>VD</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>VD</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Leak_Na_channel">
<variable name="t" units="ms" public_interface="in"/>
<variable name="VD" units="dimensionless" public_interface="in"/>
<variable name="V_Na" units="dimensionless" public_interface="in"/>
<variable name="Na_i" units="mM" public_interface="in"/>
<variable name="Na_o" units="mM" public_interface="in"/>
<variable name="I_Nalk" units="pA" public_interface="out"/>
<variable name="g_nalk" units="pA_per_mM" initial_value="0.0053"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_Nalk</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>g_nalk</ci>
<apply>
<root/>
<apply>
<times/>
<ci>Na_i</ci>
<ci>Na_o</ci>
</apply>
</apply>
<apply>
<sinh/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<minus/>
<ci>VD</ci>
<ci>V_Na</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<sinh/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>VD</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>VD</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="HCN_channel">
<variable name="t" units="ms" public_interface="in"/>
<variable name="V" units="mV" public_interface="in"/>
<variable name="VD" units="dimensionless" public_interface="in"/>
<variable name="V_Na" units="dimensionless" public_interface="in"/>
<variable name="V_K" units="dimensionless" public_interface="in"/>
<variable name="Na_i" units="mM" public_interface="in"/>
<variable name="Na_o" units="mM" public_interface="in"/>
<variable name="K_i" units="mM" public_interface="in"/>
<variable name="K_o" units="mM" public_interface="in"/>
<variable name="I_NaHCN" units="pA" public_interface="out"/>
<variable name="g_nahcn" units="pA_per_mM" initial_value="51.1"/>
<variable name="O_hcn" units="dimensionless" initial_value="0.003" public_interface="out"/>
<variable name="kf_hcn" units="per_ms"/>
<variable name="kr_hcn" units="per_ms"/>
<variable name="kf_free" units="per_ms"/>
<variable name="kr_free" units="per_ms"/>
<variable name="kf_bnd" units="per_ms"/>
<variable name="kr_bnd" units="per_ms"/>
<variable name="P_c" units="dimensionless"/>
<variable name="P_o" units="dimensionless"/>
<variable name="cAMP" units="mM" initial_value="1e-5"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>O_hcn</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>kf_hcn</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>O_hcn</ci>
</apply>
</apply>
<apply>
<times/>
<ci>kr_hcn</ci>
<ci>O_hcn</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>kf_hcn</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>kf_free</ci>
<ci>P_c</ci>
</apply>
<apply>
<times/>
<ci>kf_bnd</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_c</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>kr_hcn</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>kr_free</ci>
<ci>P_o</ci>
</apply>
<apply>
<times/>
<ci>kr_bnd</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_o</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>kf_free</ci>
<apply>
<divide/>
<cn cellml:units="per_ms">0.006</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">87.7</cn>
</apply>
<cn cellml:units="mV">6.45</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>kr_free</ci>
<apply>
<divide/>
<cn cellml:units="per_ms">0.08</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">51.7</cn>
</apply>
</apply>
<cn cellml:units="mV">7</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>kf_bnd</ci>
<apply>
<divide/>
<cn cellml:units="per_ms">0.0268</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">94.2</cn>
</apply>
<cn cellml:units="mV">13.3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>kr_bnd</ci>
<apply>
<divide/>
<cn cellml:units="per_ms">0.08</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">35.5</cn>
</apply>
</apply>
<cn cellml:units="mV">7</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_c</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>cAMP</ci>
<cn type="e-notation" cellml:units="mM">1.163 <sep/>
-3</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_o</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>cAMP</ci>
<cn type="e-notation" cellml:units="mM">1.45 <sep/>
-5</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NaHCN</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>g_nahcn</ci>
<ci>O_hcn</ci>
<apply>
<root/>
<apply>
<times/>
<ci>Na_i</ci>
<ci>Na_o</ci>
</apply>
</apply>
<apply>
<sinh/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<minus/>
<ci>VD</ci>
<ci>V_Na</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<sinh/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>VD</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>VD</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SK_K_channel">
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="K_i" units="mM" public_interface="in"/>
<variable name="K_o" units="mM" public_interface="in"/>
<variable name="VD" units="dimensionless" public_interface="in"/>
<variable name="V_K" units="dimensionless" public_interface="in"/>
<variable name="O_sk" units="dimensionless"/>
<variable name="I_Ksk" units="pA" public_interface="out"/>
<variable name="g_ksk" units="pA_per_mM" initial_value="2.2515"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>O_sk</ci>
<apply>
<divide/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">4.2</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<cn type="e-notation" cellml:units="mM">0.35 <sep/>
-3</cn>
<cn cellml:units="dimensionless">4.2</cn>
</apply>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">4.2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Ksk</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>g_ksk</ci>
<ci>O_sk</ci>
<apply>
<root/>
<apply>
<times/>
<ci>K_i</ci>
<ci>K_o</ci>
</apply>
</apply>
<apply>
<sinh/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<apply>
<minus/>
<ci>VD</ci>
<ci>V_K</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<sinh/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>VD</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>VD</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="DR_K_channel">
<variable name="t" units="ms" public_interface="in"/>
<variable name="V" units="mV" public_interface="in"/>
<variable name="V_tau" units="mV" public_interface="in"/>
<variable name="V_K" units="dimensionless" public_interface="in"/>
<variable name="O_kdr" units="dimensionless"/>
<variable name="m_kdr" units="dimensionless" initial_value="0.0932"/>
<variable name="I_Kdr" units="pA" public_interface="out"/>
<variable name="g_kdr" units="nS" initial_value="31.237"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_Kdr</ci>
<apply>
<times/>
<ci>g_kdr</ci>
<ci>O_kdr</ci>
<apply>
<minus/>
<ci>V</ci>
<apply>
<times/>
<ci>V_K</ci>
<ci>V_tau</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>O_kdr</ci>
<apply>
<power/>
<ci>m_kdr</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>m_kdr</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">25</cn>
</apply>
</apply>
<cn cellml:units="mV">12</cn>
</apply>
</apply>
</apply>
</apply>
<ci>m_kdr</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="ms">18</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">39</cn>
</apply>
<cn cellml:units="mV">8</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="ms">1</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="IR_K_channel">
<variable name="V" units="mV" public_interface="in"/>
<variable name="V_tau" units="mV" public_interface="in"/>
<variable name="V_K" units="dimensionless" public_interface="in"/>
<variable name="O_kir" units="dimensionless"/>
<variable name="I_Kir" units="pA" public_interface="out"/>
<variable name="g_kir" units="nS" initial_value="13.816"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_Kir</ci>
<apply>
<times/>
<ci>g_kir</ci>
<ci>O_kir</ci>
<apply>
<minus/>
<ci>V</ci>
<apply>
<times/>
<ci>V_K</ci>
<ci>V_tau</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>O_kir</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">85</cn>
</apply>
<cn cellml:units="mV">12.1</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="PMCA">
<variable name="t" units="ms" public_interface="in"/>
<variable name="VD" units="dimensionless" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="Ca_o" units="mM" public_interface="in"/>
<variable name="atp" units="mM" public_interface="in"/>
<variable name="CaCam" units="mM" public_interface="in"/>
<variable name="I_pmca" units="pA" public_interface="out"/>
<variable name="y_pc" units="dimensionless" initial_value="0.483"/>
<variable name="K_pmca" units="pA"/>
<variable name="k_1pc" units="per_ms"/>
<variable name="k_2pc" units="per_ms" initial_value="0.001"/>
<variable name="k_3pc" units="per_ms" initial_value="0.001"/>
<variable name="k_4pc" units="per_ms" initial_value="1"/>
<variable name="P_E1Spc" units="dimensionless"/>
<variable name="P_E2Spc" units="dimensionless"/>
<variable name="P_E1pc" units="dimensionless"/>
<variable name="P_E2pc" units="dimensionless"/>
<variable name="alpha_pc" units="per_ms"/>
<variable name="beta_pc" units="per_ms"/>
<variable name="K_pci" units="mM"/>
<variable name="K_pco" units="mM" initial_value="2"/>
<variable name="k_pmca" units="dimensionless" initial_value="2.233"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_pmca</ci>
<apply>
<times/>
<ci>K_pmca</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_1pc</ci>
<ci>P_E1Spc</ci>
<ci>y_pc</ci>
</apply>
<apply>
<times/>
<ci>k_2pc</ci>
<ci>P_E2Spc</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y_pc</ci>
</apply>
</apply>
</apply>
<cn cellml:units="ms">1</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>K_pmca</ci>
<apply>
<times/>
<ci>k_pmca</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="pA">10.56</cn>
<ci>CaCam</ci>
</apply>
<apply>
<plus/>
<ci>CaCam</ci>
<cn cellml:units="mM">0.00005</cn>
</apply>
</apply>
<cn cellml:units="pA">1.2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>y_pc</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>beta_pc</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y_pc</ci>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_pc</ci>
<ci>y_pc</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_1pc</ci>
<apply>
<divide/>
<cn cellml:units="per_ms">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<cn cellml:units="mM">0.1</cn>
<ci>atp</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_E1Spc</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>K_pci</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_E1pc</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_E1Spc</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_E2Spc</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>K_pco</ci>
<ci>Ca_o</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_E2pc</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_E2Spc</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_pc</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_2pc</ci>
<ci>P_E2Spc</ci>
</apply>
<apply>
<times/>
<ci>k_4pc</ci>
<ci>P_E2pc</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_pc</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_1pc</ci>
<ci>P_E1Spc</ci>
</apply>
<apply>
<times/>
<ci>k_3pc</ci>
<ci>P_E1pc</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>K_pci</ci>
<apply>
<times/>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">173.6</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>CaCam</ci>
<cn cellml:units="mM">0.00005</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">6.4</cn>
</apply>
<cn type="e-notation" cellml:units="mM">1 <sep/>
-5</cn>
</apply>
</apply>
</math>
</component>
<component name="NaCa">
<variable name="Na_i" units="mM" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="Na_o" units="mM" public_interface="in"/>
<variable name="Ca_o" units="mM" public_interface="in"/>
<variable name="del" units="dimensionless" initial_value="0.35"/>
<variable name="I_xm" units="pA" public_interface="out"/>
<variable name="VD" units="dimensionless" public_interface="in"/>
<variable name="k_xm" units="pA" initial_value="0.0166"/>
<variable name="Dr" units="mM4"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_xm</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>k_xm</ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<times/>
<ci>del</ci>
<ci>VD</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>Na_o</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>Ca_i</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>del</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>VD</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>Dr</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Dr</ci>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="mM4">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.001</cn>
<apply>
<plus/>
<apply>
<times/>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>Ca_o</ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>Na_o</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ca_i</ci>
<cn cellml:units="mM">0.0069</cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_pump">
<variable name="t" units="ms" public_interface="in"/>
<variable name="Na_i" units="mM" public_interface="in"/>
<variable name="Na_o" units="mM" public_interface="in"/>
<variable name="Na_eff" units="mM"/>
<variable name="K_i" units="mM" public_interface="in"/>
<variable name="K_o" units="mM" public_interface="in"/>
<variable name="atp" units="mM" public_interface="in"/>
<variable name="VD" units="dimensionless" public_interface="in"/>
<variable name="I_nk" units="pA" public_interface="out"/>
<variable name="y_nk" units="dimensionless" initial_value="0.6213"/>
<variable name="alpha_nk" units="per_ms"/>
<variable name="beta_nk" units="per_ms"/>
<variable name="P_E1Snk" units="dimensionless"/>
<variable name="P_E2Snk" units="dimensionless"/>
<variable name="P_E1Dnk" units="dimensionless"/>
<variable name="P_E2Dnk" units="dimensionless"/>
<variable name="k_1nk" units="per_ms"/>
<variable name="k_2nk" units="per_ms" initial_value="0.04"/>
<variable name="k_3nk" units="per_ms" initial_value="0.01"/>
<variable name="k_4nk" units="per_ms" initial_value="0.165"/>
<variable name="K_nknai" units="mM" initial_value="4.05"/>
<variable name="K_nknao" units="mM" initial_value="69.8"/>
<variable name="K_nkki" units="mM" initial_value="32.88"/>
<variable name="K_nkko" units="mM" initial_value="0.258"/>
<variable name="k_nk" units="pA" initial_value="1085.7"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_nk</ci>
<apply>
<times/>
<ci>k_nk</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_1nk</ci>
<ci>P_E1Snk</ci>
<ci>y_nk</ci>
</apply>
<apply>
<times/>
<ci>k_2nk</ci>
<ci>P_E2Snk</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y_nk</ci>
</apply>
</apply>
</apply>
<cn cellml:units="ms">1</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>k_1nk</ci>
<apply>
<divide/>
<cn cellml:units="per_ms">0.37</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<cn cellml:units="mM">0.094</cn>
<ci>atp</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Na_eff</ci>
<apply>
<times/>
<ci>Na_o</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.82</cn>
</apply>
<ci>VD</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>y_nk</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>beta_nk</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y_nk</ci>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_nk</ci>
<ci>y_nk</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_nk</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_2nk</ci>
<ci>P_E2Snk</ci>
</apply>
<apply>
<times/>
<ci>k_4nk</ci>
<ci>P_E2Dnk</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_nk</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_1nk</ci>
<ci>P_E1Snk</ci>
</apply>
<apply>
<times/>
<ci>k_3nk</ci>
<ci>P_E1Dnk</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_E1Snk</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>K_nknai</ci>
<ci>Na_i</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>K_i</ci>
<ci>K_nkki</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_E1Dnk</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>K_nkki</ci>
<ci>K_i</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Na_i</ci>
<ci>K_nknai</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_E2Snk</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>K_nknao</ci>
<ci>Na_eff</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>K_o</ci>
<ci>K_nkko</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_E2Dnk</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<apply>
<divide/>
<ci>K_nkko</ci>
<ci>K_o</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Na_eff</ci>
<ci>K_nknao</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="Membrane">
<component_ref component="L_type_Ca_channel"/>
<component_ref component="transient_Na_channel"/>
<component_ref component="Leak_Na_channel"/>
<component_ref component="HCN_channel"/>
<component_ref component="SK_K_channel"/>
<component_ref component="DR_K_channel"/>
<component_ref component="IR_K_channel"/>
<component_ref component="sodium_pump"/>
<component_ref component="PMCA"/>
<component_ref component="NaCa"/>
</component_ref>
<component_ref component="Cytosol">
<component_ref component="calcium_dynamics">
<component_ref component="calcium_buffer_dynamics"/>
</component_ref>
<component_ref component="sodium_dynamics"/>
<component_ref component="potassium_dynamics"/>
</component_ref>
</group>
<connection>
<map_components component_1="environment" component_2="Membrane"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="V_tau" variable_2="V_tau"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="vol_cyt" variable_2="vol_cyt"/>
<map_variables variable_1="vol_pmu" variable_2="vol_pmu"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
<map_variables variable_1="K_o" variable_2="K_o"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Cytosol"/>
<map_variables variable_1="t" variable_2="t"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="Cytosol"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
</connection>
<connection>
<map_components component_1="environment" component_2="calcium_dynamics"/>
<map_variables variable_1="vol_cyt" variable_2="vol_cyt"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="calcium_dynamics"/>
<map_variables variable_1="J_Ca" variable_2="J_Ca"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="calcium_dynamics"/>
<map_variables variable_1="J_Ca" variable_2="J_Ca"/>
</connection>
<connection>
<map_components component_1="environment" component_2="sodium_dynamics"/>
<map_variables variable_1="vol_cyt" variable_2="vol_cyt"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="sodium_dynamics"/>
<map_variables variable_1="J_Na" variable_2="J_Na"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="sodium_dynamics"/>
<map_variables variable_1="J_Na" variable_2="J_Na"/>
</connection>
<connection>
<map_components component_1="environment" component_2="potassium_dynamics"/>
<map_variables variable_1="vol_cyt" variable_2="vol_cyt"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="potassium_dynamics"/>
<map_variables variable_1="J_K" variable_2="J_K"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="potassium_dynamics"/>
<map_variables variable_1="J_K" variable_2="J_K"/>
</connection>
<connection>
<map_components component_1="environment" component_2="calcium_buffer_dynamics"/>
<map_variables variable_1="t" variable_2="t"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="calcium_buffer_dynamics"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="calcium_dynamics" component_2="calcium_buffer_dynamics"/>
<map_variables variable_1="J_calb" variable_2="J_calb"/>
<map_variables variable_1="J_cam" variable_2="J_cam"/>
</connection>
<connection>
<map_components component_1="environment" component_2="L_type_Ca_channel"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="L_type_Ca_channel"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="VD" variable_2="VD"/>
<map_variables variable_1="V_Ca" variable_2="V_Ca"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="L_type_Ca_channel"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="calcium_dynamics" component_2="L_type_Ca_channel"/>
<map_variables variable_1="I_CaL" variable_2="I_CaL"/>
</connection>
<connection>
<map_components component_1="environment" component_2="PMCA"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="PMCA"/>
<map_variables variable_1="VD" variable_2="VD"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="PMCA"/>
<map_variables variable_1="atp" variable_2="atp"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="calcium_dynamics" component_2="PMCA"/>
<map_variables variable_1="I_pmca" variable_2="I_pmca"/>
</connection>
<connection>
<map_components component_1="calcium_buffer_dynamics" component_2="PMCA"/>
<map_variables variable_1="CaCam" variable_2="CaCam"/>
</connection>
<connection>
<map_components component_1="environment" component_2="NaCa"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="NaCa"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="NaCa"/>
<map_variables variable_1="VD" variable_2="VD"/>
</connection>
<connection>
<map_components component_1="calcium_dynamics" component_2="NaCa"/>
<map_variables variable_1="I_xm" variable_2="I_xm"/>
</connection>
<connection>
<map_components component_1="sodium_dynamics" component_2="NaCa"/>
<map_variables variable_1="I_xm" variable_2="I_xm"/>
</connection>
<connection>
<map_components component_1="environment" component_2="sodium_pump"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
<map_variables variable_1="K_o" variable_2="K_o"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="sodium_pump"/>
<map_variables variable_1="VD" variable_2="VD"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="sodium_pump"/>
<map_variables variable_1="atp" variable_2="atp"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
</connection>
<connection>
<map_components component_1="sodium_dynamics" component_2="sodium_pump"/>
<map_variables variable_1="I_nk" variable_2="I_nk"/>
</connection>
<connection>
<map_components component_1="potassium_dynamics" component_2="sodium_pump"/>
<map_variables variable_1="I_nk" variable_2="I_nk"/>
</connection>
<connection>
<map_components component_1="environment" component_2="transient_Na_channel"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="transient_Na_channel"/>
<map_variables variable_1="VD" variable_2="VD"/>
<map_variables variable_1="V_Na" variable_2="V_Na"/>
</connection>
<connection>
<map_components component_1="sodium_dynamics" component_2="transient_Na_channel"/>
<map_variables variable_1="I_Na" variable_2="I_Na"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="transient_Na_channel"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Leak_Na_channel"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="Leak_Na_channel"/>
<map_variables variable_1="VD" variable_2="VD"/>
<map_variables variable_1="V_Na" variable_2="V_Na"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="Leak_Na_channel"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
</connection>
<connection>
<map_components component_1="sodium_dynamics" component_2="Leak_Na_channel"/>
<map_variables variable_1="I_Nalk" variable_2="I_Nalk"/>
</connection>
<connection>
<map_components component_1="environment" component_2="HCN_channel"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
<map_variables variable_1="K_o" variable_2="K_o"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="HCN_channel"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="VD" variable_2="VD"/>
<map_variables variable_1="V_Na" variable_2="V_Na"/>
<map_variables variable_1="V_K" variable_2="V_K"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="HCN_channel"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
</connection>
<connection>
<map_components component_1="sodium_dynamics" component_2="HCN_channel"/>
<map_variables variable_1="I_NaHCN" variable_2="I_NaHCN"/>
</connection>
<connection>
<map_components component_1="environment" component_2="SK_K_channel"/>
<map_variables variable_1="K_o" variable_2="K_o"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="SK_K_channel"/>
<map_variables variable_1="VD" variable_2="VD"/>
<map_variables variable_1="V_K" variable_2="V_K"/>
</connection>
<connection>
<map_components component_1="Cytosol" component_2="SK_K_channel"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="potassium_dynamics" component_2="SK_K_channel"/>
<map_variables variable_1="I_Ksk" variable_2="I_Ksk"/>
</connection>
<connection>
<map_components component_1="environment" component_2="DR_K_channel"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="V_tau" variable_2="V_tau"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="DR_K_channel"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="V_K" variable_2="V_K"/>
</connection>
<connection>
<map_components component_1="potassium_dynamics" component_2="DR_K_channel"/>
<map_variables variable_1="I_Kdr" variable_2="I_Kdr"/>
</connection>
<connection>
<map_components component_1="environment" component_2="IR_K_channel"/>
<map_variables variable_1="V_tau" variable_2="V_tau"/>
</connection>
<connection>
<map_components component_1="Membrane" component_2="IR_K_channel"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="V_K" variable_2="V_K"/>
</connection>
<connection>
<map_components component_1="potassium_dynamics" component_2="IR_K_channel"/>
<map_variables variable_1="I_Kir" variable_2="I_Kir"/>
</connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<rdf:Description rdf:about="">
<dc:creator rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Francis</vCard:Family>
<vCard:Given>Febe</vCard:Given>
</vCard:N>
<vCard:EMAIL rdf:parseType="Resource">
<rdf:value>Febe.Francis@nuim.ie</rdf:value>
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
</vCard:EMAIL>
<vCard:ORG rdf:parseType="Resource">
<vCard:Orgname>National University of Ireland Maynooth</vCard:Orgname>
<vCard:Orgunit>Hamilton Institute</vCard:Orgunit>
</vCard:ORG>
</dc:creator>
<dcterms:created rdf:parseType="Resource">
<dcterms:W3CDTF>2013-03-18</dcterms:W3CDTF>
</dcterms:created>
</rdf:Description>
<rdf:Description rdf:about="#aslanidi_2009">
<cmeta:species>Mouse</cmeta:species>
<cmeta:bio_entity>Substantia nigra</cmeta:bio_entity>
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>electrophysiology</rdf:li>
<rdf:li>neuron</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<bqs:reference rdf:parseType="Resource">
<bqs:JournalArticle rdf:parseType="Resource">
<dc:creator>
<rdf:Seq>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Francis</vCard:Family>
<vCard:Given>F</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Garcia</vCard:Family>
<vCard:Given>M</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Middleton</vCard:Family>
<vCard:Given>R</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
</rdf:Seq>
</dc:creator>
<dc:title>
A single Compartment model of pacemaking in dissasociated Substantia nigra neurons. Stability and energy analysis
</dc:title>
<bqs:Journal rdf:parseType="Resource">
<dc:title>Journal of Computational Neuroscience</dc:title>
</bqs:Journal>
<bqs:volume>xx</bqs:volume>
<bqs:first_page>xx</bqs:first_page>
<bqs:last_page>xx</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
</model>
