Rendering of the source text

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

CREATED :  11th July 2002

LAST MODIFIED : 27th April 2005

AUTHOR :  Catherine Lloyd
          The Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the CellML Metadata 1.0 Specification released on 16th
January, 2002.

DESCRIPTION :  This file contains a CellML description of Miftakhov et al's 1999 mathematical model of the electrical activity of a smooth muscle cell from the small intestine. 

CHANGES: 
  18/07/2002 - CML - Added more metadata
  09/04/2003 - AAC - Added publication date information.  
  27/04/2005 - PJV - Changed initial values to make them valid cellml.
--><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#" cmeta:id="miftakhov_abdusheva_christensen_1999_version01" name="miftakhov_abdusheva_christensen_1999_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Motility Patterns of the Small Bowel</title>
  <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 is the original unchecked version of the model imported from the previous
            CellML model repository, 24-Jan-2006.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
Intestinal motility describes the rhythmic, propulsive contractile activity of the distal stomach, intestine and colon.  Peristaltic waves are generated through interplay between the enteric nervous system, smooth muscle cells and interstitial cells of Cajal (ICC).  ICC are critical for initiating slow wave activity, making ICC the <emphasis>pacemakers</emphasis> of the gut.  The slow wave membrane depolarisations, generated by ICC, determine the maximum frequency and the propagation characteristics of the intestinal rhythmic contractions.
</para>

<para>
Smooth muscle contraction is the fundamental event in gastrointestinal motion.  Smooth muscle cells respond to the slow wave activity imposed by ICC pacemaker cells.  Although many of the biochemical mechanisms underlying the excitation-contraction coupling are not yet defined, it is known that cytosolic Ca<superscript>2+</superscript> is the essential component in the coupling phenomenon.  Slow wave membrane depolarisation activates smooth muscle contraction, mainly by triggering an L-type Ca<superscript>2+</superscript> current.  This is the main source of Ca<superscript>2+</superscript> for contraction.
</para>

<para>
Since the late 1960s, slow wave propagation in the gastrointestinal tract has been the subject of mathematical models.  In 1968 Nelson and Becker proposed that a chain of relaxation oscillators could simulate the electromechanical activity in the small intestine.  For the next two decades this idea persisted.  However, more recently core conductor models, which consider the ionic fluxes and the electrical activity of individual cells, have been developed.   
</para>

<para>
In 1999, Miftakhov <emphasis>et al.</emphasis> published a mathematical model, which aimed to simulate motility patterns of the small bowel.  Within the functional unit of an electromyogenic syncytium are equations that represent the dynamics of ion channels in smooth muscle cells.   These include voltage-gated Ca<superscript>2+</superscript> channels (L and T-type), a voltage-gated K<superscript>+</superscript> channel, a calcium-activated K<superscript>+</superscript> channel, and a leak chloride current.   The kinetics of these ion channels determine the electrical activity of the whole functional unit.
</para>

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

<para>
<ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&amp;_udi=B6WMD-45FSB7Y-53&amp;_coverDate=03%2F07%2F1999&amp;_alid=187124419&amp;_rdoc=1&amp;_fmt=&amp;_orig=search&amp;_qd=1&amp;_cdi=6932&amp;_sort=d&amp;view=c&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=abeef0b8b2b3eb46ad90e8680570122c">Numerical Simulation of Motility Patterns of the Small Bowel. 1. Formulation of a Mathematical Model</ulink>, R. N. Miftakhov, G. R. Abdusheva and J. Christensen, 1999, <ulink url="http://www.sciencedirect.com/science?_ob=JournalURL&amp;_cdi=6932&amp;_auth=y&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=c8e64d20ff203cda2255f4b6b7affe4f">
            <emphasis>Journal of Theoretical Biology</emphasis>
          </ulink>, 197, 89-112. (A PDF version of the article is available to subscribers on the <emphasis>Science Direct</emphasis>.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=10036210&amp;dopt=Abstract">PubMed ID: 10036210</ulink> 
</para>

<para>
The raw CellML descriptions of the model can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>.
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell_diagram</title>
    </objectinfo>
    <imagedata fileref="miftakhov_1999.png"/>
  </imageobject>
</mediaobject>
<caption>The myoelectrical activity of the gastrointestinal smooth muscle cell is governed by the dynamics of voltage-dependent Ca<superscript>2+</superscript> L and T-type channels (<emphasis>I<subscript>Ca,L</subscript> and I<subscript>Ca,T</subscript>
            </emphasis>), a voltage-gated K<superscript>+</superscript> channel (<emphasis>I<subscript>K</subscript>
            </emphasis>), a calcium-activated K<superscript>+</superscript> channel (<emphasis>I<subscript>Ca-K</subscript>
            </emphasis>), and a leak chloride current (<emphasis>I<subscript>Cl</subscript>
            </emphasis>).</caption>
</informalfigure>

</sect1>
</article>
</documentation>
  
  
  <!--
    Below, we define some additional units for association with variables and
    constants within the model. The identifiers are fairly self-explanatory.
  -->
  <units name="millisecond">
    <unit units="second" prefix="milli"/>
  </units>
  
  <units name="per_millisecond">
    <unit units="second" prefix="milli" 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="per_millivolt_millisecond">
    <unit units="millivolt" exponent="-1"/>
    <unit units="millisecond" exponent="-1"/>
  </units>
  
  <units name="milliS_per_cm2">
    <unit units="siemens" prefix="milli"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </units>
  
  <units name="microF_per_cm2">
    <unit units="farad" prefix="micro"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </units>
  
  <units name="millimolar">
    <unit units="mole" prefix="milli"/>
    <unit units="litre" exponent="-1"/>
  </units> 
  
  <units name="microA_per_cm2">
    <unit units="ampere" prefix="micro"/>
    <unit units="metre" prefix="centi" exponent="-2"/>
  </units>
  
  <component name="environment">
    <variable units="millisecond" public_interface="out" name="time"/>
  </component>
  
  <component name="model_constants">
    <variable units="dimensionless" public_interface="out" name="alpha" initial_value="0.12"/>
    <variable units="dimensionless" public_interface="out" name="lamda" initial_value="12.5"/>
  </component>
  
  <component name="membrane">
    <variable units="millivolt" public_interface="out" name="V" initial_value="-55.0"/>
    
    <variable units="microF_per_cm2" name="Cm" initial_value="2.5"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="dimensionless" public_interface="in" name="alpha"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_Ca_T"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_Ca_L"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_Ca_K"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_K"/>
    <variable units="microA_per_cm2" public_interface="in" name="i_Cl"/>
    
    <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>
            <minus/>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <ci> Cm </ci>
                <ci> alpha </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <ci> i_Ca_T </ci>
            <ci> i_Ca_L </ci>
            <ci> i_Ca_K </ci>
            <ci> i_K </ci>
            <ci> i_Cl </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="T_type_calcium_current">
    <variable units="microA_per_cm2" public_interface="out" name="i_Ca_T"/>
    <variable units="millivolt" public_interface="out" name="E_Ca" initial_value="80.0"/>
    
    <variable units="milliS_per_cm2" name="g_Ca_T" initial_value="0.51"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="dimensionless" public_interface="in" private_interface="out" name="alpha"/>
    <variable units="dimensionless" public_interface="in" private_interface="out" name="lamda"/>
    
    <variable units="dimensionless" private_interface="in" name="m"/>
    <variable units="dimensionless" private_interface="in" name="h"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Ca_T_calculation">
        <eq/>
        <ci> i_Ca_T </ci>
        <apply>
          <times/>
          <ci> g_Ca_T </ci>
          <apply>
            <power/>
            <ci> m </ci>
            <cn cellml:units="dimensionless"> 3.0 </cn>
          </apply>
          <ci> h </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="T_type_calcium_current_m_gate">
    <variable units="dimensionless" public_interface="out" name="m"/>
        
    <variable units="per_millisecond" name="alpha_m"/>
    <variable units="per_millisecond" name="beta_m"/>
        
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_m_calculation">
        <eq/>
        <ci> alpha_m </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millivolt_millisecond"> 0.1 </cn>
            <apply>
              <minus/>
              <cn cellml:units="millivolt"> 50.0 </cn>
              <ci> V </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 5.0 </cn>
                <apply>
                  <times/>
                  <ci> V </ci>
                  <cn cellml:units="dimensionless"> 0.1 </cn>
                </apply>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_m_calculation">
        <eq/>
        <ci> beta_m </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millisecond"> 4.0 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <cn cellml:units="millivolt"> 25.0 </cn>
                <ci> V </ci>
              </apply>
              <cn cellml:units="millivolt"> 18.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="dm_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> m </ci>
        </apply>
        <apply>
          <divide/>
          <ci> alpha_m </ci>
          <apply>
            <plus/>
            <ci> alpha_m </ci>
            <ci> beta_m </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="T_type_calcium_current_h_gate">
    <variable units="dimensionless" public_interface="out" name="h"/>
        
    <variable units="per_millisecond" name="alpha_h"/>
    <variable units="per_millisecond" name="beta_h"/>
    
    <variable units="dimensionless" public_interface="in" name="alpha"/>
    <variable units="dimensionless" public_interface="in" name="lamda"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_h_calculation">
        <eq/>
        <ci> alpha_h </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millisecond"> 0.07 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <cn cellml:units="millivolt"> 25.0 </cn>
                <ci> V </ci>
              </apply>
              <cn cellml:units="millivolt"> 20.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_h_calculation">
        <eq/>
        <ci> beta_h </ci>
        <apply>
          <divide/>
          <cn cellml:units="dimensionless"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="millisecond"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 5.5 </cn>
                <apply>
                  <times/>
                  <ci> V </ci>
                  <cn cellml:units="dimensionless"> 0.1 </cn>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="dh_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> h </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> alpha_h </ci>
              <apply>
                <minus/>
                <cn cellml:units="per_millivolt"> 1.0 </cn>
                <ci> h </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> beta_h </ci>
              <ci> h </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> alpha </ci>
            <ci> lamda </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_calcium_current">
    <variable units="microA_per_cm2" public_interface="out" name="i_Ca_L"/>

    <variable units="milliS_per_cm2" name="g_Ca_L" initial_value="0.004"/>
    
    <variable units="millivolt" public_interface="in" name="E_Ca"/>
    <variable units="dimensionless" public_interface="in" private_interface="out" name="alpha"/>
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    
    <variable units="dimensionless" public_interface="out" private_interface="in" name="x_Ca"/>
   
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Ca_L_calculation">
        <eq/>
        <ci> i_Ca_L </ci>
        <apply>
          <times/>
          <ci> g_Ca_L </ci>
          <ci> x_Ca </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="L_type_calcium_current_x_Ca_gate">
    <variable units="dimensionless" public_interface="out" name="x_Ca"/>
        
    <variable units="millisecond" name="tau_x_Ca" initial_value="500.0"/>
   
    <variable units="dimensionless" public_interface="in" name="alpha"/>    
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="dx_Ca_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> x_Ca </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <divide/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <plus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <exp/>
                  <apply>
                    <times/>
                    <cn cellml:units="dimensionless"> 0.15 </cn>
                    <apply>
                      <minus/>
                      <apply>
                        <minus/>
                        <ci> V </ci>
                      </apply>
                      <cn cellml:units="dimensionless"> 50.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
            <ci> x_Ca </ci>
          </apply>
          <apply>
            <times/>
            <ci> alpha </ci>
            <ci> tau_x_Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="potassium_current">
    <variable units="microA_per_cm2" public_interface="out" name="i_K"/>
    <variable units="millivolt" public_interface="out" name="E_K" initial_value="-75.0"/>
    
    <variable units="milliS_per_cm2" name="g_K" initial_value="0.3"/>
    
    <variable units="millisecond" public_interface="in" private_interface="out" name="time"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    <variable units="dimensionless" public_interface="in" private_interface="out" name="alpha"/>
    <variable units="dimensionless" public_interface="in" private_interface="out" name="lamda"/>
    
    <variable units="dimensionless" private_interface="in" name="n"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_K_calculation">
        <eq/>
        <ci> i_K </ci>
        <apply>
          <times/>
          <ci> g_K </ci>
          <apply>
            <power/>
            <ci> n </ci>
            <cn cellml:units="dimensionless"> 4.0 </cn>
          </apply>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_K </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="potassium_current_n_gate">
    <variable units="dimensionless" public_interface="out" name="n"/>
        
    <variable units="per_millisecond" name="alpha_n"/>
    <variable units="per_millisecond" name="beta_n"/>
    
    <variable units="dimensionless" public_interface="in" name="alpha"/>
    <variable units="dimensionless" public_interface="in" name="lamda"/>   
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="alpha_n_calculation">
        <eq/>
        <ci> alpha_n </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <cn cellml:units="per_millivolt_millisecond"> 0.01 </cn>
            <apply>
              <minus/>
              <cn cellml:units="millivolt"> 55.0 </cn>
              <ci> V </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt"> 55.0 </cn>
                  <ci> V </ci>
                </apply>
                <cn cellml:units="millivolt"> 10.0 </cn>
              </apply>
            </apply>
            <cn cellml:units="dimensionless"> 1.0 </cn>
          </apply>
        </apply>
      </apply>
      
      <apply id="beta_n_calculation">
        <eq/>
        <ci> beta_n </ci>
        <apply>
          <times/>
          <cn cellml:units="per_millisecond"> 0.125 </cn>
          <apply>
            <exp/>
            <apply>
              <divide/>
              <apply>
                <minus/>
                <cn cellml:units="millivolt"> 45.0 </cn>
                <ci> V </ci>
              </apply>
              <cn cellml:units="millivolt"> 80.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="dn_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> n </ci>
        </apply>
        <apply>
          <divide/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> alpha_n </ci>
              <apply>
                <minus/>
                <cn cellml:units="per_millivolt"> 1.0 </cn>
                <ci> n </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> beta_n </ci>
              <ci> n </ci>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> alpha </ci>
            <ci> lamda </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="calcium_activated_potassium_current">
    <variable units="microA_per_cm2" public_interface="out" name="i_Ca_K"/>
    
    <variable units="millimolar" name="Ca" initial_value="1E-4"/>
    <variable units="milliS_per_cm2" name="g_Ca_K" initial_value="0.03"/>
    <variable units="per_millisecond" name="rho" initial_value="0.125E3"/>
    <variable units="per_millivolt" name="K_c" initial_value="425.0E-5"/>
  
    <variable units="dimensionless" public_interface="in" name="alpha"/>
    <variable units="millivolt" public_interface="in" name="E_K"/>
    <variable units="millivolt" public_interface="in" name="E_Ca"/>
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    <variable units="dimensionless" public_interface="in" name="x_Ca"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Ca_K_calculation">
        <eq/>
        <ci> i_Ca_K </ci>
        <apply>
          <divide/>
          <apply>
            <times/>
            <ci> g_Ca_K </ci>
            <ci> Ca </ci>
            <apply>
              <minus/>
              <ci> V </ci>
              <ci> E_K </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="millimolar"> 0.5 </cn>
            <ci> Ca </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="dCa_dt">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca </ci>
        </apply>
        <apply>
          <times/>
          <apply>
            <divide/>
            <ci> rho </ci>
            <ci> alpha </ci>
          </apply>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> K_c </ci>
              <ci> x_Ca </ci>
              <apply>
                <minus/>
                <ci> E_Ca </ci>
                <ci> V </ci>
              </apply>
            </apply>
            <ci> Ca </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="leak_chloride_current">
    <variable units="microA_per_cm2" public_interface="out" name="i_Cl"/>

    <variable units="milliS_per_cm2" name="g_Cl" initial_value="0.003"/>
    <variable units="millivolt" name="E_Cl" initial_value="-40.0"/>
    
    <variable units="millisecond" public_interface="in" name="time"/>
    <variable units="millivolt" public_interface="in" name="V"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="i_Cl_calculation">
        <eq/>
        <ci> i_Cl </ci>
        <apply>
          <times/>
          <ci> g_Cl </ci>
          <apply>
            <minus/>
            <ci> V </ci>
            <ci> E_Cl </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <group>
    <relationship_ref relationship="containment"/>
    <component_ref component="membrane">
      <component_ref component="T_type_calcium_current">
        <component_ref component="T_type_calcium_current_m_gate"/>
        <component_ref component="T_type_calcium_current_h_gate"/>
      </component_ref>
      <component_ref component="L_type_calcium_current">
        <component_ref component="L_type_calcium_current_x_Ca_gate"/>
      </component_ref>
      <component_ref component="potassium_current">
        <component_ref component="potassium_current_n_gate"/>
      </component_ref>
      <component_ref component="calcium_activated_potassium_current"/>
      <component_ref component="leak_chloride_current"/>
    </component_ref>
  </group>
  
  <group>
    <relationship_ref relationship="encapsulation"/>
    <component_ref component="T_type_calcium_current">
      <component_ref component="T_type_calcium_current_m_gate"/>
      <component_ref component="T_type_calcium_current_h_gate"/>
    </component_ref>
    <component_ref component="L_type_calcium_current">
      <component_ref component="L_type_calcium_current_x_Ca_gate"/>
    </component_ref>
    <component_ref component="potassium_current">
      <component_ref component="potassium_current_n_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="T_type_calcium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="L_type_calcium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="calcium_activated_potassium_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="leak_chloride_current"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_calcium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Ca_T" variable_1="i_Ca_T"/>
  </connection>
  
  <connection>
    <map_components component_2="model_constants" component_1="membrane"/>
    <map_variables variable_2="alpha" variable_1="alpha"/> 
  </connection>
  
  <connection>
    <map_components component_2="L_type_calcium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Ca_L" variable_1="i_Ca_L"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_current" 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="calcium_activated_potassium_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Ca_K" variable_1="i_Ca_K"/>
  </connection>
  
  <connection>
    <map_components component_2="leak_chloride_current" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="i_Cl" variable_1="i_Cl"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_calcium_current" component_1="L_type_calcium_current"/>
    <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_calcium_current" component_1="calcium_activated_potassium_current"/>
    <map_variables variable_2="E_Ca" variable_1="E_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_activated_potassium_current" component_1="potassium_current"/>
    <map_variables variable_2="E_K" variable_1="E_K"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_activated_potassium_current" component_1="L_type_calcium_current"/>
    <map_variables variable_2="x_Ca" variable_1="x_Ca"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_calcium_current" component_1="model_constants"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
    <map_variables variable_2="lamda" variable_1="lamda"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_calcium_current" component_1="model_constants"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_current" component_1="model_constants"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
    <map_variables variable_2="lamda" variable_1="lamda"/>
  </connection>
  
  <connection>
    <map_components component_2="calcium_activated_potassium_current" component_1="model_constants"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_calcium_current_m_gate" component_1="T_type_calcium_current"/>
    <map_variables variable_2="m" variable_1="m"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>

  <connection>
    <map_components component_2="T_type_calcium_current_h_gate" component_1="T_type_calcium_current"/>
    <map_variables variable_2="h" variable_1="h"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
    <map_variables variable_2="lamda" variable_1="lamda"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="L_type_calcium_current_x_Ca_gate" component_1="L_type_calcium_current"/>
    <map_variables variable_2="x_Ca" variable_1="x_Ca"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_current_n_gate" component_1="potassium_current"/>
    <map_variables variable_2="n" variable_1="n"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
    <map_variables variable_2="lamda" variable_1="lamda"/>
    <map_variables variable_2="time" variable_1="time"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>


<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#0153fe15-e7aa-4e5d-95ea-4cfe44f6ac16">
    <rdf:li>gastrointestinal</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>biomechanics</rdf:li>
    <rdf:li>smooth muscle</rdf:li>
    <rdf:li>Smooth Muscle Cell</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#fd70f18d-d0f2-491a-a245-b94378c5bd6c">
    <rdf:li rdf:resource="rdf:#76b00b07-9afd-4812-83be-bf4e5990f76f"/>
    <rdf:li rdf:resource="rdf:#fcb5fe8c-662e-4a75-8485-2602f53ce02d"/>
    <rdf:li rdf:resource="rdf:#938a2fbe-f316-4029-8f79-a3f1dea9599c"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#fcb5fe8c-662e-4a75-8485-2602f53ce02d">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#3a2bc57f-5dd5-445b-a5d4-92943ee918c6"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#938a2fbe-f316-4029-8f79-a3f1dea9599c">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#a552ac2e-35f0-4eb6-88eb-b83769292350"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5628e094-3c3c-4c8f-98ca-83475e5ccd1f">
    <dcterms:W3CDTF>2002-07-11</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fa6b4d7a-3a9a-4996-a394-53894511d411">
    <vCard:N rdf:resource="rdf:#3e01ff8e-a8dd-4839-81c7-25c0ee03e1e8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d793a3d9-4fb5-4988-ae8c-b6313271e639">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="#miftakhov_abdusheva_christensen_1999_version01">
    <dc:title>
        The 1999 Miftakhov et al. Model of Smooth Muscle Electrical Activity.
      </dc:title>
    <cmeta:bio_entity>Smooth Muscle Cell</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#fdcd93ef-3a65-404a-9a79-2fab125f8bac"/>
    <bqs:reference rdf:resource="rdf:#616dab21-5f4c-4811-b347-04e9f3c8f8da"/>
    <bqs:reference rdf:resource="rdf:#a23cf77d-9692-425c-bf7f-2907218fd7d4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#64ddb1a0-1fa7-4690-8528-60c49c9453a3">
    <dc:creator rdf:resource="rdf:#fd70f18d-d0f2-491a-a245-b94378c5bd6c"/>
    <dc:title>
            Numerical Simulation of Motility Patterns of the Small Bowel. I.                 Formulation of a Mathematical Model
          </dc:title>
    <bqs:volume>197</bqs:volume>
    <bqs:first_page>89</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#390d3e54-bc68-4a7e-b3c1-15448f4047f9"/>
    <dcterms:issued rdf:resource="rdf:#df11f0bf-420d-43bd-a180-d771b658d81e"/>
    <bqs:last_page>112</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a23cf77d-9692-425c-bf7f-2907218fd7d4">
    <dc:subject rdf:resource="rdf:#6ef39430-4a39-4335-a4f2-71a1899cc1d7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f8e418aa-9efa-4de5-b688-58bc0b2bf61e">
    <vCard:N rdf:resource="rdf:#634daa5b-9cc0-482e-a4a3-b26e6034db54"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#052cca2f-1746-4d54-bb14-88b24f55db3a">
    <vCard:Given>R</vCard:Given>
    <vCard:Family>Miftakhov</vCard:Family>
    <vCard:Other>N</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f5f58381-0d91-441e-94da-83cd59e9ea89">
    <dcterms:modified rdf:resource="rdf:#ec76836b-391c-4744-832a-9952b73a93ab"/>
    <rdf:value>
          Changed initial values to make them valid cellml
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#fa6b4d7a-3a9a-4996-a394-53894511d411"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d1cc68d0-b0d3-4bf2-99e6-8305cf68efdd">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#616dab21-5f4c-4811-b347-04e9f3c8f8da">
    <bqs:Pubmed_id>10036210</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#64ddb1a0-1fa7-4690-8528-60c49c9453a3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#762f1974-8fb7-414b-8147-77bded8cde9d">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b7dd572d-f61b-4372-9949-277f85bbaf5b">
    <dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ec76836b-391c-4744-832a-9952b73a93ab">
    <dcterms:W3CDTF>2005-04-27</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#df11f0bf-420d-43bd-a180-d771b658d81e">
    <dcterms:W3CDTF>1999-03-07</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6ef39430-4a39-4335-a4f2-71a1899cc1d7">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#0153fe15-e7aa-4e5d-95ea-4cfe44f6ac16"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4481c5a5-b827-415e-ad6f-65c85a9c7796">
    <dcterms:modified rdf:resource="rdf:#4d064cae-e8b8-4300-a3ae-f56860c3645f"/>
    <rdf:value>
          Added more metadata.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#f8e418aa-9efa-4de5-b688-58bc0b2bf61e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#634daa5b-9cc0-482e-a4a3-b26e6034db54">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, The Bioengineering Institute
      </dc:publisher>
    <cmeta:modification rdf:resource="rdf:#4481c5a5-b827-415e-ad6f-65c85a9c7796"/>
    <cmeta:modification rdf:resource="rdf:#47dad6cc-9123-4644-a27e-77321b9deb35"/>
    <cmeta:modification rdf:resource="rdf:#f5f58381-0d91-441e-94da-83cd59e9ea89"/>
    <dcterms:created rdf:resource="rdf:#5628e094-3c3c-4c8f-98ca-83475e5ccd1f"/>
    <dc:creator rdf:resource="rdf:#16b03f9b-85ec-4b8f-87cb-b7969c4e3d1f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3e01ff8e-a8dd-4839-81c7-25c0ee03e1e8">
    <vCard:Given>Peter</vCard:Given>
    <vCard:Family>Villiger</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#47dad6cc-9123-4644-a27e-77321b9deb35">
    <dcterms:modified rdf:resource="rdf:#b7dd572d-f61b-4372-9949-277f85bbaf5b"/>
    <rdf:value>
          Added publication date information.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#fe4b8aca-a6e5-4245-ab8c-a9acd2e2e15b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a552ac2e-35f0-4eb6-88eb-b83769292350">
    <vCard:Given>J</vCard:Given>
    <vCard:Family>Christensen</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4d064cae-e8b8-4300-a3ae-f56860c3645f">
    <dcterms:W3CDTF>2002-07-18</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#390d3e54-bc68-4a7e-b3c1-15448f4047f9">
    <dc:title>Journal of theoretical Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3a2bc57f-5dd5-445b-a5d4-92943ee918c6">
    <vCard:Given>G</vCard:Given>
    <vCard:Family>Abdusheva</vCard:Family>
    <vCard:Other>R</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#90c92e2f-5e3a-4015-ba27-f95e5b33d833">
    <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:#fe4b8aca-a6e5-4245-ab8c-a9acd2e2e15b">
    <vCard:N rdf:resource="rdf:#d793a3d9-4fb5-4988-ae8c-b6313271e639"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#76b00b07-9afd-4812-83be-bf4e5990f76f">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#052cca2f-1746-4d54-bb14-88b24f55db3a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3fca3a36-b0fe-4fd5-8c9c-bcff9e204a44">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fdcd93ef-3a65-404a-9a79-2fab125f8bac">
    <dc:creator rdf:resource="rdf:#d1cc68d0-b0d3-4bf2-99e6-8305cf68efdd"/>
    <rdf:value>
          This is the CellML description of Miftakhov et al's 1999 mathematical 
          model of the electrical activity of a smooth muscle cell from the 
          small intestine. 
        </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#16b03f9b-85ec-4b8f-87cb-b7969c4e3d1f">
    <vCard:ORG rdf:resource="rdf:#762f1974-8fb7-414b-8147-77bded8cde9d"/>
    <vCard:EMAIL rdf:resource="rdf:#90c92e2f-5e3a-4015-ba27-f95e5b33d833"/>
    <vCard:N rdf:resource="rdf:#3fca3a36-b0fe-4fd5-8c9c-bcff9e204a44"/>
  </rdf:Description>
</rdf:RDF>
</model>