- Author:
- Shelley Fong <s.fong@auckland.ac.nz>
- Date:
- 2021-11-11 15:41:15+13:00
- Desc:
- io
- Permanent Source URI:
- https://models.cellml.org/workspace/7a8/rawfile/633993335c6e11ad2fa206a1eeeeaa1714a64233/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="PLB/BG_PLB.cellml">
<component component_ref="PLB" name="PLB"/>
</import>
<import xlink:href="SERCA/BG_SERCA.cellml">
<component component_ref="SERCA" name="SERCA"/>
</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="LCC/BG_LCC.cellml">
<component component_ref="LCC" name="LCC"/>
<component component_ref="LCC_gate" name="LCC_gate"/>
</import>
<import xlink:href="TRPN/BG_TRPN.cellml">
<component component_ref="TRPN" name="TRPN"/>
</import>
<component name="BG_parameters">
<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="2.24963e-05" name="kappa_SERCA_R1_2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="108329" name="kappa_SERCA_R2_4" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="2.77146e+06" name="kappa_SERCA_R2_2a" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="108329" name="kappa_SERCA_R4_5" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="3442.07" name="kappa_SERCA_R5_6" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="3.07427e+06" name="kappa_SERCA_R6_8" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.0306e+07" name="kappa_SERCA_R8_9" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.0306e+07" name="kappa_SERCA_R9_10" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00565005" name="kappa_SERCA_R10_1" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="3.58172e+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="8.37314" name="kappa_Ca1" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="15.6575" name="kappa_Ca2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.0929479" name="kappa_K1" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.17381" name="kappa_K2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="942.411" name="kappa_d000" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="23.2508" name="kappa_d010" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="43.4783" name="kappa_d020" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.175516" name="kappa_d001" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00433026" name="kappa_d011" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00809745" name="kappa_d021" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.73616" name="kappa_f1_000" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="8.62132" name="kappa_f1_100" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.000323344" name="kappa_f1_001" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00160565" name="kappa_f1_101" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="3.17968" name="kappa_f2_000" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="15.7895" name="kappa_f2_100" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.000592187" name="kappa_f2_001" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.00294065" name="kappa_f2_101" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="8932.14" name="kappa_f3_010" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="44354.8" name="kappa_f3_110" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.66353" name="kappa_f3_011" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="8.2607" name="kappa_f3_111" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="60096.5" name="kappa_fCa000" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="298424" name="kappa_fCa100" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1482.67" name="kappa_fCa010" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="7362.59" name="kappa_fCa110" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="2772.56" name="kappa_fCa020" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="13767.8" name="kappa_fCa120" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="223.356" name="kappa_R_TRPNCa" public_interface="out" units="fmol_per_sec"/>
<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="906.879" name="K_P1_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="2605.18" name="K_P2_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="271.664" name="K_P2a_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="6950.17" name="K_P4_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0214514" name="K_P5_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="1228.84" name="K_P6_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="76.0193" name="K_P8_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="73.055" name="K_P9_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="105.473" name="K_P10_SERCA" public_interface="out" units="per_fmol"/>
<variable initial_value="278.106" name="K_H" public_interface="out" units="per_fmol"/>
<variable initial_value="0.052177" name="K_Cai" public_interface="out" units="per_fmol"/>
<variable initial_value="0.052177" name="K_Ca_SR" public_interface="out" units="per_fmol"/>
<variable initial_value="1072.81" name="K_MgATP" public_interface="out" units="per_fmol"/>
<variable initial_value="1.34257e-05" name="K_MgADP" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0245735" name="K_P" 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="0.052177" name="K_Cao" public_interface="out" units="per_fmol"/>
<variable initial_value="0.00369486" name="K_Ki" public_interface="out" units="per_fmol"/>
<variable initial_value="0.00369486" name="K_Ko" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0150147" name="K_000_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="0.608585" name="K_010_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="0.325451" name="K_020_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="0.00302366" name="K_100_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="0.122556" name="K_110_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0655392" name="K_120_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="80.6198" name="K_001_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="3267.72" name="K_011_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="1747.47" name="K_021_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="16.2352" name="K_101_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="658.051" name="K_111_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="351.905" name="K_121_LCC" public_interface="out" units="per_fmol"/>
<variable initial_value="7.25115" name="K_TRPN" public_interface="out" units="per_fmol"/>
<variable initial_value="0.02603" name="K_Ca_TRPN" 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 initial_value="1.381e5" name="C_m" public_interface="out" units="fF"/>
<variable initial_value="-13039" name="q_membrane" public_interface="out" units="fC"/>
<variable name="V_m" public_interface="out" units="volt"/>
<variable name="v_RyR" public_interface="in" units="fmol_per_sec"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<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="6.82e-3" name="q_Cai_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"/>
<!-- var q_Ca_di_init: fmol {init: 3.09e-6};-->
<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="2.57" name="q_TRPN_init" units="fmol"/>
<variable initial_value="1e-18" name="q_Ca_TRPN_init" units="fmol"/>
<!-- Global value-->
<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_H" public_interface="out" units="fmol"/>
<variable name="q_Cai" 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"/>
<variable name="q_Ca_SR" public_interface="out" units="fmol"/>
<!-- var q_Ca_di: fmol {pub: out};-->
<!-- var q_C_RyR: fmol {pub: out};
var q_CI_RyR: fmol {pub: out};
var q_I_RyR: fmol {pub: out};
var q_O_RyR: fmol {pub: out};-->
<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_TRPN" public_interface="out" units="fmol"/>
<variable name="q_Ca_TRPN" public_interface="out" 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"/>
<!-- 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"/>
<!-- RyR imports-->
<variable name="q_Ca_SR_mRyR" public_interface="in" units="fmol"/>
<variable name="q_Cai_mRyR" public_interface="in" units="fmol"/>
<!-- RyR gate imports-->
<variable name="q_Ca_gate_complexes_mRyR" public_interface="in" units="fmol"/>
<!-- 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"/>
<!-- LCC_gate imports-->
<variable name="q_Ca_fgate_LCC" public_interface="in" units="fmol"/>
<!-- TRPN imports-->
<variable name="q_Cai_mTRPN" public_interface="in" units="fmol"/>
<variable name="q_TRPN_mTRPN" public_interface="in" units="fmol"/>
<variable name="q_Ca_TRPN_mTRPN" public_interface="in" units="fmol"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<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_H</ci>
<apply>
<plus/>
<ci>q_H_init</ci>
<ci>q_H_mSERCA</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Cai</ci>
<apply>
<plus/>
<ci>q_Cai_init</ci>
<ci>q_Cai_mSERCA</ci>
<ci>q_Cai_mLCC</ci>
<ci>q_Cai_mTRPN</ci>
<ci>q_Cai_mRyR</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Ca_SR</ci>
<apply>
<plus/>
<ci>q_Ca_SR_init</ci>
<ci>q_Ca_SR_mSERCA</ci>
<ci>q_Ca_SR_mRyR</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>
<!-- q_Ca_di = q_Ca_di_init+q_Ca_di_mRyR;-->
<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_TRPN</ci>
<apply>
<plus/>
<ci>q_TRPN_init</ci>
<ci>q_TRPN_mTRPN</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Ca_TRPN</ci>
<apply>
<plus/>
<ci>q_Ca_TRPN_init</ci>
<ci>q_Ca_TRPN_mTRPN</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="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="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="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="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"/>
</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="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="TRPN"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="q_Cai_mTRPN" variable_2="q_Cai"/>
<map_variables variable_1="q_Cai" variable_2="q_Cai_global"/>
<map_variables variable_1="q_TRPN_mTRPN" variable_2="q_TRPN"/>
<map_variables variable_1="q_TRPN" variable_2="q_TRPN_global"/>
<map_variables variable_1="q_Ca_TRPN_mTRPN" variable_2="q_Ca_TRPN"/>
<map_variables variable_1="q_Ca_TRPN" variable_2="q_Ca_TRPN_global"/>
</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="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="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="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="TRPN"/>
<map_variables variable_1="kappa_R_TRPNCa" variable_2="kappa_R_TRPNCa"/>
<map_variables variable_1="K_Cai" variable_2="K_Cai"/>
<map_variables variable_1="K_TRPN" variable_2="K_TRPN"/>
<map_variables variable_1="K_Ca_TRPN" variable_2="K_Ca_TRPN"/>
</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="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="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="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="TRPN"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
</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="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="RyR"/>
<map_variables variable_1="zCa" variable_2="zCa"/>
</connection>
</model>
