- Author:
- Shelley Fong <s.fong@auckland.ac.nz>
- Date:
- 2022-03-10 10:16:08+13:00
- Desc:
- With updated parameters following W fix, and updated RyR component
- Permanent Source URI:
- https://models.cellml.org/workspace/7a8/rawfile/2c52cc3fe8538fba448ec1255ffa3934b951ca50/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"/>
</import>
<import xlink:href="SERCA/BG_SERCA.cellml">
<component component_ref="SERCA" name="SERCA"/>
</import>
<units name="kPa">
<unit prefix="kilo" units="pascal"/>
</units>
<component name="BG_parameters">
<variable initial_value="117.289" name="kappa_R_TRPNCa" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="6.35843" name="kappa_R_BU" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.297353" name="kappa_R_UW" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00699655" name="kappa_R_WS" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="125938" name="kappa_R_SU" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="9.98335" name="kappa_Ca1" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="18.6686" name="kappa_Ca2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.0964275" name="kappa_K1" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.180317" name="kappa_K2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="826.484" name="kappa_d000" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="20.3907" name="kappa_d010" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="38.13" name="kappa_d020" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.188903" name="kappa_d001" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00466054" name="kappa_d011" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00871509" name="kappa_d021" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.52259" name="kappa_f1_000" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="7.56081" name="kappa_f1_100" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.000348007" name="kappa_f1_001" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00172812" name="kappa_f1_101" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="2.78854" name="kappa_f2_000" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="13.8472" name="kappa_f2_100" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.000637356" name="kappa_f2_001" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00316495" name="kappa_f2_101" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="7833.39" name="kappa_f3_010" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="38898.7" name="kappa_f3_110" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.79042" name="kappa_f3_011" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="8.89078" name="kappa_f3_111" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="64680.3" name="kappa_fCa000" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="321186" name="kappa_fCa100" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1595.76" name="kappa_fCa010" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="7924.17" name="kappa_fCa110" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="2984.03" name="kappa_fCa020" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="14818" 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="1.87783e+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="1.83036e-05" name="kappa_SERCA_R1_2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="108168" name="kappa_SERCA_R2_4" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="2.0993e+06" name="kappa_SERCA_R2_2a" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="108168" name="kappa_SERCA_R4_5" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="3965.46" name="kappa_SERCA_R5_6" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.35535e+06" name="kappa_SERCA_R6_8" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.43269e+07" name="kappa_SERCA_R8_9" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.43269e+07" name="kappa_SERCA_R9_10" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.0084367" name="kappa_SERCA_R10_1" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="15.2972" name="K_TRPN" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0470993" name="K_Cai" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0495694" name="K_Ca_TRPN" public_interface="out" units="per_fmol"/>
<variable initial_value="6.77815" name="K_B" public_interface="out" units="per_fmol"/>
<variable initial_value="2.5418" name="K_U" public_interface="out" units="per_fmol"/>
<variable initial_value="16.6195" name="K_W" public_interface="out" units="per_fmol"/>
<variable initial_value="4.15487e-06" name="K_S" public_interface="out" units="per_fmol"/>
<variable initial_value="0.312662" name="K_Cao" public_interface="out" units="per_fmol"/>
<variable initial_value="0.00383319" name="K_Ki" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0254461" name="K_Ko" public_interface="out" units="per_fmol"/>
<variable initial_value="0.588955" name="K_000_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="23.8718" name="K_010_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="12.7659" name="K_020_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="0.118603" name="K_100_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="4.80729" name="K_110_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="2.57078" name="K_120_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="2576.78" name="K_001_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="104443" name="K_011_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="55852.8" name="K_021_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="518.91" name="K_101_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="21032.7" name="K_111_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="11247.6" 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="2.84345" name="K_Ca_SR" public_interface="out" units="per_fmol"/>
<variable initial_value="1.8464" name="K_Ca_i" public_interface="out" units="per_fmol"/>
<variable initial_value="3962.58" name="K_C_RyR" public_interface="out" units="per_fmol"/>
<variable initial_value="39.6258" name="K_CI_RyR" public_interface="out" units="per_fmol"/>
<variable initial_value="0.06793" name="K_I_RyR" public_interface="out" units="per_fmol"/>
<variable initial_value="6.793" name="K_O_RyR" public_interface="out" units="per_fmol"/>
<variable initial_value="29236.5" name="K_P1_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="110146" name="K_P2_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="12337.4" name="K_P2a_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="239441" name="K_P4_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="0.640531" name="K_P5_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="95883.3" name="K_P6_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="1630.43" name="K_P8_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="1807.79" name="K_P9_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="2429.86" name="K_P10_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="288.622" name="K_H" public_interface="out" units="per_fmol"/>
<variable initial_value="1406.95" name="K_MgATP" public_interface="out" units="per_fmol"/>
<variable initial_value="5.13777e-06" name="K_MgADP" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0175602" 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_RyRgate_Ca_i" 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_RyRgate_Ca_i</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="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>sum_I</ci>
<apply>
<plus/>
<apply>
<minus/>
<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="t"/>
<map_variables variable_1="v_RyR" variable_2="v_RyR"/>
<map_variables variable_1="v_RyRgate_Ca_i" variable_2="v_RyRgate_Ca_i"/>
<map_variables variable_1="q_Ca_SR" variable_2="q_Ca_SR"/>
<map_variables variable_1="q_Cai" variable_2="q_Ca_i"/>
</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_Ca_i" variable_2="K_Cai"/>
<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="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>
</model>
