Location: Endresen, 1996 @ 56893a073b3d / endresen_1996.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2009-11-19 14:14:54+13:00
Desc:
Added images in ai svg and png format
Permanent Source URI:
http://models.cellml.org/workspace/endresen_1996/rawfile/56893a073b3dbb67b7f32341f4657d3c63bbb57f/endresen_1996.cellml

<?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>
In the paper described here, Lars Petter Endresen presents a very simple mathematical model of the rabbit sinoatrial action potential.  One cell is defined by two nonlinear first-order ordinary differential equations, using just ten constant parameters.  These two equations represent the slow inward calcium current, and the outward potassium current.  Despite its inherent simplicity, the model can reproduce the complex waveforms and rhythms observed in experimentally recorded action potentials.
</para>

<para>
The complete 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>