- Author:
- Catherine Lloyd <c.lloyd@auckland.ac.nz>
- Date:
- 2009-11-10 13:25:34+13:00
- Desc:
- Added image documentation ready for Hanne to create the diagram.
- Permanent Source URI:
- http://models.cellml.org/workspace/noble_noble_2001/rawfile/ec66390f5bcdcd51cb22afa68a9f17cf20185c54/noble_noble_2001.cellml
<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 02/10/2007 at 17:26:58 using:
COR (0.9.31.793)
Copyright 2002-2007 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk
CellML 1.0 was used to generate this cellular model
http://www.CellML.org/
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="noble_noble_2001" name="noble_noble_2001">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Noble Ventricular Cell Model 1998</title>
<author>
<firstname>Penny</firstname>
<surname>Noble</surname>
<affiliation>
<shortaffil>Oxford University</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This version of the model has been curated and unit checked and is known to run in PCEnv, COR and Jsim. A PCEnv session file is also available for this model.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
The Noble98 (Noble et al. 1998) model of the guinea-pig ventricular cell was a development of the Noble91 (Noble et al. 1991) model and included representation of the components of IK, (IKr and IKs), the persistent sodium current and a dyadic space between the sarcolemma and the sarcoplasmic reticulum. In several important respects this model was a major improvement on the earlier model. However on rigorous testing for new applications problems associated with the calcium dynamics equations were found in the new model.
The model demonstrated here is an attempt at 'correcting' the inability of Noble98 to reproduce a mechanism of initiation of Early After-Depolarizations (EADs), the voltage-dependence of calcium release from the sarcoplasmic reticulum (SR) and also to improve the shape of the calcium transient.
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="noble_2001.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram of the cell model.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<units name="micrometre">
<unit units="metre" prefix="micro"/>
</units>
<units name="per_second">
<unit units="second" exponent="-1"/>
</units>
<units name="millivolt">
<unit units="volt" prefix="milli"/>
</units>
<units name="per_millivolt">
<unit units="volt" prefix="milli" exponent="-1"/>
</units>
<units name="per_millivolt_second">
<unit units="millivolt" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="microS">
<unit units="siemens" prefix="micro"/>
</units>
<units name="microF">
<unit units="farad" prefix="micro"/>
</units>
<units name="nanoA">
<unit units="ampere" prefix="nano"/>
</units>
<units name="mA_nA">
<unit units="ampere" prefix="milli"/>
<unit units="nanoA" exponent="-1"/>
</units>
<units name="nanoA_per_millimolar">
<unit units="ampere" prefix="nano"/>
<unit units="millimolar" exponent="-1"/>
</units>
<units name="millimolar">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="millimolar4">
<unit units="millimolar" exponent="4"/>
</units>
<units name="millimolar_per_second">
<unit units="millimolar"/>
<unit units="second" exponent="-1"/>
</units>
<units name="per_millimolar_second">
<unit units="millimolar" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="joule_per_kilomole_kelvin">
<unit units="joule"/>
<unit units="mole" prefix="kilo" exponent="-1"/>
<unit units="kelvin" exponent="-1"/>
</units>
<units name="coulomb_per_mole">
<unit units="coulomb"/>
<unit units="mole" exponent="-1"/>
</units>
<units name="micrometre3">
<unit units="metre" prefix="micro" exponent="3"/>
</units>
<units name="litre_micrometre3">
<unit units="litre"/>
<unit units="micrometre3" exponent="-1"/>
</units>
<component name="environment">
<variable units="second" public_interface="out" name="time"/>
</component>
<component name="membrane">
<variable units="millivolt" public_interface="out" name="V" initial_value="-92.8915042"/>
<variable units="joule_per_kilomole_kelvin" public_interface="out" name="R" initial_value="8314.472"/>
<variable units="kelvin" public_interface="out" name="T" initial_value="310"/>
<variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96485.3415"/>
<variable units="microF" name="Cm" initial_value="9.5e-5"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="nanoA" public_interface="in" name="i_K1"/>
<variable units="nanoA" public_interface="in" name="i_to"/>
<variable units="nanoA" public_interface="in" name="i_Kr"/>
<variable units="nanoA" public_interface="in" name="i_Ks"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_K_cyt"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_K_ds"/>
<variable units="nanoA" public_interface="in" name="i_NaK"/>
<variable units="nanoA" public_interface="in" name="i_Na"/>
<variable units="nanoA" public_interface="in" name="i_b_Na"/>
<variable units="nanoA" public_interface="in" name="i_p_Na"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_Na_cyt"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_Na_ds"/>
<variable units="nanoA" public_interface="in" name="i_NaCa_cyt"/>
<variable units="nanoA" public_interface="in" name="i_NaCa_ds"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_Ca_cyt"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_Ca_ds"/>
<variable units="nanoA" public_interface="in" name="i_b_Ca"/>
<variable units="nanoA" name="i_Stim"/>
<variable units="second" name="stim_start" initial_value="0.1"/>
<variable units="second" name="stim_end" initial_value="10000"/>
<variable units="second" name="stim_period" initial_value="1"/>
<variable units="second" name="stim_duration" initial_value="0.002"/>
<variable units="nanoA" name="stim_amplitude" initial_value="-6"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Stim</ci>
<piecewise>
<piece>
<ci>stim_amplitude</ci>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<apply>
<leq/>
<ci>time</ci>
<ci>stim_end</ci>
</apply>
<apply>
<leq/>
<apply>
<minus/>
<apply>
<minus/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<apply>
<times/>
<apply>
<floor/>
<apply>
<divide/>
<apply>
<minus/>
<ci>time</ci>
<ci>stim_start</ci>
</apply>
<ci>stim_period</ci>
</apply>
</apply>
<ci>stim_period</ci>
</apply>
</apply>
<ci>stim_duration</ci>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="nanoA">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>V</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>Cm</ci>
</apply>
<apply>
<plus/>
<ci>i_Stim</ci>
<ci>i_K1</ci>
<ci>i_to</ci>
<ci>i_Kr</ci>
<ci>i_Ks</ci>
<ci>i_NaK</ci>
<ci>i_Na</ci>
<ci>i_b_Na</ci>
<ci>i_p_Na</ci>
<ci>i_Ca_L_Na_cyt</ci>
<ci>i_Ca_L_Na_ds</ci>
<ci>i_NaCa_cyt</ci>
<ci>i_NaCa_ds</ci>
<ci>i_Ca_L_Ca_cyt</ci>
<ci>i_Ca_L_Ca_ds</ci>
<ci>i_Ca_L_K_cyt</ci>
<ci>i_Ca_L_K_ds</ci>
<ci>i_b_Ca</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="reversal_potentials">
<variable units="millivolt" public_interface="out" name="E_Na"/>
<variable units="millivolt" public_interface="out" name="E_K"/>
<variable units="millivolt" public_interface="out" name="E_Ks"/>
<variable units="millivolt" public_interface="out" name="E_Ca"/>
<variable units="millivolt" public_interface="out" name="E_mh"/>
<variable units="dimensionless" name="P_kna" initial_value="0.03"/>
<variable units="millimolar" public_interface="in" name="K_o"/>
<variable units="millimolar" public_interface="in" name="Na_o"/>
<variable units="millimolar" public_interface="in" name="K_i"/>
<variable units="millimolar" public_interface="in" name="Na_i"/>
<variable units="millimolar" public_interface="in" name="Ca_o"/>
<variable units="millimolar" public_interface="in" name="Ca_i"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Na_o</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_o</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_Ks</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<apply>
<plus/>
<ci>K_o</ci>
<apply>
<times/>
<ci>P_kna</ci>
<ci>Na_o</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>K_i</ci>
<apply>
<times/>
<ci>P_kna</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_Ca</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Ca_o</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_mh</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<apply>
<plus/>
<ci>Na_o</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.12</cn>
<ci>K_o</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>Na_i</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.12</cn>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_independent_potassium_current">
<variable units="nanoA" public_interface="out" name="i_K1"/>
<variable units="millimolar" public_interface="out" name="K_mk1" initial_value="10"/>
<variable units="microS" name="g_K1" initial_value="0.5"/>
<variable units="millimolar" public_interface="in" name="K_o"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_K1</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>g_K1</ci>
<ci>K_o</ci>
</apply>
<apply>
<plus/>
<ci>K_o</ci>
<ci>K_mk1</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">10</cn>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">1.25</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapid_delayed_rectifier_potassium_current">
<variable units="nanoA" public_interface="out" name="i_Kr"/>
<variable units="microS" name="g_Kr1" initial_value="0.0028"/>
<variable units="microS" name="g_Kr2" initial_value="0.0017"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="dimensionless" private_interface="in" name="xr1"/>
<variable units="dimensionless" private_interface="in" name="xr2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Kr</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci>g_Kr1</ci>
<ci>xr1</ci>
</apply>
<apply>
<times/>
<ci>g_Kr2</ci>
<ci>xr2</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">9</cn>
</apply>
<cn cellml:units="millivolt">22.4</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapid_delayed_rectifier_potassium_current_xr1_gate">
<variable units="dimensionless" public_interface="out" name="xr1" initial_value="1.02e-5"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="per_second" name="alpha_xr1"/>
<variable units="per_second" name="beta_xr1"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_xr1</ci>
<apply>
<divide/>
<cn cellml:units="per_second">50</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
<cn cellml:units="millivolt">9</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_xr1</ci>
<apply>
<times/>
<cn cellml:units="per_second">0.05</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">20</cn>
</apply>
</apply>
<cn cellml:units="millivolt">15</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>xr1</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_xr1</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>xr1</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_xr1</ci>
<ci>xr1</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapid_delayed_rectifier_potassium_current_xr2_gate">
<variable units="dimensionless" public_interface="out" name="xr2" initial_value="2e-7"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="per_second" name="alpha_xr2"/>
<variable units="per_second" name="beta_xr2"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_xr2</ci>
<apply>
<divide/>
<cn cellml:units="per_second">50</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
<cn cellml:units="millivolt">9</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_xr2</ci>
<apply>
<times/>
<cn cellml:units="per_second">0.4</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">30</cn>
</apply>
<cn cellml:units="millivolt">30</cn>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>xr2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_xr2</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>xr2</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_xr2</ci>
<ci>xr2</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_delayed_rectifier_potassium_current">
<variable units="nanoA" public_interface="out" name="i_Ks"/>
<variable units="microS" name="g_Ks" initial_value="0.0016"/>
<variable units="millivolt" public_interface="in" name="E_Ks"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="dimensionless" private_interface="in" name="xs"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ks</ci>
<apply>
<times/>
<ci>g_Ks</ci>
<apply>
<power/>
<ci>xs</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Ks</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_delayed_rectifier_potassium_current_xs_gate">
<variable units="dimensionless" public_interface="out" name="xs" initial_value="0.0006469"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="per_second" name="alpha_xs"/>
<variable units="per_second" name="beta_xs"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_xs</ci>
<apply>
<divide/>
<cn cellml:units="per_second">14</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
<cn cellml:units="millivolt">9</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_xs</ci>
<apply>
<times/>
<cn cellml:units="per_second">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">45</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>xs</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_xs</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>xs</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_xs</ci>
<ci>xs</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current">
<variable units="nanoA" public_interface="out" name="i_Na"/>
<variable units="microS" name="g_Na" initial_value="0.5"/>
<variable units="millivolt" public_interface="in" name="E_mh"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="dimensionless" private_interface="in" name="m"/>
<variable units="dimensionless" private_interface="in" name="h"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Na</ci>
<apply>
<times/>
<ci>g_Na</ci>
<apply>
<power/>
<ci>m</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>h</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_mh</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current_m_gate">
<variable units="dimensionless" public_interface="out" name="m" initial_value="0.0016111"/>
<variable units="per_second" name="alpha_m"/>
<variable units="per_second" name="beta_m"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millivolt" name="delta_m" initial_value="1e-5"/>
<variable units="millivolt" name="E0_m"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E0_m</ci>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">41</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_m</ci>
<piecewise>
<piece>
<cn cellml:units="per_second">2000</cn>
<apply>
<lt/>
<apply>
<abs/>
<ci>E0_m</ci>
</apply>
<ci>delta_m</ci>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">200</cn>
<ci>E0_m</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.1</cn>
</apply>
<ci>E0_m</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_m</ci>
<apply>
<times/>
<cn cellml:units="per_second">8000</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.056</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">66</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>m</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_m</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>m</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_m</ci>
<ci>m</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current_h_gate">
<variable units="dimensionless" public_interface="out" name="h" initial_value="0.9944559"/>
<variable units="per_second" name="alpha_h"/>
<variable units="per_second" name="beta_h"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_h</ci>
<apply>
<times/>
<cn cellml:units="per_second">20</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.125</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_h</ci>
<apply>
<divide/>
<cn cellml:units="per_second">2000</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">320</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>h</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_h</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>h</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_h</ci>
<ci>h</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="persistent_sodium_current">
<variable units="nanoA" public_interface="out" name="i_p_Na"/>
<variable units="microS" name="g_pna" initial_value="0.0027"/>
<variable units="millivolt" public_interface="in" name="E_Na"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_p_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>g_pna</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">52</cn>
</apply>
</apply>
<cn cellml:units="millivolt">8</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_background_current">
<variable units="nanoA" public_interface="out" name="i_b_Na"/>
<variable units="microS" name="g_bna" initial_value="0.0006"/>
<variable units="millimolar" public_interface="in" name="Na_i"/>
<variable units="millivolt" public_interface="in" name="E_Na"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_b_Na</ci>
<apply>
<times/>
<ci>g_bna</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel">
<variable units="nanoA" public_interface="out" name="i_Ca_L_Ca_cyt"/>
<variable units="nanoA" public_interface="out" name="i_Ca_L_K_cyt"/>
<variable units="nanoA" public_interface="out" name="i_Ca_L_Na_cyt"/>
<variable units="nanoA" public_interface="out" name="i_Ca_L_Ca_ds"/>
<variable units="nanoA" public_interface="out" name="i_Ca_L_K_ds"/>
<variable units="nanoA" public_interface="out" name="i_Ca_L_Na_ds"/>
<variable units="nanoA" public_interface="out" name="i_Ca_L"/>
<variable units="nanoA_per_millimolar" name="P_Ca_L" initial_value="1"/>
<variable units="dimensionless" name="P_CaK" initial_value="0.002"/>
<variable units="dimensionless" name="P_CaNa" initial_value="0.01"/>
<variable units="millimolar" public_interface="in" name="Ca_o"/>
<variable units="millimolar" public_interface="in" private_interface="out" name="Ca_i"/>
<variable units="millimolar" public_interface="in" private_interface="out" name="Ca_ds"/>
<variable units="millimolar" public_interface="in" name="K_o"/>
<variable units="millimolar" public_interface="in" name="K_i"/>
<variable units="millimolar" public_interface="in" name="Na_o"/>
<variable units="millimolar" public_interface="in" name="Na_i"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="dimensionless" private_interface="in" name="d"/>
<variable units="dimensionless" private_interface="in" name="CaChoncyt"/>
<variable units="dimensionless" private_interface="in" name="CaChonds"/>
<variable units="millimolar" private_interface="out" name="KCaChoff" initial_value="0.01"/>
<variable units="millimolar" private_interface="out" name="Kmdsinact" initial_value="0.001"/>
<variable units="dimensionless" name="FrICa" initial_value="0.7"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ca_L_Ca_cyt</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>FrICa</ci>
</apply>
<cn cellml:units="dimensionless">4</cn>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>CaChoncyt</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">100</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_K_cyt</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>FrICa</ci>
</apply>
<ci>P_CaK</ci>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>CaChoncyt</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>K_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">50</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>K_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_Na_cyt</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>FrICa</ci>
</apply>
<ci>P_CaNa</ci>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>CaChoncyt</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Na_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">50</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Na_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_Ca_ds</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>FrICa</ci>
<cn cellml:units="dimensionless">4</cn>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>CaChonds</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Ca_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">100</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_K_ds</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>FrICa</ci>
<ci>P_CaK</ci>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>CaChonds</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>K_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">50</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>K_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L_Na_ds</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>FrICa</ci>
<ci>P_CaNa</ci>
<ci>P_Ca_L</ci>
<ci>d</ci>
<ci>CaChonds</ci>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>Na_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="millivolt">50</cn>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Na_o</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_Ca_L</ci>
<apply>
<plus/>
<ci>i_Ca_L_Ca_cyt</ci>
<ci>i_Ca_L_K_cyt</ci>
<ci>i_Ca_L_Na_cyt</ci>
<ci>i_Ca_L_Ca_ds</ci>
<ci>i_Ca_L_K_ds</ci>
<ci>i_Ca_L_Na_ds</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_d_gate">
<variable units="dimensionless" public_interface="out" name="d" initial_value="0"/>
<variable units="per_second" name="alpha_d"/>
<variable units="per_second" name="beta_d"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millivolt" name="E0_d"/>
<variable units="dimensionless" name="speed_d" initial_value="3"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E0_d</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">24</cn>
</apply>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_d</ci>
<piecewise>
<piece>
<cn cellml:units="per_second">120</cn>
<apply>
<lt/>
<apply>
<abs/>
<ci>E0_d</ci>
</apply>
<cn cellml:units="millivolt">0.0001</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">30</cn>
<ci>E0_d</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>E0_d</ci>
</apply>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_d</ci>
<piecewise>
<piece>
<cn cellml:units="per_second">120</cn>
<apply>
<lt/>
<apply>
<abs/>
<ci>E0_d</ci>
</apply>
<cn cellml:units="millivolt">0.0001</cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">12</cn>
<ci>E0_d</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>E0_d</ci>
<cn cellml:units="millivolt">10</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>d</ci>
</apply>
<apply>
<times/>
<ci>speed_d</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_d</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>d</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_d</ci>
<ci>d</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_f_gate">
<variable units="dimensionless" name="f" initial_value="0"/>
<variable units="dimensionless" public_interface="out" name="CaChoncyt"/>
<variable units="dimensionless" public_interface="out" name="CaChonds"/>
<variable units="dimensionless" name="CaChoffcyt"/>
<variable units="dimensionless" name="CaChoffds"/>
<variable units="per_second" name="alpha_f"/>
<variable units="per_second" name="beta_f"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millimolar" public_interface="in" name="Ca_i"/>
<variable units="millimolar" public_interface="in" name="Ca_ds"/>
<variable units="millimolar" public_interface="in" name="KCaChoff"/>
<variable units="millimolar" public_interface="in" name="Kmdsinact"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="dimensionless" name="speed_f" initial_value="0.3"/>
<variable units="millivolt" name="delta_f" initial_value="0.0001"/>
<variable units="millivolt" name="E0_f"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E0_f</ci>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">34</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_f</ci>
<piecewise>
<piece>
<cn cellml:units="per_second">25</cn>
<apply>
<lt/>
<apply>
<abs/>
<ci>E0_f</ci>
</apply>
<ci>delta_f</ci>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">6.25</cn>
<ci>E0_f</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>E0_f</ci>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_f</ci>
<apply>
<divide/>
<cn cellml:units="per_second">12</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">34</cn>
</apply>
</apply>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CaChoffcyt</ci>
<apply>
<divide/>
<ci>Ca_i</ci>
<apply>
<plus/>
<ci>KCaChoff</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>f</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">120</cn>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>f</ci>
</apply>
<ci>CaChoffcyt</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>f</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>CaChoffcyt</ci>
</apply>
</apply>
</apply>
<ci>speed_f</ci>
<ci>beta_f</ci>
</apply>
<apply>
<times/>
<ci>alpha_f</ci>
<ci>f</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CaChoncyt</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>f</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>CaChoffcyt</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CaChoffds</ci>
<apply>
<divide/>
<ci>Ca_ds</ci>
<apply>
<plus/>
<ci>Kmdsinact</ci>
<ci>Ca_ds</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CaChonds</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>f</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>CaChoffds</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_background_current">
<variable units="nanoA" public_interface="out" name="i_b_Ca"/>
<variable units="microS" name="g_bca" initial_value="0.00025"/>
<variable units="millivolt" public_interface="in" name="E_Ca"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_b_Ca</ci>
<apply>
<times/>
<ci>g_bca</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Ca</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="transient_outward_current">
<variable units="nanoA" public_interface="out" name="i_to"/>
<variable units="microS" name="g_to" initial_value="0.005"/>
<variable units="dimensionless" name="g_tos" initial_value="0"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="dimensionless" private_interface="in" name="s"/>
<variable units="dimensionless" private_interface="in" name="r"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_to</ci>
<apply>
<times/>
<ci>g_to</ci>
<apply>
<plus/>
<ci>g_tos</ci>
<apply>
<times/>
<ci>s</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>g_tos</ci>
</apply>
</apply>
</apply>
<ci>r</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="transient_outward_current_s_gate">
<variable units="dimensionless" public_interface="out" name="s" initial_value="0.9948645"/>
<variable units="per_second" name="alpha_s"/>
<variable units="per_second" name="beta_s"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_s</ci>
<apply>
<times/>
<cn cellml:units="per_second">0.033</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">17</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_s</ci>
<apply>
<divide/>
<cn cellml:units="per_second">33</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.125</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">10</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_s</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>s</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_s</ci>
<ci>s</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="transient_outward_current_r_gate">
<variable units="dimensionless" public_interface="out" name="r" initial_value="0"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r</ci>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">333</cn>
<apply>
<minus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
</apply>
</apply>
<ci>r</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_potassium_pump">
<variable units="nanoA" public_interface="out" name="i_NaK"/>
<variable units="nanoA" name="i_NaK_max" initial_value="0.7"/>
<variable units="millimolar" name="K_mK" initial_value="1"/>
<variable units="millimolar" name="K_mNa" initial_value="40"/>
<variable units="millimolar" public_interface="in" name="K_o"/>
<variable units="millimolar" public_interface="in" name="Na_i"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_NaK</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>i_NaK_max</ci>
<ci>K_o</ci>
</apply>
<apply>
<plus/>
<ci>K_mK</ci>
<ci>K_o</ci>
</apply>
</apply>
<ci>Na_i</ci>
</apply>
<apply>
<plus/>
<ci>K_mNa</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_calcium_exchanger">
<variable units="nanoA" public_interface="out" name="i_NaCa"/>
<variable units="nanoA" public_interface="out" name="i_NaCa_cyt"/>
<variable units="nanoA" public_interface="out" name="i_NaCa_ds"/>
<variable units="nanoA" name="k_NaCa" initial_value="0.0005"/>
<variable units="dimensionless" name="n_NaCa" initial_value="3"/>
<variable units="dimensionless" name="d_NaCa" initial_value="0"/>
<variable units="dimensionless" name="gamma" initial_value="0.2"/>
<variable units="dimensionless" name="FRiNaCa" initial_value="0.001"/>
<variable units="millimolar" public_interface="in" name="Na_i"/>
<variable units="millimolar" public_interface="in" name="Na_o"/>
<variable units="millimolar" public_interface="in" name="Ca_i"/>
<variable units="millimolar" public_interface="in" name="Ca_ds"/>
<variable units="millimolar" public_interface="in" name="Ca_o"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_NaCa_cyt</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>FRiNaCa</ci>
</apply>
<ci>k_NaCa</ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>gamma</ci>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<ci>n_NaCa</ci>
</apply>
<ci>Ca_o</ci>
</apply>
<apply>
<times/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>gamma</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci>Na_o</ci>
<ci>n_NaCa</ci>
</apply>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="millimolar4">1</cn>
<apply>
<times/>
<ci>d_NaCa</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>Ca_i</ci>
<apply>
<power/>
<ci>Na_o</ci>
<ci>n_NaCa</ci>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<power/>
<ci>Na_i</ci>
<ci>n_NaCa</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ca_i</ci>
<cn cellml:units="millimolar">0.0069</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_NaCa_ds</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>FRiNaCa</ci>
<ci>k_NaCa</ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>gamma</ci>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<ci>n_NaCa</ci>
</apply>
<ci>Ca_o</ci>
</apply>
<apply>
<times/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>gamma</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
<apply>
<minus/>
<ci>n_NaCa</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci>Na_o</ci>
<ci>n_NaCa</ci>
</apply>
<ci>Ca_ds</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="millimolar4">1</cn>
<apply>
<times/>
<ci>d_NaCa</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>Ca_ds</ci>
<apply>
<power/>
<ci>Na_o</ci>
<ci>n_NaCa</ci>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<power/>
<ci>Na_i</ci>
<ci>n_NaCa</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ca_ds</ci>
<cn cellml:units="millimolar">0.0069</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_NaCa</ci>
<apply>
<plus/>
<ci>i_NaCa_cyt</ci>
<ci>i_NaCa_ds</ci>
</apply>
</apply>
</math>
</component>
<component name="sarcoplasmic_reticulum_calcium_pump">
<variable units="millimolar_per_second" public_interface="out" name="i_up"/>
<variable units="dimensionless" name="K_1"/>
<variable units="millimolar" name="K_2"/>
<variable units="millimolar" name="K_cyca" initial_value="0.0003"/>
<variable units="dimensionless" name="K_xcs" initial_value="0.4"/>
<variable units="millimolar" name="K_srca" initial_value="0.5"/>
<variable units="millimolar_per_second" name="alpha_up" initial_value="0.4"/>
<variable units="millimolar_per_second" name="beta_up" initial_value="0.03"/>
<variable units="millimolar" public_interface="in" name="Ca_i"/>
<variable units="millimolar" public_interface="in" name="Ca_up"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>K_1</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>K_cyca</ci>
<ci>K_xcs</ci>
</apply>
<ci>K_srca</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>K_2</ci>
<apply>
<plus/>
<ci>Ca_i</ci>
<apply>
<times/>
<ci>Ca_up</ci>
<ci>K_1</ci>
</apply>
<apply>
<times/>
<ci>K_cyca</ci>
<ci>K_xcs</ci>
</apply>
<ci>K_cyca</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_up</ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<ci>Ca_i</ci>
<ci>K_2</ci>
</apply>
<ci>alpha_up</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>Ca_up</ci>
<ci>K_1</ci>
</apply>
<ci>K_2</ci>
</apply>
<ci>beta_up</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_translocation">
<variable units="millimolar_per_second" public_interface="out" name="i_trans"/>
<variable units="millimolar" public_interface="in" name="Ca_rel"/>
<variable units="millimolar" public_interface="in" name="Ca_up"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_trans</ci>
<apply>
<times/>
<cn cellml:units="per_second">50</cn>
<apply>
<minus/>
<ci>Ca_up</ci>
<ci>Ca_rel</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_release">
<variable units="millimolar_per_second" public_interface="out" name="i_rel"/>
<variable units="dimensionless" name="VoltDep"/>
<variable units="dimensionless" name="RegBindSite"/>
<variable units="dimensionless" name="CaiReg"/>
<variable units="dimensionless" name="CadsReg"/>
<variable units="per_second" name="ActRate"/>
<variable units="per_second" name="InactRate"/>
<variable units="per_second" name="K_leak_rate" initial_value="0.05"/>
<variable units="per_second" name="K_m_rel" initial_value="250"/>
<variable units="millimolar" name="K_m_Ca_cyt" initial_value="0.0005"/>
<variable units="millimolar" name="K_m_Ca_ds" initial_value="0.01"/>
<variable units="dimensionless" name="PrecFrac"/>
<variable units="dimensionless" name="ActFrac" initial_value="0.0049039"/>
<variable units="dimensionless" name="ProdFrac" initial_value="0.6950649"/>
<variable units="dimensionless" name="SpeedRel"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millimolar" public_interface="in" name="Ca_i"/>
<variable units="millimolar" public_interface="in" name="Ca_ds"/>
<variable units="millimolar" public_interface="in" name="Ca_rel"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>VoltDep</ci>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.08</cn>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CaiReg</ci>
<apply>
<divide/>
<ci>Ca_i</ci>
<apply>
<plus/>
<ci>Ca_i</ci>
<ci>K_m_Ca_cyt</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>CadsReg</ci>
<apply>
<divide/>
<ci>Ca_ds</ci>
<apply>
<plus/>
<ci>Ca_ds</ci>
<ci>K_m_Ca_ds</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>RegBindSite</ci>
<apply>
<plus/>
<ci>CaiReg</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>CaiReg</ci>
</apply>
<ci>CadsReg</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>ActRate</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="per_second">0</cn>
<ci>VoltDep</ci>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">90</cn>
<apply>
<power/>
<ci>RegBindSite</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>InactRate</ci>
<apply>
<plus/>
<cn cellml:units="per_second">60</cn>
<apply>
<times/>
<cn cellml:units="per_second">500</cn>
<apply>
<power/>
<ci>RegBindSite</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>SpeedRel</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">5</cn>
<apply>
<lt/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>PrecFrac</ci>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>ActFrac</ci>
</apply>
<ci>ProdFrac</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>ActFrac</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>PrecFrac</ci>
<ci>SpeedRel</ci>
<ci>ActRate</ci>
</apply>
<apply>
<times/>
<ci>ActFrac</ci>
<ci>SpeedRel</ci>
<ci>InactRate</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>ProdFrac</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>ActFrac</ci>
<ci>SpeedRel</ci>
<ci>InactRate</ci>
</apply>
<apply>
<times/>
<ci>SpeedRel</ci>
<cn cellml:units="per_second">1</cn>
<ci>ProdFrac</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_rel</ci>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci>ActFrac</ci>
<apply>
<plus/>
<ci>ActFrac</ci>
<cn cellml:units="dimensionless">0.25</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>K_m_rel</ci>
</apply>
<ci>K_leak_rate</ci>
</apply>
<ci>Ca_rel</ci>
</apply>
</apply>
</math>
</component>
<component name="extracellular_sodium_concentration">
<variable units="millimolar" public_interface="out" name="Na_o" initial_value="140"/>
</component>
<component name="intracellular_sodium_concentration">
<variable units="millimolar" public_interface="out" name="Na_i" initial_value="7.6713487"/>
<variable units="micrometre3" public_interface="in" name="V_i"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="nanoA" public_interface="in" name="i_Na"/>
<variable units="nanoA" public_interface="in" name="i_b_Na"/>
<variable units="nanoA" public_interface="in" name="i_p_Na"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_Na_cyt"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_Na_ds"/>
<variable units="nanoA" public_interface="in" name="i_NaK"/>
<variable units="nanoA" public_interface="in" name="i_NaCa_cyt"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Na_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="mA_nA">1</cn>
</apply>
<apply>
<times/>
<cn cellml:units="litre_micrometre3">1</cn>
<ci>V_i</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>i_Na</ci>
<ci>i_p_Na</ci>
<ci>i_b_Na</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>i_NaK</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>i_NaCa_cyt</ci>
</apply>
<ci>i_Ca_L_Na_cyt</ci>
<ci>i_Ca_L_Na_ds</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="extracellular_calcium_concentration">
<variable units="millimolar" public_interface="out" name="Ca_o" initial_value="2"/>
</component>
<component name="extracellular_potassium_concentration">
<variable units="millimolar" public_interface="out" name="K_o" initial_value="4"/>
</component>
<component name="intracellular_potassium_concentration">
<variable units="millimolar" public_interface="out" name="K_i" initial_value="136.1745362"/>
<variable units="micrometre3" public_interface="in" name="V_i"/>
<variable units="nanoA" public_interface="in" name="i_K1"/>
<variable units="nanoA" public_interface="in" name="i_Kr"/>
<variable units="nanoA" public_interface="in" name="i_Ks"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_K_cyt"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_K_ds"/>
<variable units="nanoA" public_interface="in" name="i_to"/>
<variable units="nanoA" public_interface="in" name="i_NaK"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="mA_nA">1</cn>
</apply>
<apply>
<times/>
<cn cellml:units="litre_micrometre3">1</cn>
<ci>V_i</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_K1</ci>
<ci>i_Kr</ci>
<ci>i_Ks</ci>
<ci>i_Ca_L_K_cyt</ci>
<ci>i_Ca_L_K_ds</ci>
<ci>i_to</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>i_NaK</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="intracellular_calcium_concentration">
<variable units="millimolar" public_interface="out" name="Ca_i" initial_value="1.49e-5"/>
<variable units="millimolar" public_interface="out" name="Ca_ds" initial_value="2.7e-6"/>
<variable units="millimolar" public_interface="out" name="Ca_up" initial_value="0.36963"/>
<variable units="millimolar" public_interface="out" name="Ca_rel" initial_value="0.6460487"/>
<variable units="millimolar" name="Ca_Calmod" initial_value="0.0005841"/>
<variable units="millimolar" name="Ca_Trop" initial_value="0.0003732"/>
<variable units="millimolar" name="Calmod" initial_value="0.02"/>
<variable units="millimolar" name="Trop" initial_value="0.05"/>
<variable units="per_millimolar_second" name="alpha_Calmod" initial_value="100000"/>
<variable units="per_second" name="beta_Calmod" initial_value="50"/>
<variable units="per_millimolar_second" name="alpha_Trop" initial_value="100000"/>
<variable units="per_second" name="beta_Trop" initial_value="200"/>
<variable units="micrometre" name="radius" initial_value="12"/>
<variable units="micrometre" name="length" initial_value="74"/>
<variable units="micrometre3" public_interface="out" name="V_i"/>
<variable units="micrometre3" name="V_Cell"/>
<variable units="dimensionless" name="V_i_ratio"/>
<variable units="dimensionless" name="V_ds_ratio" initial_value="0.1"/>
<variable units="dimensionless" name="V_rel_ratio" initial_value="0.1"/>
<variable units="dimensionless" name="V_e_ratio" initial_value="0.4"/>
<variable units="dimensionless" name="V_up_ratio" initial_value="0.01"/>
<variable units="per_second" name="Kdecay" initial_value="10"/>
<variable units="millimolar_per_second" public_interface="in" name="i_up"/>
<variable units="millimolar_per_second" public_interface="in" name="i_trans"/>
<variable units="millimolar_per_second" public_interface="in" name="i_rel"/>
<variable units="nanoA" public_interface="in" name="i_NaCa_cyt"/>
<variable units="nanoA" public_interface="in" name="i_NaCa_ds"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_Ca_cyt"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L_Ca_ds"/>
<variable units="nanoA" public_interface="in" name="i_b_Ca"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_Cell</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">3.141592654</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>radius</ci>
<cn cellml:units="dimensionless">1000</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>length</ci>
</apply>
<cn cellml:units="dimensionless">1000</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>V_i_ratio</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>V_e_ratio</ci>
</apply>
<ci>V_up_ratio</ci>
</apply>
<ci>V_rel_ratio</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_i</ci>
<apply>
<times/>
<ci>V_Cell</ci>
<ci>V_i_ratio</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_i</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="mA_nA">1</cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<cn cellml:units="litre_micrometre3">1</cn>
<ci>V_i</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_Ca_L_Ca_cyt</ci>
<ci>i_b_Ca</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>i_NaCa_cyt</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>i_NaCa_ds</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_ds</ci>
<ci>V_ds_ratio</ci>
<ci>Kdecay</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>i_rel</ci>
<ci>V_rel_ratio</ci>
</apply>
<ci>V_i_ratio</ci>
</apply>
</apply>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Calmod</ci>
</apply>
</apply>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Trop</ci>
</apply>
</apply>
<ci>i_up</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_ds</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="mA_nA">1</cn>
</apply>
<ci>i_Ca_L_Ca_ds</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<cn cellml:units="litre_micrometre3">1</cn>
<ci>V_ds_ratio</ci>
<ci>V_i</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_ds</ci>
<ci>Kdecay</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_up</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<ci>V_i_ratio</ci>
<ci>V_up_ratio</ci>
</apply>
<ci>i_up</ci>
</apply>
<ci>i_trans</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_rel</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<ci>V_up_ratio</ci>
<ci>V_rel_ratio</ci>
</apply>
<ci>i_trans</ci>
</apply>
<ci>i_rel</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Calmod</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_Calmod</ci>
<ci>Ca_i</ci>
<apply>
<minus/>
<ci>Calmod</ci>
<ci>Ca_Calmod</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_Calmod</ci>
<ci>Ca_Calmod</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Trop</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_Trop</ci>
<ci>Ca_i</ci>
<apply>
<minus/>
<ci>Trop</ci>
<ci>Ca_Trop</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_Trop</ci>
<ci>Ca_Trop</ci>
</apply>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="fast_sodium_current">
<component_ref component="fast_sodium_current_m_gate"/>
<component_ref component="fast_sodium_current_h_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_d_gate"/>
<component_ref component="L_type_Ca_channel_f_gate"/>
</component_ref>
<component_ref component="rapid_delayed_rectifier_potassium_current">
<component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
<component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
</component_ref>
<component_ref component="slow_delayed_rectifier_potassium_current">
<component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
</component_ref>
<component_ref component="time_independent_potassium_current"/>
<component_ref component="persistent_sodium_current"/>
<component_ref component="sodium_background_current"/>
<component_ref component="calcium_background_current"/>
<component_ref component="transient_outward_current">
<component_ref component="transient_outward_current_s_gate"/>
<component_ref component="transient_outward_current_r_gate"/>
</component_ref>
<component_ref component="sodium_potassium_pump"/>
<component_ref component="sodium_calcium_exchanger"/>
<component_ref component="sarcoplasmic_reticulum_calcium_pump"/>
<component_ref component="calcium_translocation"/>
<component_ref component="calcium_release"/>
<component_ref component="extracellular_sodium_concentration"/>
<component_ref component="intracellular_sodium_concentration"/>
<component_ref component="extracellular_potassium_concentration"/>
<component_ref component="intracellular_potassium_concentration"/>
<component_ref component="extracellular_calcium_concentration"/>
<component_ref component="intracellular_calcium_concentration"/>
<component_ref component="reversal_potentials"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="fast_sodium_current">
<component_ref component="fast_sodium_current_m_gate"/>
<component_ref component="fast_sodium_current_h_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_d_gate"/>
<component_ref component="L_type_Ca_channel_f_gate"/>
</component_ref>
<component_ref component="transient_outward_current">
<component_ref component="transient_outward_current_s_gate"/>
<component_ref component="transient_outward_current_r_gate"/>
</component_ref>
<component_ref component="rapid_delayed_rectifier_potassium_current">
<component_ref component="rapid_delayed_rectifier_potassium_current_xr1_gate"/>
<component_ref component="rapid_delayed_rectifier_potassium_current_xr2_gate"/>
</component_ref>
<component_ref component="slow_delayed_rectifier_potassium_current">
<component_ref component="slow_delayed_rectifier_potassium_current_xs_gate"/>
</component_ref>
</group>
<connection>
<map_components component_2="environment" component_1="membrane"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="reversal_potentials"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="fast_sodium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="L_type_Ca_channel"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="rapid_delayed_rectifier_potassium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="slow_delayed_rectifier_potassium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="time_independent_potassium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="persistent_sodium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sodium_background_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="calcium_background_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="transient_outward_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sodium_potassium_pump"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sodium_calcium_exchanger"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="intracellular_sodium_concentration"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="intracellular_potassium_concentration"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="intracellular_calcium_concentration"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sarcoplasmic_reticulum_calcium_pump"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="calcium_translocation"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="calcium_release"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="calcium_release"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="reversal_potentials" component_1="membrane"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Na" variable_1="i_Na"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Ca_L_Ca_cyt" variable_1="i_Ca_L_Ca_cyt"/>
<map_variables variable_2="i_Ca_L_Ca_ds" variable_1="i_Ca_L_Ca_ds"/>
<map_variables variable_2="i_Ca_L_Na_cyt" variable_1="i_Ca_L_Na_cyt"/>
<map_variables variable_2="i_Ca_L_Na_ds" variable_1="i_Ca_L_Na_ds"/>
<map_variables variable_2="i_Ca_L_K_cyt" variable_1="i_Ca_L_K_cyt"/>
<map_variables variable_2="i_Ca_L_K_ds" variable_1="i_Ca_L_K_ds"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="rapid_delayed_rectifier_potassium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Kr" variable_1="i_Kr"/>
</connection>
<connection>
<map_components component_2="slow_delayed_rectifier_potassium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Ks" variable_1="i_Ks"/>
</connection>
<connection>
<map_components component_2="transient_outward_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_to" variable_1="i_to"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_K1" variable_1="i_K1"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="persistent_sodium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_p_Na" variable_1="i_p_Na"/>
</connection>
<connection>
<map_components component_2="sodium_background_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
</connection>
<connection>
<map_components component_2="calcium_background_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_b_Ca" variable_1="i_b_Ca"/>
</connection>
<connection>
<map_components component_2="sodium_potassium_pump" component_1="membrane"/>
<map_variables variable_2="i_NaK" variable_1="i_NaK"/>
</connection>
<connection>
<map_components component_2="sodium_calcium_exchanger" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_NaCa_cyt" variable_1="i_NaCa_cyt"/>
<map_variables variable_2="i_NaCa_ds" variable_1="i_NaCa_ds"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="intracellular_sodium_concentration" component_1="membrane"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="intracellular_potassium_concentration" component_1="membrane"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="membrane"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="sodium_background_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_Na" variable_1="E_Na"/>
</connection>
<connection>
<map_components component_2="persistent_sodium_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_Na" variable_1="E_Na"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="rapid_delayed_rectifier_potassium_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="slow_delayed_rectifier_potassium_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_Ks" variable_1="E_Ks"/>
</connection>
<connection>
<map_components component_2="transient_outward_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="calcium_background_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_Ca" variable_1="E_Ca"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_mh" variable_1="E_mh"/>
</connection>
<connection>
<map_components component_2="extracellular_potassium_concentration" component_1="reversal_potentials"/>
<map_variables variable_2="K_o" variable_1="K_o"/>
</connection>
<connection>
<map_components component_2="extracellular_sodium_concentration" component_1="reversal_potentials"/>
<map_variables variable_2="Na_o" variable_1="Na_o"/>
</connection>
<connection>
<map_components component_2="intracellular_potassium_concentration" component_1="reversal_potentials"/>
<map_variables variable_2="K_i" variable_1="K_i"/>
</connection>
<connection>
<map_components component_2="intracellular_sodium_concentration" component_1="reversal_potentials"/>
<map_variables variable_2="Na_i" variable_1="Na_i"/>
</connection>
<connection>
<map_components component_2="extracellular_calcium_concentration" component_1="reversal_potentials"/>
<map_variables variable_2="Ca_o" variable_1="Ca_o"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="reversal_potentials"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
</connection>
<connection>
<map_components component_2="extracellular_potassium_concentration" component_1="time_independent_potassium_current"/>
<map_variables variable_2="K_o" variable_1="K_o"/>
</connection>
<connection>
<map_components component_2="intracellular_potassium_concentration" component_1="time_independent_potassium_current"/>
<map_variables variable_2="i_K1" variable_1="i_K1"/>
</connection>
<connection>
<map_components component_2="intracellular_potassium_concentration" component_1="rapid_delayed_rectifier_potassium_current"/>
<map_variables variable_2="i_Kr" variable_1="i_Kr"/>
</connection>
<connection>
<map_components component_2="intracellular_potassium_concentration" component_1="slow_delayed_rectifier_potassium_current"/>
<map_variables variable_2="i_Ks" variable_1="i_Ks"/>
</connection>
<connection>
<map_components component_2="intracellular_sodium_concentration" component_1="fast_sodium_current"/>
<map_variables variable_2="i_Na" variable_1="i_Na"/>
</connection>
<connection>
<map_components component_2="intracellular_sodium_concentration" component_1="sodium_background_current"/>
<map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
<map_variables variable_2="Na_i" variable_1="Na_i"/>
</connection>
<connection>
<map_components component_2="intracellular_sodium_concentration" component_1="persistent_sodium_current"/>
<map_variables variable_2="i_p_Na" variable_1="i_p_Na"/>
</connection>
<connection>
<map_components component_2="intracellular_sodium_concentration" component_1="L_type_Ca_channel"/>
<map_variables variable_2="i_Ca_L_Na_cyt" variable_1="i_Ca_L_Na_cyt"/>
<map_variables variable_2="i_Ca_L_Na_ds" variable_1="i_Ca_L_Na_ds"/>
<map_variables variable_2="Na_i" variable_1="Na_i"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="L_type_Ca_channel"/>
<map_variables variable_2="i_Ca_L_Ca_cyt" variable_1="i_Ca_L_Ca_cyt"/>
<map_variables variable_2="i_Ca_L_Ca_ds" variable_1="i_Ca_L_Ca_ds"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
<map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
</connection>
<connection>
<map_components component_2="intracellular_potassium_concentration" component_1="L_type_Ca_channel"/>
<map_variables variable_2="i_Ca_L_K_cyt" variable_1="i_Ca_L_K_cyt"/>
<map_variables variable_2="i_Ca_L_K_ds" variable_1="i_Ca_L_K_ds"/>
<map_variables variable_2="K_i" variable_1="K_i"/>
</connection>
<connection>
<map_components component_2="extracellular_potassium_concentration" component_1="L_type_Ca_channel"/>
<map_variables variable_2="K_o" variable_1="K_o"/>
</connection>
<connection>
<map_components component_2="extracellular_sodium_concentration" component_1="L_type_Ca_channel"/>
<map_variables variable_2="Na_o" variable_1="Na_o"/>
</connection>
<connection>
<map_components component_2="extracellular_calcium_concentration" component_1="L_type_Ca_channel"/>
<map_variables variable_2="Ca_o" variable_1="Ca_o"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="calcium_background_current"/>
<map_variables variable_2="i_b_Ca" variable_1="i_b_Ca"/>
</connection>
<connection>
<map_components component_2="intracellular_sodium_concentration" component_1="sodium_potassium_pump"/>
<map_variables variable_2="i_NaK" variable_1="i_NaK"/>
<map_variables variable_2="Na_i" variable_1="Na_i"/>
</connection>
<connection>
<map_components component_2="extracellular_potassium_concentration" component_1="sodium_potassium_pump"/>
<map_variables variable_2="K_o" variable_1="K_o"/>
</connection>
<connection>
<map_components component_2="intracellular_potassium_concentration" component_1="sodium_potassium_pump"/>
<map_variables variable_2="i_NaK" variable_1="i_NaK"/>
</connection>
<connection>
<map_components component_2="intracellular_potassium_concentration" component_1="transient_outward_current"/>
<map_variables variable_2="i_to" variable_1="i_to"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="sodium_calcium_exchanger"/>
<map_variables variable_2="i_NaCa_cyt" variable_1="i_NaCa_cyt"/>
<map_variables variable_2="i_NaCa_ds" variable_1="i_NaCa_ds"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
<map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
</connection>
<connection>
<map_components component_2="intracellular_sodium_concentration" component_1="sodium_calcium_exchanger"/>
<map_variables variable_2="i_NaCa_cyt" variable_1="i_NaCa_cyt"/>
<map_variables variable_2="Na_i" variable_1="Na_i"/>
</connection>
<connection>
<map_components component_2="extracellular_sodium_concentration" component_1="sodium_calcium_exchanger"/>
<map_variables variable_2="Na_o" variable_1="Na_o"/>
</connection>
<connection>
<map_components component_2="extracellular_calcium_concentration" component_1="sodium_calcium_exchanger"/>
<map_variables variable_2="Ca_o" variable_1="Ca_o"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="sarcoplasmic_reticulum_calcium_pump"/>
<map_variables variable_2="i_up" variable_1="i_up"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
<map_variables variable_2="Ca_up" variable_1="Ca_up"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="calcium_translocation"/>
<map_variables variable_2="i_trans" variable_1="i_trans"/>
<map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
<map_variables variable_2="Ca_up" variable_1="Ca_up"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="calcium_release"/>
<map_variables variable_2="i_rel" variable_1="i_rel"/>
<map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
<map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
</connection>
<connection>
<map_components component_2="rapid_delayed_rectifier_potassium_current_xr1_gate" component_1="rapid_delayed_rectifier_potassium_current"/>
<map_variables variable_2="xr1" variable_1="xr1"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="rapid_delayed_rectifier_potassium_current_xr2_gate" component_1="rapid_delayed_rectifier_potassium_current"/>
<map_variables variable_2="xr2" variable_1="xr2"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="slow_delayed_rectifier_potassium_current_xs_gate" component_1="slow_delayed_rectifier_potassium_current"/>
<map_variables variable_2="xs" variable_1="xs"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current_m_gate" component_1="fast_sodium_current"/>
<map_variables variable_2="m" variable_1="m"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current_h_gate" component_1="fast_sodium_current"/>
<map_variables variable_2="h" variable_1="h"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_d_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="d" variable_1="d"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_f_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="CaChoncyt" variable_1="CaChoncyt"/>
<map_variables variable_2="CaChonds" variable_1="CaChonds"/>
<map_variables variable_2="KCaChoff" variable_1="KCaChoff"/>
<map_variables variable_2="Kmdsinact" variable_1="Kmdsinact"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
<map_variables variable_2="Ca_ds" variable_1="Ca_ds"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="transient_outward_current_s_gate" component_1="transient_outward_current"/>
<map_variables variable_2="s" variable_1="s"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="transient_outward_current_r_gate" component_1="transient_outward_current"/>
<map_variables variable_2="r" variable_1="r"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="intracellular_sodium_concentration"/>
<map_variables variable_2="V_i" variable_1="V_i"/>
</connection>
<connection>
<map_components component_2="intracellular_calcium_concentration" component_1="intracellular_potassium_concentration"/>
<map_variables variable_2="V_i" variable_1="V_i"/>
</connection>
<rdf:RDF>
<rdf:Seq rdf:about="rdf:#32f313ee-3a26-4cd5-b916-7f25624cc712">
<rdf:li rdf:resource="rdf:#c067f7bf-13da-4fee-897c-76b499face0f"/>
<rdf:li rdf:resource="rdf:#aa1e031c-f91a-4f2b-bdb0-9935efb2fb88"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#2db40e56-79b9-4436-8f11-428fbeefd30d">
<dc:title>Journal of Physiology</dc:title>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>Department of Physiology, Anatomy & Genetics, University of Oxford</dc:publisher>
<cmeta:comment rdf:resource="rdf:#6acd60b6-13a2-48be-b0e4-400fc9f66731"/>
<dcterms:created rdf:resource="rdf:#da292423-b0bd-45d4-b3cd-e353cf43d806"/>
<dc:creator rdf:resource="rdf:#fe5a656e-2eb6-45ac-ae0a-243fb2ba9e5e"/>
<cmeta:modification rdf:resource="rdf:#5d7eec37-9b3f-4639-92ba-04ca1157d258"/>
<cmeta:modification rdf:resource="rdf:#308faf4e-60ac-47bc-9e31-fcb57dfe50a5"/>
<cmeta:modification rdf:resource="rdf:#ce0f379b-9ef1-46bb-b80d-11ce4710161d"/>
</rdf:Description>
<rdf:Description rdf:about="#noble_noble_2001">
<ns7:simulation rdf:resource="rdf:#$TOyx12"/>
<bqs:reference rdf:resource="rdf:#baae06ff-7a2e-492c-83e0-97c33444f1c9"/>
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>calcium dynamics</rdf:li>
<rdf:li>cardiac</rdf:li>
<rdf:li>electrophysiology</rdf:li>
<rdf:li>ventricular myocyte</rdf:li>
<rdf:li>guinea pig</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<cmeta:comment rdf:resource="rdf:#d9dbe043-322d-4a5b-8141-3bb322ab15c2"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#76b620ed-c827-4594-b1d0-c567267b0e50">
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
<rdf:value>penny.noble@dpag.ox.ac.uk</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#64d5a12b-304f-487a-b048-ffc31c0482f4">
<vCard:N rdf:resource="rdf:#f344d2c4-cd63-4827-9532-65e8568e90b0"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#45171d91-e70f-4aa9-a29e-ec87aacdcea3">
<vCard:N rdf:resource="rdf:#86f82cf5-520a-4107-8e1f-3e3e1c162e0b"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#89c922db-fae5-41e0-aaf0-87dbfe639ba1">
<dcterms:W3CDTF>2007-09-06T10:33:50+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6acd60b6-13a2-48be-b0e4-400fc9f66731">
<dc:creator rdf:resource="rdf:#efc1b103-a476-48a3-9297-f64897cec6e9"/>
<rdf:value>This version of the model has been curated and unit checked and is known to run in PCEnv, COR and Jsim. A PCEnv session file is also available for this model.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f344d2c4-cd63-4827-9532-65e8568e90b0">
<vCard:Given>James</vCard:Given>
<vCard:Family>Lawson</vCard:Family>
<vCard:Other>Richard</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fe5a656e-2eb6-45ac-ae0a-243fb2ba9e5e">
<vCard:ORG rdf:resource="rdf:#6308a8ff-044a-4ccb-97e1-bec0b0487b7b"/>
<vCard:EMAIL rdf:resource="rdf:#76b620ed-c827-4594-b1d0-c567267b0e50"/>
<vCard:N rdf:resource="rdf:#e383ffe6-cb7f-41ed-9763-2a53f0c171c5"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$VOyx12">
<ns7:nonstandard-pointDensity>100000</ns7:nonstandard-pointDensity>
<ns7:maximumStepSize>0.001</ns7:maximumStepSize>
</rdf:Description>
<rdf:Description rdf:about="rdf:#baae06ff-7a2e-492c-83e0-97c33444f1c9">
<bqs:JournalArticle rdf:resource="rdf:#ad36d3f7-e884-4125-99fc-36cd188399d1"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#abb451a0-aa82-4405-9e72-f8b0c01d92b5">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
<vCard:Other/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#da292423-b0bd-45d4-b3cd-e353cf43d806">
<dcterms:W3CDTF>2007-05-17T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#308faf4e-60ac-47bc-9e31-fcb57dfe50a5">
<dcterms:modified rdf:resource="rdf:#89c922db-fae5-41e0-aaf0-87dbfe639ba1"/>
<rdf:value>INaCads was not included in the dCads/dt ODE and hence the system was out
of balance. This issue has been fixed in version 05</rdf:value>
<cmeta:modifier rdf:resource="rdf:#64d5a12b-304f-487a-b048-ffc31c0482f4"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c067f7bf-13da-4fee-897c-76b499face0f">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#e1bc35af-e3aa-47a9-b10c-e9b8f51ad0fd"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8996b364-9f31-4ff6-ac94-2ff65799b92b">
<dcterms:W3CDTF>2001-01-01</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ce0f379b-9ef1-46bb-b80d-11ce4710161d">
<dcterms:modified rdf:resource="rdf:#fa8cc3a7-9cad-4cac-9914-69ee14ac6e27"/>
<rdf:value>Units checked, curated.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#45171d91-e70f-4aa9-a29e-ec87aacdcea3"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#089609f2-eee9-4e27-8b6f-73963d094ac8">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d9dbe043-322d-4a5b-8141-3bb322ab15c2">
<dc:creator rdf:resource="rdf:#e96b3e9a-a8b4-4ca5-94e8-d2a2a7cb951a"/>
<rdf:value>The Noble98 (Noble et al. 1998) model of the guinea-pig ventricular cell was a development of the Noble91 (Noble et al. 1991) model and included representation of the components of IK, (IKr and IKs), the persistent sodium current and a dyadic space between the sarcolemma and the sarcoplasmic reticulum. In several important respects this model was a major improvement on the earlier model. However on rigorous testing for new applications problems associated with the calcium dynamics equations were found in the new model.
The model demonstrated here is an attempt at 'correcting' the inability of Noble98 to reproduce a mechanism of initiation of Early After-Depolarizations (EADs), the voltage-dependence of calcium release from the sarcoplasmic reticulum (SR) and also to improve the shape of the calcium transient.
</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#aa1e031c-f91a-4f2b-bdb0-9935efb2fb88">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#089609f2-eee9-4e27-8b6f-73963d094ac8"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e96b3e9a-a8b4-4ca5-94e8-d2a2a7cb951a">
<vCard:FN>James Lawson</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e383ffe6-cb7f-41ed-9763-2a53f0c171c5">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5d7eec37-9b3f-4639-92ba-04ca1157d258">
<dcterms:modified rdf:resource="rdf:#60453a7d-4dc0-43e5-8386-ec2ee9bf0c7c"/>
<rdf:value>Note that the missing i_NCX_ds component was added to the dCai/dt rather than the dCads/dt equation in order to achieve balance in the model.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#463d4ed0-1fc1-4070-82d6-16f894fc341f"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fa8cc3a7-9cad-4cac-9914-69ee14ac6e27">
<dcterms:W3CDTF>2007-12-04T15:19:56+13:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e1bc35af-e3aa-47a9-b10c-e9b8f51ad0fd">
<vCard:Given>Dennis</vCard:Given>
<vCard:Family>Noble</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ad36d3f7-e884-4125-99fc-36cd188399d1">
<dc:title>Remodelling of Calcium Dynamics in Guinea-Pig Ventricular Cells</dc:title>
<dcterms:issued rdf:resource="rdf:#8996b364-9f31-4ff6-ac94-2ff65799b92b"/>
<bqs:Journal rdf:resource="rdf:#2db40e56-79b9-4436-8f11-428fbeefd30d"/>
<dc:creator rdf:resource="rdf:#32f313ee-3a26-4cd5-b916-7f25624cc712"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#efc1b103-a476-48a3-9297-f64897cec6e9">
<vCard:FN>James Lawson</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#60453a7d-4dc0-43e5-8386-ec2ee9bf0c7c">
<dcterms:W3CDTF>2008-05-14T03:19:46+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#463d4ed0-1fc1-4070-82d6-16f894fc341f">
<vCard:N rdf:resource="rdf:#abb451a0-aa82-4405-9e72-f8b0c01d92b5"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#86f82cf5-520a-4107-8e1f-3e3e1c162e0b">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$UOyx12">
<rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
<rdf:first rdf:resource="rdf:#$VOyx12"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$TOyx12">
<ns7:boundIntervals rdf:resource="rdf:#$UOyx12"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6308a8ff-044a-4ccb-97e1-bec0b0487b7b">
<vCard:Orgname>Oxford University</vCard:Orgname>
<vCard:Orgunit>Department of Physiology, Anatomy & Genetics</vCard:Orgunit>
</rdf:Description>
</rdf:RDF>
</model>
