<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 13/06/2008 at 1:52:31 at p.m. using:
COR (0.9.31.981)
Copyright 2002-2008 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk
CellML 1.0 was used to generate this 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#" name="pandit_clark_giles_demir_2001" cmeta:id="pandit_clark_giles_demir_2001">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Adult Rat Left Ventricular Myocyte Model, 2001</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This CellML version of the model represents the epicardial cell. In addition to the formulation given by the paper and the author's later Corrections document, the IStim current has been adjusted to produce 1 Hz stimulations for 10 seconds.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: Mathematical models were developed to reconstruct the action potentials (AP) recorded in epicardial and endocardial myocytes isolated from the adult rat left ventricle. The main goal was to obtain additional insight into the ionic mechanisms responsible for the transmural AP heterogeneity. The simulation results support the hypothesis that the smaller density and the slower reactivation kinetics of the Ca(2+)-independent transient outward K(+) current (I(t)) in the endocardial myocytes can account for the longer action potential duration (APD), and more prominent rate dependence in that cell type. The larger density of the Na(+) current (I(Na)) in the endocardial myocytes results in a faster upstroke (dV/dt(max)). This, in addition to the smaller magnitude of I(t), is responsible for the larger peak overshoot of the simulated endocardial AP. The prolonged APD in the endocardial cell also leads to an enhanced amplitude of the sustained K(+) current (I(ss)), and a larger influx of Ca(2+) ions via the L-type Ca(2+) current (I(CaL)). The latter results in an increased sarcoplasmic reticulum (SR) load, which is mainly responsible for the higher peak systolic value of the Ca(2+) transient [Ca(2+)](i), and the resultant increase in the Na(+)-Ca(2+) exchanger (I(NaCa)) activity, associated with the simulated endocardial AP. In combination, these calculations provide novel, quantitative insights into the repolarization process and its naturally occurring transmural variations in the rat left ventricle.
</para>
<para>
The original paper reference is cited below:
</para>
<para>
A Mathematical Model of Action Potential Heterogeneity in Adult Rat Left Ventricular Myocytes, Sandeep V. Pandit, Robert B. Clark, Wayne R. Giles and Semahat S. Demir, 2001, <emphasis>Biophysical Journal</emphasis>, 81, 3029-3051.<ulink url="http://www.ncbi.nlm.nih.gov/pubmed/11720973">PubMed ID: 11720973</ulink>
</para>
<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>fluid compartment model of the rat epicardial/endocardial ventricular cell</title>
</objectinfo>
<imagedata fileref="pandit_2001b.png"/>
</imageobject>
</mediaobject>
<caption>A fluid compartment model of the rat epicardial/endocardial ventricular cell.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<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_millivolt2">
<unit units="volt" prefix="milli" exponent="-2"/>
</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="millimolar">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="millimolar_per_second">
<unit units="millimolar"/>
<unit units="second" exponent="-1"/>
</units>
<units name="per_millimolar_per_second">
<unit units="millimolar" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="per_millimolar3_per_second">
<unit units="millimolar" exponent="-3"/>
<unit units="second" exponent="-1"/>
</units>
<units name="per_millimolar4">
<unit units="millimolar" exponent="-4"/>
</units>
<units name="per_millimolar4_per_second">
<unit units="per_millimolar4"/>
<unit units="second" exponent="-1"/>
</units>
<units name="nanoA_per_millimolar4">
<unit units="nanoA"/>
<unit units="per_millimolar4"/>
</units>
<units name="millijoule_per_mole_kelvin">
<unit units="joule" prefix="milli"/>
<unit units="kelvin" exponent="-1"/>
<unit units="mole" exponent="-1"/>
</units>
<units name="coulomb_per_mole">
<unit units="coulomb"/>
<unit units="mole" exponent="-1"/>
</units>
<units name="micro_litre">
<unit units="litre" prefix="micro"/>
</units>
<component name="environment">
<variable units="second" public_interface="out" name="time" cmeta:id="environment_time"/>
</component>
<component name="membrane">
<variable units="millivolt" public_interface="out" cmeta:id="membrane_V" name="V" initial_value="-80.50146"/>
<variable units="millijoule_per_mole_kelvin" public_interface="out" name="R" initial_value="8314.5"/>
<variable units="kelvin" public_interface="out" name="T" initial_value="295"/>
<variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96487"/>
<variable units="microF" name="Cm" initial_value="0.0001"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="nanoA" public_interface="in" name="i_Na"/>
<variable units="nanoA" public_interface="in" name="i_Ca_L"/>
<variable units="nanoA" public_interface="in" name="i_t"/>
<variable units="nanoA" public_interface="in" name="i_ss"/>
<variable units="nanoA" public_interface="in" name="i_f"/>
<variable units="nanoA" public_interface="in" name="i_K1"/>
<variable units="nanoA" public_interface="in" name="i_B"/>
<variable units="nanoA" public_interface="in" name="i_NaK"/>
<variable units="nanoA" public_interface="in" name="i_NaCa"/>
<variable units="nanoA" public_interface="in" name="i_Ca_P"/>
<variable units="nanoA" name="i_Stim"/>
<variable units="second" name="stim_start" initial_value="0.1"/>
<variable units="second" name="stim_end" initial_value="10"/>
<variable units="second" name="stim_period" initial_value="1"/>
<variable units="second" name="stim_duration" initial_value="0.005"/>
<variable units="nanoA" name="stim_amplitude" initial_value="-0.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>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_Na</ci>
<ci>i_Ca_L</ci>
<ci>i_t</ci>
<ci>i_ss</ci>
<ci>i_f</ci>
<ci>i_K1</ci>
<ci>i_B</ci>
<ci>i_NaK</ci>
<ci>i_NaCa</ci>
<ci>i_Ca_P</ci>
<ci>i_Stim</ci>
</apply>
</apply>
<ci>Cm</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_current">
<variable units="nanoA" public_interface="out" name="i_Na" cmeta:id="sodium_current_i_Na"/>
<variable units="millivolt" public_interface="out" name="E_Na"/>
<variable units="microS" name="g_Na" initial_value="0.8"/>
<variable units="millimolar" public_interface="in" name="Na_i"/>
<variable units="millimolar" public_interface="in" name="Na_o"/>
<variable units="millijoule_per_mole_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="m"/>
<variable units="dimensionless" private_interface="in" name="h"/>
<variable units="dimensionless" private_interface="in" name="j"/>
<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>i_Na</ci>
<apply>
<times/>
<ci>g_Na</ci>
<apply>
<power/>
<ci>m</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>h</ci>
<ci>j</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_m_gate">
<variable units="dimensionless" public_interface="out" name="m" initial_value="0.004164108"/>
<variable units="dimensionless" name="m_max"/>
<variable units="second" name="tau_m"/>
<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>m_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">45</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">6.5</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_m</ci>
<apply>
<divide/>
<cn cellml:units="second">0.00136</cn>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.32</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">47.13</cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</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">47.13</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.08</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">11</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>m</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>m_max</ci>
<ci>m</ci>
</apply>
<ci>tau_m</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_h_gate">
<variable units="dimensionless" public_interface="out" name="h" initial_value="0.6735613"/>
<variable units="dimensionless" name="h_max"/>
<variable units="second" name="tau_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>h_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">76.1</cn>
</apply>
<cn cellml:units="millivolt">6.07</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_h</ci>
<piecewise>
<piece>
<apply>
<times/>
<cn cellml:units="second">0.0004537</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">10.66</cn>
</apply>
</apply>
<cn cellml:units="millivolt">11.1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<geq/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<cn cellml:units="second">0.00349</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.135</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">80</cn>
</apply>
</apply>
<cn cellml:units="millivolt">6.8</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">3.56</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.079</cn>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">310000</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.35</cn>
<ci>V</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>h</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>h_max</ci>
<ci>h</ci>
</apply>
<ci>tau_h</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_j_gate">
<variable units="dimensionless" public_interface="out" name="j" initial_value="0.6729362"/>
<variable units="dimensionless" name="j_max"/>
<variable units="second" name="tau_j"/>
<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>j_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">76.1</cn>
</apply>
<cn cellml:units="millivolt">6.07</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_j</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="second">0.01163</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</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">32</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.0000002535</cn>
</apply>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<geq/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<cn cellml:units="second">0.00349</cn>
<apply>
<plus/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">37.78</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.311</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">79.23</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">127140</cn>
</apply>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.2444</cn>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.00003474</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.04391</cn>
</apply>
<ci>V</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.1212</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.01052</cn>
</apply>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.1378</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">40.14</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>j</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>j_max</ci>
<ci>j</ci>
</apply>
<ci>tau_j</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel">
<variable units="nanoA" public_interface="out" name="i_Ca_L" cmeta:id="L_type_Ca_channel_i_Ca_L"/>
<variable units="microS" name="g_Ca_L" initial_value="0.031"/>
<variable units="millivolt" name="E_Ca_L" initial_value="65"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="millimolar" public_interface="in" private_interface="out" name="Ca_ss"/>
<variable units="dimensionless" private_interface="in" name="d"/>
<variable units="dimensionless" private_interface="in" name="f_11"/>
<variable units="dimensionless" private_interface="in" name="f_12"/>
<variable units="dimensionless" private_interface="in" name="Ca_inact"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ca_L</ci>
<apply>
<times/>
<ci>g_Ca_L</ci>
<ci>d</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">0.9</cn>
<apply>
<divide/>
<ci>Ca_inact</ci>
<cn cellml:units="dimensionless">10</cn>
</apply>
</apply>
<ci>f_11</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.1</cn>
<apply>
<divide/>
<ci>Ca_inact</ci>
<cn cellml:units="dimensionless">10</cn>
</apply>
</apply>
<ci>f_12</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Ca_L</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_d_gate">
<variable units="dimensionless" public_interface="out" name="d" initial_value="0.000002171081"/>
<variable units="dimensionless" name="d_max"/>
<variable units="second" name="tau_d"/>
<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>d_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">15.3</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_d</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="second">0.00305</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt2">0.0045</cn>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">7</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="second">0.00105</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt2">0.002</cn>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">18</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.00025</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>d</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>d_max</ci>
<ci>d</ci>
</apply>
<ci>tau_d</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_f_11_gate">
<variable units="dimensionless" public_interface="out" name="f_11" initial_value="0.9999529"/>
<variable units="dimensionless" name="f_11_max"/>
<variable units="second" name="tau_f_11"/>
<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>f_11_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">26.7</cn>
</apply>
<cn cellml:units="millivolt">5.4</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_f_11</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="second">0.105</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">45</cn>
</apply>
<cn cellml:units="millivolt">12</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="second">0.04</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">25</cn>
</apply>
<cn cellml:units="millivolt">25</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="second">0.015</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
<cn cellml:units="millivolt">25</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.0017</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>f_11</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>f_11_max</ci>
<ci>f_11</ci>
</apply>
<ci>tau_f_11</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_f_12_gate">
<variable units="dimensionless" public_interface="out" name="f_12" initial_value="0.9999529"/>
<variable units="dimensionless" name="f_12_max"/>
<variable units="second" name="tau_f_12"/>
<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>f_12_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">26.7</cn>
</apply>
<cn cellml:units="millivolt">5.4</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_f_12</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="second">0.041</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">47</cn>
</apply>
<cn cellml:units="millivolt">12</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="second">0.08</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">55</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="second">0.015</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
<cn cellml:units="millivolt">25</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.0017</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>f_12</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>f_12_max</ci>
<ci>f_12</ci>
</apply>
<ci>tau_f_12</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_Ca_inact_gate">
<variable units="dimensionless" public_interface="out" name="Ca_inact" initial_value="0.9913102"/>
<variable units="second" name="tau_Ca_inact" initial_value="0.009"/>
<variable units="dimensionless" name="Ca_inact_max"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millimolar" public_interface="in" name="Ca_ss"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Ca_inact_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>Ca_ss</ci>
<cn cellml:units="millimolar">0.01</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_inact</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>Ca_inact_max</ci>
<ci>Ca_inact</ci>
</apply>
<ci>tau_Ca_inact</ci>
</apply>
</apply>
</math>
</component>
<component name="Ca_independent_transient_outward_K_current">
<variable units="nanoA" public_interface="out" name="i_t" cmeta:id="Ca_independent_transient_outward_K_current_i_t"/>
<variable units="millivolt" public_interface="out" name="E_K"/>
<variable units="microS" name="g_t" initial_value="0.035"/>
<variable units="dimensionless" name="a" initial_value="0.886"/>
<variable units="dimensionless" name="b" initial_value="0.114"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="millijoule_per_mole_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="millimolar" public_interface="in" name="K_o"/>
<variable units="millimolar" public_interface="in" name="K_i"/>
<variable units="dimensionless" private_interface="in" name="r"/>
<variable units="dimensionless" private_interface="in" name="s"/>
<variable units="dimensionless" private_interface="in" name="s_slow"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<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>i_t</ci>
<apply>
<times/>
<ci>g_t</ci>
<ci>r</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>a</ci>
<ci>s</ci>
</apply>
<apply>
<times/>
<ci>b</ci>
<ci>s_slow</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Ca_independent_transient_outward_K_current_r_gate">
<variable units="dimensionless" public_interface="out" name="r" initial_value="0.002191519"/>
<variable units="second" name="tau_r"/>
<variable units="dimensionless" name="r_max"/>
<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>r_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">10.6</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">11.42</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_r</ci>
<apply>
<divide/>
<cn cellml:units="second">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">45.16</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.03577</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">98.9</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">38</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>r_max</ci>
<ci>r</ci>
</apply>
<ci>tau_r</ci>
</apply>
</apply>
</math>
</component>
<component name="Ca_independent_transient_outward_K_current_s_gate">
<variable units="dimensionless" public_interface="out" name="s" initial_value="0.9842542"/>
<variable units="second" name="tau_s"/>
<variable units="dimensionless" name="s_max"/>
<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>s_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">45.3</cn>
</apply>
<cn cellml:units="millivolt">6.8841</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_s</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="second">0.35</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">70</cn>
</apply>
<cn cellml:units="millivolt">15</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.035</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>s_max</ci>
<ci>s</ci>
</apply>
<ci>tau_s</ci>
</apply>
</apply>
</math>
</component>
<component name="Ca_independent_transient_outward_K_current_s_slow_gate">
<variable units="dimensionless" public_interface="out" name="s_slow" initial_value="0.6421196"/>
<variable units="second" name="tau_s_slow"/>
<variable units="dimensionless" name="s_slow_max"/>
<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>s_slow_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">45.3</cn>
</apply>
<cn cellml:units="millivolt">6.8841</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_s_slow</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="second">3.7</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">70</cn>
</apply>
<cn cellml:units="millivolt">30</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.035</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s_slow</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>s_slow_max</ci>
<ci>s_slow</ci>
</apply>
<ci>tau_s_slow</ci>
</apply>
</apply>
</math>
</component>
<component name="steady_state_outward_K_current">
<variable units="nanoA" public_interface="out" name="i_ss" cmeta:id="steady_state_outward_K_current_i_ss"/>
<variable units="microS" name="g_ss" initial_value="0.007"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<variable units="dimensionless" private_interface="in" name="r_ss"/>
<variable units="dimensionless" private_interface="in" name="s_ss"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_ss</ci>
<apply>
<times/>
<ci>g_ss</ci>
<ci>r_ss</ci>
<ci>s_ss</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="steady_state_outward_K_current_r_ss_gate">
<variable units="dimensionless" public_interface="out" name="r_ss" initial_value="0.002907171"/>
<variable units="second" name="tau_r_ss"/>
<variable units="dimensionless" name="r_ss_max"/>
<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>r_ss_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">11.5</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">11.82</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_r_ss</ci>
<apply>
<divide/>
<cn cellml:units="second">10</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">45.16</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.03577</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">98.9</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">38</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r_ss</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>r_ss_max</ci>
<ci>r_ss</ci>
</apply>
<ci>tau_r_ss</ci>
</apply>
</apply>
</math>
</component>
<component name="steady_state_outward_K_current_s_ss_gate">
<variable units="dimensionless" public_interface="out" name="s_ss" initial_value="0.3142767"/>
<variable units="second" name="tau_s_ss"/>
<variable units="dimensionless" name="s_ss_max"/>
<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>s_ss_max</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">87.5</cn>
</apply>
<cn cellml:units="millivolt">10.3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_s_ss</ci>
<cn cellml:units="second">2.1</cn>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s_ss</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>s_ss_max</ci>
<ci>s_ss</ci>
</apply>
<ci>tau_s_ss</ci>
</apply>
</apply>
</math>
</component>
<component name="inward_rectifier">
<variable units="nanoA" public_interface="out" name="i_K1" cmeta:id="inward_rectifier_i_K1"/>
<variable units="microS" name="g_K1" initial_value="0.024"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millijoule_per_mole_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="millimolar" public_interface="in" name="K_o"/>
<variable units="millimolar" public_interface="in" name="K_i"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_K1</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="nanoA">48</cn>
<apply>
<plus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">37</cn>
</apply>
<cn cellml:units="millivolt">25</cn>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">37</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">25</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="nanoA">10</cn>
</apply>
<cn cellml:units="dimensionless">0.001</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
<apply>
<plus/>
<ci>E_K</ci>
<cn cellml:units="millivolt">76.77</cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">17</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>g_K1</ci>
<apply>
<minus/>
<ci>V</ci>
<apply>
<plus/>
<ci>E_K</ci>
<cn cellml:units="millivolt">1.73</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1.613</cn>
<ci>F</ci>
<apply>
<minus/>
<ci>V</ci>
<apply>
<plus/>
<ci>E_K</ci>
<cn cellml:units="millivolt">1.73</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>K_o</ci>
<cn cellml:units="millimolar">0.9988</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millimolar">0.124</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="hyperpolarisation_activated_current">
<variable units="nanoA" public_interface="out" name="i_f" cmeta:id="hyperpolarisation_activated_current_i_f"/>
<variable units="nanoA" public_interface="out" name="i_f_Na"/>
<variable units="nanoA" public_interface="out" name="i_f_K"/>
<variable units="microS" name="g_f" initial_value="0.00145"/>
<variable units="dimensionless" name="f_Na" initial_value="0.2"/>
<variable units="dimensionless" name="f_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="millivolt" public_interface="in" name="E_K"/>
<variable units="millivolt" public_interface="in" name="E_Na"/>
<variable units="dimensionless" private_interface="in" name="y"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>f_K</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>f_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_f_Na</ci>
<apply>
<times/>
<ci>g_f</ci>
<ci>y</ci>
<ci>f_Na</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_f_K</ci>
<apply>
<times/>
<ci>g_f</ci>
<ci>y</ci>
<ci>f_K</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_f</ci>
<apply>
<plus/>
<ci>i_f_Na</ci>
<ci>i_f_K</ci>
</apply>
</apply>
</math>
</component>
<component name="hyperpolarisation_activated_current_y_gate">
<variable units="dimensionless" public_interface="out" name="y" initial_value="0.003578708"/>
<variable units="second" name="tau_y"/>
<variable units="dimensionless" name="y_infinity"/>
<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>y_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">138.6</cn>
</apply>
<cn cellml:units="millivolt">10.48</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_y</ci>
<apply>
<divide/>
<cn cellml:units="second">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.11885</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">80</cn>
</apply>
<cn cellml:units="millivolt">28.37</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5623</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">80</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">14.19</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>y</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>y_infinity</ci>
<ci>y</ci>
</apply>
<ci>tau_y</ci>
</apply>
</apply>
</math>
</component>
<component name="background_currents">
<variable units="nanoA" public_interface="out" name="i_B" cmeta:id="background_currents_i_B"/>
<variable units="nanoA" public_interface="out" name="i_B_Na"/>
<variable units="nanoA" public_interface="out" name="i_B_Ca"/>
<variable units="nanoA" public_interface="out" name="i_B_K"/>
<variable units="microS" name="g_B_Na" initial_value="0.00008015"/>
<variable units="microS" name="g_B_Ca" initial_value="0.0000324"/>
<variable units="microS" name="g_B_K" initial_value="0.000138"/>
<variable units="millivolt" name="E_Ca" initial_value="65"/>
<variable units="millivolt" public_interface="in" name="E_Na"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<variable units="millimolar" public_interface="in" name="Ca_o"/>
<variable units="millimolar" public_interface="in" name="Ca_i"/>
<variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<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_B_Na</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_B_Ca</ci>
<apply>
<times/>
<ci>g_B_Ca</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Ca</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_B_K</ci>
<apply>
<times/>
<ci>g_B_K</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_B</ci>
<apply>
<plus/>
<ci>i_B_Na</ci>
<ci>i_B_Ca</ci>
<ci>i_B_K</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_potassium_pump">
<variable units="nanoA" public_interface="out" name="i_NaK" cmeta:id="sodium_potassium_pump_i_NaK"/>
<variable units="nanoA" name="i_NaK_max" initial_value="0.08"/>
<variable units="millimolar" name="K_m_K" initial_value="1.5"/>
<variable units="millimolar" name="K_m_Na" initial_value="10"/>
<variable units="dimensionless" name="sigma"/>
<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="Na_i"/>
<variable units="millijoule_per_mole_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"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>sigma</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Na_o</ci>
<cn cellml:units="millimolar">67.3</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
<cn cellml:units="dimensionless">7</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>i_NaK</ci>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>i_NaK_max</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.1245</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.1</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.0365</cn>
<ci>sigma</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>K_o</ci>
</apply>
<apply>
<plus/>
<ci>K_o</ci>
<ci>K_m_K</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>K_m_Na</ci>
<ci>Na_i</ci>
</apply>
<cn cellml:units="dimensionless">1.5</cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sarcolemmal_calcium_pump_current">
<variable units="nanoA" public_interface="out" name="i_Ca_P" cmeta:id="sarcolemmal_calcium_pump_current_i_Ca_P"/>
<variable units="nanoA" name="i_Ca_P_max" initial_value="0.004"/>
<variable units="millimolar" public_interface="in" name="Ca_i"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ca_P</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>i_Ca_P_max</ci>
<ci>Ca_i</ci>
</apply>
<apply>
<plus/>
<ci>Ca_i</ci>
<cn cellml:units="millimolar">0.0004</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Na_Ca_ion_exchanger_current">
<variable units="nanoA" public_interface="out" name="i_NaCa" cmeta:id="Na_Ca_ion_exchanger_current_i_NaCa"/>
<variable units="nanoA_per_millimolar4" name="K_NaCa" initial_value="0.000009984"/>
<variable units="per_millimolar4" name="d_NaCa" initial_value="0.0001"/>
<variable units="dimensionless" name="gamma_NaCa" initial_value="0.5"/>
<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_o"/>
<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</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>K_NaCa</ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>Ca_o</ci>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.03743</cn>
<ci>V</ci>
<ci>gamma_NaCa</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>Na_o</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>Ca_i</ci>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.03743</cn>
<ci>V</ci>
<apply>
<minus/>
<ci>gamma_NaCa</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>d_NaCa</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>Ca_i</ci>
<apply>
<power/>
<ci>Na_o</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_o</ci>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SR_Ca_release_channel">
<variable units="millimolar_per_second" public_interface="out" name="J_rel" cmeta:id="SR_Ca_release_channel_J_rel"/>
<variable units="per_second" name="v1" initial_value="1.8e3"/>
<variable units="per_millimolar4_per_second" name="k_a_plus" initial_value="12.15e12"/>
<variable units="per_second" name="k_a_minus" initial_value="576"/>
<variable units="per_millimolar3_per_second" name="k_b_plus" initial_value="4.05e9"/>
<variable units="per_second" name="k_b_minus" initial_value="1930"/>
<variable units="per_second" name="k_c_plus" initial_value="100"/>
<variable units="per_second" name="k_c_minus" initial_value="0.8"/>
<variable units="dimensionless" name="P_O1" initial_value="0.0004327548"/>
<variable units="dimensionless" name="P_O2" initial_value="0.000000000606254"/>
<variable units="dimensionless" name="P_C1" initial_value="0.6348229"/>
<variable units="dimensionless" name="P_C2" initial_value="0.3647471"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millimolar" public_interface="in" name="Ca_ss"/>
<variable units="millimolar" public_interface="in" name="Ca_JSR"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P_C1</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<ci>k_a_plus</ci>
</apply>
<apply>
<power/>
<ci>Ca_ss</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<ci>P_C1</ci>
</apply>
<apply>
<times/>
<ci>k_a_minus</ci>
<ci>P_O1</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P_O1</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>k_a_plus</ci>
<apply>
<power/>
<ci>Ca_ss</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
<ci>P_C1</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>k_a_minus</ci>
<ci>P_O1</ci>
</apply>
<apply>
<times/>
<ci>k_b_plus</ci>
<apply>
<power/>
<ci>Ca_ss</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>P_O1</ci>
</apply>
<apply>
<times/>
<ci>k_c_plus</ci>
<ci>P_O1</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>k_b_minus</ci>
<ci>P_O2</ci>
</apply>
<apply>
<times/>
<ci>k_c_minus</ci>
<ci>P_C2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P_O2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>k_b_plus</ci>
<apply>
<power/>
<ci>Ca_ss</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>P_O1</ci>
</apply>
<apply>
<times/>
<ci>k_b_minus</ci>
<ci>P_O2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P_C2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>k_c_plus</ci>
<ci>P_O1</ci>
</apply>
<apply>
<times/>
<ci>k_c_minus</ci>
<ci>P_C2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_rel</ci>
<apply>
<times/>
<ci>v1</ci>
<apply>
<plus/>
<ci>P_O1</ci>
<ci>P_O2</ci>
</apply>
<apply>
<minus/>
<ci>Ca_JSR</ci>
<ci>Ca_ss</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SERCA2a_pump">
<variable units="millimolar_per_second" public_interface="out" name="J_up" cmeta:id="SERCA2a_pump_J_up"/>
<variable units="millimolar" name="K_fb" initial_value="0.000168"/>
<variable units="millimolar" name="K_rb" initial_value="3.29"/>
<variable units="dimensionless" name="fb"/>
<variable units="dimensionless" name="rb"/>
<variable units="millimolar_per_second" name="Vmaxf" initial_value="0.04"/>
<variable units="millimolar_per_second" name="Vmaxr" initial_value="0.9"/>
<variable units="dimensionless" name="K_SR" initial_value="1"/>
<variable units="dimensionless" name="N_fb" initial_value="1.2"/>
<variable units="dimensionless" name="N_rb" initial_value="1"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millimolar" public_interface="in" name="Ca_i"/>
<variable units="millimolar" public_interface="in" name="Ca_NSR"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>fb</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>Ca_i</ci>
<ci>K_fb</ci>
</apply>
<ci>N_fb</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>rb</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>Ca_NSR</ci>
<ci>K_rb</ci>
</apply>
<ci>N_rb</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_up</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>K_SR</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>Vmaxf</ci>
<ci>fb</ci>
</apply>
<apply>
<times/>
<ci>Vmaxr</ci>
<ci>rb</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>fb</ci>
<ci>rb</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="intracellular_and_SR_Ca_fluxes">
<variable units="millimolar_per_second" public_interface="out" name="J_tr"/>
<variable units="millimolar_per_second" public_interface="out" name="J_xfer" cmeta:id="intracellular_and_SR_Ca_fluxes_J_xfer"/>
<variable units="millimolar_per_second" public_interface="out" name="J_trpn"/>
<variable units="second" name="tau_tr" initial_value="0.0005747"/>
<variable units="second" name="tau_xfer" initial_value="0.0267"/>
<variable units="millimolar" name="HTRPNCa" initial_value="1.394301e-1"/>
<variable units="millimolar" name="LTRPNCa" initial_value="5.1619e-3"/>
<variable units="millimolar_per_second" name="J_HTRPNCa"/>
<variable units="millimolar_per_second" name="J_LTRPNCa"/>
<variable units="millimolar" name="HTRPN_tot" initial_value="0.14"/>
<variable units="millimolar" name="LTRPN_tot" initial_value="0.07"/>
<variable units="per_millimolar_per_second" name="k_htrpn_plus" initial_value="200000"/>
<variable units="per_second" name="k_htrpn_minus" initial_value="0.066"/>
<variable units="per_millimolar_per_second" name="k_ltrpn_plus" initial_value="40000"/>
<variable units="per_second" name="k_ltrpn_minus" initial_value="40"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millimolar" public_interface="in" name="Ca_ss"/>
<variable units="millimolar" public_interface="in" name="Ca_i"/>
<variable units="millimolar" public_interface="in" name="Ca_NSR"/>
<variable units="millimolar" public_interface="in" name="Ca_JSR"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_tr</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Ca_NSR</ci>
<ci>Ca_JSR</ci>
</apply>
<ci>tau_tr</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_xfer</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Ca_ss</ci>
<ci>Ca_i</ci>
</apply>
<ci>tau_xfer</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_HTRPNCa</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_htrpn_plus</ci>
<ci>Ca_i</ci>
<apply>
<minus/>
<ci>HTRPN_tot</ci>
<ci>HTRPNCa</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k_htrpn_minus</ci>
<ci>HTRPNCa</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>HTRPNCa</ci>
</apply>
<ci>J_HTRPNCa</ci>
</apply>
<apply>
<eq/>
<ci>J_LTRPNCa</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>k_ltrpn_plus</ci>
<ci>Ca_i</ci>
<apply>
<minus/>
<ci>LTRPN_tot</ci>
<ci>LTRPNCa</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k_ltrpn_minus</ci>
<ci>LTRPNCa</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>LTRPNCa</ci>
</apply>
<ci>J_LTRPNCa</ci>
</apply>
<apply>
<eq/>
<ci>J_trpn</ci>
<apply>
<plus/>
<ci>J_HTRPNCa</ci>
<ci>J_LTRPNCa</ci>
</apply>
</apply>
</math>
</component>
<component name="intracellular_ion_concentrations">
<variable units="millimolar" public_interface="out" name="Na_i" initial_value="10.73519"/>
<variable units="millimolar" public_interface="out" name="Ca_i" initial_value="0.00007901351"/>
<variable units="millimolar" public_interface="out" name="K_i" initial_value="139.2751"/>
<variable units="millimolar" public_interface="out" name="Ca_ss" initial_value="0.00008737212"/>
<variable units="millimolar" public_interface="out" name="Ca_JSR" initial_value="0.06607948"/>
<variable units="millimolar" public_interface="out" name="Ca_NSR" initial_value="0.06600742"/>
<variable units="micro_litre" name="V_myo" initial_value="9.36e-6"/>
<variable units="micro_litre" name="V_JSR" initial_value="5.6e-8"/>
<variable units="micro_litre" name="V_NSR" initial_value="5.04e-7"/>
<variable units="micro_litre" name="V_SS" initial_value="1.2e-9"/>
<variable units="millimolar" name="K_mCMDN" initial_value="0.00238"/>
<variable units="millimolar" name="K_mCSQN" initial_value="0.8"/>
<variable units="millimolar" name="K_mEGTA" initial_value="0.00015"/>
<variable units="millimolar" name="CMDN_tot" initial_value="0.05"/>
<variable units="millimolar" name="CSQN_tot" initial_value="15"/>
<variable units="millimolar" name="EGTA_tot" initial_value="10"/>
<variable units="dimensionless" name="beta_i"/>
<variable units="dimensionless" name="beta_SS"/>
<variable units="dimensionless" name="beta_JSR"/>
<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_Ca_L"/>
<variable units="nanoA" public_interface="in" name="i_B_Na"/>
<variable units="nanoA" public_interface="in" name="i_NaCa"/>
<variable units="nanoA" public_interface="in" name="i_NaK"/>
<variable units="nanoA" public_interface="in" name="i_f_Na"/>
<variable units="nanoA" public_interface="in" name="i_f_K"/>
<variable units="nanoA" public_interface="in" name="i_B_K"/>
<variable units="nanoA" public_interface="in" name="i_K1"/>
<variable units="nanoA" public_interface="in" name="i_t"/>
<variable units="nanoA" public_interface="in" name="i_ss"/>
<variable units="nanoA" public_interface="in" name="i_Ca_P"/>
<variable units="nanoA" public_interface="in" name="i_B_Ca"/>
<variable units="millimolar_per_second" public_interface="in" name="J_up"/>
<variable units="millimolar_per_second" public_interface="in" name="J_rel"/>
<variable units="millimolar_per_second" public_interface="in" name="J_xfer"/>
<variable units="millimolar_per_second" public_interface="in" name="J_trpn"/>
<variable units="millimolar_per_second" public_interface="in" name="J_tr"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>beta_i</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>CMDN_tot</ci>
<ci>K_mCMDN</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>K_mCMDN</ci>
<ci>Ca_i</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>EGTA_tot</ci>
<ci>K_mEGTA</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>K_mEGTA</ci>
<ci>Ca_i</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_SS</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>CMDN_tot</ci>
<ci>K_mCMDN</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>K_mCMDN</ci>
<ci>Ca_ss</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_JSR</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>CSQN_tot</ci>
<ci>K_mCSQN</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>K_mCSQN</ci>
<ci>Ca_JSR</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<ci>beta_i</ci>
<apply>
<minus/>
<ci>J_xfer</ci>
<apply>
<plus/>
<ci>J_up</ci>
<ci>J_trpn</ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<minus/>
<ci>i_B_Ca</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>i_NaCa</ci>
</apply>
</apply>
<ci>i_Ca_P</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>V_myo</ci>
<ci>F</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Na_i</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_Na</ci>
<ci>i_B_Na</ci>
<apply>
<times/>
<ci>i_NaCa</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<apply>
<times/>
<ci>i_NaK</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>i_f_Na</ci>
</apply>
</apply>
<apply>
<times/>
<ci>V_myo</ci>
<ci>F</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_i</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>i_ss</ci>
<ci>i_B_K</ci>
<ci>i_t</ci>
<ci>i_K1</ci>
<ci>i_f_K</ci>
<apply>
<times/>
<ci>i_NaK</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>V_myo</ci>
<ci>F</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_ss</ci>
</apply>
<apply>
<times/>
<ci>beta_SS</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci>J_rel</ci>
<ci>V_JSR</ci>
</apply>
<ci>V_SS</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>J_xfer</ci>
<ci>V_myo</ci>
</apply>
<ci>V_SS</ci>
</apply>
</apply>
<apply>
<divide/>
<ci>i_Ca_L</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>V_SS</ci>
<ci>F</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_JSR</ci>
</apply>
<apply>
<times/>
<ci>beta_JSR</ci>
<apply>
<minus/>
<ci>J_tr</ci>
<ci>J_rel</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_NSR</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci>J_up</ci>
<ci>V_myo</ci>
</apply>
<ci>V_NSR</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>J_tr</ci>
<ci>V_JSR</ci>
</apply>
<ci>V_NSR</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="standard_ionic_concentrations">
<variable units="millimolar" public_interface="out" name="Na_o" initial_value="140"/>
<variable units="millimolar" public_interface="out" name="Ca_o" initial_value="1.2"/>
<variable units="millimolar" public_interface="out" name="K_o" initial_value="5.4"/>
</component>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="sodium_current">
<component_ref component="sodium_current_m_gate"/>
<component_ref component="sodium_current_h_gate"/>
<component_ref component="sodium_current_j_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_11_gate"/>
<component_ref component="L_type_Ca_channel_f_12_gate"/>
<component_ref component="L_type_Ca_channel_Ca_inact_gate"/>
</component_ref>
<component_ref component="Ca_independent_transient_outward_K_current">
<component_ref component="Ca_independent_transient_outward_K_current_r_gate"/>
<component_ref component="Ca_independent_transient_outward_K_current_s_gate"/>
<component_ref component="Ca_independent_transient_outward_K_current_s_slow_gate"/>
</component_ref>
<component_ref component="steady_state_outward_K_current">
<component_ref component="steady_state_outward_K_current_r_ss_gate"/>
<component_ref component="steady_state_outward_K_current_s_ss_gate"/>
</component_ref>
<component_ref component="hyperpolarisation_activated_current">
<component_ref component="hyperpolarisation_activated_current_y_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="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="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="steady_state_outward_K_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="inward_rectifier"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="hyperpolarisation_activated_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sarcolemmal_calcium_pump_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="Na_Ca_ion_exchanger_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="SR_Ca_release_channel"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="SERCA2a_pump"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="intracellular_and_SR_Ca_fluxes"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="intracellular_ion_concentrations"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="sodium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
<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" variable_1="i_Ca_L"/>
</connection>
<connection>
<map_components component_2="Ca_independent_transient_outward_K_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_t" variable_1="i_t"/>
<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="steady_state_outward_K_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_ss" variable_1="i_ss"/>
</connection>
<connection>
<map_components component_2="inward_rectifier" 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="hyperpolarisation_activated_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_f" variable_1="i_f"/>
</connection>
<connection>
<map_components component_2="background_currents" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_B" variable_1="i_B"/>
<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="sodium_potassium_pump" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_NaK" variable_1="i_NaK"/>
<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="sarcolemmal_calcium_pump_current" component_1="membrane"/>
<map_variables variable_2="i_Ca_P" variable_1="i_Ca_P"/>
</connection>
<connection>
<map_components component_2="Na_Ca_ion_exchanger_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="membrane"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="sodium_current"/>
<map_variables variable_2="i_Na" variable_1="i_Na"/>
<map_variables variable_2="Na_i" variable_1="Na_i"/>
</connection>
<connection>
<map_components component_2="background_currents" component_1="sodium_current"/>
<map_variables variable_2="E_Na" variable_1="E_Na"/>
</connection>
<connection>
<map_components component_2="hyperpolarisation_activated_current" component_1="sodium_current"/>
<map_variables variable_2="E_Na" variable_1="E_Na"/>
</connection>
<connection>
<map_components component_2="standard_ionic_concentrations" component_1="sodium_current"/>
<map_variables variable_2="Na_o" variable_1="Na_o"/>
</connection>
<connection>
<map_components component_2="standard_ionic_concentrations" component_1="background_currents"/>
<map_variables variable_2="Ca_o" variable_1="Ca_o"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="L_type_Ca_channel"/>
<map_variables variable_2="i_Ca_L" variable_1="i_Ca_L"/>
<map_variables variable_2="Ca_ss" variable_1="Ca_ss"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="i_t" variable_1="i_t"/>
<map_variables variable_2="K_i" variable_1="K_i"/>
</connection>
<connection>
<map_components component_2="steady_state_outward_K_current" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="inward_rectifier" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="hyperpolarisation_activated_current" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="background_currents" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="standard_ionic_concentrations" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="K_o" variable_1="K_o"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="steady_state_outward_K_current"/>
<map_variables variable_2="i_ss" variable_1="i_ss"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="inward_rectifier"/>
<map_variables variable_2="i_K1" variable_1="i_K1"/>
<map_variables variable_2="K_i" variable_1="K_i"/>
</connection>
<connection>
<map_components component_2="standard_ionic_concentrations" component_1="inward_rectifier"/>
<map_variables variable_2="K_o" variable_1="K_o"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="hyperpolarisation_activated_current"/>
<map_variables variable_2="i_f_Na" variable_1="i_f_Na"/>
<map_variables variable_2="i_f_K" variable_1="i_f_K"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="background_currents"/>
<map_variables variable_2="i_B_Na" variable_1="i_B_Na"/>
<map_variables variable_2="i_B_K" variable_1="i_B_K"/>
<map_variables variable_2="i_B_Ca" variable_1="i_B_Ca"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" 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="standard_ionic_concentrations" component_1="sodium_potassium_pump"/>
<map_variables variable_2="K_o" variable_1="K_o"/>
<map_variables variable_2="Na_o" variable_1="Na_o"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="sarcolemmal_calcium_pump_current"/>
<map_variables variable_2="i_Ca_P" variable_1="i_Ca_P"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="Na_Ca_ion_exchanger_current"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
<map_variables variable_2="Na_i" variable_1="Na_i"/>
<map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
</connection>
<connection>
<map_components component_2="standard_ionic_concentrations" component_1="Na_Ca_ion_exchanger_current"/>
<map_variables variable_2="Ca_o" variable_1="Ca_o"/>
<map_variables variable_2="Na_o" variable_1="Na_o"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="SR_Ca_release_channel"/>
<map_variables variable_2="J_rel" variable_1="J_rel"/>
<map_variables variable_2="Ca_ss" variable_1="Ca_ss"/>
<map_variables variable_2="Ca_JSR" variable_1="Ca_JSR"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="SERCA2a_pump"/>
<map_variables variable_2="J_up" variable_1="J_up"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
<map_variables variable_2="Ca_NSR" variable_1="Ca_NSR"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="intracellular_and_SR_Ca_fluxes"/>
<map_variables variable_2="J_tr" variable_1="J_tr"/>
<map_variables variable_2="J_xfer" variable_1="J_xfer"/>
<map_variables variable_2="J_trpn" variable_1="J_trpn"/>
<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
<map_variables variable_2="Ca_NSR" variable_1="Ca_NSR"/>
<map_variables variable_2="Ca_ss" variable_1="Ca_ss"/>
<map_variables variable_2="Ca_JSR" variable_1="Ca_JSR"/>
</connection>
<connection>
<map_components component_2="sodium_current_m_gate" component_1="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="sodium_current_h_gate" component_1="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="sodium_current_j_gate" component_1="sodium_current"/>
<map_variables variable_2="j" variable_1="j"/>
<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_11_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="f_11" variable_1="f_11"/>
<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_12_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="f_12" variable_1="f_12"/>
<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_Ca_inact_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="Ca_inact" variable_1="Ca_inact"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="Ca_ss" variable_1="Ca_ss"/>
</connection>
<connection>
<map_components component_2="Ca_independent_transient_outward_K_current_r_gate" component_1="Ca_independent_transient_outward_K_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="Ca_independent_transient_outward_K_current_s_gate" component_1="Ca_independent_transient_outward_K_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="Ca_independent_transient_outward_K_current_s_slow_gate" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="s_slow" variable_1="s_slow"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="steady_state_outward_K_current_r_ss_gate" component_1="steady_state_outward_K_current"/>
<map_variables variable_2="r_ss" variable_1="r_ss"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="steady_state_outward_K_current_s_ss_gate" component_1="steady_state_outward_K_current"/>
<map_variables variable_2="s_ss" variable_1="s_ss"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="hyperpolarisation_activated_current_y_gate" component_1="hyperpolarisation_activated_current"/>
<map_variables variable_2="y" variable_1="y"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<rdf:RDF>
<rdf:Seq rdf:about="rdf:#cc3ea2c4-6178-48e5-af65-7906d40d1fa4">
<rdf:li rdf:resource="rdf:#952a6e88-d239-4cd5-8db7-a0438bcc0557"/>
<rdf:li rdf:resource="rdf:#b720b6a7-d649-4704-b58c-3e03c078a323"/>
<rdf:li rdf:resource="rdf:#71546bf6-ff74-4c92-8272-878e3c7ff18b"/>
<rdf:li rdf:resource="rdf:#7c633a24-2630-4ed1-87fc-b33b8f35b18f"/>
</rdf:Seq>
<rdf:Seq rdf:about="rdf:#1316a37c-6c41-4c5f-8968-444b36ee3de5">
<rdf:li rdf:resource="rdf:#51eaaf89-1fe5-454e-81d4-7e5561c13e5c"/>
<rdf:li rdf:resource="rdf:#d13ad4d7-2b77-45c4-a989-9a2c3ccc204e"/>
<rdf:li rdf:resource="rdf:#359aa759-ca02-411f-ae3f-3da81bc0f064"/>
<rdf:li rdf:resource="rdf:#3a2362e4-af59-47e3-81e2-ff8a464cb4c3"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#91960679-9437-40bf-9215-85f10e927a1b">
<vCard:N rdf:resource="rdf:#a00fa64d-2c70-4560-8de5-182db60b50ae"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9996133d-52a8-465f-8c28-1da0b94d2322">
<vCard:Given>W</vCard:Given>
<vCard:Family>Giles</vCard:Family>
<vCard:Other>R</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3d34ba23-dc1f-4b23-b33d-578a6d2d8c57">
<bqs:Pubmed_id>11720973</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#d213868e-80c0-4c8a-bd5b-9e7f68be57fe"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d79d8c97-5b38-4c39-9e96-0aad2089274a">
<vCard:N rdf:resource="rdf:#19d6d867-3aae-40f9-b8ee-22c24fad5dc0"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$QvprG">
<ns7:boundIntervals rdf:resource="rdf:#$RvprG"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5fa88c97-91e4-486a-a985-229f233b984e">
<vCard:Given>S</vCard:Given>
<vCard:Family>Demir</vCard:Family>
<vCard:Other>S</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a89fcb96-f39b-4075-847b-e46f395266af">
<vCard:N rdf:resource="rdf:#819c7015-b101-47cf-a2a7-ec1ff211e163"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9dd15c9e-7d73-41d7-a27f-5a4585bbdad1">
<dc:title>Biophysical Journal</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#589a05f4-5e47-4059-90d3-e577f2fb8755">
<vCard:N rdf:resource="rdf:#c2b4910b-52ce-4333-9f1b-d570efb93892"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#dfcc85e4-f2cd-44b9-81f0-c1a1b38e1193">
<vCard:Given>Michael</vCard:Given>
<vCard:Family>Cooling</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2ecf70be-c03a-409e-ab2d-0f76aa6536c5">
<vCard:Given>James</vCard:Given>
<vCard:Family>Lawson</vCard:Family>
<vCard:Other>Richard</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7a7ac380-fabd-4cc9-b603-061dce055eea">
<dc:creator rdf:resource="rdf:#1316a37c-6c41-4c5f-8968-444b36ee3de5"/>
<dc:title>A mathematical model of action potential heterogeneity in adult rat left ventricular myocytes</dc:title>
<bqs:volume>81 6</bqs:volume>
<bqs:first_page>3029</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#e2d39fc6-49fd-4c39-93ec-3a3519b1d8ff"/>
<dcterms:issued rdf:resource="rdf:#d71f5fd8-902a-4360-b367-ee3d45f19635"/>
<bqs:last_page>3051</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#19d6d867-3aae-40f9-b8ee-22c24fad5dc0">
<vCard:Given>James</vCard:Given>
<vCard:Family>Lawson</vCard:Family>
<vCard:Other>Richard</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f41e4fd4-bf58-44e8-8f42-bf58bf1aae93">
<vCard:N rdf:resource="rdf:#1a06d1a1-b9f2-47f8-a61a-fac50ece9828"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8a9f2967-4170-417e-9653-efdad57d4d0c">
<vCard:Given>S</vCard:Given>
<vCard:Family>Demir</vCard:Family>
<vCard:Other>S</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:title>Rat ventricular myocyte model from the original Pandit 2001 paper: epicardial cell</dc:title>
<dc:publisher>Department of Physiology, Anatomy & Genetics, University of Oxford</dc:publisher>
<cmeta:comment rdf:resource="rdf:#59cb9d3e-3dca-45f1-9857-f9e5a07de4f0"/>
<dcterms:created rdf:resource="rdf:#c9681df6-408f-4fae-ae6f-5a3ff09c0475"/>
<dc:creator rdf:resource="rdf:#c4600241-ec5b-47ee-8294-30cb2854a60e"/>
<cmeta:modification rdf:resource="rdf:#4e120d40-0111-469a-8f20-e5bfb5ede2e4"/>
<cmeta:modification rdf:resource="rdf:#136ced93-18a7-4f9b-9628-1b29057940bc"/>
<cmeta:modification rdf:resource="rdf:#2f0fc38a-3c1c-4bfb-a830-47c3131b56c3"/>
<cmeta:modification rdf:resource="rdf:#4377fe6c-add1-493b-84bb-f38bb0c075c7"/>
<cmeta:modification rdf:resource="rdf:#43a014e2-42c1-44fe-82ac-55a21c5f0ed8"/>
<cmeta:modification rdf:resource="rdf:#6cf4bb45-eef1-4cad-965f-93ee42368494"/>
<cmeta:modification rdf:resource="rdf:#847270ad-ac16-478c-ba30-577f1e0483fe"/>
<cmeta:modification rdf:resource="rdf:#8ea99a18-ee92-4fda-a003-f64ae2947ee2"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2908c9a7-4e05-48a8-a4da-8f324e079c93">
<vCard:Given>R</vCard:Given>
<vCard:Family>Clark</vCard:Family>
<vCard:Other>B</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6324c2a1-51e4-4a74-8693-1021f04bc9b1">
<vCard:Orgname>Oxford University</vCard:Orgname>
<vCard:Orgunit>Department of Physiology, Anatomy & Genetics</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="#pandit_clark_giles_demir_2001">
<ns7:simulation rdf:resource="rdf:#$QvprG"/>
<cmeta:comment rdf:resource="rdf:#aa39cce2-2692-4f0c-9409-68cd0b6580f7"/>
<bqs:reference rdf:resource="rdf:#3d34ba23-dc1f-4b23-b33d-578a6d2d8c57"/>
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>electrophysiology</rdf:li>
<rdf:li>ventricular myocyte</rdf:li>
<rdf:li>cardiac</rdf:li>
<rdf:li>rat</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<bqs:reference rdf:resource="rdf:#c5725b6b-4fe6-4323-b8b5-c46dbd8b8635"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#010abb91-8c98-4c4b-9125-9611c9229f62">
<vCard:FN/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3a2362e4-af59-47e3-81e2-ff8a464cb4c3">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#5fa88c97-91e4-486a-a985-229f233b984e"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#48f139cb-d71c-4c43-b9fc-30a96302b9f5">
<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:#819c7015-b101-47cf-a2a7-ec1ff211e163">
<vCard:Given>Michael</vCard:Given>
<vCard:Family>Cooling</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b720b6a7-d649-4704-b58c-3e03c078a323">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#2908c9a7-4e05-48a8-a4da-8f324e079c93"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#51eaaf89-1fe5-454e-81d4-7e5561c13e5c">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#7702b560-a156-4e0a-b769-51c3ae395590"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8ea99a18-ee92-4fda-a003-f64ae2947ee2">
<dcterms:modified rdf:resource="rdf:#338a6b4b-01d1-4ad6-9f4e-4bd24d16edd1"/>
<rdf:value>Updated documentation image to full cell diagram.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#3801bb39-63d2-4ea7-b5a5-26d159ddbb04"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1b9f3170-1975-4598-acb3-6c6499b9c9dc">
<dcterms:W3CDTF>2008-05-26T11:07:32+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$SvprG">
<ns7:nonstandard-pointDensity>100000</ns7:nonstandard-pointDensity>
<ns7:maximumStepSize>0.001</ns7:maximumStepSize>
</rdf:Description>
<rdf:Description rdf:about="rdf:#43a014e2-42c1-44fe-82ac-55a21c5f0ed8">
<dcterms:modified rdf:resource="rdf:#70008440-24f7-4d2c-a98f-b41499e81135"/>
<rdf:value>This model has been curated by Penny Noble of Oxford University</rdf:value>
<cmeta:modifier rdf:resource="rdf:#f41e4fd4-bf58-44e8-8f42-bf58bf1aae93"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8b76b0bc-1b43-4263-bef7-9ccda50484ab">
<vCard:FN>Unknown Person</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6691363a-783d-4a5a-a783-51c9684f1fb1">
<dcterms:W3CDTF>2001-12-00 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2f0fc38a-3c1c-4bfb-a830-47c3131b56c3">
<dcterms:modified rdf:resource="rdf:#1b9f3170-1975-4598-acb3-6c6499b9c9dc"/>
<cmeta:modifier rdf:resource="rdf:#f45fe86a-a280-460b-a9a9-ff3aedc56221"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4e120d40-0111-469a-8f20-e5bfb5ede2e4">
<dcterms:modified rdf:resource="rdf:#4547f559-c84e-4c50-a77b-57e61d1ca6c7"/>
<rdf:value>Added some cmeta:id's to current and flux variables</rdf:value>
<cmeta:modifier rdf:resource="rdf:#d79d8c97-5b38-4c39-9e96-0aad2089274a"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#359aa759-ca02-411f-ae3f-3da81bc0f064">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#1c7efc4b-1d4c-4068-a851-0d7610c1e695"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4377fe6c-add1-493b-84bb-f38bb0c075c7">
<dcterms:modified rdf:resource="rdf:#343cfb12-4cd6-4689-9fe9-ae348c1f9e50"/>
<rdf:value>Summary of code changes from v07:
1. The inward rectifying current (Ik1) had a 0.0001 term in the numerator which has been changed to 0.001 to match the paper and corrections document.
2. Mmax, hmax, jmax, dmax, f11max, f12max, Caactmax, rmax, smax sslowmax, rssmax and sssmax were all called minfinity, hinfinity, jinfinity, dinfinity...etc. These have been changed to ...max variables for consistency internally and with the paper.
3. The Ryanodine equations had odd dimensionless units and terms which divide by 1 to make the units work. Here variables 'n' and 'm' were removed and the units set as per the original model of Keizer, J. and Levine, L. Ryanodine receptor adaptation and Ca2+(-)induced Ca2+ release-dependent Ca2+ oscillations. Biophys. J. 71,6, 3477-3487 (1996). This was after discussion with both Catherine Lloyd and Penny Noble. The units are now correct.
4. Units millimolar_4 was misnamed and is now per_millimolar4.
5. Units nanoA_millimolar_4 was misnamed and is now nanoA_per_millimolar4.
6. Not a change but just to note that this implementation has a 1Hz stimulus for 10 seconds built in.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#a89fcb96-f39b-4075-847b-e46f395266af"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c9681df6-408f-4fae-ae6f-5a3ff09c0475">
<dcterms:W3CDTF>2007-01-16T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#68041856-fc2b-41b9-b67e-434a1f2af75f">
<vCard:Given>James</vCard:Given>
<vCard:Family>Lawson</vCard:Family>
<vCard:Other>Richard</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#847270ad-ac16-478c-ba30-577f1e0483fe">
<dcterms:modified rdf:resource="rdf:#d33019dc-19a2-4363-9795-a23e24743685"/>
<rdf:value>Several variables were given cmeta:id's to allow creation of a PCEnv session file.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#5701be65-dbe4-4143-aa31-f1fc094c2d2c"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#70008440-24f7-4d2c-a98f-b41499e81135">
<dcterms:W3CDTF>2007-01-16T14:20:00+13:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f45fe86a-a280-460b-a9a9-ff3aedc56221">
<vCard:N rdf:resource="rdf:#dfcc85e4-f2cd-44b9-81f0-c1a1b38e1193"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#aa39cce2-2692-4f0c-9409-68cd0b6580f7">
<dc:creator rdf:resource="rdf:#8b76b0bc-1b43-4263-bef7-9ccda50484ab"/>
<rdf:value>This version was created by Penny Noble of Oxford University and represents the EPICARDIAL CELL VARIANT. The model has been checked in COR.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#338a6b4b-01d1-4ad6-9f4e-4bd24d16edd1">
<dcterms:W3CDTF>2008-07-01T16:01:54+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d33019dc-19a2-4363-9795-a23e24743685">
<dcterms:W3CDTF>2007-06-22T14:17:56+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#$RvprG">
<rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
<rdf:first rdf:resource="rdf:#$SvprG"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5701be65-dbe4-4143-aa31-f1fc094c2d2c">
<vCard:N rdf:resource="rdf:#2ecf70be-c03a-409e-ab2d-0f76aa6536c5"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6cf4bb45-eef1-4cad-965f-93ee42368494">
<dcterms:modified rdf:resource="rdf:#2ee4b150-2e99-4617-8c15-fb11b9a9cb12"/>
<rdf:value>Added two differential equations to calculate HTRPNCa and LTRPNCa (the last two equations on page 3047 of the original paper). In doing so fixed the model to recreate the published results. Frank Sachse has checked this model in JSim and it also runs in this tool.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#91960679-9437-40bf-9215-85f10e927a1b"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d13ad4d7-2b77-45c4-a989-9a2c3ccc204e">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#abdd93f9-81ec-4d28-b70d-a0ae4eec7a39"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a00fa64d-2c70-4560-8de5-182db60b50ae">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7439b4d0-6082-44e6-ad00-3c6a36c12aec">
<vCard:Given>S</vCard:Given>
<vCard:Family>Pandit</vCard:Family>
<vCard:Other>V</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#27c55eb9-83e5-42a1-be0e-6ce759afb405">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#59cb9d3e-3dca-45f1-9857-f9e5a07de4f0">
<dc:creator rdf:resource="rdf:#010abb91-8c98-4c4b-9125-9611c9229f62"/>
<rdf:value/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4547f559-c84e-4c50-a77b-57e61d1ca6c7">
<dcterms:W3CDTF>2008-07-31T10:37:18+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c4600241-ec5b-47ee-8294-30cb2854a60e">
<vCard:ORG rdf:resource="rdf:#6324c2a1-51e4-4a74-8693-1021f04bc9b1"/>
<vCard:EMAIL rdf:resource="rdf:#48f139cb-d71c-4c43-b9fc-30a96302b9f5"/>
<vCard:N rdf:resource="rdf:#27c55eb9-83e5-42a1-be0e-6ce759afb405"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#952a6e88-d239-4cd5-8db7-a0438bcc0557">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#7439b4d0-6082-44e6-ad00-3c6a36c12aec"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d213868e-80c0-4c8a-bd5b-9e7f68be57fe">
<dc:creator rdf:resource="rdf:#cc3ea2c4-6178-48e5-af65-7906d40d1fa4"/>
<dc:title>A Mathematical Model of Action Potential Heterogeneity in Adult Rat Left Ventricular Myocytes</dc:title>
<bqs:volume>81</bqs:volume>
<bqs:first_page>3029</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#9dd15c9e-7d73-41d7-a27f-5a4585bbdad1"/>
<dcterms:issued rdf:resource="rdf:#6691363a-783d-4a5a-a783-51c9684f1fb1"/>
<bqs:last_page>3051</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3801bb39-63d2-4ea7-b5a5-26d159ddbb04">
<vCard:N rdf:resource="rdf:#68041856-fc2b-41b9-b67e-434a1f2af75f"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c2b4910b-52ce-4333-9f1b-d570efb93892">
<vCard:Given>Mike</vCard:Given>
<vCard:Family>Cooling</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1a06d1a1-b9f2-47f8-a61a-fac50ece9828">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c5725b6b-4fe6-4323-b8b5-c46dbd8b8635">
<bqs:Pubmed_id>11720973</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#7a7ac380-fabd-4cc9-b603-061dce055eea"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7c633a24-2630-4ed1-87fc-b33b8f35b18f">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#8a9f2967-4170-417e-9653-efdad57d4d0c"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d71f5fd8-902a-4360-b367-ee3d45f19635">
<dcterms:W3CDTF>2001-12-01 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#560f4835-5400-4f6c-8379-e837c2f3987d">
<dcterms:W3CDTF>2008-06-13T13:56:13+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e2d39fc6-49fd-4c39-93ec-3a3519b1d8ff">
<dc:title>Biophysical Journal</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#71546bf6-ff74-4c92-8272-878e3c7ff18b">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#9996133d-52a8-465f-8c28-1da0b94d2322"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2ee4b150-2e99-4617-8c15-fb11b9a9cb12">
<dcterms:W3CDTF>2007-11-01T12:18:01+13:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1c7efc4b-1d4c-4068-a851-0d7610c1e695">
<vCard:Given>W</vCard:Given>
<vCard:Family>Giles</vCard:Family>
<vCard:Other>R</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#136ced93-18a7-4f9b-9628-1b29057940bc">
<dcterms:modified rdf:resource="rdf:#560f4835-5400-4f6c-8379-e837c2f3987d"/>
<rdf:value>Noticed the volumes for V_NSR and V_SS were wrong and have corrected them. They were wrong in version 7.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#589a05f4-5e47-4059-90d3-e577f2fb8755"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#abdd93f9-81ec-4d28-b70d-a0ae4eec7a39">
<vCard:Given>R</vCard:Given>
<vCard:Family>Clark</vCard:Family>
<vCard:Other>B</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7702b560-a156-4e0a-b769-51c3ae395590">
<vCard:Given>S</vCard:Given>
<vCard:Family>Pandit</vCard:Family>
<vCard:Other>V</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#343cfb12-4cd6-4689-9fe9-ae348c1f9e50">
<dcterms:W3CDTF>2008-05-26T11:07:32+12:00</dcterms:W3CDTF>
</rdf:Description>
</rdf:RDF>
</model>
