Location: FCU_ExcitationContraction_Coupling @ 8bb03dd23f61 / FCU_EC_coupling.cellml

Author:
Shelley Fong <s.fong@auckland.ac.nz>
Date:
2022-02-18 14:17:34+13:00
Desc:
Updating way cellml is written
Permanent Source URI:
http://models.cellml.org/workspace/7a8/rawfile/8bb03dd23f61cba2ac9583321a845c02ff6f4155/FCU_EC_coupling.cellml

<?xml version='1.0' encoding='UTF-8'?>
<model name="FCU_EC_coupling" 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="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"/>
        <units name="per_sec" units_ref="per_sec"/>
        <units name="J_per_K_per_mol" units_ref="J_per_K_per_mol"/>
        <units name="fmol_per_L" units_ref="fmol_per_L"/>
        <units name="fmol_per_L_per_sec" units_ref="fmol_per_L_per_sec"/>
        <units name="per_sec_per_fmol_per_L" units_ref="per_sec_per_fmol_per_L"/>
        <units name="uM" units_ref="uM"/>
        <units name="mM_per_sec" units_ref="mM_per_sec"/>
        <units name="uM_per_sec" units_ref="uM_per_sec"/>
        <units name="pL" units_ref="pL"/>
        <units name="m_to_u" units_ref="m_to_u"/>
    </import>
    <import xlink:href="units_and_constants/constants_BG.cellml">
        <component component_ref="constants" name="constants"/>
    </import>
    <import xlink:href="units_and_constants/ion_valency.cellml">
        <component component_ref="ion_valences" name="ion_valences"/>
    </import>
    <import xlink:href="crossbridge_TRPN/BG_crossbridge_TRPN.cellml">
        <component component_ref="crossbridge_TRPN" name="crossbridge_TRPN"/>
    </import>
    <import xlink:href="LCC/BG_LCC.cellml">
        <component component_ref="LCC" name="LCC"/>
        <component component_ref="LCC_gate" name="LCC_gate"/>
    </import>
    <import xlink:href="PLB/BG_PLB.cellml">
        <component component_ref="PLB" name="PLB"/>
    </import>
    <import xlink:href="RyR/BG_RyR.cellml">
        <component component_ref="RyR" name="RyR"/>
        <component component_ref="RyR_gate" name="RyR_gate"/>
    </import>
    <import xlink:href="SERCA/BG_SERCA.cellml">
        <component component_ref="SERCA" name="SERCA"/>
    </import>
    <component name="BG_parameters">
        <variable initial_value="114.037" name="kappa_R_TRPNCa" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6.18867" name="kappa_R_BU" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.298505" name="kappa_R_UW" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00702365" name="kappa_R_WS" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="126426" name="kappa_R_SU" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="8.20611" name="kappa_Ca1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="15.3452" name="kappa_Ca2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.0925573" name="kappa_K1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.17308" name="kappa_K2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="956.69" name="kappa_d000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="23.603" name="kappa_d010" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="44.137" name="kappa_d020" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.174044" name="kappa_d001" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00429394" name="kappa_d011" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00802955" name="kappa_d021" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.76246" name="kappa_f1_000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="8.75195" name="kappa_f1_100" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.000320633" name="kappa_f1_001" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00159218" name="kappa_f1_101" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="3.22785" name="kappa_f2_000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="16.0287" name="kappa_f2_100" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.000587221" name="kappa_f2_001" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00291599" name="kappa_f2_101" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="9067.48" name="kappa_f3_010" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="45026.8" name="kappa_f3_110" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.64958" name="kappa_f3_011" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="8.19142" name="kappa_f3_111" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="59592.5" name="kappa_fCa000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="295921" name="kappa_fCa100" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1470.24" name="kappa_fCa010" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="7300.85" name="kappa_fCa110" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="2749.31" name="kappa_fCa020" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="13652.4" name="kappa_fCa120" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="45.5263" name="kappa_PLBph1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6.55904" name="kappa_PLBph2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.386674" name="kappa_PLBd1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.21269" name="kappa_PLBd2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="431.435" name="kappa_Inh" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="3.53995e+06" name="kappa_RyR" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="8.83262" name="kappa_OC" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.12618" name="kappa_CCI" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="883.262" name="kappa_CII" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="73.6052" name="kappa_IO" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="2.27229e-05" name="kappa_SERCA_R1_2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="106884" name="kappa_SERCA_R2_4" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="2.78115e+06" name="kappa_SERCA_R2_2a" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="106884" name="kappa_SERCA_R4_5" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="3440.8" name="kappa_SERCA_R5_6" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="3.04364e+06" name="kappa_SERCA_R6_8" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.03099e+07" name="kappa_SERCA_R8_9" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.03099e+07" name="kappa_SERCA_R9_10" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00568924" name="kappa_SERCA_R10_1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="14.0366" name="K_TRPN" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0527926" name="K_Cai" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0509829" name="K_Ca_TRPN" public_interface="out" units="per_fmol"/>
        <variable initial_value="6.75199" name="K_B" public_interface="out" units="per_fmol"/>
        <variable initial_value="2.532" name="K_U" public_interface="out" units="per_fmol"/>
        <variable initial_value="16.5554" name="K_W" public_interface="out" units="per_fmol"/>
        <variable initial_value="4.13884e-06" name="K_S" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0527926" name="K_Cao" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.00367934" name="K_Ki" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.00367934" name="K_Ko" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0147906" name="K_000_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.599501" name="K_010_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.320594" name="K_020_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.00297853" name="K_100_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.120727" name="K_110_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.064561" name="K_120_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="81.3015" name="K_001_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="3295.35" name="K_011_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="1762.25" name="K_021_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="16.3725" name="K_101_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="663.616" name="K_111_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="354.881" name="K_121_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.00235741" name="K_PLB" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.395191" name="K_PKACI" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.638527" name="K_PLB_PKACI" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.361989" name="K_PP1" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.203757" name="K_PLBp_PP1" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.326014" name="K_PLBp" public_interface="out" units="per_fmol"/>
        <variable initial_value="5.41093" name="K_Ip" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0673793" name="K_Ip_PP1" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0527926" name="K_Ca_SR" public_interface="out" units="per_fmol"/>
        <variable initial_value="115.191" name="K_C_RyR" public_interface="out" units="per_fmol"/>
        <variable initial_value="1.15191" name="K_CI_RyR" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.00197471" name="K_I_RyR" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.197471" name="K_O_RyR" public_interface="out" units="per_fmol"/>
        <variable initial_value="901.774" name="K_P1_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="2579.2" name="K_P2_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="270.718" name="K_P2a_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="7044.18" name="K_P4_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0214593" name="K_P5_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="1241.21" name="K_P6_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="75.0041" name="K_P8_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="73.0274" name="K_P9_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="104.747" name="K_P10_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="279.929" name="K_H" public_interface="out" units="per_fmol"/>
        <variable initial_value="1068.12" name="K_MgATP" public_interface="out" units="per_fmol"/>
        <variable initial_value="1.32969e-05" name="K_MgADP" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0245424" name="K_P" public_interface="out" units="per_fmol"/>
        <variable initial_value="2" name="z_rCa" public_interface="out" units="dimensionless"/>
        <variable initial_value="2.1404" name="z_fd" public_interface="out" units="dimensionless"/>
        <variable initial_value="-1.1495" name="z_ff1" public_interface="out" units="dimensionless"/>
        <variable initial_value="0.72162" name="z_ff2" public_interface="out" units="dimensionless"/>
        <variable initial_value="4.2933" name="z_ff3" public_interface="out" units="dimensionless"/>
        <variable initial_value="-2.1404" name="z_rd" public_interface="out" units="dimensionless"/>
        <variable initial_value="1.8993" name="z_rf1" public_interface="out" units="dimensionless"/>
        <variable initial_value="-0.52288" name="z_rf2" public_interface="out" units="dimensionless"/>
        <variable initial_value="0" name="z_rf3" public_interface="out" units="dimensionless"/>
    </component>
    <component name="environment">
        <variable name="time" public_interface="out" units="second"/>
        <variable initial_value="34.4" name="vol_myo" public_interface="out" units="pL"/>
        <variable initial_value="500" name="freq" units="dimensionless"/>
        <variable initial_value="1.381e5" name="C_m" public_interface="out" units="fF"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable initial_value="-13039" name="q_membrane" public_interface="out" units="fC"/>
        <variable name="V_m" public_interface="out" units="volt"/>
        <variable initial_value="2.57" name="q_TRPN" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_Ca_TRPN" public_interface="out" units="fmol"/>
        <variable initial_value="6.82e-3" name="q_Cai" public_interface="out" units="fmol"/>
        <variable initial_value="1" name="q_B" public_interface="out" units="fmol"/>
        <variable initial_value="0" name="q_U" public_interface="out" units="fmol"/>
        <variable initial_value="0" name="q_W" public_interface="out" units="fmol"/>
        <variable initial_value="0" name="q_S" public_interface="out" units="fmol"/>
        <variable initial_value="6.84" name="q_Cao" public_interface="out" units="fmol"/>
        <variable initial_value="5.51E+03" name="q_Ki" public_interface="out" units="fmol"/>
        <variable initial_value="2.05E+02" name="q_Ko" public_interface="out" units="fmol"/>
        <variable initial_value="4.028E+00" name="q_PLB" public_interface="out" units="fmol"/>
        <variable initial_value="2.234E-03" name="q_PKACI" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_PLB_PKACI" public_interface="out" units="fmol"/>
        <variable initial_value="3.382E-02" name="q_PP1" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_PLBp_PP1" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_PLBp" public_interface="out" units="fmol"/>
        <variable initial_value="1.999E-03" name="q_Ip" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_Ip_PP1" public_interface="out" units="fmol"/>
        <variable initial_value="0.641" name="q_Ca_SR" public_interface="out" units="fmol"/>
        <variable initial_value="0.004028" name="q_H" public_interface="out" units="fmol"/>
        <variable initial_value="570" name="q_P" public_interface="out" units="fmol"/>
        <variable initial_value="1.3794" name="q_MgADP" public_interface="out" units="fmol"/>
        <variable initial_value="3.8" name="q_MgATP" public_interface="out" units="fmol"/>
        <!-- crossbridge_TRPN imports-->
        <variable name="v_R_TRPNCa" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_R_BU" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_R_UW" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_R_WS" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_R_SU" public_interface="in" units="fmol_per_sec"/>
        <variable name="Tension" public_interface="in" units="kPa"/>
        <variable name="SL" public_interface="in" units="metre"/>
        <!-- LCC imports-->
        <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"/>
        <!-- PLB imports-->
        <variable name="v_PLBph1" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_PLBph2" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_PLBd1" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_PLBd2" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Inh" public_interface="in" units="fmol_per_sec"/>
        <!-- RyR imports-->
        <variable name="v_RyR" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_gate_Cai" public_interface="in" units="fmol_per_sec"/>
        <!-- SERCA imports-->
        <variable name="v_Cai_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Ca_SR_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_H_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_MgADP_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_MgATP_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_P_SERCA" public_interface="in" units="fmol_per_sec"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Cai</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>v_R_TRPNCa</ci>
                    </apply>
                    <ci>v_CaGHK_i</ci>
                    <ci>v_RyR</ci>
                    <ci>v_gate_Cai</ci>
                    <ci>v_Cai_SERCA</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_TRPN</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_R_TRPNCa</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ca_TRPN</ci>
                </apply>
                <ci>v_R_TRPNCa</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_B</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_R_BU</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_U</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>v_R_BU</ci>
                        <ci>v_R_UW</ci>
                    </apply>
                    <ci>v_R_SU</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_W</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_R_UW</ci>
                    <ci>v_R_WS</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_R_WS</ci>
                    <ci>v_R_SU</ci>
                </apply>
            </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>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PLB</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBd2</ci>
                    <ci>v_PLBph1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PKACI</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBph2</ci>
                    <ci>v_PLBph1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PLB_PKACI</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBph1</ci>
                    <ci>v_PLBph2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PP1</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <minus/>
                        <ci>v_PLBd2</ci>
                        <ci>v_PLBd1</ci>
                    </apply>
                    <ci>v_Inh</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PLBp_PP1</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBd1</ci>
                    <ci>v_PLBd2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PLBp</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBph2</ci>
                    <ci>v_PLBd1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ip</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_Inh</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ip_PP1</ci>
                </apply>
                <ci>v_Inh</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ca_SR</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>v_RyR</ci>
                    </apply>
                    <ci>v_Ca_SR_SERCA</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_H</ci>
                </apply>
                <ci>v_H_SERCA</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_MgADP</ci>
                </apply>
                <ci>v_MgADP_SERCA</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_MgATP</ci>
                </apply>
                <ci>v_MgATP_SERCA</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_P</ci>
                </apply>
                <ci>v_P_SERCA</ci>
            </apply>
        </math>
        <!-- membrane potential change-->
        <variable name="I_pulse" public_interface="out" units="fA"/>
        <variable initial_value="4e-2" name="pulse_start" units="second"/>
        <variable initial_value="4.1e-2" name="pulse_end" units="second"/>
        <variable initial_value="1e8" name="pulseMag" units="fA"/>
        <variable initial_value="0" name="pulseHolding" units="fA"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_pulse</ci>
                <piecewise>
                    <piece>
                        <ci>pulseMag</ci>
                        <apply>
                            <and/>
                            <apply>
                                <gt/>
                                <ci>time</ci>
                                <ci>pulse_start</ci>
                            </apply>
                            <apply>
                                <lt/>
                                <ci>time</ci>
                                <ci>pulse_end</ci>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <ci>pulseHolding</ci>
                    </otherwise>
                </piecewise>
            </apply>
        </math>
        <variable name="v_E_RyR" units="fA"/>
        <variable name="sum_I" units="fA"/>
        <variable name="zCa" public_interface="in" units="dimensionless"/>
        <variable name="zK" public_interface="in" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>v_E_RyR</ci>
                <apply>
                    <times/>
                    <ci>zCa</ci>
                    <ci>v_RyR</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>sum_I</ci>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <ci>v_E_RyR</ci>
                        </apply>
                        <ci>v_E_LCC</ci>
                    </apply>
                    <ci>I_pulse</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_membrane</ci>
                </apply>
                <ci>sum_I</ci>
            </apply>
            <apply>
                <eq/>
                <ci>V_m</ci>
                <apply>
                    <divide/>
                    <ci>q_membrane</ci>
                    <ci>C_m</ci>
                </apply>
            </apply>
        </math>
        <!-- compute conserved quantities-->
        <!--        var Ca_T: fmol;-->
        <variable name="PLB_T" units="fmol"/>
        <variable name="PKACI_T" units="fmol"/>
        <variable name="Ip_T" units="fmol"/>
        <variable name="K_i_T" units="fmol"/>
        <variable name="K_o_T" units="fmol"/>
        <!--        Ca_T = q_Cai+q_Ca_TRPN+q_Ca_SR+q_Cao+q_Ca_fgate_LCC+q_SERCA_Ca_complexes+q_Ca_gate_complexes_mRyR;-->
        <!--+ q_Ca_di-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>PLB_T</ci>
                <apply>
                    <plus/>
                    <ci>q_PLB</ci>
                    <ci>q_PLB_PKACI</ci>
                    <ci>q_PLBp</ci>
                    <ci>q_PLBp_PP1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>PKACI_T</ci>
                <apply>
                    <plus/>
                    <ci>q_PKACI</ci>
                    <ci>q_PLB_PKACI</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Ip_T</ci>
                <apply>
                    <plus/>
                    <ci>q_Ip</ci>
                    <ci>q_Ip_PP1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K_i_T</ci>
                <ci>q_Ki</ci>
            </apply>
            <apply>
                <eq/>
                <ci>K_o_T</ci>
                <ci>q_Ko</ci>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="environment" component_2="crossbridge_TRPN"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="q_TRPN" variable_2="q_TRPN"/>
        <map_variables variable_1="q_Cai" variable_2="q_Cai"/>
        <map_variables variable_1="q_Ca_TRPN" variable_2="q_Ca_TRPN"/>
        <map_variables variable_1="q_B" variable_2="q_B"/>
        <map_variables variable_1="q_U" variable_2="q_U"/>
        <map_variables variable_1="q_W" variable_2="q_W"/>
        <map_variables variable_1="q_S" variable_2="q_S"/>
        <map_variables variable_1="v_R_TRPNCa" variable_2="v_R_TRPNCa"/>
        <map_variables variable_1="v_R_BU" variable_2="v_R_BU"/>
        <map_variables variable_1="v_R_UW" variable_2="v_R_UW"/>
        <map_variables variable_1="v_R_WS" variable_2="v_R_WS"/>
        <map_variables variable_1="v_R_SU" variable_2="v_R_SU"/>
        <map_variables variable_1="Tension" variable_2="T_total"/>
        <map_variables variable_1="SL" variable_2="SL"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="LCC"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="q_membrane" variable_2="q_E"/>
        <map_variables variable_1="q_Cai" variable_2="q_Cai"/>
        <map_variables variable_1="q_Cao" variable_2="q_Cao"/>
        <map_variables variable_1="q_Ki" variable_2="q_Ki"/>
        <map_variables variable_1="q_Ko" variable_2="q_Ko"/>
        <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="v_E_LCC" variable_2="v_E_LCC"/>
        <map_variables variable_1="I_pulse" variable_2="I_pulse"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="PLB"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="q_PLB" variable_2="q_PLB"/>
        <map_variables variable_1="q_PKACI" variable_2="q_PKACI"/>
        <map_variables variable_1="q_PLB_PKACI" variable_2="q_PLB_PKACI"/>
        <map_variables variable_1="q_PP1" variable_2="q_PP1"/>
        <map_variables variable_1="q_PLBp_PP1" variable_2="q_PLBp_PP1"/>
        <map_variables variable_1="q_PLBp" variable_2="q_PLBp"/>
        <map_variables variable_1="q_Ip" variable_2="q_Ip"/>
        <map_variables variable_1="q_Ip_PP1" variable_2="q_Ip_PP1"/>
        <map_variables variable_1="v_PLBph1" variable_2="v_PLBph1"/>
        <map_variables variable_1="v_PLBph2" variable_2="v_PLBph2"/>
        <map_variables variable_1="v_PLBd1" variable_2="v_PLBd1"/>
        <map_variables variable_1="v_PLBd2" variable_2="v_PLBd2"/>
        <map_variables variable_1="v_Inh" variable_2="v_Inh"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="RyR"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="V_m" variable_2="V_m"/>
        <map_variables variable_1="v_RyR" variable_2="v_RyR"/>
        <map_variables variable_1="v_gate_Cai" variable_2="v_gate_Cai"/>
        <map_variables variable_1="q_Ca_SR" variable_2="q_Ca_SR"/>
        <map_variables variable_1="q_Cai" variable_2="q_Cai"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="RyR_gate"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="SERCA"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="q_H" variable_2="q_H"/>
        <map_variables variable_1="q_Cai" variable_2="q_Cai"/>
        <map_variables variable_1="q_Ca_SR" variable_2="q_Ca_SR"/>
        <map_variables variable_1="q_MgATP" variable_2="q_MgATP"/>
        <map_variables variable_1="q_MgADP" variable_2="q_MgADP"/>
        <map_variables variable_1="q_P" variable_2="q_P"/>
        <map_variables variable_1="v_Cai_SERCA" variable_2="v_Cai"/>
        <map_variables variable_1="v_Ca_SR_SERCA" variable_2="v_Ca_SR"/>
        <map_variables variable_1="v_H_SERCA" variable_2="v_H"/>
        <map_variables variable_1="v_MgADP_SERCA" variable_2="v_MgADP"/>
        <map_variables variable_1="v_MgATP_SERCA" variable_2="v_MgATP"/>
        <map_variables variable_1="v_P_SERCA" variable_2="v_P"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="crossbridge_TRPN"/>
        <map_variables variable_1="kappa_R_TRPNCa" variable_2="kappa_R_TRPNCa"/>
        <map_variables variable_1="kappa_R_BU" variable_2="kappa_R_BU"/>
        <map_variables variable_1="kappa_R_UW" variable_2="kappa_R_UW"/>
        <map_variables variable_1="kappa_R_WS" variable_2="kappa_R_WS"/>
        <map_variables variable_1="kappa_R_SU" variable_2="kappa_R_SU"/>
        <map_variables variable_1="K_TRPN" variable_2="K_TRPN"/>
        <map_variables variable_1="K_Cai" variable_2="K_Cai"/>
        <map_variables variable_1="K_Ca_TRPN" variable_2="K_Ca_TRPN"/>
        <map_variables variable_1="K_B" variable_2="K_B"/>
        <map_variables variable_1="K_U" variable_2="K_U"/>
        <map_variables variable_1="K_W" variable_2="K_W"/>
        <map_variables variable_1="K_S" variable_2="K_S"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="LCC"/>
        <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"/>
        <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"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="LCC_gate"/>
        <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"/>
        <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="z_rCa" variable_2="z_rCa"/>
        <map_variables variable_1="z_fd" variable_2="z_fd"/>
        <map_variables variable_1="z_ff1" variable_2="z_ff1"/>
        <map_variables variable_1="z_ff2" variable_2="z_ff2"/>
        <map_variables variable_1="z_ff3" variable_2="z_ff3"/>
        <map_variables variable_1="z_rd" variable_2="z_rd"/>
        <map_variables variable_1="z_rf1" variable_2="z_rf1"/>
        <map_variables variable_1="z_rf2" variable_2="z_rf2"/>
        <map_variables variable_1="z_rf3" variable_2="z_rf3"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="PLB"/>
        <map_variables variable_1="kappa_PLBph1" variable_2="kappa_PLBph1"/>
        <map_variables variable_1="kappa_PLBph2" variable_2="kappa_PLBph2"/>
        <map_variables variable_1="kappa_PLBd1" variable_2="kappa_PLBd1"/>
        <map_variables variable_1="kappa_PLBd2" variable_2="kappa_PLBd2"/>
        <map_variables variable_1="kappa_Inh" variable_2="kappa_Inh"/>
        <map_variables variable_1="K_PLB" variable_2="K_PLB"/>
        <map_variables variable_1="K_PKACI" variable_2="K_PKACI"/>
        <map_variables variable_1="K_PLB_PKACI" variable_2="K_PLB_PKACI"/>
        <map_variables variable_1="K_PP1" variable_2="K_PP1"/>
        <map_variables variable_1="K_PLBp_PP1" variable_2="K_PLBp_PP1"/>
        <map_variables variable_1="K_PLBp" variable_2="K_PLBp"/>
        <map_variables variable_1="K_Ip" variable_2="K_Ip"/>
        <map_variables variable_1="K_Ip_PP1" variable_2="K_Ip_PP1"/>
    </connection>
    <connection>
        <map_components component_1="RyR" component_2="BG_parameters"/>
        <map_variables variable_1="kappa_RyR" variable_2="kappa_RyR"/>
        <map_variables variable_1="K_Ca_SR" variable_2="K_Ca_SR"/>
        <map_variables variable_1="K_Cai" variable_2="K_Cai"/>
    </connection>
    <connection>
        <map_components component_1="RyR_gate" component_2="BG_parameters"/>
        <map_variables variable_1="kappa_OC" variable_2="kappa_OC"/>
        <map_variables variable_1="kappa_CCI" variable_2="kappa_CCI"/>
        <map_variables variable_1="kappa_CII" variable_2="kappa_CII"/>
        <map_variables variable_1="kappa_IO" variable_2="kappa_IO"/>
        <map_variables variable_1="K_C_RyR" variable_2="K_C_RyR"/>
        <map_variables variable_1="K_CI_RyR" variable_2="K_CI_RyR"/>
        <map_variables variable_1="K_I_RyR" variable_2="K_I_RyR"/>
        <map_variables variable_1="K_O_RyR" variable_2="K_O_RyR"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="SERCA"/>
        <map_variables variable_1="kappa_SERCA_R1_2" variable_2="kappa_SERCA_R1_2"/>
        <map_variables variable_1="kappa_SERCA_R2_4" variable_2="kappa_SERCA_R2_4"/>
        <map_variables variable_1="kappa_SERCA_R2_2a" variable_2="kappa_SERCA_R2_2a"/>
        <map_variables variable_1="kappa_SERCA_R4_5" variable_2="kappa_SERCA_R4_5"/>
        <map_variables variable_1="kappa_SERCA_R5_6" variable_2="kappa_SERCA_R5_6"/>
        <map_variables variable_1="kappa_SERCA_R6_8" variable_2="kappa_SERCA_R6_8"/>
        <map_variables variable_1="kappa_SERCA_R8_9" variable_2="kappa_SERCA_R8_9"/>
        <map_variables variable_1="kappa_SERCA_R9_10" variable_2="kappa_SERCA_R9_10"/>
        <map_variables variable_1="kappa_SERCA_R10_1" variable_2="kappa_SERCA_R10_1"/>
        <map_variables variable_1="K_P1_SERCA" variable_2="K_P1_SERCA"/>
        <map_variables variable_1="K_P2_SERCA" variable_2="K_P2_SERCA"/>
        <map_variables variable_1="K_P2a_SERCA" variable_2="K_P2a_SERCA"/>
        <map_variables variable_1="K_P4_SERCA" variable_2="K_P4_SERCA"/>
        <map_variables variable_1="K_P5_SERCA" variable_2="K_P5_SERCA"/>
        <map_variables variable_1="K_P6_SERCA" variable_2="K_P6_SERCA"/>
        <map_variables variable_1="K_P8_SERCA" variable_2="K_P8_SERCA"/>
        <map_variables variable_1="K_P9_SERCA" variable_2="K_P9_SERCA"/>
        <map_variables variable_1="K_P10_SERCA" variable_2="K_P10_SERCA"/>
        <map_variables variable_1="K_H" variable_2="K_H"/>
        <map_variables variable_1="K_Cai" variable_2="K_Cai"/>
        <map_variables variable_1="K_Ca_SR" variable_2="K_Ca_SR"/>
        <map_variables variable_1="K_MgATP" variable_2="K_MgATP"/>
        <map_variables variable_1="K_MgADP" variable_2="K_MgADP"/>
        <map_variables variable_1="K_P" variable_2="K_P"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="crossbridge_TRPN"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="LCC"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="PLB"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="RyR"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="RyR_gate"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="SERCA"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="ion_valences" component_2="environment"/>
        <map_variables variable_1="zCa" variable_2="zCa"/>
        <map_variables variable_1="zK" variable_2="zK"/>
    </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="constants" component_2="environment"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="ion_valences" component_2="LCC"/>
        <map_variables variable_1="zCa" variable_2="zCa"/>
        <map_variables variable_1="zK" variable_2="zK"/>
    </connection>
    <connection>
        <map_components component_1="ion_valences" component_2="RyR"/>
        <map_variables variable_1="zCa" variable_2="zCa"/>
    </connection>
</model>