Rendering of the source text

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 17/07/2007 at 16:19:18 using:

COR (0.9.31.723)
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#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:ns9="http://www.cellml.org/metadata/simulation/1.0#" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="endresen_1997" name="endresen_1997">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Mathematical modeling and qualitative analysis of insulin therapies</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
        This CellML version of the model has been checked in COR and OpenCell and is known to be valid CellML. The units are consistent and the model runs to recreate the published results.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: A model of the rabbit sinoatrial action potential is introduced, based on a model by Morris and Lecar. One cell is described by two nonlinear first-order ordinary differential equations, with ten constant parameters. The model is much simpler than most other models in use, but can reproduce perfectly experimentally recorded action potentials. The dynamics of two coupled cells, with and without the presence of periodic acetylcholine pulses, shows examples of bifurcations and strange attractors, mathematical phenomena characterizing chaotic motion. It remains to be clarified whether such dynamics is actually observed, for example in the small irregular variations of the normal heart rate.
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="endresen_1996.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the cell model</caption>
</informalfigure>

<para>
The original paper reference is cited below:
</para>

<para>
Chaos in weakly-coupled pacemaker cells, Lars Petter Endresen, 1997, <emphasis>Journal of Theoretical Biology</emphasis>, 184, 41-50. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=9039399&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 9039399</ulink>
</para>

</sect1>
</article>
</documentation> 


   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="per_millivolt">
      <unit units="volt" prefix="milli" exponent="-1"/>
   </units>
   <units name="picoS">
      <unit units="siemens" prefix="pico"/>
   </units>
   <units name="picoF">
      <unit units="farad" prefix="pico"/>
   </units>
   <units name="femtoA">
      <unit units="ampere" prefix="femto"/>
   </units>
   <units name="molar">
      <unit units="mole"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <component name="environment">
      <variable units="second" public_interface="out" name="time"/>
   </component>
   <component name="membrane">
      <variable units="millivolt" public_interface="out" name="V" initial_value="-52.07606"/>
      <variable units="picoF" name="Cm" initial_value="60"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="femtoA" public_interface="in" name="i_s"/>
      <variable units="femtoA" public_interface="in" name="i_K"/>
      <variable units="femtoA" public_interface="in" name="i_K_ACh"/>
      <variable units="femtoA" public_interface="in" name="i_j"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_s</ci>
                     <ci>i_K</ci>
                     <ci>i_K_ACh</ci>
                     <ci>i_j</ci>
                  </apply>
               </apply>
               <ci>Cm</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="calcium_channel">
      <variable units="femtoA" public_interface="out" name="i_s"/>
      <variable units="picoS" name="g_s" initial_value="382.9118"/>
      <variable units="millivolt" name="V_s" initial_value="214.1429"/>
      <variable units="millivolt" name="V_1" initial_value="-35.9358"/>
      <variable units="millivolt" name="V_2" initial_value="7.8589"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_s</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <cn cellml:units="dimensionless">1</cn>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <ci>g_s</ci>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <tanh/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <ci>V_1</ci>
                        </apply>
                        <ci>V_2</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_s</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="potassium_channel">
      <variable units="femtoA" public_interface="out" name="i_K"/>
      <variable units="picoS" name="g_K" initial_value="536.1093"/>
      <variable units="millivolt" name="V_K" initial_value="-259.0783"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="w"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K</ci>
            <apply>
               <times/>
               <ci>g_K</ci>
               <ci>w</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="potassium_channel_w_gate">
      <variable units="dimensionless" public_interface="out" name="w" initial_value="0.0008971"/>
      <variable units="per_second" name="lambda_w" initial_value="20.7796"/>
      <variable units="millivolt" name="V_3" initial_value="-27.9375"/>
      <variable units="millivolt" name="V_4" initial_value="6.321"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>w</ci>
            </apply>
            <apply>
               <times/>
               <ci>lambda_w</ci>
               <apply>
                  <cosh/>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <ci>V</ci>
                        <ci>V_3</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>V_4</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <tanh/>
                           <apply>
                              <divide/>
                              <apply>
                                 <minus/>
                                 <ci>V</ci>
                                 <ci>V_3</ci>
                              </apply>
                              <ci>V_4</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <ci>w</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="acetyl_choline_activated_potassium_channel">
      <variable units="femtoA" public_interface="out" name="i_K_ACh"/>
      <variable units="second" public_interface="in" private_interface="out" name="time"/>
      <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
      <variable units="dimensionless" private_interface="in" name="u"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K_ACh</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">1</cn>
               <cn cellml:units="picoS">0.27</cn>
               <ci>u</ci>
               <apply>
                  <plus/>
                  <ci>V</ci>
                  <cn cellml:units="millivolt">90</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="acetyl_choline_activated_potassium_channel_u_gate">
      <variable units="dimensionless" public_interface="out" name="u" initial_value="0.2344555"/>
      <variable units="per_second" name="alpha"/>
      <variable units="per_second" name="beta"/>
      <variable units="molar" name="ACh" initial_value="1e-6"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="second" public_interface="in" name="time"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>alpha</ci>
            <apply>
               <divide/>
               <cn cellml:units="per_second">0.012332</cn>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <cn cellml:units="molar" type="e-notation">4.2<sep/>-6</cn>
                     <ci>ACh</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>beta</ci>
            <apply>
               <times/>
               <cn cellml:units="per_second">0.01</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <cn cellml:units="per_millivolt">0.0133</cn>
                     <apply>
                        <plus/>
                        <ci>V</ci>
                        <cn cellml:units="millivolt">40</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>u</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>alpha</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>u</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta</ci>
                  <ci>u</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="coupling_current">
      <variable units="femtoA" public_interface="out" name="i_j"/>
      <variable units="picoS" name="g_j" initial_value="0"/>
      <variable units="millivolt" public_interface="in" name="V"/>
      <variable units="millivolt" name="V_B" initial_value="-50"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_j</ci>
            <apply>
               <times/>
               <ci>g_j</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_B</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="calcium_channel"/>
         <component_ref component="potassium_channel">
            <component_ref component="potassium_channel_w_gate"/>
         </component_ref>
         <component_ref component="acetyl_choline_activated_potassium_channel">
            <component_ref component="acetyl_choline_activated_potassium_channel_u_gate"/>
         </component_ref>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="potassium_channel">
         <component_ref component="potassium_channel_w_gate"/>
      </component_ref>
      <component_ref component="acetyl_choline_activated_potassium_channel">
         <component_ref component="acetyl_choline_activated_potassium_channel_u_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="potassium_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="acetyl_choline_activated_potassium_channel"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="calcium_channel" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_s" variable_1="i_s"/>
   </connection>
   <connection>
      <map_components component_2="potassium_channel" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K" variable_1="i_K"/>
   </connection>
   <connection>
      <map_components component_2="acetyl_choline_activated_potassium_channel" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_K_ACh" variable_1="i_K_ACh"/>
   </connection>
   <connection>
      <map_components component_2="potassium_channel_w_gate" component_1="potassium_channel"/>
      <map_variables variable_2="w" variable_1="w"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="acetyl_choline_activated_potassium_channel_u_gate" component_1="acetyl_choline_activated_potassium_channel"/>
      <map_variables variable_2="u" variable_1="u"/>
      <map_variables variable_2="time" variable_1="time"/>
      <map_variables variable_2="V" variable_1="V"/>
   </connection>
   <connection>
      <map_components component_2="coupling_current" component_1="membrane"/>
      <map_variables variable_2="V" variable_1="V"/>
      <map_variables variable_2="i_j" variable_1="i_j"/>
   </connection>

<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#681639cc-a2ef-42ea-a526-bc882f863c8b">
    <rdf:li rdf:resource="rdf:#09699ae2-00bd-46f7-938d-33197d0eeff6"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#$cFMBC2">
    <ns7:boundIntervals rdf:resource="rdf:#$dFMBC2"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#251aeac3-44bd-4bbb-ae3e-23b01226b981">
    <vCard:FN>James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>Oxford University</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#2c80736c-e859-421a-bea6-25ae4cf08209"/>
    <dcterms:created rdf:resource="rdf:#1d03292c-8dbc-4762-b560-4d8d86a43832"/>
    <dc:creator rdf:resource="rdf:#ffda8d00-d5dd-4d3e-bc89-317dc176d460"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#596283fa-cef8-49b7-b36a-a836d28b607e">
    <dcterms:W3CDTF>1996-01-07 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c7b7b9f9-9a5c-4a68-ad5f-965d6005886e">
    <vCard:Orgname>Oxford University</vCard:Orgname>
    <vCard:Orgunit>Dept of Physiology, Anatomy &amp; Genetics</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#eb58e65a-c1ac-4367-96da-1f07380a294f">
    <dc:title>Journal of Theoretical Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1d03292c-8dbc-4762-b560-4d8d86a43832">
    <dcterms:W3CDTF>2007-09-11T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ffda8d00-d5dd-4d3e-bc89-317dc176d460">
    <vCard:ORG rdf:resource="rdf:#c7b7b9f9-9a5c-4a68-ad5f-965d6005886e"/>
    <vCard:EMAIL rdf:resource="rdf:#366ce1a5-ad3d-4ac5-8e0c-c6df1f1961f6"/>
    <vCard:N rdf:resource="rdf:#37220945-17e2-4486-8f5b-aab7e12c651c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#37220945-17e2-4486-8f5b-aab7e12c651c">
    <vCard:Given>Penny</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$eFMBC2">
    <ns7:nonstandard-algorithm>am112</ns7:nonstandard-algorithm>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#50da750b-df24-407d-9c83-0b3dc29c9768">
    <bqs:Pubmed_id>9039399</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#703fc683-0ae4-402e-88e5-b1225749cc90"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#366ce1a5-ad3d-4ac5-8e0c-c6df1f1961f6">
    <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:#809b8eed-69f0-4f13-920c-dea2a96a1898">
    <vCard:FN>James Lawson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#703fc683-0ae4-402e-88e5-b1225749cc90">
    <dc:creator rdf:resource="rdf:#681639cc-a2ef-42ea-a526-bc882f863c8b"/>
    <dc:title>Chaos in Weakly-coupled Pacemaker Cells</dc:title>
    <bqs:volume>184(1)</bqs:volume>
    <bqs:first_page>41</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#eb58e65a-c1ac-4367-96da-1f07380a294f"/>
    <dcterms:issued rdf:resource="rdf:#596283fa-cef8-49b7-b36a-a836d28b607e"/>
    <bqs:last_page>50</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#51aea0e7-d1f8-478c-b157-e3bef4323451">
    <dc:creator rdf:resource="rdf:#809b8eed-69f0-4f13-920c-dea2a96a1898"/>
    <rdf:value>This file is a CellML description of Endresen's 1996 model of the rabbit sinoatrial action potential from the paper &quot;Chaos in Weakly-coupled Pacemaker Cells&quot;</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#endresen_1997">
    <ns7:simulation rdf:resource="rdf:#$cFMBC2"/>
    <bqs:reference rdf:resource="rdf:#50da750b-df24-407d-9c83-0b3dc29c9768"/>
<bqs:reference rdf:parseType="Resource">
  <dc:subject rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>electrophysiology</rdf:li>
        <rdf:li>cardiac</rdf:li>
        <rdf:li>rabbit</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
    <cmeta:comment rdf:resource="rdf:#51aea0e7-d1f8-478c-b157-e3bef4323451"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2c80736c-e859-421a-bea6-25ae4cf08209">
    <dc:creator rdf:resource="rdf:#251aeac3-44bd-4bbb-ae3e-23b01226b981"/>
    <rdf:value>This model is known to run in COR and PCEnv 0.2 and reproduces the published results.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#09699ae2-00bd-46f7-938d-33197d0eeff6">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#bd8bd4f8-4d61-4311-8c10-4364c1b8d375"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bd8bd4f8-4d61-4311-8c10-4364c1b8d375">
    <vCard:Given>Lars</vCard:Given>
    <vCard:Family>Endresen</vCard:Family>
    <vCard:Other>P</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#$dFMBC2">
    <rdf:rest rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
    <rdf:first rdf:resource="rdf:#$eFMBC2"/>
  </rdf:Description>
</rdf:RDF>
</model>