Location: Hunt, Tilunaite, Bass, Soeller, Roderick, Rajagopal, Crampin @ 348331e8cfcf / Hunt_et_al_ECC_ETC_calcium_dynamics_cytosol.cellml
- Author:
- Hilary Hunt <hhunt1@6200L-124424-M.local>
- Date:
- 2020-05-06 15:54:20+10:00
- Desc:
- readded cellml
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/5ee/rawfile/348331e8cfcf170ad73d8831fa6ecb1f69454ef6/Hunt_et_al_ECC_ETC_calcium_dynamics_cytosol.cellml
<?xml version='1.0'?>
<model cmeta:id="Hinch_et_al_2004" name="Hinch_et_al_2004" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<units name="per_ms">
<unit exponent="-1" prefix="milli" units="second"/>
</units>
<units name="ms">
<unit prefix="milli" units="second"/>
</units>
<units name="mV">
<unit prefix="milli" units="volt"/>
</units>
<units name="mM">
<unit prefix="milli" units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="mM_per_ms">
<unit units="mM"/>
<unit exponent="-1" prefix="milli" units="second"/>
</units>
<units name="per_mM_ms">
<unit exponent="-1" units="mM"/>
<unit exponent="-1" prefix="milli" units="second"/>
</units>
<units name="mM_per_mV_ms">
<unit units="mM"/>
<unit exponent="-1" prefix="milli" units="second"/>
<unit exponent="-1" prefix="milli" units="volt"/>
</units>
<units name="per_ms3">
<unit exponent="-3" prefix="milli" units="second"/>
</units>
<units name="um2">
<unit exponent="2" prefix="micro" units="meter"/>
</units>
<units name="um3">
<unit exponent="3" prefix="micro" units="meter"/>
</units>
<units name="um3_per_ms">
<unit exponent="3" prefix="micro" units="meter"/>
<unit exponent="-1" prefix="milli" units="second"/>
</units>
<units name="um3_mM_per_ms">
<unit exponent="3" prefix="micro" units="meter"/>
<unit units="mM"/>
<unit exponent="-1" prefix="milli" units="second"/>
</units>
<units name="mJ_per_mole_K">
<unit prefix="milli" units="joule"/>
<unit exponent="-1" units="mole"/>
<unit exponent="-1" units="kelvin"/>
</units>
<units name="C_per_mole">
<unit units="coulomb"/>
<unit exponent="-1" units="mole"/>
</units>
<units name="first_order_rate_constant">
<unit exponent="-1" units="second"/>
</units>
<units name="second_order_rate_constant">
<unit exponent="-1" units="mM"/>
<unit exponent="-1" units="second"/>
</units>
<component name="environment">
<variable name="time" public_interface="out" units="ms"/>
</component>
<component name="cell_geometry">
<variable initial_value="25.84e3" name="V_myo" public_interface="out" units="um3"/>
<variable initial_value="2.098e3" name="V_SR" public_interface="out" units="um3"/>
<variable initial_value="1.534e4" name="A_cap" public_interface="out" units="um2"/>
</component>
<component name="membrane">
<variable name="V" public_interface="out" units="mV"/>
<variable initial_value="8314.5" name="R" public_interface="out" units="mJ_per_mole_K"/>
<variable initial_value="295" name="T" public_interface="out" units="kelvin"/>
<variable initial_value="96487" name="F" public_interface="out" units="C_per_mole"/>
<variable name="time" public_interface="in" units="ms"/>
<variable name="FVRT" public_interface="out" units="dimensionless"/>
<variable name="FVRT_Ca" public_interface="out" units="dimensionless"/>
<variable initial_value="1000" name="beat_period" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>FVRT</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>F</ci>
<ci>V</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>FVRT_Ca</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>FVRT</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V</ci>
<piecewise>
<piece>
<apply>
<minus/>
<cn cellml:units="mV">35</cn>
<apply>
<times/>
<cn cellml:units="mV">1.4375</cn>
<apply>
<rem/>
<ci>time</ci>
<ci>beat_period</ci>
</apply>
</apply>
</apply>
<apply>
<and/>
<apply>
<geq/>
<apply>
<rem/>
<ci>time</ci>
<ci>beat_period</ci>
</apply>
<cn cellml:units="ms">0</cn>
</apply>
<apply>
<leq/>
<apply>
<rem/>
<ci>time</ci>
<ci>beat_period</ci>
</apply>
<cn cellml:units="ms">80</cn>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<minus/>
<cn cellml:units="mV">80</cn>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="CaRU">
<variable initial_value="0.065" name="g_D" private_interface="out" units="um3_per_ms"/>
<variable initial_value="0.02" name="J_R" private_interface="out" units="um3_per_ms"/>
<variable initial_value="9.13e-4" name="J_L" private_interface="out" units="um3_per_ms"/>
<variable name="V_myo" private_interface="out" public_interface="in" units="um3"/>
<variable initial_value="50000" name="N" private_interface="out" units="dimensionless"/>
<variable name="V" private_interface="out" public_interface="in" units="mV"/>
<variable name="Ca_i" private_interface="out" public_interface="in" units="mM"/>
<variable name="Ca_o" private_interface="out" public_interface="in" units="mM"/>
<variable name="Ca_SR" private_interface="out" public_interface="in" units="mM"/>
<variable name="time" private_interface="out" public_interface="in" units="ms"/>
<variable name="FVRT_Ca" private_interface="out" public_interface="in" units="dimensionless"/>
<variable name="I_RyR" private_interface="in" public_interface="out" units="mM_per_ms"/>
<variable name="I_LCC" private_interface="in" public_interface="out" units="mM_per_ms"/>
</component>
<component name="CaRU_Transitions">
<variable name="V" public_interface="in" units="mV"/>
<variable name="C_oc" public_interface="in" units="mM"/>
<variable name="C_co" public_interface="in" units="mM"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable initial_value="-2" name="V_L" units="mV"/>
<variable initial_value="7" name="del_VL" units="mV"/>
<variable initial_value="2.35" name="phi_L" units="dimensionless"/>
<variable initial_value="1" name="t_L" units="ms"/>
<variable initial_value="650" name="tau_L" units="ms"/>
<variable name="t_R" units="ms"/>
<variable initial_value="2.43" name="tau_R" units="ms"/>
<variable initial_value="0.05" name="phi_R" units="dimensionless"/>
<variable initial_value="0.012" name="theta_R" units="dimensionless"/>
<variable initial_value="41e-3" name="K_RyR" units="mM"/>
<variable initial_value="0.22e-3" name="K_L" units="mM"/>
<variable initial_value="0.0625" name="a" units="dimensionless"/>
<variable initial_value="14" name="b" units="dimensionless"/>
<variable initial_value="0.01" name="c" units="dimensionless"/>
<variable initial_value="100" name="d" units="dimensionless"/>
<variable name="expVL" units="dimensionless"/>
<variable name="alpha_p" public_interface="out" units="per_ms"/>
<variable name="alpha_m" public_interface="out" units="per_ms"/>
<variable name="beta_poc" public_interface="out" units="per_ms"/>
<variable name="beta_pcc" public_interface="out" units="per_ms"/>
<variable name="beta_m" public_interface="out" units="per_ms"/>
<variable name="epsilon_pco" public_interface="out" units="per_ms"/>
<variable name="epsilon_pcc" public_interface="out" units="per_ms"/>
<variable name="epsilon_m" public_interface="out" units="per_ms"/>
<variable name="mu_poc" public_interface="out" units="per_ms"/>
<variable name="mu_pcc" public_interface="out" units="per_ms"/>
<variable name="mu_moc" public_interface="out" units="per_ms"/>
<variable name="mu_mcc" public_interface="out" units="per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>expVL</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
<ci>V_L</ci>
</apply>
<ci>del_VL</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>t_R</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">1.17</cn>
<ci>t_L</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_p</ci>
<apply>
<divide/>
<ci>expVL</ci>
<apply>
<times/>
<ci>t_L</ci>
<apply>
<plus/>
<ci>expVL</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_m</ci>
<apply>
<divide/>
<ci>phi_L</ci>
<ci>t_L</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_poc</ci>
<apply>
<divide/>
<apply>
<power/>
<ci>C_oc</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>t_R</ci>
<apply>
<plus/>
<apply>
<power/>
<ci>C_oc</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_pcc</ci>
<apply>
<divide/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>t_R</ci>
<apply>
<plus/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_m</ci>
<apply>
<divide/>
<ci>phi_R</ci>
<ci>t_R</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>epsilon_pco</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>C_co</ci>
<apply>
<plus/>
<ci>expVL</ci>
<ci>a</ci>
</apply>
</apply>
<apply>
<times/>
<ci>tau_L</ci>
<ci>K_L</ci>
<apply>
<plus/>
<ci>expVL</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>epsilon_pcc</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>Ca_i</ci>
<apply>
<plus/>
<ci>expVL</ci>
<ci>a</ci>
</apply>
</apply>
<apply>
<times/>
<ci>tau_L</ci>
<ci>K_L</ci>
<apply>
<plus/>
<ci>expVL</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>epsilon_m</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>b</ci>
<apply>
<plus/>
<ci>expVL</ci>
<ci>a</ci>
</apply>
</apply>
<apply>
<times/>
<ci>tau_L</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>b</ci>
<ci>expVL</ci>
</apply>
<ci>a</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_poc</ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<power/>
<ci>C_oc</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>c</ci>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>tau_R</ci>
<apply>
<plus/>
<apply>
<power/>
<ci>C_oc</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_pcc</ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>c</ci>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>tau_R</ci>
<apply>
<plus/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_moc</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>theta_R</ci>
<ci>d</ci>
<apply>
<plus/>
<apply>
<power/>
<ci>C_oc</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>c</ci>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>tau_R</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>d</ci>
<apply>
<power/>
<ci>C_oc</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<ci>c</ci>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_mcc</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>theta_R</ci>
<ci>d</ci>
<apply>
<plus/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>c</ci>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>tau_R</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>d</ci>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<ci>c</ci>
<apply>
<power/>
<ci>K_RyR</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="DS_Calcium_Concentrations">
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable name="Ca_SR" public_interface="in" units="mM"/>
<variable name="Ca_o" public_interface="in" units="mM"/>
<variable name="FVRT_Ca" public_interface="in" units="dimensionless"/>
<variable name="g_D" public_interface="in" units="um3_per_ms"/>
<variable name="J_R" public_interface="in" units="um3_per_ms"/>
<variable name="J_L" public_interface="in" units="um3_per_ms"/>
<variable name="C_cc" public_interface="out" units="mM"/>
<variable name="C_oc" public_interface="out" units="mM"/>
<variable name="C_co" public_interface="out" units="mM"/>
<variable name="C_oo" public_interface="out" units="mM"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>C_cc</ci>
<ci>Ca_i</ci>
</apply>
<apply>
<eq/>
<ci>C_co</ci>
<apply>
<divide/>
<apply>
<plus/>
<ci>Ca_i</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<ci>Ca_SR</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>C_oc</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<plus/>
<ci>Ca_i</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>Ca_o</ci>
<ci>FVRT_Ca</ci>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>FVRT_Ca</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<gt/>
<apply>
<abs/>
<ci>FVRT_Ca</ci>
</apply>
<cn cellml:units="dimensionless">0.000000001</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<plus/>
<ci>Ca_i</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>Ca_o</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>C_oo</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<plus/>
<ci>Ca_i</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<ci>Ca_SR</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>Ca_o</ci>
<ci>FVRT_Ca</ci>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>FVRT_Ca</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<gt/>
<apply>
<abs/>
<ci>FVRT_Ca</ci>
</apply>
<cn cellml:units="dimensionless">0.000000001</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<plus/>
<ci>Ca_i</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<ci>Ca_SR</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>Ca_o</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="LCC_and_RyR_fluxes">
<variable name="g_D" public_interface="in" units="um3_per_ms"/>
<variable name="J_R" public_interface="in" units="um3_per_ms"/>
<variable name="J_L" public_interface="in" units="um3_per_ms"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable name="Ca_SR" public_interface="in" units="mM"/>
<variable name="Ca_o" public_interface="in" units="mM"/>
<variable name="FVRT_Ca" public_interface="in" units="dimensionless"/>
<variable name="J_Loo" public_interface="out" units="um3_mM_per_ms"/>
<variable name="J_Loc" public_interface="out" units="um3_mM_per_ms"/>
<variable name="J_Rco" public_interface="out" units="um3_mM_per_ms"/>
<variable name="J_Roo" public_interface="out" units="um3_mM_per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_Rco</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>J_R</ci>
<apply>
<minus/>
<ci>Ca_SR</ci>
<ci>Ca_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Roo</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<times/>
<ci>J_R</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>Ca_SR</ci>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>FVRT_Ca</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>Ca_SR</ci>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>FVRT_Ca</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<gt/>
<apply>
<abs/>
<ci>FVRT_Ca</ci>
</apply>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<ci>J_R</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>Ca_SR</ci>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>Ca_SR</ci>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>J_Loc</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>J_L</ci>
<ci>FVRT_Ca</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
<ci>Ca_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>FVRT_Ca</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<gt/>
<apply>
<abs/>
<ci>FVRT_Ca</ci>
</apply>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>J_L</ci>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
<ci>Ca_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>J_Loo</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>J_L</ci>
<ci>FVRT_Ca</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<minus/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
<ci>Ca_SR</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<ci>FVRT_Ca</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<ci>FVRT_Ca</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<gt/>
<apply>
<abs/>
<ci>FVRT_Ca</ci>
</apply>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>J_L</ci>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
<ci>Ca_SR</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>J_R</ci>
<ci>g_D</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>J_L</ci>
<ci>g_D</ci>
</apply>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.00001</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="CaRU_states">
<variable name="alpha_p" public_interface="in" units="per_ms"/>
<variable name="alpha_m" public_interface="in" units="per_ms"/>
<variable name="beta_pcc" public_interface="in" units="per_ms"/>
<variable name="beta_poc" public_interface="in" units="per_ms"/>
<variable name="beta_m" public_interface="in" units="per_ms"/>
<variable name="denom" units="per_ms3"/>
<variable name="y_oc" public_interface="out" units="dimensionless"/>
<variable name="y_co" public_interface="out" units="dimensionless"/>
<variable name="y_oo" public_interface="out" units="dimensionless"/>
<variable name="y_cc" public_interface="out" units="dimensionless"/>
<variable name="y_ci" public_interface="out" units="dimensionless"/>
<variable name="y_oi" public_interface="out" units="dimensionless"/>
<variable name="y_ic" public_interface="out" units="dimensionless"/>
<variable name="y_io" public_interface="out" units="dimensionless"/>
<variable name="y_ii" public_interface="out" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>denom</ci>
<apply>
<times/>
<apply>
<plus/>
<ci>alpha_p</ci>
<ci>alpha_m</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<ci>alpha_m</ci>
<ci>beta_m</ci>
<ci>beta_poc</ci>
</apply>
<apply>
<plus/>
<ci>beta_m</ci>
<ci>beta_pcc</ci>
</apply>
</apply>
<apply>
<times/>
<ci>alpha_p</ci>
<apply>
<plus/>
<ci>beta_m</ci>
<ci>beta_poc</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>y_oc</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>alpha_p</ci>
<ci>beta_m</ci>
<apply>
<plus/>
<ci>alpha_p</ci>
<ci>alpha_m</ci>
<ci>beta_m</ci>
<ci>beta_pcc</ci>
</apply>
</apply>
<ci>denom</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>y_co</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>alpha_m</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>beta_pcc</ci>
<apply>
<plus/>
<ci>alpha_m</ci>
<ci>beta_m</ci>
<ci>beta_poc</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_poc</ci>
<ci>alpha_p</ci>
</apply>
</apply>
</apply>
<ci>denom</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>y_oo</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>alpha_p</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>beta_poc</ci>
<apply>
<plus/>
<ci>alpha_p</ci>
<ci>beta_m</ci>
<ci>beta_pcc</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_pcc</ci>
<ci>alpha_m</ci>
</apply>
</apply>
</apply>
<ci>denom</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>y_cc</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>alpha_m</ci>
<ci>beta_m</ci>
<apply>
<plus/>
<ci>alpha_m</ci>
<ci>alpha_p</ci>
<ci>beta_m</ci>
<ci>beta_poc</ci>
</apply>
</apply>
<ci>denom</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>y_ci</ci>
<apply>
<divide/>
<ci>alpha_m</ci>
<apply>
<plus/>
<ci>alpha_p</ci>
<ci>alpha_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>y_oi</ci>
<apply>
<divide/>
<ci>alpha_p</ci>
<apply>
<plus/>
<ci>alpha_p</ci>
<ci>alpha_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>y_ic</ci>
<apply>
<divide/>
<ci>beta_m</ci>
<apply>
<plus/>
<ci>beta_pcc</ci>
<ci>beta_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>y_io</ci>
<apply>
<divide/>
<ci>beta_pcc</ci>
<apply>
<plus/>
<ci>beta_pcc</ci>
<ci>beta_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>y_ii</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>y_oc</ci>
</apply>
<ci>y_co</ci>
</apply>
<ci>y_oo</ci>
</apply>
<ci>y_cc</ci>
</apply>
<ci>y_ci</ci>
</apply>
<ci>y_ic</ci>
</apply>
<ci>y_oi</ci>
</apply>
<ci>y_io</ci>
</apply>
</apply>
</math>
</component>
<component name="CaRU_reduced_states">
<variable name="y_oc" public_interface="in" units="dimensionless"/>
<variable name="y_cc" public_interface="in" units="dimensionless"/>
<variable name="y_co" public_interface="in" units="dimensionless"/>
<variable name="alpha_p" public_interface="in" units="per_ms"/>
<variable name="alpha_m" public_interface="in" units="per_ms"/>
<variable name="beta_pcc" public_interface="in" units="per_ms"/>
<variable name="beta_m" public_interface="in" units="per_ms"/>
<variable name="epsilon_pco" public_interface="in" units="per_ms"/>
<variable name="epsilon_pcc" public_interface="in" units="per_ms"/>
<variable name="epsilon_m" public_interface="in" units="per_ms"/>
<variable name="mu_poc" public_interface="in" units="per_ms"/>
<variable name="mu_pcc" public_interface="in" units="per_ms"/>
<variable name="mu_moc" public_interface="in" units="per_ms"/>
<variable name="mu_mcc" public_interface="in" units="per_ms"/>
<variable name="r_1" public_interface="out" units="per_ms"/>
<variable name="r_2" public_interface="out" units="per_ms"/>
<variable name="r_3" public_interface="out" units="per_ms"/>
<variable name="r_4" public_interface="out" units="per_ms"/>
<variable name="r_5" public_interface="out" units="per_ms"/>
<variable name="r_6" public_interface="out" units="per_ms"/>
<variable name="r_7" public_interface="out" units="per_ms"/>
<variable name="r_8" public_interface="out" units="per_ms"/>
<variable initial_value="0.988494715189347" name="z_1" public_interface="out" units="dimensionless"/>
<variable initial_value="0.00875264158667971" name="z_2" public_interface="out" units="dimensionless"/>
<variable initial_value="0.00272874220458962" name="z_3" public_interface="out" units="dimensionless"/>
<variable name="z_4" public_interface="out" units="dimensionless"/>
<variable name="time" public_interface="in" units="ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>r_1</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>y_oc</ci>
<ci>mu_poc</ci>
</apply>
<apply>
<times/>
<ci>y_cc</ci>
<ci>mu_pcc</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_2</ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_p</ci>
<ci>mu_moc</ci>
</apply>
<apply>
<times/>
<ci>alpha_m</ci>
<ci>mu_mcc</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>alpha_p</ci>
<ci>alpha_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_3</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>beta_m</ci>
<ci>mu_pcc</ci>
</apply>
<apply>
<plus/>
<ci>beta_m</ci>
<ci>beta_pcc</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_4</ci>
<ci>mu_mcc</ci>
</apply>
<apply>
<eq/>
<ci>r_5</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>y_co</ci>
<ci>epsilon_pco</ci>
</apply>
<apply>
<times/>
<ci>y_cc</ci>
<ci>epsilon_pcc</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_6</ci>
<ci>epsilon_m</ci>
</apply>
<apply>
<eq/>
<ci>r_7</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>alpha_m</ci>
<ci>epsilon_pcc</ci>
</apply>
<apply>
<plus/>
<ci>alpha_p</ci>
<ci>alpha_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_8</ci>
<ci>epsilon_m</ci>
</apply>
<apply>
<eq/>
<ci>z_4</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>z_1</ci>
</apply>
<ci>z_2</ci>
</apply>
<ci>z_3</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>z_1</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<plus/>
<ci>r_1</ci>
<ci>r_5</ci>
</apply>
</apply>
<ci>z_1</ci>
</apply>
<apply>
<times/>
<ci>r_2</ci>
<ci>z_2</ci>
</apply>
<apply>
<times/>
<ci>r_6</ci>
<ci>z_3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>z_2</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>r_1</ci>
<ci>z_1</ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci>r_2</ci>
<ci>r_7</ci>
</apply>
<ci>z_2</ci>
</apply>
</apply>
<apply>
<times/>
<ci>r_8</ci>
<ci>z_4</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>z_3</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>r_5</ci>
<ci>z_1</ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci>r_6</ci>
<ci>r_3</ci>
</apply>
<ci>z_3</ci>
</apply>
</apply>
<apply>
<times/>
<ci>r_4</ci>
<ci>z_4</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="RyR_current">
<variable name="y_oo" public_interface="in" units="dimensionless"/>
<variable name="y_co" public_interface="in" units="dimensionless"/>
<variable name="y_oc" public_interface="in" units="dimensionless"/>
<variable name="beta_pcc" public_interface="in" units="per_ms"/>
<variable name="beta_m" public_interface="in" units="per_ms"/>
<variable name="J_Roo" public_interface="in" units="um3_mM_per_ms"/>
<variable name="J_Rco" public_interface="in" units="um3_mM_per_ms"/>
<variable name="z_1" public_interface="in" units="dimensionless"/>
<variable name="z_3" public_interface="in" units="dimensionless"/>
<variable name="V_myo" public_interface="in" units="um3"/>
<variable name="N" public_interface="in" units="dimensionless"/>
<variable name="J_R1" public_interface="out" units="um3_mM_per_ms"/>
<variable name="J_R3" public_interface="out" units="um3_mM_per_ms"/>
<variable name="I_RyR" public_interface="out" units="mM_per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_R1</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>y_oo</ci>
<ci>J_Roo</ci>
</apply>
<apply>
<times/>
<ci>J_Rco</ci>
<ci>y_co</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_R3</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>J_Rco</ci>
<ci>beta_pcc</ci>
</apply>
<apply>
<plus/>
<ci>beta_m</ci>
<ci>beta_pcc</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_RyR</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci>z_1</ci>
<ci>J_R1</ci>
</apply>
<apply>
<times/>
<ci>z_3</ci>
<ci>J_R3</ci>
</apply>
</apply>
<ci>N</ci>
</apply>
<ci>V_myo</ci>
</apply>
</apply>
</math>
</component>
<component name="LCC_current">
<variable name="y_oo" public_interface="in" units="dimensionless"/>
<variable name="y_oc" public_interface="in" units="dimensionless"/>
<variable name="alpha_p" public_interface="in" units="per_ms"/>
<variable name="alpha_m" public_interface="in" units="per_ms"/>
<variable name="J_Loo" public_interface="in" units="um3_mM_per_ms"/>
<variable name="J_Loc" public_interface="in" units="um3_mM_per_ms"/>
<variable name="z_1" public_interface="in" units="dimensionless"/>
<variable name="z_2" public_interface="in" units="dimensionless"/>
<variable name="V_myo" public_interface="in" units="um3"/>
<variable name="N" public_interface="in" units="dimensionless"/>
<variable name="J_L1" public_interface="out" units="um3_mM_per_ms"/>
<variable name="J_L2" public_interface="out" units="um3_mM_per_ms"/>
<variable name="I_LCC" public_interface="out" units="mM_per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_L1</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>J_Loo</ci>
<ci>y_oo</ci>
</apply>
<apply>
<times/>
<ci>J_Loc</ci>
<ci>y_oc</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_L2</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>J_Loc</ci>
<ci>alpha_p</ci>
</apply>
<apply>
<plus/>
<ci>alpha_p</ci>
<ci>alpha_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_LCC</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci>z_1</ci>
<ci>J_L1</ci>
</apply>
<apply>
<times/>
<ci>z_2</ci>
<ci>J_L2</ci>
</apply>
</apply>
<ci>N</ci>
</apply>
<ci>V_myo</ci>
</apply>
</apply>
</math>
</component>
<component name="Na_Ca_Exchanger">
<variable name="FVRT" public_interface="in" units="dimensionless"/>
<variable initial_value="87.5" name="K_mNa" units="mM"/>
<variable initial_value="1.38" name="K_mCa" units="mM"/>
<variable initial_value="0.35" name="eta" units="dimensionless"/>
<variable initial_value="0.1" name="k_sat" units="dimensionless"/>
<variable initial_value="38.5e-3" name="g_NCX" units="mM_per_ms"/>
<variable name="Na_i" public_interface="in" units="mM"/>
<variable name="Na_o" public_interface="in" units="mM"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable name="Ca_o" public_interface="in" units="mM"/>
<variable name="I_NaCa" public_interface="out" units="mM_per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_NaCa</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>g_NCX</ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<exp/>
<apply>
<times/>
<ci>eta</ci>
<ci>FVRT</ci>
</apply>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>Ca_o</ci>
</apply>
<apply>
<times/>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>eta</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>FVRT</ci>
</apply>
</apply>
<apply>
<power/>
<ci>Na_o</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<power/>
<ci>Na_o</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<apply>
<power/>
<ci>K_mNa</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
<apply>
<plus/>
<ci>Ca_o</ci>
<ci>K_mCa</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>k_sat</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>eta</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>FVRT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SERCA">
<variable initial_value="0.45e-3" name="g_SERCA" units="mM_per_ms"/>
<variable initial_value="0.5e-3" name="K_SERCA" units="mM"/>
<variable name="I_SERCA" public_interface="out" units="mM_per_ms"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_SERCA</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>g_SERCA</ci>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>K_SERCA</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Sarcolemmal_Ca_pump">
<variable initial_value="0.0035e-3" name="g_pCa" units="mM_per_ms"/>
<variable initial_value="0.5e-3" name="K_pCa" units="mM"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable name="I_pCa" public_interface="out" units="mM_per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_pCa</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>g_pCa</ci>
<ci>Ca_i</ci>
</apply>
<apply>
<plus/>
<ci>K_pCa</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Background_Ca_current">
<variable name="V" public_interface="in" units="mV"/>
<variable name="R" public_interface="in" units="mJ_per_mole_K"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="F" public_interface="in" units="C_per_mole"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable name="Ca_o" public_interface="in" units="mM"/>
<variable name="E_Ca" public_interface="out" units="mV"/>
<variable initial_value="2.6875e-8" name="g_CaB" units="mM_per_mV_ms"/>
<variable name="I_CaB" public_interface="out" units="mM_per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E_Ca</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Ca_o</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_CaB</ci>
<apply>
<times/>
<ci>g_CaB</ci>
<apply>
<minus/>
<ci>E_Ca</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SR_Ca_leak_current">
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable name="Ca_SR" public_interface="in" units="mM"/>
<variable initial_value="1.8951e-5" name="g_SRl" units="per_ms"/>
<variable name="I_SR" public_interface="out" units="mM_per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_SR</ci>
<apply>
<times/>
<ci>g_SRl</ci>
<apply>
<minus/>
<ci>Ca_SR</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="troponin_Ca_buffer">
<variable initial_value="0.04" name="k_m_TRPN" units="per_ms"/>
<variable initial_value="0.04e3" name="k_p_TRPN" units="per_mM_ms"/>
<variable initial_value="70e-3" name="B_TRPN" units="mM"/>
<variable name="TRPN" public_interface="in" units="mM"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable name="I_TRPN" public_interface="out" units="mM_per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_TRPN</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_m_TRPN</ci>
<apply>
<minus/>
<ci>B_TRPN</ci>
<ci>TRPN</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k_p_TRPN</ci>
<ci>TRPN</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calmodulin_Ca_buffer">
<variable initial_value="2.382e-3" name="k_CMDN" units="mM"/>
<variable initial_value="50e-3" name="B_CMDN" units="mM"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable name="beta_CMDN" public_interface="out" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>beta_CMDN</ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>k_CMDN</ci>
<ci>B_CMDN</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>k_CMDN</ci>
<ci>Ca_i</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fluo4AM_buffer">
<variable initial_value="1" name="k_fluo" units="mM"/>
<!--From Table 3 Thomas2000 found in HeLa cells-->
<variable initial_value="1e-3" name="B_fluo" units="mM"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable name="beta_fluo" public_interface="out" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>beta_fluo</ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>k_fluo</ci>
<ci>B_fluo</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>k_fluo</ci>
<ci>Ca_i</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="intracell">
<variable initial_value="10" name="Na_i" public_interface="out" units="mM"/>
<variable initial_value="0.0001" name="Ca_i" public_interface="out" units="mM"/>
<variable initial_value="0.76" name="Ca_SR" public_interface="out" units="mM"/>
<variable initial_value="0.06" name="TRPN" public_interface="out" units="mM"/>
<variable name="time" public_interface="in" units="ms"/>
<variable name="V_myo" public_interface="in" units="um3"/>
<variable name="V_SR" public_interface="in" units="um3"/>
<variable name="beta_CMDN" public_interface="in" units="dimensionless"/>
<variable name="beta_fluo" public_interface="in" units="dimensionless"/>
<variable name="I_SERCA" public_interface="in" units="mM_per_ms"/>
<variable name="I_pCa" public_interface="in" units="mM_per_ms"/>
<variable name="I_CaB" public_interface="in" units="mM_per_ms"/>
<variable name="I_TRPN" public_interface="in" units="mM_per_ms"/>
<variable name="I_SR" public_interface="in" units="mM_per_ms"/>
<variable name="I_NaCa" public_interface="in" units="mM_per_ms"/>
<variable name="I_RyR" public_interface="in" units="mM_per_ms"/>
<variable name="I_LCC" public_interface="in" units="mM_per_ms"/>
<variable name="CaSR_plot" units="mM"/>
<variable name="open_probability" public_interface="in" units="dimensionless"/>
<variable name="I_IP3R" units="mM_per_ms"/>
<variable initial_value="0.45" name="k_f" units="um3_per_ms"/>
<variable name="I_other" units="mM"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_other</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_LCC</ci>
<ci>I_SR</ci>
<ci>I_NaCa</ci>
</apply>
<ci>I_pCa</ci>
</apply>
<ci>I_CaB</ci>
<ci>I_TRPN</ci>
</apply>
</apply>
</math>
<variable initial_value="1000" name="N_IP3R" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_IP3R</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>k_f</ci>
<ci>N_IP3R</ci>
<ci>open_probability</ci>
<apply>
<minus/>
<ci>Ca_SR</ci>
<ci>Ca_i</ci>
</apply>
</apply>
<ci>V_myo</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>TRPN</ci>
</apply>
<ci>I_TRPN</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<ci>beta_fluo</ci>
<ci>beta_CMDN</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>I_RyR</ci>
<ci>I_SERCA</ci>
</apply>
<ci>I_IP3R</ci>
<ci>I_other</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_SR</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>V_myo</ci>
<ci>V_SR</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>I_RyR</ci>
</apply>
<ci>I_SERCA</ci>
</apply>
<ci>I_SR</ci>
</apply>
<ci>I_IP3R</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CaSR_plot</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>Ca_SR</ci>
<ci>V_SR</ci>
</apply>
<ci>V_myo</ci>
</apply>
</apply>
</math>
</component>
<component name="extracellular_ion_concentrations">
<variable initial_value="140" name="Na_o" public_interface="out" units="mM"/>
<variable initial_value="1" name="Ca_o" public_interface="out" units="mM"/>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="intracell"/>
<component_ref component="extracellular_ion_concentrations"/>
<component_ref component="cell_geometry"/>
<component_ref component="CaRU">
<component_ref component="CaRU_Transitions"/>
<component_ref component="CaRU_states"/>
<component_ref component="CaRU_reduced_states"/>
<component_ref component="DS_Calcium_Concentrations"/>
<component_ref component="LCC_and_RyR_fluxes"/>
<component_ref component="RyR_current"/>
<component_ref component="LCC_current"/>
</component_ref>
<component_ref component="Na_Ca_Exchanger"/>
<component_ref component="SERCA"/>
<component_ref component="Sarcolemmal_Ca_pump"/>
<component_ref component="Background_Ca_current"/>
<component_ref component="SR_Ca_leak_current"/>
<component_ref component="troponin_Ca_buffer"/>
<component_ref component="calmodulin_Ca_buffer"/>
<component_ref component="fluo4AM_buffer"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="CaRU">
<component_ref component="CaRU_Transitions"/>
<component_ref component="CaRU_states"/>
<component_ref component="CaRU_reduced_states"/>
<component_ref component="DS_Calcium_Concentrations"/>
<component_ref component="LCC_and_RyR_fluxes"/>
<component_ref component="RyR_current"/>
<component_ref component="LCC_current"/>
</component_ref>
</group>
<connection>
<map_components component_1="environment" component_2="membrane"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="CaRU"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="intracell"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="cell_geometry" component_2="CaRU"/>
<map_variables variable_1="V_myo" variable_2="V_myo"/>
</connection>
<connection>
<map_components component_1="CaRU" component_2="CaRU_reduced_states"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="CaRU" component_2="LCC_and_RyR_fluxes"/>
<map_variables variable_1="g_D" variable_2="g_D"/>
<map_variables variable_1="J_R" variable_2="J_R"/>
<map_variables variable_1="J_L" variable_2="J_L"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Ca_SR" variable_2="Ca_SR"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
<map_variables variable_1="FVRT_Ca" variable_2="FVRT_Ca"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="CaRU"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="FVRT_Ca" variable_2="FVRT_Ca"/>
</connection>
<connection>
<map_components component_1="CaRU" component_2="CaRU_Transitions"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="CaRU_Transitions" component_2="CaRU_states"/>
<map_variables variable_1="alpha_p" variable_2="alpha_p"/>
<map_variables variable_1="alpha_m" variable_2="alpha_m"/>
<map_variables variable_1="beta_pcc" variable_2="beta_pcc"/>
<map_variables variable_1="beta_poc" variable_2="beta_poc"/>
<map_variables variable_1="beta_m" variable_2="beta_m"/>
</connection>
<connection>
<map_components component_1="CaRU_Transitions" component_2="RyR_current"/>
<map_variables variable_1="beta_pcc" variable_2="beta_pcc"/>
<map_variables variable_1="beta_m" variable_2="beta_m"/>
</connection>
<connection>
<map_components component_1="CaRU_Transitions" component_2="LCC_current"/>
<map_variables variable_1="alpha_p" variable_2="alpha_p"/>
<map_variables variable_1="alpha_m" variable_2="alpha_m"/>
</connection>
<connection>
<map_components component_1="CaRU_states" component_2="RyR_current"/>
<map_variables variable_1="y_oo" variable_2="y_oo"/>
<map_variables variable_1="y_co" variable_2="y_co"/>
<map_variables variable_1="y_oc" variable_2="y_oc"/>
</connection>
<connection>
<map_components component_1="CaRU_states" component_2="LCC_current"/>
<map_variables variable_1="y_oo" variable_2="y_oo"/>
<map_variables variable_1="y_oc" variable_2="y_oc"/>
</connection>
<connection>
<map_components component_1="LCC_and_RyR_fluxes" component_2="RyR_current"/>
<map_variables variable_1="J_Roo" variable_2="J_Roo"/>
<map_variables variable_1="J_Rco" variable_2="J_Rco"/>
</connection>
<connection>
<map_components component_1="LCC_and_RyR_fluxes" component_2="LCC_current"/>
<map_variables variable_1="J_Loo" variable_2="J_Loo"/>
<map_variables variable_1="J_Loc" variable_2="J_Loc"/>
</connection>
<connection>
<map_components component_1="CaRU_reduced_states" component_2="RyR_current"/>
<map_variables variable_1="z_1" variable_2="z_1"/>
<map_variables variable_1="z_3" variable_2="z_3"/>
</connection>
<connection>
<map_components component_1="CaRU_reduced_states" component_2="LCC_current"/>
<map_variables variable_1="z_1" variable_2="z_1"/>
<map_variables variable_1="z_2" variable_2="z_2"/>
</connection>
<connection>
<map_components component_1="DS_Calcium_Concentrations" component_2="CaRU_Transitions"/>
<map_variables variable_1="C_oc" variable_2="C_oc"/>
<map_variables variable_1="C_co" variable_2="C_co"/>
</connection>
<connection>
<map_components component_1="CaRU_Transitions" component_2="CaRU_reduced_states"/>
<map_variables variable_1="alpha_p" variable_2="alpha_p"/>
<map_variables variable_1="alpha_m" variable_2="alpha_m"/>
<map_variables variable_1="beta_pcc" variable_2="beta_pcc"/>
<map_variables variable_1="beta_m" variable_2="beta_m"/>
<map_variables variable_1="epsilon_pco" variable_2="epsilon_pco"/>
<map_variables variable_1="epsilon_pcc" variable_2="epsilon_pcc"/>
<map_variables variable_1="epsilon_m" variable_2="epsilon_m"/>
<map_variables variable_1="mu_poc" variable_2="mu_poc"/>
<map_variables variable_1="mu_pcc" variable_2="mu_pcc"/>
<map_variables variable_1="mu_moc" variable_2="mu_moc"/>
<map_variables variable_1="mu_mcc" variable_2="mu_mcc"/>
</connection>
<connection>
<map_components component_1="CaRU_states" component_2="CaRU_reduced_states"/>
<map_variables variable_1="y_oc" variable_2="y_oc"/>
<map_variables variable_1="y_cc" variable_2="y_cc"/>
<map_variables variable_1="y_co" variable_2="y_co"/>
</connection>
<connection>
<map_components component_1="CaRU" component_2="DS_Calcium_Concentrations"/>
<map_variables variable_1="g_D" variable_2="g_D"/>
<map_variables variable_1="J_R" variable_2="J_R"/>
<map_variables variable_1="J_L" variable_2="J_L"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
<map_variables variable_1="Ca_SR" variable_2="Ca_SR"/>
<map_variables variable_1="FVRT_Ca" variable_2="FVRT_Ca"/>
</connection>
<connection>
<map_components component_1="CaRU" component_2="RyR_current"/>
<map_variables variable_1="V_myo" variable_2="V_myo"/>
<map_variables variable_1="N" variable_2="N"/>
<map_variables variable_1="I_RyR" variable_2="I_RyR"/>
</connection>
<connection>
<map_components component_1="CaRU" component_2="LCC_current"/>
<map_variables variable_1="V_myo" variable_2="V_myo"/>
<map_variables variable_1="N" variable_2="N"/>
<map_variables variable_1="I_LCC" variable_2="I_LCC"/>
</connection>
<connection>
<map_components component_1="intracell" component_2="CaRU"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Ca_SR" variable_2="Ca_SR"/>
<map_variables variable_1="I_RyR" variable_2="I_RyR"/>
<map_variables variable_1="I_LCC" variable_2="I_LCC"/>
</connection>
<connection>
<map_components component_1="extracellular_ion_concentrations" component_2="CaRU"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
</connection>
<connection>
<map_components component_1="intracell" component_2="Na_Ca_Exchanger"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="I_NaCa" variable_2="I_NaCa"/>
</connection>
<connection>
<map_components component_1="extracellular_ion_concentrations" component_2="Na_Ca_Exchanger"/>
<map_variables variable_1="Na_o" variable_2="Na_o"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="Na_Ca_Exchanger"/>
<map_variables variable_1="FVRT" variable_2="FVRT"/>
</connection>
<connection>
<map_components component_1="intracell" component_2="SERCA"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="I_SERCA" variable_2="I_SERCA"/>
</connection>
<connection>
<map_components component_1="intracell" component_2="Sarcolemmal_Ca_pump"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="I_pCa" variable_2="I_pCa"/>
</connection>
<connection>
<map_components component_1="intracell" component_2="Background_Ca_current"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="I_CaB" variable_2="I_CaB"/>
</connection>
<connection>
<map_components component_1="extracellular_ion_concentrations" component_2="Background_Ca_current"/>
<map_variables variable_1="Ca_o" variable_2="Ca_o"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="Background_Ca_current"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="intracell" component_2="SR_Ca_leak_current"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Ca_SR" variable_2="Ca_SR"/>
<map_variables variable_1="I_SR" variable_2="I_SR"/>
</connection>
<connection>
<map_components component_1="intracell" component_2="troponin_Ca_buffer"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="TRPN" variable_2="TRPN"/>
<map_variables variable_1="I_TRPN" variable_2="I_TRPN"/>
</connection>
<connection>
<map_components component_1="intracell" component_2="calmodulin_Ca_buffer"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="beta_CMDN" variable_2="beta_CMDN"/>
</connection>
<connection>
<map_components component_1="intracell" component_2="fluo4AM_buffer"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="beta_fluo" variable_2="beta_fluo"/>
</connection>
<connection>
<map_components component_1="cell_geometry" component_2="intracell"/>
<map_variables variable_1="V_myo" variable_2="V_myo"/>
<map_variables variable_1="V_SR" variable_2="V_SR"/>
</connection>
<!-- *********-->
<component name="open_probability">
<variable name="open_probability" public_interface="out" units="dimensionless"/>
<variable name="beta" public_interface="in" units="dimensionless"/>
<variable name="alpha" public_interface="in" units="dimensionless"/>
<variable initial_value="0.324" name="k_beta" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>open_probability</ci>
<apply>
<divide/>
<ci>beta</ci>
<apply>
<plus/>
<ci>beta</ci>
<apply>
<times/>
<ci>k_beta</ci>
<apply>
<plus/>
<ci>beta</ci>
<ci>alpha</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="alphabeta">
<variable name="alpha" public_interface="out" units="dimensionless"/>
<variable name="beta" public_interface="out" units="dimensionless"/>
<variable name="B" units="dimensionless"/>
<variable name="m" units="dimensionless"/>
<variable initial_value="0.99" name="h" units="dimensionless"/>
<variable name="h_alpha" units="dimensionless"/>
<variable name="Ca_i" public_interface="in" units="mM"/>
<variable initial_value="8e-3" name="K_c" units="mM"/>
<variable initial_value="2.2e-3" name="K_h" units="mM"/>
<variable initial_value="2e-4" name="K_p" units="mM"/>
<variable name="IP_3" units="mM"/>
<variable name="time" public_interface="in" units="ms"/>
<variable initial_value="10" name="t_max" units="per_ms"/>
<variable initial_value="1e-3" name="K_t" units="mM"/>
<variable initial_value="2.7e-5" name="t_p" units="per_ms"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>B</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>m</ci>
<ci>h_alpha</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta</ci>
<apply>
<times/>
<ci>B</ci>
<ci>m</ci>
<ci>h</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>m</ci>
<apply>
<divide/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>K_c</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>h_alpha</ci>
<apply>
<divide/>
<apply>
<power/>
<ci>K_h</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>K_h</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>B</ci>
<apply>
<divide/>
<apply>
<power/>
<ci>IP_3</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>K_p</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>IP_3</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>h</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<ci>h_alpha</ci>
<ci>h</ci>
</apply>
<ci>t_max</ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>K_t</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>t_max</ci>
<apply>
<power/>
<ci>K_t</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>IP_3</ci>
<cn cellml:units="mM" type="e-notation">1<sep/>-2</cn>
</apply>
</math>
<!-- According to Remus2006, this is the amount of IP3 released by ET-1-->
</component>
<connection>
<map_components component_1="open_probability" component_2="alphabeta"/>
<map_variables variable_1="alpha" variable_2="alpha"/>
<map_variables variable_1="beta" variable_2="beta"/>
</connection>
<connection>
<map_components component_1="open_probability" component_2="intracell"/>
<map_variables variable_1="open_probability" variable_2="open_probability"/>
</connection>
<connection>
<map_components component_1="alphabeta" component_2="intracell"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="alphabeta" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
</model>