- Author:
- pmr2.import <nobody@models.cellml.org>
- Date:
- 2009-06-17 15:43:04+12:00
- Desc:
- committing version02 of rice_jafri_winslow_2000
- Permanent Source URI:
- https://models.cellml.org/workspace/rice_jafri_winslow_2000/rawfile/4533913da7a3d2c499b34891dc5df5d240f99ab1/rice_jafri_winslow_2000.cellml
<?xml version='1.0' encoding='utf-8'?>
<!-- FILE : rice_model2_2000_raw.xml
CREATED : 24th June 2002
LAST MODIFIED : 9th April 2003
AUTHOR : Catherine Lloyd
The Bioengineering Institute
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 Rice et al's 2000
detailed model of ventricular cardiac myocytes. It simulates
action potentials, Ca2+-handling mechanisms and isometric force
generation by the myofilaments.
CHANGES:
18/07/2002 - CML - Added more metadata.
09/04/2003 - AAC - Added publication date information.
--><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:pathway_editor="http://www.physiome.com/pathway_editor/1.0#" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" pathway_editor:rendering_config_file="rice_model2_2000_CellMLrender.xml" cmeta:id="rice_jafri_winslow_2000_version02" name="rice_jafri_winslow_2000_version02">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Modelling Interval-Force Relations in Cardiac Muscle</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>
Short-term interval-force (I-F) relations describe the dependence of cardiac muscle contraction strength for short interbeat intervals. Two phenomena that characterise these short-term I-F relations are restitution and postextrasystolic potentiation. In their 2000 paper, J. Jeremy Rice, M. Saleet Jafri and Raimond L. Winslow use two modelling approaches to study short-term interval-force relations in cardiac muscle. Their first approach employs a relatively simple (or minimal), discrete-time model of excitation-contraction coupling. This model explains restitution and postextrasystolic potentiation in terms of the total amount of Ca<superscript>2+</superscript> stored in the sarcoplasmic reticulum (SR) and a feature in which Ca<superscript>2+</superscript> stored in the SR slowly becomes available for release into the sarcoplasm (see <xref linkend="fig_simple_diagram"/> below).
</para>
<informalfigure float="0" id="fig_simple_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>simple_model</title>
</objectinfo>
<imagedata fileref="rice_2000a.png"/>
</imageobject>
<caption>
Schematic diagram of a simple, discrete-time model of Ca<superscript>2+</superscript>-handling in restitution and postextrasystolic potentiation. In<subscript>SC</subscript> and Out<subscript>SC</subscript> represent Ca<superscript>2+</superscript> influx and efflux across the sarcolemma, and In<subscript>SR</subscript> and Out<subscript>SR</subscript> represent sarcoplasmic reticulum Ca<superscript>2+</superscript> release and uptake.
</caption>
</mediaobject>
</informalfigure>
<para>
In their second modelling approach, J. Jeremy Rice, M. Saleet Jafri and Raimond L. Winslow develop a more detailed ventricular cell model that is able to simulate action potentials, Ca<superscript>2+</superscript>-handling mechanisms, and isometric force generation by the myofilaments. This single comprehensive (or maximal) model combines two previously developed models (see <ulink url="${HTML_EXMPL_JRW_MODEL}"> The Jafri-Rice-Winslow Ventricular Model, 1998</ulink> and <ulink url="${HTML_EXMPL_RICE_MODEL}"> Cooperative Mechanisms in Cardiac Muscle, Rice <emphasis>et al.</emphasis>, 1999</ulink>). The first ventricular cell model describes membrane currents and Ca<superscript>2+</superscript>-handling, and the second model of myofilaments provides the force generation of muscle contraction (see <xref linkend="fig_detailed_diagram"/> below). The ventricular cell model provides the Ca<superscript>2+</superscript> needed to drive the myofilament contraction. This is a feedforward pathway in which Ca<superscript>2+</superscript> binds to troponin and initiates force generation by triggering myofilament contraction. A feedback pathway also links the two models, as the affinity of troponin for Ca<superscript>2+</superscript> is a function of developed force.
</para>
<para>
<ulink url="http://ajpheart.physiology.org/cgi/content/abstract/278/3/H913">Modeling short-term interval-force relations in cardiac muscle</ulink> J. Jeremy Rice, M. Saleet Jafri and Raimond L. Winslow, 2000, <ulink url="http://ajpheart.physiology.org/">
<emphasis>American Journal of Physiology</emphasis>
</ulink>, 278, H913-H931. (<ulink url="http://ajpheart.physiology.org/cgi/content/full/278/3/H913">Full text</ulink> and <ulink url="http://ajpheart.physiology.org/cgi/reprint/278/3/H913.pdf">PDF</ulink> versions of the article are available for Journal Members on the American Journal of Physiology website.) <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=10710361&dopt=Abstract">PubMed ID: 10710361</ulink>
</para>
<informalfigure float="0" id="fig_detailed_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>detailed_model</title>
</objectinfo>
<imagedata fileref="rice_2000b.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram of a detailed, single ventricular cell model which is able to simulate action potentials, Ca<superscript>2+</superscript>-handling and isometric force generation by myofilaments.</caption>
</informalfigure>
<para>
Simulation results generated from the model support experimental data. They suggest that short-term interval-force relations result mainly from the interplay between ryanodine receptor adaptation and SR Ca<superscript>2+</superscript> loading with additional contributions from membrane currents and myofilament activation.
</para>
</sect1>
</article>
</documentation>
<!--
Below, we define some additional units for association with variables and
constants within the model. The identifiers are fairly self-explanatory.
-->
<units name="millisecond">
<unit units="second" prefix="milli"/>
</units>
<units name="per_millisecond">
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<units name="per_second">
<unit units="second" exponent="-1"/>
</units>
<units name="micro_litre">
<unit units="litre" prefix="micro"/>
</units>
<units name="micrometre">
<unit units="metre" prefix="micro"/>
</units>
<units name="cm2">
<unit units="metre" prefix="centi" exponent="2"/>
</units>
<units name="cm_per_second">
<unit units="metre" prefix="centi"/>
<unit units="second" exponent="-1"/>
</units>
<units name="cm_per_millisecond">
<unit units="metre" prefix="centi"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="millivolt">
<unit units="volt" prefix="milli"/>
</units>
<units name="per_millivolt">
<unit units="millivolt" exponent="-1"/>
</units>
<units name="per_millivolt_millisecond">
<unit units="millivolt" exponent="-1"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="millimolar">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="micromolar">
<unit units="mole" prefix="micro"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="micromolar_per_second">
<unit units="micromolar"/>
<unit units="second" exponent="-1"/>
</units>
<units name="micromolar_per_millisecond">
<unit units="micromolar"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="per_micromolar_per_second">
<unit units="micromolar" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="per_micromolar3_per_millisecond">
<unit units="micromolar" exponent="-3"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="per_micromolar4_per_millisecond">
<unit units="micromolar" exponent="-4"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units name="microF_per_cm2">
<unit units="farad" prefix="micro"/>
<unit units="metre" prefix="centi" exponent="-2"/>
</units>
<units name="microA_per_microF">
<unit units="ampere" prefix="micro"/>
<unit units="farad" prefix="micro" exponent="-1"/>
</units>
<units name="milliS_per_microF">
<unit units="siemens" prefix="milli"/>
<unit units="farad" prefix="micro" 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_millimole">
<unit units="coulomb"/>
<unit units="mole" prefix="milli" exponent="-1"/>
</units>
<units name="N_per_mm2">
<unit units="newton"/>
<unit units="metre" prefix="milli" exponent="-2"/>
</units>
<!--
The "environment" component is used to declare variables that are used by
all or most of the other components, in this case just "time".
-->
<component name="environment">
<variable units="millisecond" public_interface="out" name="time"/>
</component>
<component name="cell_geometry">
<variable units="cm2" public_interface="out" name="A_cap" initial_value="1.534E-4"/>
<variable units="micro_litre" public_interface="out" name="V_myo" initial_value="25.84E-6"/>
<variable units="micro_litre" public_interface="out" name="V_JSR" initial_value="0.16E-6"/>
<variable units="micro_litre" public_interface="out" name="V_NSR" initial_value="1.40E-6"/>
<variable units="micro_litre" public_interface="out" name="V_SS" initial_value="0.495E-6"/>
</component>
<!--
The "membrane" component is really the `root' node of our model.
It defines the action potential variable "V" among other things.
-->
<component name="membrane">
<variable units="millivolt" public_interface="out" name="V" initial_value="-84.1638"/>
<variable units="joule_per_mole_kelvin" public_interface="out" name="R" initial_value="8.314"/>
<variable units="kelvin" public_interface="out" name="T" initial_value="310.0"/>
<variable units="coulomb_per_millimole" public_interface="out" name="F" initial_value="96.5"/>
<variable units="microF_per_cm2" name="Cm" initial_value="1.0"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="microA_per_microF" public_interface="in" name="i_Na"/>
<variable units="microA_per_microF" public_interface="in" name="i_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="i_Ca_K"/>
<variable units="microA_per_microF" public_interface="in" name="i_K"/>
<variable units="microA_per_microF" public_interface="in" name="i_K1"/>
<variable units="microA_per_microF" public_interface="in" name="i_NaCa"/>
<variable units="microA_per_microF" public_interface="in" name="i_Kp"/>
<variable units="microA_per_microF" public_interface="in" name="i_p_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="i_Na_b"/>
<variable units="microA_per_microF" public_interface="in" name="i_Ca_b"/>
<variable units="microA_per_microF" public_interface="in" name="i_NaK"/>
<variable units="microA_per_microF" public_interface="in" name="i_ns_Ca"/>
<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>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless"> -1.0 </cn>
<ci> Cm </ci>
</apply>
<apply>
<plus/>
<ci> i_Na </ci>
<ci> i_Ca </ci>
<ci> i_Ca_K </ci>
<ci> i_K </ci>
<ci> i_K1 </ci>
<ci> i_Kp </ci>
<ci> i_NaCa </ci>
<ci> i_NaK </ci>
<ci> i_ns_Ca </ci>
<ci> i_p_Ca </ci>
<ci> i_Ca_b </ci>
<ci> i_Na_b </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current">
<variable units="microA_per_microF" public_interface="out" name="i_Na"/>
<variable units="millivolt" public_interface="out" name="E_Na"/>
<variable units="milliS_per_microF" name="g_Na" initial_value="12.8"/>
<variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_millimole" public_interface="in" name="F"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="dimensionless" private_interface="in" name="m"/>
<variable units="dimensionless" private_interface="in" name="h"/>
<variable units="dimensionless" private_interface="in" name="j"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_Na_calculation">
<eq/>
<ci> i_Na </ci>
<apply>
<times/>
<ci> g_Na </ci>
<apply>
<power/>
<ci> m </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> h </ci>
<ci> j </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_Na </ci>
</apply>
</apply>
</apply>
<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>
</math>
</component>
<component name="fast_sodium_current_m_gate">
<variable units="dimensionless" public_interface="out" name="m" initial_value="0.0326174"/>
<variable units="per_millisecond" name="alpha_m"/>
<variable units="per_millisecond" name="beta_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="alpha_m_calculation">
<eq/>
<ci> alpha_m </ci>
<apply>
<times/>
<cn cellml:units="per_millivolt_millisecond"> 0.32 </cn>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 47.13 </cn>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> -0.1 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 47.13 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_m_calculation">
<eq/>
<ci> beta_m </ci>
<apply>
<times/>
<cn cellml:units="per_millisecond"> 0.08 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 11.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="dm_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> m </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> alpha_m </ci>
<apply>
<minus/>
<cn cellml:units="per_millivolt"> 1.0 </cn>
<ci> m </ci>
</apply>
</apply>
<apply>
<times/>
<ci> beta_m </ci>
<ci> m </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current_h_gate">
<variable units="dimensionless" public_interface="out" name="h" initial_value="0.962369"/>
<variable units="per_millisecond" name="alpha_h"/>
<variable units="per_millisecond" name="beta_h"/>
<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="alpha_h_calculation">
<eq/>
<ci> alpha_h </ci>
<piecewise>
<piece>
<apply>
<times/>
<cn cellml:units="per_millisecond"> 0.135 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<cn cellml:units="millivolt"> 80.0 </cn>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> -6.8 </cn>
</apply>
</apply>
</apply>
<apply>
<lt/>
<ci> V </ci>
<cn cellml:units="millivolt"> -40.0 </cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="per_millisecond"> 0.0 </cn>
</otherwise>
</piecewise>
</apply>
<apply id="beta_h_calculation">
<eq/>
<ci> beta_h </ci>
<piecewise>
<piece>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="per_millisecond"> 3.56 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="millivolt"> 0.079 </cn>
<ci> V </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_millisecond"> 310000.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> 0.35 </cn>
<ci> V </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<lt/>
<ci> V </ci>
<cn cellml:units="millivolt"> -40.0 </cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<cn cellml:units="millisecond"> 0.13 </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"> 10.66 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 11.1 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply id="dh_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> h </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> alpha_h </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> h </ci>
</apply>
</apply>
<apply>
<times/>
<ci> beta_h </ci>
<ci> h </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current_j_gate">
<variable units="dimensionless" public_interface="out" name="j" initial_value="0.975658"/>
<variable units="per_millisecond" name="alpha_j"/>
<variable units="per_millisecond" name="beta_j"/>
<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="alpha_j_calculation">
<eq/>
<ci> alpha_j </ci>
<piecewise>
<piece>
<apply>
<times/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_millivolt_millisecond"> -127140.0
</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> 0.2444 </cn>
<ci> V </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_millivolt_millisecond"> 0.00003474 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> -0.04391 </cn>
<ci> V </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 37.78 </cn>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> 0.311 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 79.23 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<lt/>
<ci> V </ci>
<cn cellml:units="millivolt"> -40.0 </cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="per_millisecond"> 0.0 </cn>
</otherwise>
</piecewise>
</apply>
<apply id="beta_j_calculation">
<eq/>
<ci> beta_j </ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millisecond"> 0.1212 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> -0.01052 </cn>
<ci> V </ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> -0.1378 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 40.14 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<lt/>
<ci> V </ci>
<cn cellml:units="millivolt"> -40.0 </cn>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millisecond"> 0.3 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> -0.0000002535 </cn>
<ci> V </ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> -0.1 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 32.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply id="dj_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> j </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> alpha_j </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> j </ci>
</apply>
</apply>
<apply>
<times/>
<ci> beta_j </ci>
<ci> j </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_dependent_potassium_current">
<variable units="microA_per_microF" public_interface="out" name="i_K"/>
<variable units="milliS_per_microF" name="g_K"/>
<variable units="millivolt" name="E_K"/>
<variable units="cm_per_second" name="PNa_K" initial_value="0.01833"/>
<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="Ko"/>
<variable units="millimolar" public_interface="in" name="Ki"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_millimole" public_interface="in" name="F"/>
<variable units="dimensionless" private_interface="in" name="X"/>
<variable units="dimensionless" private_interface="in" name="Xi"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="g_K_calculation">
<eq/>
<ci> g_K </ci>
<apply>
<times/>
<cn cellml:units="milliS_per_microF"> 0.1128 </cn>
<apply>
<root/>
<apply>
<divide/>
<ci> Ko </ci>
<cn cellml:units="millimolar"> 5.4 </cn>
</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/>
<apply>
<plus/>
<ci> Ko </ci>
<apply>
<times/>
<ci> PNa_K </ci>
<ci> Nao </ci>
</apply>
</apply>
<apply>
<plus/>
<ci> Ki </ci>
<apply>
<times/>
<ci> PNa_K </ci>
<ci> Nai </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="i_K_calculation">
<eq/>
<ci> i_K </ci>
<apply>
<times/>
<ci> g_K </ci>
<ci> Xi </ci>
<apply>
<power/>
<ci> X </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_dependent_potassium_current_X_gate">
<variable units="dimensionless" public_interface="out" name="X" initial_value="0.000337175"/>
<variable units="per_millisecond" name="alpha_X"/>
<variable units="per_millisecond" name="beta_X"/>
<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="alpha_X_calculation">
<eq/>
<ci> alpha_X </ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_millisecond"> 0.0000719 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 30.0 </cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> -0.148 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 30.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_X_calculation">
<eq/>
<ci> beta_X </ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_millisecond"> 0.000131 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 30.0 </cn>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> -1.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> 0.0687 </cn>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 30.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="dX_dt">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> X </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> alpha_X </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> X </ci>
</apply>
</apply>
<apply>
<times/>
<ci> beta_X </ci>
<ci> X </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_dependent_potassium_current_Xi_gate">
<variable units="dimensionless" public_interface="out" name="Xi"/>
<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="Xi_calculation">
<eq/>
<ci> Xi </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/>
<ci> V </ci>
<cn cellml:units="millivolt"> 40.0 </cn>
</apply>
<cn cellml:units="millivolt"> 40.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_independent_potassium_current">
<variable units="microA_per_microF" public_interface="out" name="i_K1"/>
<variable units="millivolt" public_interface="out" private_interface="out" name="E_K1"/>
<variable units="milliS_per_microF" name="g_K1"/>
<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="Ko"/>
<variable units="millimolar" public_interface="in" name="Ki"/>
<variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_millimole" public_interface="in" name="F"/>
<variable units="dimensionless" private_interface="in" name="K1_infinity"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="g_K1_calculation">
<eq/>
<ci> g_K1 </ci>
<apply>
<times/>
<cn cellml:units="milliS_per_microF"> 0.75 </cn>
<apply>
<root/>
<apply>
<divide/>
<ci> Ko </ci>
<cn cellml:units="millimolar"> 5.4 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="E_K1_calculation">
<eq/>
<ci> E_K1 </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>
<apply id="i_K1_calculation">
<eq/>
<ci> i_K1 </ci>
<apply>
<times/>
<ci> g_K1 </ci>
<ci> K1_infinity </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K1 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_independent_potassium_current_K1_gate">
<variable units="dimensionless" public_interface="out" name="K1_infinity"/>
<variable units="per_millisecond" name="alpha_K1"/>
<variable units="per_millisecond" name="beta_K1"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="E_K1"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="alpha_K1_calculation">
<eq/>
<ci> alpha_K1 </ci>
<apply>
<divide/>
<cn cellml:units="per_millisecond"> 1.02 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> 0.2385 </cn>
<apply>
<minus/>
<ci> V </ci>
<apply>
<plus/>
<ci> E_K1 </ci>
<cn cellml:units="millivolt"> 59.215 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_K1_calculation">
<eq/>
<ci> beta_K1 </ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="per_millisecond"> 0.49124 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.08032 </cn>
<apply>
<plus/>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K1 </ci>
</apply>
<cn cellml:units="millivolt"> 5.476 </cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.06175 </cn>
<apply>
<minus/>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K1 </ci>
</apply>
<cn cellml:units="millivolt"> 594.31 </cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt"> -0.5143 </cn>
<apply>
<plus/>
<apply>
<minus/>
<ci> V </ci>
<ci> E_K1 </ci>
</apply>
<cn cellml:units="millivolt"> 4.753 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="K1_infinity_calculation">
<eq/>
<ci> K1_infinity </ci>
<apply>
<divide/>
<ci> alpha_K1 </ci>
<apply>
<plus/>
<ci> alpha_K1 </ci>
<ci> beta_K1 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="plateau_potassium_current">
<variable units="microA_per_microF" public_interface="out" name="i_Kp"/>
<variable units="millivolt" name="E_Kp"/>
<variable units="milliS_per_microF" name="g_Kp" initial_value="0.00828"/>
<variable units="dimensionless" name="Kp"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millivolt" public_interface="in" name="E_K1"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="E_Kp_calculation">
<eq/>
<ci> E_Kp </ci>
<ci> E_K1 </ci>
</apply>
<apply id="Kp_calculation">
<eq/>
<ci> Kp </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/>
<cn cellml:units="millivolt"> 7.488 </cn>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 5.98 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="i_Kp_calculation">
<eq/>
<ci> i_Kp </ci>
<apply>
<times/>
<ci> g_Kp </ci>
<ci> Kp </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_Kp </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Na_Ca_exchanger">
<variable units="microA_per_microF" public_interface="out" name="i_NaCa"/>
<variable units="microA_per_microF" name="k_NaCa" initial_value="9000.0"/>
<variable units="millimolar" name="K_mNa" initial_value="87.5"/>
<variable units="millimolar" name="K_mCa" initial_value="1.38"/>
<variable units="dimensionless" name="k_sat" initial_value="0.1"/>
<variable units="dimensionless" name="eta" initial_value="0.35"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_millimole" public_interface="in" name="F"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="micromolar" public_interface="in" name="Cai"/>
<variable units="millimolar" public_interface="in" name="Cao"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="Na_Ca_exchanger">
<eq/>
<ci> i_NaCa </ci>
<apply>
<times/>
<ci> k_NaCa </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<apply>
<power/>
<ci> K_mNa </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<apply>
<power/>
<ci> Nao </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<ci> K_mCa </ci>
<ci> Cao </ci>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<ci> k_sat </ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci> eta </ci>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<exp/>
<apply>
<times/>
<ci> eta </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci> Nai </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> Cao </ci>
</apply>
<apply>
<times/>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci> eta </ci>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci> Nao </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> Cai </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_potassium_pump">
<variable units="microA_per_microF" public_interface="out" name="i_NaK"/>
<variable units="microA_per_microF" name="I_NaK" initial_value="2.47"/>
<variable units="dimensionless" name="f_NaK"/>
<variable units="millimolar" name="K_mNai" initial_value="10.0"/>
<variable units="millimolar" name="K_mKo" initial_value="1.5"/>
<variable units="dimensionless" name="sigma"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_millimole" 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="Ko"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="f_NaK_calculation">
<eq/>
<ci> f_NaK </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.1245 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> -0.1 </cn>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.0365 </cn>
<ci> sigma </ci>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="sigma_calculation">
<eq/>
<ci> sigma </ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<cn cellml:units="dimensionless"> 7.0 </cn>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci> Nao </ci>
<cn cellml:units="dimensionless"> 67.3 </cn>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
</apply>
<apply id="i_NaK_calculation">
<eq/>
<ci> i_NaK </ci>
<apply>
<times/>
<ci> I_NaK </ci>
<ci> f_NaK </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<power/>
<apply>
<divide/>
<ci> K_mNai </ci>
<ci> Nai </ci>
</apply>
<cn cellml:units="dimensionless"> 1.5 </cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<ci> Ko </ci>
<apply>
<plus/>
<ci> Ko </ci>
<ci> K_mKo </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="non_specific_calcium_activated_current">
<variable units="microA_per_microF" public_interface="out" name="i_ns_Ca"/>
<variable units="microA_per_microF" public_interface="out" name="i_ns_Na"/>
<variable units="microA_per_microF" public_interface="out" name="i_ns_K"/>
<variable units="microA_per_microF" name="I_ns_Ca"/>
<variable units="microA_per_microF" name="I_ns_Na"/>
<variable units="microA_per_microF" name="I_ns_K"/>
<variable units="micromolar" name="K_m_ns_Ca" initial_value="1.2"/>
<variable units="cm_per_second" name="Pns_Na" initial_value="0.0"/>
<variable units="cm_per_second" name="Pns_K" initial_value="0.0"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="micromolar" public_interface="in" name="Cai"/>
<variable units="micromolar" public_interface="in" name="Nai"/>
<variable units="micromolar" public_interface="in" name="Ki"/>
<variable units="micromolar" public_interface="in" name="Nao"/>
<variable units="micromolar" public_interface="in" name="Ko"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_millimole" public_interface="in" name="F"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_ns_Na_calculation">
<eq/>
<ci> i_ns_Na </ci>
<apply>
<times/>
<ci> I_ns_Na </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<power/>
<apply>
<divide/>
<ci> K_m_ns_Ca </ci>
<ci> Cai </ci>
</apply>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="i_ns_K_calculation">
<eq/>
<ci> i_ns_K </ci>
<apply>
<times/>
<ci> I_ns_K </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<power/>
<apply>
<divide/>
<ci> K_m_ns_Ca </ci>
<ci> Cai </ci>
</apply>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="i_ns_Ca_calculation">
<eq/>
<ci> i_ns_Ca </ci>
<apply>
<plus/>
<ci> i_ns_Na </ci>
<ci> i_ns_K </ci>
</apply>
</apply>
<apply id="I_ns_Na_calculation">
<eq/>
<ci> I_ns_Na </ci>
<apply>
<times/>
<ci> Pns_Na </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<apply>
<power/>
<ci> F </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.75 </cn>
<ci> Nai </ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.75 </cn>
<ci> Nao </ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="I_ns_K_calculation">
<eq/>
<ci> I_ns_K </ci>
<apply>
<times/>
<ci> Pns_K </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<apply>
<power/>
<ci> F </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.75 </cn>
<ci> Ki </ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.75 </cn>
<ci> Ko </ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sarcolemmal_calcium_pump">
<variable units="microA_per_microF" public_interface="out" name="i_p_Ca"/>
<variable units="micromolar" name="K_mpCa" initial_value="0.5"/>
<variable units="microA_per_microF" name="I_pCa" initial_value="0.575"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="micromolar" public_interface="in" name="Cai"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_p_Ca_calculation">
<eq/>
<ci> i_p_Ca </ci>
<apply>
<times/>
<ci> I_pCa </ci>
<apply>
<divide/>
<ci> Cai </ci>
<apply>
<plus/>
<ci> K_mpCa </ci>
<ci> Cai </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_background_current">
<variable units="microA_per_microF" public_interface="out" name="i_Ca_b"/>
<variable units="milliS_per_microF" name="g_Cab" initial_value="0.006032"/>
<variable units="millivolt" name="E_CaN"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_millimole" public_interface="in" name="F"/>
<variable units="micromolar" public_interface="in" name="Cai"/>
<variable units="millimolar" public_interface="in" name="Cao"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="E_CaN_calculation">
<eq/>
<ci> E_CaN </ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> F </ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci> Cao </ci>
<ci> Cai </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="i_Ca_b_calculation">
<eq/>
<ci> i_Ca_b </ci>
<apply>
<times/>
<ci> g_Cab </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_CaN </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_background_current">
<variable units="microA_per_microF" public_interface="out" name="i_Na_b"/>
<variable units="milliS_per_microF" name="g_Nab" initial_value="0.00141"/>
<variable units="millivolt" name="E_NaN"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millivolt" public_interface="in" name="E_Na"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="E_NaN_calculation">
<eq/>
<ci> E_NaN </ci>
<ci> E_Na </ci>
</apply>
<apply id="i_Na_b_calculation">
<eq/>
<ci> i_Na_b </ci>
<apply>
<times/>
<ci> g_Nab </ci>
<apply>
<minus/>
<ci> V </ci>
<ci> E_NaN </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel">
<variable units="microA_per_microF" public_interface="out" name="i_Ca"/>
<variable units="microA_per_microF" public_interface="out" name="i_Ca_K"/>
<variable units="cm_per_second" name="P_Ca" initial_value="4.32E-3"/>
<variable units="cm_per_second" name="P_K" initial_value="1.11E-11"/>
<variable units="cm_per_second" name="p_k"/>
<variable units="microA_per_microF" name="i_Ca_half" initial_value="-2.16"/>
<variable units="microA_per_microF" name="i_Ca_max"/>
<variable units="dimensionless" name="O" initial_value="4.21370E-19"/>
<variable units="dimensionless" name="O_Ca" initial_value="0.0"/>
<variable units="per_millisecond" name="alpha"/>
<variable units="per_millisecond" name="beta"/>
<variable units="per_millisecond" name="gamma"/>
<variable units="per_millisecond" name="alpha_a"/>
<variable units="per_millisecond" name="beta_b"/>
<variable units="dimensionless" name="a" initial_value="2.0"/>
<variable units="dimensionless" name="b" initial_value="2.0"/>
<variable units="per_millisecond" name="g" initial_value="2.0"/>
<variable units="per_millisecond" name="f" initial_value="0.3"/>
<variable units="per_millisecond" name="g_" initial_value="0.0"/>
<variable units="per_millisecond" name="f_" initial_value="0.0"/>
<variable units="per_millisecond" name="omega" initial_value="0.01"/>
<variable units="dimensionless" name="C0" initial_value="0.998666"/>
<variable units="dimensionless" name="C1" initial_value="1.63594E-4"/>
<variable units="dimensionless" name="C2" initial_value="1.00495E-8"/>
<variable units="dimensionless" name="C3" initial_value="2.74371E-13"/>
<variable units="dimensionless" name="C4" initial_value="2.80908E-18"/>
<variable units="dimensionless" name="C_Ca0" initial_value="0.00272826"/>
<variable units="dimensionless" name="C_Ca1" initial_value="7.67674E-7"/>
<variable units="dimensionless" name="C_Ca2" initial_value="1.88635E-10"/>
<variable units="dimensionless" name="C_Ca3" initial_value="2.06007E-14"/>
<variable units="dimensionless" name="C_Ca4" initial_value="8.43672E-19"/>
<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="Ca_SS"/>
<variable units="millimolar" public_interface="in" name="Cao"/>
<variable units="millimolar" public_interface="in" name="Ko"/>
<variable units="millimolar" public_interface="in" name="Ki"/>
<variable units="joule_per_mole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_millimole" public_interface="in" name="F"/>
<variable units="dimensionless" private_interface="in" name="y"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_Ca_calculation">
<eq/>
<ci> i_Ca </ci>
<apply>
<times/>
<ci> i_Ca_max </ci>
<ci> y </ci>
<apply>
<plus/>
<ci> O </ci>
<ci> O_Ca </ci>
</apply>
</apply>
</apply>
<apply id="i_Ca_max_calculation">
<eq/>
<ci> i_Ca_max </ci>
<apply>
<times/>
<apply>
<power/>
<ci> P_Ca </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<apply>
<power/>
<ci> F </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.001 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.341 </cn>
<ci> Cao </ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="i_Ca_K_calculation">
<eq/>
<ci> i_Ca_K </ci>
<apply>
<times/>
<ci> p_k </ci>
<ci> y </ci>
<apply>
<plus/>
<ci> O </ci>
<ci> O_Ca </ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<apply>
<power/>
<ci> F </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci> Ki </ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<ci> Ko </ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="p_k_calculation">
<eq/>
<ci> p_k </ci>
<apply>
<divide/>
<ci> P_K </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> i_Ca_max </ci>
<ci> i_Ca_half </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="alpha_calculation">
<eq/>
<ci> alpha </ci>
<apply>
<times/>
<cn cellml:units="per_millisecond"> 0.4 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 12.0 </cn>
</apply>
<cn cellml:units="millivolt"> 10.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_calculation">
<eq/>
<ci> beta </ci>
<apply>
<times/>
<cn cellml:units="per_millisecond"> 0.05 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 12.0 </cn>
</apply>
</apply>
<cn cellml:units="millivolt"> 13.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="alpha_a_calculation">
<eq/>
<ci> alpha_a </ci>
<apply>
<times/>
<ci> alpha </ci>
<ci> a </ci>
</apply>
</apply>
<apply id="beta_b_calculation">
<eq/>
<ci> beta_b </ci>
<apply>
<divide/>
<ci> beta </ci>
<ci> b </ci>
</apply>
</apply>
<apply id="gamma_calculation">
<eq/>
<ci> gamma </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.5625 </cn>
<ci> Ca_SS </ci>
</apply>
</apply>
<apply id="C0_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C0 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta </ci>
<ci> C1 </ci>
</apply>
<apply>
<times/>
<ci> omega </ci>
<ci> C_Ca0 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> alpha </ci>
</apply>
<ci> gamma </ci>
</apply>
<ci> C0 </ci>
</apply>
</apply>
</apply>
<apply id="C1_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C1 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> alpha </ci>
<ci> C0 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> beta </ci>
<ci> C2 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega </ci>
<ci> b </ci>
</apply>
<ci> C_Ca1 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci> beta </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> alpha </ci>
</apply>
<apply>
<times/>
<ci> gamma </ci>
<ci> a </ci>
</apply>
</apply>
<ci> C1 </ci>
</apply>
</apply>
</apply>
<apply id="C2_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C2 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> alpha </ci>
<ci> C1 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta </ci>
<ci> C3 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega </ci>
<apply>
<power/>
<ci> b </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<ci> C_Ca2 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha </ci>
</apply>
<apply>
<times/>
<ci> gamma </ci>
<apply>
<power/>
<ci> a </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<ci> C2 </ci>
</apply>
</apply>
</apply>
<apply id="C3_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C3 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha </ci>
<ci> C2 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> beta </ci>
<ci> C4 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega </ci>
<apply>
<power/>
<ci> b </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
<ci> C_Ca3 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> alpha </ci>
<apply>
<times/>
<ci> gamma </ci>
<apply>
<power/>
<ci> a </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
<ci> C3 </ci>
</apply>
</apply>
</apply>
<apply id="C4_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C4 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> alpha </ci>
<ci> C3 </ci>
</apply>
<apply>
<times/>
<ci> g </ci>
<ci> O </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega </ci>
<apply>
<power/>
<ci> b </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
<ci> C_Ca4 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> f </ci>
<apply>
<times/>
<ci> gamma </ci>
<apply>
<power/>
<ci> a </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
</apply>
<ci> C4 </ci>
</apply>
</apply>
</apply>
<apply id="O_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> O </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> f </ci>
<ci> C4 </ci>
</apply>
<apply>
<times/>
<ci> g </ci>
<ci> O </ci>
</apply>
</apply>
</apply>
<apply id="C_Ca0_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C_Ca0 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta_b </ci>
<ci> C_Ca1 </ci>
</apply>
<apply>
<times/>
<ci> gamma </ci>
<ci> C_Ca0 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> alpha_a </ci>
</apply>
<ci> omega </ci>
</apply>
<ci> C_Ca0 </ci>
</apply>
</apply>
</apply>
<apply id="C_Ca1_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C_Ca1 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> alpha_a </ci>
<ci> C_Ca0 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> beta_b </ci>
<ci> C_Ca2 </ci>
</apply>
<apply>
<times/>
<ci> gamma </ci>
<ci> a </ci>
<ci> C1 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci> beta_b </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> alpha_a </ci>
</apply>
<apply>
<divide/>
<ci> omega </ci>
<ci> b </ci>
</apply>
</apply>
<ci> C_Ca1 </ci>
</apply>
</apply>
</apply>
<apply id="C_Ca2_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C_Ca2 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> alpha_a </ci>
<ci> C_Ca1 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta_b </ci>
<ci> C_Ca3 </ci>
</apply>
<apply>
<times/>
<ci> gamma </ci>
<apply>
<power/>
<ci> a </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> C2 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta_b </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_a </ci>
</apply>
<apply>
<divide/>
<ci> omega </ci>
<apply>
<power/>
<ci> b </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<ci> C_Ca2 </ci>
</apply>
</apply>
</apply>
<apply id="C_Ca3_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C_Ca3 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_a </ci>
<ci> C_Ca2 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> beta_b </ci>
<ci> C_Ca4 </ci>
</apply>
<apply>
<times/>
<ci> gamma </ci>
<apply>
<power/>
<ci> a </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> C3 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta_b </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> alpha_a </ci>
<apply>
<divide/>
<ci> omega </ci>
<apply>
<power/>
<ci> b </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
<ci> C_Ca3 </ci>
</apply>
</apply>
</apply>
<apply id="C_Ca4_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> C_Ca4 </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> alpha_a </ci>
<ci> C_Ca3 </ci>
</apply>
<apply>
<times/>
<ci> g_ </ci>
<ci> O_Ca </ci>
</apply>
<apply>
<times/>
<ci> gamma </ci>
<apply>
<power/>
<ci> a </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> C4 </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta_b </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> f_ </ci>
<apply>
<divide/>
<ci> omega </ci>
<apply>
<power/>
<ci> b </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
</apply>
<ci> C_Ca4 </ci>
</apply>
</apply>
</apply>
<apply id="O_Ca_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> O_Ca </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> f_ </ci>
<ci> C_Ca4 </ci>
</apply>
<apply>
<times/>
<ci> g_ </ci>
<ci> O_Ca </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_y_gate">
<variable units="dimensionless" public_interface="out" name="y" initial_value="0.862010"/>
<variable units="dimensionless" name="y_infinity"/>
<variable units="dimensionless" 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_diff_eq">
<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>
<apply id="y_infinity_calculation">
<eq/>
<ci> y_infinity </ci>
<apply>
<plus/>
<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"> 55.0 </cn>
</apply>
<cn cellml:units="millivolt"> 7.5 </cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 0.1 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 21.0 </cn>
</apply>
<cn cellml:units="millivolt"> 6.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_y_calculation">
<eq/>
<ci> tau_y </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 20.0 </cn>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 600.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"> 9.5 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="RyR_channel_states">
<variable units="dimensionless" public_interface="out" name="P_O1" initial_value="0.428917E-3"/>
<variable units="dimensionless" public_interface="out" name="P_O2" initial_value="0.487455E-9"/>
<variable units="dimensionless" name="P_C1" initial_value="0.637931"/>
<variable units="dimensionless" name="P_C2" initial_value="0.361643"/>
<variable units="per_micromolar4_per_millisecond" name="k_a_plus" initial_value="0.01215"/>
<variable units="per_millisecond" name="k_a_minus" initial_value="0.1425"/>
<variable units="per_micromolar3_per_millisecond" name="k_b_plus" initial_value="0.00405"/>
<variable units="per_millisecond" name="k_b_minus" initial_value="1.930"/>
<variable units="per_millisecond" name="k_c_plus" initial_value="0.010"/>
<variable units="per_millisecond" name="k_c_minus" initial_value="0.0008"/>
<variable units="dimensionless" name="n" initial_value="4.0"/>
<variable units="dimensionless" name="m" initial_value="3.0"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="millimolar" public_interface="in" name="Ca_SS"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="P_C1_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> P_C1 </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<ci> k_a_plus </ci>
</apply>
<apply>
<power/>
<ci> Ca_SS </ci>
<ci> n </ci>
</apply>
<ci> P_C1 </ci>
</apply>
<apply>
<times/>
<ci> k_a_minus </ci>
<ci> P_O1 </ci>
</apply>
</apply>
</apply>
<apply id="P_O1_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> P_O1 </ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci> k_a_plus </ci>
<apply>
<power/>
<ci> Ca_SS </ci>
<ci> n </ci>
</apply>
<ci> P_C1 </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> k_a_minus </ci>
<ci> P_O1 </ci>
</apply>
<apply>
<times/>
<ci> k_b_plus </ci>
<apply>
<power/>
<ci> Ca_SS </ci>
<ci> m </ci>
</apply>
<ci> P_O1 </ci>
</apply>
<apply>
<times/>
<ci> k_c_plus </ci>
<ci> P_O1 </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> k_b_minus </ci>
<ci> P_O2 </ci>
</apply>
<apply>
<times/>
<ci> k_c_minus </ci>
<ci> P_C2 </ci>
</apply>
</apply>
</apply>
<apply id="P_O2_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> P_O2 </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_b_plus </ci>
<apply>
<power/>
<ci> Ca_SS </ci>
<ci> m </ci>
</apply>
<ci> P_O1 </ci>
</apply>
<apply>
<times/>
<ci> k_b_minus </ci>
<ci> P_O2 </ci>
</apply>
</apply>
</apply>
<apply id="P_C2_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> P_C2 </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_c_plus </ci>
<ci> P_O1 </ci>
</apply>
<apply>
<times/>
<ci> k_c_minus </ci>
<ci> P_C2 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="SERCA_pump">
<variable units="micromolar_per_millisecond" public_interface="out" name="J_up"/>
<variable units="micromolar" name="k_fb" initial_value="0.24"/>
<variable units="millimolar" name="k_rb" initial_value="1.64"/>
<variable units="micromolar_per_second" name="fb" initial_value="0.292"/>
<variable units="micromolar_per_second" name="rb" initial_value="0.311"/>
<variable units="micromolar_per_second" name="Vmax_f"/>
<variable units="micromolar_per_second" name="Vmax_r"/>
<variable units="dimensionless" name="n_fb" initial_value="1.4"/>
<variable units="dimensionless" name="n_rb"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="micromolar" public_interface="in" name="Ca_NSR"/>
<variable units="micromolar" public_interface="in" name="Cai"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="J_up_calculation">
<eq/>
<ci> J_up </ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci> Vmax_f </ci>
<ci> fb </ci>
</apply>
<apply>
<times/>
<ci> Vmax_r </ci>
<ci> rb </ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> fb </ci>
<ci> rb </ci>
</apply>
</apply>
</apply>
<apply id="fb_calculation">
<eq/>
<ci> fb </ci>
<apply>
<power/>
<apply>
<divide/>
<ci> Cai </ci>
<ci> k_fb </ci>
</apply>
<ci> n_fb </ci>
</apply>
</apply>
<apply id="rb_calculation">
<eq/>
<ci> rb </ci>
<apply>
<power/>
<apply>
<divide/>
<ci> Ca_NSR </ci>
<ci> k_rb </ci>
</apply>
<ci> n_rb </ci>
</apply>
</apply>
</math>
</component>
<component name="intracellular_Ca_fluxes">
<variable units="micromolar_per_millisecond" public_interface="out" name="J_rel"/>
<variable units="micromolar_per_millisecond" public_interface="out" name="J_trpn"/>
<variable units="micromolar_per_millisecond" public_interface="out" name="J_tr"/>
<variable units="micromolar_per_millisecond" public_interface="out" name="J_xfer"/>
<variable units="per_millisecond" name="v1" initial_value="3.6"/>
<variable units="millisecond" name="tau_tr" initial_value="0.5747"/>
<variable units="millisecond" name="tau_xfer" initial_value="9.09"/>
<variable units="micromolar" public_interface="in" name="HTRPN_tot"/>
<variable units="micromolar" public_interface="in" name="LTRPN_tot"/>
<variable units="millimolar" public_interface="in" name="HTRPNCa"/>
<variable units="millimolar" public_interface="in" name="LTRPNCa"/>
<variable units="per_micromolar_per_second" public_interface="in" name="k_htrpn_plus"/>
<variable units="per_second" public_interface="in" name="k_htrpn_minus"/>
<variable units="per_micromolar_per_second" public_interface="in" name="k_ltrpn_plus"/>
<variable units="per_micromolar_per_second" public_interface="in" name="k_ltrpn_minus"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="dimensionless" public_interface="in" name="P_O1"/>
<variable units="dimensionless" public_interface="in" name="P_O2"/>
<variable units="micromolar" public_interface="in" name="Cai"/>
<variable units="millimolar" public_interface="in" name="Ca_SS"/>
<variable units="micromolar" public_interface="in" name="Ca_NSR"/>
<variable units="millimolar" public_interface="in" name="Ca_JSR"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="J_rel_calculation">
<eq/>
<ci> J_rel </ci>
<apply>
<times/>
<ci> v1 </ci>
<apply>
<plus/>
<ci> P_O1 </ci>
<ci> P_O2 </ci>
</apply>
<apply>
<minus/>
<ci> Ca_JSR </ci>
<ci> Ca_SS </ci>
</apply>
</apply>
</apply>
<apply id="J_tr_calculation">
<eq/>
<ci> J_tr </ci>
<apply>
<divide/>
<apply>
<minus/>
<ci> Ca_NSR </ci>
<ci> Ca_JSR </ci>
</apply>
<ci> tau_tr </ci>
</apply>
</apply>
<apply id="J_xfer_calculation">
<eq/>
<ci> J_xfer </ci>
<apply>
<divide/>
<apply>
<minus/>
<ci> Ca_SS </ci>
<ci> Cai </ci>
</apply>
<ci> tau_xfer </ci>
</apply>
</apply>
<apply id="J_trpn_calculation">
<eq/>
<ci> J_trpn </ci>
<apply>
<plus/>
<apply>
<times/>
<ci> k_htrpn_plus </ci>
<ci> Cai </ci>
<apply>
<minus/>
<apply>
<minus/>
<ci> HTRPN_tot </ci>
<ci> HTRPNCa </ci>
</apply>
<apply>
<times/>
<ci> k_htrpn_minus </ci>
<ci> HTRPNCa </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> k_ltrpn_plus </ci>
<ci> Cai </ci>
<apply>
<minus/>
<apply>
<minus/>
<ci> LTRPN_tot </ci>
<ci> LTRPNCa </ci>
</apply>
<apply>
<times/>
<ci> k_ltrpn_minus </ci>
<ci> LTRPNCa </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="intracellular_ionic_concentrations">
<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="Cai"/>
<variable units="millimolar" public_interface="out" name="Cao" initial_value="1.8"/>
<variable units="millimolar" public_interface="out" name="Ca_SS"/>
<variable units="millimolar" public_interface="out" name="Ca_NSR"/>
<variable units="millimolar" public_interface="out" name="Ca_JSR"/>
<variable units="dimensionless" name="betai"/>
<variable units="dimensionless" name="beta_SS"/>
<variable units="dimensionless" name="beta_JSR"/>
<variable units="micromolar" name="K_mCMDN" initial_value="2.38"/>
<variable units="millimolar" name="K_mCSQN" initial_value="0.8"/>
<variable units="millimolar" name="CSQN_tot" initial_value="15.0"/>
<variable units="micromolar" name="CMDN_tot" initial_value="50.0"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="micromolar_per_millisecond" public_interface="in" name="J_rel"/>
<variable units="micromolar_per_millisecond" public_interface="in" name="J_up"/>
<variable units="micromolar_per_millisecond" public_interface="in" name="J_tr"/>
<variable units="micromolar_per_millisecond" public_interface="in" name="J_xfer"/>
<variable units="micromolar_per_millisecond" public_interface="in" name="J_trpn"/>
<variable units="coulomb_per_millimole" public_interface="in" name="F"/>
<variable units="microA_per_microF" public_interface="in" name="i_Na"/>
<variable units="microA_per_microF" public_interface="in" name="i_Na_b"/>
<variable units="microA_per_microF" public_interface="in" name="i_ns_Na"/>
<variable units="microA_per_microF" public_interface="in" name="i_NaCa"/>
<variable units="microA_per_microF" public_interface="in" name="i_NaK"/>
<variable units="microA_per_microF" public_interface="in" name="i_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="i_p_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="i_Ca_b"/>
<variable units="microA_per_microF" public_interface="in" name="i_Ca_K"/>
<variable units="microA_per_microF" public_interface="in" name="i_K"/>
<variable units="microA_per_microF" public_interface="in" name="i_K1"/>
<variable units="microA_per_microF" public_interface="in" name="i_Kp"/>
<variable units="microA_per_microF" public_interface="in" name="i_ns_K"/>
<variable units="cm2" public_interface="in" name="A_cap"/>
<variable units="micro_litre" public_interface="in" name="V_myo"/>
<variable units="micro_litre" public_interface="in" name="V_JSR"/>
<variable units="micro_litre" public_interface="in" name="V_NSR"/>
<variable units="micro_litre" public_interface="in" name="V_SS"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="betai_calculation">
<eq/>
<ci> betai </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<apply>
<times/>
<ci> CMDN_tot </ci>
<ci> K_mCMDN </ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci> K_mCMDN </ci>
<ci> Cai </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_SS_calculation">
<eq/>
<ci> beta_SS </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<apply>
<times/>
<ci> CMDN_tot </ci>
<ci> K_mCMDN </ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci> K_mCMDN </ci>
<ci> Ca_SS </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_JSR_calculation">
<eq/>
<ci> beta_JSR </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<apply>
<times/>
<ci> CSQN_tot </ci>
<ci> K_mCSQN </ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci> K_mCSQN </ci>
<ci> Ca_JSR </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="calcium_internal_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Cai </ci>
</apply>
<apply>
<times/>
<ci> betai </ci>
<apply>
<minus/>
<apply>
<plus/>
<ci> J_xfer </ci>
<ci> J_up </ci>
<ci> J_trpn </ci>
<apply>
<times/>
<apply>
<plus/>
<ci> i_Ca_b </ci>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> i_NaCa </ci>
</apply>
</apply>
<ci> i_p_Ca </ci>
</apply>
<apply>
<divide/>
<ci> A_cap </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> V_myo </ci>
<ci> F </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="calcium_subspace_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Ca_SS </ci>
</apply>
<apply>
<times/>
<ci> beta_SS </ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<ci> J_rel </ci>
<apply>
<divide/>
<ci> V_JSR </ci>
<ci> V_SS </ci>
</apply>
</apply>
<apply>
<times/>
<ci> J_xfer </ci>
<apply>
<divide/>
<ci> V_myo </ci>
<ci> V_SS </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> i_Ca </ci>
<apply>
<divide/>
<ci> A_cap </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> V_SS </ci>
<ci> F </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="calcium_JSR_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Ca_JSR </ci>
</apply>
<apply>
<times/>
<ci> beta_JSR </ci>
<apply>
<minus/>
<ci> J_tr </ci>
<ci> J_rel </ci>
</apply>
</apply>
</apply>
<apply id="calcium_NSR_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Ca_NSR </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> J_up </ci>
<apply>
<divide/>
<ci> V_myo </ci>
<ci> V_NSR </ci>
</apply>
</apply>
<apply>
<times/>
<ci> J_tr </ci>
<apply>
<divide/>
<ci> V_JSR </ci>
<ci> V_NSR </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="sodium_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Nai </ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<apply>
<plus/>
<ci> i_Na </ci>
<ci> i_Na_b </ci>
<ci> i_ns_Na </ci>
<apply>
<times/>
<ci> i_NaCa </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<apply>
<times/>
<ci> i_NaK </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<ci> A_cap </ci>
<apply>
<times/>
<ci> V_myo </ci>
<ci> F </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="potassium_internal_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Ki </ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<apply>
<plus/>
<ci> i_Ca_K </ci>
<ci> i_K </ci>
<ci> i_K1 </ci>
<ci> i_Kp </ci>
<ci> i_ns_K </ci>
<apply>
<minus/>
<apply>
<times/>
<ci> i_NaK </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<ci> A_cap </ci>
<apply>
<times/>
<ci> V_myo </ci>
<ci> F </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="troponin">
<variable units="per_micromolar_per_second" public_interface="out" name="k_htrpn_plus" initial_value="1.0E8"/>
<variable units="per_second" public_interface="out" name="k_htrpn_minus" initial_value="0.33"/>
<variable units="per_micromolar_per_second" public_interface="out" name="k_ltrpn_plus" initial_value="4.0E7"/>
<variable units="per_micromolar_per_second" public_interface="out" name="k_ltrpn_minus" initial_value="40.0"/>
<variable units="micromolar" public_interface="out" name="HTRPN_tot" initial_value="140.0"/>
<variable units="micromolar" public_interface="out" name="LTRPN_tot" initial_value="70.0"/>
<variable units="millimolar" public_interface="out" name="HTRPNCa" initial_value="0.132944"/>
<variable units="millimolar" public_interface="out" name="LTRPNCa" initial_value="0.4068232"/>
<variable units="N_per_mm2" public_interface="in" name="Force_norm"/>
<variable units="millimolar" public_interface="in" name="Cai"/>
<variable units="millisecond" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="HTRPNCa_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> HTRPNCa </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_htrpn_plus </ci>
<ci> Cai </ci>
<apply>
<minus/>
<ci> HTRPN_tot </ci>
<ci> HTRPNCa </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k_htrpn_minus </ci>
<ci> HTRPNCa </ci>
</apply>
</apply>
</apply>
<apply id="LTRPNCa_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> LTRPNCa </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_ltrpn_plus </ci>
<ci> Cai </ci>
<apply>
<minus/>
<ci> LTRPN_tot </ci>
<ci> LTRPNCa </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k_ltrpn_minus </ci>
<ci> LTRPNCa </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 0.333 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.667 </cn>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> Force_norm </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="tropomyosin_cross_bridges">
<variable units="per_second" public_interface="out" name="f_01"/>
<variable units="per_second" public_interface="out" name="f_12"/>
<variable units="per_second" public_interface="out" name="f_23"/>
<variable units="per_second" public_interface="out" name="g_XB" initial_value="30.0"/>
<variable units="micrometre" public_interface="out" name="SL" initial_value="2.15"/>
<variable units="dimensionless" public_interface="out" name="P1" initial_value="1.47026E-6"/>
<variable units="dimensionless" public_interface="out" name="P2" initial_value="1.83527E-6"/>
<variable units="dimensionless" public_interface="out" name="P3" initial_value="1.06786E-6"/>
<variable units="dimensionless" public_interface="out" name="N1" initial_value="3.88506E-6"/>
<variable units="per_second" name="f_XB" initial_value="10.0"/>
<variable units="per_second" name="g_XB_SL"/>
<variable units="per_second" name="g_01_SL"/>
<variable units="per_second" name="g_12_SL"/>
<variable units="per_second" name="g_23_SL"/>
<variable units="dimensionless" name="N0" initial_value="0.99999027"/>
<variable units="dimensionless" name="P0" initial_value="1.46758E-6"/>
<variable units="micrometre" name="SL_norm"/>
<variable units="micrometre" name="N_trop"/>
<variable units="micromolar" name="K_trop"/>
<variable units="micromolar" name="k_trop_pn"/>
<variable units="micromolar" name="k_trop_np"/>
<variable units="micromolar" name="K_trop_half"/>
<variable units="micromolar" name="K_trop_Ca"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="per_micromolar_per_second" public_interface="in" name="k_ltrpn_plus"/>
<variable units="per_micromolar_per_second" public_interface="in" name="k_ltrpn_minus"/>
<variable units="micromolar" public_interface="in" name="LTRPN_tot"/>
<variable units="millimolar" public_interface="in" name="LTRPNCa"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="f_01_calculation">
<eq/>
<ci> f_01 </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> f_XB </ci>
</apply>
</apply>
<apply id="f_12_calculation">
<eq/>
<ci> f_12 </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 10.0 </cn>
<ci> f_XB </ci>
</apply>
</apply>
<apply id="f_23_calculation">
<eq/>
<ci> f_23 </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 7.0 </cn>
<ci> f_XB </ci>
</apply>
</apply>
<apply id="g_01_SL_calculation">
<eq/>
<ci> g_01_SL </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> g_XB_SL </ci>
</apply>
</apply>
<apply id="g_12_SL_calculation">
<eq/>
<ci> g_12_SL </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> g_XB_SL </ci>
</apply>
</apply>
<apply id="g_23_SL_calculation">
<eq/>
<ci> g_23_SL </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB_SL </ci>
</apply>
</apply>
<apply id="g_XB_SL_calculation">
<eq/>
<ci> g_XB_SL </ci>
<apply>
<times/>
<ci> g_XB </ci>
<apply>
<plus/>
<cn cellml:units="micrometre"> 1.0 </cn>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="micrometre"> 1.0 </cn>
<ci> SL_norm </ci>
</apply>
<cn cellml:units="dimensionless"> 1.6 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="SL_norm_calculation">
<eq/>
<ci> SL_norm </ci>
<apply>
<divide/>
<apply>
<minus/>
<ci> SL </ci>
<cn cellml:units="micrometre"> 1.7 </cn>
</apply>
<cn cellml:units="micrometre"> 0.7 </cn>
</apply>
</apply>
<apply id="k_trop_np_calculation">
<eq/>
<ci> k_trop_np </ci>
<apply>
<times/>
<ci> k_trop_pn </ci>
<apply>
<power/>
<apply>
<divide/>
<apply>
<divide/>
<ci> LTRPNCa </ci>
<ci> LTRPN_tot </ci>
</apply>
<ci> K_trop_half </ci>
</apply>
<ci> N_trop </ci>
</apply>
</apply>
</apply>
<apply id="N_trop_calculation">
<eq/>
<ci> N_trop </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 3.5 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.5 </cn>
<ci> SL_norm </ci>
</apply>
</apply>
</apply>
<apply id="K_trop_half_calculation">
<eq/>
<ci> K_trop_half </ci>
<apply>
<power/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> K_trop_Ca </ci>
<apply>
<minus/>
<cn cellml:units="micromolar"> 1.7 </cn>
<apply>
<times/>
<cn cellml:units="micromolar"> 0.9 </cn>
<ci> SL_norm </ci>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> -1.0 </cn>
</apply>
</apply>
<apply id="K_trop_Ca_calculation">
<eq/>
<ci> K_trop_Ca </ci>
<apply>
<divide/>
<ci> k_ltrpn_minus </ci>
<ci> k_ltrpn_plus </ci>
</apply>
</apply>
<apply id="N0_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> N0 </ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci> k_trop_pn </ci>
<ci> P0 </ci>
</apply>
<apply>
<times/>
<ci> k_trop_np </ci>
<ci> N0 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> g_01_SL </ci>
<ci> N1 </ci>
</apply>
</apply>
</apply>
<apply id="P0_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> P0 </ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<apply>
<plus/>
<ci> k_trop_pn </ci>
<ci> f_01 </ci>
</apply>
<ci> P0 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k_trop_np </ci>
<ci> N0 </ci>
</apply>
<apply>
<times/>
<ci> g_01_SL </ci>
<ci> P1 </ci>
</apply>
</apply>
</apply>
<apply id="P1_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> P1 </ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<apply>
<plus/>
<ci> k_trop_pn </ci>
<ci> f_12 </ci>
<ci> g_01_SL </ci>
</apply>
<ci> P1 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k_trop_np </ci>
<ci> N1 </ci>
</apply>
<apply>
<times/>
<ci> k_trop_np </ci>
<ci> N1 </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<ci> P0 </ci>
</apply>
<apply>
<times/>
<ci> g_12_SL </ci>
<ci> P2 </ci>
</apply>
</apply>
</apply>
<apply id="P2_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> P2 </ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<apply>
<plus/>
<ci> f_23 </ci>
<ci> g_12_SL </ci>
</apply>
<ci> P2 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> f_12 </ci>
<ci> P1 </ci>
</apply>
<apply>
<times/>
<ci> g_23_SL </ci>
<ci> P3 </ci>
</apply>
</apply>
</apply>
<apply id="P3_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> P3 </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<ci> g_23_SL </ci>
</apply>
<ci> P3 </ci>
</apply>
<apply>
<times/>
<ci> f_23 </ci>
<ci> P2 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="force_computation">
<variable units="N_per_mm2" public_interface="out" name="Force_norm"/>
<variable units="N_per_mm2" name="Force"/>
<variable units="N_per_mm2" name="Force_max"/>
<variable units="N_per_mm2" name="zeta" initial_value="0.1"/>
<variable units="dimensionless" name="phi_SL"/>
<variable units="per_second" name="P1_max"/>
<variable units="per_second" name="P2_max"/>
<variable units="per_second" name="P3_max"/>
<variable units="millisecond" public_interface="in" name="time"/>
<variable units="micrometre" public_interface="in" name="SL"/>
<variable units="dimensionless" public_interface="in" name="P1"/>
<variable units="dimensionless" public_interface="in" name="P2"/>
<variable units="dimensionless" public_interface="in" name="P3"/>
<variable units="dimensionless" public_interface="in" name="N1"/>
<variable units="per_second" public_interface="in" name="f_01"/>
<variable units="per_second" public_interface="in" name="f_12"/>
<variable units="per_second" public_interface="in" name="f_23"/>
<variable units="per_second" public_interface="in" name="g_XB"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="Force_calculation">
<eq/>
<ci> Force </ci>
<apply>
<times/>
<ci> zeta </ci>
<ci> Force_norm </ci>
</apply>
</apply>
<apply id="Force_norm_calculation">
<eq/>
<ci> Force_norm </ci>
<apply>
<times/>
<ci> phi_SL </ci>
<apply>
<divide/>
<apply>
<plus/>
<ci> P1 </ci>
<ci> N1 </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> P2 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> P3 </ci>
</apply>
</apply>
<ci> Force_max </ci>
</apply>
</apply>
</apply>
<apply id="Force_max_calculation">
<eq/>
<ci> Force_max </ci>
<apply>
<plus/>
<ci> P1_max </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> P2_max </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> P3_max </ci>
</apply>
</apply>
</apply>
<apply id="phi_SL_calculation">
<eq/>
<ci> phi_SL </ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<minus/>
<ci> SL </ci>
<cn cellml:units="micrometre"> 0.6 </cn>
</apply>
<cn cellml:units="micrometre"> 1.4 </cn>
</apply>
<apply>
<and/>
<apply>
<lt/>
<ci> SL </ci>
<cn cellml:units="micrometre"> 2.0 </cn>
</apply>
<apply>
<gt/>
<ci> SL </ci>
<cn cellml:units="micrometre"> 1.7 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<and/>
<apply>
<lt/>
<ci> SL </ci>
<cn cellml:units="micrometre"> 2.2 </cn>
</apply>
<apply>
<gt/>
<ci> SL </ci>
<cn cellml:units="micrometre"> 2.0 </cn>
</apply>
</apply>
</piece>
<piece>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="micrometre"> 3.6 </cn>
<ci> SL </ci>
</apply>
<cn cellml:units="micrometre"> 1.4 </cn>
</apply>
<apply>
<and/>
<apply>
<lt/>
<ci> SL </ci>
<cn cellml:units="micrometre"> 2.3 </cn>
</apply>
<apply>
<gt/>
<ci> SL </ci>
<cn cellml:units="micrometre"> 2.2 </cn>
</apply>
</apply>
</piece>
</piecewise>
</apply>
<apply id="P1_max_calculation">
<eq/>
<ci> P1_max </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> f_01 </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<ci> f_12 </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<ci> f_12 </ci>
<ci> f_23 </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="P2_max_calculation">
<eq/>
<ci> P2_max </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> f_01 </ci>
<ci> f_12 </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<ci> f_12 </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<ci> f_12 </ci>
<ci> f_23 </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="P3_max_calculation">
<eq/>
<ci> P3_max </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> f_01 </ci>
<ci> f_12 </ci>
<ci> f_23 </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> g_XB </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<ci> f_12 </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> g_XB </ci>
</apply>
<apply>
<times/>
<ci> f_01 </ci>
<ci> f_12 </ci>
<ci> f_23 </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="fast_sodium_current">
<component_ref component="fast_sodium_current_m_gate"/>
<component_ref component="fast_sodium_current_h_gate"/>
<component_ref component="fast_sodium_current_j_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_y_gate"/>
</component_ref>
<component_ref component="time_dependent_potassium_current">
<component_ref component="time_dependent_potassium_current_X_gate"/>
<component_ref component="time_dependent_potassium_current_Xi_gate"/>
</component_ref>
<component_ref component="Na_Ca_exchanger"/>
<component_ref component="time_independent_potassium_current">
<component_ref component="time_independent_potassium_current_K1_gate"/>
</component_ref>
<component_ref component="plateau_potassium_current"/>
<component_ref component="sarcolemmal_calcium_pump"/>
<component_ref component="sodium_background_current"/>
<component_ref component="calcium_background_current"/>
<component_ref component="sodium_potassium_pump"/>
<component_ref component="non_specific_calcium_activated_current"/>
<component_ref component="intracellular_Ca_fluxes"/>
<component_ref component="SERCA_pump"/>
<component_ref component="RyR_channel_states"/>
<component_ref component="intracellular_ionic_concentrations"/>
<component_ref component="troponin"/>
<component_ref component="tropomyosin_cross_bridges"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="fast_sodium_current">
<component_ref component="fast_sodium_current_m_gate"/>
<component_ref component="fast_sodium_current_h_gate"/>
<component_ref component="fast_sodium_current_j_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_y_gate"/>
</component_ref>
<component_ref component="time_dependent_potassium_current">
<component_ref component="time_dependent_potassium_current_X_gate"/>
<component_ref component="time_dependent_potassium_current_Xi_gate"/>
</component_ref>
<component_ref component="time_independent_potassium_current">
<component_ref component="time_independent_potassium_current_K1_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="fast_sodium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="L_type_Ca_channel"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="time_dependent_potassium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="Na_Ca_exchanger"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="time_independent_potassium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="plateau_potassium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sarcolemmal_calcium_pump"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sodium_background_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="calcium_background_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="sodium_potassium_pump"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="non_specific_calcium_activated_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="RyR_channel_states"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="SERCA_pump"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="intracellular_Ca_fluxes"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="intracellular_ionic_concentrations"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="troponin"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="tropomyosin_cross_bridges"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="force_computation"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<!--
Several variables are passed between the "membrane" and its sub-components.
-->
<connection>
<map_components component_2="fast_sodium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Na" variable_1="i_Na"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Ca" variable_1="i_Ca"/>
<map_variables variable_2="i_Ca_K" variable_1="i_Ca_K"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="time_dependent_potassium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_K" variable_1="i_K"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="Na_Ca_exchanger" component_1="membrane"/>
<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="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_K1" variable_1="i_K1"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="plateau_potassium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Kp" variable_1="i_Kp"/>
</connection>
<connection>
<map_components component_2="sarcolemmal_calcium_pump" component_1="membrane"/>
<map_variables variable_2="i_p_Ca" variable_1="i_p_Ca"/>
</connection>
<connection>
<map_components component_2="sodium_background_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Na_b" variable_1="i_Na_b"/>
</connection>
<connection>
<map_components component_2="calcium_background_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_Ca_b" variable_1="i_Ca_b"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="sodium_potassium_pump" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_NaK" variable_1="i_NaK"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="non_specific_calcium_activated_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="i_ns_Ca" variable_1="i_ns_Ca"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="membrane"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<!-- Several variables are passed between the sibling components. -->
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="fast_sodium_current"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="i_Na" variable_1="i_Na"/>
</connection>
<connection>
<map_components component_2="sodium_background_current" component_1="fast_sodium_current"/>
<map_variables variable_2="E_Na" variable_1="E_Na"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="L_type_Ca_channel"/>
<map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
<map_variables variable_2="i_Ca" variable_1="i_Ca"/>
<map_variables variable_2="i_Ca_K" variable_1="i_Ca_K"/>
<map_variables variable_2="Ki" variable_1="Ki"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
<map_variables variable_2="Cao" variable_1="Cao"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="time_dependent_potassium_current"/>
<map_variables variable_2="Ki" variable_1="Ki"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="i_K" variable_1="i_K"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="Na_Ca_exchanger"/>
<map_variables variable_2="Cao" variable_1="Cao"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
<map_variables variable_2="i_NaCa" variable_1="i_NaCa"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="time_independent_potassium_current"/>
<map_variables variable_2="Ki" variable_1="Ki"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
<map_variables variable_2="i_K1" variable_1="i_K1"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current" component_1="plateau_potassium_current"/>
<map_variables variable_2="E_K1" variable_1="E_K1"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="plateau_potassium_current"/>
<map_variables variable_2="i_Kp" variable_1="i_Kp"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="sodium_background_current"/>
<map_variables variable_2="i_Na_b" variable_1="i_Na_b"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="sarcolemmal_calcium_pump"/>
<map_variables variable_2="i_p_Ca" variable_1="i_p_Ca"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="calcium_background_current"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
<map_variables variable_2="Cao" variable_1="Cao"/>
<map_variables variable_2="i_Ca_b" variable_1="i_Ca_b"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="sodium_potassium_pump"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<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_ionic_concentrations" component_1="non_specific_calcium_activated_current"/>
<map_variables variable_2="i_ns_Na" variable_1="i_ns_Na"/>
<map_variables variable_2="i_ns_K" variable_1="i_ns_K"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
<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="cell_geometry" component_1="intracellular_ionic_concentrations"/>
<map_variables variable_2="A_cap" variable_1="A_cap"/>
<map_variables variable_2="V_myo" variable_1="V_myo"/>
<map_variables variable_2="V_SS" variable_1="V_SS"/>
<map_variables variable_2="V_JSR" variable_1="V_JSR"/>
<map_variables variable_2="V_NSR" variable_1="V_NSR"/>
</connection>
<connection>
<map_components component_2="intracellular_Ca_fluxes" component_1="RyR_channel_states"/>
<map_variables variable_2="P_O1" variable_1="P_O1"/>
<map_variables variable_2="P_O2" variable_1="P_O2"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="RyR_channel_states"/>
<map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
</connection>
<connection>
<map_components component_2="SERCA_pump" component_1="intracellular_ionic_concentrations"/>
<map_variables variable_2="J_up" variable_1="J_up"/>
<map_variables variable_2="Ca_NSR" variable_1="Ca_NSR"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
</connection>
<connection>
<map_components component_2="intracellular_ionic_concentrations" component_1="intracellular_Ca_fluxes"/>
<map_variables variable_2="J_rel" variable_1="J_rel"/>
<map_variables variable_2="J_trpn" variable_1="J_trpn"/>
<map_variables variable_2="J_tr" variable_1="J_tr"/>
<map_variables variable_2="J_xfer" variable_1="J_xfer"/>
<map_variables variable_2="Ca_SS" variable_1="Ca_SS"/>
<map_variables variable_2="Ca_JSR" variable_1="Ca_JSR"/>
<map_variables variable_2="Ca_NSR" variable_1="Ca_NSR"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
</connection>
<connection>
<map_components component_2="troponin" component_1="intracellular_Ca_fluxes"/>
<map_variables variable_2="HTRPNCa" variable_1="HTRPNCa"/>
<map_variables variable_2="HTRPN_tot" variable_1="HTRPN_tot"/>
<map_variables variable_2="LTRPNCa" variable_1="LTRPNCa"/>
<map_variables variable_2="LTRPN_tot" variable_1="LTRPN_tot"/>
<map_variables variable_2="k_htrpn_minus" variable_1="k_htrpn_minus"/>
<map_variables variable_2="k_htrpn_plus" variable_1="k_htrpn_plus"/>
<map_variables variable_2="k_ltrpn_minus" variable_1="k_ltrpn_minus"/>
<map_variables variable_2="k_ltrpn_plus" variable_1="k_ltrpn_plus"/>
</connection>
<connection>
<map_components component_2="troponin" component_1="force_computation"/>
<map_variables variable_2="Force_norm" variable_1="Force_norm"/>
</connection>
<connection>
<map_components component_2="troponin" component_1="intracellular_ionic_concentrations"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
</connection>
<connection>
<map_components component_2="tropomyosin_cross_bridges" component_1="troponin"/>
<map_variables variable_2="k_ltrpn_minus" variable_1="k_ltrpn_minus"/>
<map_variables variable_2="k_ltrpn_plus" variable_1="k_ltrpn_plus"/>
<map_variables variable_2="LTRPNCa" variable_1="LTRPNCa"/>
<map_variables variable_2="LTRPN_tot" variable_1="LTRPN_tot"/>
</connection>
<connection>
<map_components component_2="tropomyosin_cross_bridges" component_1="force_computation"/>
<map_variables variable_2="f_01" variable_1="f_01"/>
<map_variables variable_2="f_12" variable_1="f_12"/>
<map_variables variable_2="f_23" variable_1="f_23"/>
<map_variables variable_2="g_XB" variable_1="g_XB"/>
<map_variables variable_2="SL" variable_1="SL"/>
<map_variables variable_2="P1" variable_1="P1"/>
<map_variables variable_2="P2" variable_1="P2"/>
<map_variables variable_2="N1" variable_1="N1"/>
<map_variables variable_2="P3" variable_1="P3"/>
</connection>
<!--
Various variables are passed between parent components and their
encapsulated gates.
-->
<connection>
<map_components component_2="fast_sodium_current_m_gate" component_1="fast_sodium_current"/>
<map_variables variable_2="m" variable_1="m"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current_h_gate" component_1="fast_sodium_current"/>
<map_variables variable_2="h" variable_1="h"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current_j_gate" component_1="fast_sodium_current"/>
<map_variables variable_2="j" variable_1="j"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_y_gate" component_1="L_type_Ca_channel"/>
<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="time_dependent_potassium_current_X_gate" component_1="time_dependent_potassium_current"/>
<map_variables variable_2="X" variable_1="X"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="time_dependent_potassium_current_Xi_gate" component_1="time_dependent_potassium_current"/>
<map_variables variable_2="Xi" variable_1="Xi"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current_K1_gate" component_1="time_independent_potassium_current"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="E_K1" variable_1="E_K1"/>
<map_variables variable_2="K1_infinity" variable_1="K1_infinity"/>
</connection>
<rdf:RDF>
<rdf:Bag rdf:about="rdf:#3327930a-1af9-4c82-a86d-70774454de3d">
<rdf:li>calcium dynamics</rdf:li>
<rdf:li>cardiac myocyte</rdf:li>
<rdf:li>Cardiac Myocyte</rdf:li>
<rdf:li>electrophysiology</rdf:li>
<rdf:li>myofilament mechanics</rdf:li>
<rdf:li>cardiac</rdf:li>
</rdf:Bag>
<rdf:Seq rdf:about="rdf:#citationAuthorsSeq">
<rdf:li rdf:resource="rdf:#author1Vcard"/>
<rdf:li rdf:resource="rdf:#author2Vcard"/>
<rdf:li rdf:resource="rdf:#author3Vcard"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#cf4e6901-4481-4416-b461-cc264e0a9d59">
<dcterms:W3CDTF>2002-06-24</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8dc79e2d-d733-46bf-9823-6cba8b2c2588">
<dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#cfa41065-f3e9-41d7-b2ae-a47cbc211110">
<dcterms:modified rdf:resource="rdf:#735565b1-a521-4d7c-b198-09bd6d8f0b7a"/>
<rdf:value>
Added more metadata.
</rdf:value>
<cmeta:modifier rdf:resource="rdf:#a165ac20-3a7c-4713-8edf-bcb47e57a145"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7d005e9c-1569-4f07-b235-bc637c2c3012">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value rdf:resource="rdf:#3327930a-1af9-4c82-a86d-70774454de3d"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6de2c6e3-a9e5-4924-9298-8704b91b2c4e">
<dc:title>American Journal of Physiology: Heart and Circulatory Physiology</dc:title>
</rdf:Description>
<rdf:Description rdf:about="#rice_jafri_winslow_2000_version02">
<dc:title>
The Rice et al's 2000 detailed model of ventricular cardiac
myocytes
</dc:title>
<cmeta:bio_entity>Cardiac Myocyte</cmeta:bio_entity>
<cmeta:comment rdf:resource="rdf:#fdac2647-6378-4d56-96e0-1b7906026c31"/>
<bqs:reference rdf:resource="rdf:#b3a3df16-9fbf-4039-abc9-e2b61a621801"/>
<bqs:reference rdf:resource="rdf:#dd6972d1-9747-47f0-9bb1-26d908bbb8e2"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author3VcardN">
<vCard:Given>Raimond</vCard:Given>
<vCard:Family>Winslow</vCard:Family>
<vCard:Other>L</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>
The University of Auckland, Bioengineering Institute
</dc:publisher>
<cmeta:modification rdf:resource="rdf:#24eb0b3f-f381-4ccf-ac5d-70644a657622"/>
<cmeta:modification rdf:resource="rdf:#cfa41065-f3e9-41d7-b2ae-a47cbc211110"/>
<dcterms:created rdf:resource="rdf:#cf4e6901-4481-4416-b461-cc264e0a9d59"/>
<dc:creator rdf:resource="rdf:#2f43be39-8e3c-4605-9fde-02dbe2610013"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4ca4b792-d074-4e5c-b3f9-041f9bd45d57">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2f43be39-8e3c-4605-9fde-02dbe2610013">
<vCard:ORG rdf:resource="rdf:#14d2d51a-cf53-434a-b830-54a20a4b8ac3"/>
<vCard:EMAIL rdf:resource="rdf:#9260fb1a-dd92-4fb1-b17e-d6545f37316d"/>
<vCard:N rdf:resource="rdf:#9b89d4f9-1535-466b-8e1e-45d105a5df0e"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author2Vcard">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#author2VcardN"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#dd6972d1-9747-47f0-9bb1-26d908bbb8e2">
<dc:subject rdf:resource="rdf:#7d005e9c-1569-4f07-b235-bc637c2c3012"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author2VcardN">
<vCard:Given>M</vCard:Given>
<vCard:Family>Jafri</vCard:Family>
<vCard:Other>Saleet</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#14d2d51a-cf53-434a-b830-54a20a4b8ac3">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1be093f0-0715-463d-8a1b-938f856c7139">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#045b84e5-d335-4460-a7e8-db2f57cffc3e">
<dcterms:W3CDTF>2000-03-01</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#74b49ad5-99d9-486c-a023-ed1073ab0991">
<vCard:N rdf:resource="rdf:#59c07677-3e2a-4a4d-ada1-48f44959a4f4"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a165ac20-3a7c-4713-8edf-bcb47e57a145">
<vCard:N rdf:resource="rdf:#4ca4b792-d074-4e5c-b3f9-041f9bd45d57"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fb91e637-d3e7-42cc-979f-b922d3b981c0">
<dc:creator rdf:resource="rdf:#citationAuthorsSeq"/>
<dc:title>Modeling short-term interval-force relations in cardiac muscle</dc:title>
<bqs:volume>278</bqs:volume>
<bqs:first_page>H913</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#6de2c6e3-a9e5-4924-9298-8704b91b2c4e"/>
<dcterms:issued rdf:resource="rdf:#045b84e5-d335-4460-a7e8-db2f57cffc3e"/>
<bqs:last_page>H931</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#59c07677-3e2a-4a4d-ada1-48f44959a4f4">
<vCard:Given>Autumn</vCard:Given>
<vCard:Family>Cuellar</vCard:Family>
<vCard:Other>A</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b3a3df16-9fbf-4039-abc9-e2b61a621801">
<bqs:Pubmed_id>10710361</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#fb91e637-d3e7-42cc-979f-b922d3b981c0"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fdac2647-6378-4d56-96e0-1b7906026c31">
<dc:creator rdf:resource="rdf:#1be093f0-0715-463d-8a1b-938f856c7139"/>
<rdf:value>This is the CellML description of Rice et al's 2000 detailed model of ventricular cardiac myocytes. It simulates action potentials, Ca2+-handling mechanisms and isometric force generation by the myofilaments.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author1Vcard">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#author1VcardN"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#author1VcardN">
<vCard:Given>John</vCard:Given>
<vCard:Family>Rice</vCard:Family>
<vCard:Other>Jeremy</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9260fb1a-dd92-4fb1-b17e-d6545f37316d">
<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:#author3Vcard">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#author3VcardN"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9b89d4f9-1535-466b-8e1e-45d105a5df0e">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#735565b1-a521-4d7c-b198-09bd6d8f0b7a">
<dcterms:W3CDTF>2002-07-18</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#24eb0b3f-f381-4ccf-ac5d-70644a657622">
<dcterms:modified rdf:resource="rdf:#8dc79e2d-d733-46bf-9823-6cba8b2c2588"/>
<rdf:value>
Added publication date information.
</rdf:value>
<cmeta:modifier rdf:resource="rdf:#74b49ad5-99d9-486c-a023-ed1073ab0991"/>
</rdf:Description>
</rdf:RDF>
</model>
