Location: Roux, Noble, Noble, Marhl, 2006 @ f3105f9ab946 / roux_noble_noble_marhl_2006.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2008-03-03 22:44:45+13:00
Desc:
committing version01 of roux_noble_noble_marhl_2006
Permanent Source URI:
https://staging.physiomeproject.org/workspace/roux_noble_noble_marhl_2006/rawfile/f3105f9ab9465303708a7e4910457ba8cde05295/roux_noble_noble_marhl_2006.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  roux_model_2006.xml

CREATED :  29th February 2008

AUTHOR :  Catherine Lloyd
          The Bioengineering Institute
          The University of Auckland

DESCRIPTION :  This file contains a CellML description of Roux et al.'s 2006 mathematical model of calcium handling in airway myocytes.

--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="roux_model_2006" cmeta:id="roux_model_2006">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
        This CellML model can be opened in both COR and PCEnv however the model quickly crashes.  The units have been checked and are consistent.  The original model author is working with us on the process of curation and we soon hope to have the CellML model running to replicate the published results.
          </para>
  </section>
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
Abstract:  Airway myocytes are the primary effectors of airway reactivity which modulates airway resistance and hence ventilation. Stimulation of airway myocytes results in an increase in the cytosolic Ca2+ concentration ([Ca2+]i) and the subsequent activation of the contractile apparatus. Many contractile agonists, including acetylcholine, induce [Ca2+]i increase via Ca2+ release from the sarcoplasmic reticulum through InsP3 receptors. Several models have been developed to explain the characteristics of InsP3-induced [Ca2+]i responses, in particular Ca2+ oscillations. The article reviews the modelling of the major structures implicated in intracellular Ca2+ handling, i.e., InsP3 receptors, SERCAs, mitochondria and Ca2+-binding cytosolic proteins. We developed theoretical models specifically dedicated to the airway myocyte which include the major mechanisms responsible for intracellular Ca2+ handling identified in these cells. These biocomputations pointed out the importance of the relative proportion of InsP3 receptor isoforms and the respective role of the different mechanisms responsible for cytosolic Ca2+ clearance in the pattern of [Ca2+]i variations. We have developed a theoretical model of membrane conductances that predicts the variations in membrane potential and extracellular Ca2+ influx. Stimulation of this model by simulated increase in [Ca2+]i predicts membrane depolarisation, but not great enough to trigger a significant opening of voltage-dependant Ca2+ channels. This may explain why airway contraction induced by cholinergic stimulation does not greatly depend on extracellular calcium. The development of such models of airway myocytes is important for the understanding of the cellular mechanisms of airway reactivity and their possible modulation by pharmacological agents.
</para>

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

<para>
<ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&amp;_udi=B6TBN-4GBWFDJ-2&amp;_user=140507&amp;_rdoc=1&amp;_fmt=&amp;_orig=search&amp;_sort=d&amp;view=c&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=fb58937e9ac3f82ee7e81e9b28c139de">Modelling of calcium handling in airway myocytes</ulink>, Etienne Roux, Penelope J. Noble, Denis Noble and Marko Marhl, 2006, <ulink url="http://www.sciencedirect.com/science/journal/00796107">
            <emphasis>Progress in Biophysics and Molecular Biology</emphasis>
          </ulink>, 15, 556-565.  (<ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&amp;_udi=B6TBN-4GBWFDJ-2&amp;_user=140507&amp;_rdoc=1&amp;_fmt=&amp;_orig=search&amp;_sort=d&amp;view=c&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=fb58937e9ac3f82ee7e81e9b28c139de#">Full text</ulink> and PDF versions of the article are available to journal subscribers on the <emphasis>Progress in Biophysics and Molecular Biology</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=15982722&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 15982722</ulink>
</para>

<informalfigure float="0" id="cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram</title>
    </objectinfo>
    <imagedata fileref="roux_2006.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the cell model. The model includes L-type voltage-operated Ca2+ current (ICaL), Ca2+-activated Cl- current (IClCa), Ca2+-activated K+ current (IKCa), delayed rectifier K+ current (IKdr), non-specific cationic current (Icationic), basal currents for Ca2+, Na+ and K+, the plasma membrane Ca2+ ATPase (PMCA) and the Na+-K+ATPase (NKA). Extracellular and intracellular concentrations of Na+, K+, Ca2+ and Cl- are fixed or forced.</caption>
</informalfigure>

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



  
  
  <units name="per_second">
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="per_millimolar_per_second">
    <unit units="millimolar" exponent="-1"/>
    <unit units="second" exponent="-1"/>
  </units>

  <units name="millivolt">
    <unit units="volt" prefix="milli"/>
  </units>

  <units name="per_millivolt">
    <unit units="millivolt" exponent="-1"/>
  </units>
  
  <units name="per_millivolt_per_second">
    <unit units="millivolt" exponent="-1"/>
    <unit units="second" exponent="-1"/>
  </units>

  <units name="per_millivolt_second">
    <unit units="millivolt" exponent="-1"/>
    <unit units="second"/>
  </units>
  
  <units name="microS">
    <unit units="siemens" prefix="micro"/>
  </units>
  
  <units name="microF">
    <unit units="farad" prefix="micro"/>
  </units>

  <units name="picoA">
    <unit units="ampere" prefix="pico"/>
  </units>

  <units name="nanoA">
    <unit units="ampere" prefix="nano"/>
  </units>
  
  <units name="nanoA_per_millimolar">
    <unit units="nanoA"/>
    <unit units="millimolar" exponent="-1"/>
  </units>

  <units name="nanoA_per_millimolar_mole_per_C">
    <unit units="nanoA"/>
    <unit units="millimolar" exponent="-1"/>
    <unit units="mole"/>
    <unit units="coulomb" exponent="-1"/>
  </units>
  
  <units name="millimolar">
    <unit units="mole" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="nanomolar">
    <unit units="mole" prefix="nano"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="millijoule_per_mole_kelvin">
    <unit units="joule" prefix="milli"/>
    <unit units="mole" exponent="-1"/>
    <unit units="kelvin" exponent="-1"/>
  </units>
  
  <units name="coulomb_per_mole">
    <unit units="coulomb"/>
    <unit units="mole" exponent="-1"/>
  </units>
 
  
  
  <component name="environment">
    <variable units="second" public_interface="out" name="time"/>
  </component>
  
  <component name="reversal_potentials">
    <variable units="millivolt" public_interface="out" name="E_Na"/>
    <variable units="millivolt" public_interface="out" name="E_Ca"/>
    <variable units="millivolt" public_interface="out" name="E_K"/>
    <variable units="millivolt" public_interface="out" name="E_Cl"/>
    
    <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    <variable units="second" public_interface="in" name="time"/>
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    <variable units="millimolar" public_interface="in" name="Ki"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="millimolar" public_interface="in" name="Cli"/>
    <variable units="millimolar" public_interface="in" name="Clo"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="E_Na_calculation">
        <eq/>
        <ci> E_Na </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>                
            <apply>
              <divide/>
              <ci> Nao </ci>
              <ci> Nai </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="E_Ca_calculation">
        <eq/>
        <ci> E_Ca </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 2.0 </cn>
              <ci> F </ci>
            </apply>
          </apply>
          <apply>
            <ln/>                
            <apply>
              <divide/>
              <ci> Cao </ci>
              <ci> Cai </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="E_K_calculation">
        <eq/>
        <ci> E_K </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>                
            <apply>
              <divide/>
              <ci> Ko </ci>
              <ci> Ki </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="E_Cl_calculation">
        <eq/>
        <ci> E_Cl </ci>
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
            <ci> F </ci>
          </apply>
          <apply>
            <ln/>                
            <apply>
              <divide/>
              <ci> Cli </ci>
              <ci> Clo </ci>
            </apply>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>

  <component name="delayed_rectifier_K_current"> 
    <variable units="nanoA" public_interface="out" name="i_Kdr"/>          
    
    <variable units="microS" name="g_Kdr" initial_value="0.035"/>
    <variable units="dimensionless" name="gamma_KSS" initial_value="0.15"/>
    
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>   
    <variable units="millivolt" public_interface="in" name="E_K"/>
    
    <variable units="dimensionless" private_interface="in" name="xa"/>
    <variable units="dimensionless" private_interface="in" name="xi1"/>
    <variable units="dimensionless" private_interface="in" name="xi2"/>
      
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Kdr_calculation">
        <eq/>
        <ci> i_Kdr </ci>
        <apply>
          <times/>
          <ci> g_Kdr </ci>
          <apply>
            <plus/>
            <ci> gamma_KSS </ci>
            <apply>
              <times/>
              <apply>
                <plus/>
                <ci> xi1 </ci>
                <ci> xi2 </ci>
              </apply>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> gamma_KSS </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <power/>
            <ci> xa </ci>
            <cn cellml:units="dimensionless"> 2.0 </cn>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="delayed_rectifier_K_current_xa_gate">
    <variable units="dimensionless" public_interface="out" name="xa" initial_value="0.001"/>
   
    <variable units="dimensionless" name="xa_infinity"/>
    <variable units="second" name="tau_xa"/>
    
    <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 id="xa_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> xa </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> xa_infinity </ci>
            <ci> xa </ci>
          </apply>
          <ci> tau_xa </ci>
        </apply>
      </apply>
      
      <apply id="xa_infinity_calculation">
        <eq/>
        <ci> xa_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt"> 5.5 </cn>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="millivolt"> 6.0 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply id="tau_xa_calculation">
        <eq/>
        <ci> tau_xa </ci>
        <apply>
          <minus/>
          <cn cellml:units="second"> 0.003262 </cn>
          <apply>
            <times/>
            <cn cellml:units="per_millivolt_second"> 0.00003552 </cn>
            <ci> V </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="delayed_rectifier_K_current_xi1_gate">
    <variable units="dimensionless" public_interface="out" name="xi1" initial_value="0.0"/>
   
    <variable units="dimensionless" name="xi1_infinity"/>
    <variable units="second" name="tau_xi1" initial_value="0.25"/>
    
    <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 id="xi1_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> xi1 </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> xi1_infinity </ci>
            <ci> xi1 </ci>
          </apply>
          <ci> tau_xi1 </ci>
        </apply>
      </apply>
      
      <apply id="xi1_infinity_calculation">
        <eq/>
        <ci> xi1_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 4.3 </cn>
                </apply>
                <cn cellml:units="millivolt"> 7.5 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
    
  <component name="delayed_rectifier_K_current_xi2_gate">
    <variable units="dimensionless" public_interface="out" name="xi2" initial_value="0.0"/>
   
    <variable units="dimensionless" name="xi2_infinity"/>
    <variable units="second" name="tau_xi2" initial_value="1.98"/>
    
    <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 id="xi2_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> xi2 </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> xi2_infinity </ci>
            <ci> xi2 </ci>
          </apply>
          <ci> tau_xi2 </ci>
        </apply>
      </apply>
      
      <apply id="xi2_infinity_calculation">
        <eq/>
        <ci> xi2_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 4.3 </cn>
                </apply>
                <cn cellml:units="millivolt"> 7.5 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="L_type_Ca_current">
    <variable units="nanoA" public_interface="out" name="i_CaL"/>
    
    <variable units="nanoA_per_millimolar" name="P_CaL" initial_value="0.0030"/>
    <variable units="millivolt" name="V_surf" initial_value="150.0"/>
    
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    
    <variable units="dimensionless" public_interface="out" private_interface="in" name="d"/>
    <variable units="dimensionless" private_interface="in" name="f"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">    
      <apply id="i_CaL_calculation">
        <eq/>
        <ci> i_CaL </ci>
        <apply>
          <times/>
          <cn cellml:units="dimensionless"> 4.0 </cn>
          <ci> P_CaL </ci>
          <ci> d </ci>
          <ci> f </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> V_surf </ci>
          </apply>
          <apply>
            <divide/>
            <apply>
              <divide/>
              <ci> F </ci>
              <apply>
                <times/>
                <ci> R </ci>
                <ci> T </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <apply>
                        <minus/>
                        <ci> V </ci>
                        <ci> V_surf </ci>
                      </apply>
                    </apply>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> Cai </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <apply>
                        <minus/>
                        <ci> V </ci>
                        <ci> V_surf </ci>
                      </apply>
                    </apply>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> Cao </ci>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <apply>
                      <minus/>
                      <apply>
                        <minus/>
                        <ci> V </ci>
                        <ci> V_surf </ci>
                      </apply>
                    </apply>
                    <cn cellml:units="dimensionless"> 2.0 </cn>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_Ca_current_d_gate">
    <variable units="dimensionless" public_interface="out" name="d" initial_value="0.0"/> 

    <variable units="per_second" name="alpha_d"/>
    <variable units="per_second" name="beta_d"/>
    
    <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 id="alpha_d_calculation">
        <eq/>
        <ci> alpha_d </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millivolt_per_second"> 30.0 </cn>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 18.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="per_millivolt"> 0.25 </cn>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 18.0 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply> 
      
      <apply id="beta_d_calculation">
        <eq/>
        <ci> beta_d </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millivolt_per_second"> 12.0 </cn>
            <apply>
              <plus/>
              <ci> V </ci>
              <cn cellml:units="millivolt"> 18.0 </cn>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <times/>
                <cn cellml:units="per_millivolt"> 0.1 </cn>
                <apply>
                  <plus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 18.0 </cn>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply> 
         
      <apply id="dd_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> d </ci>
        </apply>
        <apply>
          <minus/>
          <ci> alpha_d </ci>
          <apply>
            <times/>
            <apply>
              <plus/>
              <ci> alpha_d </ci>
              <ci> beta_d </ci>
            </apply>
            <ci> d </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_Ca_current_f_gate">
    <variable units="dimensionless" public_interface="out" name="f" initial_value="1.0"/> 

    <variable units="second" name="tau_f" initial_value="0.0173"/>
    <variable units="dimensionless" name="f_infinity"/>
    
    <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 id="f_infinity_calculation">
        <eq/>
        <ci> f_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <plus/>
                  <cn cellml:units="millivolt"> 23.0 </cn>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="millivolt"> 6.6 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="df_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> f </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> f_infinity </ci>
            <ci> f </ci>
          </apply>
          <ci> tau_f </ci>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="Ca_activated_K_current"> 
    <variable units="nanoA" public_interface="out" name="i_KCa"/>          
    
    <variable units="microS" name="g_KCa" initial_value="2.45"/>
    
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/> 
    <variable units="millimolar" public_interface="in" private_interface="out" name="Cai"/>     
    <variable units="millivolt" public_interface="in" name="E_K"/>
    <variable units="millijoule_per_mole_kelvin" public_interface="in" private_interface="out" name="R"/>
    <variable units="kelvin" public_interface="in" private_interface="out" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" private_interface="out" name="F"/>
    
    <variable units="dimensionless" private_interface="in" name="xCa1"/>
    <variable units="dimensionless" private_interface="in" name="B"/>
      
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_KCa_calculation">
        <eq/>
        <ci> i_KCa </ci>
        <apply>
          <times/>
          <ci> g_KCa </ci>
          <ci> xCa1 </ci>
          <ci> B </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>
  
  <component name="Ca_activated_K_current_xCa1_gate">
    <variable units="dimensionless" public_interface="out" name="xCa1" initial_value="1.0"/>
   
    <variable units="dimensionless" name="xCa1_infinity"/>
    <variable units="second" name="tau_xCa1" initial_value="1.0"/>
    <variable units="millimolar" name="K2"/>
    <variable units="millimolar" name="K4"/>
    <variable units="per_second" name="alpha" initial_value="280.0"/>
    <variable units="per_second" name="beta" initial_value="480.0"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="second" public_interface="in" name="time"/>
    <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="xCa1_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> xCa1 </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> xCa1_infinity </ci>
            <ci> xCa1 </ci>
          </apply>
          <ci> tau_xCa1 </ci>
        </apply>
      </apply>
      
      <apply id="xCa1_infinity_calculation">
        <eq/>
        <ci> xCa1_infinity </ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <apply>
              <power/>
              <ci> Cai </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <times/>
              <ci> K4 </ci>
              <ci> Cai </ci>
            </apply>
          </apply> 
          <apply>
            <plus/>
            <apply>
              <power/>
              <ci> Cai </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
            <apply>
              <times/>
              <ci> K4 </ci>
              <ci> Cai </ci>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <divide/>
                  <ci> alpha </ci>
                  <ci> beta </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> K4 </ci>
              <ci> K2 </ci>
              <apply>
                <divide/>
                <ci> alpha </ci>
                <ci> beta </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply id="K4_calculation">
        <eq/>
        <ci> K4 </ci>
        <apply>
          <times/>
          <cn cellml:units="millimolar"> 0.0000125 </cn>
          <apply>
            <exp/>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.99 </cn>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> V </ci>
                  <ci> F </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> R </ci>
                  <ci> T </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply id="K2_calculation">
        <eq/>
        <ci> K2 </ci>
        <apply>
          <times/>
          <cn cellml:units="millimolar"> 0.000275 </cn>
          <apply>
            <exp/>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.51 </cn>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <times/>
                  <ci> V </ci>
                  <ci> F </ci>
                </apply>
                <apply>
                  <times/>
                  <ci> R </ci>
                  <ci> T </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="Ca_activated_K_current_B_gate">
    <variable units="dimensionless" public_interface="out" name="B" initial_value="0.001"/> 

    <variable units="per_millimolar_per_second" name="K1"/>
    <variable units="per_second" name="K_1"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="dimensionless" public_interface="in" name="xCa1"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="K1_calculation">
        <eq/>
        <ci> K1 </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millimolar_per_second"> 0.85 </cn>
          <apply>
            <exp/>
            <apply>
              <times/>
              <cn cellml:units="per_millivolt"> 0.04 </cn>
              <ci> V </ci>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply id="K_1_calculation">
        <eq/>
        <ci> K_1 </ci>
        <apply>
          <times/>
          <cn cellml:units="per_second"> 0.24 </cn>
          <apply>
            <exp/>
            <apply>
              <times/>
              <cn cellml:units="per_millivolt"> -0.012 </cn>
              <ci> V </ci>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="dB_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> B </ci>
        </apply>
        <piecewise>
          <piece>
            <apply>
              <minus/>
              <apply>
                <times/>
                <ci> K1 </ci>
                <ci> Cai </ci>
                <ci> xCa1 </ci>
              </apply>
              <apply>
                <times/>
                <ci> K_1 </ci>
                <ci> B </ci>
              </apply>
            </apply>
            <apply>
              <geq/>
              <ci> Cai </ci>
              <cn cellml:units="millimolar"> 0.1 </cn>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="per_second"> 1.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
    </math>
  </component>

  <component name="Ca_activated_Cl_current"> 
    <variable units="nanoA" public_interface="out" name="i_ClCa"/>          
    
    <variable units="microS" name="g_Cl" initial_value="0.01"/>
    <variable units="millimolar" name="CaCT" initial_value="0.0005"/>
    <variable units="dimensionless" name="h" initial_value="3"/>
    
    <variable units="millivolt" public_interface="in" name="V"/> 
    <variable units="millimolar" public_interface="in" name="Cai"/>     
    <variable units="millivolt" public_interface="in" name="E_Cl"/>
      
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_ClCa_calculation">
        <eq/>
        <ci> i_ClCa </ci>
        <apply>
          <times/>
          <ci> g_Cl </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Cl </ci>
          </apply>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <power/>
                <apply>
                  <divide/>
                  <ci> CaCT </ci>
                  <ci> Cai </ci>
                </apply>
                <ci> h </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>

  <component name="nonspecific_cationic_current">
    <variable units="nanoA" public_interface="out" name="i_cationic"/>
    
	<variable units="nanoA" name="i_nsNa"/>
	<variable units="nanoA" name="i_nsK"/>
	<variable units="nanoA" name="i_nsCa"/>
    <variable units="nanoA_per_millimolar_mole_per_C" name="P_nsCa" initial_value="0.0000000175"/>
    <variable units="dimensionless" name="alpha_Ca"/>
    <variable units="millimolar" name="Km_nsCa" initial_value="0.0012"/>
    
    <variable units="second" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="millimolar" public_interface="in" name="Cai"/>
    <variable units="millimolar" public_interface="in" name="Cao"/>
    <variable units="millimolar" public_interface="in" private_interface="out" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Nao"/>
    <variable units="millimolar" public_interface="in" private_interface="out" name="Ki"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    <variable units="millijoule_per_mole_kelvin" public_interface="in" name="R"/>
    <variable units="kelvin" public_interface="in" name="T"/>
    <variable units="coulomb_per_mole" public_interface="in" name="F"/>
    
    <variable units="dimensionless" private_interface="in" name="ns"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">    
      <apply id="i_cationic_calculation">
        <eq/>
        <ci> i_cationic </ci>
        <apply>
          <plus/>
          <ci> i_nsNa </ci>
          <ci> i_nsK </ci>
          <ci> i_nsCa </ci>
        </apply>
      </apply>

      <apply id="i_nsNa_calculation">
        <eq/>
        <ci> i_nsNa </ci>
        <apply>
          <times/>
          <ci> ns </ci>
          <ci> alpha_Ca </ci>
          <ci> P_nsCa </ci>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> V </ci>
              <apply>
                <power/>
                <ci> F </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.75 </cn>
                <ci> Nai </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.75 </cn>
                <ci> Nao </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="i_nsK_calculation">
        <eq/>
        <ci> i_nsK </ci>
        <apply>
          <times/>
          <ci> ns </ci>
          <ci> alpha_Ca </ci>
          <ci> P_nsCa </ci>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> V </ci>
              <apply>
                <power/>
                <ci> F </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.75 </cn>
                <ci> Ki </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 0.75 </cn>
                <ci> Ko </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="i_nsCa_calculation">
        <eq/>
        <ci> i_nsCa </ci>
        <apply>
          <times/>
          <ci> ns </ci>
          <ci> alpha_Ca </ci>
          <ci> P_nsCa </ci>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> V </ci>
              <apply>
                <power/>
                <ci> F </ci>
                <cn cellml:units="dimensionless"> 2.0 </cn>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> R </ci>
              <ci> T </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1.5 </cn>
                <ci> Cai </ci>
                <apply>
                  <exp/>
                  <apply>
                    <divide/>
                    <apply>
                      <times/>
                      <ci> V </ci>
                      <ci> F </ci>
                    </apply>
                    <apply>
                      <times/>
                      <ci> R </ci>
                      <ci> T </ci>
                    </apply>
                  </apply>
                </apply>
              </apply>
              <apply>
                <times/>
                <cn cellml:units="dimensionless"> 1.5 </cn>
                <ci> Cao </ci>
              </apply>
            </apply>
            <apply>
              <minus/>
              <apply>
                <exp/>
                <apply>
                  <divide/>
                  <apply>
                    <times/>
                    <ci> V </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> R </ci>
                    <ci> T </ci>
                  </apply>
                </apply>
              </apply>
              <cn cellml:units="dimensionless"> 1.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
       
      <apply id="alpha_Ca_calculation">
        <eq/>
        <ci> alpha_Ca </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci> Km_nsCa </ci>
                <ci> Cai </ci>
              </apply>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="nonspecific_cationic_current_ns_gate">
    <variable units="dimensionless" public_interface="out" name="ns" initial_value="1.0"/> 

    <variable units="second" name="tau_ns"/>
    <variable units="dimensionless" name="ns_infinity"/>
    
    <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 id="ns_infinity_calculation">
        <eq/>
        <ci> ns_infinity </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <ci> V </ci>
                  <cn cellml:units="millivolt"> 69.8 </cn>
                </apply>
                <cn cellml:units="millivolt"> -11.9 </cn>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="dns_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> ns </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <ci> ns_infinity </ci>
            <ci> ns </ci>
          </apply>
          <ci> tau_ns </ci>
        </apply>
      </apply>

      <apply id="tau_ns_calculation">
        <eq/>
        <ci> tau_ns </ci>
        <apply>
          <minus/>
          <cn cellml:units="second"> 0.20025 </cn>
          <apply>
            <times/>
            <cn cellml:units="per_millivolt_second"> 0.00089875 </cn>
            <ci> V </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="Ca_pump_current"> 
    <variable units="nanoA" public_interface="out" name="i_pCa"/>          
    
    <variable units="nanoA" name="i_pCamax" initial_value="1.15"/>
    <variable units="millimolar" name="Kmp_Ca" initial_value="0.05"/>
    
    <variable units="millimolar" public_interface="in" name="Cai"/>
      
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_pCa_calculation">
        <eq/>
        <ci> i_pCa </ci>
        <apply>
          <times/>
          <ci> i_pCamax </ci>
          <apply>
            <divide/>
            <ci> Cai </ci>
            <apply>
              <plus/>
              <ci> Kmp_Ca </ci>
              <ci> Cai </ci>
            </apply>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>

  <component name="Na_K_pump_current">
    <variable units="nanoA" public_interface="out" name="i_NaK"/>
    
    <variable units="nanoA" name="i_NaKmax" initial_value="0.7"/>
    <variable units="millimolar" name="Km_Na" initial_value="40.0"/>
    <variable units="millimolar" name="Km_K" initial_value="1.0"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millimolar" public_interface="in" name="Nai"/>
    <variable units="millimolar" public_interface="in" name="Ko"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_NaK_calculation">
        <eq/>
        <ci> i_NaK </ci>
        <apply>
          <times/>
          <ci> i_NaKmax </ci>
          <apply>
            <divide/>
            <ci> Ko </ci>
            <apply>
              <plus/>
              <ci> Ko </ci>
              <ci> Km_K </ci>
            </apply>
          </apply>
          <apply>
            <divide/>
            <ci> Nai </ci>
            <apply>
              <plus/>
              <ci> Nai </ci>
              <ci> Km_Na </ci>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="Ca_background_current"> 
    <variable units="nanoA" public_interface="out" name="i_bCa"/>          
    <variable units="microS" name="g_bCa" initial_value="1.194E-5"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>   
    <variable units="millivolt" public_interface="in" name="E_Ca"/>         
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_bCa_calculation">
        <eq/>
        <ci> i_bCa </ci>
        <apply>
          <times/>
          <ci> g_bCa </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Ca </ci>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>

  <component name="K_background_current"> 
    <variable units="nanoA" public_interface="out" name="i_bK"/>          
    <variable units="microS" name="g_bK" initial_value="0.008729"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>   
    <variable units="millivolt" public_interface="in" name="E_K"/>         
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_bK_calculation">
        <eq/>
        <ci> i_bK </ci>
        <apply>
          <times/>
          <ci> g_bK </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>

  <component name="Na_background_current"> 
    <variable units="nanoA" public_interface="out" name="i_bNa"/>          
    <variable units="microS" name="g_bNa" initial_value="0.003263"/>
    
    <variable units="millivolt" public_interface="in" name="V"/>   
    <variable units="millivolt" public_interface="in" name="E_Na"/>         
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_bNa_calculation">
        <eq/>
        <ci> i_bNa </ci>
        <apply>
          <times/>
          <ci> g_bNa </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Na </ci>
          </apply>
        </apply>
      </apply>  
    </math>
  </component>
 
  <component name="membrane">
    <variable units="millivolt" public_interface="out" name="V" initial_value="-60.0"/>
    
    <variable units="microF" name="cap" initial_value="0.00002"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="nanoA" public_interface="in" name="i_CaL"/>
    <variable units="nanoA" public_interface="in" name="i_Kdr"/>
    <variable units="nanoA" public_interface="in" name="i_KCa"/>
    <variable units="nanoA" public_interface="in" name="i_ClCa"/>
    <variable units="nanoA" public_interface="in" name="i_cationic"/>
    <variable units="nanoA" public_interface="in" name="i_NaK"/>
    <variable units="nanoA" public_interface="in" name="i_pCa"/>
    <variable units="nanoA" public_interface="in" name="i_bCa"/>
    <variable units="nanoA" public_interface="in" name="i_bK"/>
    <variable units="nanoA" public_interface="in" name="i_bNa"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="membrane_voltage_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> V </ci>
        </apply>
        <apply>
          <times/>
          <apply>
            <divide/>
            <cn cellml:units="dimensionless"> -1.0 </cn>
            <ci> cap </ci>
          </apply>
          <apply>
            <plus/>
            <ci> i_CaL </ci>
            <ci> i_Kdr </ci>
            <ci> i_KCa </ci>
            <ci> i_ClCa </ci>
            <ci> i_cationic </ci>
            <ci> i_NaK </ci>
            <ci> i_pCa </ci>
            <ci> i_bCa </ci>
            <ci> i_bK </ci>
            <ci> i_bNa </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="model_parameters">
    <variable units="millijoule_per_mole_kelvin" public_interface="out" name="R" initial_value="834.472"/>
    <variable units="kelvin" public_interface="out" name="T" initial_value="310.0"/>
    <variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96485.3415"/>
    <variable units="millimolar" public_interface="out" name="Nai" initial_value="12.0"/>
    <variable units="millimolar" public_interface="out" name="Nao" initial_value="130.0"/>
    <variable units="millimolar" public_interface="out" name="Cai" initial_value="0.1"/>
    <variable units="millimolar" public_interface="out" name="Cao" initial_value="2.0"/>
    <variable units="millimolar" public_interface="out" name="Ki" initial_value="150.0"/>
    <variable units="millimolar" public_interface="out" name="Ko" initial_value="5.0"/>
    <variable units="millimolar" public_interface="out" name="Cli" initial_value="55.0"/>
    <variable units="millimolar" public_interface="out" name="Clo" initial_value="140.0"/>
  </component>
  
  
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="L_type_Ca_current">
        <component_ref component="L_type_Ca_current_d_gate"/>
        <component_ref component="L_type_Ca_current_f_gate"/>
      </component_ref>
      <component_ref component="delayed_rectifier_K_current">
        <component_ref component="delayed_rectifier_K_current_xa_gate"/>
        <component_ref component="delayed_rectifier_K_current_xi1_gate"/>
        <component_ref component="delayed_rectifier_K_current_xi2_gate"/>
      </component_ref>
      <component_ref component="Ca_activated_K_current">
        <component_ref component="Ca_activated_K_current_xCa1_gate"/>
        <component_ref component="Ca_activated_K_current_B_gate"/>
      </component_ref>
      <component_ref component="nonspecific_cationic_current">
        <component_ref component="nonspecific_cationic_current_ns_gate"/>
      </component_ref>
      <component_ref component="Ca_activated_Cl_current"/>
      <component_ref component="Na_K_pump_current"/>
      <component_ref component="Ca_pump_current"/>
      <component_ref component="Na_background_current"/>
      <component_ref component="Ca_background_current"/>
      <component_ref component="K_background_current"/>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="encapsulation"/>
    <component_ref component="L_type_Ca_current">
      <component_ref component="L_type_Ca_current_d_gate"/>
      <component_ref component="L_type_Ca_current_f_gate"/>
    </component_ref>
    <component_ref component="delayed_rectifier_K_current">
      <component_ref component="delayed_rectifier_K_current_xa_gate"/>
      <component_ref component="delayed_rectifier_K_current_xi1_gate"/>
      <component_ref component="delayed_rectifier_K_current_xi2_gate"/>
    </component_ref>
    <component_ref component="Ca_activated_K_current">
      <component_ref component="Ca_activated_K_current_xCa1_gate"/>
      <component_ref component="Ca_activated_K_current_B_gate"/>
    </component_ref>
      <component_ref component="nonspecific_cationic_current">
      <component_ref component="nonspecific_cationic_current_ns_gate"/>
    </component_ref>
  </group>
  
  
  <connection>
    <map_components component_2="environment" component_1="reversal_potentials"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <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="L_type_Ca_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Ca_activated_K_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="delayed_rectifier_K_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
 
  <connection>
    <map_components component_2="environment" component_1="nonspecific_cationic_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="delayed_rectifier_K_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Kdr" variable_1="i_Kdr"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_CaL" variable_1="i_CaL"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_activated_K_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_KCa" variable_1="i_KCa"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_activated_Cl_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_ClCa" variable_1="i_ClCa"/>
  </connection>
  
  <connection>
    <map_components component_2="nonspecific_cationic_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_cationic" variable_1="i_cationic"/>
  </connection>
  
  <connection>
    <map_components component_2="Na_K_pump_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_NaK" variable_1="i_NaK"/>
  </connection>
  
  <connection>
    <map_components component_2="K_background_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_bK" variable_1="i_bK"/>
  </connection>
  
  <connection>
    <map_components component_2="Na_background_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_bNa" variable_1="i_bNa"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_background_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_bCa" variable_1="i_bCa"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_pump_current" component_1="membrane"/>
    <map_variables variable_2="i_pCa" variable_1="i_pCa"/>
  </connection>
  
  <connection>
    <map_components component_2="delayed_rectifier_K_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>

  <connection>
    <map_components component_2="Ca_activated_K_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="K_background_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
    
  <connection>
    <map_components component_2="Na_background_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Na" variable_1="E_Na"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_background_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_activated_Cl_current" component_1="reversal_potentials"/>
    <map_variables variable_2="E_Cl" variable_1="E_Cl"/>
  </connection>
  
  <connection>
    <map_components component_2="reversal_potentials" component_1="model_parameters"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="Cao" variable_1="Cao"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="Clo" variable_1="Clo"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
    <map_variables variable_2="Cli" variable_1="Cli"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_current" component_1="model_parameters"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="Cao" variable_1="Cao"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_activated_K_current" component_1="model_parameters"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>

  <connection>
    <map_components component_2="Ca_activated_Cl_current" component_1="model_parameters"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="nonspecific_cationic_current" component_1="model_parameters"/>
    <map_variables variable_2="Nao" variable_1="Nao"/>
    <map_variables variable_2="Cao" variable_1="Cao"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="Ki" variable_1="Ki"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="Ca_pump_current" component_1="model_parameters"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>

  <connection>
    <map_components component_2="Na_K_pump_current" component_1="model_parameters"/>
    <map_variables variable_2="Ko" variable_1="Ko"/>
    <map_variables variable_2="Nai" variable_1="Nai"/>
  </connection>
  
  <connection>
    <map_components component_2="delayed_rectifier_K_current_xa_gate" component_1="delayed_rectifier_K_current"/>
    <map_variables variable_2="xa" variable_1="xa"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>

  <connection>
    <map_components component_2="delayed_rectifier_K_current_xi1_gate" component_1="delayed_rectifier_K_current"/>
    <map_variables variable_2="xi1" variable_1="xi1"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
 
  <connection>
    <map_components component_2="delayed_rectifier_K_current_xi2_gate" component_1="delayed_rectifier_K_current"/>
    <map_variables variable_2="xi2" variable_1="xi2"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_current_d_gate" component_1="L_type_Ca_current"/>
    <map_variables variable_2="d" variable_1="d"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_Ca_current_f_gate" component_1="L_type_Ca_current"/>
    <map_variables variable_2="f" variable_1="f"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_activated_K_current_xCa1_gate" component_1="Ca_activated_K_current"/>
    <map_variables variable_2="xCa1" variable_1="xCa1"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="T" variable_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="Ca_activated_K_current_B_gate" component_1="Ca_activated_K_current"/>
    <map_variables variable_2="B" variable_1="B"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="Cai" variable_1="Cai"/>
  </connection>
  
  <connection>
    <map_components component_2="Ca_activated_K_current_xCa1_gate" component_1="Ca_activated_K_current_B_gate"/>
    <map_variables variable_2="xCa1" variable_1="xCa1"/>
  </connection>
  
  <connection>
    <map_components component_2="nonspecific_cationic_current_ns_gate" component_1="nonspecific_cationic_current"/>
    <map_variables variable_2="ns" variable_1="ns"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#90e29205-0196-415c-b442-c23739901f41">
    <rdf:li>airway myocyte</rdf:li>
    <rdf:li>smooth muscle</rdf:li>
    <rdf:li>calcium</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#c61e82c6-0348-4747-b08b-f5817e29ca18">
    <rdf:li rdf:resource="rdf:#bf783065-6f8f-4fe2-8590-8cfce358463c"/>
    <rdf:li rdf:resource="rdf:#2cf2a889-08da-4ca5-94c1-9f1e81500ed9"/>
    <rdf:li rdf:resource="rdf:#a72436db-66c4-4141-86d7-4e4d4a72fa91"/>
    <rdf:li rdf:resource="rdf:#f2af3133-665e-4d5c-84b7-704478e4ed59"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#f2af3133-665e-4d5c-84b7-704478e4ed59">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#8ed48a09-82a2-4a6b-85c5-867bb4def2f9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2cf2a889-08da-4ca5-94c1-9f1e81500ed9">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#91709f85-997f-4711-80ee-b2bf8f9af12f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3f7323d2-71d8-4dda-a731-04b88752a91a">
    <vCard:Given>D</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, The Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#19e6c756-e3ce-4713-be5b-cce1809d9fa1"/>
    <dcterms:created rdf:resource="rdf:#53ce3854-e7e7-4e58-b660-db04e045688d"/>
    <dc:creator rdf:resource="rdf:#e3c11c76-0d44-4641-8ae9-55bd5a445e7a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bf783065-6f8f-4fe2-8590-8cfce358463c">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#e6ff530d-dc3b-47f7-8076-7e07afdbdf09"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#53ce3854-e7e7-4e58-b660-db04e045688d">
    <dcterms:W3CDTF>2008-02-29T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#310194c4-a4f6-4271-a0ad-f81b0ecd29d1">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8ed48a09-82a2-4a6b-85c5-867bb4def2f9">
    <vCard:Given>Marko</vCard:Given>
    <vCard:Family>Marhl</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6e8755e9-c009-47db-b21b-f2d2274cdd32">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3a836de6-629a-4b08-afe4-404d9abef616">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4f7b087a-26cb-4bb4-8fc9-7028dc53ddc9">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#90e29205-0196-415c-b442-c23739901f41"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d740a9be-9f95-45ad-a4e0-49b77a606b4d">
    <dc:subject rdf:resource="rdf:#4f7b087a-26cb-4bb4-8fc9-7028dc53ddc9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e6ff530d-dc3b-47f7-8076-7e07afdbdf09">
    <vCard:Given>Etienne</vCard:Given>
    <vCard:Family>Roux</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3f9dbf34-815d-4984-9d3a-e9f9f31ddc2e">
    <dcterms:W3CDTF>2006</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5fa49291-829a-480d-8c25-53c2c9da49a8">
    <bqs:Pubmed_id>15982722</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#a0698944-5d2e-43ce-8dfb-6d5891a3b9d3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e3c11c76-0d44-4641-8ae9-55bd5a445e7a">
    <vCard:ORG rdf:resource="rdf:#6e8755e9-c009-47db-b21b-f2d2274cdd32"/>
    <vCard:EMAIL rdf:resource="rdf:#3a836de6-629a-4b08-afe4-404d9abef616"/>
    <vCard:N rdf:resource="rdf:#a3dc09c3-e63a-4bc0-99ee-98406de45ea5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a3dc09c3-e63a-4bc0-99ee-98406de45ea5">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#92803c4a-2862-4ca8-b743-a9ed2d312fcb">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a72436db-66c4-4141-86d7-4e4d4a72fa91">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#3f7323d2-71d8-4dda-a731-04b88752a91a"/>
  </rdf:Description>
  <rdf:Description rdf:about="#roux_model_2006">
    <dc:title>
        Roux et al.'s 2006 mathematical model of calcium handling in airway myocytes.
      </dc:title>
    <cmeta:bio_entity>airway myocyte</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#b4b9c3c9-8f26-4534-a7b4-14733bc17303"/>
    <bqs:reference rdf:resource="rdf:#5fa49291-829a-480d-8c25-53c2c9da49a8"/>
    <bqs:reference rdf:resource="rdf:#d740a9be-9f95-45ad-a4e0-49b77a606b4d"/>
    <cmeta:species>Homo sapiens</cmeta:species>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b4b9c3c9-8f26-4534-a7b4-14733bc17303">
    <dc:creator rdf:resource="rdf:#310194c4-a4f6-4271-a0ad-f81b0ecd29d1"/>
    <rdf:value>This is the CellML description of Roux et al.'s 2006 mathematical model of calcium handling in airway myocytes.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#19e6c756-e3ce-4713-be5b-cce1809d9fa1">
    <dc:creator rdf:resource="rdf:#92803c4a-2862-4ca8-b743-a9ed2d312fcb"/>
    <rdf:value>This CellML model can be opened in both COR and PCEnv however the model quickly crashes.  The units have been checked and are consistent.  The original model author is working with us on the process of curation and we soon hope to have the CellML model running to replicate the published results.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2c0db3ea-b819-4bc0-bd76-b6bbd67160c3">
    <dc:title>Progress in Biophysics and Molecular Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#91709f85-997f-4711-80ee-b2bf8f9af12f">
    <vCard:Given>P</vCard:Given>
    <vCard:Family>Noble</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a0698944-5d2e-43ce-8dfb-6d5891a3b9d3">
    <dc:creator rdf:resource="rdf:#c61e82c6-0348-4747-b08b-f5817e29ca18"/>
    <dc:title>Modelling of calcium handling in airway myocytes</dc:title>
    <bqs:volume>90</bqs:volume>
    <bqs:first_page>64</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#2c0db3ea-b819-4bc0-bd76-b6bbd67160c3"/>
    <dcterms:issued rdf:resource="rdf:#3f9dbf34-815d-4984-9d3a-e9f9f31ddc2e"/>
    <bqs:last_page>87</bqs:last_page>
  </rdf:Description>
</rdf:RDF>
</model>