Location: Fenton, Karma, 1998 @ 5925f43454dd / fenton_karma_1998_MLR1.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2010-04-30 13:43:03+12:00
Desc:
Remade session file with integrated clickable SVG
Permanent Source URI:
https://models.cellml.org/workspace/5bb/rawfile/5925f43454dd91b8e055bee3b59391ad10631382/fenton_karma_1998_MLR1.cellml

<?xml version="1.0"?>
<!--
This CellML file was generated on 30/10/2007 at 13:57:00 using:

COR (0.9.31.801)
Copyright 2002-2007 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#" xml:base="" cmeta:id="fenton_karma_1998_MLR-1" name="fenton_karma_1998_MLR-1">
  <documentation xmlns="http://cellml.org/tmp-documentation">
    <article>
      <articleinfo>
        <title>A Simplified Ventricular Myocyte Model</title>
        <author>
          <firstname>Penny</firstname>
          <surname>Noble</surname>
          <affiliation>
            <shortaffil>Oxford University</shortaffil>
          </affiliation>
        </author>
      </articleinfo>
      <section id="sec_status">
        <title>Model Status</title>
        <para>
            This version of this model has had a stimulus protocol added (James Lawson, 15/06/07) to allow it to simulate trains of action potentials. The previous version was created by Penny Noble of Oxford University, and is known to read in COR and PCEnv. The parameter values used in this variant (MLR-1) of the Fenton-Karma model are consistent with the modified (1994) Luo-Rudy model (see Table 1 of the 1998 model errata).  Simulations of this CellML model can be run using CMISS.
          </para>
      </section>
      <sect1 id="sec_structure">
        <title>Model Structure</title>
        <para>
Ventricular fibrillation (VF) is a disorganised electrical wave activity that disrupts the regular and synchronised contraction of the ventricular muscle and consequently destroys the main pumping function of the heart.  Research to date has suggested that electrical vortices are the mechanism underlying VF.  These are manisfested in the 2D heart as spiral waves, and in the 3D heart as scroll waves of action potential.  Compared with the extensive knowledge available on vortex filament behaviour in an isotropic medium, relatively little is known about how the vortices behave in the anisotropic ventricular muscle.  This 1998 publication by Fenton and Karma (fully referenced below) explores the dynamics of the vortex filament in continuous myocardium via numerical simulation.    </para>
        <para>
Over the past 25 to 30 years, mathematical models that describe ventricular action potential have become increasingly complex as new experimental data has become available and has been incorporated into the mathematical equations.  Although these complex models are more realistic, they are also computationally expensive to run, and isolating subsets of essential parameters from the model is difficult.  One traditional method for avoiding this complexity is to use simplified models such as <ulink url="${HTML_EXMPL_FN_SIMPLE}">the FitzHugh-Nagumo model, 1961</ulink>.  However, these simplified models have been criticised as being too simple and as not having the capacity to fully capture certain important features of the ventricular action potential.
</para>
        <para>
The modelling approach that Fenton and Karma take is to use a simplified ionic model of ventricular action potential with three membrane currents.  This model retains enough detail to quantitatively reproduce the behaviour of the ventricular action potential captured by the more complex ionic models of cardiac action potential (such as <ulink url="${HTML_EXMPL_BR_MODEL}">the Beeler-Reuter 1977 model</ulink>, and <ulink url="${HTML_EXMPL_LR_I_MODEL}">the original Luo-Rudy 1991 model</ulink>), but it is less computationally expensive than these other models.
</para>
        <para>
The three currents in the Fenton-Karma model are: 
</para>
        <itemizedlist>
          <listitem>
            <para>
              <emphasis>I<subscript>fi</subscript>
              </emphasis>, a fast inward current which corresponds to the <emphasis>I<subscript>Na</subscript>
              </emphasis> current</para>
          </listitem>;
  <listitem>
            <para>
              <emphasis>I<subscript>so</subscript>
              </emphasis>, a slow outward current which corresponds to the <emphasis>I<subscript>K</subscript>
              </emphasis> current</para>
          </listitem>; and     
  <listitem>
            <para>
              <emphasis>I<subscript>si</subscript>
              </emphasis>, a slow inward current which corresponds to the <emphasis>I<subscript>Ca</subscript>
              </emphasis> current</para>
          </listitem>.
</itemizedlist>
        <para>
(see <xref linkend="fig_cell_diagram"/> below).
The authors choose to use these labels as opposed to Na, K and Ca, as a reminder that <emphasis>I<subscript>fi</subscript>
          </emphasis>, <emphasis>I<subscript>so</subscript>
          </emphasis>, and  <emphasis>I<subscript>si</subscript>
          </emphasis> do not actually represent measured currents, but only their activation, inactivation and reactivation dynamics which are needed to quantitatively reproduce restitution properties.
</para>
        <para>
Model validation indicates that this simplified model is able to faithfully reproduce the 2D patterns of reentry of the more complex models.
</para>
        <para>
The complete original paper reference is cited below:
</para>
        <para>
          Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation, Flavio Fenton and Alain Karma, 1998, <emphasis>Chaos</emphasis>,8, 20-47.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=12779708&amp;dopt=Abstract">PubMed ID: 12779708</ulink>
        </para>
        <informalfigure float="0" id="fig_cell_diagram">
          <mediaobject>
            <imageobject>
              <objectinfo>
                <title>cell diagram</title>
              </objectinfo>
              <imagedata fileref="fenton_1998a.png"/>
            </imageobject>
          </mediaobject>
          <caption>A schematic diagram of the three ionic currents described by the Fenton-Karma model of a ventricular myocyte.</caption>
        </informalfigure>
        <para>
The Fenton-Karma model has been described here in CellML (the raw CellML description of the Fenton-Karma model can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>).          
</para>
      </sect1>
    </article>
  </documentation>
   <units name="ms">
      <unit prefix="milli" units="second"/>
   </units>
   <units name="per_ms">
      <unit exponent="-1" prefix="milli" units="second"/>
   </units>
   <units name="mV">
      <unit prefix="milli" units="volt"/>
   </units>
   <units name="per_mV">
      <unit exponent="-1" prefix="milli" units="volt"/>
   </units>
   <units name="per_mV_ms">
      <unit exponent="-1" units="mV"/>
      <unit exponent="-1" units="ms"/>
   </units>
   <units name="mS_per_cm2">
      <unit prefix="milli" units="siemens"/>
      <unit exponent="-2" prefix="centi" units="metre"/>
   </units>
   <units name="uF_per_cm2">
      <unit prefix="micro" units="farad"/>
      <unit exponent="-2" prefix="centi" units="metre"/>
   </units>
   <units name="uA_per_cm2">
      <unit prefix="micro" units="ampere"/>
      <unit exponent="-2" prefix="centi" units="metre"/>
   </units>
   <units name="concentration_units">
      <unit prefix="nano" units="mole"/>
      <unit exponent="-3" prefix="milli" units="metre"/>
   </units>
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="ms"/>
   </component>

 <component name="stimulus_protocol">
    <variable cmeta:id="stimulus_protocol_Istim" name="Istim" public_interface="out" units="per_ms"/>
    <variable cmeta:id="stimulus_protocol_IstimStart" initial_value="100" name="IstimStart" units="ms"/>
    <variable cmeta:id="stimulus_protocol_IstimEnd" initial_value="50000" name="IstimEnd" units="ms"/>
    <variable cmeta:id="stimulus_protocol_IstimAmplitude" initial_value="-0.2" name="IstimAmplitude" units="per_ms"/>
    <variable cmeta:id="stimulus_protocol_IstimPeriod" initial_value="1000" name="IstimPeriod" units="ms"/>
    <variable cmeta:id="stimulus_protocol_IstimPulseDuration" initial_value="1" name="IstimPulseDuration" units="ms"/>
    <variable name="time" public_interface="in" units="ms"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML" id="Stimulus_Current_Equation">
      <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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_ms">0</cn>
          </otherwise>
        </piecewise>
      </apply>
    </math>
    
  </component>

<connection>
    <map_components component_1="membrane" component_2="stimulus_protocol"/>
    <map_variables variable_1="J_stim" variable_2="Istim"/>
  </connection>
  <connection>
    <map_components component_1="environment" component_2="stimulus_protocol"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>

   <component name="membrane">
      <variable initial_value="0" name="u" public_interface="out" units="dimensionless"/>
      <variable initial_value="1" name="Cm" public_interface="out" units="uF_per_cm2"/>
      <variable cmeta:id="membrane_Vm" name="Vm" units="mV"/>
      <variable initial_value="-85" name="V_0" units="mV"/>
      <variable initial_value="15" name="V_fi" units="mV"/>
      <variable name="time" public_interface="in" units="ms"/>
      <variable name="J_fi" public_interface="in" units="per_ms"/>
      <variable name="J_so" public_interface="in" units="per_ms"/>
      <variable name="J_si" public_interface="in" units="per_ms"/>
      <variable name="J_stim" public_interface="in" units="per_ms"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         
         <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="p">
      <variable name="p" public_interface="out" units="dimensionless"/>
      <variable initial_value="0.13" name="u_c" public_interface="out" units="dimensionless"/>
      <variable name="u" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="ms"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>p</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
                  <apply>
                     <lt/>
                     <ci>u</ci>
                     <ci>u_c</ci>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="q">
      <variable name="q" public_interface="out" units="dimensionless"/>
      <variable initial_value="0" name="u_v" units="dimensionless"/>
      <variable name="u" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="ms"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>q</ci>
            <piecewise>
               <piece>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
                  <apply>
                     <lt/>
                     <ci>u</ci>
                     <ci>u_v</ci>
                  </apply>
               </piece>
               <otherwise>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="fast_inward_current">
      <variable cmeta:id="fast_inward_current_J_fi" name="J_fi" public_interface="out" units="per_ms"/>
      <variable name="tau_d" units="ms"/>
      <variable initial_value="5.8" name="g_fi_max" units="mS_per_cm2"/>
      <variable name="Cm" public_interface="in" units="uF_per_cm2"/>
      <variable name="p" private_interface="out" public_interface="in" units="dimensionless"/>
      <variable name="q" private_interface="out" public_interface="in" units="dimensionless"/>
      <variable name="u_c" public_interface="in" units="dimensionless"/>
      <variable name="time" private_interface="out" public_interface="in" units="ms"/>
      <variable name="u" private_interface="out" public_interface="in" units="dimensionless"/>
      <variable name="v" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau_d</ci>
            <apply>
               <divide/>
               <ci>Cm</ci>
               <ci>g_fi_max</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_fi</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>v</ci>
                  </apply>
                  <ci>p</ci>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>u</ci>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>u</ci>
                     <ci>u_c</ci>
                  </apply>
               </apply>
               <ci>tau_d</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_inward_current_v_gate">
      <variable initial_value="1" name="v" public_interface="out" units="dimensionless"/>
      <variable name="tau_v_minus" units="ms"/>
      <variable initial_value="18.2" name="tau_v1_minus" units="ms"/>
      <variable initial_value="18.2" name="tau_v2_minus" units="ms"/>
      <variable initial_value="10" name="tau_v_plus" units="ms"/>
      <variable name="q" public_interface="in" units="dimensionless"/>
      <variable name="p" public_interface="in" units="dimensionless"/>
      <variable name="u" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="ms"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau_v_minus</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <ci>q</ci>
                  <ci>tau_v1_minus</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <ci>q</ci>
                  </apply>
                  <ci>tau_v2_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 xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <ci>p</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <ci>v</ci>
                     </apply>
                  </apply>
                  <ci>tau_v_minus</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>p</ci>
                     <ci>v</ci>
                  </apply>
                  <ci>tau_v_plus</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_outward_current">
      <variable cmeta:id="slow_outward_current_J_so" name="J_so" public_interface="out" units="per_ms"/>
      <variable initial_value="12.5" name="tau_0" units="ms"/>
      <variable initial_value="130" name="tau_r" units="ms"/>
      <variable name="p" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="ms"/>
      <variable name="u" public_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_so</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>u</ci>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <ci>p</ci>
                     </apply>
                  </apply>
                  <ci>tau_0</ci>
               </apply>
               <apply>
                  <divide/>
                  <ci>p</ci>
                  <ci>tau_r</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_inward_current">
      <variable cmeta:id="slow_inward_current_J_si" name="J_si" public_interface="out" units="per_ms"/>
      <variable initial_value="127" name="tau_si" units="ms"/>
      <variable initial_value="0.85" name="u_csi" units="dimensionless"/>
      <variable initial_value="10" name="k" units="dimensionless"/>
      <variable name="p" private_interface="out" public_interface="in" units="dimensionless"/>
      <variable name="time" private_interface="out" public_interface="in" units="ms"/>
      <variable name="u" private_interface="out" public_interface="in" units="dimensionless"/>
      <variable name="w" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_si</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>w</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <tanh/>
                        <apply>
                           <times/>
                           <ci>k</ci>
                           <apply>
                              <minus/>
                              <ci>u</ci>
                              <ci>u_csi</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  <ci>tau_si</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_inward_current_w_gate">
      <variable initial_value="1" name="w" public_interface="out" units="dimensionless"/>
      <variable initial_value="80" name="tau_w_minus" units="ms"/>
      <variable initial_value="1020" name="tau_w_plus" units="ms"/>
      <variable name="p" public_interface="in" units="dimensionless"/>
      <variable name="u" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="ms"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>w</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <ci>p</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                        <ci>w</ci>
                     </apply>
                  </apply>
                  <ci>tau_w_minus</ci>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>p</ci>
                     <ci>w</ci>
                  </apply>
                  <ci>tau_w_plus</ci>
               </apply>
            </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_ref component="p"/>
         <component_ref component="q"/>
      </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>
   </group>
   <connection>
      <map_components component_1="membrane" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="p" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="q" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="fast_inward_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="slow_outward_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="membrane" component_2="fast_inward_current"/>
      <map_variables variable_1="u" variable_2="u"/>
      <map_variables variable_1="J_fi" variable_2="J_fi"/>
      <map_variables variable_1="Cm" variable_2="Cm"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="slow_outward_current"/>
      <map_variables variable_1="u" variable_2="u"/>
      <map_variables variable_1="J_so" variable_2="J_so"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="slow_inward_current"/>
      <map_variables variable_1="u" variable_2="u"/>
      <map_variables variable_1="J_si" variable_2="J_si"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="p"/>
      <map_variables variable_1="u" variable_2="u"/>
   </connection>
   <connection>
      <map_components component_1="membrane" component_2="q"/>
      <map_variables variable_1="u" variable_2="u"/>
   </connection>
   <connection>
      <map_components component_1="p" component_2="fast_inward_current"/>
      <map_variables variable_1="p" variable_2="p"/>
      <map_variables variable_1="u_c" variable_2="u_c"/>
   </connection>
   <connection>
      <map_components component_1="p" component_2="slow_outward_current"/>
      <map_variables variable_1="p" variable_2="p"/>
   </connection>
   <connection>
      <map_components component_1="p" component_2="slow_inward_current"/>
      <map_variables variable_1="p" variable_2="p"/>
   </connection>
   <connection>
      <map_components component_1="q" component_2="fast_inward_current"/>
      <map_variables variable_1="q" variable_2="q"/>
   </connection>
   <connection>
      <map_components component_1="fast_inward_current" component_2="fast_inward_current_v_gate"/>
      <map_variables variable_1="v" variable_2="v"/>
      <map_variables variable_1="p" variable_2="p"/>
      <map_variables variable_1="q" variable_2="q"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="u" variable_2="u"/>
   </connection>
   <connection>
      <map_components component_1="slow_inward_current" component_2="slow_inward_current_w_gate"/>
      <map_variables variable_1="w" variable_2="w"/>
      <map_variables variable_1="p" variable_2="p"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="u" variable_2="u"/>
   </connection>



<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <rdf:Seq rdf:about="rdf:#e1a6af36-d357-4fa2-b42e-2f8fff4d3d2b">
    <rdf:li rdf:resource="rdf:#37c5369e-8152-4b87-9fc5-50cb63352246"/>
    <rdf:li rdf:resource="rdf:#48710d37-a154-47e1-ab64-b5cf01241a50"/>
  </rdf:Seq>
  <rdf:Description rdf:about="">
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation (MLR-1 Model)</dc:title>
    <dc:publisher xmlns:dc="http://purl.org/dc/elements/1.1/">Department of Physiology, Anatomy &amp; Genetics, University of Oxford</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#9527679f-e896-4978-a838-d117106b14da"/>
    <dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#7d59717d-b702-435a-8fdd-a8f6511df88a"/>
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#51687ec5-6bae-481b-a7be-db061ae06540"/>
    <cmeta:modification rdf:resource="rdf:#8cda4f63-bf5d-46ec-9649-234812d8b89c"/>
    <cmeta:modification rdf:resource="rdf:#07da05f9-d914-481f-a55e-6a95e47ede81"/>
    <cmeta:modification rdf:resource="rdf:#a5dccb73-684b-41a7-a2c9-ac3edbfea151"/>
    <cmeta:modification rdf:resource="rdf:#db07fa98-6fcf-4205-ab91-ca0f05d8aba1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3f034e7c-4448-437d-93aa-5b46d9e96c26">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-06-20T14:47:44+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a5dccb73-684b-41a7-a2c9-ac3edbfea151">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#4f4084db-8927-46ca-8e3b-3179730f95c5"/>
    <rdf:value>This file has been rebuilt and validated by Penny Noble of Oxford University</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#7952f2dd-a5ea-48a1-8a9a-699a729b35c1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#21cf87eb-1048-4cba-ba19-fc65a60492e5">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lawson</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Richard</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#db07fa98-6fcf-4205-ab91-ca0f05d8aba1">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#3f034e7c-4448-437d-93aa-5b46d9e96c26"/>
    <rdf:value>Some cmeta:id's were added to several variables to allow creation of a PCEnv session file.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#6f9413dd-4364-4ac8-8b3f-9356f0ce9a35"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7952f2dd-a5ea-48a1-8a9a-699a729b35c1">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#92d3b138-0307-4a3c-bf95-27dc11a42608"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a6f6a513-d9d4-450d-9d8f-c49079452a18">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8a60798b-5b2c-451c-bd7d-2493bf59fad1">
    <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:#92d3b138-0307-4a3c-bf95-27dc11a42608">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lawson</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">R</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#55f66b44-9410-453e-a73a-5c28a7eee008">
    <vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6f9413dd-4364-4ac8-8b3f-9356f0ce9a35">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#521a22ca-2fd1-4697-8ca8-4614ac566c61"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#48710d37-a154-47e1-ab64-b5cf01241a50">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#0e0f2633-8c14-4205-9a2a-d04cf8ddfdaf"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#521a22ca-2fd1-4697-8ca8-4614ac566c61">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lawson</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Richard</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7d59717d-b702-435a-8fdd-a8f6511df88a">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2006-05-10T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#0816f6c9-7db8-477e-a104-8cfe4b1d7b79"/>
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#c3b52c74-8a66-49c1-b04c-d897edfeb748"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0e0f2633-8c14-4205-9a2a-d04cf8ddfdaf">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">A</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Karma</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#48db750a-59d2-429b-a88e-65fb7ce6c511">
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Chaos</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fa164a73-6a68-44ec-8f3f-f784c68664e9">
    <vCard:Orgname xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Oxford University</vCard:Orgname>
    <vCard:Orgunit xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Department of Physiology, Anatomy &amp; Genetics</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$O7UWt">
    <ns7:endingValue xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#">10000</ns7:endingValue>
    <ns7:nonstandard-pointDensity xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#">100000</ns7:nonstandard-pointDensity>
    <ns7:maximumStepSize xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#">0.1</ns7:maximumStepSize>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#62c4c0d3-764a-45c7-b177-d0d6954d30a2">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>unknown</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9527679f-e896-4978-a838-d117106b14da">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#4a89c100-bcbe-4a8c-887f-211e53c61633"/>
    <rdf:value>This version has been modified from version 03 to add a stimulus protocol to allow simulation of trains of action potentials. Version 03 was created and curated by Penny Noble of Oxford University.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0ac312cd-a0ec-4236-9f20-87f3632c492d">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#21cf87eb-1048-4cba-ba19-fc65a60492e5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4a89c100-bcbe-4a8c-887f-211e53c61633">
    <vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$M7UWt">
    <ns7:boundIntervals xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" rdf:resource="rdf:#$N7UWt"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3ee5c5ee-409b-4522-881a-4ef5fbc830a1">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Penny</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Noble</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0816f6c9-7db8-477e-a104-8cfe4b1d7b79">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2006-01-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#801bdca1-3eec-4250-a525-97486b03e091">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">F</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Fenton</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#37c5369e-8152-4b87-9fc5-50cb63352246">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#801bdca1-3eec-4250-a525-97486b03e091"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8cda4f63-bf5d-46ec-9649-234812d8b89c">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#2ef2a14b-9ad3-4513-8c77-b712cc17c9b8"/>
    <rdf:value>Units checked, curated.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#083f577e-689b-47b6-bd95-f0f2e81ad61f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c3b52c74-8a66-49c1-b04c-d897edfeb748">
    <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#da7c67aa-e852-48ab-a6f5-25c3a34aa1d6"/>
    <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#62c4c0d3-764a-45c7-b177-d0d6954d30a2"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#a6f6a513-d9d4-450d-9d8f-c49079452a18"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#85a5d21a-fb00-4c97-9dc4-29540bda0159">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Penny</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Noble</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e1ec28a9-bd5c-449e-a0b4-35cb2936625b">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">1998-03-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ea1ed633-e446-411e-a3db-f1f7de108c35">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#55f66b44-9410-453e-a73a-5c28a7eee008"/>
    <rdf:value>This version &amp; variant was created by Penny Noble of Oxford University and uses the parameters of the modified Luo-Rudy 1994 model.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#fenton_karma_1998_MLR-1">
    <ns7:simulation xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" rdf:resource="rdf:#$M7UWt"/>
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">unknown</dc:title>
    <cmeta:comment rdf:resource="rdf:#ea1ed633-e446-411e-a3db-f1f7de108c35"/>
    <bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#8d2a6af5-7074-45f9-ab67-2030458352ee"/>
<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
  <dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>Ventricular Myocyte</rdf:li>
        <rdf:li>myocardium</rdf:li>
        <rdf:li>electrophysiology</rdf:li>
        <rdf:li>vortex dynamics</rdf:li>
        <rdf:li>cardiac</rdf:li>
        <rdf:li>3d</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#07da05f9-d914-481f-a55e-6a95e47ede81">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#6de984eb-012e-4c41-94eb-5d1823a334a6"/>
    <rdf:value>A stimulus protocol was added to this model to allow simulation of a train of action potentials. The amplitude (-0.2 per second) and duration (1ms) of the stimulus was taken from the single stimulus in Penny's file.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#0ac312cd-a0ec-4236-9f20-87f3632c492d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4f4084db-8927-46ca-8e3b-3179730f95c5">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2006-05-10T14:40:00+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#83e582e9-20a9-4361-952e-6f34c72b1494">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#e1a6af36-d357-4fa2-b42e-2f8fff4d3d2b"/>
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation</dc:title>
    <bqs:volume xmlns:bqs="http://www.cellml.org/bqs/1.0#">8 1</bqs:volume>
    <bqs:first_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">20</bqs:first_page>
    <bqs:Journal xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#48db750a-59d2-429b-a88e-65fb7ce6c511"/>
    <dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#e1ec28a9-bd5c-449e-a0b4-35cb2936625b"/>
    <bqs:last_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">47</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8d2a6af5-7074-45f9-ab67-2030458352ee">
    <bqs:Pubmed_id xmlns:bqs="http://www.cellml.org/bqs/1.0#">12779708</bqs:Pubmed_id>
    <bqs:JournalArticle xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#83e582e9-20a9-4361-952e-6f34c72b1494"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$N7UWt">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$O7UWt"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#da7c67aa-e852-48ab-a6f5-25c3a34aa1d6">
    <vCard:Orgname xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Orgname>
    <vCard:Orgunit xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">unknown</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#51687ec5-6bae-481b-a7be-db061ae06540">
    <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#fa164a73-6a68-44ec-8f3f-f784c68664e9"/>
    <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#8a60798b-5b2c-451c-bd7d-2493bf59fad1"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#85a5d21a-fb00-4c97-9dc4-29540bda0159"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2ef2a14b-9ad3-4513-8c77-b712cc17c9b8">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-11-29T10:19:30+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#083f577e-689b-47b6-bd95-f0f2e81ad61f">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#3ee5c5ee-409b-4522-881a-4ef5fbc830a1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6de984eb-012e-4c41-94eb-5d1823a334a6">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-06-15T11:54:51+12:00</dcterms:W3CDTF>
  </rdf:Description>
</rdf:RDF>
</model>