<?xml version='1.0' encoding='utf-8'?>
<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#" cmeta:id="demir_clark_giles_murphey_1994_version01" name="demir_clark_giles_murphey_1994_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>A Mathematical Model Of A Rabbit Sinoatrial Node Cell, 1994</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 version of the model has been curated by Penny Noble and is known to run in COR and PCEnv to replicate the published results. Units have been checked and they are balanced.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: A mathematical model for the electrophysiological responses of a rabbit sinoatrial node cell that is based on whole cell recordings from enzymatically isolated single pacemaker cells at 37 degrees C has been developed. The ion channels, Na(+)-K+ and Ca2+ pumps, and Na(+)-Ca2+ exchanger in the surface membrane (sarcolemma) are described using equations for these known currents in mammalian pacemaker cells. The extracellular environment is treated as a diffusion-limited space, and the myoplasm contains Ca(2+)-binding proteins (calmodulin and troponin). Original features of this model include 1) new equations for the hyperpolarization-activated inward current, 2) assessment of the role of the transient-type Ca2+ current during pacemaker depolarization, 3) inclusion of an Na+ current based on recent experimental data, and 4) demonstration of the possible influence of pump and exchanger currents and background currents on the pacemaker rate. This model provides acceptable fits to voltage-clamp and action potential data and can be used to seek biophysically based explanations of the electrophysiological activity in the rabbit sinoatrial node cell.
</para>
<para>
The original paper reference is cited below:
</para>
<para>
A mathematical model of a rabbit sinoatrial node cell, S.S. Demir, J.W. Clark, C.R. Murphey and W.R. Giles, 1994,
<emphasis>American Journal of Physiology</emphasis>, 266, C832-C852. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=8166247&dopt=Abstract">PubMed ID: 8166247</ulink>
</para>
<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>cell diagram of the Demir et al SAN model showing ionic currents, pumps and exchangers within the sarcolemma</title>
</objectinfo>
<imagedata fileref="demir_1994.png"/>
</imageobject>
</mediaobject>
<caption>A schematic diagram describing the current flows across the cell membrane that are captured in the Demir <emphasis>et al</emphasis> 1994 model of the electrophysiological activity in a SA node 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_millivolt_second">
<unit units="millivolt" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="microS">
<unit units="siemens" prefix="micro"/>
</units>
<units name="microF">
<unit units="farad" prefix="micro"/>
</units>
<units name="nanoA">
<unit units="ampere" prefix="nano"/>
</units>
<units name="millimolar">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="per_millimolar">
<unit units="millimolar" exponent="-1"/>
</units>
<units name="per_millimolar_second">
<unit units="millimolar" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="millimolar_per_second">
<unit units="millimolar"/>
<unit units="second" exponent="-1"/>
</units>
<units name="microS_per_millimolar">
<unit units="microS"/>
<unit units="millimolar" exponent="-1"/>
</units>
<units name="nanoA_per_millimolar">
<unit units="nanoA"/>
<unit units="millimolar" exponent="-1"/>
</units>
<units name="millimolar4">
<unit units="millimolar" exponent="4"/>
</units>
<units name="joule_per_kilomole_kelvin">
<unit units="joule"/>
<unit units="mole" prefix="kilo" exponent="-1"/>
<unit units="kelvin" exponent="-1"/>
</units>
<units name="coulomb_per_mole">
<unit units="coulomb"/>
<unit units="mole" exponent="-1"/>
</units>
<units name="mul_per_second">
<unit units="litre" prefix="micro"/>
<unit units="second" exponent="-1"/>
</units>
<units name="microLitre">
<unit units="litre" prefix="micro"/>
</units>
<component name="environment">
<variable name="time" units="second" public_interface="out"/>
</component>
<component name="membrane">
<variable name="V" units="millivolt" initial_value="-49.54105" public_interface="out"/>
<variable name="R" units="joule_per_kilomole_kelvin" initial_value="8314.472" public_interface="out"/>
<variable name="T" units="kelvin" initial_value="310" public_interface="out"/>
<variable name="F" units="coulomb_per_mole" initial_value="96485.3415" public_interface="out"/>
<variable name="Cm" units="microF" initial_value="5.5e-5"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="i_Na" units="nanoA" public_interface="in"/>
<variable name="i_Ca_T" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L" units="nanoA" public_interface="in"/>
<variable name="i_K" units="nanoA" public_interface="in"/>
<variable name="i_f" units="nanoA" public_interface="in"/>
<variable name="i_B" units="nanoA" public_interface="in"/>
<variable name="i_NaK" units="nanoA" public_interface="in"/>
<variable name="i_NaCa" units="nanoA" public_interface="in"/>
<variable name="i_Ca_P" units="nanoA" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<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_T</ci>
<ci>i_Ca_L</ci>
<ci>i_K</ci>
<ci>i_f</ci>
<ci>i_B</ci>
<ci>i_NaK</ci>
<ci>i_NaCa</ci>
<ci>i_Ca_P</ci>
</apply>
</apply>
<ci>Cm</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_current">
<variable name="i_Na" units="nanoA" public_interface="out"/>
<variable name="P_Na" units="mul_per_second" initial_value="0.00344"/>
<variable name="E_Na" units="millivolt" public_interface="in"/>
<variable name="Na_c" units="millimolar" public_interface="in"/>
<variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="m" units="dimensionless" private_interface="in"/>
<variable name="h1" units="dimensionless" private_interface="in"/>
<variable name="h2" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Na</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>P_Na</ci>
<apply>
<power/>
<ci>m</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>h1</ci>
<ci>h2</ci>
<ci>Na_c</ci>
<ci>V</ci>
<apply>
<power/>
<ci>F</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_m_gate">
<variable name="m" units="dimensionless" initial_value="0.250113" public_interface="out"/>
<variable name="m_infinity" units="dimensionless"/>
<variable name="tau_m" units="second"/>
<variable name="alpha_m" units="per_second"/>
<variable name="beta_m" units="per_second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_m</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt_second">824</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">51.9</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">51.9</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">8.9</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_m</ci>
<apply>
<times/>
<cn cellml:units="per_second">32960</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">51.9</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">8.9</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>m_infinity</ci>
<apply>
<divide/>
<ci>alpha_m</ci>
<apply>
<plus/>
<ci>alpha_m</ci>
<ci>beta_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_m</ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_m</ci>
<ci>beta_m</ci>
</apply>
</apply>
<cn cellml:units="second">0.000015</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>m</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>m_infinity</ci>
<ci>m</ci>
</apply>
<ci>tau_m</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_h_gate">
<variable name="h1" units="dimensionless" initial_value="0.001386897" public_interface="out"/>
<variable name="h2" units="dimensionless" initial_value="0.002065463" public_interface="out"/>
<variable name="h1_infinity" units="dimensionless"/>
<variable name="h2_infinity" units="dimensionless"/>
<variable name="tau_h1" units="second"/>
<variable name="tau_h2" units="second"/>
<variable name="alpha_h1" units="per_second"/>
<variable name="beta_h1" units="per_second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_h1</ci>
<apply>
<times/>
<cn cellml:units="per_second">165</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">101.3</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">12.6</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_h1</ci>
<apply>
<divide/>
<cn cellml:units="per_second">12360</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">320</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">101.3</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">12.6</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>h1_infinity</ci>
<apply>
<divide/>
<ci>alpha_h1</ci>
<apply>
<plus/>
<ci>alpha_h1</ci>
<ci>beta_h1</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_h1</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_h1</ci>
<ci>beta_h1</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>h2_infinity</ci>
<ci>h1_infinity</ci>
</apply>
<apply>
<eq/>
<ci>tau_h2</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">20</cn>
<ci>tau_h1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>h1</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>h1_infinity</ci>
<ci>h1</ci>
</apply>
<ci>tau_h1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>h2</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>h2_infinity</ci>
<ci>h2</ci>
</apply>
<ci>tau_h2</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel">
<variable name="i_Ca_L" units="nanoA" public_interface="out"/>
<variable name="g_Ca_L" units="microS" initial_value="0.02115"/>
<variable name="E_Ca_L" units="millivolt" initial_value="46.4"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="d_L" units="dimensionless" private_interface="in"/>
<variable name="d_L_infinity" units="dimensionless" private_interface="in"/>
<variable name="f_L" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ca_L</ci>
<apply>
<times/>
<ci>g_Ca_L</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_L</ci>
<ci>d_L</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.095</cn>
<ci>d_L_infinity</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 name="d_L" units="dimensionless" initial_value="0.002572773" public_interface="out"/>
<variable name="d_L_infinity" units="dimensionless" public_interface="out"/>
<variable name="alpha_d_L" units="per_second"/>
<variable name="beta_d_L" units="per_second"/>
<variable name="tau_d_L" units="second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_d_L</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt_second">28.39</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">35</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">35</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">2.5</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt_second">84.9</cn>
</apply>
<ci>V</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.208</cn>
</apply>
<ci>V</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_d_L</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">11.43</cn>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.4</cn>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_d_L</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_d_L</ci>
<ci>beta_d_L</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>d_L_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">14.1</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">6</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>d_L</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>d_L_infinity</ci>
<ci>d_L</ci>
</apply>
<ci>tau_d_L</ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_f_gate">
<variable name="f_L" units="dimensionless" initial_value="0.98651" public_interface="out"/>
<variable name="alpha_f_L" units="per_second"/>
<variable name="beta_f_L" units="per_second"/>
<variable name="f_L_infinity" units="dimensionless"/>
<variable name="tau_f_L" units="second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_f_L</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">3.75</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">28</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">28</cn>
</apply>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_f_L</ci>
<apply>
<divide/>
<cn cellml:units="per_second">30</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">28</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_f_L</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_f_L</ci>
<ci>beta_f_L</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>f_L_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">30</cn>
</apply>
<cn cellml:units="millivolt">5</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>f_L</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>f_L_infinity</ci>
<ci>f_L</ci>
</apply>
<ci>tau_f_L</ci>
</apply>
</apply>
</math>
</component>
<component name="T_type_Ca_channel">
<variable name="i_Ca_T" units="nanoA" public_interface="out"/>
<variable name="g_Ca_T" units="microS" initial_value="0.02521"/>
<variable name="E_Ca_T" units="millivolt" initial_value="45"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="d_T" units="dimensionless" private_interface="in"/>
<variable name="f_T" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Ca_T</ci>
<apply>
<times/>
<ci>g_Ca_T</ci>
<ci>d_T</ci>
<ci>f_T</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Ca_T</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="T_type_Ca_channel_d_gate">
<variable name="d_T" units="dimensionless" initial_value="0.02012114" public_interface="out"/>
<variable name="alpha_d_T" units="per_second"/>
<variable name="beta_d_T" units="per_second"/>
<variable name="d_T_infinity" units="dimensionless"/>
<variable name="tau_d_T" units="second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_d_T</ci>
<apply>
<times/>
<cn cellml:units="per_second">1068</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">26.3</cn>
</apply>
<cn cellml:units="millivolt">30</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_d_T</ci>
<apply>
<times/>
<cn cellml:units="per_second">1068</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">26.3</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">30</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_d_T</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_d_T</ci>
<ci>beta_d_T</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>d_T_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">26.3</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">6</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>d_T</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>d_T_infinity</ci>
<ci>d_T</ci>
</apply>
<ci>tau_d_T</ci>
</apply>
</apply>
</math>
</component>
<component name="T_type_Ca_channel_f_gate">
<variable name="f_T" units="dimensionless" initial_value="0.1945111" public_interface="out"/>
<variable name="alpha_f_T" units="per_second"/>
<variable name="beta_f_T" units="per_second"/>
<variable name="f_T_infinity" units="dimensionless"/>
<variable name="tau_f_T" units="second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_f_T</ci>
<apply>
<times/>
<cn cellml:units="per_second">15.3</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">61.7</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">83.3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_f_T</ci>
<apply>
<times/>
<cn cellml:units="per_second">15</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">61.7</cn>
</apply>
<cn cellml:units="millivolt">15.38</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_f_T</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>alpha_f_T</ci>
<ci>beta_f_T</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>f_T_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">61.7</cn>
</apply>
<cn cellml:units="millivolt">5.6</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>f_T</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>f_T_infinity</ci>
<ci>f_T</ci>
</apply>
<ci>tau_f_T</ci>
</apply>
</apply>
</math>
</component>
<component name="delayed_rectifying_potassium_current">
<variable name="i_K" units="nanoA" public_interface="out"/>
<variable name="g_K" units="microS"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="K_b" units="millimolar" public_interface="in"/>
<variable name="P_a" units="dimensionless" private_interface="in"/>
<variable name="P_i" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>g_K</ci>
<apply>
<times/>
<cn cellml:units="microS">0.00693</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>K_b</ci>
<cn cellml:units="millimolar">1</cn>
</apply>
<cn cellml:units="dimensionless">0.59</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_K</ci>
<apply>
<times/>
<ci>g_K</ci>
<ci>P_a</ci>
<ci>P_i</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="delayed_rectifying_potassium_current_P_a_gate">
<variable name="P_a" units="dimensionless" initial_value="0.02302278" public_interface="out"/>
<variable name="tau_P_a" units="second"/>
<variable name="P_a_infinity" units="dimensionless"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>P_a_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">5.1</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">7.4</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_P_a</ci>
<apply>
<divide/>
<cn cellml:units="second">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">17</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.0398</cn>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2.11</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.051</cn>
</apply>
<ci>V</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P_a</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>P_a_infinity</ci>
<ci>P_a</ci>
</apply>
<ci>tau_P_a</ci>
</apply>
</apply>
</math>
</component>
<component name="delayed_rectifying_potassium_current_P_i_gate">
<variable name="P_i" units="dimensionless" initial_value="0.3777728" public_interface="out"/>
<variable name="alpha_P_i" units="per_second"/>
<variable name="beta_P_i" units="per_second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_P_i</ci>
<apply>
<times/>
<cn cellml:units="per_second">100</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.0183</cn>
</apply>
<ci>V</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_P_i</ci>
<apply>
<times/>
<cn cellml:units="per_second">656</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.00942</cn>
<ci>V</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P_i</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_P_i</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_i</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_P_i</ci>
<ci>P_i</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="linear_background_current">
<variable name="i_B" units="nanoA" public_interface="out"/>
<variable name="i_B_Na" units="nanoA" public_interface="out"/>
<variable name="i_B_Ca" units="nanoA" public_interface="out"/>
<variable name="i_B_K" units="nanoA" public_interface="out"/>
<variable name="g_B_Na" units="microS" initial_value="0.00016"/>
<variable name="g_B_Ca" units="microS" initial_value="0.0000364"/>
<variable name="g_B_K" units="microS" initial_value="0.0000694"/>
<variable name="E_Na" units="millivolt" public_interface="in"/>
<variable name="E_Ca" units="millivolt" public_interface="in"/>
<variable name="E_K" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<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="hyperpolarisation_activated_current">
<variable name="i_f" units="nanoA" public_interface="out"/>
<variable name="i_f_Na" units="nanoA" public_interface="out"/>
<variable name="i_f_K" units="nanoA" public_interface="out"/>
<variable name="g_f_Na" units="microS" initial_value="0.0067478"/>
<variable name="g_f_K" units="microS" initial_value="0.0128821"/>
<variable name="time" units="second" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="y" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_f_Na</ci>
<apply>
<times/>
<ci>g_f_Na</ci>
<apply>
<power/>
<ci>y</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">75</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_f_K</ci>
<apply>
<times/>
<ci>g_f_K</ci>
<apply>
<power/>
<ci>y</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">85</cn>
</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 name="y" units="dimensionless" initial_value="0.09227776" public_interface="out"/>
<variable name="y_infinity" units="dimensionless"/>
<variable name="tau_y" units="second"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<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">72.2</cn>
</apply>
<cn cellml:units="millivolt">9</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">1.6483</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">54.06</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">24.33</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">14.01055</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">0.7</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">60</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">5.5</cn>
</apply>
</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="sodium_potassium_pump">
<variable name="i_NaK" units="nanoA" public_interface="out"/>
<variable name="K_m_Na" units="millimolar" initial_value="5.46"/>
<variable name="K_m_K" units="millimolar" initial_value="0.621"/>
<variable name="i_NaK_max" units="nanoA" initial_value="0.2192"/>
<variable name="Na_i" units="millimolar" public_interface="in"/>
<variable name="K_c" units="millimolar" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_NaK</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>i_NaK_max</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>Na_i</ci>
<apply>
<plus/>
<ci>K_m_Na</ci>
<ci>Na_i</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">3</cn>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>K_c</ci>
<apply>
<plus/>
<ci>K_m_K</ci>
<ci>K_c</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<cn cellml:units="dimensionless">1.6</cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1.5</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">60</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_pump_current">
<variable name="i_Ca_P" units="nanoA" public_interface="out"/>
<variable name="i_Ca_P_max" units="nanoA" initial_value="0.02869"/>
<variable name="Ca_i" units="millimolar" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<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="sodium_calcium_pump">
<variable name="i_NaCa" units="nanoA" public_interface="out"/>
<variable name="K_NaCa" units="nanoA" initial_value="0.00001248"/>
<variable name="d_NaCa" units="dimensionless" initial_value="0.0001"/>
<variable name="gamma" units="dimensionless" initial_value="0.5"/>
<variable name="Na_i" units="millimolar" public_interface="in"/>
<variable name="Na_c" units="millimolar" public_interface="in"/>
<variable name="Ca_i" units="millimolar" public_interface="in"/>
<variable name="Ca_c" units="millimolar" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<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_c</ci>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.03743</cn>
<ci>V</ci>
<ci>gamma</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>Na_c</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</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="millimolar4">1</cn>
<apply>
<times/>
<ci>d_NaCa</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>Ca_i</ci>
<apply>
<power/>
<ci>Na_c</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_c</ci>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="intracellular_concentrations_and_buffer_equations">
<variable name="Na_i" units="millimolar" initial_value="9.701621" public_interface="out"/>
<variable name="Ca_i" units="millimolar" initial_value="3.787018e-4" public_interface="out"/>
<variable name="K_i" units="millimolar" initial_value="1.407347e2" public_interface="out"/>
<variable name="Ca_Calmod" units="dimensionless" initial_value="0.1411678"/>
<variable name="Ca_Trop" units="dimensionless" initial_value="0.07331396"/>
<variable name="Ca_Mg_Trop" units="dimensionless" initial_value="0.7618549"/>
<variable name="Mg_Mg_Trop" units="dimensionless" initial_value="0.2097049"/>
<variable name="phi_C" units="per_second"/>
<variable name="phi_TC" units="per_second"/>
<variable name="phi_TMgC" units="per_second"/>
<variable name="phi_TMgM" units="per_second"/>
<variable name="phi_B" units="millimolar_per_second"/>
<variable name="Mg_i" units="millimolar" initial_value="2.5"/>
<variable name="F_C" units="millimolar_per_second"/>
<variable name="F_TC" units="millimolar_per_second"/>
<variable name="F_TMgC" units="millimolar_per_second"/>
<variable name="Vol" units="microLitre" public_interface="in"/>
<variable name="V_i" units="microLitre" public_interface="out"/>
<variable name="i_NaK" units="nanoA" public_interface="in"/>
<variable name="i_NaCa" units="nanoA" public_interface="in"/>
<variable name="i_B_Na" units="nanoA" public_interface="in"/>
<variable name="i_f_Na" units="nanoA" public_interface="in"/>
<variable name="i_Na" units="nanoA" public_interface="in"/>
<variable name="i_K" units="nanoA" public_interface="in"/>
<variable name="i_B_K" units="nanoA" public_interface="in"/>
<variable name="i_f_K" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L" units="nanoA" public_interface="in"/>
<variable name="i_Ca_T" units="nanoA" public_interface="in"/>
<variable name="i_Ca_P" units="nanoA" public_interface="in"/>
<variable name="i_B_Ca" units="nanoA" public_interface="in"/>
<variable name="i_up" units="nanoA" public_interface="in"/>
<variable name="i_rel" units="nanoA" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_i</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.465</cn>
<ci>Vol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>phi_C</ci>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_millimolar_second">129000</cn>
<ci>Ca_i</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>Ca_Calmod</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">307</cn>
<ci>Ca_Calmod</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Calmod</ci>
</apply>
<ci>phi_C</ci>
</apply>
<apply>
<eq/>
<ci>phi_TC</ci>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_millimolar_second">50500</cn>
<ci>Ca_i</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>Ca_Trop</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">252</cn>
<ci>Ca_Trop</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Trop</ci>
</apply>
<ci>phi_TC</ci>
</apply>
<apply>
<eq/>
<ci>phi_TMgC</ci>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_millimolar_second">129000</cn>
<ci>Ca_i</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>Ca_Mg_Trop</ci>
<ci>Mg_Mg_Trop</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">4.25</cn>
<ci>Ca_Mg_Trop</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Mg_Trop</ci>
</apply>
<ci>phi_TMgC</ci>
</apply>
<apply>
<eq/>
<ci>phi_TMgM</ci>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_millimolar_second">1290</cn>
<ci>Mg_i</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>Ca_Mg_Trop</ci>
<ci>Mg_Mg_Trop</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">429</cn>
<ci>Mg_Mg_Trop</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Mg_Mg_Trop</ci>
</apply>
<ci>phi_TMgM</ci>
</apply>
<apply>
<eq/>
<ci>F_C</ci>
<apply>
<times/>
<cn cellml:units="millimolar">0.09</cn>
<ci>phi_C</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>F_TC</ci>
<apply>
<times/>
<cn cellml:units="millimolar">0.031</cn>
<ci>phi_TC</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>F_TMgC</ci>
<apply>
<times/>
<cn cellml:units="millimolar">0.062</cn>
<ci>phi_TMgC</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>phi_B</ci>
<apply>
<plus/>
<ci>F_C</ci>
<ci>F_TC</ci>
<ci>F_TMgC</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Na_i</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>i_NaK</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>i_NaCa</ci>
</apply>
<ci>i_B_Na</ci>
<ci>i_f_Na</ci>
<ci>i_Na</ci>
</apply>
</apply>
<apply>
<times/>
<ci>F</ci>
<ci>V_i</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_i</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>i_NaK</ci>
</apply>
<apply>
<plus/>
<ci>i_K</ci>
<ci>i_f_K</ci>
<ci>i_B_K</ci>
</apply>
</apply>
<apply>
<times/>
<ci>F</ci>
<ci>V_i</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_i</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>i_NaCa</ci>
</apply>
<ci>i_rel</ci>
</apply>
<apply>
<plus/>
<ci>i_Ca_L</ci>
<ci>i_Ca_T</ci>
<ci>i_Ca_P</ci>
<ci>i_B_Ca</ci>
<ci>i_up</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>V_i</ci>
<ci>F</ci>
</apply>
</apply>
<ci>phi_B</ci>
</apply>
</apply>
</math>
</component>
<component name="cleft_space_equations">
<variable name="Na_c" units="millimolar" initial_value="139.9988" public_interface="out"/>
<variable name="Ca_c" units="millimolar" initial_value="2.00474" public_interface="out"/>
<variable name="K_c" units="millimolar" initial_value="5.389014" public_interface="out"/>
<variable name="K_b" units="millimolar" initial_value="5.4" public_interface="out"/>
<variable name="Na_b" units="millimolar" initial_value="140"/>
<variable name="Ca_b" units="millimolar" initial_value="2"/>
<variable name="Vol" units="microLitre" initial_value="3.497e-6" public_interface="out"/>
<variable name="V_c" units="microLitre"/>
<variable name="tau_p" units="second" initial_value="0.01"/>
<variable name="i_NaK" units="nanoA" public_interface="in"/>
<variable name="i_NaCa" units="nanoA" public_interface="in"/>
<variable name="i_B_Na" units="nanoA" public_interface="in"/>
<variable name="i_f_Na" units="nanoA" public_interface="in"/>
<variable name="i_Na" units="nanoA" public_interface="in"/>
<variable name="i_K" units="nanoA" public_interface="in"/>
<variable name="i_B_K" units="nanoA" public_interface="in"/>
<variable name="i_f_K" units="nanoA" public_interface="in"/>
<variable name="i_Ca_L" units="nanoA" public_interface="in"/>
<variable name="i_Ca_T" units="nanoA" public_interface="in"/>
<variable name="i_Ca_P" units="nanoA" public_interface="in"/>
<variable name="i_B_Ca" units="nanoA" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_c</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.136</cn>
<ci>Vol</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Na_c</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Na_b</ci>
<ci>Na_c</ci>
</apply>
<ci>tau_p</ci>
</apply>
<apply>
<divide/>
<apply>
<plus/>
<ci>i_Na</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>i_NaCa</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>i_NaK</ci>
</apply>
<ci>i_B_Na</ci>
<ci>i_f_Na</ci>
</apply>
<apply>
<times/>
<ci>F</ci>
<ci>V_c</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_c</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<ci>K_b</ci>
<ci>K_c</ci>
</apply>
<ci>tau_p</ci>
</apply>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>i_NaK</ci>
</apply>
<ci>i_K</ci>
<ci>i_B_K</ci>
<ci>i_f_K</ci>
</apply>
<apply>
<times/>
<ci>F</ci>
<ci>V_c</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_c</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Ca_b</ci>
<ci>Ca_c</ci>
</apply>
<ci>tau_p</ci>
</apply>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>i_NaCa</ci>
</apply>
<ci>i_Ca_L</ci>
<ci>i_Ca_T</ci>
<ci>i_Ca_P</ci>
<ci>i_B_Ca</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
<ci>V_c</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SR_Ca_uptake_and_release">
<variable name="i_rel" units="nanoA" public_interface="out"/>
<variable name="i_up" units="nanoA" public_interface="out"/>
<variable name="i_tr" units="nanoA"/>
<variable name="Ca_up" units="millimolar" initial_value="16.95311"/>
<variable name="Ca_rel" units="millimolar" initial_value="16.85024"/>
<variable name="Ca_i" units="millimolar" public_interface="in"/>
<variable name="alpha_up" units="nanoA" initial_value="0.08"/>
<variable name="beta_up" units="nanoA" initial_value="0.072"/>
<variable name="alpha_rel" units="nanoA_per_millimolar" initial_value="0.5"/>
<variable name="K1" units="dimensionless"/>
<variable name="K2" units="millimolar"/>
<variable name="k_cyca" units="millimolar" initial_value="0.00005"/>
<variable name="k_xcs" units="dimensionless" initial_value="0.9"/>
<variable name="k_SRCa" units="millimolar" initial_value="22"/>
<variable name="k_rel" units="millimolar" initial_value="0.004"/>
<variable name="r_act" units="per_second"/>
<variable name="r_inact" units="per_second"/>
<variable name="Ca_Calse" units="dimensionless" initial_value="0.9528726"/>
<variable name="phi_Calse" units="per_second"/>
<variable name="F1" units="dimensionless" initial_value="0.1133251"/>
<variable name="F2" units="dimensionless" initial_value="0.0007594214"/>
<variable name="F3" units="dimensionless" initial_value="0.8859153"/>
<variable name="V_i" units="microLitre" public_interface="in"/>
<variable name="V_up" units="microLitre"/>
<variable name="V_rel" units="microLitre"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_up</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.01166</cn>
<ci>V_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_rel</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">0.001296</cn>
<ci>V_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>r_act</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="per_second">240</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="millivolt">40</cn>
</apply>
<cn cellml:units="per_millivolt">0.08</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">240</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Ca_i</ci>
<apply>
<plus/>
<ci>Ca_i</ci>
<ci>k_rel</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r_inact</ci>
<apply>
<plus/>
<cn cellml:units="per_second">40</cn>
<apply>
<times/>
<cn cellml:units="per_second">240</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>Ca_i</ci>
<apply>
<plus/>
<ci>Ca_i</ci>
<ci>k_rel</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>phi_Calse</ci>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_millimolar_second">770</cn>
<ci>Ca_rel</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>Ca_Calse</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">641</cn>
<ci>Ca_Calse</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_Calse</ci>
</apply>
<ci>phi_Calse</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>F1</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_second">0.96</cn>
<ci>F3</ci>
</apply>
<apply>
<times/>
<ci>r_act</ci>
<ci>F1</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>F2</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>r_act</ci>
<ci>F1</ci>
</apply>
<apply>
<times/>
<ci>r_inact</ci>
<ci>F2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>F3</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>r_inact</ci>
<ci>F2</ci>
</apply>
<apply>
<times/>
<cn cellml:units="per_second">0.96</cn>
<ci>F3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>K1</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>k_cyca</ci>
<ci>k_xcs</ci>
</apply>
<ci>k_SRCa</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>K2</ci>
<apply>
<plus/>
<ci>Ca_i</ci>
<apply>
<times/>
<ci>Ca_up</ci>
<ci>K1</ci>
</apply>
<apply>
<times/>
<ci>k_cyca</ci>
<ci>k_xcs</ci>
</apply>
<ci>k_cyca</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_up</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_up</ci>
<ci>Ca_i</ci>
</apply>
<apply>
<times/>
<ci>beta_up</ci>
<ci>Ca_up</ci>
<ci>K1</ci>
</apply>
</apply>
<ci>K2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_rel</ci>
<apply>
<times/>
<ci>alpha_rel</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>F2</ci>
<apply>
<plus/>
<ci>F2</ci>
<cn cellml:units="dimensionless">0.25</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>Ca_rel</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>i_tr</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>Ca_up</ci>
<ci>Ca_rel</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
<ci>V_up</ci>
</apply>
<cn cellml:units="second">0.06418</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_up</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>i_up</ci>
<ci>i_tr</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>V_up</ci>
<ci>F</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_rel</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<minus/>
<ci>i_tr</ci>
<ci>i_rel</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>V_rel</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millimolar">11.48</cn>
<ci>phi_Calse</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="reversal_potentials">
<variable name="E_Na" units="millivolt" public_interface="out"/>
<variable name="E_K" units="millivolt" public_interface="out"/>
<variable name="E_Ca" units="millivolt" public_interface="out"/>
<variable name="K_c" units="millimolar" public_interface="in"/>
<variable name="Na_c" units="millimolar" public_interface="in"/>
<variable name="K_i" units="millimolar" public_interface="in"/>
<variable name="Na_i" units="millimolar" public_interface="in"/>
<variable name="Ca_c" units="millimolar" public_interface="in"/>
<variable name="Ca_i" units="millimolar" public_interface="in"/>
<variable name="R" units="joule_per_kilomole_kelvin" public_interface="in"/>
<variable name="T" units="kelvin" public_interface="in"/>
<variable name="F" units="coulomb_per_mole" public_interface="in"/>
<variable name="time" units="second" public_interface="in"/>
<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_c</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_c</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_Ca</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.5</cn>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Ca_c</ci>
<ci>Ca_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="sodium_current">
<component_ref component="sodium_current_m_gate"/>
<component_ref component="sodium_current_h_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_d_gate"/>
<component_ref component="L_type_Ca_channel_f_gate"/>
</component_ref>
<component_ref component="T_type_Ca_channel">
<component_ref component="T_type_Ca_channel_d_gate"/>
<component_ref component="T_type_Ca_channel_f_gate"/>
</component_ref>
<component_ref component="delayed_rectifying_potassium_current">
<component_ref component="delayed_rectifying_potassium_current_P_a_gate"/>
<component_ref component="delayed_rectifying_potassium_current_P_i_gate"/>
</component_ref>
<component_ref component="linear_background_current"/>
<component_ref component="hyperpolarisation_activated_current">
<component_ref component="hyperpolarisation_activated_current_y_gate"/>
</component_ref>
<component_ref component="sodium_potassium_pump"/>
<component_ref component="calcium_pump_current"/>
<component_ref component="sodium_calcium_pump"/>
<component_ref component="intracellular_concentrations_and_buffer_equations"/>
<component_ref component="cleft_space_equations"/>
<component_ref component="SR_Ca_uptake_and_release"/>
<component_ref component="reversal_potentials"/>
</component_ref>
</group>
<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_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_d_gate"/>
<component_ref component="L_type_Ca_channel_f_gate"/>
</component_ref>
<component_ref component="T_type_Ca_channel">
<component_ref component="T_type_Ca_channel_d_gate"/>
<component_ref component="T_type_Ca_channel_f_gate"/>
</component_ref>
<component_ref component="delayed_rectifying_potassium_current">
<component_ref component="delayed_rectifying_potassium_current_P_a_gate"/>
<component_ref component="delayed_rectifying_potassium_current_P_i_gate"/>
</component_ref>
<component_ref component="hyperpolarisation_activated_current">
<component_ref component="hyperpolarisation_activated_current_y_gate"/>
</component_ref>
</group>
<connection>
<map_components component_1="membrane" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="T_type_Ca_channel" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="delayed_rectifying_potassium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="linear_background_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="hyperpolarisation_activated_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_potassium_pump" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="calcium_pump_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_pump" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="intracellular_concentrations_and_buffer_equations" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="cleft_space_equations" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="SR_Ca_uptake_and_release" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_current"/>
<map_variables variable_1="i_Na" variable_2="i_Na"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="L_type_Ca_channel"/>
<map_variables variable_1="i_Ca_L" variable_2="i_Ca_L"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="T_type_Ca_channel"/>
<map_variables variable_1="i_Ca_T" variable_2="i_Ca_T"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="delayed_rectifying_potassium_current"/>
<map_variables variable_1="i_K" variable_2="i_K"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="linear_background_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_B" variable_2="i_B"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="hyperpolarisation_activated_current"/>
<map_variables variable_1="i_f" variable_2="i_f"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="calcium_pump_current" component_2="membrane"/>
<map_variables variable_1="i_Ca_P" variable_2="i_Ca_P"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_potassium_pump"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_calcium_pump"/>
<map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="intracellular_concentrations_and_buffer_equations"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="cleft_space_equations"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="SR_Ca_uptake_and_release"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="reversal_potentials"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="sodium_current"/>
<map_variables variable_1="E_Na" variable_2="E_Na"/>
</connection>
<connection>
<map_components component_1="intracellular_concentrations_and_buffer_equations" component_2="sodium_current"/>
<map_variables variable_1="i_Na" variable_2="i_Na"/>
</connection>
<connection>
<map_components component_1="cleft_space_equations" component_2="sodium_current"/>
<map_variables variable_1="i_Na" variable_2="i_Na"/>
<map_variables variable_1="Na_c" variable_2="Na_c"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="cleft_space_equations"/>
<map_variables variable_1="i_Ca_L" variable_2="i_Ca_L"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="intracellular_concentrations_and_buffer_equations"/>
<map_variables variable_1="i_Ca_L" variable_2="i_Ca_L"/>
</connection>
<connection>
<map_components component_1="T_type_Ca_channel" component_2="cleft_space_equations"/>
<map_variables variable_1="i_Ca_T" variable_2="i_Ca_T"/>
</connection>
<connection>
<map_components component_1="T_type_Ca_channel" component_2="intracellular_concentrations_and_buffer_equations"/>
<map_variables variable_1="i_Ca_T" variable_2="i_Ca_T"/>
</connection>
<connection>
<map_components component_1="cleft_space_equations" component_2="delayed_rectifying_potassium_current"/>
<map_variables variable_1="i_K" variable_2="i_K"/>
<map_variables variable_1="K_b" variable_2="K_b"/>
</connection>
<connection>
<map_components component_1="intracellular_concentrations_and_buffer_equations" component_2="delayed_rectifying_potassium_current"/>
<map_variables variable_1="i_K" variable_2="i_K"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="delayed_rectifying_potassium_current"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="linear_background_current"/>
<map_variables variable_1="E_Na" variable_2="E_Na"/>
<map_variables variable_1="E_Ca" variable_2="E_Ca"/>
<map_variables variable_1="E_K" variable_2="E_K"/>
</connection>
<connection>
<map_components component_1="cleft_space_equations" component_2="linear_background_current"/>
<map_variables variable_1="i_B_Na" variable_2="i_B_Na"/>
<map_variables variable_1="i_B_Ca" variable_2="i_B_Ca"/>
<map_variables variable_1="i_B_K" variable_2="i_B_K"/>
</connection>
<connection>
<map_components component_1="intracellular_concentrations_and_buffer_equations" component_2="linear_background_current"/>
<map_variables variable_1="i_B_Na" variable_2="i_B_Na"/>
<map_variables variable_1="i_B_Ca" variable_2="i_B_Ca"/>
<map_variables variable_1="i_B_K" variable_2="i_B_K"/>
</connection>
<connection>
<map_components component_1="intracellular_concentrations_and_buffer_equations" component_2="hyperpolarisation_activated_current"/>
<map_variables variable_1="i_f_K" variable_2="i_f_K"/>
<map_variables variable_1="i_f_Na" variable_2="i_f_Na"/>
</connection>
<connection>
<map_components component_1="cleft_space_equations" component_2="hyperpolarisation_activated_current"/>
<map_variables variable_1="i_f_K" variable_2="i_f_K"/>
<map_variables variable_1="i_f_Na" variable_2="i_f_Na"/>
</connection>
<connection>
<map_components component_1="sodium_potassium_pump" component_2="intracellular_concentrations_and_buffer_equations"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
</connection>
<connection>
<map_components component_1="sodium_potassium_pump" component_2="cleft_space_equations"/>
<map_variables variable_1="K_c" variable_2="K_c"/>
<map_variables variable_1="i_NaK" variable_2="i_NaK"/>
</connection>
<connection>
<map_components component_1="calcium_pump_current" component_2="cleft_space_equations"/>
<map_variables variable_1="i_Ca_P" variable_2="i_Ca_P"/>
</connection>
<connection>
<map_components component_1="calcium_pump_current" component_2="intracellular_concentrations_and_buffer_equations"/>
<map_variables variable_1="i_Ca_P" variable_2="i_Ca_P"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_pump" component_2="intracellular_concentrations_and_buffer_equations"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
</connection>
<connection>
<map_components component_1="sodium_calcium_pump" component_2="cleft_space_equations"/>
<map_variables variable_1="Ca_c" variable_2="Ca_c"/>
<map_variables variable_1="Na_c" variable_2="Na_c"/>
<map_variables variable_1="i_NaCa" variable_2="i_NaCa"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="intracellular_concentrations_and_buffer_equations"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection>
<map_components component_1="reversal_potentials" component_2="cleft_space_equations"/>
<map_variables variable_1="K_c" variable_2="K_c"/>
<map_variables variable_1="Na_c" variable_2="Na_c"/>
<map_variables variable_1="Ca_c" variable_2="Ca_c"/>
</connection>
<connection>
<map_components component_1="intracellular_concentrations_and_buffer_equations" component_2="cleft_space_equations"/>
<map_variables variable_1="Vol" variable_2="Vol"/>
</connection>
<connection>
<map_components component_1="intracellular_concentrations_and_buffer_equations" component_2="SR_Ca_uptake_and_release"/>
<map_variables variable_1="i_rel" variable_2="i_rel"/>
<map_variables variable_1="i_up" variable_2="i_up"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="V_i" variable_2="V_i"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="sodium_current_m_gate"/>
<map_variables variable_1="m" variable_2="m"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="sodium_current_h_gate"/>
<map_variables variable_1="h1" variable_2="h1"/>
<map_variables variable_1="h2" variable_2="h2"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_d_gate"/>
<map_variables variable_1="d_L" variable_2="d_L"/>
<map_variables variable_1="d_L_infinity" variable_2="d_L_infinity"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="L_type_Ca_channel" component_2="L_type_Ca_channel_f_gate"/>
<map_variables variable_1="f_L" variable_2="f_L"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="T_type_Ca_channel" component_2="T_type_Ca_channel_d_gate"/>
<map_variables variable_1="d_T" variable_2="d_T"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="T_type_Ca_channel" component_2="T_type_Ca_channel_f_gate"/>
<map_variables variable_1="f_T" variable_2="f_T"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="delayed_rectifying_potassium_current" component_2="delayed_rectifying_potassium_current_P_a_gate"/>
<map_variables variable_1="P_a" variable_2="P_a"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="delayed_rectifying_potassium_current" component_2="delayed_rectifying_potassium_current_P_i_gate"/>
<map_variables variable_1="P_i" variable_2="P_i"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="hyperpolarisation_activated_current" component_2="hyperpolarisation_activated_current_y_gate"/>
<map_variables variable_1="y" variable_2="y"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<rdf:RDF>
<rdf:Seq rdf:about="rdf:#citationAuthorsSeq">
<rdf:li rdf:resource="rdf:#author1Vcard"/>
<rdf:li rdf:resource="rdf:#author2Vcard"/>
<rdf:li rdf:resource="rdf:#author3Vcard"/>
<rdf:li rdf:resource="rdf:#author4Vcard"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#be39a29b-0829-4d71-8ca3-5cb8e71d6220">
<bqs:JournalArticle rdf:resource="rdf:#ccca7b13-1d58-4de0-880b-bb8e5b3a0e84"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1f7e39e4-c589-424e-b1e1-087042a2918e">
<vCard:Given>Peter</vCard:Given>
<vCard:Family>Villiger</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author4Vcard">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#author4VcardN"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a047a17c-984a-42c7-92b8-b1fce8634695">
<vCard:N rdf:resource="rdf:#1f7e39e4-c589-424e-b1e1-087042a2918e"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author1Vcard">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#author1VcardN"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ff75bcdc-f5c1-470d-83a6-6e441fd6628c">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
<vCard:Other>unknown</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="">
<cmeta:modification rdf:resource="rdf:#fff5b863-1b9a-4086-97e4-2706cd17995b"/>
<dcterms:created rdf:resource="rdf:#44b0544d-7a05-4a72-8615-0a3569d6b566"/>
<dc:creator rdf:resource="rdf:#e5a3ad1c-9241-4f04-93fc-23e48f5cdc9d"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author3VcardN">
<vCard:Given>W</vCard:Given>
<vCard:Family>Giles</vCard:Family>
<vCard:Other>R</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author2VcardN">
<vCard:Given>J</vCard:Given>
<vCard:Family>Clark</vCard:Family>
<vCard:Other>W</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author2Vcard">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#author2VcardN"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b859390d-732f-4fc4-bf64-c23b19a7e169">
<vCard:Orgname>Oxford University</vCard:Orgname>
<vCard:Orgunit>unknown</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#81a6687f-782b-4061-879d-088b6ffd132e">
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
<rdf:value>penny.noble@physiol.ox.ac.uk</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ccca7b13-1d58-4de0-880b-bb8e5b3a0e84">
<dc:creator rdf:resource="rdf:#citationAuthorsSeq"/>
<dc:title>A Mathematical Model of a Rabbit Sinoatrial Node Cell</dc:title>
<bqs:volume>266</bqs:volume>
<bqs:first_page>C832</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#44eae6fb-5245-45f2-96ae-2582f6506a4b"/>
<dcterms:issued rdf:resource="rdf:#618443dc-d57f-4b90-931c-dd75f55f6ff5"/>
<bqs:last_page>C852</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#18f15d86-eb6b-44d5-966b-0204df5e20a1">
<dcterms:W3CDTF>2006-03-31</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#618443dc-d57f-4b90-931c-dd75f55f6ff5">
<dcterms:W3CDTF>1994-01-01</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#44b0544d-7a05-4a72-8615-0a3569d6b566">
<dcterms:W3CDTF>2002-01-30</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author3Vcard">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#author3VcardN"/>
</rdf:Description>
<rdf:Description rdf:about="#demir_clark_giles_murphey_1994_version01">
<bqs:reference rdf:resource="rdf:#be39a29b-0829-4d71-8ca3-5cb8e71d6220"/>
<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>sinoatrial node</rdf:li>
<rdf:li>rabbit</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author1VcardN">
<vCard:Given>S</vCard:Given>
<vCard:Family>Demir</vCard:Family>
<vCard:Other>S</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fff5b863-1b9a-4086-97e4-2706cd17995b">
<dcterms:modified rdf:resource="rdf:#18f15d86-eb6b-44d5-966b-0204df5e20a1"/>
<rdf:value>
added metadata
</rdf:value>
<cmeta:modifier rdf:resource="rdf:#a047a17c-984a-42c7-92b8-b1fce8634695"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e5a3ad1c-9241-4f04-93fc-23e48f5cdc9d">
<vCard:ORG rdf:resource="rdf:#b859390d-732f-4fc4-bf64-c23b19a7e169"/>
<vCard:EMAIL rdf:resource="rdf:#81a6687f-782b-4061-879d-088b6ffd132e"/>
<vCard:N rdf:resource="rdf:#ff75bcdc-f5c1-470d-83a6-6e441fd6628c"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author4VcardN">
<vCard:Given>C</vCard:Given>
<vCard:Family>Murphey</vCard:Family>
<vCard:Other>R</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#44eae6fb-5245-45f2-96ae-2582f6506a4b">
<dc:title>American Journal of Physiology</dc:title>
</rdf:Description>
</rdf:RDF>
</model>
