Location: BG_LCC @ 2e4d8c4eb984 / BG_LCC.cellml

Author:
Shelley Fong <s.fong@auckland.ac.nz>
Date:
2021-11-08 11:15:30+13:00
Desc:
Change zCa to z_rCa
Permanent Source URI:
https://models.physiomeproject.org/workspace/6d7/rawfile/2e4d8c4eb9847def4a75dd06f6d1b0a9ffdc6146/BG_LCC.cellml

<?xml version='1.0' encoding='UTF-8'?>
<!-- L_type Ca channel with gates
 same as Pan and LRd
 flux of charge (dq_membrane/dt) = sum of all currents-->
<!-- 5 Aug repurposing this code to use as general, modular LCC channel
 6 Aug updating parameters with values for an isolated LCC system (no other channels present when finding BG parameters)-->
<!-- 4 Nov adding pulse to represent false action potential as a stimulus
 adding membrane potential contribution from gate variables-->
<model name="LCC_module" 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="units_and_constants/units_BG.cellml">
        <units name="mM" units_ref="mM"/>
        <units name="fmol" units_ref="fmol"/>
        <units name="per_fmol" units_ref="per_fmol"/>
        <units name="J_per_mol" units_ref="J_per_mol"/>
        <units name="J_per_K_per_mol" units_ref="J_per_K_per_mol"/>
        <units name="fmol_per_sec" units_ref="fmol_per_sec"/>
        <units name="C_per_mol" units_ref="C_per_mol"/>
        <units name="J_per_C" units_ref="J_per_C"/>
        <units name="microm3" units_ref="microm3"/>
        <units name="fF" units_ref="fF"/>
        <units name="fC" units_ref="fC"/>
        <units name="fA" units_ref="fA"/>
        <units name="per_second" units_ref="per_second"/>
        <units name="millivolt" units_ref="millivolt"/>
    </import>
    <import xlink:href="units_and_constants/constants_BG.cellml">
        <component component_ref="constants" name="constants"/>
    </import>
    <component name="environment">
        <variable name="time" public_interface="out" units="second"/>
        <variable name="Vmem" units="volt"/>
        <variable initial_value="-13039" name="q_membrane" public_interface="out" units="fC"/>
        <variable initial_value="1.381e5" name="C_m" public_interface="out" units="fF"/>
        <!--        var R: J_per_K_per_mol {pub: in};
        var T: kelvin {pub: in};
        var RT: J_per_mol {pub: out};-->
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable initial_value="9.3276" name="q_Cao_init" units="fmol"/>
        <variable initial_value="0.00456" name="q_Cai_init" units="fmol"/>
        <variable initial_value="27.9828" name="q_Ko_init" units="fmol"/>
        <variable initial_value="5510" name="q_Ki_init" units="fmol"/>
        <variable name="v_CaGHK_i" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_CaGHK_o" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_KGHK_i" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_KGHK_o" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_E_LCC" public_interface="in" units="fA"/>
        <variable initial_value="2" name="zCa" public_interface="out" units="dimensionless"/>
        <variable initial_value="1" name="zK" public_interface="out" units="dimensionless"/>
        <!--        ode(q_Cai, time) = v_CaGHK_i;
        ode(q_Cao, time) = v_CaGHK_o;
        ode(q_Ki, time) = v_KGHK_i;
        ode(q_Ko, time) = v_KGHK_o;-->
        <variable name="I_pulse" public_interface="out" units="fA"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_pulse</ci>
                <piecewise>
                    <piece>
                        <cn cellml:units="fA" type="e-notation">4<sep/>6</cn>
                        <apply>
                            <and/>
                            <apply>
                                <gt/>
                                <ci>time</ci>
                                <cn cellml:units="second" type="e-notation">1<sep/>-4</cn>
                            </apply>
                            <apply>
                                <lt/>
                                <ci>time</ci>
                                <cn cellml:units="second" type="e-notation">2<sep/>-4</cn>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <cn cellml:units="fA">0</cn>
                    </otherwise>
                </piecewise>
            </apply>
            <!--            case (time > 2e-4{second}) and (time < 2.1e-4{second}):
                4e6{fA};-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_membrane</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>v_E_LCC</ci>
                    </apply>
                    <ci>I_pulse</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Vmem</ci>
                <apply>
                    <divide/>
                    <ci>q_membrane</ci>
                    <ci>C_m</ci>
                </apply>
            </apply>
        </math>
        <!-- global value-->
        <variable name="q_Cao" public_interface="out" units="fmol"/>
        <variable name="q_Cai" public_interface="out" units="fmol"/>
        <variable name="q_Ko" public_interface="out" units="fmol"/>
        <variable name="q_Ki" public_interface="out" units="fmol"/>
        <!-- from submodule-->
        <variable name="q_Cao_mLCC" public_interface="in" units="fmol"/>
        <variable name="q_Cai_mLCC" public_interface="in" units="fmol"/>
        <variable name="q_Ko_mLCC" public_interface="in" units="fmol"/>
        <variable name="q_Ki_mLCC" public_interface="in" units="fmol"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>q_Cao</ci>
                <apply>
                    <plus/>
                    <ci>q_Cao_mLCC</ci>
                    <ci>q_Cao_init</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>q_Cai</ci>
                <apply>
                    <plus/>
                    <ci>q_Cai_mLCC</ci>
                    <ci>q_Cai_init</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>q_Ko</ci>
                <apply>
                    <plus/>
                    <ci>q_Ko_mLCC</ci>
                    <ci>q_Ko_init</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>q_Ki</ci>
                <apply>
                    <plus/>
                    <ci>q_Ki_mLCC</ci>
                    <ci>q_Ki_init</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="LCC_parameters">
        <variable initial_value="115.689" name="kappa_Ca1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="216.336" name="kappa_Ca2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.160891" name="kappa_K1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.300862" name="kappa_K2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="132.798" name="kappa_d000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="3.27634" name="kappa_d010" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6.12666" name="kappa_d020" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.525901" name="kappa_d001" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.0129748" name="kappa_d011" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.0242625" name="kappa_d021" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.244648" name="kappa_f1_000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.21486" name="kappa_f1_100" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.000968841" name="kappa_f1_001" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00481102" name="kappa_f1_101" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.448058" name="kappa_f2_000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="2.22495" name="kappa_f2_100" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00177438" name="kappa_f2_001" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00881111" name="kappa_f2_101" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1258.66" name="kappa_f3_010" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6250.17" name="kappa_f3_110" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="4.98447" name="kappa_f3_011" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="24.7516" name="kappa_f3_111" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="180068" name="kappa_fCa000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="894171" name="kappa_fCa100" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="4442.55" name="kappa_fCa010" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="22060.6" name="kappa_fCa110" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="8307.45" name="kappa_fCa020" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="41252.7" name="kappa_fCa120" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.0102432" name="K_Cai" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0751143" name="K_Cao" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.00578984" name="K_Ki" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0424574" name="K_Ko" public_interface="out" units="per_fmol"/>
        <variable initial_value="3.66543" name="K_000_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="148.569" name="K_010_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="79.4498" name="K_020_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.738142" name="K_100_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="29.9187" name="K_110_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="15.9996" name="K_120_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="925.578" name="K_001_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="37516" name="K_011_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="20062.3" name="K_021_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="186.392" name="K_101_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="7554.95" name="K_111_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="4040.14" name="K_121_LCC" public_interface="out" units="per_fmol"/>
    </component>
    <component name="LCC">
        <variable name="time" public_interface="in" units="second"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <variable name="RT" public_interface="out" units="J_per_mol"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable name="C_m" public_interface="in" units="fF"/>
        <variable name="K_Cai" public_interface="in" units="per_fmol"/>
        <variable name="K_Cao" public_interface="in" units="per_fmol"/>
        <variable name="K_Ki" public_interface="in" units="per_fmol"/>
        <variable name="K_Ko" public_interface="in" units="per_fmol"/>
        <variable name="kappa_Ca1" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_Ca2" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_K1" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_K2" public_interface="in" units="fmol_per_sec"/>
        <variable name="zCa" public_interface="in" units="dimensionless"/>
        <variable name="zK" public_interface="in" units="dimensionless"/>
        <variable name="q_Cai_global" public_interface="in" units="fmol"/>
        <variable name="q_Cao_global" public_interface="in" units="fmol"/>
        <variable name="q_Ki_global" public_interface="in" units="fmol"/>
        <variable name="q_Ko_global" public_interface="in" units="fmol"/>
        <variable initial_value="1e-16" name="q_Cai" public_interface="out" units="fmol"/>
        <variable initial_value="1e-16" name="q_Cao" public_interface="out" units="fmol"/>
        <variable initial_value="1e-16" name="q_Ki" public_interface="out" units="fmol"/>
        <variable initial_value="1e-16" name="q_Ko" public_interface="out" units="fmol"/>
        <variable name="q_E" public_interface="in" units="fC"/>
        <variable name="mu_E" public_interface="out" units="J_per_C"/>
        <variable name="mu_mem_Ca" units="J_per_mol"/>
        <variable name="mu_mem_K" units="J_per_mol"/>
        <variable name="v_CaGHK_i" units="fmol_per_sec"/>
        <variable name="v_CaGHK_o" units="fmol_per_sec"/>
        <variable name="v_KGHK_i" units="fmol_per_sec"/>
        <variable name="v_KGHK_o" units="fmol_per_sec"/>
        <variable name="mu_Cai" public_interface="out" units="J_per_mol"/>
        <variable name="mu_Cao" units="J_per_mol"/>
        <variable name="mu_Ki" units="J_per_mol"/>
        <variable name="mu_Ko" units="J_per_mol"/>
        <variable name="mu_S111" public_interface="in" units="J_per_mol"/>
        <variable name="mu_S121" public_interface="in" units="J_per_mol"/>
        <variable name="Af_Ca1" units="J_per_mol"/>
        <variable name="Ar_Ca1" units="J_per_mol"/>
        <variable name="v_Ca1" public_interface="out" units="fmol_per_sec"/>
        <variable name="v_Ca2" public_interface="out" units="fmol_per_sec"/>
        <variable name="v_K1" public_interface="out" units="fmol_per_sec"/>
        <variable name="v_K2" public_interface="out" units="fmol_per_sec"/>
        <variable name="Af_Ca2" units="J_per_mol"/>
        <variable name="Ar_Ca2" units="J_per_mol"/>
        <variable name="Af_K1" units="J_per_mol"/>
        <variable name="Ar_K1" units="J_per_mol"/>
        <variable name="Af_K2" units="J_per_mol"/>
        <variable name="Ar_K2" units="J_per_mol"/>
        <variable name="v_fCa00" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_fCa01" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_fCa02" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_fCa10" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_fCa11" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_fCa12" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_E_LCC" public_interface="out" units="fA"/>
        <variable name="v_gate" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Ca_fgate" public_interface="out" units="fmol_per_sec"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>RT</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                </apply>
            </apply>
            <!--        v_E_LCC = zCa*F*(v_CaGHK_o+v_KGHK_o);-->
            <apply>
                <eq/>
                <ci>v_E_LCC</ci>
                <apply>
                    <times/>
                    <ci>F</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <ci>v_gate</ci>
                            <apply>
                                <times/>
                                <ci>zCa</ci>
                                <apply>
                                    <plus/>
                                    <ci>v_Ca1</ci>
                                    <ci>v_Ca2</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <ci>zK</ci>
                            <apply>
                                <plus/>
                                <ci>v_K1</ci>
                                <ci>v_K2</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!-- potential energies-->
            <apply>
                <eq/>
                <ci>mu_E</ci>
                <apply>
                    <divide/>
                    <ci>q_E</ci>
                    <ci>C_m</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_Cai</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_Cai</ci>
                            <ci>q_Cai_global</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_Cao</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_Cao</ci>
                            <ci>q_Cao_global</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_Ki</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_Ki</ci>
                            <ci>q_Ki_global</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_Ko</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_Ko</ci>
                            <ci>q_Ko_global</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!-- Constitutive relations - Convention: positive flow from inside to outside-->
            <apply>
                <eq/>
                <ci>mu_mem_Ca</ci>
                <apply>
                    <times/>
                    <ci>zCa</ci>
                    <ci>F</ci>
                    <ci>mu_E</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_mem_K</ci>
                <apply>
                    <times/>
                    <ci>zK</ci>
                    <ci>F</ci>
                    <ci>mu_E</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Af_Ca1</ci>
                <apply>
                    <plus/>
                    <ci>mu_Cai</ci>
                    <ci>mu_mem_Ca</ci>
                    <ci>mu_S111</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Ar_Ca1</ci>
                <apply>
                    <plus/>
                    <ci>mu_Cao</ci>
                    <ci>mu_S111</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Af_Ca2</ci>
                <apply>
                    <plus/>
                    <ci>mu_Cai</ci>
                    <ci>mu_mem_Ca</ci>
                    <ci>mu_S121</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Ar_Ca2</ci>
                <apply>
                    <plus/>
                    <ci>mu_Cao</ci>
                    <ci>mu_S121</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Af_K1</ci>
                <apply>
                    <plus/>
                    <ci>mu_Ki</ci>
                    <ci>mu_mem_K</ci>
                    <ci>mu_S111</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Ar_K1</ci>
                <apply>
                    <plus/>
                    <ci>mu_Ko</ci>
                    <ci>mu_S111</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Af_K2</ci>
                <apply>
                    <plus/>
                    <ci>mu_Ki</ci>
                    <ci>mu_mem_K</ci>
                    <ci>mu_S121</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Ar_K2</ci>
                <apply>
                    <plus/>
                    <ci>mu_Ko</ci>
                    <ci>mu_S121</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_Ca1</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <ci>kappa_Ca1</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Af_Ca1</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Ar_Ca1</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <eq/>
                            <ci>mu_mem_Ca</ci>
                            <cn cellml:units="J_per_mol">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>kappa_Ca1</ci>
                                        <ci>mu_mem_Ca</ci>
                                    </apply>
                                    <ci>RT</ci>
                                </apply>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <ci>mu_mem_Ca</ci>
                                            <ci>RT</ci>
                                        </apply>
                                    </apply>
                                    <cn cellml:units="dimensionless">1</cn>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Af_Ca1</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Ar_Ca1</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>v_Ca2</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <ci>kappa_Ca2</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Af_Ca2</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Ar_Ca2</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <eq/>
                            <ci>mu_mem_Ca</ci>
                            <cn cellml:units="J_per_mol">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>kappa_Ca2</ci>
                                        <ci>mu_mem_Ca</ci>
                                    </apply>
                                    <ci>RT</ci>
                                </apply>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <ci>mu_mem_Ca</ci>
                                            <ci>RT</ci>
                                        </apply>
                                    </apply>
                                    <cn cellml:units="dimensionless">1</cn>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Af_Ca2</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Ar_Ca2</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>v_K1</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <ci>kappa_K1</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Af_K1</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Ar_K1</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <eq/>
                            <ci>mu_mem_K</ci>
                            <cn cellml:units="J_per_mol">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>kappa_K1</ci>
                                        <ci>mu_mem_K</ci>
                                    </apply>
                                    <ci>RT</ci>
                                </apply>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <ci>mu_mem_K</ci>
                                            <ci>RT</ci>
                                        </apply>
                                    </apply>
                                    <cn cellml:units="dimensionless">1</cn>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Af_K1</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Ar_K1</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>v_K2</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <ci>kappa_K2</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Af_K2</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Ar_K2</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <eq/>
                            <ci>mu_mem_K</ci>
                            <cn cellml:units="J_per_mol">0</cn>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>kappa_K2</ci>
                                        <ci>mu_mem_K</ci>
                                    </apply>
                                    <ci>RT</ci>
                                </apply>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <ci>mu_mem_K</ci>
                                            <ci>RT</ci>
                                        </apply>
                                    </apply>
                                    <cn cellml:units="dimensionless">1</cn>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Af_K2</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <ci>Ar_K2</ci>
                                        <ci>RT</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>v_Ca_fgate</ci>
                <apply>
                    <times/>
                    <ci>zCa</ci>
                    <apply>
                        <plus/>
                        <ci>v_fCa10</ci>
                        <ci>v_fCa11</ci>
                        <ci>v_fCa02</ci>
                        <ci>v_fCa01</ci>
                        <ci>v_fCa00</ci>
                        <ci>v_fCa12</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_CaGHK_i</ci>
                <apply>
                    <minus/>
                    <apply>
                        <minus/>
                        <ci>v_Ca_fgate</ci>
                        <ci>v_Ca1</ci>
                    </apply>
                    <ci>v_Ca2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_CaGHK_o</ci>
                <apply>
                    <plus/>
                    <ci>v_Ca2</ci>
                    <ci>v_Ca1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_KGHK_i</ci>
                <apply>
                    <minus/>
                    <apply>
                        <minus/>
                        <ci>v_K1</ci>
                    </apply>
                    <ci>v_K2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_KGHK_o</ci>
                <apply>
                    <plus/>
                    <ci>v_K2</ci>
                    <ci>v_K1</ci>
                </apply>
            </apply>
        </math>
        <variable name="I_pulse" public_interface="in" units="fA"/>
        <variable name="v_pulse" units="fmol_per_sec"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>v_pulse</ci>
                <apply>
                    <divide/>
                    <ci>I_pulse</ci>
                    <apply>
                        <times/>
                        <ci>zCa</ci>
                        <ci>F</ci>
                    </apply>
                </apply>
            </apply>
            <!--        ode(q_Cai, time) = v_CaGHK_i+v_pulse;
        ode(q_Cao, time) = v_CaGHK_o-v_pulse;-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Cai</ci>
                </apply>
                <ci>v_CaGHK_i</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Cao</ci>
                </apply>
                <ci>v_CaGHK_o</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ki</ci>
                </apply>
                <ci>v_KGHK_i</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ko</ci>
                </apply>
                <ci>v_KGHK_o</ci>
            </apply>
        </math>
    </component>
    <component name="LCC_gate">
        <variable name="time" public_interface="in" units="second"/>
        <variable name="Vmem" public_interface="in" units="volt"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable name="RT" public_interface="in" units="J_per_mol"/>
        <!-- gating parameters-->
        <variable name="kappa_d000" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_d010" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_d020" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_d001" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_d011" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_d021" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f1_000" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f1_100" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f1_001" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f1_101" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f2_000" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f2_100" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f2_001" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f2_101" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f3_010" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f3_110" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f3_011" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_f3_111" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_fCa000" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_fCa100" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_fCa010" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_fCa110" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_fCa020" public_interface="in" units="fmol_per_sec"/>
        <variable name="kappa_fCa120" public_interface="in" units="fmol_per_sec"/>
        <variable name="K_000" public_interface="in" units="per_fmol"/>
        <variable name="K_010" public_interface="in" units="per_fmol"/>
        <variable name="K_020" public_interface="in" units="per_fmol"/>
        <variable name="K_100" public_interface="in" units="per_fmol"/>
        <variable name="K_110" public_interface="in" units="per_fmol"/>
        <variable name="K_120" public_interface="in" units="per_fmol"/>
        <variable name="K_001" public_interface="in" units="per_fmol"/>
        <variable name="K_011" public_interface="in" units="per_fmol"/>
        <variable name="K_021" public_interface="in" units="per_fmol"/>
        <variable name="K_101" public_interface="in" units="per_fmol"/>
        <variable name="K_111" public_interface="in" units="per_fmol"/>
        <variable name="K_121" public_interface="in" units="per_fmol"/>
        <!-- BG mol quantities  -->
        <variable initial_value="3.6989e-07" name="q_S000" units="fmol"/>
        <variable initial_value="7.3239e-05" name="q_S010" units="fmol"/>
        <variable initial_value="3.6989e-07" name="q_S020" units="fmol"/>
        <variable initial_value="3.7363e-09" name="q_S100" units="fmol"/>
        <variable initial_value="7.3979e-07" name="q_S110" units="fmol"/>
        <variable initial_value="3.7363e-09" name="q_S120" units="fmol"/>
        <variable initial_value="4.1099e-08" name="q_S001" units="fmol"/>
        <variable initial_value="8.1377e-06" name="q_S011" units="fmol"/>
        <variable initial_value="4.1099e-08" name="q_S021" units="fmol"/>
        <variable initial_value="4.1514e-10" name="q_S101" units="fmol"/>
        <variable initial_value="8.2199e-08" name="q_S111" units="fmol"/>
        <variable initial_value="4.1514e-10" name="q_S121" units="fmol"/>
        <variable name="q_Ca_fgate" public_interface="out" units="fmol"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>q_Ca_fgate</ci>
                <apply>
                    <times/>
                    <ci>z_rCa</ci>
                    <apply>
                        <plus/>
                        <ci>q_S000</ci>
                        <ci>q_S010</ci>
                        <ci>q_S020</ci>
                        <ci>q_S100</ci>
                        <ci>q_S110</ci>
                        <ci>q_S120</ci>
                    </apply>
                </apply>
            </apply>
        </math>
        <!-- gating variable potential energy-->
        <variable name="mu_Cai" public_interface="in" units="J_per_mol"/>
        <variable name="mu_S000" units="J_per_mol"/>
        <variable name="mu_S010" units="J_per_mol"/>
        <variable name="mu_S020" units="J_per_mol"/>
        <variable name="mu_S100" units="J_per_mol"/>
        <variable name="mu_S110" units="J_per_mol"/>
        <variable name="mu_S120" units="J_per_mol"/>
        <variable name="mu_S001" units="J_per_mol"/>
        <variable name="mu_S011" units="J_per_mol"/>
        <variable name="mu_S021" units="J_per_mol"/>
        <variable name="mu_S101" units="J_per_mol"/>
        <variable name="mu_S111" public_interface="out" units="J_per_mol"/>
        <variable name="mu_S121" public_interface="out" units="J_per_mol"/>
        <variable name="v_fCa00" public_interface="out" units="fmol_per_sec"/>
        <variable name="v_fCa01" public_interface="out" units="fmol_per_sec"/>
        <variable name="v_fCa02" public_interface="out" units="fmol_per_sec"/>
        <variable name="v_fCa10" public_interface="out" units="fmol_per_sec"/>
        <variable name="v_fCa11" public_interface="out" units="fmol_per_sec"/>
        <variable name="v_fCa12" public_interface="out" units="fmol_per_sec"/>
        <variable name="v_f1_000" units="fmol_per_sec"/>
        <variable name="v_f2_000" units="fmol_per_sec"/>
        <variable name="v_f3_010" units="fmol_per_sec"/>
        <variable name="v_f1_100" units="fmol_per_sec"/>
        <variable name="v_f2_100" units="fmol_per_sec"/>
        <variable name="v_f3_110" units="fmol_per_sec"/>
        <variable name="v_d000" units="fmol_per_sec"/>
        <variable name="v_d010" units="fmol_per_sec"/>
        <variable name="v_d020" units="fmol_per_sec"/>
        <variable name="v_d001" units="fmol_per_sec"/>
        <variable name="v_d011" units="fmol_per_sec"/>
        <variable name="v_d021" units="fmol_per_sec"/>
        <variable name="v_f1_001" units="fmol_per_sec"/>
        <variable name="v_f2_001" units="fmol_per_sec"/>
        <variable name="v_f3_011" units="fmol_per_sec"/>
        <variable name="v_f1_101" units="fmol_per_sec"/>
        <variable name="v_f2_101" units="fmol_per_sec"/>
        <variable name="v_f3_111" units="fmol_per_sec"/>
        <variable name="v_S000" units="fmol_per_sec"/>
        <variable name="v_S010" units="fmol_per_sec"/>
        <variable name="v_S020" units="fmol_per_sec"/>
        <variable name="v_S100" units="fmol_per_sec"/>
        <variable name="v_S110" units="fmol_per_sec"/>
        <variable name="v_S120" units="fmol_per_sec"/>
        <variable name="v_S001" units="fmol_per_sec"/>
        <variable name="v_S011" units="fmol_per_sec"/>
        <variable name="v_S021" units="fmol_per_sec"/>
        <variable name="v_S101" units="fmol_per_sec"/>
        <variable name="v_S111" units="fmol_per_sec"/>
        <variable name="v_S121" units="fmol_per_sec"/>
        <variable initial_value="2" name="z_rCa" units="dimensionless"/>
        <variable initial_value="2.1404" name="z_fd" units="dimensionless"/>
        <variable initial_value="-1.1495" name="z_ff1" units="dimensionless"/>
        <variable initial_value="0.72162" name="z_ff2" units="dimensionless"/>
        <variable initial_value="4.2933" name="z_ff3" units="dimensionless"/>
        <variable initial_value="-2.1404" name="z_rd" units="dimensionless"/>
        <variable initial_value="1.8993" name="z_rf1" units="dimensionless"/>
        <variable initial_value="-0.52288" name="z_rf2" units="dimensionless"/>
        <variable initial_value="0" name="z_rf3" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>mu_S000</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_000</ci>
                            <ci>q_S000</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S001</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_001</ci>
                            <ci>q_S001</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S010</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_010</ci>
                            <ci>q_S010</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S011</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_011</ci>
                            <ci>q_S011</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S020</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_020</ci>
                            <ci>q_S020</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S021</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_021</ci>
                            <ci>q_S021</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S100</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_100</ci>
                            <ci>q_S100</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S101</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_101</ci>
                            <ci>q_S101</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S110</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_110</ci>
                            <ci>q_S110</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S111</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_111</ci>
                            <ci>q_S111</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S120</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_120</ci>
                            <ci>q_S120</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_S121</ci>
                <apply>
                    <times/>
                    <ci>RT</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_121</ci>
                            <ci>q_S121</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!-- individual gate transition rates-->
            <apply>
                <eq/>
                <ci>v_fCa00</ci>
                <apply>
                    <times/>
                    <ci>kappa_fCa000</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_S000</ci>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S001</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rCa</ci>
                                        <ci>mu_Cai</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_fCa01</ci>
                <apply>
                    <times/>
                    <ci>kappa_fCa010</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_S010</ci>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S011</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rCa</ci>
                                        <ci>mu_Cai</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_fCa02</ci>
                <apply>
                    <times/>
                    <ci>kappa_fCa020</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_S020</ci>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S021</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rCa</ci>
                                        <ci>mu_Cai</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_fCa10</ci>
                <apply>
                    <times/>
                    <ci>kappa_fCa100</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_S100</ci>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S101</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rCa</ci>
                                        <ci>mu_Cai</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_fCa11</ci>
                <apply>
                    <times/>
                    <ci>kappa_fCa110</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_S110</ci>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S111</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rCa</ci>
                                        <ci>mu_Cai</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_fCa12</ci>
                <apply>
                    <times/>
                    <ci>kappa_fCa120</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_S120</ci>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S121</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rCa</ci>
                                        <ci>mu_Cai</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_d000</ci>
                <apply>
                    <times/>
                    <ci>kappa_d000</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S000</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_fd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S100</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_d010</ci>
                <apply>
                    <times/>
                    <ci>kappa_d010</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S010</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_fd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S110</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_d020</ci>
                <apply>
                    <times/>
                    <ci>kappa_d020</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S020</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_fd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S120</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_d001</ci>
                <apply>
                    <times/>
                    <ci>kappa_d001</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S001</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_fd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S101</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_d011</ci>
                <apply>
                    <times/>
                    <ci>kappa_d011</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S011</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_fd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S111</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_d021</ci>
                <apply>
                    <times/>
                    <ci>kappa_d021</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S021</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_fd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S121</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rd</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f1_000</ci>
                <apply>
                    <times/>
                    <ci>kappa_f1_000</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S000</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff1</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S010</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf1</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f2_000</ci>
                <apply>
                    <times/>
                    <ci>kappa_f2_000</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S000</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff2</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S020</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf2</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f3_010</ci>
                <apply>
                    <times/>
                    <ci>kappa_f3_010</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S010</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff3</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S020</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf3</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f1_100</ci>
                <apply>
                    <times/>
                    <ci>kappa_f1_100</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S100</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff1</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S110</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf1</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f2_100</ci>
                <apply>
                    <times/>
                    <ci>kappa_f2_100</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S100</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff2</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S120</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf2</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f3_110</ci>
                <apply>
                    <times/>
                    <ci>kappa_f3_110</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S110</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff3</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S120</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf3</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f1_001</ci>
                <apply>
                    <times/>
                    <ci>kappa_f1_001</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S001</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff1</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S011</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf1</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f2_001</ci>
                <apply>
                    <times/>
                    <ci>kappa_f2_001</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S001</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff2</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S021</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf2</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f3_011</ci>
                <apply>
                    <times/>
                    <ci>kappa_f3_011</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S011</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff3</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S021</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf3</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f1_101</ci>
                <apply>
                    <times/>
                    <ci>kappa_f1_101</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S101</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff1</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S111</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf1</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f2_101</ci>
                <apply>
                    <times/>
                    <ci>kappa_f2_101</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S101</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff2</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S121</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf2</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_f3_111</ci>
                <apply>
                    <times/>
                    <ci>kappa_f3_111</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S111</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_ff3</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_S121</ci>
                                    <apply>
                                        <times/>
                                        <ci>z_rf3</ci>
                                        <ci>F</ci>
                                        <ci>Vmem</ci>
                                    </apply>
                                </apply>
                                <ci>RT</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!-- whole gate transition rates-->
            <apply>
                <eq/>
                <ci>v_S000</ci>
                <apply>
                    <minus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <apply>
                                <minus/>
                                <ci>v_d000</ci>
                            </apply>
                            <ci>v_f1_000</ci>
                        </apply>
                        <ci>v_f2_000</ci>
                    </apply>
                    <ci>v_fCa00</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S001</ci>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <apply>
                                <minus/>
                                <ci>v_d001</ci>
                            </apply>
                            <ci>v_f1_001</ci>
                        </apply>
                        <ci>v_f2_001</ci>
                    </apply>
                    <ci>v_fCa00</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S010</ci>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <apply>
                            <minus/>
                            <apply>
                                <minus/>
                                <ci>v_d010</ci>
                            </apply>
                            <ci>v_f3_010</ci>
                        </apply>
                        <ci>v_f1_100</ci>
                    </apply>
                    <ci>v_fCa01</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S011</ci>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <ci>v_d011</ci>
                        </apply>
                        <ci>v_f3_011</ci>
                    </apply>
                    <ci>v_f1_101</ci>
                    <ci>v_fCa01</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S020</ci>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <apply>
                            <minus/>
                            <ci>v_d020</ci>
                        </apply>
                        <ci>v_f3_010</ci>
                        <ci>v_f2_000</ci>
                    </apply>
                    <ci>v_fCa02</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S021</ci>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>v_d021</ci>
                    </apply>
                    <ci>v_f3_011</ci>
                    <ci>v_f2_001</ci>
                    <ci>v_fCa02</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S100</ci>
                <apply>
                    <minus/>
                    <apply>
                        <minus/>
                        <apply>
                            <plus/>
                            <apply>
                                <minus/>
                                <ci>v_f2_100</ci>
                            </apply>
                            <ci>v_d000</ci>
                        </apply>
                        <ci>v_f1_100</ci>
                    </apply>
                    <ci>v_fCa10</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S101</ci>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <plus/>
                            <apply>
                                <minus/>
                                <ci>v_f2_101</ci>
                            </apply>
                            <ci>v_d001</ci>
                        </apply>
                        <ci>v_f1_101</ci>
                    </apply>
                    <ci>v_fCa10</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S110</ci>
                <apply>
                    <minus/>
                    <apply>
                        <minus/>
                        <apply>
                            <plus/>
                            <ci>v_f1_100</ci>
                            <ci>v_d010</ci>
                        </apply>
                        <ci>v_fCa11</ci>
                    </apply>
                    <ci>v_f3_110</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S111</ci>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <ci>v_f1_101</ci>
                        <ci>v_d011</ci>
                        <ci>v_fCa11</ci>
                    </apply>
                    <ci>v_f3_111</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S120</ci>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <plus/>
                            <ci>v_f2_100</ci>
                            <ci>v_d020</ci>
                        </apply>
                        <ci>v_fCa12</ci>
                    </apply>
                    <ci>v_f3_110</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v_S121</ci>
                <apply>
                    <plus/>
                    <ci>v_f2_101</ci>
                    <ci>v_d021</ci>
                    <ci>v_fCa12</ci>
                    <ci>v_f3_111</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S000</ci>
                </apply>
                <ci>v_S000</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S001</ci>
                </apply>
                <ci>v_S001</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S010</ci>
                </apply>
                <ci>v_S010</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S011</ci>
                </apply>
                <ci>v_S011</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S020</ci>
                </apply>
                <ci>v_S020</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S021</ci>
                </apply>
                <ci>v_S021</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S100</ci>
                </apply>
                <ci>v_S100</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S101</ci>
                </apply>
                <ci>v_S101</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S110</ci>
                </apply>
                <ci>v_S110</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S111</ci>
                </apply>
                <ci>v_S111</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S120</ci>
                </apply>
                <ci>v_S120</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S121</ci>
                </apply>
                <ci>v_S121</ci>
            </apply>
        </math>
        <variable name="v_gate" public_interface="out" units="fmol_per_sec"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>v_gate</ci>
                <apply>
                    <plus/>
                    <apply>
                        <times/>
                        <ci>v_f2_001</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf2</ci>
                            <ci>z_ff2</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f3_011</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf3</ci>
                            <ci>z_ff3</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f1_101</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf1</ci>
                            <ci>z_ff1</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f2_101</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf2</ci>
                            <ci>z_ff2</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f3_111</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf3</ci>
                            <ci>z_ff3</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f1_001</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf1</ci>
                            <ci>z_ff1</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_d021</ci>
                        <apply>
                            <minus/>
                            <ci>z_rd</ci>
                            <ci>z_fd</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_d011</ci>
                        <apply>
                            <minus/>
                            <ci>z_rd</ci>
                            <ci>z_fd</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_d001</ci>
                        <apply>
                            <minus/>
                            <ci>z_rd</ci>
                            <ci>z_fd</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f3_110</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf3</ci>
                            <ci>z_ff3</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f2_100</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf2</ci>
                            <ci>z_ff2</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f1_100</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf1</ci>
                            <ci>z_ff1</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f3_010</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf3</ci>
                            <ci>z_ff3</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f2_000</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf2</ci>
                            <ci>z_ff2</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_f1_000</ci>
                        <apply>
                            <minus/>
                            <ci>z_rf1</ci>
                            <ci>z_ff1</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_d020</ci>
                        <apply>
                            <minus/>
                            <ci>z_rd</ci>
                            <ci>z_fd</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_d010</ci>
                        <apply>
                            <minus/>
                            <ci>z_rd</ci>
                            <ci>z_fd</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>v_d000</ci>
                        <apply>
                            <minus/>
                            <ci>z_rd</ci>
                            <ci>z_fd</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="environment" component_2="LCC"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="q_Cai_mLCC" variable_2="q_Cai"/>
        <map_variables variable_1="q_Cai" variable_2="q_Cai_global"/>
        <map_variables variable_1="q_Cao_mLCC" variable_2="q_Cao"/>
        <map_variables variable_1="q_Cao" variable_2="q_Cao_global"/>
        <map_variables variable_1="q_Ki_mLCC" variable_2="q_Ki"/>
        <map_variables variable_1="q_Ki" variable_2="q_Ki_global"/>
        <map_variables variable_1="q_Ko_mLCC" variable_2="q_Ko"/>
        <map_variables variable_1="q_Ko" variable_2="q_Ko_global"/>
        <map_variables variable_1="q_membrane" variable_2="q_E"/>
        <map_variables variable_1="C_m" variable_2="C_m"/>
        <map_variables variable_1="v_CaGHK_i" variable_2="v_CaGHK_i"/>
        <map_variables variable_1="v_CaGHK_o" variable_2="v_CaGHK_o"/>
        <map_variables variable_1="v_KGHK_i" variable_2="v_KGHK_i"/>
        <map_variables variable_1="v_KGHK_o" variable_2="v_KGHK_o"/>
        <map_variables variable_1="zCa" variable_2="zCa"/>
        <map_variables variable_1="zK" variable_2="zK"/>
        <map_variables variable_1="v_E_LCC" variable_2="v_E_LCC"/>
        <map_variables variable_1="I_pulse" variable_2="I_pulse"/>
        <!--        vars v_LCC_Ca1 and v_Ca1;
        vars v_LCC_Ca2 and v_Ca2;
        vars v_LCC_K1 and v_K1;
        vars v_LCC_K2 and v_K2;-->
    </connection>
    <connection>
        <map_components component_1="LCC" component_2="LCC_gate"/>
        <map_variables variable_1="RT" variable_2="RT"/>
        <map_variables variable_1="mu_E" variable_2="Vmem"/>
        <map_variables variable_1="mu_Cai" variable_2="mu_Cai"/>
        <map_variables variable_1="mu_S111" variable_2="mu_S111"/>
        <map_variables variable_1="mu_S121" variable_2="mu_S121"/>
        <map_variables variable_1="v_fCa00" variable_2="v_fCa00"/>
        <map_variables variable_1="v_fCa01" variable_2="v_fCa01"/>
        <map_variables variable_1="v_fCa02" variable_2="v_fCa02"/>
        <map_variables variable_1="v_fCa10" variable_2="v_fCa10"/>
        <map_variables variable_1="v_fCa11" variable_2="v_fCa11"/>
        <map_variables variable_1="v_fCa12" variable_2="v_fCa12"/>
        <map_variables variable_1="v_gate" variable_2="v_gate"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="LCC_gate"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="LCC_gate"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="LCC_parameters" component_2="LCC"/>
        <map_variables variable_1="K_Cai" variable_2="K_Cai"/>
        <map_variables variable_1="K_Cao" variable_2="K_Cao"/>
        <map_variables variable_1="K_Ki" variable_2="K_Ki"/>
        <map_variables variable_1="K_Ko" variable_2="K_Ko"/>
        <map_variables variable_1="kappa_Ca1" variable_2="kappa_Ca1"/>
        <map_variables variable_1="kappa_Ca2" variable_2="kappa_Ca2"/>
        <map_variables variable_1="kappa_K1" variable_2="kappa_K1"/>
        <map_variables variable_1="kappa_K2" variable_2="kappa_K2"/>
    </connection>
    <connection>
        <map_components component_1="LCC_parameters" component_2="LCC_gate"/>
        <map_variables variable_1="kappa_d000" variable_2="kappa_d000"/>
        <map_variables variable_1="kappa_d010" variable_2="kappa_d010"/>
        <map_variables variable_1="kappa_d020" variable_2="kappa_d020"/>
        <map_variables variable_1="kappa_d001" variable_2="kappa_d001"/>
        <map_variables variable_1="kappa_d011" variable_2="kappa_d011"/>
        <map_variables variable_1="kappa_d021" variable_2="kappa_d021"/>
        <map_variables variable_1="kappa_f1_000" variable_2="kappa_f1_000"/>
        <map_variables variable_1="kappa_f1_100" variable_2="kappa_f1_100"/>
        <map_variables variable_1="kappa_f1_001" variable_2="kappa_f1_001"/>
        <map_variables variable_1="kappa_f1_101" variable_2="kappa_f1_101"/>
        <map_variables variable_1="kappa_f2_000" variable_2="kappa_f2_000"/>
        <map_variables variable_1="kappa_f2_100" variable_2="kappa_f2_100"/>
        <map_variables variable_1="kappa_f2_001" variable_2="kappa_f2_001"/>
        <map_variables variable_1="kappa_f2_101" variable_2="kappa_f2_101"/>
        <map_variables variable_1="kappa_f3_010" variable_2="kappa_f3_010"/>
        <map_variables variable_1="kappa_f3_110" variable_2="kappa_f3_110"/>
        <map_variables variable_1="kappa_f3_011" variable_2="kappa_f3_011"/>
        <map_variables variable_1="kappa_f3_111" variable_2="kappa_f3_111"/>
        <map_variables variable_1="kappa_fCa000" variable_2="kappa_fCa000"/>
        <map_variables variable_1="kappa_fCa100" variable_2="kappa_fCa100"/>
        <map_variables variable_1="kappa_fCa010" variable_2="kappa_fCa010"/>
        <map_variables variable_1="kappa_fCa110" variable_2="kappa_fCa110"/>
        <map_variables variable_1="kappa_fCa020" variable_2="kappa_fCa020"/>
        <map_variables variable_1="kappa_fCa120" variable_2="kappa_fCa120"/>
        <map_variables variable_1="K_000_LCC" variable_2="K_000"/>
        <map_variables variable_1="K_010_LCC" variable_2="K_010"/>
        <map_variables variable_1="K_020_LCC" variable_2="K_020"/>
        <map_variables variable_1="K_100_LCC" variable_2="K_100"/>
        <map_variables variable_1="K_110_LCC" variable_2="K_110"/>
        <map_variables variable_1="K_120_LCC" variable_2="K_120"/>
        <map_variables variable_1="K_001_LCC" variable_2="K_001"/>
        <map_variables variable_1="K_011_LCC" variable_2="K_011"/>
        <map_variables variable_1="K_021_LCC" variable_2="K_021"/>
        <map_variables variable_1="K_101_LCC" variable_2="K_101"/>
        <map_variables variable_1="K_111_LCC" variable_2="K_111"/>
        <map_variables variable_1="K_121_LCC" variable_2="K_121"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="environment"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="LCC"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="R" variable_2="R"/>
    </connection>
</model>