- Author:
- pmr2.import <nobody@models.cellml.org>
- Date:
- 2008-05-07 15:22:00+12:00
- Desc:
- committing version06 of beeler_reuter_1977
- Permanent Source URI:
- https://models.cellml.org/workspace/beeler_reuter_1977/rawfile/575ef44b878c3de8946da61b270ba65564aad462/beeler_reuter_1977.cellml
<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 07/05/2008 at 16:14:45 using:
COR (0.9.31.927)
Copyright 2002-2008 Dr Alan Garny
http://COR.physiol.ox.ac.uk/ - COR@physiol.ox.ac.uk
CellML 1.0 was used to generate this cellular model
http://www.CellML.org/
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="beeler_reuter_1977_version06" name="beeler_reuter_1977_version06">
<!--
Below, we define some additional units for association with variables and
constants within the model. The identifiers are fairly self-explanatory.
-->
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Beeler-Reuter Mammalian Ventricular Model 1977</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 model has been curated and is known to run and reproduce the published results in PCEnv and COR. A PCEnv session file is also associated with this model.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
In contrast to the earlier Purkinje fibre ionic current models of <ulink url="${HTML_EXMPL_N_MODEL}">D. Noble (1962)</ulink> and <ulink url="${HTML_EXMPL_MNT_MODEL}">R.E. McAllister, D. Noble and R.W. Tsien (1975)</ulink>, the G.W. Beeler and H. Reuter 1977 model was developed to describe the mammalian ventricular action potential. Not all the ionic currents of the Purkinje fibre model are present in ventricular tissue; therefore, this model is simpler than the MNT model. The total ionic flux is divided into only four discrete, individual ionic currents (see <xref linkend="fig_cell_diagram"/> below). The main additional feature of the Beeler-Reuter ionic current model is a representation of the intracellular calcium ion concentration.
</para>
<para>
The complete original paper reference is cited below:
</para>
<para>
Reconstruction of the action potential of ventricular myocardial fibres, Beeler, G.W. and Reuter, H. 1977 <ulink url="http://www.jphysiol.org/">
<emphasis>Journal of Physiology</emphasis>
</ulink>, 268, 177-210. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=874889&dopt=Abstract">PubMed ID: 874889</ulink>
</para>
<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>cell diagram of the Beeler-Reuter model showing ionic currents across the cell surface membrane</title>
</objectinfo>
<imagedata fileref="beeler_reuter_1977.png"/>
</imageobject>
</mediaobject>
<caption>A schematic diagram describing the current flows across the cell membrane that are captured in the BR model.</caption>
</informalfigure>
<informalfigure float="0" id="fig_cellml_rendering">
<mediaobject>
<imageobject>
<objectinfo>
<title>the cellml rendering of the Beeler-Reuter model</title>
</objectinfo>
<imagedata fileref="cellml_rendering.gif"/>
</imageobject>
</mediaobject>
<caption>The network defined in the CellML description of the Beeler-Reuter model. A key describing the significance of the shapes of the components and the colours of the connections between them is in the <ulink url="http://www.cellml.org/tutorial/notation/?searchterm=notation%20guide">notation guide</ulink>. For simplicity, not all the variables are shown.</caption>
</informalfigure>
<para>
The membrane physically contains the currents as indicated by the blue arrows in <xref linkend="fig_cellml_rendering"/>. The currents act independently and are not connected to each other. Several of the channels encapsulate <emphasis>and</emphasis> contain further components which represent activation and inactivation gates. The addition of an encapsulation relationship informs modellers and processing software that the gates are important parts of the current model. It also prevents any other components that aren't also encapsulated by the parent component from connecting to its gates, effectively hiding them from the rest of the model.
</para>
<para>
The breakdown of the model into components and the definition of encapsulation and containment relationships between them is somewhat arbitrary. When considering how a model should be broken into components, modellers are encouraged to consider which parts of a model might be re-used and how the physiological elements of the system being modelled are naturally bounded. Containment relationships should be used to provide simple rendering information for processing software (ideally, this will correspond to the layout of the physical system), and encapsulation should be used to group sets of components into sub-models.
</para>
</sect1>
</article>
</documentation>
<units name="ms">
<unit units="second" prefix="milli"/>
</units>
<units name="per_ms">
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<units name="mV">
<unit units="volt" prefix="milli"/>
</units>
<units name="per_mV">
<unit units="volt" prefix="milli" exponent="-1"/>
</units>
<units name="per_mV_ms">
<unit units="mV" exponent="-1"/>
<unit units="ms" exponent="-1"/>
</units>
<units name="mS_per_mm2">
<unit units="siemens" prefix="milli"/>
<unit units="metre" prefix="milli" exponent="-2"/>
</units>
<units name="uF_per_mm2">
<unit units="farad" prefix="micro"/>
<unit units="metre" prefix="milli" exponent="-2"/>
</units>
<units name="uA_per_mm2">
<unit units="ampere" prefix="micro"/>
<unit units="metre" prefix="milli" exponent="-2"/>
</units>
<units name="concentration_units">
<unit units="mole" prefix="nano"/>
<unit units="metre" prefix="milli" exponent="-3"/>
</units>
<units name="per_concentration_units">
<unit units="concentration_units" exponent="-1"/>
</units>
<units name="coulomb_per_mole">
<unit units="coulomb"/>
<unit units="mole" exponent="-1"/>
</units>
<units name="per_mm">
<unit units="metre" prefix="-3" exponent="-1"/>
</units>
<component name="environment">
<variable name="time" units="ms" public_interface="out"/>
</component>
<component name="membrane">
<variable name="V" units="mV" initial_value="-84.624" public_interface="out"/>
<variable name="C" units="uF_per_mm2" initial_value="0.01"/>
<variable name="time" units="ms" public_interface="in"/>
<variable name="i_Na" units="uA_per_mm2" public_interface="in"/>
<variable name="i_s" units="uA_per_mm2" public_interface="in"/>
<variable name="i_x1" units="uA_per_mm2" public_interface="in"/>
<variable name="i_K1" units="uA_per_mm2" public_interface="in"/>
<variable name="Istim" units="uA_per_mm2" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>V</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>Istim</ci>
<apply>
<plus/>
<ci>i_Na</ci>
<ci>i_s</ci>
<ci>i_x1</ci>
<ci>i_K1</ci>
</apply>
</apply>
<ci>C</ci>
</apply>
</apply>
</math>
</component>
<component name="sodium_current">
<variable name="i_Na" units="uA_per_mm2" public_interface="out"/>
<variable name="g_Na" units="mS_per_mm2" initial_value="4e-2"/>
<variable name="E_Na" units="mV" initial_value="50"/>
<variable name="g_Nac" units="mS_per_mm2" initial_value="3e-5"/>
<variable name="time" units="ms" public_interface="in" private_interface="out"/>
<variable name="V" units="mV" public_interface="in" private_interface="out"/>
<variable name="m" units="dimensionless" private_interface="in"/>
<variable name="h" units="dimensionless" private_interface="in"/>
<variable name="j" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_Na</ci>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<ci>g_Na</ci>
<apply>
<power/>
<ci>m</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>h</ci>
<ci>j</ci>
</apply>
<ci>g_Nac</ci>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_m_gate">
<variable name="m" units="dimensionless" initial_value="0.011" public_interface="out"/>
<variable name="alpha_m" units="per_ms"/>
<variable name="beta_m" units="per_ms"/>
<variable name="V" units="mV" public_interface="in"/>
<variable name="time" units="ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_m</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mV_ms">1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">47</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mV">0.1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">47</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_m</ci>
<apply>
<times/>
<cn cellml:units="per_ms">40</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mV">0.056</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">72</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>m</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_m</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>m</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_m</ci>
<ci>m</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_h_gate">
<variable name="h" units="dimensionless" initial_value="0.988" public_interface="out"/>
<variable name="alpha_h" units="per_ms"/>
<variable name="beta_h" units="per_ms"/>
<variable name="V" units="mV" public_interface="in"/>
<variable name="time" units="ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_h</ci>
<apply>
<times/>
<cn cellml:units="per_ms">0.126</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mV">0.25</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">77</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_h</ci>
<apply>
<divide/>
<cn cellml:units="per_ms">1.7</cn>
<apply>
<plus/>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mV">0.082</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">22.5</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>h</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_h</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>h</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_h</ci>
<ci>h</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_current_j_gate">
<variable name="j" units="dimensionless" initial_value="0.975" public_interface="out"/>
<variable name="alpha_j" units="per_ms"/>
<variable name="beta_j" units="per_ms"/>
<variable name="V" units="mV" public_interface="in"/>
<variable name="time" units="ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_j</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_ms">0.055</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mV">0.25</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">78</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mV">0.2</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">78</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_j</ci>
<apply>
<divide/>
<cn cellml:units="per_ms">0.3</cn>
<apply>
<plus/>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mV">0.1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">32</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<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</cn>
<ci>j</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_j</ci>
<ci>j</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_inward_current">
<variable name="i_s" units="uA_per_mm2" public_interface="out"/>
<variable name="g_s" units="mS_per_mm2" initial_value="9e-4"/>
<variable name="E_s" units="mV"/>
<variable name="Cai" units="concentration_units" initial_value="1e-4"/>
<variable name="time" units="ms" public_interface="in" private_interface="out"/>
<variable name="V" units="mV" public_interface="in" private_interface="out"/>
<variable name="d" units="dimensionless" private_interface="in"/>
<variable name="f" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E_s</ci>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="mV">82.3</cn>
</apply>
<apply>
<times/>
<cn cellml:units="mV">13.0287</cn>
<apply>
<ln/>
<apply>
<times/>
<ci>Cai</ci>
<cn cellml:units="per_concentration_units">0.001</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>i_s</ci>
<apply>
<times/>
<ci>g_s</ci>
<ci>d</ci>
<ci>f</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_s</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Cai</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mm">0.01</cn>
</apply>
<ci>i_s</ci>
</apply>
<cn cellml:units="coulomb_per_mole">1</cn>
</apply>
<apply>
<times/>
<cn cellml:units="per_ms">0.07</cn>
<apply>
<minus/>
<cn cellml:units="concentration_units">0.0001</cn>
<ci>Cai</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_inward_current_d_gate">
<variable name="d" units="dimensionless" initial_value="0.003" public_interface="out"/>
<variable name="alpha_d" units="per_ms"/>
<variable name="beta_d" units="per_ms"/>
<variable name="V" units="mV" public_interface="in"/>
<variable name="time" units="ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_d</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_ms">0.095</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="mV">5</cn>
</apply>
</apply>
<cn cellml:units="mV">100</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<cn cellml:units="mV">5</cn>
</apply>
</apply>
<cn cellml:units="mV">13.89</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_d</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_ms">0.07</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">44</cn>
</apply>
</apply>
<cn cellml:units="mV">59</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">44</cn>
</apply>
<cn cellml:units="mV">20</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>d</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_d</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>d</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_d</ci>
<ci>d</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_inward_current_f_gate">
<variable name="f" units="dimensionless" initial_value="0.994" public_interface="out"/>
<variable name="alpha_f" units="per_ms"/>
<variable name="beta_f" units="per_ms"/>
<variable name="V" units="mV" public_interface="in"/>
<variable name="time" units="ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_f</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_ms">0.012</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">28</cn>
</apply>
</apply>
<cn cellml:units="mV">125</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">28</cn>
</apply>
<cn cellml:units="mV">6.67</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_f</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_ms">0.0065</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">30</cn>
</apply>
</apply>
<cn cellml:units="mV">50</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">30</cn>
</apply>
</apply>
<cn cellml:units="mV">5</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>f</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_f</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>f</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_f</ci>
<ci>f</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_dependent_outward_current">
<variable name="i_x1" units="uA_per_mm2" public_interface="out"/>
<variable name="time" units="ms" public_interface="in" private_interface="out"/>
<variable name="V" units="mV" public_interface="in" private_interface="out"/>
<variable name="x1" units="dimensionless" private_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_x1</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x1</ci>
<cn cellml:units="uA_per_mm2">8e-3</cn>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_mV">0.04</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">77</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_mV">0.04</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">35</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_dependent_outward_current_x1_gate">
<variable name="x1" units="dimensionless" initial_value="0.0001" public_interface="out"/>
<variable name="alpha_x1" units="per_ms"/>
<variable name="beta_x1" units="per_ms"/>
<variable name="V" units="mV" public_interface="in"/>
<variable name="time" units="ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_x1</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_ms">5e-4</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">50</cn>
</apply>
<cn cellml:units="mV">12.1</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">50</cn>
</apply>
<cn cellml:units="mV">17.5</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_x1</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_ms">0.0013</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">20</cn>
</apply>
</apply>
<cn cellml:units="mV">16.67</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">20</cn>
</apply>
</apply>
<cn cellml:units="mV">25</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>x1</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_x1</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>x1</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_x1</ci>
<ci>x1</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_independent_outward_current">
<variable name="i_K1" units="uA_per_mm2" public_interface="out"/>
<variable name="time" units="ms" public_interface="in"/>
<variable name="V" units="mV" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>i_K1</ci>
<apply>
<times/>
<cn cellml:units="uA_per_mm2">0.0035</cn>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">4</cn>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_mV">0.04</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">85</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<plus/>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_mV">0.08</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">53</cn>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_mV">0.04</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">53</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_mV">0.2</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">23</cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_mV">0.04</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="mV">23</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="stimulus_protocol">
<variable name="Istim" units="uA_per_mm2" public_interface="out"/>
<variable name="IstimStart" units="ms" initial_value="10"/>
<variable name="IstimEnd" units="ms" initial_value="50000"/>
<variable name="IstimAmplitude" units="uA_per_mm2" initial_value="0.5"/>
<variable name="IstimPeriod" units="ms" initial_value="1000"/>
<variable name="IstimPulseDuration" units="ms" initial_value="1"/>
<variable name="time" units="ms" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Istim</ci>
<piecewise>
<piece>
<ci>IstimAmplitude</ci>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<ci>IstimStart</ci>
</apply>
<apply>
<leq/>
<ci>time</ci>
<ci>IstimEnd</ci>
</apply>
<apply>
<leq/>
<apply>
<minus/>
<apply>
<minus/>
<ci>time</ci>
<ci>IstimStart</ci>
</apply>
<apply>
<times/>
<apply>
<floor/>
<apply>
<divide/>
<apply>
<minus/>
<ci>time</ci>
<ci>IstimStart</ci>
</apply>
<ci>IstimPeriod</ci>
</apply>
</apply>
<ci>IstimPeriod</ci>
</apply>
</apply>
<ci>IstimPulseDuration</ci>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="uA_per_mm2">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="sodium_current">
<component_ref component="sodium_current_m_gate"/>
<component_ref component="sodium_current_h_gate"/>
<component_ref component="sodium_current_j_gate"/>
</component_ref>
<component_ref component="slow_inward_current">
<component_ref component="slow_inward_current_d_gate"/>
<component_ref component="slow_inward_current_f_gate"/>
</component_ref>
<component_ref component="time_dependent_outward_current">
<component_ref component="time_dependent_outward_current_x1_gate"/>
</component_ref>
<component_ref component="time_independent_outward_current"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="sodium_current">
<component_ref component="sodium_current_m_gate"/>
<component_ref component="sodium_current_h_gate"/>
<component_ref component="sodium_current_j_gate"/>
</component_ref>
<component_ref component="slow_inward_current">
<component_ref component="slow_inward_current_d_gate"/>
<component_ref component="slow_inward_current_f_gate"/>
</component_ref>
<component_ref component="time_dependent_outward_current">
<component_ref component="time_dependent_outward_current_x1_gate"/>
</component_ref>
</group>
<connection>
<map_components component_1="membrane" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="slow_inward_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="time_dependent_outward_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="time_independent_outward_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="sodium_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_Na" variable_2="i_Na"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="slow_inward_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_s" variable_2="i_s"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="time_dependent_outward_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_x1" variable_2="i_x1"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="time_independent_outward_current"/>
<map_variables variable_1="V" variable_2="V"/>
<map_variables variable_1="i_K1" variable_2="i_K1"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="sodium_current_m_gate"/>
<map_variables variable_1="m" variable_2="m"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="sodium_current_h_gate"/>
<map_variables variable_1="h" variable_2="h"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="sodium_current" component_2="sodium_current_j_gate"/>
<map_variables variable_1="j" variable_2="j"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="time_dependent_outward_current" component_2="time_dependent_outward_current_x1_gate"/>
<map_variables variable_1="x1" variable_2="x1"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="slow_inward_current" component_2="slow_inward_current_d_gate"/>
<map_variables variable_1="d" variable_2="d"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="slow_inward_current" component_2="slow_inward_current_f_gate"/>
<map_variables variable_1="f" variable_2="f"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection>
<map_components component_1="membrane" component_2="stimulus_protocol"/>
<map_variables variable_1="Istim" variable_2="Istim"/>
</connection>
<connection>
<map_components component_1="environment" component_2="stimulus_protocol"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<rdf:RDF>
<rdf:Seq rdf:about="rdf:#8424fc29-10dd-4d52-9d8b-5ed3584e2708">
<rdf:li rdf:resource="rdf:#85ef4ed4-88b8-482b-80da-a37edce16c95"/>
<rdf:li rdf:resource="rdf:#34cf9f7a-d252-464a-bf07-f9d0e989acfe"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#302ce03c-5af0-45f0-89ba-d5518cb6d91c">
<dc:creator rdf:resource="rdf:#8424fc29-10dd-4d52-9d8b-5ed3584e2708"/>
<dc:title>Reconstruction of the action potential of ventricular myocardial fibres</dc:title>
<bqs:volume>268</bqs:volume>
<bqs:first_page>177</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#ec81b107-def6-4277-87e8-0a89f837d856"/>
<dcterms:issued rdf:resource="rdf:#38a3a9b5-d63a-4736-8cbb-6c6e1edd1027"/>
<bqs:last_page>210</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="">
<cmeta:comment rdf:resource="rdf:#cdb74574-9908-4a67-9c3b-f9d5bc24039c"/>
<dcterms:created rdf:resource="rdf:#e8ff124c-9968-4230-afd2-da5e24ce4e74"/>
<dc:creator rdf:resource="rdf:#0df426b2-7ef4-40f2-baf4-4d05eac8b3b5"/>
<cmeta:modification rdf:resource="rdf:#77f98dda-c979-49f1-9bcf-da4944906f91"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ef9ea29b-ae5e-4db0-975e-42fe8f0092ed">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#77f98dda-c979-49f1-9bcf-da4944906f91">
<dcterms:modified rdf:resource="rdf:#e7527bce-f615-4263-afc0-c075d70b6e23"/>
<rdf:value>Added an intial value for X1 to enable the model to run.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#baef8605-6d17-47bc-b214-261f23273934"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#13b8a84a-b5f2-4bae-ae2d-2b1153145132">
<vCard:Given>H</vCard:Given>
<vCard:Family>Reuter</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e8ff124c-9968-4230-afd2-da5e24ce4e74">
<dcterms:W3CDTF>2008-05-08T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0df426b2-7ef4-40f2-baf4-4d05eac8b3b5">
<vCard:ORG rdf:resource="rdf:#e289c80e-007b-4ac0-b716-546e0926483d"/>
<vCard:EMAIL rdf:resource="rdf:#806732fe-3bdb-4bac-92e6-8286ab04a46b"/>
<vCard:N rdf:resource="rdf:#ee66d65b-4cbe-4af4-810a-a99bbfb17ad1"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2324bd00-6aaf-4cf6-a159-2eb626f467a7">
<dc:creator rdf:resource="rdf:#a31de35b-050f-40be-a915-76b474997067"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3277b835-42e2-42b1-9729-b5884be39914">
<vCard:Given>G</vCard:Given>
<vCard:Family>Beeler</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ec81b107-def6-4277-87e8-0a89f837d856">
<dc:title>Journal of Physiology</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e7527bce-f615-4263-afc0-c075d70b6e23">
<dcterms:W3CDTF>2008-05-08T03:15:26+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#85ef4ed4-88b8-482b-80da-a37edce16c95">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#3277b835-42e2-42b1-9729-b5884be39914"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#34cf9f7a-d252-464a-bf07-f9d0e989acfe">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#13b8a84a-b5f2-4bae-ae2d-2b1153145132"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e289c80e-007b-4ac0-b716-546e0926483d">
<vCard:Orgname>Bioengineering Institute, University of Auckland</vCard:Orgname>
</rdf:Description>
<rdf:Description rdf:about="rdf:#38a3a9b5-d63a-4736-8cbb-6c6e1edd1027">
<dcterms:W3CDTF>1977-06-00 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#806732fe-3bdb-4bac-92e6-8286ab04a46b">
<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="#beeler_reuter_1977_version06">
<bqs:reference rdf:resource="rdf:#9585afad-8063-410a-843f-15a21ac4da7d"/>
<cmeta:comment rdf:resource="rdf:#2324bd00-6aaf-4cf6-a159-2eb626f467a7"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#baef8605-6d17-47bc-b214-261f23273934">
<vCard:N rdf:resource="rdf:#ef9ea29b-ae5e-4db0-975e-42fe8f0092ed"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#cdb74574-9908-4a67-9c3b-f9d5bc24039c">
<dc:creator rdf:resource="rdf:#0212d7cf-9ee6-41a5-b633-aaddb12a73cd"/>
<rdf:value>In contrast to the earlier Purkinje fibre ionic current models of <ulink url="${HTML_EXMPL_N_MODEL}">D. Noble (1962)</ulink> and <ulink url="${HTML_EXMPL_MNT_MODEL}">R.E. McAllister, D. Noble and R.W. Tsien (1975)</ulink>, the G.W. Beeler and H. Reuter 1977 model was developed to describe the mammalian ventricular action potential. Not all the ionic currents of the Purkinje fibre model are present in ventricular tissue; therefore, this model is simpler than the MNT model. The total ionic flux is divided into only four discrete, individual ionic currents (see <xref linkend="fig_cell_diagram"/> below). The main additional feature of the Beeler-Reuter ionic current model is a representation of the intracellular calcium ion concentration.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ee66d65b-4cbe-4af4-810a-a99bbfb17ad1">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0212d7cf-9ee6-41a5-b633-aaddb12a73cd">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9585afad-8063-410a-843f-15a21ac4da7d">
<bqs:Pubmed_id>874889</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#302ce03c-5af0-45f0-89ba-d5518cb6d91c"/>
</rdf:Description>
</rdf:RDF>
</model>
