- Author:
- James Lawson <j.lawson@auckland.ac.nz>
- Date:
- 2009-08-27 16:21:39+12:00
- Desc:
- added 1.1 model
- Permanent Source URI:
- http://models.cellml.org/workspace/terkildsen_niederer_crampin_hunter_smith_2008/rawfile/6149bed591084b9ad24439be12b5aa8fe4dbf387/Hinch_et_al_2004.cellml
<?xml version="1.0" encoding="utf-8"?>
<model name="Hinch_et_al_2004"
cmeta:id="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 units="second" prefix="milli" exponent="-1"/>
</units>
<units name="ms">
<unit units="second" prefix="milli"/>
</units>
<units name="mV">
<unit units="volt" prefix="milli"/>
</units>
<units name="mM">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="mM_per_ms">
<unit units="mM"/>
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<units name="per_mM_ms">
<unit units="mM" exponent="-1"/>
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<units name="mM_per_mV_ms">
<unit units="mM"/>
<unit units="second" prefix="milli" exponent="-1"/>
<unit units="volt" prefix="milli" exponent="-1"/>
</units>
<units name="per_ms3">
<unit units="second" prefix="milli" exponent="-3"/>
</units>
<units name="um2">
<unit units="meter" prefix="micro" exponent="2"/>
</units>
<units name="um3">
<unit units="meter" prefix="micro" exponent="3"/>
</units>
<units name="um3_per_ms">
<unit units="meter" prefix="micro" exponent="3"/>
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<units name="um3_mM_per_ms">
<unit units="meter" prefix="micro" exponent="3"/>
<unit units="mM"/>
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<units name="mJ_per_mole_K">
<unit units="joule" prefix="milli"/>
<unit units="mole" exponent="-1"/>
<unit units="kelvin" exponent="-1"/>
</units>
<units name="C_per_mole">
<unit units="coulomb"/>
<unit units="mole" exponent="-1"/>
</units>
<component name="environment">
<variable name="time" units="ms" public_interface="out"/>
</component>
<component name="cell_geometry">
<variable name="V_myo" units="um3" initial_value="25.84e3" public_interface="out"/>
<variable name="V_SR" units="um3" initial_value="2.098e3" public_interface="out"/>
<variable name="A_cap" units="um2" initial_value="1.534e4" public_interface="out"/>
</component>
<component name="membrane">
<variable name="V" units="mV" public_interface="out"/>
<variable name="R" units="mJ_per_mole_K" initial_value="8314.5" public_interface="out"/>
<variable name="T" units="kelvin" initial_value="295" public_interface="out"/>
<variable name="F" units="C_per_mole" initial_value="96487" public_interface="out"/>
<variable name="time" units="ms" public_interface="in"/>
<variable name="FVRT" units="dimensionless" public_interface="out"/>
<variable name="FVRT_Ca" units="dimensionless" public_interface="out"/>
<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>
<cn cellml:units="mV">0</cn>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<cn cellml:units="ms">0</cn>
</apply>
<apply>
<leq/>
<ci>time</ci>
<cn cellml:units="ms">200</cn>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<minus/>
<cn cellml:units="mV">80</cn>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="CaRU">
<variable name="g_D" units="um3_per_ms" initial_value="0.065" private_interface="out"/>
<variable name="J_R" units="um3_per_ms" initial_value="0.02" private_interface="out"/>
<variable name="J_L" units="um3_per_ms" initial_value="9.13e-4" private_interface="out"/>
<variable name="V_myo" units="um3" public_interface="in" private_interface="out"/>
<variable name="N" units="dimensionless" initial_value="50000" private_interface="out"/>
<variable name="V" units="mV" public_interface="in" private_interface="out"/>
<variable name="Ca_i" units="mM" public_interface="in" private_interface="out"/>
<variable name="Ca_o" units="mM" public_interface="in" private_interface="out"/>
<variable name="Ca_SR" units="mM" public_interface="in" private_interface="out"/>
<variable name="time" units="ms" public_interface="in" private_interface="out"/>
<variable name="FVRT_Ca" units="dimensionless" public_interface="in" private_interface="out"/>
<variable name="I_RyR" units="mM_per_ms" public_interface="out" private_interface="in"/>
<variable name="I_LCC" units="mM_per_ms" public_interface="out" private_interface="in"/>
</component>
<component name="CaRU_Transitions">
<variable name="V" units="mV" public_interface="in"/>
<variable name="C_oc" units="mM" public_interface="in"/>
<variable name="C_co" units="mM" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="V_L" units="mV" initial_value="-2"/>
<variable name="del_VL" units="mV" initial_value="7"/>
<variable name="phi_L" units="dimensionless" initial_value="2.35"/>
<variable name="t_L" units="ms" initial_value="1"/>
<variable name="tau_L" units="ms" initial_value="650"/>
<variable name="t_R" units="ms"/>
<variable name="tau_R" units="ms" initial_value="2.43"/>
<variable name="phi_R" units="dimensionless" initial_value="0.05"/>
<variable name="theta_R" units="dimensionless" initial_value="0.012"/>
<variable name="K_RyR" units="mM" initial_value="41e-3"/>
<variable name="K_L" units="mM" initial_value="0.22e-3"/>
<variable name="a" units="dimensionless" initial_value="0.0625"/>
<variable name="b" units="dimensionless" initial_value="14"/>
<variable name="c" units="dimensionless" initial_value="0.01"/>
<variable name="d" units="dimensionless" initial_value="100"/>
<variable name="expVL" units="dimensionless"/>
<variable name="alpha_p" units="per_ms" public_interface="out"/>
<variable name="alpha_m" units="per_ms" public_interface="out"/>
<variable name="beta_poc" units="per_ms" public_interface="out"/>
<variable name="beta_pcc" units="per_ms" public_interface="out"/>
<variable name="beta_m" units="per_ms" public_interface="out"/>
<variable name="epsilon_pco" units="per_ms" public_interface="out"/>
<variable name="epsilon_pcc" units="per_ms" public_interface="out"/>
<variable name="epsilon_m" units="per_ms" public_interface="out"/>
<variable name="mu_poc" units="per_ms" public_interface="out"/>
<variable name="mu_pcc" units="per_ms" public_interface="out"/>
<variable name="mu_moc" units="per_ms" public_interface="out"/>
<variable name="mu_mcc" units="per_ms" public_interface="out"/>
<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" units="mM" public_interface="in"/>
<variable name="Ca_SR" units="mM" public_interface="in"/>
<variable name="Ca_o" units="mM" public_interface="in"/>
<variable name="FVRT_Ca" units="dimensionless" public_interface="in"/>
<variable name="g_D" units="um3_per_ms" public_interface="in"/>
<variable name="J_R" units="um3_per_ms" public_interface="in"/>
<variable name="J_L" units="um3_per_ms" public_interface="in"/>
<variable name="C_cc" units="mM" public_interface="out"/>
<variable name="C_oc" units="mM" public_interface="out"/>
<variable name="C_co" units="mM" public_interface="out"/>
<variable name="C_oo" units="mM" public_interface="out"/>
<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" units="um3_per_ms" public_interface="in"/>
<variable name="J_R" units="um3_per_ms" public_interface="in"/>
<variable name="J_L" units="um3_per_ms" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="Ca_SR" units="mM" public_interface="in"/>
<variable name="Ca_o" units="mM" public_interface="in"/>
<variable name="FVRT_Ca" units="dimensionless" public_interface="in"/>
<variable name="J_Loo" units="um3_mM_per_ms" public_interface="out"/>
<variable name="J_Loc" units="um3_mM_per_ms" public_interface="out"/>
<variable name="J_Rco" units="um3_mM_per_ms" public_interface="out"/>
<variable name="J_Roo" units="um3_mM_per_ms" public_interface="out"/>
<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" units="per_ms" public_interface="in"/>
<variable name="alpha_m" units="per_ms" public_interface="in"/>
<variable name="beta_pcc" units="per_ms" public_interface="in"/>
<variable name="beta_poc" units="per_ms" public_interface="in"/>
<variable name="beta_m" units="per_ms" public_interface="in"/>
<variable name="denom" units="per_ms3"/>
<variable name="y_oc" units="dimensionless" public_interface="out"/>
<variable name="y_co" units="dimensionless" public_interface="out"/>
<variable name="y_oo" units="dimensionless" public_interface="out"/>
<variable name="y_cc" units="dimensionless" public_interface="out"/>
<variable name="y_ci" units="dimensionless" public_interface="out"/>
<variable name="y_oi" units="dimensionless" public_interface="out"/>
<variable name="y_ic" units="dimensionless" public_interface="out"/>
<variable name="y_io" units="dimensionless" public_interface="out"/>
<variable name="y_ii" units="dimensionless" public_interface="out"/>
<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" units="dimensionless" public_interface="in"/>
<variable name="y_cc" units="dimensionless" public_interface="in"/>
<variable name="y_co" units="dimensionless" public_interface="in"/>
<variable name="alpha_p" units="per_ms" public_interface="in"/>
<variable name="alpha_m" units="per_ms" public_interface="in"/>
<variable name="beta_pcc" units="per_ms" public_interface="in"/>
<variable name="beta_m" units="per_ms" public_interface="in"/>
<variable name="epsilon_pco" units="per_ms" public_interface="in"/>
<variable name="epsilon_pcc" units="per_ms" public_interface="in"/>
<variable name="epsilon_m" units="per_ms" public_interface="in"/>
<variable name="mu_poc" units="per_ms" public_interface="in"/>
<variable name="mu_pcc" units="per_ms" public_interface="in"/>
<variable name="mu_moc" units="per_ms" public_interface="in"/>
<variable name="mu_mcc" units="per_ms" public_interface="in"/>
<variable name="r_1" units="per_ms" public_interface="out"/>
<variable name="r_2" units="per_ms" public_interface="out"/>
<variable name="r_3" units="per_ms" public_interface="out"/>
<variable name="r_4" units="per_ms" public_interface="out"/>
<variable name="r_5" units="per_ms" public_interface="out"/>
<variable name="r_6" units="per_ms" public_interface="out"/>
<variable name="r_7" units="per_ms" public_interface="out"/>
<variable name="r_8" units="per_ms" public_interface="out"/>
<variable name="z_1" units="dimensionless" initial_value="0.98859" public_interface="out"/>
<variable name="z_2" units="dimensionless" initial_value="0.0087302" public_interface="out"/>
<variable name="z_3" units="dimensionless" initial_value="0.0026566" public_interface="out"/>
<variable name="z_4" units="dimensionless" public_interface="out"/>
<variable name="time" units="ms" public_interface="in"/>
<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" units="dimensionless" public_interface="in"/>
<variable name="y_co" units="dimensionless" public_interface="in"/>
<variable name="y_oc" units="dimensionless" public_interface="in"/>
<variable name="beta_pcc" units="per_ms" public_interface="in"/>
<variable name="beta_m" units="per_ms" public_interface="in"/>
<variable name="J_Roo" units="um3_mM_per_ms" public_interface="in"/>
<variable name="J_Rco" units="um3_mM_per_ms" public_interface="in"/>
<variable name="z_1" units="dimensionless" public_interface="in"/>
<variable name="z_3" units="dimensionless" public_interface="in"/>
<variable name="V_myo" units="um3" public_interface="in"/>
<variable name="N" units="dimensionless" public_interface="in"/>
<variable name="J_R1" units="um3_mM_per_ms" public_interface="out"/>
<variable name="J_R3" units="um3_mM_per_ms" public_interface="out"/>
<variable name="I_RyR" units="mM_per_ms" public_interface="out"/>
<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" units="dimensionless" public_interface="in"/>
<variable name="y_oc" units="dimensionless" public_interface="in"/>
<variable name="alpha_p" units="per_ms" public_interface="in"/>
<variable name="alpha_m" units="per_ms" public_interface="in"/>
<variable name="J_Loo" units="um3_mM_per_ms" public_interface="in"/>
<variable name="J_Loc" units="um3_mM_per_ms" public_interface="in"/>
<variable name="z_1" units="dimensionless" public_interface="in"/>
<variable name="z_2" units="dimensionless" public_interface="in"/>
<variable name="V_myo" units="um3" public_interface="in"/>
<variable name="N" units="dimensionless" public_interface="in"/>
<variable name="J_L1" units="um3_mM_per_ms" public_interface="out"/>
<variable name="J_L2" units="um3_mM_per_ms" public_interface="out"/>
<variable name="I_LCC" units="mM_per_ms" public_interface="out"/>
<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" units="dimensionless" public_interface="in"/>
<variable name="K_mNa" units="mM" initial_value="87.5"/>
<variable name="K_mCa" units="mM" initial_value="1.38"/>
<variable name="eta" units="dimensionless" initial_value="0.35"/>
<variable name="k_sat" units="dimensionless" initial_value="0.1"/>
<variable name="g_NCX" units="mM_per_ms" initial_value="38.5e-3"/>
<variable name="Na_i" units="mM" public_interface="in"/>
<variable name="Na_o" units="mM" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="Ca_o" units="mM" public_interface="in"/>
<variable name="I_NaCa" units="mM_per_ms" public_interface="out"/>
<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 name="g_SERCA" units="mM_per_ms" initial_value="0.45e-3"/>
<variable name="K_SERCA" units="mM" initial_value="0.5e-3"/>
<variable name="I_SERCA" units="mM_per_ms" public_interface="out"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<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 name="g_pCa" units="mM_per_ms" initial_value="0.0035e-3"/>
<variable name="K_mpCa" units="mM" initial_value="0.5e-3"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="I_pCa" units="mM_per_ms" public_interface="out"/>
<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_mpCa</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Background_Ca_current">
<variable name="V" units="mV" public_interface="in"/>
<variable name="R" units="mJ_per_mole_K" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="F" units="C_per_mole" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="Ca_o" units="mM" public_interface="in"/>
<variable name="E_Ca" units="mV" public_interface="out"/>
<variable name="g_CaB" units="mM_per_mV_ms" initial_value="2.6875e-8"/>
<variable name="I_CaB" units="mM_per_ms" public_interface="out"/>
<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" units="mM" public_interface="in"/>
<variable name="Ca_SR" units="mM" public_interface="in"/>
<variable name="g_SRl" units="per_ms" initial_value="1.8951e-5"/>
<variable name="I_SR" units="mM_per_ms" public_interface="out"/>
<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 name="k_m_TRPN" units="per_ms" initial_value="0.04"/>
<variable name="k_p_TRPN" units="per_mM_ms" initial_value="0.04e3"/>
<variable name="B_TRPN" units="mM" initial_value="70e-3"/>
<variable name="TRPN" units="mM" public_interface="in"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="I_TRPN" units="mM_per_ms" public_interface="out"/>
<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 name="k_CMDN" units="mM" initial_value="2.382e-3"/>
<variable name="B_CMDN" units="mM" initial_value="50e-3"/>
<variable name="Ca_i" units="mM" public_interface="in"/>
<variable name="beta_CMDN" units="dimensionless" public_interface="out"/>
<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="intracellular_ion_concentrations">
<variable name="Na_i" units="mM" initial_value="10" public_interface="out"/>
<variable name="Ca_i" units="mM" initial_value="0.0001" public_interface="out"/>
<variable name="Ca_SR" units="mM" initial_value="700e-3" public_interface="out"/>
<variable name="TRPN" units="mM" initial_value="0.0636364" public_interface="out"/>
<variable name="time" units="ms" public_interface="in"/>
<variable name="V_myo" units="um3" public_interface="in"/>
<variable name="V_SR" units="um3" public_interface="in"/>
<variable name="beta_CMDN" units="dimensionless" public_interface="in"/>
<variable name="I_SERCA" units="mM_per_ms" public_interface="in"/>
<variable name="I_pCa" units="mM_per_ms" public_interface="in"/>
<variable name="I_CaB" units="mM_per_ms" public_interface="in"/>
<variable name="I_TRPN" units="mM_per_ms" public_interface="in"/>
<variable name="I_SR" units="mM_per_ms" public_interface="in"/>
<variable name="I_NaCa" units="mM_per_ms" public_interface="in"/>
<variable name="I_RyR" units="mM_per_ms" public_interface="in"/>
<variable name="I_LCC" units="mM_per_ms" public_interface="in"/>
<variable name="CaSR_plot" units="mM"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<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_CMDN</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_LCC</ci>
<ci>I_RyR</ci>
</apply>
<ci>I_SERCA</ci>
</apply>
<ci>I_SR</ci>
<ci>I_NaCa</ci>
</apply>
<ci>I_pCa</ci>
</apply>
<ci>I_CaB</ci>
<ci>I_TRPN</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>
<plus/>
<apply>
<minus/>
<ci>I_RyR</ci>
</apply>
<ci>I_SERCA</ci>
</apply>
<ci>I_SR</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 name="Na_o" units="mM" initial_value="140" public_interface="out"/>
<variable name="Ca_o" units="mM" initial_value="1" public_interface="out"/>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="intracellular_ion_concentrations"/>
<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>
</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="intracellular_ion_concentrations"/>
<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="intracellular_ion_concentrations" 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="intracellular_ion_concentrations" 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="intracellular_ion_concentrations" 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="intracellular_ion_concentrations" 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="intracellular_ion_concentrations" 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="intracellular_ion_concentrations" 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="intracellular_ion_concentrations" 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="intracellular_ion_concentrations" 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="cell_geometry" component_2="intracellular_ion_concentrations"/>
<map_variables variable_1="V_myo" variable_2="V_myo"/>
<map_variables variable_1="V_SR" variable_2="V_SR"/>
</connection>
</model>
