- Author:
- Shelley Fong <s.fong@auckland.ac.nz>
- Date:
- 2022-01-28 09:31:09+13:00
- Desc:
- Adding working MATLAB file
- Permanent Source URI:
- https://models.cellml.org/workspace/8a0/rawfile/de0f4010f8be6cebc02cb3ea514727a337f2f779/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"/>
</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 name="v_RyR" public_interface="in" units="fmol_per_sec"/>
<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_init" units="fmol"/>
<variable initial_value="1e-18" name="q_Ca_TRPN_init" units="fmol"/>
<variable initial_value="6.82e-3" name="q_Cai_init" units="fmol"/>
<variable initial_value="1" name="q_B_init" units="fmol"/>
<variable initial_value="0" name="q_U_init" units="fmol"/>
<variable initial_value="0" name="q_W_init" units="fmol"/>
<variable initial_value="0" name="q_S_init" units="fmol"/>
<variable initial_value="6.84" name="q_Cao_init" units="fmol"/>
<variable initial_value="5.51E+03" name="q_Ki_init" units="fmol"/>
<variable initial_value="2.05E+02" name="q_Ko_init" units="fmol"/>
<variable initial_value="4.028E+00" name="q_PLB_init" units="fmol"/>
<variable initial_value="2.234E-03" name="q_PKACI_init" units="fmol"/>
<variable initial_value="1e-18" name="q_PLB_PKACI_init" units="fmol"/>
<variable initial_value="3.382E-02" name="q_PP1_init" units="fmol"/>
<variable initial_value="1e-18" name="q_PLBp_PP1_init" units="fmol"/>
<variable initial_value="1e-18" name="q_PLBp_init" units="fmol"/>
<variable initial_value="1.999E-03" name="q_Ip_init" units="fmol"/>
<variable initial_value="1e-18" name="q_Ip_PP1_init" units="fmol"/>
<variable initial_value="0.641" name="q_Ca_SR_init" units="fmol"/>
<variable initial_value="0.004028" name="q_H_init" units="fmol"/>
<variable initial_value="570" name="q_P_init" units="fmol"/>
<variable initial_value="1.3794" name="q_MgADP_init" units="fmol"/>
<variable initial_value="3.8" name="q_MgATP_init" units="fmol"/>
<!-- mass conservation checks-->
<!-- Global value-->
<variable name="q_TRPN" public_interface="out" units="fmol"/>
<variable name="q_Cai" public_interface="out" units="fmol"/>
<variable name="q_Ca_TRPN" public_interface="out" units="fmol"/>
<variable name="q_B" public_interface="out" units="fmol"/>
<variable name="q_U" public_interface="out" units="fmol"/>
<variable name="q_W" public_interface="out" units="fmol"/>
<variable name="q_S" public_interface="out" units="fmol"/>
<variable name="q_Cao" public_interface="out" units="fmol"/>
<variable name="q_Ki" public_interface="out" units="fmol"/>
<variable name="q_Ko" public_interface="out" units="fmol"/>
<variable name="q_PLB" public_interface="out" units="fmol"/>
<variable name="q_PKACI" public_interface="out" units="fmol"/>
<variable name="q_PLB_PKACI" public_interface="out" units="fmol"/>
<variable name="q_PP1" public_interface="out" units="fmol"/>
<variable name="q_PLBp_PP1" public_interface="out" units="fmol"/>
<variable name="q_PLBp" public_interface="out" units="fmol"/>
<variable name="q_Ip" public_interface="out" units="fmol"/>
<variable name="q_Ip_PP1" public_interface="out" units="fmol"/>
<variable name="q_Ca_SR" public_interface="out" units="fmol"/>
<variable name="q_H" public_interface="out" units="fmol"/>
<variable name="q_MgATP" public_interface="out" units="fmol"/>
<variable name="q_MgADP" public_interface="out" units="fmol"/>
<variable name="q_P" public_interface="out" units="fmol"/>
<!-- crossbridge_TRPN imports-->
<variable name="q_TRPN_mcrossbridge_TRPN" public_interface="in" units="fmol"/>
<variable name="q_Cai_mcrossbridge_TRPN" public_interface="in" units="fmol"/>
<variable name="q_Ca_TRPN_mcrossbridge_TRPN" public_interface="in" units="fmol"/>
<variable name="q_B_mcrossbridge_TRPN" public_interface="in" units="fmol"/>
<variable name="q_U_mcrossbridge_TRPN" public_interface="in" units="fmol"/>
<variable name="q_W_mcrossbridge_TRPN" public_interface="in" units="fmol"/>
<variable name="q_S_mcrossbridge_TRPN" public_interface="in" units="fmol"/>
<variable name="Tension" public_interface="in" units="kPa"/>
<variable name="SL" public_interface="in" units="metre"/>
<!-- LCC imports-->
<variable name="q_Cai_mLCC" public_interface="in" units="fmol"/>
<variable name="q_Cao_mLCC" public_interface="in" units="fmol"/>
<variable name="q_Ki_mLCC" public_interface="in" units="fmol"/>
<variable name="q_Ko_mLCC" public_interface="in" units="fmol"/>
<variable name="v_Ca_fgate_mLCC" public_interface="in" units="fmol_per_sec"/>
<variable name="q_Ca_fgate_LCC" public_interface="in" units="fmol"/>
<!-- PLB imports-->
<variable name="q_PLB_mPLB" public_interface="in" units="fmol"/>
<variable name="q_PKACI_mPLB" public_interface="in" units="fmol"/>
<variable name="q_PLB_PKACI_mPLB" public_interface="in" units="fmol"/>
<variable name="q_PP1_mPLB" public_interface="in" units="fmol"/>
<variable name="q_PLBp_PP1_mPLB" public_interface="in" units="fmol"/>
<variable name="q_PLBp_mPLB" public_interface="in" units="fmol"/>
<variable name="q_Ip_mPLB" public_interface="in" units="fmol"/>
<variable name="q_Ip_PP1_mPLB" public_interface="in" units="fmol"/>
<!-- RyR imports-->
<variable name="q_Ca_SR_mRyR" public_interface="in" units="fmol"/>
<variable name="q_Cai_mRyR" public_interface="in" units="fmol"/>
<variable name="q_Ca_gate_complexes_mRyR" public_interface="in" units="fmol"/>
<!-- SERCA imports-->
<variable name="q_H_mSERCA" public_interface="in" units="fmol"/>
<variable name="q_Cai_mSERCA" public_interface="in" units="fmol"/>
<variable name="q_Ca_SR_mSERCA" public_interface="in" units="fmol"/>
<variable name="q_MgATP_mSERCA" public_interface="in" units="fmol"/>
<variable name="q_MgADP_mSERCA" public_interface="in" units="fmol"/>
<variable name="q_P_mSERCA" public_interface="in" units="fmol"/>
<variable name="q_SERCA_Ca_complexes" public_interface="in" units="fmol"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>q_TRPN</ci>
<apply>
<plus/>
<ci>q_TRPN_init</ci>
<ci>q_TRPN_mcrossbridge_TRPN</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Cai</ci>
<apply>
<plus/>
<ci>q_Cai_init</ci>
<ci>q_Cai_mcrossbridge_TRPN</ci>
<ci>q_Cai_mLCC</ci>
<ci>q_Cai_mRyR</ci>
<ci>q_Cai_mSERCA</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Ca_TRPN</ci>
<apply>
<plus/>
<ci>q_Ca_TRPN_init</ci>
<ci>q_Ca_TRPN_mcrossbridge_TRPN</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_B</ci>
<apply>
<plus/>
<ci>q_B_init</ci>
<ci>q_B_mcrossbridge_TRPN</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_U</ci>
<apply>
<plus/>
<ci>q_U_init</ci>
<ci>q_U_mcrossbridge_TRPN</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_W</ci>
<apply>
<plus/>
<ci>q_W_init</ci>
<ci>q_W_mcrossbridge_TRPN</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_S</ci>
<apply>
<plus/>
<ci>q_S_init</ci>
<ci>q_S_mcrossbridge_TRPN</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Cao</ci>
<apply>
<plus/>
<ci>q_Cao_init</ci>
<ci>q_Cao_mLCC</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Ki</ci>
<apply>
<plus/>
<ci>q_Ki_init</ci>
<ci>q_Ki_mLCC</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Ko</ci>
<apply>
<plus/>
<ci>q_Ko_init</ci>
<ci>q_Ko_mLCC</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_PLB</ci>
<apply>
<plus/>
<ci>q_PLB_init</ci>
<ci>q_PLB_mPLB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_PKACI</ci>
<apply>
<plus/>
<ci>q_PKACI_init</ci>
<ci>q_PKACI_mPLB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_PLB_PKACI</ci>
<apply>
<plus/>
<ci>q_PLB_PKACI_init</ci>
<ci>q_PLB_PKACI_mPLB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_PP1</ci>
<apply>
<plus/>
<ci>q_PP1_init</ci>
<ci>q_PP1_mPLB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_PLBp_PP1</ci>
<apply>
<plus/>
<ci>q_PLBp_PP1_init</ci>
<ci>q_PLBp_PP1_mPLB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_PLBp</ci>
<apply>
<plus/>
<ci>q_PLBp_init</ci>
<ci>q_PLBp_mPLB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Ip</ci>
<apply>
<plus/>
<ci>q_Ip_init</ci>
<ci>q_Ip_mPLB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Ip_PP1</ci>
<apply>
<plus/>
<ci>q_Ip_PP1_init</ci>
<ci>q_Ip_PP1_mPLB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Ca_SR</ci>
<apply>
<plus/>
<ci>q_Ca_SR_init</ci>
<ci>q_Ca_SR_mRyR</ci>
<ci>q_Ca_SR_mSERCA</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H</ci>
<apply>
<plus/>
<ci>q_H_init</ci>
<ci>q_H_mSERCA</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_MgATP</ci>
<apply>
<plus/>
<ci>q_MgATP_init</ci>
<ci>q_MgATP_mSERCA</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_MgADP</ci>
<apply>
<plus/>
<ci>q_MgADP_init</ci>
<ci>q_MgADP_mSERCA</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_P</ci>
<apply>
<plus/>
<ci>q_P_init</ci>
<ci>q_P_mSERCA</ci>
</apply>
</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="v_E_LCC" public_interface="in" 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-->
<variable name="Ca_T" units="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"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Ca_T</ci>
<apply>
<plus/>
<ci>q_Cai</ci>
<ci>q_Ca_TRPN</ci>
<ci>q_Ca_SR</ci>
<ci>q_Cao</ci>
<ci>q_Ca_fgate_LCC</ci>
<ci>q_SERCA_Ca_complexes</ci>
<ci>q_Ca_gate_complexes_mRyR</ci>
</apply>
</apply>
<!--+ q_Ca_di-->
<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_mcrossbridge_TRPN" variable_2="q_TRPN"/>
<map_variables variable_1="q_TRPN" variable_2="q_TRPN_global"/>
<map_variables variable_1="q_Cai_mcrossbridge_TRPN" variable_2="q_Cai"/>
<map_variables variable_1="q_Cai" variable_2="q_Cai_global"/>
<map_variables variable_1="q_Ca_TRPN_mcrossbridge_TRPN" variable_2="q_Ca_TRPN"/>
<map_variables variable_1="q_Ca_TRPN" variable_2="q_Ca_TRPN_global"/>
<map_variables variable_1="q_B_mcrossbridge_TRPN" variable_2="q_B"/>
<map_variables variable_1="q_B" variable_2="q_B_global"/>
<map_variables variable_1="q_U_mcrossbridge_TRPN" variable_2="q_U"/>
<map_variables variable_1="q_U" variable_2="q_U_global"/>
<map_variables variable_1="q_W_mcrossbridge_TRPN" variable_2="q_W"/>
<map_variables variable_1="q_W" variable_2="q_W_global"/>
<map_variables variable_1="q_S_mcrossbridge_TRPN" variable_2="q_S"/>
<map_variables variable_1="q_S" variable_2="q_S_global"/>
<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="C_m" variable_2="C_m"/>
<map_variables variable_1="q_membrane" variable_2="q_E"/>
<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="v_E_LCC" variable_2="v_E_LCC"/>
<map_variables variable_1="I_pulse" variable_2="I_pulse"/>
<map_variables variable_1="v_Ca_fgate_mLCC" variable_2="v_Ca_fgate"/>
</connection>
<connection>
<map_components component_1="environment" component_2="PLB"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="q_PLB_mPLB" variable_2="q_PLB"/>
<map_variables variable_1="q_PLB" variable_2="q_PLB_global"/>
<map_variables variable_1="q_PKACI_mPLB" variable_2="q_PKACI"/>
<map_variables variable_1="q_PKACI" variable_2="q_PKACI_global"/>
<map_variables variable_1="q_PLB_PKACI_mPLB" variable_2="q_PLB_PKACI"/>
<map_variables variable_1="q_PLB_PKACI" variable_2="q_PLB_PKACI_global"/>
<map_variables variable_1="q_PP1_mPLB" variable_2="q_PP1"/>
<map_variables variable_1="q_PP1" variable_2="q_PP1_global"/>
<map_variables variable_1="q_PLBp_PP1_mPLB" variable_2="q_PLBp_PP1"/>
<map_variables variable_1="q_PLBp_PP1" variable_2="q_PLBp_PP1_global"/>
<map_variables variable_1="q_PLBp_mPLB" variable_2="q_PLBp"/>
<map_variables variable_1="q_PLBp" variable_2="q_PLBp_global"/>
<map_variables variable_1="q_Ip_mPLB" variable_2="q_Ip"/>
<map_variables variable_1="q_Ip" variable_2="q_Ip_global"/>
<map_variables variable_1="q_Ip_PP1_mPLB" variable_2="q_Ip_PP1"/>
<map_variables variable_1="q_Ip_PP1" variable_2="q_Ip_PP1_global"/>
</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="q_Ca_SR_mRyR" variable_2="q_Ca_SR"/>
<map_variables variable_1="q_Ca_SR" variable_2="q_Ca_SR_global"/>
<map_variables variable_1="q_Cai_mRyR" variable_2="q_Cai"/>
<map_variables variable_1="q_Cai" variable_2="q_Cai_global"/>
<map_variables variable_1="v_RyR" variable_2="v_RyR"/>
</connection>
<connection>
<map_components component_1="environment" component_2="RyR_gate"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="q_Ca_gate_complexes_mRyR" variable_2="q_Ca_gate_complexes"/>
</connection>
<connection>
<map_components component_1="environment" component_2="SERCA"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="q_H_mSERCA" variable_2="q_H"/>
<map_variables variable_1="q_H" variable_2="q_H_global"/>
<map_variables variable_1="q_Cai_mSERCA" variable_2="q_Cai"/>
<map_variables variable_1="q_Cai" variable_2="q_Cai_global"/>
<map_variables variable_1="q_Ca_SR_mSERCA" variable_2="q_Ca_SR"/>
<map_variables variable_1="q_Ca_SR" variable_2="q_Ca_SR_global"/>
<map_variables variable_1="q_MgATP_mSERCA" variable_2="q_MgATP"/>
<map_variables variable_1="q_MgATP" variable_2="q_MgATP_global"/>
<map_variables variable_1="q_MgADP_mSERCA" variable_2="q_MgADP"/>
<map_variables variable_1="q_MgADP" variable_2="q_MgADP_global"/>
<map_variables variable_1="q_P_mSERCA" variable_2="q_P"/>
<map_variables variable_1="q_P" variable_2="q_P_global"/>
<map_variables variable_1="q_SERCA_Ca_complexes" variable_2="q_SERCA_Ca_complexes"/>
</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"/>
</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"/>
<map_variables variable_1="q_Ca_fgate_LCC" variable_2="q_Ca_fgate"/>
</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>
