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_1999" name="miftakhov_1999">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Numerical simulation of motility patterns of the small bowel. 1. formulation of a mathematical model</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
           This CellML model represents the smooth muscle cell described in the published paper (equations 20-27). The model runs in both OpenCell and COR and the units are consistent. However the model does not replicate the results in figure 2 of the paper - this is most likely to be due to a lack of a stimulus input. The original paper describes several cell types which are linked together, and a mechanical stimulus which is detected by a neuron (eq. 31). Therefore to get this model to run correctly the other cell models will have to be coded up and linked together. CellML 1.1 is suited for this purpose.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: A complete mathematical model of the periodic myoelectrical activity of a functional unit of the small intestine is presented. Based on real morphological and electrophysiological data, the model assumes that: the functional unit is an electromyogenic syncytium; the kinetics of L-type Ca2+, T-type Ca2+, Ca2+-activated K+, voltage dependent K+and Cl-channels determine the electrical activity of the functional unit; the enteric nervous system is satisfactorily represented by an efferent cholinergic neuron that provides an excitatory input to the functional unit through receptor-linked L-type Ca2+channels and by an afferent pathway composed of the primary and secondary sensory neurons; the dynamics of propagation of the wave of depolarization along the unmyelinated nerve axons satisfy the Hodgkin-Huxley model; the electrical activity of the neural soma reflects the interaction of N-type Ca2+channels, Ca2+-activated K+and voltage dependent Na+, K+and Cl-channels; the smooth muscle syncytium of the locus is a null-dimensional contractile system. With the proposed model the dynamics of active force generation are determined entirely by the concentration of cytosolic calcium. The model describes: the mechanical excitation of the free nerve endings of the mechanoreceptor of the receptive field of the pathway; the electrical processes of the propagation of excitation along the afferent and efferent neural circuits; the chemical mechanisms of nerve-pulse transmission at the synaptic zones; the slow wave and bursting type electrical activity; cytosolic calcium concentration; the dynamics of active force generation. Numerical simulations have shown that the model can display different electrical patterns and mechanical responses of the locus. The results show good qualitative and quantitative agreement with the results of experiments conducted on the small intestine.
</para>

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

<para>
Numerical Simulation of Motility Patterns of the Small Bowel. 1. Formulation of a Mathematical Model, R. N. Miftakhov, G. R. Abdusheva and J. Christensen, 1999, <emphasis>Journal of Theoretical Biology</emphasis>, 197, 89-112. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/10036210">PubMed ID: 10036210</ulink> 
</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 Ca2+ L and T-type channels (ICa,L and ICa,T), a voltage-gated K+ channel (IK), a calcium-activated K+ channel (ICa-K), and a leak chloride current (ICl).</caption>
</informalfigure>

</sect1>
</article>
</documentation>
  
  
  <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="millimolar_per_millivolt">
    <unit units="millimolar"/>
    <unit units="millivolt" 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="gate_voltage">
    <variable units="millivolt" public_interface="out" name="V_tilde"/>
	
	<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
    
	<math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="V_tilde_calculation">
        <eq/>
        <ci> V_tilde </ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
			<apply>
              <times/>
			  <cn cellml:units="dimensionless"> 127.0 </cn>
              <ci> V </ci>
            </apply>
			<cn cellml:units="millivolt"> 8265.0 </cn>
		  </apply>
          <cn cellml:units="dimensionless"> 105.0 </cn>
        </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" name="V"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V_tilde"/>
    <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_tilde"/>
    
    <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_tilde </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 5.0 </cn>
                <apply>
                  <times/>
                  <ci> V_tilde </ci>
                  <cn cellml:units="per_millivolt"> 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_tilde </ci>
              </apply>
              <cn cellml:units="millivolt"> 18.0 </cn>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="m_equation">
        <eq/>
        <ci> m </ci>
        <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" initial_value="0.01" 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_tilde"/>
    <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_tilde </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="per_millisecond"> 1.0 </cn>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 5.5 </cn>
                <apply>
                  <times/>
                  <ci> V_tilde </ci>
                  <cn cellml:units="per_millivolt"> 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="dimensionless"> 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" name="V"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V_tilde"/>
    
    <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" initial_value="0.01" 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_tilde"/>
    <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="per_millivolt"> 0.15 </cn>
                    <apply>
                      <minus/>
                      <apply>
                        <minus/>
                        <ci> V_tilde </ci>
                      </apply>
                      <cn cellml:units="millivolt"> 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" name="V"/>
    <variable units="millivolt" public_interface="in" private_interface="out" name="V_tilde"/>
    <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" initial_value="0.01" 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_tilde"/>
    <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_tilde </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <minus/>
                  <cn cellml:units="millivolt"> 55.0 </cn>
                  <ci> V_tilde </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_tilde </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="dimensionless"> 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="millimolar_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="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="gate_voltage" component_1="membrane"/>
    <map_variables variable_2="V" variable_1="V"/>
  </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="L_type_calcium_current" component_1="gate_voltage"/>
    <map_variables variable_2="V_tilde" variable_1="V_tilde"/>
  </connection>
  
  <connection>
    <map_components component_2="potassium_current" component_1="gate_voltage"/>
    <map_variables variable_2="V_tilde" variable_1="V_tilde"/>
  </connection>
  
  <connection>
    <map_components component_2="T_type_calcium_current" component_1="gate_voltage"/>
    <map_variables variable_2="V_tilde" variable_1="V_tilde"/>
  </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="V_tilde" variable_1="V_tilde"/>
  </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_tilde" variable_1="V_tilde"/>
  </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_tilde" variable_1="V_tilde"/>
  </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_tilde" variable_1="V_tilde"/>
  </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_1999">
    <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>Auckland 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>
        Auckland Bioengineering Institute, The University of Auckland
      </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>