- Author:
- pmr2.import <nobody@models.cellml.org>
- Date:
- 2006-07-09 08:01:51+12:00
- Desc:
- committing version01 of kurata_hisatome_imanishi_shibamoto_2002
- Permanent Source URI:
- https://models.cellml.org/workspace/kurata_hisatome_imanishi_shibamoto_2002/rawfile/08d216a5a63620e6e2fbc5e9fc0048929edef390/kurata_hisatome_imanishi_shibamoto_2002.cellml
<?xml version='1.0' encoding='utf-8'?>
<!-- FILE : kurata_model_2002.xml
CREATED : 7th October 2003
LAST MODIFIED : 20th April 2005
AUTHOR : Catherine Lloyd
Bioengineering Instute
The University of Auckland
MODEL STATUS : This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/01/2002 CellML Metadata 1.0 Specification.
DESCRIPTION : This file contains a CellML description of Kurata et al.'s 2002 mathematical model of the primary pacemaker cell.
CHANGES:
20/04/2005 - PJV - Made MathML id's unique
--><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="kurata_hisatome_imanishi_shibamoto_2002_version01" name="kurata_hisatome_imanishi_shibamoto_2002_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Improved Mathematical Model for the Primary Pacemaker Cell</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 is the original unchecked version of the model imported from the previous
CellML model repository, 24-Jan-2006.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
The pacemaker activity of the sinoatrial (SA) node initiates the spontaneous beating of the heart. Over the past decade a large body of data on the ionic currents underlying this pacemaker activity has been elucidated in patch-clamp experiments. On the basis of this data, several mathematical models describing the pace maker activity of a single rabbit SA node cell have been developed. These include:
<itemizedlist>
<listitem>
<para>
<ulink url="${HTML_EXMPL_D_SAN_MODEL}">Demir <emphasis>et al.</emphasis> Sinoatrial Node Model, 1994</ulink>
</para>
</listitem>;
<listitem>
<para>
<ulink url="${HTML_EXMPL_DOKOS_MODEL}">Modelling the Ion Currents Underlying Sinoatrial Node Pacemaker Activity, Dokos <emphasis>et al.</emphasis>, 1996</ulink>
</para>
</listitem>;
<listitem>
<para>
<ulink url="${HTML_EXMPL_DOKOS_MODEL_II}">A Model of Sinoatrial Node Vagal Control, Dokos <emphasis>et al.</emphasis>, 1996</ulink>
</para>
</listitem>;
<listitem>
<para>
<ulink url="${HTML_EXMPL_D99_SAN_MODEL}">Demir <emphasis>et al.</emphasis> Sinoatrial Node Model, 1999</ulink>
</para>
</listitem>;
<listitem>
<para>
<ulink url="${HTML_EXMPL_Z_SAN_MODEL}">Zhang <emphasis>et al.</emphasis> Sinoatrial Node Model, 2000</ulink>
</para>
</listitem>; and
<listitem>
<para>
<ulink url="${HTML_EXMPL_B_SAN_MODEL}">Boyett <emphasis>et al.</emphasis> Sinoatrial Node Model, 2001</ulink>
</para>
</listitem>.
</itemizedlist>
</para>
<para>
With the elucidation of more experimental data, these models have become increasingly complex, with equations that describe calcium buffering, intracellular compartmentation, and additional ionic currents through channels, pumps and exchangers. In the study by Kurata <emphasis>et al.</emphasis> described here, the aim is to develop an improved mathematical model of a single primary pacemaker cell of the rabbit SA node, in order to investigate the mechanisms underlying pacemaker generation. On the basis of the most recent experimental data, Kurata <emphasis>et al.</emphasis> updated the previously published models (see list above) in several ways:
</para>
<itemizedlist>
<listitem>
<para>A sustained inward current has been added (<emphasis>I<subscript>st</subscript>
</emphasis>).</para>
</listitem>
<listitem>
<para>The voltage- and Ca<superscript>2+</superscript>-dependent activation kinetics of the L-type calcium channel current (<emphasis>I<subscript>Ca,L</subscript>
</emphasis>) have been reformulated.</para>
</listitem>
<listitem>
<para>Expressions for the activation kinetics of the rapidly activating delayed rectifier potassium current (<emphasis>I<subscript>Kr</subscript>
</emphasis>) have been updated.</para>
</listitem>
<listitem>
<para>Revised kinetic equations for two 4-AP sensitive currents (<emphasis>I<subscript>to</subscript>
</emphasis> and <emphasis>I<subscript>sus</subscript>
</emphasis>) have been incorporated.</para>
</listitem>
<listitem>
<para>Voltage- and concentration-dependent kinetics of the Na<superscript>+</superscript>-K<superscript>+</superscript> pump current (<emphasis>I<subscript>NaK</subscript>
</emphasis>) have been reformulated.</para>
</listitem>
<listitem>
<para>The subsarcolemmal space (see <xref linkend="fig_cell_diagram"/> below) as a diffusion barrier for intracellular Ca<superscript>2+</superscript> has been added.</para>
</listitem>
</itemizedlist>
<para>
The model was validated by comparing simulation results with experimental data, and also by comparing them with those of previous mathematical models. The authors concluded that their model represented a significant improvement over previous models because it can:
<itemizedlist>
<listitem>
<para>Simulate whole cell voltage-clamp data for <emphasis>I<subscript>Ca,L</subscript>
</emphasis>, <emphasis>I<subscript>Kr</subscript>
</emphasis>, and <emphasis>I<subscript>st</subscript>
</emphasis>;</para>
</listitem>
<listitem>
<para>Reproduce the wave shapes of spontaneous action potentials and ionic currents during action potential clamp recordings; and</para>
</listitem>
<listitem>
<para>Mimic the effects of channel blockers or Ca <superscript>2+</superscript> buffers on pacemaker activity more accurately than the previous models.</para>
</listitem>
</itemizedlist>
</para>
<para>
The model has been described here in CellML (the raw CellML description of the Kurata <emphasis>et al.</emphasis> 2002 model can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>).
</para>
<para>
The complete original paper reference is cited below:
</para>
<para>
<ulink url="http://ajpheart.physiology.org/cgi/content/abstract/283/5/H2074">Dynamical description of sinoatrial node pacemaking: improved mathematical model for primary pacemaker cell</ulink>, Yasutaka Kurata, Ichiro Hisatome, Sunao Imanishi, and Toshishige Shibamoto, 2002, <ulink url="http://ajpheart.physiology.org/">
<emphasis>American Journal of Physiology</emphasis>
</ulink>, 283, H2074-H2101. (<ulink url="http://ajpheart.physiology.org/cgi/content/full/283/5/H2074">Full text (HTML)</ulink> and <ulink url="http://ajpheart.physiology.org/cgi/reprint/283/5/H2074.pdf">PDF</ulink> versions of the article are available on the <emphasis>American Journal of Physiology</emphasis> website.) <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=12384487&dopt=Abstract">PubMed ID: 12384487</ulink>
</para>
<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>reaction diagram</title>
</objectinfo>
<imagedata fileref="kurata_2002.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram depicting the intracellular compartments for Ca<superscript>2+</superscript>.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<units name="first_order_rate_constant">
<unit units="millisecond" exponent="-1"/>
</units>
<units name="second_order_rate_constant">
<unit units="millimolar" exponent="-1"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="millisecond">
<unit units="second" prefix="milli"/>
</units>
<units name="millivolt">
<unit units="volt" prefix="milli"/>
</units>
<units name="flux">
<unit units="millimolar"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="per_millivolt_millisecond">
<unit units="millivolt" exponent="-1"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="nanoS_per_picoF">
<unit units="siemens" prefix="nano"/>
<unit units="farad" prefix="pico" exponent="-1"/>
</units>
<units name="picoF">
<unit units="farad" prefix="pico"/>
</units>
<units name="per_picoF">
<unit units="farad" prefix="pico" exponent="-1"/>
</units>
<units name="picoA_per_millimolar_picoF">
<unit units="ampere" prefix="pico"/>
<unit units="millimolar" exponent="-1"/>
<unit units="farad" prefix="pico" exponent="-1"/>
</units>
<units name="picoA">
<unit units="ampere" prefix="pico"/>
</units>
<units name="millimolar">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="joule_per_mole_kelvin">
<unit units="joule"/>
<unit units="mole" exponent="-1"/>
<unit units="kelvin" exponent="-1"/>
</units>
<units name="coulomb_per_mole">
<unit units="coulomb" exponent="-1"/>
<unit units="mole" prefix="milli"/>
</units>
<units name="picol">
<unit units="litre" prefix="pico"/>
</units>
<component name="environment">
<variable units="millisecond" public_interface="out" name="time"/>
</component>
<component name="membrane">
<variable units="millivolt" public_interface="out" name="V"/>
<variable units="joule_per_mole_kelvin" public_interface="out" name="R" initial_value="8.3144"/>
<variable units="kelvin" public_interface="out" name="T" initial_value="310.15"/>
<variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96485"/>
<variable units="picoF" name="Cm" initial_value="32.0"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="picoA" public_interface="in" name="i_Na"/>
<variable units="picoA" public_interface="in" name="i_CaT"/>
<variable units="picoA" public_interface="in" name="i_CaL"/>
<variable units="picoA" public_interface="in" name="i_Kr"/>
<variable units="picoA" public_interface="in" name="i_Ks"/>
<variable units="picoA" public_interface="in" name="i_to"/>
<variable units="picoA" public_interface="in" name="i_sus"/>
<variable units="picoA" public_interface="in" name="i_h"/>
<variable units="picoA" public_interface="in" name="i_st"/>
<variable units="picoA" public_interface="in" name="i_b_Na"/>
<variable units="picoA" public_interface="in" name="i_NaCa"/>
<variable units="picoA" public_interface="in" name="i_NaK"/>
<variable units="picoA" public_interface="in" name="i_K_ACh"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="membrane_voltage_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> V </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> i_CaL </ci>
<ci> i_CaT </ci>
<ci> i_Kr </ci>
<ci> i_Ks </ci>
<ci> i_to </ci>
<ci> i_sus </ci>
<ci> i_h </ci>
<ci> i_st </ci>
<ci> i_Na </ci>
<ci> i_b_Na </ci>
<ci> i_K_ACh </ci>
<ci> i_NaK </ci>
<ci> i_NaCa </ci>
</apply>
</apply>
<ci> Cm </ci>
</apply>
</apply>
</math>
</component>
<component name="reversal_potentials">
<variable units="millivolt" public_interface="out" name="E_Na"/>
<variable units="millivolt" public_interface="out" name="E_K"/>
<variable units="joule_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="millimolar" public_interface="in" name="Nai"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="millimolar" public_interface="in" name="Ki"/>
<variable units="millimolar" public_interface="in" name="Ko"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="E_Na_calculation">
<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> Nao </ci>
<ci> Nai </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="E_K_calculation">
<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> Ko </ci>
<ci> Ki </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_calcium_channel_current">
<variable units="picoA" public_interface="out" name="i_CaL"/>
<variable units="nanoS_per_picoF" name="g_CaL" initial_value="0.58"/>
<variable units="millivolt" name="E_CaL" initial_value="45.0"/>
<variable units="millisecond" 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_sub"/>
<variable units="dimensionless" private_interface="in" name="d"/>
<variable units="dimensionless" private_interface="in" name="f"/>
<variable units="dimensionless" private_interface="in" name="fCa"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_CaL_calculation">
<eq/>
<ci> i_CaL </ci>
<apply>
<times/>
<ci> g_CaL </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_CaL </ci>
</apply>
<ci> d </ci>
<ci> f </ci>
<ci> fCa </ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_calcium_channel_current_d_gate">
<variable units="dimensionless" public_interface="out" name="d"/>
<variable units="dimensionless" name="d_infinity"/>
<variable units="millisecond" name="tau_d"/>
<variable units="first_order_rate_constant" name="alpha_d"/>
<variable units="first_order_rate_constant" name="beta_d"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="L_type_calcium_channel_current_d_gate_d_infinity_calculation">
<eq/>
<ci> d_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 14.1 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 6.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="L_type_calcium_channel_current_d_gate_tau_d_calculation">
<eq/>
<ci> tau_d </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<ci> alpha_d </ci>
<ci> beta_d </ci>
</apply>
</apply>
</apply>
<apply id="alpha_d_calculation">
<eq/>
<ci> alpha_d </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_millisecond"> -0.02839 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 35.0 </cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 35.0 </cn>
</apply>
<cn cellml:units="millivolt"> 2.5 </cn>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_millivolt_millisecond"> 0.0849 </cn>
<ci> V </ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 4.8 </cn>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_d_calculation">
<eq/>
<ci> beta_d </ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_millisecond"> 0.01143 </cn>
<apply>
<minus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 5.0 </cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 5.0 </cn>
</apply>
<cn cellml:units="millivolt"> 2.5 </cn>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
</apply>
<apply id="L_type_calcium_channel_current_d_gate_tau_dd_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> d </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> d_infinity </ci>
<ci> d </ci>
</apply>
<ci> tau_d </ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_calcium_channel_current_f_gate">
<variable units="dimensionless" public_interface="out" name="f"/>
<variable units="dimensionless" name="f_infinity"/>
<variable units="millisecond" name="tau_f"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="L_type_calcium_channel_current_f_gate_infinity_calculation">
<eq/>
<ci> f_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 30.0 </cn>
</apply>
<cn cellml:units="millivolt"> 5.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="L_type_calcium_channel_current_f_gate_tau_calculation">
<eq/>
<ci> tau_f </ci>
<apply>
<plus/>
<cn cellml:units="millisecond"> 44.3 </cn>
<apply>
<times/>
<cn cellml:units="millisecond"> 257.1 </cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 32.5 </cn>
</apply>
<cn cellml:units="millivolt"> 13.9 </cn>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="L_type_calcium_channel_current_f_gate_df_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> f </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> f_infinity </ci>
<ci> f </ci>
</apply>
<ci> tau_f </ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_calcium_channel_current_fCa_gate">
<variable units="dimensionless" public_interface="out" name="fCa"/>
<variable units="first_order_rate_constant" name="alpha_fCa" initial_value="0.035"/>
<variable units="second_order_rate_constant" name="beta_fCa" initial_value="60.0"/>
<variable units="dimensionless" name="fCa_infinity"/>
<variable units="millisecond" name="tau_fCa"/>
<variable units="millimolar" name="Km_fCa" initial_value="0.00035"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millimolar" public_interface="in" name="Ca_sub"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="fCa_infinity_calculation">
<eq/>
<ci> fCa_infinity </ci>
<apply>
<divide/>
<ci> Km_fCa </ci>
<apply>
<plus/>
<ci> Km_fCa </ci>
<ci> Ca_sub </ci>
</apply>
</apply>
</apply>
<apply id="tau_fCa_calculation">
<eq/>
<ci> tau_fCa </ci>
<apply>
<divide/>
<ci> fCa_infinity </ci>
<ci> alpha_fCa </ci>
</apply>
</apply>
<apply id="dfCa_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> fCa </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> fCa_infinity </ci>
<ci> fCa </ci>
</apply>
<ci> tau_fCa </ci>
</apply>
</apply>
</math>
</component>
<component name="T_type_calcium_channel_current">
<variable units="picoA" public_interface="out" name="i_CaT"/>
<variable units="nanoS_per_picoF" name="g_CaT" initial_value="0.458"/>
<variable units="millivolt" name="E_CaT" initial_value="45.0"/>
<variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="dimensionless" private_interface="in" name="d"/>
<variable units="dimensionless" private_interface="in" name="f"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_CaT_calculation">
<eq/>
<ci> i_CaT </ci>
<apply>
<times/>
<ci> g_CaT </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_CaT </ci>
</apply>
<ci> d </ci>
<ci> f </ci>
</apply>
</apply>
</math>
</component>
<component name="T_type_calcium_channel_current_d_gate">
<variable units="dimensionless" public_interface="out" name="d"/>
<variable units="dimensionless" name="d_infinity"/>
<variable units="millisecond" name="tau_d"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="T_type_calcium_channel_current_d_gate_d_infinity_calculation">
<eq/>
<ci> d_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 26.3 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 6.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="T_type_calcium_channel_current_d_gate_tau_d_calculation">
<eq/>
<ci> tau_d </ci>
<apply>
<divide/>
<cn cellml:units="millisecond"> 1.0 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="millisecond"> 1.068 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 26.3 </cn>
</apply>
<cn cellml:units="millivolt"> 30.0 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millisecond"> 1.068 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 26.3 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 30.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="T_type_calcium_channel_current_d_gate_dd_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> d </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> d_infinity </ci>
<ci> d </ci>
</apply>
<ci> tau_d </ci>
</apply>
</apply>
</math>
</component>
<component name="T_type_calcium_channel_current_f_gate">
<variable units="dimensionless" public_interface="out" name="f"/>
<variable units="dimensionless" name="f_infinity"/>
<variable units="millisecond" name="tau_f"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="T_type_calcium_channel_current_f_gate_infinity_calculation">
<eq/>
<ci> f_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </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>
<apply id="T_type_calcium_channel_current_f_gate_tau_calculation">
<eq/>
<ci> tau_f </ci>
<apply>
<divide/>
<cn cellml:units="millisecond"> 1.0 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.0153 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 61.7 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 83.3 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.015 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 26.3 </cn>
</apply>
<cn cellml:units="millivolt"> 15.380 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="T_type_calcium_channel_current_f_gate_df_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> f </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> f_infinity </ci>
<ci> f </ci>
</apply>
<ci> tau_f </ci>
</apply>
</apply>
</math>
</component>
<component name="rapidly_activating_delayed_rectifier_potassium_current">
<variable units="picoA" public_interface="out" name="i_Kr"/>
<variable units="nanoS_per_picoF" name="g_Kr"/>
<variable units="millisecond" 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="millimolar" public_interface="in" name="Ko"/>
<variable units="dimensionless" private_interface="in" name="paS"/>
<variable units="dimensionless" private_interface="in" name="paF"/>
<variable units="dimensionless" private_interface="in" name="pi"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_Kr_calculation">
<eq/>
<ci> i_Kr </ci>
<apply>
<times/>
<ci> g_Kr </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.6 </cn>
<ci> paF </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.4 </cn>
<ci> paS </ci>
</apply>
<ci> pi </ci>
</apply>
</apply>
<apply id="g_Kr_calculation">
<eq/>
<ci> g_Kr </ci>
<apply>
<times/>
<cn cellml:units="nanoS_per_picoF"> 0.025 </cn>
<apply>
<power/>
<ci> Ko </ci>
<cn cellml:units="dimensionless"> 0.59 </cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="rapidly_activating_delayed_rectifier_potassium_current_pa_gate">
<variable units="dimensionless" public_interface="out" name="paS"/>
<variable units="dimensionless" public_interface="out" name="paF"/>
<variable units="dimensionless" name="pa_infinity"/>
<variable units="millisecond" name="tau_paS"/>
<variable units="millisecond" name="tau_paF"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="pa_infinity_calculation">
<eq/>
<ci> pa_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 23.2 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 10.6 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_paS_calculation">
<eq/>
<ci> tau_paS </ci>
<apply>
<divide/>
<cn cellml:units="millisecond"> 0.84655 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.0042 </cn>
<apply>
<exp/>
<apply>
<divide/>
<ci> V </ci>
<cn cellml:units="millivolt"> 17.0 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.00015 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 21.6 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_paF_calculation">
<eq/>
<ci> tau_paF </ci>
<apply>
<divide/>
<cn cellml:units="millisecond"> 0.84655 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.0372 </cn>
<apply>
<exp/>
<apply>
<divide/>
<ci> V </ci>
<cn cellml:units="millivolt"> 15.9 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.00096 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 22.5 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="dpaS_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> paS </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> pa_infinity </ci>
<ci> paS </ci>
</apply>
<ci> tau_paS </ci>
</apply>
</apply>
<apply id="dpaF_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> paF </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> pa_infinity </ci>
<ci> paF </ci>
</apply>
<ci> tau_paF </ci>
</apply>
</apply>
</math>
</component>
<component name="rapidly_activating_delayed_rectifier_potassium_current_pi_gate">
<variable units="dimensionless" public_interface="out" name="pi"/>
<variable units="dimensionless" name="pi_infinity"/>
<variable units="millisecond" name="tau_pi"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="pi_infinity_calculation">
<eq/>
<ci> pi_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 28.6 </cn>
</apply>
<cn cellml:units="millivolt"> 17.1 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_pi_calculation">
<eq/>
<ci> tau_pi </ci>
<apply>
<divide/>
<cn cellml:units="millisecond"> 1.0 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.1 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 54.645 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.656 </cn>
<apply>
<exp/>
<apply>
<divide/>
<ci> V </ci>
<cn cellml:units="millivolt"> 106.157 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="dpi_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> pi </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> pi_infinity </ci>
<ci> pi </ci>
</apply>
<ci> tau_pi </ci>
</apply>
</apply>
</math>
</component>
<component name="slowly_activating_delayed_rectifier_potassium_current">
<variable units="picoA" public_interface="out" name="i_Ks"/>
<variable units="nanoS_per_picoF" name="g_Ks" initial_value="0.0259"/>
<variable units="millivolt" name="E_Ks"/>
<variable units="millisecond" 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" name="Ki"/>
<variable units="millimolar" public_interface="in" name="Ko"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="joule_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="dimensionless" private_interface="in" name="n"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_Ks_calculation">
<eq/>
<ci> i_Ks </ci>
<apply>
<times/>
<ci> g_Ks </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_Ks </ci>
</apply>
<apply>
<power/>
<ci> n </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<apply id="E_Ks_calculation">
<eq/>
<ci> E_Ks </ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
<ci> F </ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<apply>
<plus/>
<ci> Ko </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.12 </cn>
<ci> Nao </ci>
</apply>
</apply>
<apply>
<plus/>
<ci> Ki </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.12 </cn>
<ci> Nai </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slowly_activating_delayed_rectifier_potassium_current_n_gate">
<variable units="dimensionless" public_interface="out" name="n"/>
<variable units="dimensionless" name="n_infinity"/>
<variable units="millisecond" name="tau_n"/>
<variable units="first_order_rate_constant" name="alpha_n"/>
<variable units="first_order_rate_constant" name="beta_n"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="n_infinity_calculation">
<eq/>
<ci> n_infinity </ci>
<apply>
<divide/>
<ci> alpha_n </ci>
<apply>
<plus/>
<ci> alpha_n </ci>
<ci> beta_n </ci>
</apply>
</apply>
</apply>
<apply id="tau_n_calculation">
<eq/>
<ci> tau_n </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<ci> alpha_n </ci>
<ci> beta_n </ci>
</apply>
</apply>
</apply>
<apply id="alpha_n_calculation">
<eq/>
<ci> alpha_n </ci>
<apply>
<divide/>
<cn cellml:units="first_order_rate_constant"> 0.014 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 40.0 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 9.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_n_calculation">
<eq/>
<ci> beta_n </ci>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant"> 0.001 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 45.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="dn_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> n </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> n_infinity </ci>
<ci> n </ci>
</apply>
<ci> tau_n </ci>
</apply>
</apply>
</math>
</component>
<component name="AP_sensitive_currents">
<variable units="picoA" public_interface="out" name="i_to"/>
<variable units="picoA" public_interface="out" name="i_sus"/>
<variable units="nanoS_per_picoF" name="g_to" initial_value="0.18"/>
<variable units="nanoS_per_picoF" name="g_sus" initial_value="0.02"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<variable units="millisecond" 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="q"/>
<variable units="dimensionless" private_interface="in" name="r"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_to_calculation">
<eq/>
<ci> i_to </ci>
<apply>
<times/>
<ci> g_to </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K </ci>
</apply>
<ci> q </ci>
<ci> r </ci>
</apply>
</apply>
<apply id="i_sus_calculation">
<eq/>
<ci> i_sus </ci>
<apply>
<times/>
<ci> g_sus </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K </ci>
</apply>
<ci> r </ci>
</apply>
</apply>
</math>
</component>
<component name="AP_sensitive_currents_q_gate">
<variable units="dimensionless" public_interface="out" name="q"/>
<variable units="dimensionless" name="q_infinity"/>
<variable units="millisecond" name="tau_q"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="q_infinity_calculation">
<eq/>
<ci> q_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 49.0 </cn>
</apply>
<cn cellml:units="millivolt"> 13.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_q_calculation">
<eq/>
<ci> tau_q </ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="millisecond"> 65.17 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.57 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> -0.08 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 44.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.065 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.1 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 45.93 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="millisecond"> 10.1 </cn>
</apply>
</apply>
<apply id="dq_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> q </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> q_infinity </ci>
<ci> q </ci>
</apply>
<ci> tau_q </ci>
</apply>
</apply>
</math>
</component>
<component name="AP_sensitive_currents_r_gate">
<variable units="dimensionless" public_interface="out" name="r"/>
<variable units="dimensionless" name="r_infinity"/>
<variable units="millisecond" name="tau_r"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="r_infinity_calculation">
<eq/>
<ci> r_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 19.3 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 15.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_r_calculation">
<eq/>
<ci> tau_r </ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="millisecond"> 21.826 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 1.037 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.09 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 30.61 </cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.369 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> -0.12 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 23.84 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="millisecond"> 4.172 </cn>
</apply>
</apply>
<apply id="dr_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> r </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> r_infinity </ci>
<ci> r </ci>
</apply>
<ci> tau_r </ci>
</apply>
</apply>
</math>
</component>
<component name="hyperpolarisation_activated_current">
<variable units="picoA" public_interface="out" name="i_h"/>
<variable units="picoA" public_interface="out" name="i_h_Na"/>
<variable units="picoA" public_interface="out" name="i_h_K"/>
<variable units="nanoS_per_picoF" name="g_h_Na" initial_value="0.1368381"/>
<variable units="nanoS_per_picoF" name="g_h_K" initial_value="0.2201619"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<variable units="millivolt" public_interface="in" name="E_Na"/>
<variable units="millisecond" 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="y"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_h_calculation">
<eq/>
<ci> i_h </ci>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci> g_h_Na </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_Na </ci>
</apply>
</apply>
<apply>
<times/>
<ci> g_h_K </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K </ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci> y </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<apply id="i_h_Na_calculation">
<eq/>
<ci> i_h_Na </ci>
<apply>
<minus/>
<ci> i_h </ci>
<ci> i_h_K </ci>
</apply>
</apply>
<apply id="i_h_K_calculation">
<eq/>
<ci> i_h_K </ci>
<apply>
<minus/>
<ci> i_h </ci>
<ci> i_h_Na </ci>
</apply>
</apply>
</math>
</component>
<component name="hyperpolarisation_activated_current_y_gate">
<variable units="dimensionless" public_interface="out" name="y"/>
<variable units="dimensionless" name="y_infinity"/>
<variable units="millisecond" name="tau_y"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="y_infinity_calculation">
<eq/>
<ci> y_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 64.0 </cn>
</apply>
<cn cellml:units="millivolt"> 13.5 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_y_calculation">
<eq/>
<ci> tau_y </ci>
<apply>
<divide/>
<cn cellml:units="millisecond"> 0.71665 </cn>
<apply>
<plus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 386.9 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 45.3 </cn>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 73.08 </cn>
</apply>
<cn cellml:units="millivolt"> 19.23 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="dy_dt">
<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="sustained_inward_current">
<variable units="picoA" public_interface="out" name="i_st"/>
<variable units="nanoS_per_picoF" name="g_st" initial_value="0.015"/>
<variable units="millivolt" name="E_st" initial_value="37.4"/>
<variable units="millisecond" 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="qa"/>
<variable units="dimensionless" private_interface="in" name="qi"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_st_calculation">
<eq/>
<ci> i_st </ci>
<apply>
<times/>
<ci> g_st </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_st </ci>
</apply>
<ci> qa </ci>
<ci> qi </ci>
</apply>
</apply>
</math>
</component>
<component name="sustained_inward_current_qa_gate">
<variable units="dimensionless" public_interface="out" name="qa"/>
<variable units="dimensionless" name="qa_infinity"/>
<variable units="millisecond" name="tau_qa"/>
<variable units="first_order_rate_constant" name="alpha_qa"/>
<variable units="first_order_rate_constant" name="beta_qa"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="qa_infinity_calculation">
<eq/>
<ci> qa_infinity </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 57.0 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 5.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_qa_calculation">
<eq/>
<ci> tau_qa </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<ci> alpha_qa </ci>
<ci> beta_qa </ci>
</apply>
</apply>
</apply>
<apply id="alpha_qa_calculation">
<eq/>
<ci> alpha_qa </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.15 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 11.0 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.2 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 700.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_qa_calculation">
<eq/>
<ci> beta_qa </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="millisecond"> 16.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<ci> V </ci>
<cn cellml:units="millivolt"> 8.0 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millisecond"> 15.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<ci> V </ci>
<cn cellml:units="millivolt"> 50.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="dqa_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> qa </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> qa_infinity </ci>
<ci> qa </ci>
</apply>
<ci> tau_qa </ci>
</apply>
</apply>
</math>
</component>
<component name="sustained_inward_current_qi_gate">
<variable units="dimensionless" public_interface="out" name="qi"/>
<variable units="dimensionless" name="qi_infinity"/>
<variable units="millisecond" name="tau_qi"/>
<variable units="first_order_rate_constant" name="alpha_qi"/>
<variable units="first_order_rate_constant" name="beta_qi"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="qi_infinity_calculation">
<eq/>
<ci> qi_infinity </ci>
<apply>
<divide/>
<ci> alpha_qi </ci>
<apply>
<plus/>
<ci> alpha_qi </ci>
<ci> beta_qi </ci>
</apply>
</apply>
</apply>
<apply id="tau_qi_calculation">
<eq/>
<ci> tau_qi </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<ci> alpha_qi </ci>
<ci> beta_qi </ci>
</apply>
</apply>
</apply>
<apply id="alpha_qi_calculation">
<eq/>
<ci> alpha_qi </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 0.1504 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="millisecond"> 3100.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<ci> V </ci>
<cn cellml:units="millivolt"> 13.0 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millisecond"> 700.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<ci> V </ci>
<cn cellml:units="millivolt"> 70.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_qi_calculation">
<eq/>
<ci> beta_qi </ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 0.1504 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="millisecond"> 95.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 10.0 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="millisecond"> 50.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 700.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 0.000229 </cn>
<apply>
<plus/>
<cn cellml:units="millisecond"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 5.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="dqi_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> qi </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> qi_infinity </ci>
<ci> qi </ci>
</apply>
<ci> tau_qi </ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_channel_current">
<variable units="picoA" public_interface="out" name="i_Na"/>
<variable units="nanoS_per_picoF" name="g_Na" initial_value="0.00"/>
<variable units="millivolt" name="E_mh"/>
<variable units="dimensionless" name="F_Na"/>
<variable units="millisecond" 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" name="Ki"/>
<variable units="millimolar" public_interface="in" name="Ko"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="joule_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="dimensionless" private_interface="in" name="m"/>
<variable units="dimensionless" private_interface="in" name="hs"/>
<variable units="dimensionless" private_interface="in" name="hf"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_Na_calculation">
<eq/>
<ci> i_Na </ci>
<apply>
<times/>
<ci> g_Na </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_mh </ci>
</apply>
<apply>
<power/>
<ci> m </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> F_Na </ci>
</apply>
<ci> hf </ci>
</apply>
<apply>
<times/>
<ci> F_Na </ci>
<ci> hs </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="E_mh_calculation">
<eq/>
<ci> E_mh </ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
<ci> F </ci>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<apply>
<plus/>
<ci> Nao </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.12 </cn>
<ci> Ko </ci>
</apply>
</apply>
<apply>
<plus/>
<ci> Nai </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.12 </cn>
<ci> Ki </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="F_Na_calculation">
<eq/>
<ci> F_Na </ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.0952 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> -0.063 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 14.4 </cn>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 1.66 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> -0.225 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 43.7 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> 0.0869 </cn>
</apply>
</apply>
</math>
</component>
<component name="sodium_channel_current_m_gate">
<variable units="dimensionless" public_interface="out" name="m"/>
<variable units="dimensionless" name="m_infinity"/>
<variable units="millisecond" name="tau_m"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="m_infinitm_calculation">
<eq/>
<ci> m_infinity </ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 21.4 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 6.0 </cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless"> -1.0 </cn>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
<apply id="tau_m_calculation">
<eq/>
<ci> tau_m </ci>
<apply>
<divide/>
<cn cellml:units="millisecond"> 0.6247 </cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.832 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 36.7 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 2.985 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.627 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 45.01 </cn>
</apply>
<cn cellml:units="millivolt"> 12.195 </cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> 0.04 </cn>
</apply>
</apply>
</apply>
<apply id="dm_dt">
<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_channel_current_h_gate">
<variable units="dimensionless" public_interface="out" name="hs"/>
<variable units="dimensionless" public_interface="out" name="hf"/>
<variable units="dimensionless" name="h_infinity"/>
<variable units="millisecond" name="tau_hs"/>
<variable units="millisecond" name="tau_hf"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="h_infinith_calculation">
<eq/>
<ci> h_infinity </ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 70.6 </cn>
</apply>
<cn cellml:units="millivolt"> 6.4 </cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
</apply>
<apply id="tau_hf_calculation">
<eq/>
<ci> tau_hf </ci>
<apply>
<divide/>
<cn cellml:units="millisecond"> 0.4757 </cn>
<apply>
<plus/>
<apply>
<power/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 15600 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 20.0 </cn>
</apply>
<cn cellml:units="millivolt"> 16.949 </cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 0.92 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.037 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 20.0 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 11.494 </cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_hs_calculation">
<eq/>
<ci> tau_hs </ci>
<apply>
<divide/>
<cn cellml:units="millisecond"> 0.4757 </cn>
<apply>
<plus/>
<apply>
<power/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 168000 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 20.0 </cn>
</apply>
<cn cellml:units="millivolt"> 18.182 </cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 2.58 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.94 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 20.0 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 15.873 </cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="dhf_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> hf </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> h_infinity </ci>
<ci> hf </ci>
</apply>
<ci> tau_hf </ci>
</apply>
</apply>
<apply id="dhs_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> hs </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> h_infinity </ci>
<ci> hs </ci>
</apply>
<ci> tau_hs </ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_dependent_background_current">
<variable units="picoA" public_interface="out" name="i_b_Na"/>
<variable units="nanoS_per_picoF" name="g_b_Na" initial_value="0.0054"/>
<variable units="millivolt" public_interface="in" name="E_Na"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_b_Na_calculation">
<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>
</math>
</component>
<component name="background_muscarinic_potassium_channel_current">
<variable units="picoA" public_interface="out" name="i_K_ACh"/>
<variable units="picoA_per_millimolar_picoF" name="g_K_ACh"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millimolar" public_interface="in" name="Ki"/>
<variable units="millimolar" public_interface="in" name="Ko"/>
<variable units="joule_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"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_K_ACh_calculation">
<eq/>
<ci> i_K_ACh </ci>
<apply>
<times/>
<ci> g_K_ACh </ci>
<apply>
<minus/>
<ci> Ki </ci>
<apply>
<times/>
<ci> Ko </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>
</apply>
<apply id="g_K_ACh_calculation">
<eq/>
<ci> g_K_ACh </ci>
<apply>
<times/>
<cn cellml:units="picoA_per_millimolar_picoF"> 0.0011 </cn>
<apply>
<power/>
<ci> Ko </ci>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_potassium_pump_current">
<variable units="picoA" public_interface="out" name="i_NaK"/>
<variable units="millimolar" name="Km_Kp" initial_value="1.4"/>
<variable units="millimolar" name="Km_Nap" initial_value="14.0"/>
<variable units="picoA" name="i_NaK_max" initial_value="3.6"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="millimolar" public_interface="in" name="Ko"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millivolt" public_interface="in" name="E_Na"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_NaK_calculation">
<eq/>
<ci> i_NaK </ci>
<apply>
<times/>
<ci> i_NaK_max </ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<power/>
<apply>
<divide/>
<ci> Km_Kp </ci>
<ci> Ko </ci>
</apply>
<cn cellml:units="dimensionless"> 1.2 </cn>
</apply>
</apply>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<power/>
<apply>
<divide/>
<ci> Km_Nap </ci>
<ci> Nai </ci>
</apply>
<cn cellml:units="dimensionless"> 1.3 </cn>
</apply>
</apply>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci> V </ci>
<ci> E_Na </ci>
</apply>
<cn cellml:units="millivolt"> 120.0 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 30.0 </cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_calcium_exchange_current">
<variable units="picoA" public_interface="out" name="i_NaCa"/>
<variable units="per_picoF" name="kNaCa" initial_value="125.0"/>
<variable units="dimensionless" name="x1"/>
<variable units="dimensionless" name="x2"/>
<variable units="dimensionless" name="x3"/>
<variable units="dimensionless" name="x4"/>
<variable units="dimensionless" name="k41"/>
<variable units="dimensionless" name="k34"/>
<variable units="dimensionless" name="k23"/>
<variable units="dimensionless" name="k21"/>
<variable units="dimensionless" name="k32"/>
<variable units="dimensionless" name="k43"/>
<variable units="dimensionless" name="k12"/>
<variable units="dimensionless" name="k14"/>
<variable units="dimensionless" name="Qci" initial_value="0.1369"/>
<variable units="dimensionless" name="Qn" initial_value="0.4315"/>
<variable units="dimensionless" name="Qco" initial_value="0.0"/>
<variable units="millimolar" name="K3ni" initial_value="26.44"/>
<variable units="millimolar" name="Kci" initial_value="0.0207"/>
<variable units="millimolar" name="K1ni" initial_value="395.3"/>
<variable units="millimolar" name="K2ni" initial_value="2.289"/>
<variable units="millimolar" name="Kcni" initial_value="26.44"/>
<variable units="millimolar" name="K3no" initial_value="4.663"/>
<variable units="millimolar" name="K1no" initial_value="1628.0"/>
<variable units="millimolar" name="K2no" initial_value="561.4"/>
<variable units="millimolar" name="Kco" initial_value="3.663"/>
<variable units="dimensionless" name="do"/>
<variable units="dimensionless" name="di"/>
<variable units="millimolar" public_interface="in" name="Cao"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="millimolar" public_interface="in" name="Ca_sub"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="joule_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 id="i_NaCa_calculation">
<eq/>
<ci> i_NaCa </ci>
<apply>
<times/>
<ci> kNaCa </ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci> x2 </ci>
<ci> k21 </ci>
</apply>
<apply>
<times/>
<ci> x1 </ci>
<ci> k12 </ci>
</apply>
</apply>
<apply>
<plus/>
<ci> x1 </ci>
<ci> x2 </ci>
<ci> x3 </ci>
<ci> x4 </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="x1_calculation">
<eq/>
<ci> x1 </ci>
<apply>
<plus/>
<apply>
<times/>
<ci> k41 </ci>
<ci> k34 </ci>
<apply>
<plus/>
<ci> k23 </ci>
<ci> k21 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k21 </ci>
<ci> k32 </ci>
<apply>
<plus/>
<ci> k43 </ci>
<ci> k41 </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="x2_calculation">
<eq/>
<ci> x2 </ci>
<apply>
<plus/>
<apply>
<times/>
<ci> k32 </ci>
<ci> k43 </ci>
<apply>
<plus/>
<ci> k14 </ci>
<ci> k12 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k41 </ci>
<ci> k12 </ci>
<apply>
<plus/>
<ci> k34 </ci>
<ci> k32 </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="x3_calculation">
<eq/>
<ci> x3 </ci>
<apply>
<plus/>
<apply>
<times/>
<ci> k14 </ci>
<ci> k43 </ci>
<apply>
<plus/>
<ci> k23 </ci>
<ci> k21 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k12 </ci>
<ci> k23 </ci>
<apply>
<plus/>
<ci> k43 </ci>
<ci> k41 </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="x4_calculation">
<eq/>
<ci> x4 </ci>
<apply>
<plus/>
<apply>
<times/>
<ci> k23 </ci>
<ci> k34 </ci>
<apply>
<plus/>
<ci> k14 </ci>
<ci> k12 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k14 </ci>
<ci> k21 </ci>
<apply>
<plus/>
<ci> k34 </ci>
<ci> k32 </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="k43_calculation">
<eq/>
<ci> k43 </ci>
<apply>
<divide/>
<ci> Nai </ci>
<apply>
<plus/>
<ci> K3ni </ci>
<ci> Nai </ci>
</apply>
</apply>
</apply>
<apply id="k12_calculation">
<eq/>
<ci> k12 </ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci> Ca_sub </ci>
<ci> Kci </ci>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci> Qci </ci>
</apply>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<ci> di </ci>
</apply>
</apply>
<apply id="k14_calculation">
<eq/>
<ci> k14 </ci>
<apply>
<times/>
<apply>
<divide/>
<ci> Nai </ci>
<ci> K1ni </ci>
</apply>
<apply>
<divide/>
<ci> Nai </ci>
<ci> K2ni </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> Nai </ci>
<ci> K3ni </ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci> Qn </ci>
<ci> V </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
<ci> di </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="k41_calculation">
<eq/>
<ci> k41 </ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci> Qn </ci>
</apply>
<ci> V </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="di_calculation">
<eq/>
<ci> di </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<apply>
<divide/>
<ci> Ca_sub </ci>
<ci> Kci </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci> Qci </ci>
</apply>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<ci> Nai </ci>
<ci> Kcni </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> Nai </ci>
<ci> K1ni </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<apply>
<divide/>
<ci> Nai </ci>
<ci> K2ni </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> Nai </ci>
<ci> K3ni </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="k34_calculation">
<eq/>
<ci> k34 </ci>
<apply>
<divide/>
<ci> Nao </ci>
<apply>
<plus/>
<ci> K3no </ci>
<ci> Nao </ci>
</apply>
</apply>
</apply>
<apply id="k21_calculation">
<eq/>
<ci> k21 </ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci> Cao </ci>
<ci> Kco </ci>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> Qco </ci>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<ci> do </ci>
</apply>
</apply>
<apply id="k23_calculation">
<eq/>
<ci> k23 </ci>
<apply>
<times/>
<apply>
<divide/>
<ci> Nao </ci>
<ci> K1no </ci>
</apply>
<apply>
<divide/>
<ci> Nao </ci>
<ci> K2no </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> Nao </ci>
<ci> K3no </ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> Qn </ci>
<ci> V </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
<ci> do </ci>
</apply>
</apply>
<apply id="k32_calculation">
<eq/>
<ci> k32 </ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> Qn </ci>
<ci> V </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="do_calculation">
<eq/>
<ci> do </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<apply>
<divide/>
<ci> Cao </ci>
<ci> Kco </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci> Qco </ci>
</apply>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> Nao </ci>
<ci> K1no </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> Nao </ci>
<ci> K2no </ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> Nao </ci>
<ci> K3no </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="intracellular_calcium_dynamics">
<variable units="flux" public_interface="out" name="j_Ca_dif"/>
<variable units="flux" public_interface="out" name="j_rel"/>
<variable units="flux" public_interface="out" name="j_up"/>
<variable units="flux" public_interface="out" name="j_tr"/>
<variable units="millisecond" name="tau_dif_Ca" initial_value="0.04"/>
<variable units="millisecond" name="tau_tr" initial_value="60.0"/>
<variable units="millimolar" name="K_rel" initial_value="0.0012"/>
<variable units="flux" name="P_up" initial_value="0.005"/>
<variable units="first_order_rate_constant" name="P_rel" initial_value="5.0"/>
<variable units="millimolar" name="K_up" initial_value="0.0006"/>
<variable units="millimolar" public_interface="in" name="Ca_up"/>
<variable units="millimolar" public_interface="in" name="Cai"/>
<variable units="millimolar" public_interface="in" name="Ca_sub"/>
<variable units="millimolar" public_interface="in" name="Ca_rel"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="j_Ca_dif_calculation">
<eq/>
<ci> j_Ca_dif </ci>
<apply>
<divide/>
<apply>
<minus/>
<ci> Ca_sub </ci>
<ci> Cai </ci>
</apply>
<ci> tau_dif_Ca </ci>
</apply>
</apply>
<apply id="j_rel_calculation">
<eq/>
<ci> j_rel </ci>
<apply>
<times/>
<ci> P_rel </ci>
<apply>
<divide/>
<apply>
<minus/>
<ci> Ca_rel </ci>
<ci> Ca_sub </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<power/>
<apply>
<divide/>
<ci> K_rel </ci>
<ci> Ca_sub </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="j_up_calculation">
<eq/>
<ci> j_up </ci>
<apply>
<times/>
<ci> P_up </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> K_up </ci>
<ci> Cai </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="j_tr_calculation">
<eq/>
<ci> j_tr </ci>
<apply>
<divide/>
<apply>
<minus/>
<ci> Ca_up </ci>
<ci> Ca_rel </ci>
</apply>
<ci> tau_tr </ci>
</apply>
</apply>
</math>
</component>
<component name="intracellular_ion_concentrations">
<variable units="millimolar" public_interface="out" name="Cai"/>
<variable units="millimolar" public_interface="out" name="Cao" initial_value="2.0"/>
<variable units="millimolar" public_interface="out" name="Ca_up"/>
<variable units="millimolar" public_interface="out" name="Ca_rel"/>
<variable units="millimolar" public_interface="out" name="Ca_sub"/>
<variable units="millimolar" public_interface="out" name="Nai"/>
<variable units="millimolar" public_interface="out" name="Nao" initial_value="140.0"/>
<variable units="millimolar" public_interface="out" name="Ki"/>
<variable units="millimolar" public_interface="out" name="Ko" initial_value="5.4"/>
<variable units="millimolar" public_interface="out" name="Mgi" initial_value="2.5"/>
<variable units="picol" name="V_i" initial_value="1.5939"/>
<variable units="picol" name="V_rel" initial_value="0.0042"/>
<variable units="picol" name="V_up" initial_value="0.0406"/>
<variable units="picol" name="V_sub" initial_value="0.035"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="picoA" public_interface="in" name="i_b_Na"/>
<variable units="picoA" public_interface="in" name="i_CaL"/>
<variable units="picoA" public_interface="in" name="i_CaT"/>
<variable units="picoA" public_interface="in" name="i_NaCa"/>
<variable units="flux" public_interface="in" name="j_rel"/>
<variable units="flux" public_interface="in" name="j_Ca_dif"/>
<variable units="flux" public_interface="in" name="j_tr"/>
<variable units="flux" public_interface="in" name="j_up"/>
<variable units="picoA" public_interface="in" name="i_h_Na"/>
<variable units="picoA" public_interface="in" name="i_st"/>
<variable units="picoA" public_interface="in" name="i_NaK"/>
<variable units="picoA" public_interface="in" name="i_Na"/>
<variable units="picoA" public_interface="in" name="i_Kr"/>
<variable units="picoA" public_interface="in" name="i_Ks"/>
<variable units="picoA" public_interface="in" name="i_to"/>
<variable units="picoA" public_interface="in" name="i_sus"/>
<variable units="picoA" public_interface="in" name="i_h_K"/>
<variable units="picoA" public_interface="in" name="i_K_ACh"/>
<variable units="millimolar" public_interface="in" name="TMC_tot"/>
<variable units="millimolar" public_interface="in" name="CM_tot"/>
<variable units="millimolar" public_interface="in" name="TC_tot"/>
<variable units="millimolar" public_interface="in" name="CQ_tot"/>
<variable units="flux" public_interface="in" name="delta_fTMC"/>
<variable units="flux" public_interface="in" name="delta_fCMi"/>
<variable units="flux" public_interface="in" name="delta_fCMs"/>
<variable units="flux" public_interface="in" name="delta_fTC"/>
<variable units="flux" public_interface="in" name="delta_fCQ"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="dNai_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Nai </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> i_h_Na </ci>
<ci> i_st </ci>
<ci> i_b_Na </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> i_NaK </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> i_NaCa </ci>
</apply>
<ci> i_Na </ci>
</apply>
</apply>
<apply>
<times/>
<ci> F </ci>
<ci> V_i </ci>
</apply>
</apply>
</apply>
<apply id="dKi_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Ki </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> i_Kr </ci>
<ci> i_Ks </ci>
<ci> i_to </ci>
<ci> i_sus </ci>
<ci> i_h_K </ci>
<ci> i_K_ACh </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> -2.0 </cn>
<ci> i_NaK </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> F </ci>
<ci> V_i </ci>
</apply>
</apply>
</apply>
<apply id="dCai_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Cai </ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci> j_Ca_dif </ci>
<ci> V_sub </ci>
</apply>
<apply>
<times/>
<ci> j_up </ci>
<ci> V_up </ci>
</apply>
</apply>
<ci> V_i </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> CM_tot </ci>
<ci> delta_fCMi </ci>
</apply>
<apply>
<times/>
<ci> TC_tot </ci>
<ci> delta_fTC </ci>
</apply>
<apply>
<times/>
<ci> TMC_tot </ci>
<ci> delta_fTMC </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="dCa_sub_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Ca_sub </ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci> i_CaL </ci>
<ci> i_CaT </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> i_NaCa </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> F </ci>
</apply>
</apply>
<apply>
<times/>
<ci> j_rel </ci>
<ci> V_rel </ci>
</apply>
</apply>
<ci> V_sub </ci>
</apply>
<apply>
<plus/>
<ci> j_Ca_dif </ci>
<apply>
<times/>
<ci> CM_tot </ci>
<ci> delta_fCMs </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="Ca_up_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Ca_up </ci>
</apply>
<apply>
<minus/>
<ci> j_up </ci>
<apply>
<times/>
<ci> j_tr </ci>
<apply>
<divide/>
<ci> V_rel </ci>
<ci> V_up </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="Ca_rel_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Ca_rel </ci>
</apply>
<apply>
<minus/>
<ci> j_tr </ci>
<apply>
<plus/>
<ci> j_rel </ci>
<apply>
<times/>
<ci> CQ_tot </ci>
<ci> delta_fCQ </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_buffering">
<variable units="millimolar" public_interface="out" name="TC_tot" initial_value="0.031"/>
<variable units="millimolar" public_interface="out" name="TMC_tot" initial_value="0.062"/>
<variable units="millimolar" public_interface="out" name="CM_tot" initial_value="0.045"/>
<variable units="millimolar" public_interface="out" name="CQ_tot" initial_value="10.0"/>
<variable units="flux" public_interface="out" name="delta_fTC"/>
<variable units="flux" public_interface="out" name="delta_fTMC"/>
<variable units="flux" public_interface="out" name="delta_fCMs"/>
<variable units="flux" public_interface="out" name="delta_fCMi"/>
<variable units="flux" public_interface="out" name="delta_fCQ"/>
<variable units="flux" name="delta_fTMM"/>
<variable units="millimolar" name="fTMM"/>
<variable units="millimolar" name="fCMi"/>
<variable units="millimolar" name="fCMs"/>
<variable units="millimolar" name="fTC"/>
<variable units="millimolar" name="fTMC"/>
<variable units="millimolar" name="fCQ"/>
<variable units="second_order_rate_constant" name="kf_TC" initial_value="88.8"/>
<variable units="second_order_rate_constant" name="kf_TMM" initial_value="2.277"/>
<variable units="second_order_rate_constant" name="kf_TMC" initial_value="227.7"/>
<variable units="second_order_rate_constant" name="kf_CM" initial_value="227.7"/>
<variable units="second_order_rate_constant" name="kf_CQ" initial_value="0.534"/>
<variable units="first_order_rate_constant" name="kb_TC" initial_value="0.446"/>
<variable units="first_order_rate_constant" name="kb_TMC" initial_value="0.00751"/>
<variable units="first_order_rate_constant" name="kb_TMM" initial_value="0.751"/>
<variable units="first_order_rate_constant" name="kb_CM" initial_value="0.542"/>
<variable units="first_order_rate_constant" name="kb_CQ" initial_value="0.445"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millimolar" public_interface="in" name="Mgi"/>
<variable units="millimolar" public_interface="in" name="Cai"/>
<variable units="millimolar" public_interface="in" name="Ca_sub"/>
<variable units="millimolar" public_interface="in" name="Ca_rel"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="dfTC_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> fTC </ci>
</apply>
<ci> delta_fTC </ci>
</apply>
<apply id="delta_fTC_calculation">
<eq/>
<ci> delta_fTC </ci>
<apply>
<minus/>
<apply>
<times/>
<ci> kf_TC </ci>
<ci> Cai </ci>
<apply>
<minus/>
<cn cellml:units="millimolar"> 1.0 </cn>
<ci> fTC </ci>
</apply>
</apply>
<apply>
<times/>
<ci> kb_TC </ci>
<ci> fTC </ci>
</apply>
</apply>
</apply>
<apply id="dfTMC_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> fTMC </ci>
</apply>
<ci> delta_fTMC </ci>
</apply>
<apply id="delta_fTMC_calculation">
<eq/>
<ci> delta_fTMC </ci>
<apply>
<minus/>
<apply>
<times/>
<ci> kf_TMC </ci>
<ci> Cai </ci>
<apply>
<minus/>
<cn cellml:units="millimolar"> 1.0 </cn>
<apply>
<plus/>
<ci> fTMC </ci>
<ci> fTMM </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kb_TMC </ci>
<ci> fTMC </ci>
</apply>
</apply>
</apply>
<apply id="dfTMM_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> fTMM </ci>
</apply>
<ci> delta_fTMM </ci>
</apply>
<apply id="delta_fTMM_calculation">
<eq/>
<ci> delta_fTMM </ci>
<apply>
<minus/>
<apply>
<times/>
<ci> kf_TMM </ci>
<ci> Mgi </ci>
<apply>
<minus/>
<cn cellml:units="millimolar"> 1.0 </cn>
<apply>
<plus/>
<ci> fTMC </ci>
<ci> fTMM </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kb_TMM </ci>
<ci> fTMM </ci>
</apply>
</apply>
</apply>
<apply id="dfCMi_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> fCMi </ci>
</apply>
<ci> delta_fCMi </ci>
</apply>
<apply id="delta_fCMi_calculation">
<eq/>
<ci> delta_fCMi </ci>
<apply>
<minus/>
<apply>
<times/>
<ci> kf_CM </ci>
<ci> Cai </ci>
<apply>
<minus/>
<cn cellml:units="millimolar"> 1.0 </cn>
<ci> fCMi </ci>
</apply>
</apply>
<apply>
<times/>
<ci> kb_CM </ci>
<ci> fCMi </ci>
</apply>
</apply>
</apply>
<apply id="dfCMs_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> fCMs </ci>
</apply>
<ci> delta_fCMs </ci>
</apply>
<apply id="delta_fCMs_calculation">
<eq/>
<ci> delta_fCMs </ci>
<apply>
<minus/>
<apply>
<times/>
<ci> kf_CM </ci>
<ci> Ca_sub </ci>
<apply>
<minus/>
<cn cellml:units="millimolar"> 1.0 </cn>
<ci> fCMs </ci>
</apply>
</apply>
<apply>
<times/>
<ci> kb_CM </ci>
<ci> fCMs </ci>
</apply>
</apply>
</apply>
<apply id="dfCQ_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> fCQ </ci>
</apply>
<ci> delta_fCQ </ci>
</apply>
<apply id="delta_fCQ_calculation">
<eq/>
<ci> delta_fCQ </ci>
<apply>
<minus/>
<apply>
<times/>
<ci> kf_CQ </ci>
<ci> Ca_rel </ci>
<apply>
<minus/>
<cn cellml:units="millimolar"> 1.0 </cn>
<ci> fCQ </ci>
</apply>
</apply>
<apply>
<times/>
<ci> kb_CQ </ci>
<ci> fCQ </ci>
</apply>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="L_type_calcium_channel_current">
<component_ref component="L_type_calcium_channel_current_d_gate"/>
<component_ref component="L_type_calcium_channel_current_f_gate"/>
<component_ref component="L_type_calcium_channel_current_fCa_gate"/>
</component_ref>
<component_ref component="T_type_calcium_channel_current">
<component_ref component="T_type_calcium_channel_current_d_gate"/>
<component_ref component="T_type_calcium_channel_current_f_gate"/>
</component_ref>
<component_ref component="rapidly_activating_delayed_rectifier_potassium_current">
<component_ref component="rapidly_activating_delayed_rectifier_potassium_current_pa_gate"/>
<component_ref component="rapidly_activating_delayed_rectifier_potassium_current_pi_gate"/>
</component_ref>
<component_ref component="slowly_activating_delayed_rectifier_potassium_current">
<component_ref component="slowly_activating_delayed_rectifier_potassium_current_n_gate"/>
</component_ref>
<component_ref component="AP_sensitive_currents">
<component_ref component="AP_sensitive_currents_q_gate"/>
<component_ref component="AP_sensitive_currents_r_gate"/>
</component_ref>
<component_ref component="hyperpolarisation_activated_current">
<component_ref component="hyperpolarisation_activated_current_y_gate"/>
</component_ref>
<component_ref component="sustained_inward_current">
<component_ref component="sustained_inward_current_qa_gate"/>
<component_ref component="sustained_inward_current_qi_gate"/>
</component_ref>
<component_ref component="sodium_channel_current">
<component_ref component="sodium_channel_current_m_gate"/>
<component_ref component="sodium_channel_current_h_gate"/>
</component_ref>
<component_ref component="sodium_dependent_background_current"/>
<component_ref component="background_muscarinic_potassium_channel_current"/>
<component_ref component="sodium_potassium_pump_current"/>
<component_ref component="sodium_calcium_exchange_current"/>
<component_ref component="intracellular_calcium_dynamics"/>
<component_ref component="intracellular_ion_concentrations"/>
<component_ref component="calcium_buffering"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="L_type_calcium_channel_current">
<component_ref component="L_type_calcium_channel_current_d_gate"/>
<component_ref component="L_type_calcium_channel_current_f_gate"/>
<component_ref component="L_type_calcium_channel_current_fCa_gate"/>
</component_ref>
<component_ref component="T_type_calcium_channel_current">
<component_ref component="T_type_calcium_channel_current_d_gate"/>
<component_ref component="T_type_calcium_channel_current_f_gate"/>
</component_ref>
<component_ref component="rapidly_activating_delayed_rectifier_potassium_current">
<component_ref component="rapidly_activating_delayed_rectifier_potassium_current_pa_gate"/>
<component_ref component="rapidly_activating_delayed_rectifier_potassium_current_pi_gate"/>
</component_ref>
<component_ref component="slowly_activating_delayed_rectifier_potassium_current">
<component_ref component="slowly_activating_delayed_rectifier_potassium_current_n_gate"/>
</component_ref>
<component_ref component="AP_sensitive_currents">
<component_ref component="AP_sensitive_currents_q_gate"/>
<component_ref component="AP_sensitive_currents_r_gate"/>
</component_ref>
<component_ref component="hyperpolarisation_activated_current">
<component_ref component="hyperpolarisation_activated_current_y_gate"/>
</component_ref>
<component_ref component="sustained_inward_current">
<component_ref component="sustained_inward_current_qa_gate"/>
<component_ref component="sustained_inward_current_qi_gate"/>
</component_ref>
<component_ref component="sodium_channel_current">
<component_ref component="sodium_channel_current_m_gate"/>
<component_ref component="sodium_channel_current_h_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="L_type_calcium_channel_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="T_type_calcium_channel_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="slowly_activating_delayed_rectifier_potassium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="AP_sensitive_currents"/>
<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="sustained_inward_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sodium_channel_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sodium_dependent_background_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="background_muscarinic_potassium_channel_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sodium_calcium_exchange_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sodium_potassium_pump_current"/>
<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="environment" component_1="calcium_buffering"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="reversal_potentials"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="F" variable_1="F"/>
<map_variables variable_2="T" variable_1="T"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="L_type_calcium_channel_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_CaL" variable_1="i_CaL"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="T_type_calcium_channel_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_CaT" variable_1="i_CaT"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Kr" variable_1="i_Kr"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="slowly_activating_delayed_rectifier_potassium_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Ks" variable_1="i_Ks"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="F" variable_1="F"/>
<map_variables variable_2="T" variable_1="T"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="AP_sensitive_currents"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_to" variable_1="i_to"/>
<map_variables variable_2="i_sus" variable_1="i_sus"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="hyperpolarisation_activated_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_h" variable_1="i_h"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="sustained_inward_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_st" variable_1="i_st"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="sodium_channel_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Na" variable_1="i_Na"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="F" variable_1="F"/>
<map_variables variable_2="T" variable_1="T"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="sodium_dependent_background_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="background_muscarinic_potassium_channel_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_K_ACh" variable_1="i_K_ACh"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="F" variable_1="F"/>
<map_variables variable_2="T" variable_1="T"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="sodium_potassium_pump_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_NaK" variable_1="i_NaK"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="sodium_calcium_exchange_current"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="F" variable_1="F"/>
<map_variables variable_2="T" variable_1="T"/>
</connection>
<connection>
<map_components component_2="membrane" component_1="intracellular_ion_concentrations"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="reversal_potentials"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="Ki" variable_1="Ki"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
</connection>
<connection>
<map_components component_2="rapidly_activating_delayed_rectifier_potassium_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="AP_sensitive_currents" component_1="reversal_potentials"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="reversal_potentials" component_1="hyperpolarisation_activated_current"/>
<map_variables variable_2="E_Na" variable_1="E_Na"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="sodium_dependent_background_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_Na" variable_1="E_Na"/>
</connection>
<connection>
<map_components component_2="sodium_potassium_pump_current" component_1="reversal_potentials"/>
<map_variables variable_2="E_Na" variable_1="E_Na"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="L_type_calcium_channel_current"/>
<map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
<map_variables variable_2="i_CaL" variable_1="i_CaL"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="T_type_calcium_channel_current"/>
<map_variables variable_2="i_CaT" variable_1="i_CaT"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
<map_variables variable_2="i_Kr" variable_1="i_Kr"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="slowly_activating_delayed_rectifier_potassium_current"/>
<map_variables variable_2="i_Ks" variable_1="i_Ks"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="Ki" variable_1="Ki"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="AP_sensitive_currents"/>
<map_variables variable_2="i_to" variable_1="i_to"/>
<map_variables variable_2="i_sus" variable_1="i_sus"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="hyperpolarisation_activated_current"/>
<map_variables variable_2="i_h_Na" variable_1="i_h_Na"/>
<map_variables variable_2="i_h_K" variable_1="i_h_K"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="sustained_inward_current"/>
<map_variables variable_2="i_st" variable_1="i_st"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="sodium_channel_current"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="Ki" variable_1="Ki"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
<map_variables variable_2="i_Na" variable_1="i_Na"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="sodium_dependent_background_current"/>
<map_variables variable_2="i_b_Na" variable_1="i_b_Na"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="background_muscarinic_potassium_channel_current"/>
<map_variables variable_2="Ki" variable_1="Ki"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
<map_variables variable_2="i_K_ACh" variable_1="i_K_ACh"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="sodium_potassium_pump_current"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
<map_variables variable_2="i_NaK" variable_1="i_NaK"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="sodium_calcium_exchange_current"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="Cao" variable_1="Cao"/>
<map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
<map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="intracellular_calcium_dynamics"/>
<map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
<map_variables variable_2="Ca_up" variable_1="Ca_up"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
<map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
<map_variables variable_2="j_Ca_dif" variable_1="j_Ca_dif"/>
<map_variables variable_2="j_up" variable_1="j_up"/>
<map_variables variable_2="j_rel" variable_1="j_rel"/>
<map_variables variable_2="j_tr" variable_1="j_tr"/>
</connection>
<connection>
<map_components component_2="intracellular_ion_concentrations" component_1="calcium_buffering"/>
<map_variables variable_2="Ca_rel" variable_1="Ca_rel"/>
<map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
<map_variables variable_2="Mgi" variable_1="Mgi"/>
<map_variables variable_2="TC_tot" variable_1="TC_tot"/>
<map_variables variable_2="TMC_tot" variable_1="TMC_tot"/>
<map_variables variable_2="CM_tot" variable_1="CM_tot"/>
<map_variables variable_2="CQ_tot" variable_1="CQ_tot"/>
<map_variables variable_2="delta_fCMi" variable_1="delta_fCMi"/>
<map_variables variable_2="delta_fCMs" variable_1="delta_fCMs"/>
<map_variables variable_2="delta_fTMC" variable_1="delta_fTMC"/>
<map_variables variable_2="delta_fTC" variable_1="delta_fTC"/>
<map_variables variable_2="delta_fCQ" variable_1="delta_fCQ"/>
</connection>
<connection>
<map_components component_2="L_type_calcium_channel_current_d_gate" component_1="L_type_calcium_channel_current"/>
<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_calcium_channel_current_f_gate" component_1="L_type_calcium_channel_current"/>
<map_variables variable_2="f" variable_1="f"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="L_type_calcium_channel_current_fCa_gate" component_1="L_type_calcium_channel_current"/>
<map_variables variable_2="fCa" variable_1="fCa"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="Ca_sub" variable_1="Ca_sub"/>
</connection>
<connection>
<map_components component_2="T_type_calcium_channel_current_d_gate" component_1="T_type_calcium_channel_current"/>
<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="T_type_calcium_channel_current_f_gate" component_1="T_type_calcium_channel_current"/>
<map_variables variable_2="f" variable_1="f"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="sodium_channel_current_m_gate" component_1="sodium_channel_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_channel_current_h_gate" component_1="sodium_channel_current"/>
<map_variables variable_2="hs" variable_1="hs"/>
<map_variables variable_2="hf" variable_1="hf"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="rapidly_activating_delayed_rectifier_potassium_current_pa_gate" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
<map_variables variable_2="paF" variable_1="paF"/>
<map_variables variable_2="paS" variable_1="paS"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="rapidly_activating_delayed_rectifier_potassium_current_pi_gate" component_1="rapidly_activating_delayed_rectifier_potassium_current"/>
<map_variables variable_2="pi" variable_1="pi"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="slowly_activating_delayed_rectifier_potassium_current_n_gate" component_1="slowly_activating_delayed_rectifier_potassium_current"/>
<map_variables variable_2="n" variable_1="n"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="AP_sensitive_currents_q_gate" component_1="AP_sensitive_currents"/>
<map_variables variable_2="q" variable_1="q"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="AP_sensitive_currents_r_gate" component_1="AP_sensitive_currents"/>
<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="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>
<connection>
<map_components component_2="sustained_inward_current_qa_gate" component_1="sustained_inward_current"/>
<map_variables variable_2="qa" variable_1="qa"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="sustained_inward_current_qi_gate" component_1="sustained_inward_current"/>
<map_variables variable_2="qi" variable_1="qi"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<rdf:RDF>
<rdf:Bag rdf:about="rdf:#ede04a5f-ec70-4ac7-8127-4bc9c4bfe447">
<rdf:li>Sinoatrial Node</rdf:li>
<rdf:li>cardiac</rdf:li>
<rdf:li>electrophysiology</rdf:li>
</rdf:Bag>
<rdf:Seq rdf:about="rdf:#420607cd-30df-4b97-b70c-c6b596d78eaa">
<rdf:li rdf:resource="rdf:#76cbf67f-6f27-4e74-868d-8e0e893d2cef"/>
<rdf:li rdf:resource="rdf:#a94b1b64-0c2d-4d98-8c82-1b17e7d55b73"/>
<rdf:li rdf:resource="rdf:#861c060a-1e3b-4392-af38-f5522b50a239"/>
<rdf:li rdf:resource="rdf:#6582fbfd-4194-4d3c-8c9e-0c8be7b39034"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#6582fbfd-4194-4d3c-8c9e-0c8be7b39034">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#2219fdd6-6055-42f9-a78c-332dd455fc1c"/>
</rdf:Description>
<rdf:Description rdf:about="#kurata_hisatome_imanishi_shibamoto_2002_version01">
<dc:title>
Kurata et al.'s 2002 mathematical model of the primary pacemaker cell.
</dc:title>
<cmeta:bio_entity>Sinoatrial Node Cell</cmeta:bio_entity>
<cmeta:comment rdf:resource="rdf:#4f9e49d0-204a-4940-a819-f6e99086ee2d"/>
<bqs:reference rdf:resource="rdf:#2dbeaae8-f794-4c1e-9a66-2d44551f4ef2"/>
<bqs:reference rdf:resource="rdf:#fab7a55c-e7e2-410c-80fc-7595d2e07e68"/>
<cmeta:species>Rabbit</cmeta:species>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>
The University of Auckland, Bioengineering Institute
</dc:publisher>
<dcterms:created rdf:resource="rdf:#29923b8a-32eb-4fec-9d8e-6c4acc475ee4"/>
<dc:creator rdf:resource="rdf:#f0aa26ee-dd73-4831-b2af-15b5cdcd42d9"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2219fdd6-6055-42f9-a78c-332dd455fc1c">
<vCard:Given>toshishige</vCard:Given>
<vCard:Family>Shibamoto</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#26f94198-fbca-49b5-b3c3-de0f8d7a5d67">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a94b1b64-0c2d-4d98-8c82-1b17e7d55b73">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#b50ddf4e-2bf2-440f-84c8-e858e55883c3"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#76cbf67f-6f27-4e74-868d-8e0e893d2cef">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#6eb02013-bb77-49ea-a31d-96f946a74d09"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6eb02013-bb77-49ea-a31d-96f946a74d09">
<vCard:Given>Yasutaka</vCard:Given>
<vCard:Family>Kurata</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b50ddf4e-2bf2-440f-84c8-e858e55883c3">
<vCard:Given>Ichiro</vCard:Given>
<vCard:Family>Hisatome</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e43d81fb-4e2f-42a9-8562-41bad3507813">
<vCard:Given>Sunao</vCard:Given>
<vCard:Family>Imanishi</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b809ed12-7ecc-4193-9790-42f82bef66f1">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#84f830bb-62eb-4b72-a93a-63cde4a2fad9">
<dc:creator rdf:resource="rdf:#420607cd-30df-4b97-b70c-c6b596d78eaa"/>
<dc:title>
Dynamical description of sinoatrial node pacemaking: improved
mathematical model for primary pacemaker cell
</dc:title>
<bqs:volume>283</bqs:volume>
<bqs:first_page>H2074</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#e6ce91f0-d15b-4786-bff3-c102316c4619"/>
<dcterms:issued rdf:resource="rdf:#5f0feec3-80f9-49a7-8c88-0bf812e1e834"/>
<bqs:last_page>H2101</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e6ce91f0-d15b-4786-bff3-c102316c4619">
<dc:title>American Journal of Physiology</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f0aa26ee-dd73-4831-b2af-15b5cdcd42d9">
<vCard:ORG rdf:resource="rdf:#5349cdd3-c828-4560-991f-c21f94a4db30"/>
<vCard:EMAIL rdf:resource="rdf:#a71aae75-9c13-44f1-9e99-e47ba286b3eb"/>
<vCard:N rdf:resource="rdf:#b809ed12-7ecc-4193-9790-42f82bef66f1"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4f9e49d0-204a-4940-a819-f6e99086ee2d">
<dc:creator rdf:resource="rdf:#26f94198-fbca-49b5-b3c3-de0f8d7a5d67"/>
<rdf:value>
This is the CellML description of Kurata et al.'s 2002 mathematical
model of the primary pacemaker cell.
</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#21d2ce40-45c8-41ee-85a3-3097a7859a46">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value rdf:resource="rdf:#ede04a5f-ec70-4ac7-8127-4bc9c4bfe447"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fab7a55c-e7e2-410c-80fc-7595d2e07e68">
<dc:subject rdf:resource="rdf:#21d2ce40-45c8-41ee-85a3-3097a7859a46"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5f0feec3-80f9-49a7-8c88-0bf812e1e834">
<dcterms:W3CDTF>2002-11</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#29923b8a-32eb-4fec-9d8e-6c4acc475ee4">
<dcterms:W3CDTF>2003-10-07</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#861c060a-1e3b-4392-af38-f5522b50a239">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#e43d81fb-4e2f-42a9-8562-41bad3507813"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2dbeaae8-f794-4c1e-9a66-2d44551f4ef2">
<bqs:Pubmed_id>12384487</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#84f830bb-62eb-4b72-a93a-63cde4a2fad9"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a71aae75-9c13-44f1-9e99-e47ba286b3eb">
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
<rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5349cdd3-c828-4560-991f-c21f94a4db30">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
</rdf:RDF>
</model>
