Location: Fenton, Karma, 1998 @ 129998a1e151 / fenton_karma_1998_GP.cellml

Author:
AnandR <a.rampadarath@auckland.ac.nz>
Date:
2019-11-14 11:57:30+13:00
Desc:
replaced MLR-1 with MLR_1 in .cellml file
Permanent Source URI:
https://models.cellml.org/workspace/5bb/rawfile/129998a1e15191a21f08ce3f93b0d96e3c0b1ce0/fenton_karma_1998_GP.cellml

<?xml version="1.0"?>
<!--
This CellML file was generated on 30/10/2007 at 13:53:52 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#" cmeta:id="fenton_karma_1998_GP" name="fenton_karma_1998_GP">
  <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 CellML model runs in OpenCell and COR to reproduce the published results. This particular version of the CellML model has had a stimulus protocol added to allow it to simulate trains of action potentials. The parameter values used in this variant (GP) of the Fenton-Karma model are the GP values - the experimental steady-state curves extracted by Girouard et al. (Girouard et al., Circulation 93, pp. 603 to 613, 1996) from optical recordings of membrane potentials on the epicardial surface of the left ventricle of a guinea pig during plane wave pacing at fixed cycle lengths - 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>
ABSTRACT: Wave propagation in ventricular muscle is rendered highly anisotropic by the intramural rotation of the fiber. This rotational anisotropy is especially important because it can produce a twist of electrical vortices, which measures the rate of rotation (in degree/mm) of activation wavefronts in successive planes perpendicular to a line of phase singularity, or filament. This twist can then significantly alter the dynamics of the filament. This paper explores this dynamics via numerical simulation. After a review of the literature, we present modeling tools that include: (i) a simplified ionic model with three membrane currents that approximates well the restitution properties and spiral wave behavior of more complex ionic models of cardiac action potential (Beeler-Reuter and others), and (ii) a semi-implicit algorithm for the fast solution of monodomain cable equations with rotational anisotropy. We then discuss selected results of a simulation study of vortex dynamics in a parallelepipedal slab of ventricular muscle of varying wall thickness (S) and fiber rotation rate (theta(z)). The main finding is that rotational anisotropy generates a sufficiently large twist to destabilize a single transmural filament and cause a transition to a wave turbulent state characterized by a high density of chaotically moving filaments. This instability is manifested by the propagation of localized disturbances along the filament and has no previously known analog in isotropic excitable media. These disturbances correspond to highly twisted and distorted regions of filament, or "twistons," that create vortex rings when colliding with the natural boundaries of the ventricle. Moreover, when sufficiently twisted, these rings expand and create additional filaments by further colliding with boundaries. This instability mechanism is distinct from the commonly invoked patchy failure or wave breakup that is not observed here during the initial instability. For modified Beeler-Reuter-like kinetics with stable reentry in two dimensions, decay into turbulence occurs in the left ventricle in about one second above a critical wall thickness in the range of 4-6 mm that matches experiment. However this decay is suppressed by uniformly decreasing excitability. Specific experiments to test these results, and a method to characterize the filament density during fibrillation are discussed. Results are contrasted with other mechanisms of fibrillation and future prospects are summarized.
</para>
        
        <para>
The 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>
       
      </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="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.025" 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="8.7" 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="333" name="tau_v1_minus" units="ms"/>
      <variable initial_value="40" 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="25" 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="22.22" 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="65" name="tau_w_minus" units="ms"/>
      <variable initial_value="1000" 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>

  <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="10" 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"/>
    <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>
    <variable name="time" public_interface="in" units="ms"/>
  </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>



<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <rdf:Seq rdf:about="rdf:#f98bb645-1c00-4fb4-bc6a-c94ec8a151e4">
    <rdf:li rdf:resource="rdf:#4cc2a6e8-fe5d-4dff-b9cf-79240e6f68f1"/>
    <rdf:li rdf:resource="rdf:#f59bb867-c3ce-4a0f-ad8a-a7a2d6699623"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#b958a859-0963-402a-aa3e-2732447e9f0f">
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Chaos</dc:title>
  </rdf:Description>
  <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 (GP 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:#aaf0d05e-816b-4181-b469-ddcefbb160b1"/>
    <dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#1d9fb65c-7777-4566-9f8e-35a95a8e8301"/>
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#59e7ebd3-b0e2-408e-9aec-9cadacbc792c"/>
    <cmeta:modification rdf:resource="rdf:#8862e9e4-cd61-47e9-b5e7-5409066df28f"/>
    <cmeta:modification rdf:resource="rdf:#e6795066-26ba-484e-b0ce-095ffd5ba671"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#600fd408-69bb-43c8-bda2-68634c5a611d">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#38467f4c-fe9b-4f0c-bbd0-374f90c2669b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5220f008-a8c3-4fe8-9b3c-be74cc6a002b">
    <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:#949b9cea-614c-4e13-9ae8-e4873709d439"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$JGITr1">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$KGITr1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#92413ee5-bf1f-4b7d-a8f8-7b1155412c08">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#80b334ae-4929-4249-9c30-be799106431f"/>
    <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="rdf:#5c3d1fd1-82f8-49a2-b836-75076e88b21e">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#6aa5ef66-5683-4358-87cc-ab829bc372d8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#75928e86-db26-44fa-bdd3-3458f7b22e0e">
    <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:#0b41b465-d499-4d95-9c9e-4bd429578afd">
    <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:#3dde52e9-f7bb-45c1-b876-f7a10d53fe99">
    <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:#5e153a46-f6af-4507-a62b-b59f36d41e08">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-11-29T10:23:13+13:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1d9fb65c-7777-4566-9f8e-35a95a8e8301">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2006-05-10T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8862e9e4-cd61-47e9-b5e7-5409066df28f">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#5e153a46-f6af-4507-a62b-b59f36d41e08"/>
    <rdf:value>Units checked, curated.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#600fd408-69bb-43c8-bda2-68634c5a611d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d6a5a5ae-53c8-49be-8846-d274b30df388">
    <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:#7b02b72a-bdae-4e09-85d5-09e301b954f4">
    <vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#480dfac4-232b-4f29-96f6-f28e5e16adc5">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>unknown</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#26277ed1-a67f-4e6b-a690-0ad2d7fcf844"/>
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#e6c49398-f338-4b0b-8d9e-6e2d14c2e9e9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#949b9cea-614c-4e13-9ae8-e4873709d439">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#f98bb645-1c00-4fb4-bc6a-c94ec8a151e4"/>
    <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:#b958a859-0963-402a-aa3e-2732447e9f0f"/>
    <dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#bb6a26ce-0c58-4738-8477-9dd978132cac"/>
    <bqs:last_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">47</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#59e7ebd3-b0e2-408e-9aec-9cadacbc792c">
    <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#93de0b0a-21c9-4fce-8862-17484a4458d6"/>
    <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#0b41b465-d499-4d95-9c9e-4bd429578afd"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#01f1677b-6b4a-4b0c-b913-94a641fa0504"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$KGITr1">
    <ns7:endingValue xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#">10000</ns7:endingValue>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#80b334ae-4929-4249-9c30-be799106431f">
    <vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#93de0b0a-21c9-4fce-8862-17484a4458d6">
    <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:#4cc2a6e8-fe5d-4dff-b9cf-79240e6f68f1">
    <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:#3dde52e9-f7bb-45c1-b876-f7a10d53fe99"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f59bb867-c3ce-4a0f-ad8a-a7a2d6699623">
    <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:#75928e86-db26-44fa-bdd3-3458f7b22e0e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aaf0d05e-816b-4181-b469-ddcefbb160b1">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#7b02b72a-bdae-4e09-85d5-09e301b954f4"/>
    <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:#e6795066-26ba-484e-b0ce-095ffd5ba671">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#c4db9901-6076-4c8a-8942-c404735abbcc"/>
    <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:#5c3d1fd1-82f8-49a2-b836-75076e88b21e"/>
  </rdf:Description>
  <rdf:Description rdf:about="#fenton_karma_1998_GP">
    <ns7:simulation xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" rdf:resource="rdf:#$IGITr1"/>
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">unknown</dc:title>
    <cmeta:comment rdf:resource="rdf:#92413ee5-bf1f-4b7d-a8f8-7b1155412c08"/>
    <bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#5220f008-a8c3-4fe8-9b3c-be74cc6a002b"/>
<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:#d4e849d3-2fc4-495e-8714-8a3dacd846f7">
    <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:#01f1677b-6b4a-4b0c-b913-94a641fa0504">
    <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:#e6c49398-f338-4b0b-8d9e-6e2d14c2e9e9">
    <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#d4e849d3-2fc4-495e-8714-8a3dacd846f7"/>
    <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#480dfac4-232b-4f29-96f6-f28e5e16adc5"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#d6a5a5ae-53c8-49be-8846-d274b30df388"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#38467f4c-fe9b-4f0c-bbd0-374f90c2669b">
    <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:#6aa5ef66-5683-4358-87cc-ab829bc372d8">
    <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:#26277ed1-a67f-4e6b-a690-0ad2d7fcf844">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2006-01-01</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$IGITr1">
    <ns7:boundIntervals xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" rdf:resource="rdf:#$JGITr1"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c4db9901-6076-4c8a-8942-c404735abbcc">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2007-06-15T12:03:12+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bb6a26ce-0c58-4738-8477-9dd978132cac">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">1998-03-01</dcterms:W3CDTF>
  </rdf:Description>
</rdf:RDF>
</model>