- Author:
- pmr2.import <nobody@models.cellml.org>
- Date:
- 2009-05-28 04:10:56+12:00
- Desc:
- committing version02 of bueno_2007
- Permanent Source URI:
- http://models.cellml.org/workspace/bueno_2007/rawfile/b1ca358bc87d1dfd65ec0500fba923913bc615dd/bueno_2007.cellml
<?xml version='1.0' encoding='utf-8'?>
<!--
This CellML file was generated on 25/02/2008 at 13:56:30 using:
COR (0.9.31.873)
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#" cmeta:id="bueno_2007" name="bueno_2007">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This CellML model is known to run in both COR and PCEnv to recreate the results of the original model.
</para>
<para>
ValidateCellML verifies this model as valid CellML with fully consistent units.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
Modeling the dynamics of electrical wave propagation in cardiac tissue requires the use of appropriate ionic models capable to reproduce the key characteristics of this type of systems, suitable fitted for the specie under study. In the current study Alfonso Bueno presents a new minimal model for human ventricular action potentials which is able to accurately reproduce a range of experimentally measured characteristics of cardiac tissue, including threshold for excitation, action potential morphology, rate of rise, and rate dependence properties (action potential duration and conduction velocity) for epicardial, endocardial and midmyocardial cells. Since the number of parameters involved is significantly smaller than in other models, it is easier to understand the role that each parameter plays in each of the above mentioned characteristics, therefore making substantially easier to obtain fittings for different data. Besides, in terms of computational efficiency the new model can be simulated several orders of magnitude faster than other human ventricular models, which makes the new model especially suitable for large scale and whole organ simulations.
</para>
<para>
The applicability and usefulness of the new model have been proved by addressing the study of the Brugada syndrome, an important hereditary genetic disease able to cause sudden cardiac death in human patients even without identifiable structural abnormalities. Despite not incorporating independent descriptions for each of the known currents measured through the cell membrane, the model is however capable to correctly reproduce both the healthy and pathological conditions of the epicardial tissue layer, the latter being characterized by the presence of a much deeper notch in the cardiac action potential as a result of a loss of functionality of sodium channels which leaves unopposed the transient outward potassium current.
</para>
<para>
Despite being able of reproducing a great amount of available experimental data, one important limitation of the proposed minimal model is that it lacks of a complete description of intracellular calcium dynamics. Hence, the model cannot be used to study particular conditions of importance in the generation and development of arrhythmias, such as calcium overload, spontaneous calcium release or calcium-induced alternans. This limitation could be overcome however by coupling our simplified model of the cardiac transmembrane potential to a reduced model of the intracellular calcium cycling of the ventricular myocyte, such as the one suggested by Shiferaw et al. (2003), therefore yielding the first implementation of a second-generation simplified ionic model of the cardiac myocyte.
</para>
<!--
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="../images/bueno_model_2007/bueno_2007.png"/>
</imageobject>
</mediaobject>
<caption>A schematic diagram of ...</caption>
</informalfigure>
-->
<para>
This model was taken from the PhD thesis of Alfonso Bueno:
</para>
<para>
Mathematical modeling and spectral simulation of genetic diseases in the human heart, Alfonso Bueno, Department of Mathematics, University of Castilla-La Mancha, Ciudad Real, May 2007.
</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_cm2">
<unit units="siemens" prefix="milli"/>
<unit units="metre" prefix="centi" exponent="-2"/>
</units>
<units name="uF_per_cm2">
<unit units="farad" prefix="micro"/>
<unit units="metre" prefix="centi" exponent="-2"/>
</units>
<units name="uA_per_cm2">
<unit units="ampere" prefix="micro"/>
<unit units="metre" prefix="centi" exponent="-2"/>
</units>
<units name="concentration_units">
<unit units="mole" prefix="nano"/>
<unit units="metre" prefix="milli" exponent="-3"/>
</units>
<component name="environment">
<variable units="ms" public_interface="out" name="time"/>
<variable units="dimensionless" public_interface="out" name="epi" initial_value="1"/>
<variable units="dimensionless" public_interface="out" name="endo" initial_value="0"/>
<variable units="dimensionless" public_interface="out" name="mcell" initial_value="0"/>
</component>
<component name="membrane">
<variable units="dimensionless" public_interface="out" name="u" initial_value="0"/>
<variable units="mV" name="Vm"/>
<variable units="mV" name="V_0" initial_value="-83"/>
<variable units="mV" name="V_fi" initial_value="2.7"/>
<variable units="ms" public_interface="in" name="time"/>
<variable units="per_ms" public_interface="in" name="J_fi"/>
<variable units="per_ms" public_interface="in" name="J_so"/>
<variable units="per_ms" public_interface="in" name="J_si"/>
<variable units="per_ms" name="J_stim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_stim</ci>
<piecewise>
<piece>
<apply>
<minus/>
<cn cellml:units="per_ms">1</cn>
</apply>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<cn cellml:units="ms">100</cn>
</apply>
<apply>
<leq/>
<ci>time</ci>
<cn cellml:units="ms">101</cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="per_ms">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>u</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>J_fi</ci>
<ci>J_so</ci>
<ci>J_si</ci>
<ci>J_stim</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Vm</ci>
<apply>
<plus/>
<ci>V_0</ci>
<apply>
<times/>
<ci>u</ci>
<apply>
<minus/>
<ci>V_fi</ci>
<ci>V_0</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="m">
<variable units="dimensionless" public_interface="out" name="m"/>
<variable units="dimensionless" public_interface="out" name="u_m" initial_value="0.3"/>
<variable units="dimensionless" public_interface="in" name="u"/>
<variable units="ms" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>m</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0</cn>
<apply>
<lt/>
<ci>u</ci>
<ci>u_m</ci>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="p">
<variable units="dimensionless" public_interface="out" name="p"/>
<variable units="dimensionless" public_interface="out" name="u_p" initial_value="0.13"/>
<variable units="dimensionless" public_interface="in" name="u"/>
<variable units="ms" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>p</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0</cn>
<apply>
<lt/>
<ci>u</ci>
<ci>u_p</ci>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="q">
<variable units="dimensionless" public_interface="out" name="q"/>
<variable units="dimensionless" public_interface="in" name="epi"/>
<variable units="dimensionless" public_interface="in" name="endo"/>
<variable units="dimensionless" public_interface="in" name="mcell"/>
<variable units="dimensionless" public_interface="out" name="u_q"/>
<variable units="dimensionless" public_interface="in" name="u"/>
<variable units="ms" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>u_q</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0.006</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">0.024</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0.1</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>q</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0</cn>
<apply>
<lt/>
<ci>u</ci>
<ci>u_q</ci>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="r">
<variable units="dimensionless" public_interface="out" name="r"/>
<variable units="dimensionless" public_interface="in" name="epi"/>
<variable units="dimensionless" public_interface="in" name="endo"/>
<variable units="dimensionless" public_interface="in" name="mcell"/>
<variable units="dimensionless" name="u_r"/>
<variable units="dimensionless" public_interface="in" name="u"/>
<variable units="ms" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>u_r</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0.006</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">0.006</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0.005</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>r</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0</cn>
<apply>
<lt/>
<ci>u</ci>
<ci>u_r</ci>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1</cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="fast_inward_current">
<variable units="per_ms" public_interface="out" name="J_fi"/>
<variable units="ms" name="tau_fi"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="epi"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="endo"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="mcell"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="m"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="q"/>
<variable units="dimensionless" public_interface="in" name="u_p"/>
<variable units="dimensionless" name="u_u"/>
<variable units="dimensionless" public_interface="in" name="u_m"/>
<variable units="ms" public_interface="in" private_interface="out" name="time"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="u"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="u_q"/>
<variable units="dimensionless" private_interface="in" name="v"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>u_u</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">1.55</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">1.56</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">1.61</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_fi</ci>
<piecewise>
<piece>
<cn cellml:units="ms">0.11</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">0.104</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">0.078</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>J_fi</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>m</ci>
</apply>
<ci>v</ci>
<apply>
<minus/>
<ci>u</ci>
<ci>u_m</ci>
</apply>
<apply>
<minus/>
<ci>u_u</ci>
<ci>u</ci>
</apply>
</apply>
<ci>tau_fi</ci>
</apply>
</apply>
</math>
</component>
<component name="fast_inward_current_v_gate">
<variable units="dimensionless" public_interface="out" name="v" initial_value="1"/>
<variable units="dimensionless" name="v_inf"/>
<variable units="dimensionless" public_interface="in" name="epi"/>
<variable units="dimensionless" public_interface="in" name="endo"/>
<variable units="dimensionless" public_interface="in" name="mcell"/>
<variable units="ms" name="tau_v_minus"/>
<variable units="ms" name="tau_v1_minus"/>
<variable units="ms" name="tau_v2_minus"/>
<variable units="ms" name="tau_v_plus" initial_value="1.45"/>
<variable units="dimensionless" public_interface="in" name="q"/>
<variable units="dimensionless" public_interface="in" name="m"/>
<variable units="dimensionless" public_interface="in" name="u"/>
<variable units="dimensionless" public_interface="in" name="u_q"/>
<variable units="ms" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>tau_v2_minus</ci>
<piecewise>
<piece>
<cn cellml:units="ms">1150</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">10</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">1.45</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_v1_minus</ci>
<piecewise>
<piece>
<cn cellml:units="ms">60</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">75</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">80</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>v_inf</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">1</cn>
<apply>
<lt/>
<ci>u</ci>
<ci>u_q</ci>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_v_minus</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>q</ci>
<ci>tau_v2_minus</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>q</ci>
</apply>
<ci>tau_v1_minus</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>m</ci>
</apply>
<apply>
<minus/>
<ci>v_inf</ci>
<ci>v</ci>
</apply>
</apply>
<ci>tau_v_minus</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>m</ci>
<ci>v</ci>
</apply>
<ci>tau_v_plus</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_outward_current">
<variable units="per_ms" public_interface="out" name="J_so"/>
<variable units="dimensionless" public_interface="in" name="epi"/>
<variable units="dimensionless" public_interface="in" name="endo"/>
<variable units="dimensionless" public_interface="in" name="mcell"/>
<variable units="ms" name="tau_o"/>
<variable units="ms" name="tau_o1"/>
<variable units="ms" name="tau_o2"/>
<variable units="ms" name="tau_so"/>
<variable units="ms" name="tau_so1"/>
<variable units="ms" name="tau_so2"/>
<variable units="dimensionless" name="k_so"/>
<variable units="dimensionless" name="u_so"/>
<variable units="dimensionless" public_interface="in" name="r"/>
<variable units="dimensionless" public_interface="in" name="p"/>
<variable units="ms" public_interface="in" name="time"/>
<variable units="dimensionless" public_interface="in" name="u"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>tau_o1</ci>
<piecewise>
<piece>
<cn cellml:units="ms">400</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">470</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">410</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_o2</ci>
<piecewise>
<piece>
<cn cellml:units="ms">6</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">6</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">7</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_so1</ci>
<piecewise>
<piece>
<cn cellml:units="ms">30.02</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">40</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">91</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_so2</ci>
<piecewise>
<piece>
<cn cellml:units="ms">0.996</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">1.2</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">0.8</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>k_so</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">2.046</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">2</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">2.1</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>u_so</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0.65</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">0.65</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0.6</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_o</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>r</ci>
</apply>
<ci>tau_o1</ci>
</apply>
<apply>
<times/>
<ci>r</ci>
<ci>tau_o2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_so</ci>
<apply>
<plus/>
<ci>tau_so1</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>tau_so2</ci>
<ci>tau_so1</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<tanh/>
<apply>
<times/>
<ci>k_so</ci>
<apply>
<minus/>
<ci>u</ci>
<ci>u_so</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_so</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>u</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>p</ci>
</apply>
</apply>
<ci>tau_o</ci>
</apply>
<apply>
<divide/>
<ci>p</ci>
<ci>tau_so</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_inward_current">
<variable units="per_ms" public_interface="out" name="J_si"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="epi"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="endo"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="mcell"/>
<variable units="ms" name="tau_si"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="p"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="r"/>
<variable units="ms" public_interface="in" private_interface="out" name="time"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="u"/>
<variable units="dimensionless" private_interface="in" name="w"/>
<variable units="dimensionless" private_interface="in" name="s"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>tau_si</ci>
<piecewise>
<piece>
<cn cellml:units="ms">1.8875</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">2.9013</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">3.3849</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>J_si</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>p</ci>
</apply>
<ci>w</ci>
<ci>s</ci>
</apply>
<ci>tau_si</ci>
</apply>
</apply>
</math>
</component>
<component name="slow_inward_current_w_gate">
<variable units="dimensionless" public_interface="out" name="w" initial_value="1"/>
<variable units="dimensionless" public_interface="in" name="epi"/>
<variable units="dimensionless" public_interface="in" name="endo"/>
<variable units="dimensionless" public_interface="in" name="mcell"/>
<variable units="dimensionless" name="w_inf"/>
<variable units="ms" name="tau_winf"/>
<variable units="dimensionless" name="wstar_inf"/>
<variable units="ms" name="tau_w_minus"/>
<variable units="ms" name="tau_w1_minus"/>
<variable units="ms" name="tau_w2_minus"/>
<variable units="dimensionless" name="k_w_minus"/>
<variable units="dimensionless" name="u_w_minus"/>
<variable units="ms" name="tau_w_plus"/>
<variable units="dimensionless" public_interface="in" name="r"/>
<variable units="dimensionless" public_interface="in" name="u"/>
<variable units="ms" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>tau_winf</ci>
<piecewise>
<piece>
<cn cellml:units="ms">0.07</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">0.0273</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">0.01</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>wstar_inf</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0.94</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">0.78</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0.5</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_w1_minus</ci>
<piecewise>
<piece>
<cn cellml:units="ms">60</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">6</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">70</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_w2_minus</ci>
<piecewise>
<piece>
<cn cellml:units="ms">15</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">140</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">8</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>k_w_minus</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">65</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">200</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">200</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>u_w_minus</ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless">0.03</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="dimensionless">0.016</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless">0.016</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_w_plus</ci>
<piecewise>
<piece>
<cn cellml:units="ms">200</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">280</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">280</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>w_inf</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>r</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>u</ci>
<cn cellml:units="ms">1</cn>
</apply>
<ci>tau_winf</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>r</ci>
<ci>wstar_inf</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_w_minus</ci>
<apply>
<plus/>
<ci>tau_w1_minus</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>tau_w2_minus</ci>
<ci>tau_w1_minus</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<tanh/>
<apply>
<times/>
<ci>k_w_minus</ci>
<apply>
<minus/>
<ci>u</ci>
<ci>u_w_minus</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>w</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>r</ci>
</apply>
<apply>
<minus/>
<ci>w_inf</ci>
<ci>w</ci>
</apply>
</apply>
<ci>tau_w_minus</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>r</ci>
<ci>w</ci>
</apply>
<ci>tau_w_plus</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_inward_current_s_gate">
<variable units="dimensionless" public_interface="out" name="s" initial_value="0"/>
<variable units="dimensionless" public_interface="in" name="epi"/>
<variable units="dimensionless" public_interface="in" name="endo"/>
<variable units="dimensionless" public_interface="in" name="mcell"/>
<variable units="ms" name="tau_s"/>
<variable units="ms" name="tau_s1" initial_value="2.7342"/>
<variable units="ms" name="tau_s2"/>
<variable units="dimensionless" name="k_s" initial_value="2.0994"/>
<variable units="dimensionless" name="u_s" initial_value="0.9087"/>
<variable units="dimensionless" public_interface="in" name="p"/>
<variable units="dimensionless" public_interface="in" name="u"/>
<variable units="ms" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>tau_s2</ci>
<piecewise>
<piece>
<cn cellml:units="ms">16</cn>
<apply>
<eq/>
<ci>epi</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<piece>
<cn cellml:units="ms">2</cn>
<apply>
<eq/>
<ci>endo</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="ms">4</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>tau_s</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>p</ci>
</apply>
<ci>tau_s1</ci>
</apply>
<apply>
<times/>
<ci>p</ci>
<ci>tau_s2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<tanh/>
<apply>
<times/>
<ci>k_s</ci>
<apply>
<minus/>
<ci>u</ci>
<ci>u_s</ci>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>s</ci>
</apply>
<ci>tau_s</ci>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="fast_inward_current">
<component_ref component="fast_inward_current_v_gate"/>
</component_ref>
<component_ref component="slow_outward_current"/>
<component_ref component="slow_inward_current">
<component_ref component="slow_inward_current_w_gate"/>
<component_ref component="slow_inward_current_s_gate"/>
</component_ref>
<component_ref component="m"/>
<component_ref component="p"/>
<component_ref component="q"/>
<component_ref component="r"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="fast_inward_current">
<component_ref component="fast_inward_current_v_gate"/>
</component_ref>
<component_ref component="slow_inward_current">
<component_ref component="slow_inward_current_w_gate"/>
<component_ref component="slow_inward_current_s_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="m"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="p"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="q"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="epi" variable_1="epi"/>
<map_variables variable_2="endo" variable_1="endo"/>
<map_variables variable_2="mcell" variable_1="mcell"/>
</connection>
<connection>
<map_components component_2="environment" component_1="r"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="epi" variable_1="epi"/>
<map_variables variable_2="endo" variable_1="endo"/>
<map_variables variable_2="mcell" variable_1="mcell"/>
</connection>
<connection>
<map_components component_2="environment" component_1="fast_inward_current"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="epi" variable_1="epi"/>
<map_variables variable_2="endo" variable_1="endo"/>
<map_variables variable_2="mcell" variable_1="mcell"/>
</connection>
<connection>
<map_components component_2="environment" component_1="slow_outward_current"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="epi" variable_1="epi"/>
<map_variables variable_2="endo" variable_1="endo"/>
<map_variables variable_2="mcell" variable_1="mcell"/>
</connection>
<connection>
<map_components component_2="environment" component_1="slow_inward_current"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="epi" variable_1="epi"/>
<map_variables variable_2="endo" variable_1="endo"/>
<map_variables variable_2="mcell" variable_1="mcell"/>
</connection>
<connection>
<map_components component_2="fast_inward_current" component_1="membrane"/>
<map_variables variable_2="u" variable_1="u"/>
<map_variables variable_2="J_fi" variable_1="J_fi"/>
</connection>
<connection>
<map_components component_2="slow_outward_current" component_1="membrane"/>
<map_variables variable_2="u" variable_1="u"/>
<map_variables variable_2="J_so" variable_1="J_so"/>
</connection>
<connection>
<map_components component_2="slow_inward_current" component_1="membrane"/>
<map_variables variable_2="u" variable_1="u"/>
<map_variables variable_2="J_si" variable_1="J_si"/>
</connection>
<connection>
<map_components component_2="m" component_1="membrane"/>
<map_variables variable_2="u" variable_1="u"/>
</connection>
<connection>
<map_components component_2="p" component_1="membrane"/>
<map_variables variable_2="u" variable_1="u"/>
</connection>
<connection>
<map_components component_2="q" component_1="membrane"/>
<map_variables variable_2="u" variable_1="u"/>
</connection>
<connection>
<map_components component_2="r" component_1="membrane"/>
<map_variables variable_2="u" variable_1="u"/>
</connection>
<connection>
<map_components component_2="fast_inward_current" component_1="m"/>
<map_variables variable_2="m" variable_1="m"/>
<map_variables variable_2="u_m" variable_1="u_m"/>
</connection>
<connection>
<map_components component_2="fast_inward_current" component_1="p"/>
<map_variables variable_2="u_p" variable_1="u_p"/>
</connection>
<connection>
<map_components component_2="slow_outward_current" component_1="p"/>
<map_variables variable_2="p" variable_1="p"/>
</connection>
<connection>
<map_components component_2="slow_outward_current" component_1="r"/>
<map_variables variable_2="r" variable_1="r"/>
</connection>
<connection>
<map_components component_2="slow_inward_current" component_1="p"/>
<map_variables variable_2="p" variable_1="p"/>
</connection>
<connection>
<map_components component_2="slow_inward_current" component_1="r"/>
<map_variables variable_2="r" variable_1="r"/>
</connection>
<connection>
<map_components component_2="fast_inward_current" component_1="q"/>
<map_variables variable_2="q" variable_1="q"/>
<map_variables variable_2="u_q" variable_1="u_q"/>
</connection>
<connection>
<map_components component_2="fast_inward_current_v_gate" component_1="fast_inward_current"/>
<map_variables variable_2="v" variable_1="v"/>
<map_variables variable_2="q" variable_1="q"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="epi" variable_1="epi"/>
<map_variables variable_2="endo" variable_1="endo"/>
<map_variables variable_2="mcell" variable_1="mcell"/>
<map_variables variable_2="m" variable_1="m"/>
<map_variables variable_2="u" variable_1="u"/>
<map_variables variable_2="u_q" variable_1="u_q"/>
</connection>
<connection>
<map_components component_2="slow_inward_current_w_gate" component_1="slow_inward_current"/>
<map_variables variable_2="w" variable_1="w"/>
<map_variables variable_2="r" variable_1="r"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="epi" variable_1="epi"/>
<map_variables variable_2="endo" variable_1="endo"/>
<map_variables variable_2="mcell" variable_1="mcell"/>
<map_variables variable_2="u" variable_1="u"/>
</connection>
<connection>
<map_components component_2="slow_inward_current_s_gate" component_1="slow_inward_current"/>
<map_variables variable_2="s" variable_1="s"/>
<map_variables variable_2="p" variable_1="p"/>
<map_variables variable_2="u" variable_1="u"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="epi" variable_1="epi"/>
<map_variables variable_2="endo" variable_1="endo"/>
<map_variables variable_2="mcell" variable_1="mcell"/>
</connection>
<rdf:RDF>
<rdf:Seq rdf:about="rdf:#679f10ee-12d7-4474-ad1d-8685b5aab3e9">
<rdf:li rdf:resource="rdf:#927faa0c-8772-4c18-baec-973b868c6ce4"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#f3061df5-2fd2-4f4c-a840-4c8cbf3c99b4">
<dcterms:modified rdf:resource="rdf:#acd5ae7b-56ad-43e5-841e-3e1b69b05f73"/>
<rdf:value>updated curation status</rdf:value>
<cmeta:modifier rdf:resource="rdf:#1eed0db6-d2ea-43dc-948b-e0f048014981"/>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher/>
<cmeta:comment rdf:resource="rdf:#30daf8a5-05a4-4492-9dea-3cf509e5072d"/>
<dcterms:created rdf:resource="rdf:#2ab16259-22dd-467a-8d88-da85a853cb66"/>
<dc:creator rdf:resource="rdf:#b3ec4e5a-01ac-4f03-a847-efbdc8809aaa"/>
<cmeta:modification rdf:resource="rdf:#f3061df5-2fd2-4f4c-a840-4c8cbf3c99b4"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fa433c77-a70b-477b-81b6-0121d11471b1">
<dc:creator rdf:resource="rdf:#b21c1c2f-13f1-4e89-b2da-d0ece71736da"/>
<rdf:value>This model of a ventricular cell was taken from Alfonso Bueno's 2007 PhD thesis, which can be found at http://matematicas.uclm.es/abueno/webpdfs/abueno_thesis.pdf (as of 2008-06-16)</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b21c1c2f-13f1-4e89-b2da-d0ece71736da">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4c0b675b-b50c-40de-8026-022a54b4ecf7">
<vCard:Given>Alfonso</vCard:Given>
<vCard:Family>Bueno</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#acd5ae7b-56ad-43e5-841e-3e1b69b05f73">
<dcterms:W3CDTF>2009-05-28T16:10:40+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5db67e3d-fa76-42da-bc41-f89649c9aa62">
<bqs:Pubmed_id/>
<bqs:JournalArticle rdf:resource="rdf:#f2ae2bd3-7f9f-4812-9700-23139002c009"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#48ebd61a-83c5-44df-9400-4c2b3a81fb6b">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b4f1cc2b-e1e9-4652-8ddf-2cca4ab30114">
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
<rdf:value>penny.noble@dpag.ox.ac.uk</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1eed0db6-d2ea-43dc-948b-e0f048014981">
<vCard:N rdf:resource="rdf:#833e0cc1-5ecb-4216-9086-0ba193c74ade"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0915ccb0-e15d-49b2-bb73-e058ea84d8ae">
<dc:title/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2ab16259-22dd-467a-8d88-da85a853cb66">
<dcterms:W3CDTF>2008-02-28T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="#bueno_2007">
<bqs:reference rdf:resource="rdf:#5db67e3d-fa76-42da-bc41-f89649c9aa62"/>
<cmeta:comment rdf:resource="rdf:#fa433c77-a70b-477b-81b6-0121d11471b1"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#927faa0c-8772-4c18-baec-973b868c6ce4">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#4c0b675b-b50c-40de-8026-022a54b4ecf7"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#833e0cc1-5ecb-4216-9086-0ba193c74ade">
<vCard:Given>James</vCard:Given>
<vCard:Family>Lawson</vCard:Family>
<vCard:Other>Richard</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f2ae2bd3-7f9f-4812-9700-23139002c009">
<dc:creator rdf:resource="rdf:#679f10ee-12d7-4474-ad1d-8685b5aab3e9"/>
<dc:title>Mathematical modeling and spectral simulation of genetic diseases in the human heart</dc:title>
<bqs:volume/>
<bqs:first_page/>
<bqs:Journal rdf:resource="rdf:#0915ccb0-e15d-49b2-bb73-e058ea84d8ae"/>
<dcterms:issued rdf:resource="rdf:#0e6ae837-14b7-41fc-bea4-c1bb9aa16225"/>
<bqs:last_page/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4b1baa5d-3259-40f1-b3c4-3b25006396d9">
<vCard:Orgname>The University of Oxford</vCard:Orgname>
<vCard:Orgunit>Department of Physiology, Anatomy & Genetics</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0b09c387-f86a-4586-8e47-8f493a62f324">
<vCard:Given>Penny</vCard:Given>
<vCard:Family>Noble</vCard:Family>
<vCard:Other/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b3ec4e5a-01ac-4f03-a847-efbdc8809aaa">
<vCard:ORG rdf:resource="rdf:#4b1baa5d-3259-40f1-b3c4-3b25006396d9"/>
<vCard:EMAIL rdf:resource="rdf:#b4f1cc2b-e1e9-4652-8ddf-2cca4ab30114"/>
<vCard:N rdf:resource="rdf:#0b09c387-f86a-4586-8e47-8f493a62f324"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0e6ae837-14b7-41fc-bea4-c1bb9aa16225">
<dcterms:W3CDTF>2007-05-00 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#30daf8a5-05a4-4492-9dea-3cf509e5072d">
<dc:creator rdf:resource="rdf:#48ebd61a-83c5-44df-9400-4c2b3a81fb6b"/>
<rdf:value>This CellML model is known to run in both COR and PCEnv to recreate the results of the original model. The units have been checked and are consistent. Model has been curated using Bueno 2007 thesis.</rdf:value>
</rdf:Description>
</rdf:RDF>
</model>
