Location: Corrias, Buist, 2007 @ 6ecf1843c9a0 / corrias_buist_2007.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2009-12-07 12:14:03+13:00
Desc:
Added images in ai and svg format, removed non pub med references
Permanent Source URI:
https://models.cellml.org/workspace/corrias_buist_2007/rawfile/6ecf1843c9a0c605e502f56c2096d481c72fea43/corrias_buist_2007.cellml

<?xml version='1.0' encoding='utf-8'?>
<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:ns8="http://www.cellml.org/metadata/graphs/1.0#" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" name="SM_model" cmeta:id="corrias_model_2007">
	
  <documentation xmlns="http://cellml.org/tmp-documentation">
    <article>
      <articleinfo>
        <title>A quantitative model of gastric smooth muscle cellular activation</title>
        <author>
          <firstname>Alberto</firstname>
          <surname>Corrias</surname>
          <affiliation>
            <shortaffil>Division of Bioengineering, National University of Singapore</shortaffil>
          </affiliation>
        </author>
      </articleinfo>
      <section id="sec_status">
        <title>Model Status</title>
        <para>
          This CellML version of the model has been encoded by the model author, Alberto Corrias, and is known to produce correct results in JSim and CellMLSimulator. It also passes the online PyCml validation tests with no errors or warnings. This model includes the specification of four pseudo ICC stimulus pulses with a 20s period. More information and simulation results can be found at the <ulink url="http://www.bioeng.nus.edu.sg/compbiolab/projects/gi-project.html#sm-cell-ephys">Computational Bioengineering Laboratory</ulink>.
        </para>
      </section>
      <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
The coordinated activity of the enteric nervous system (ENS), interstitial cells of Cajal (ICC), and smooth muscle (SM) cells is required to generate gastric motility.  While the ICC provide a constant source of electrical activity, and can be regarded as the <emphasis>pacemaker</emphasis> of the gastrointestinal (GI) system, the ENS adds an extra level of control.  In combination the ICC and ENS act as a stimulus, causing the SM cells to contract and generate slow wave activity in the GI.  These regular, autonomous depolarisations are then propagated via gap junctions to the neighbouring SM cells in the gut wall. 
</para>

<para>
SM cells themselves are unable to generate or propagate slow waves autonomously.  Instead they response to the electrical stimuli generated by neighbouring ICC as voltage-dependent ion channels embedded within their cell surface membrane become activated and other extracellular ligands and intracellular second messengers (such as calcium ions) also modify the response of the SM cells.  As such, the inherent complexity of such a system quickly becomes apparent, and the role of mathematical models in dissecting the system and characterising the different functions of each individual component becomes clear.  However, while computational models of cardiac electrophysiology have advanced considerably over the past few years, to date, models of gastric electrophysiology have been limited by the absence of single cell models which form the foundation for the more complex, large scale, integrated models.  
</para>

<para>
In the present study Alberto Corrias and Martin Buist address this issue by presenting an electrophysiolgical model of a single gastric SM cell.  The model parameters are based on real, experimentally recorded, physiological data, and the model has also been kept relatively simple to make it computationally efficient, and therefore suitable for embedding within multicellular simulations.
</para>

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

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

<para>
A quantitative model of gastric smooth muscle cellular activation, Alberto Corrias and Martin L. Buist, 2007, <emphasis>Annals of Biomedical Engineering</emphasis>, 35, 1595-1607. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=17486452&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 17486452</ulink>
</para>

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

  
	
	<!--UNIT DECLARATION-->
	<units name="conductance_units">
		<unit units="siemens" prefix="nano"/>
	</units>
	<units name="millifarads">
		<unit units="farad" prefix="milli"/>
	</units>
	<units name="voltage_units">
		<unit units="volt" prefix="milli"/>
	</units>
	<units name="volume_units">
		<unit units="meter" prefix="micro" exponent="3"/>
	</units>
	<units name="Inverse_Voltage_units">
		<unit units="voltage_units" exponent="-1"/>
	</units>
	<units name="Inverse_Voltage_Units_per_time_units">
		<unit units="Inverse_Voltage_units"/> 
		<unit units="time_units" exponent="-1"/>
	</units>    
	<units name="time_units">
		<unit units="second" prefix="milli"/>
	</units>
	<units name="current_units">
		<unit units="voltage_units"/>
		<unit units="conductance_units"/>
	</units>
	<units name="rate_constants_units">
		<unit units="time_units" exponent="-1"/>
	</units>
	<units name="capacitance_units">
		<unit units="farad" prefix="pico"/>
	</units>
	<units name="Temperature_units">
		<unit units="kelvin"/>
	</units>
	<units name="Inverse_Temperature_units_times_conductance">
		<unit units="Temperature_units" exponent="-1"/>
		<unit units="conductance_units"/>
	</units>
	<units name="F_units">
		<unit units="mole" exponent="-1"/>
		<unit units="coulomb"/>
	</units>
	<units name="R_units">
		<unit units="joule" prefix="milli"/>
		<unit units="mole" exponent="-1"/>
		<unit units="kelvin" exponent="-1"/>
	</units>
	<units name="millimolar">
		<unit units="mole" prefix="milli"/>
		<unit units="litre" exponent="-1"/>
	</units>
	<units name="per_millimolar">
		<unit units="litre"/>
		<unit units="mole" prefix="milli" exponent="-1"/>
	</units>
	<units name="millimolar_per_millisecond">
		<unit units="mole" prefix="milli"/>
		<unit units="litre" exponent="-1"/>
		<unit units="second" prefix="milli" exponent="-1"/>
	</units>
	<units name="millimolar_per_millisecond_per_millivolt">
		<unit units="mole" prefix="milli"/>
		<unit units="litre" exponent="-1"/>
		<unit units="second" exponent="-1"/>
		<unit units="volt" prefix="milli" exponent="-1"/>
	</units>
	<units name="mole_to_millimole">
		<unit units="mole"/>
		<unit units="mole" prefix="milli" exponent="-1"/>
	</units>
	
	<!--END OF UNIT DECLARATION-->
	
	<component name="Time">
		<variable units="time_units" public_interface="out" private_interface="out" name="time" cmeta:id="time"/>
	</component> <!-- closure of time component-->	
	
	<component name="Neural_input">
		<variable units="millimolar" public_interface="out" private_interface="out" name="Ach" initial_value="0.00001"/>
	</component> <!-- closure of neural input component-->	
	
	<component name="Gap_junction">
		<variable units="conductance_units" public_interface="out" private_interface="out" name="Gcouple" initial_value="1.3"/>
	</component>
	
	<component name="Environment">
		<variable units="Temperature_units" private_interface="out" name="T" initial_value="310"/>
		<variable units="Temperature_units" private_interface="out" name="T_exp" initial_value="297"/>
		<variable units="F_units" public_interface="out" private_interface="out" name="F" initial_value="96486"/>
		<variable units="R_units" private_interface="out" name="R" initial_value="8314.4"/>
		<variable units="dimensionless" private_interface="out" name="Q10Ca" initial_value="2.1"/>
		<variable units="dimensionless" private_interface="out" name="Q10K" initial_value="1.365"/>
		<variable units="dimensionless" private_interface="out" name="Q10Na" initial_value="2.45"/>
		<variable units="millimolar" public_interface="out" private_interface="out" name="Ca_o" initial_value="2.5"/>
		<variable units="millimolar" public_interface="out" private_interface="out" name="Na_o" initial_value="137"/>
		<variable units="millimolar" public_interface="out" private_interface="out" name="K_o" initial_value="5.9"/>
		<variable units="millimolar" public_interface="out" private_interface="out" name="Cl_o" initial_value="134"/>
		<variable units="dimensionless" public_interface="out" private_interface="out" name="T_correction_Na"/>	
		<variable units="dimensionless" public_interface="out" private_interface="out" name="T_correction_K"/>
		<variable units="dimensionless" public_interface="out" private_interface="out" name="T_correction_Ca"/>
		<variable units="conductance_units" public_interface="out" private_interface="out" name="T_correction_BK"/>
		
		<variable units="Inverse_Voltage_units" public_interface="out" private_interface="out" name="FoRT"/>
		<variable units="voltage_units" public_interface="out" private_interface="out" name="RToF"/>
		
		<!-- equation that calculates FoRT-->	
		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply>
        <eq/>		
				<ci> FoRT </ci>
				<apply>
          <divide/>
					<ci> F </ci>
					<apply>
            <times/>
						<ci> R </ci>
						<ci> T </ci>
					</apply>
				</apply>		
			</apply>
      <!-- closure of equation that calculates FoRT-->	
		</math>
		
		<!-- equation that calculates RToF-->	
		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply>
        <eq/>		
				<ci> RToF </ci>
				<apply>
          <divide/>
					<apply>
            <times/>
						<ci> R </ci>
						<ci> T </ci>
					</apply>
					<ci> F </ci>
				</apply>		
			</apply>
      <!-- closure of equation that calculates FoRT-->	
		</math>
		
		<!-- equation that calculates Tcorrection_Ca-->
		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply>
        <eq/>		
				<ci>T_correction_Ca</ci>
				<apply>
          <power/>
					<ci>Q10Ca </ci>
					<apply>
            <divide/>
						<apply>
              <minus/>
							<ci>T </ci>
							<ci>T_exp </ci>
						</apply>
						<cn cellml:units="Temperature_units">10</cn>
					</apply>
				</apply>	
			</apply>
      <!-- closure of equation that calculates TcorrectionCa-->	
		</math>
		
		<!-- equation that calculates Tcorrection_K-->
		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply>
        <eq/>		
				<ci>T_correction_K</ci>
				<apply>
          <power/>
					<ci>Q10K </ci>
					<apply>
            <divide/>
						<apply>
              <minus/>
							<ci>T </ci>
							<ci>T_exp </ci>
						</apply>
						<cn cellml:units="Temperature_units">10</cn>
					</apply>
				</apply>	
			</apply>
      <!-- closure of equation that calculates TcorrectionK-->	
		</math>
		
		<!-- equation that calculates Tcorrection_Na-->
		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply>
        <eq/>		
				<ci>T_correction_Na</ci>
				<apply>
          <power/>
					<ci>Q10Na </ci>
					<apply>
            <divide/>
						<apply>
              <minus/>
							<ci>T </ci>
							<ci>T_exp </ci>
						</apply>
						<cn cellml:units="Temperature_units">10</cn>
					</apply>
				</apply>	
			</apply>
      <!-- closure of equation that calculates TcorrectionNa-->	
		</math>
		<!-- equation that calculates Tcorrection_Na-->
		<math xmlns="http://www.w3.org/1998/Math/MathML">
			<apply>
        <eq/>		
				<ci>T_correction_BK</ci>
				<apply>
          <times/>
					<cn cellml:units="Inverse_Temperature_units_times_conductance">1.1</cn>
					<apply>
            <minus/>
						<ci>T </ci>
						<ci>T_exp </ci>
					</apply>
				</apply>	
			</apply>
      <!-- closure of equation that calculates TcorrectionNa-->	
		</math>
	</component> <!-- closure of environment component-->
	
	<component name="SM_Membrane">
		
		<variable units="capacitance_units" name="Cm_SM" initial_value="77"/>
		<variable units="volume_units" public_interface="out" private_interface="out" name="Vol_SM" initial_value="3500"/>
		<variable units="voltage_units" private_interface="out" name="Vm_SM" cmeta:id="Vm" initial_value="-69.75"/>
		<variable units="time_units" public_interface="in" private_interface="out" name="time"/>
		<variable units="millimolar" private_interface="out" name="Ca_i" cmeta:id="Cai" initial_value="0.00008"/>
		<variable units="millimolar" public_interface="out" private_interface="out" name="Na_i" initial_value="10"/>
		<variable units="millimolar" private_interface="out" name="K_i" initial_value="164"/>
		<variable units="millimolar" public_interface="in" private_interface="out" name="Ca_o"/>
		<variable units="millimolar" public_interface="in" private_interface="out" name="Ach"/>
		<variable units="millimolar" public_interface="in" private_interface="out" name="Na_o"/>
		<variable units="millimolar" public_interface="in" private_interface="out" name="K_o"/>
		<variable units="millimolar" public_interface="in" private_interface="out" name="Cl_o"/>
		<variable units="dimensionless" public_interface="in" private_interface="out" name="T_correction_Ca"/>
		<variable units="dimensionless" public_interface="in" private_interface="out" name="T_correction_Na"/>
		<variable units="dimensionless" public_interface="in" private_interface="out" name="T_correction_K"/>
		<variable units="conductance_units" public_interface="in" private_interface="out" name="T_correction_BK"/>
		<variable units="voltage_units" public_interface="in" private_interface="out" name="RToF"/>
		<variable units="Inverse_Voltage_units" public_interface="in" private_interface="out" name="FoRT"/>
		<variable units="F_units" public_interface="in" private_interface="out" name="F"/>
		
		
		<variable units="current_units" private_interface="in" name="I_Na_SM" cmeta:id="INa"/>
		<variable units="current_units" private_interface="in" name="I_Ltype_SM" cmeta:id="ICaL"/>
		<variable units="current_units" private_interface="in" name="I_LVA_SM" cmeta:id="ILVA"/>
		<variable units="current_units" private_interface="in" name="I_kr_SM" cmeta:id="IKr"/>
		<variable units="current_units" private_interface="in" name="I_ka_SM" cmeta:id="IKA"/>
		<variable units="current_units" private_interface="in" name="I_BK_SM" cmeta:id="IBK"/>
		<variable units="current_units" private_interface="in" name="I_NSCC_SM"/>
		<variable units="current_units" private_interface="in" name="I_bk_SM"/>
		<variable units="millimolar_per_millisecond" private_interface="in" name="J_CaSR_SM"/>
		<variable units="current_units" public_interface="in" name="I_stim"/>
		

		<!-- equation that calculates dVm_SM/dt-->	     	
		<math xmlns="http://www.w3.org/1998/Math/MathML">   	
			<apply>
        <eq/>		
				<apply>
          <diff/>
					<bvar>
            <ci>time </ci>
          </bvar>
					<ci> Vm_SM </ci>
				</apply>
				<apply>
          <times/>
					<cn cellml:units="dimensionless">-1.0</cn>
					<apply>
            <divide/>
						<cn cellml:units="dimensionless">1.0</cn>
						<ci>Cm_SM</ci>
					</apply>
					<apply>
            <plus/>   			
						<ci> I_Na_SM </ci>
						<ci> I_Ltype_SM </ci>
						<ci> I_LVA_SM </ci>
						<ci> I_kr_SM </ci>
						<ci> I_ka_SM </ci>
						<ci> I_BK_SM </ci>
						<ci> I_NSCC_SM </ci>
						<ci> I_bk_SM </ci>
						<apply>
              <times/>
							<cn cellml:units="dimensionless">-1.0</cn>
							<ci> I_stim </ci>
						</apply>
					</apply>
				</apply>
        <!-- closure of equation dVm_SM/dt= (-1/Cm)*(sum of ionic currents)-->	
			</apply>	
		</math> 
		
		<!-- equation that calculates dCai/dt-->
		<math xmlns="http://www.w3.org/1998/Math/MathML">   	
			<apply>
        <eq/>		
				<apply>
          <diff/>
					<bvar>
            <ci>time </ci>
          </bvar>
					<ci> Ca_i </ci>
				</apply>
				<apply>
          <plus/>
					<apply>
            <divide/>
						<apply>
              <plus/>
							<apply>
                <times/>
								<cn cellml:units="dimensionless">-1</cn>
								<ci> I_Ltype_SM </ci>
							</apply>
							<apply>
                <times/>
								<cn cellml:units="dimensionless">-1</cn>
								<ci> I_LVA_SM </ci>
							</apply>
						</apply>
						<apply>
              <times/>
							<cn cellml:units="dimensionless">2.0</cn>
							<cn cellml:units="mole_to_millimole">0.001</cn>
							<ci> F </ci>
							<ci> Vol_SM </ci>
						</apply>
					</apply>
					<apply>
            <times/>
						<cn cellml:units="dimensionless">-1</cn>
						<ci>J_CaSR_SM</ci>
					</apply>
				</apply>	
			</apply>
      <!--closure of equation Ca_i-->	
		</math>
		

</component> <!-- closure of membrane component-->

<component cmeta:id="ICC_stimulus" name="I_stim">
	
        <variable units="time_units" public_interface="in" name="time"/>
	<variable units="time_units" name="local_time" cmeta:id="local_time"/>
        <variable units="time_units" name="period" initial_value="20000"/>
        <variable units="time_units" name="stim_start" cmeta:id="stim_start"/>
	<variable units="voltage_units" name="delta_VICC" initial_value="59"/>
	<variable units="conductance_units" public_interface="in" name="Gcouple"/>
	<variable units="time_units" name="t_ICCpeak" initial_value="98"/>
	<variable units="time_units" name="t_ICCplateau" initial_value="7582"/>
	<variable units="time_units" name="t_ICC_stimulus" initial_value="10000"/>
	<variable units="voltage_units" name="V_decay" initial_value="37.25"/>
        <variable units="current_units" public_interface="out" name="I_stim" cmeta:id="I_stim"/>
        
	
        <!-- set the current stimulus start time -->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
          <apply>
        <eq/>
            <ci>stim_start</ci>
            <piecewise>
              <piece>
                <apply>
              <times/>
                  <ci>period</ci>
                  <cn cellml:units="dimensionless">1.0</cn>
                </apply>
                <apply>
              <and/>
                  <apply>
                <gt/>
                    <ci>time</ci>
                    <apply>
                  <times/>
                      <ci>period</ci>
                      <cn cellml:units="dimensionless">1.0</cn>
                    </apply>
                  </apply>
                  <apply>
                <leq/>
                    <ci>time</ci>
                    <apply>
                  <times/>
                      <ci>period</ci>
                      <cn cellml:units="dimensionless">2.0</cn>
                    </apply>
                  </apply>
                </apply>
              </piece>
              <piece>
                <apply>
              <times/>
                  <ci>period</ci>
                  <cn cellml:units="dimensionless">2.0</cn>
                </apply>
                <apply>
              <and/>
                  <apply>
                <gt/>
                    <ci>time</ci>
                    <apply>
                  <times/>
                      <ci>period</ci>
                      <cn cellml:units="dimensionless">2.0</cn>
                    </apply>
                  </apply>
                  <apply>
                <leq/>
                    <ci>time</ci>
                    <apply>
                  <times/>
                      <ci>period</ci>
                      <cn cellml:units="dimensionless">3.0</cn>
                    </apply>
                  </apply>
                </apply>
              </piece>
              <piece>
                <apply>
              <times/>
                  <ci>period</ci>
                  <cn cellml:units="dimensionless">3.0</cn>
                </apply>
                <apply>
              <and/>
                  <apply>
                <gt/>
                    <ci>time</ci>
                    <apply>
                  <times/>
                      <ci>period</ci>
                      <cn cellml:units="dimensionless">3.0</cn>
                    </apply>
                  </apply>
                  <apply>
                <leq/>
                    <ci>time</ci>
                    <apply>
                  <times/>
                      <ci>period</ci>
                      <cn cellml:units="dimensionless">4.0</cn>
                    </apply>
                  </apply>
                </apply>
              </piece>
              <piece>
                <apply>
              <times/>
                  <ci>period</ci>
                  <cn cellml:units="dimensionless">4.0</cn>
                </apply>
                <apply>
              <and/>
                  <apply>
                <gt/>
                    <ci>time</ci>
                    <apply>
                  <times/>
                      <ci>period</ci>
                      <cn cellml:units="dimensionless">4.0</cn>
                    </apply>
                  </apply>
                  <apply>
                <leq/>
                    <ci>time</ci>
                    <apply>
                  <times/>
                      <ci>period</ci>
                      <cn cellml:units="dimensionless">5.0</cn>
                    </apply>
                  </apply>
                </apply>
              </piece>
              <otherwise>
                <cn cellml:units="time_units">0.0</cn>
              </otherwise>
            </piecewise>
          </apply>
        </math>
        <!-- calculate the local time to use for the stimulus protocol calculation -->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
          <apply>
        <eq/>
            <ci>local_time</ci>
            <apply>
          <minus/>
              <ci>time</ci>
              <apply>
            <plus/>
                <ci>stim_start</ci>
                <ci>t_ICCpeak</ci>
              </apply>
            </apply>
          </apply>
        </math>
<!-- equation that calculates I_stim-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> I_stim </ci>
		<piecewise>
			<piece>
				<apply>
              <times/>
					<ci> Gcouple </ci>
					<ci> delta_VICC </ci>
            			</apply>
				<apply>
              <lt/>
					<ci>local_time</ci>
					<ci>t_ICCpeak</ci>
				</apply>
			</piece>
			<piece>
				<apply>
              <times/>
					<ci> Gcouple </ci>
					<ci> delta_VICC </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> local_time </ci>
										<cn cellml:units="time_units"> 8000 </cn>
									</apply>
									<cn cellml:units="time_units"> 1000 </cn>
								</apply>
							</apply>
						</apply>
					</apply>
				</apply>
				<apply>
              <and/>
					<apply>
                <geq/>
						<ci> local_time </ci>
						<ci> t_ICCpeak </ci>
					</apply>
					<apply>
                <leq/>
						<ci> local_time </ci>
						<ci> t_ICCplateau </ci>
					</apply>
            			</apply>
			</piece>	
			<piece>
				<apply>
              <times/>
					<ci> Gcouple </ci>
					<ci> V_decay </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> local_time </ci>
										<cn cellml:units="time_units"> 8000 </cn>
									</apply>
									<cn cellml:units="time_units"> 150 </cn>
								</apply>
							</apply>
						</apply>
					</apply>
				</apply>
				<apply>
              <and/>
					<apply>
                <gt/>
						<ci> local_time </ci>
						<ci> t_ICCplateau  </ci>
					</apply>
					<apply>
                <lt/>
						<ci> local_time </ci>
						<ci> t_ICC_stimulus</ci>
					</apply>
				</apply>
			</piece>
			<otherwise>
				<cn cellml:units="current_units"> 0.0 </cn>
			</otherwise>
		</piecewise>
              </apply>
	</math>
</component>
<component name="d_Ltype_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="d_inf_Ltype_SM"/>
	<variable units="time_units" name="tau_d_Ltype_SM"/>
	<variable units="dimensionless" public_interface="out" name="d_Ltype_SM" initial_value="0.0"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_Ca"/>
	
	<!-- equation that calculates d_inf_Ltype-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> d_inf_Ltype_SM </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> Vm_SM </ci>
								<cn cellml:units="voltage_units">17</cn>
							</apply>
							<cn cellml:units="voltage_units">-4.3</cn>
						</apply>
					</apply>
				</apply>
			</apply>		
		</apply>
      <!-- closure of equation that calculates d_inf_Ltype-->	
	</math>
	
	<!-- equation that calculates tau_d_Ltype-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_d_Ltype_SM </ci>
			<apply>
          <times/>
				<ci>T_correction_Ca</ci>
				<cn cellml:units="time_units">0.47</cn>
			</apply>	
		</apply>
      <!-- closure of equation that calculates tau_d_Ltype-->	
	</math>
	<!-- equation that calculates d_Ltype-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> d_Ltype_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> d_inf_Ltype_SM </ci>
					<ci>d_Ltype_SM</ci>
				</apply>
				<ci> tau_d_Ltype_SM </ci>
			</apply>
		</apply>
      <!-- closure of equation d(d_Ltype)/dt=(d_inf_Ltype-d_Ltype)/tau_d_Ltype-->	
	</math> 
</component> <!-- closure of d_Ltype-->

<component name="f_Ltype_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="f_inf_Ltype_SM"/>
	<variable units="time_units" name="tau_f_Ltype_SM"/>
	<variable units="dimensionless" public_interface="out" name="f_Ltype_SM" initial_value="0.95"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_Ca"/>
	
	<!-- equation that calculates f_inf_Ltype-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> f_inf_Ltype_SM </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> Vm_SM </ci>
								<cn cellml:units="voltage_units">43</cn>
							</apply>
							<cn cellml:units="voltage_units">8.9</cn>
						</apply>
					</apply>
				</apply>
			</apply>		
		</apply>
      <!-- closure of equation that calculates f_inf_Ltype-->	
	</math>
	
	<!-- equation that calculates tau_f_Ltype-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_f_Ltype_SM </ci>
			<apply>
          <times/>
				<ci>T_correction_Ca</ci>
				<cn cellml:units="time_units">86</cn>	
			</apply>
        <!-- closure of equation that calculates tau_f_Ltype-->	
		</apply>
	</math>
	<!-- equation that calculates f_Ltype-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> f_Ltype_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> f_inf_Ltype_SM </ci>
					<ci>f_Ltype_SM</ci>
				</apply>
				<ci> tau_f_Ltype_SM </ci>
			</apply>
		</apply>
      <!-- closure of equation d(f_Ltype)/dt=(f_inf_Ltype-f_Ltype)/tau_f_Ltype-->	
	</math> 
</component> <!-- closure of f_Ltype-->

<component name="f_ca_Ltype_SM">
	<variable units="millimolar" public_interface="in" name="Ca_i"/>
	<variable units="dimensionless" name="f_ca_inf_Ltype_SM"/>
	<variable units="time_units" name="tau_f_ca_Ltype_SM"/>
	<variable units="dimensionless" public_interface="out" name="f_ca_Ltype_SM" initial_value="1.0"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_Ca"/>
	
	<!-- equation that calculates f_inf_ca Ltype-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> f_ca_inf_Ltype_SM </ci>
			<apply>
          <minus/>
				<cn cellml:units="dimensionless">1.0</cn>
				<apply>
            <divide/>
					<cn cellml:units="dimensionless">1.0</cn>
					<apply>
              <plus/>
						<cn cellml:units="dimensionless">1.0</cn>
						<apply>
                <exp/>
							<apply>
                  <divide/>
								<apply>
                    <minus/>
									<apply>
                      <minus/>
										<ci> Ca_i </ci>
										<cn cellml:units="millimolar">0.00008999</cn>
									</apply>
									<cn cellml:units="millimolar">0.000214</cn>
								</apply>
								<cn cellml:units="millimolar">-0.0000131</cn>
							</apply>
						</apply>
					</apply>
				</apply>
			</apply>		
		</apply>
      <!-- closure of equation that calculates f_ca_inf_Ltype-->	
	</math>
	
	<!-- equation that calculates tau_f_ca_Ltype-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_f_ca_Ltype_SM </ci>
			<apply>
          <times/>
				<ci>T_correction_Ca</ci>
				<cn cellml:units="time_units">2</cn>	
			</apply>
        <!-- closure of equation that calculates tau_f_ca_Ltype-->	
		</apply>
	</math>
	<!-- equation that calculates f_ca_Ltype-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> f_ca_Ltype_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> f_ca_inf_Ltype_SM </ci>
					<ci>f_ca_Ltype_SM</ci>
				</apply>
				<ci> tau_f_ca_Ltype_SM </ci>
			</apply>
		</apply>
      <!-- closure of equation d(f_ca_Ltype)/dt=(f_ca_inf_Ltype-f_ca_Ltype)/tau_f_ca_Ltype-->	
	</math> 
</component> <!-- closure of f_Ltype-->
<component name="I_Ltype_SM">
	
	<variable units="current_units" public_interface="out" name="I_Ltype_SM"/>
	<variable units="voltage_units" name="E_Ca"/>
	<variable units="conductance_units" name="G_max_Ltype" initial_value="65"/>
	<variable units="voltage_units" public_interface="in" private_interface="out" name="Vm_SM"/>
	<variable units="dimensionless" private_interface="in" name="f_Ltype_SM"/>
	<variable units="dimensionless" private_interface="in" name="d_Ltype_SM"/>
	<variable units="dimensionless" private_interface="in" name="f_ca_Ltype_SM"/>
	<variable units="voltage_units" public_interface="in" name="RToF"/>
	<variable units="millimolar" public_interface="in" name="Ca_o"/>
	<variable units="millimolar" public_interface="in" private_interface="out" name="Ca_i"/>
	<variable units="dimensionless" public_interface="in" private_interface="out" name="T_correction_Ca"/>
	<variable units="time_units" public_interface="in" private_interface="out" name="time"/>
	
	<!-- EQUATION OF E_Ca-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>		
			<ci>E_Ca </ci>
			<apply>
          <times/>
				<cn cellml:units="dimensionless">0.5</cn>
				<ci> RToF </ci>
				<apply>
            <ln/>
					<apply>
              <divide/>
						<ci> Ca_o </ci>
						<ci> Ca_i </ci>
					</apply>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation E_Ca -->
	<!-- EQUATION OF I_Ltype-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>I_Ltype_SM </ci>
			<apply>
          <times/>     			
				<ci> G_max_Ltype </ci>
				<apply>
            <times/>
					<ci> f_Ltype_SM </ci>
					<ci> d_Ltype_SM </ci>
					<ci> f_ca_Ltype_SM </ci>
				</apply>
				<apply>
            <minus/>
					<ci> Vm_SM </ci>
					<ci> E_Ca </ci>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation I_Ltype -->
</component> <!-- closure of I_Ltype-->

<component name="J_CaSR_SM">
	
	<variable units="millimolar_per_millisecond" public_interface="out" name="J_CaSR_SM"/>
	<variable units="millimolar_per_millisecond" name="J_max_CaSR" initial_value="0.31705"/>
	<variable units="millimolar" public_interface="in" name="Ca_i"/>
	
	<!-- EQUATION OF J_CaSR-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>J_CaSR_SM</ci>
			<apply>
          <times/>
				<ci> J_max_CaSR </ci>
				<apply>
            <power/>
					<apply>
              <times/>
						<ci>Ca_i</ci>
						<cn cellml:units="per_millimolar">1.0</cn>
					</apply>
					<cn cellml:units="dimensionless">1.34</cn>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation J_CaSR -->
</component> <!-- closure of J_CaSR-->

<component name="d_LVA_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="d_inf_LVA_SM"/>
	<variable units="time_units" name="tau_d_LVA_SM"/>
	<variable units="dimensionless" public_interface="out" name="d_LVA_SM" initial_value="0.02"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_Ca"/>
	
	<!-- equation that calculates d_inf_LVA-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> d_inf_LVA_SM </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> Vm_SM </ci>
								<cn cellml:units="voltage_units">27.5</cn>
							</apply>
							<cn cellml:units="voltage_units">-10.9</cn>
						</apply>
					</apply>
				</apply>
			</apply>		
		</apply>
      <!-- closure of equation that calculates d_inf_LVA-->	
	</math>
	
	<!-- equation that calculates tau_d_LVA-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_d_LVA_SM </ci>
			<apply>
          <times/>
				<ci>T_correction_Ca</ci>
				<cn cellml:units="time_units">3.0</cn>
			</apply>	
		</apply>
      <!-- closure of equation that calculates tau_d_LVA-->	
	</math>
	<!-- equation that calculates d_LVA-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> d_LVA_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> d_inf_LVA_SM </ci>
					<ci>d_LVA_SM</ci>
				</apply>
				<ci> tau_d_LVA_SM </ci>
			</apply>
		</apply>
      <!-- closure of equation d(d_LVA)/dt=(d_inf_LVA-d_LVA)/tau_d_LVA-->	
	</math> 
</component> <!-- closure of d_LVA-->

<component name="f_LVA_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="f_inf_LVA_SM"/>
	<variable units="time_units" name="tau_f_LVA_SM"/>
	<variable units="dimensionless" public_interface="out" name="f_LVA_SM" initial_value="0.99"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_Ca"/>
	
	<!-- equation that calculates f_inf_LVA-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> f_inf_LVA_SM </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> Vm_SM </ci>
								<cn cellml:units="voltage_units">15.8</cn>
							</apply>
							<cn cellml:units="voltage_units">7</cn>
						</apply>
					</apply>
				</apply>
			</apply>		
		</apply>
      <!-- closure of equation that calculates f_inf_LVA-->	
	</math>
	
	<!-- equation that calculates tau_f_LVA-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_f_LVA_SM </ci>
			<apply>
          <times/>
				<ci>T_correction_Ca</ci>
				<apply>
            <times/>
					<cn cellml:units="time_units">7.58</cn>	
					<apply>
              <exp/>
						<apply>
                <times/>
							<ci>Vm_SM</ci>
							<cn cellml:units="Inverse_Voltage_units">0.00817</cn>
						</apply>
					</apply>
				</apply>
			</apply>
        <!-- closure of equation that calculates tau_f_LVA-->	
		</apply>
	</math>
	<!-- equation that calculates f_LVA-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> f_LVA_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> f_inf_LVA_SM </ci>
					<ci>f_LVA_SM</ci>
				</apply>
				<ci> tau_f_LVA_SM </ci>
			</apply>
		</apply>
      <!-- closure of equation d(f_LVA)/dt=(f_inf_LVA-f_LVA)/tau_f_LVA-->	
	</math> 
</component> <!-- closure of f_LVA-->

<component name="I_LVA_SM">
	
	<variable units="current_units" public_interface="out" name="I_LVA_SM"/>
	<variable units="voltage_units" name="E_Ca"/>
	<variable units="conductance_units" name="G_max_LVA" initial_value="0.18"/>
	<variable units="voltage_units" public_interface="in" private_interface="out" name="Vm_SM"/>
	<variable units="dimensionless" private_interface="in" name="f_LVA_SM"/>
	<variable units="dimensionless" private_interface="in" name="d_LVA_SM"/>
	<variable units="voltage_units" public_interface="in" name="RToF"/>
	<variable units="millimolar" public_interface="in" name="Ca_o"/>
	<variable units="millimolar" public_interface="in" name="Ca_i"/>
	<variable units="dimensionless" public_interface="in" private_interface="out" name="T_correction_Ca"/>
	<variable units="time_units" public_interface="in" private_interface="out" name="time"/>
	
	<!-- EQUATION OF E_Ca-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>		
			<ci>E_Ca </ci>
			<apply>
          <times/>
				<cn cellml:units="dimensionless">0.5</cn>
				<ci> RToF </ci>
				<apply>
            <ln/>
					<apply>
              <divide/>
						<ci> Ca_o </ci>
						<ci> Ca_i </ci>
					</apply>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation E_Ca -->
	<!-- EQUATION OF I_LVA-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>I_LVA_SM </ci>
			<apply>
          <times/>     			
				<ci> G_max_LVA </ci>
				<apply>
            <times/>
					<ci> f_LVA_SM </ci>
					<ci> d_LVA_SM </ci>
				</apply>
				<apply>
            <minus/>
					<ci> Vm_SM </ci>
					<ci> E_Ca </ci>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation I_LVA -->
</component> <!-- closure of I_LVA-->



<component name="d_BK_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="millimolar" public_interface="in" name="Ca_i"/>
	<variable units="dimensionless" public_interface="out" name="d_BK_SM"/>
	<variable units="time_units" public_interface="in" name="time"/>
	
	<!-- equation that calculates d_BK-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> d_BK_SM </ci>
			<apply>
          <divide/>
				<cn cellml:units="dimensionless">1.0</cn>
				<apply>
            <plus/>
					<cn cellml:units="dimensionless">1.0</cn>
					<apply>
              <exp/>
						<apply>
                <minus/>
							<apply>
                  <divide/>
								<ci> Vm_SM </ci>
								<cn cellml:units="voltage_units">-17</cn>
							</apply>
							<apply>
                  <times/>
								<cn cellml:units="dimensionless">2.0</cn>
								<apply>
                    <ln/>
									<apply>
                      <divide/>
										<ci> Ca_i </ci>
										<cn cellml:units="millimolar">0.001</cn>
									</apply>
								</apply>
							</apply>
						</apply>
					</apply>
				</apply>
			</apply>
		</apply>
      <!-- closure of equation that calculates d_inf_BK-->	
	</math>
</component> <!-- closure of d_BK-->

<component name="I_BK_SM">
	
	<variable units="current_units" public_interface="out" name="I_BK_SM"/>
	<variable units="millimolar" public_interface="in" private_interface="out" name="Ca_i"/>
	<variable units="voltage_units" name="E_K"/>
	<variable units="conductance_units" name="G_max_BK" initial_value="45.7"/>
	<variable units="voltage_units" public_interface="in" private_interface="out" name="Vm_SM"/>
	<variable units="dimensionless" private_interface="in" name="d_BK_SM"/>
	<variable units="voltage_units" public_interface="in" name="RToF"/>
	<variable units="millimolar" public_interface="in" name="K_o"/>
	<variable units="millimolar" public_interface="in" name="K_i"/>
	<variable units="conductance_units" public_interface="in" name="T_correction_BK"/>
	<variable units="time_units" public_interface="in" private_interface="out" name="time"/>
	<!-- EQUATION OF E_K-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>		
			<ci>E_K </ci>
			<apply>
          <times/>     			
				<ci> RToF </ci>
				<apply>
            <ln/>
					<apply>
              <divide/>
						<ci> K_o </ci>
						<ci> K_i </ci>
					</apply>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation E_K -->
	<!-- EQUATION OF I_BK-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>I_BK_SM </ci>
			<apply>
          <times/>
				<apply>
            <plus/>
					<ci> G_max_BK </ci>
					<ci> T_correction_BK </ci>
				</apply>
				<ci> d_BK_SM </ci>
				<apply>
            <minus/>
					<ci> Vm_SM </ci>
					<ci> E_K </ci>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation I_BK -->
</component> <!-- closure of I_BK-->

<component name="I_bk_SM">
	
	<variable units="current_units" public_interface="out" name="I_bk_SM"/>
	<variable units="voltage_units" name="E_K"/>
	<variable units="conductance_units" name="G_max_bk" initial_value="0.0144"/>
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="voltage_units" public_interface="in" name="RToF"/>
	<variable units="millimolar" public_interface="in" name="K_o"/>
	<variable units="millimolar" public_interface="in" name="K_i"/>
	<!-- EQUATION OF E_K-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>		
			<ci>E_K </ci>
			<apply>
          <times/>     			
				<ci> RToF </ci>
				<apply>
            <ln/>
					<apply>
              <divide/>
						<ci> K_o </ci>
						<ci> K_i </ci>
					</apply>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation E_K -->
	<!-- EQUATION OF I_bk-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>I_bk_SM </ci>
			<apply>
          <times/>
				<ci> G_max_bk </ci>
				<apply>
            <minus/>
					<ci> Vm_SM </ci>
					<ci> E_K </ci>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation I_bk -->
</component> <!-- closure of I_bk-->



<component name="xr1_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="xr1_inf_SM"/>
	<variable units="time_units" name="tau_xr1_SM"/>
	<variable units="dimensionless" public_interface="out" name="xr1_SM" initial_value="0.0"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_K"/>
	
	<!-- equation that calculates xr1-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> xr1_inf_SM </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> Vm_SM </ci>
								<cn cellml:units="voltage_units">27</cn>
							</apply>
							<cn cellml:units="voltage_units">-5.0</cn>
						</apply>
					</apply>
				</apply>
			</apply>		
		</apply>
      <!-- closure of equation that calculates xr1-->	
	</math>
	
	<!-- equation that calculates tau_xr1-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_xr1_SM </ci>
			<apply>
          <times/>
				<ci>T_correction_K</ci>
				<cn cellml:units="time_units">80</cn>
			</apply>	
		</apply>
      <!-- closure of equation that calculates xr1-->	
	</math>
	<!-- equation that calculates xr1-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> xr1_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> xr1_inf_SM</ci>
					<ci>xr1_SM</ci>
				</apply>
				<ci> tau_xr1_SM</ci>
			</apply>
		</apply>
      <!-- closure of equation d(xr1)/dt=(xr1_inf-xr1)/tau_xr1-->	
	</math> 
</component> <!-- closure of xr1-->

<component name="xr2_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="xr2_inf_SM"/>
	<variable units="time_units" name="tau_xr2_SM"/>
	<variable units="dimensionless" public_interface="out" name="xr2_SM" initial_value="0.82"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_K"/>
	
	<!-- equation that calculates xr2_inf-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> xr2_inf_SM </ci>
			<apply>
          <plus/>
				<cn cellml:units="dimensionless">0.2</cn>
				<apply>
            <divide/>
					<cn cellml:units="dimensionless">0.8</cn>
					<apply>
              <plus/>
						<cn cellml:units="dimensionless">1.0</cn>
						<apply>
                <exp/>
							<apply>
                  <divide/>
								<apply>
                    <plus/>
									<ci> Vm_SM </ci>
									<cn cellml:units="voltage_units">58</cn>
								</apply>
								<cn cellml:units="voltage_units">10</cn>
							</apply>
						</apply>
					</apply>
				</apply>	
			</apply>
		</apply>
      <!-- closure of equation that calculates xr2_inf-->	
	</math>
	
	<!-- equation that calculates tau_xr2-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_xr2_SM </ci>
			<apply>
          <times/>
				<ci>T_correction_K</ci>
				<apply>
            <plus/>
					<cn cellml:units="time_units">-707</cn>
					<apply>
              <times/>
						<cn cellml:units="time_units">1481</cn>
						<apply>
                <exp/>
							<apply>
                  <divide/>
								<apply>
                    <plus/>
									<ci> Vm_SM </ci>
									<cn cellml:units="voltage_units">36</cn>
								</apply>
								<cn cellml:units="voltage_units">95</cn>
							</apply>
						</apply>
					</apply>
				</apply>
			</apply>	
		</apply>
      <!-- closure of equation that calculates tau_xr2-->	
	</math>
	<!-- equation that calculates d_kr-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> xr2_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> xr2_inf_SM</ci>
					<ci>xr2_SM</ci>
				</apply>
				<ci> tau_xr2_SM</ci>
			</apply>
		</apply>
      <!-- closure of equation d(xr2)/dt=(xr2_inf-xr2)/tau_xr2-->	
	</math> 
</component> <!-- closure of xr2-->

<component name="I_kr_SM">
	
	<variable units="current_units" public_interface="out" name="I_kr_SM"/>
	<variable units="voltage_units" name="E_K"/>
	<variable units="conductance_units" name="G_max_kr_SM" initial_value="35"/>
	<variable units="voltage_units" public_interface="in" private_interface="out" name="Vm_SM"/>
	<variable units="dimensionless" private_interface="in" name="xr1_SM"/>
	<variable units="dimensionless" private_interface="in" name="xr2_SM"/>
	<variable units="voltage_units" public_interface="in" name="RToF"/>
	<variable units="millimolar" public_interface="in" name="K_o"/>
	<variable units="millimolar" public_interface="in" name="K_i"/>
	<variable units="dimensionless" public_interface="in" private_interface="out" name="T_correction_K"/>
	<variable units="time_units" public_interface="in" private_interface="out" name="time"/>
	
	<!-- EQUATION OF E_K-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>		
			<ci>E_K </ci>
			<apply>
          <times/>     			
				<ci> RToF </ci>
				<apply>
            <ln/>
					<apply>
              <divide/>
						<ci> K_o </ci>
						<ci> K_i </ci>
					</apply>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation E_K -->
	<!-- EQUATION OF I_kr-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>I_kr_SM </ci>
			<apply>
          <times/>     			
				<ci> G_max_kr_SM </ci>
				<apply>
            <times/>
					<ci> xr1_SM </ci>
					<ci> xr2_SM </ci>
				</apply>
				<apply>
            <minus/>
					<ci> Vm_SM </ci>
					<ci> E_K </ci>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation I_kr -->
</component> <!-- closure of I_kr-->

<component name="m_Na_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="m_inf_Na"/>
	<variable units="time_units" name="tau_m_Na"/>
	<variable units="dimensionless" public_interface="out" name="m_Na_SM" initial_value="0.005"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_Na"/>
	
	<!-- equation that calculates m_Na-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> m_inf_Na </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> Vm_SM </ci>
								<cn cellml:units="voltage_units">47</cn>
							</apply>
							<cn cellml:units="voltage_units">-4.8</cn>
						</apply>
					</apply>
				</apply>
			</apply>		
		</apply>
      <!-- closure of equation that calculates m_inf_Na-->	
	</math>
	
	<!-- equation that calculates tau_m_Na-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_m_Na </ci>
			<apply>
          <times/>
				<ci>T_correction_Na</ci>
				<apply>
            <plus/>
					<apply>
              <times/>
						<ci>Vm_SM</ci>
						<cn cellml:units="Inverse_Voltage_units">-0.017</cn>
						<cn cellml:units="time_units">1</cn>
					</apply>
					<cn cellml:units="time_units">0.440</cn>
				</apply>
			</apply>	
		</apply>
      <!-- closure of equation that calculates tau_m_Na-->	
	</math>
	<!-- equation that calculates m_Na-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> m_Na_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> m_inf_Na </ci>
					<ci>m_Na_SM</ci>
				</apply>
				<ci> tau_m_Na </ci>
			</apply>
		</apply>
      <!-- closure of equation d(m_Na)/dt=(m_inf_Na-m_Na)/tau_m_Na-->	
	</math> 
</component> <!-- closure of m_Na-->

<component name="h_Na_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="h_inf_Na"/>
	<variable units="time_units" name="tau_h_Na"/>
	<variable units="dimensionless" public_interface="out" name="h_Na_SM" initial_value="0.05787"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_Na"/>
	
	<!-- equation that calculates f_inf_Na-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> h_inf_Na </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> Vm_SM </ci>
								<cn cellml:units="voltage_units">78</cn>
							</apply>
							<cn cellml:units="voltage_units">3</cn>
						</apply>
					</apply>
				</apply>
			</apply>		
		</apply>
      <!-- closure of equation that calculates h_inf_Na-->	
	</math>
	
	<!-- equation that calculates tau_h_Na-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_h_Na </ci>
			<apply>
          <times/>
				<ci>T_correction_Na</ci>
				<apply>
            <plus/>
					<apply>
              <times/>
						<ci>Vm_SM</ci>
						<cn cellml:units="Inverse_Voltage_units">-0.25</cn>
						<cn cellml:units="time_units">1</cn>
					</apply>
					<cn cellml:units="time_units">5.50</cn>
				</apply>
			</apply>	
		</apply>
      <!-- closure of equation that calculates tau_h_Na-->	
	</math>
	<!-- equation that calculates f_Na-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> h_Na_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> h_inf_Na </ci>
					<ci>h_Na_SM</ci>
				</apply>
				<ci> tau_h_Na </ci>
			</apply>
		</apply>
      <!-- closure of equation d(h_Na)/dt=(h_inf_Na-f_Na)/tau_h_Na-->	
	</math> 
</component> <!-- closure of hq_Na-->

<component name="I_Na_SM">
	
	<variable units="current_units" public_interface="out" name="I_Na_SM"/>
	<variable units="voltage_units" name="E_Na"/>
	<variable units="conductance_units" name="G_max_Na_SM" initial_value="3"/>
	<variable units="voltage_units" public_interface="in" private_interface="out" name="Vm_SM"/>
	<variable units="dimensionless" private_interface="in" name="h_Na_SM"/>
	<variable units="dimensionless" private_interface="in" name="m_Na_SM"/>
	<variable units="voltage_units" public_interface="in" name="RToF"/>
	<variable units="millimolar" public_interface="in" name="Na_o"/>
	<variable units="millimolar" public_interface="in" name="Na_i"/>
	<variable units="dimensionless" public_interface="in" private_interface="out" name="T_correction_Na"/>
	<variable units="time_units" public_interface="in" private_interface="out" name="time"/>
	
	<!-- EQUATION OF E_Na-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>E_Na </ci>
			<apply>
          <times/>     			
				<ci> RToF </ci>
				<apply>
            <ln/>
					<apply>
              <divide/>
						<ci> Na_o </ci>
						<ci> Na_i </ci>
					</apply>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation E_Na -->
	<!-- EQUATION OF I_Na-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>I_Na_SM </ci>
			<apply>
          <times/>     			
				<ci> G_max_Na_SM </ci>
				<apply>
            <times/>
					<ci> h_Na_SM </ci>
					<ci> m_Na_SM </ci>
				</apply>
				<apply>
            <minus/>
					<ci> Vm_SM </ci>
					<ci> E_Na </ci>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation I_Na -->
</component> <!-- closure of I_Na-->

<component name="xa1_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="xa1_inf_SM"/>
	<variable units="time_units" name="tau_xa1_SM"/>
	<variable units="dimensionless" public_interface="out" name="xa1_SM" initial_value="0.00414"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_K"/>
	
	<!-- equation that calculates xa1_inf-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> xa1_inf_SM </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> Vm_SM </ci>
								<cn cellml:units="voltage_units">26.5</cn>
							</apply>
							<cn cellml:units="voltage_units">-7.9</cn>
						</apply>
					</apply>
				</apply>
			</apply>		
		</apply>
      <!-- closure of equation that calculates xa1_inf-->	
	</math>
	
	<!-- equation that calculates tau_xa1-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_xa1_SM </ci>
			<apply>
          <times/>
				<ci>T_correction_K</ci>
				<apply>
            <plus/>
					<cn cellml:units="time_units">31.8</cn>
					<apply>
              <times/>
						<cn cellml:units="time_units">175</cn>	
						<apply>
                <exp/>
							<apply>
                  <times/>
								<cn cellml:units="dimensionless">-0.5</cn>
								<apply>
                    <power/>
									<apply>
                      <divide/>
										<apply>
                        <plus/>
											<ci>Vm_SM</ci>
											<cn cellml:units="voltage_units">44.4</cn>
										</apply>
										<cn cellml:units="voltage_units">22.3</cn>
									</apply>
									<cn cellml:units="dimensionless">2</cn>
								</apply>
							</apply>
						</apply>
					</apply>
				</apply>
			</apply>	
		</apply>
      <!-- closure of equation that calculates tau_xa1-->	
	</math>
	<!-- equation that calculates xa1-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> xa1_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> xa1_inf_SM</ci>
					<ci>xa1_SM</ci>
				</apply>
				<ci> tau_xa1_SM</ci>
			</apply>
		</apply>
      <!-- closure of equation d(xa1)/dt=(xa1_inf-xa1)/tau_xa1-->	
	</math> 
</component> <!-- closure of xa1-->

<component name="xa2_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="xa2_inf_SM"/>
	<variable units="time_units" name="tau_xa2_SM"/>
	<variable units="dimensionless" public_interface="out" name="xa2_SM" initial_value="0.72"/>
	<variable units="time_units" public_interface="in" name="time"/>
	<variable units="dimensionless" public_interface="in" name="T_correction_K"/>
	
	<!-- equation that calculates xa2_inf-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> xa2_inf_SM </ci>
			<apply>
          <plus/>
				<cn cellml:units="dimensionless">0.1</cn>
				<apply>
            <divide/>
					<cn cellml:units="dimensionless">0.9</cn>
					<apply>
              <plus/>
						<cn cellml:units="dimensionless">1.0</cn>
						<apply>
                <exp/>
							<apply>
                  <divide/>
								<apply>
                    <plus/>
									<ci> Vm_SM </ci>
									<cn cellml:units="voltage_units">65</cn>
								</apply>
								<cn cellml:units="voltage_units">6.2</cn>
							</apply>
						</apply>
					</apply>
				</apply>	
			</apply>
		</apply>
      <!-- closure of equation that calculates xa2_inf-->	
	</math>
	
	<!-- equation that calculates tau_xa2-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_xa2_SM </ci>
			<apply>
          <times/>
				<ci>T_correction_K</ci>
				<cn cellml:units="time_units">90</cn>
			</apply>	
		</apply>
      <!-- closure of equation that calculates tau_xa2-->	
	</math>
	<!-- equation that calculates d_ka-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> xa2_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> xa2_inf_SM</ci>
					<ci>xa2_SM</ci>
				</apply>
				<ci> tau_xa2_SM</ci>
			</apply>
		</apply>
      <!-- closure of equation d(xa2)/dt=(xa2_inf-d_ka)/tau_xa2-->	
	</math> 
</component> <!-- closure of xa2-->

<component name="I_ka_SM">
	
	<variable units="current_units" public_interface="out" name="I_ka_SM"/>
	<variable units="voltage_units" name="E_K"/>
	<variable units="conductance_units" name="G_max_ka_SM" initial_value="9"/>
	<variable units="voltage_units" public_interface="in" private_interface="out" name="Vm_SM"/>
	<variable units="dimensionless" private_interface="in" name="xa1_SM"/>
	<variable units="dimensionless" private_interface="in" name="xa2_SM"/>
	<variable units="voltage_units" public_interface="in" name="RToF"/>
	<variable units="millimolar" public_interface="in" name="K_o"/>
	<variable units="millimolar" public_interface="in" name="K_i"/>
	<variable units="dimensionless" public_interface="in" private_interface="out" name="T_correction_K"/>
	<variable units="time_units" public_interface="in" private_interface="out" name="time"/>
	
	<!-- EQUATION OF E_K-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>		
			<ci>E_K </ci>
			<apply>
          <times/>     			
				<ci> RToF </ci>
				<apply>
            <ln/>
					<apply>
              <divide/>
						<ci> K_o </ci>
						<ci> K_i </ci>
					</apply>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation E_K -->
	<!-- EQUATION OF I_ka-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>I_ka_SM </ci>
			<apply>
          <times/>     			
				<ci> G_max_ka_SM </ci>
				<apply>
            <times/>
					<ci> xa1_SM </ci>
					<ci> xa2_SM </ci>
				</apply>
				<apply>
            <minus/>
					<ci> Vm_SM </ci>
					<ci> E_K </ci>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation I_ka -->
</component> <!-- closure of I_ka-->


<component name="m_NSCC_SM">
	<variable units="voltage_units" public_interface="in" name="Vm_SM"/>
	<variable units="dimensionless" name="m_inf_NSCC_SM"/>
	<variable units="time_units" name="tau_m_NSCC_SM"/>
	<variable units="dimensionless" public_interface="out" name="m_NSCC_SM" initial_value="0.0"/>
	<variable units="time_units" public_interface="in" name="time"/>
	
	<!-- equation that calculates m_inf_NSCC-->	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> m_inf_NSCC_SM </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> Vm_SM </ci>
								<cn cellml:units="voltage_units">25.0</cn>
							</apply>
							<cn cellml:units="voltage_units">-20</cn>
						</apply>
					</apply>
				</apply>
			</apply>
		</apply>
      <!-- closure of equation that calculates m_inf_NSCC-->	
	</math>
	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> tau_m_NSCC_SM </ci>
			<apply>
          <times/>
				<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>Vm_SM</ci>
									<cn cellml:units="voltage_units">66</cn>
								</apply>
								<cn cellml:units="voltage_units">-26</cn>
							</apply>
						</apply>
					</apply>
				</apply>
				<cn cellml:units="time_units">150</cn>
			</apply>	
		</apply>
      <!-- closure of equation that calculates tau_m_NSCC-->	
	</math>
	
	<!-- equation that calculates d_NSCC-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">	
		<apply>
        <eq/>	
			<apply>
          <diff/>
				<bvar>
            <ci>time </ci>
          </bvar>
				<ci> m_NSCC_SM </ci>
			</apply>
			<apply>
          <divide/>	
				<apply>
            <minus/>	
					<ci> m_inf_NSCC_SM </ci>
					<ci>m_NSCC_SM</ci>
				</apply>
				<ci> tau_m_NSCC_SM </ci>
			</apply>
		</apply>
      <!-- closure of equation d(m_NSCC)/dt=(m_inf_NSCC-m_NSCC)/tau_m_NSCC-->	
	</math> 
</component> <!-- closure of d_NSCC-->

<component name="I_NSCC_SM">
	
	<variable units="current_units" public_interface="out" name="I_NSCC_SM"/>
	<variable units="voltage_units" name="E_NSCC" initial_value="-28"/>
	<variable units="conductance_units" name="G_max_NSCC_SM" initial_value="50"/>
	<variable units="voltage_units" public_interface="in" private_interface="out" name="Vm_SM"/>
	<variable units="millimolar" public_interface="in" private_interface="out" name="Ca_i"/>
	<variable units="millimolar" public_interface="in" name="Ach"/>
	<variable units="dimensionless" private_interface="in" name="m_NSCC_SM"/>
	<variable units="dimensionless" name="f_ca_NSCC_SM"/>
	<variable units="dimensionless" name="rach_NSCC_SM"/>
	<variable units="voltage_units" public_interface="in" name="RToF"/>
	<variable units="millimolar" public_interface="in" name="Na_o"/>
	<variable units="millimolar" public_interface="in" name="Na_i"/>
	<variable units="millimolar" public_interface="in" name="K_o"/>
	<variable units="millimolar" public_interface="in" name="K_i"/>
	<variable units="time_units" public_interface="in" private_interface="out" name="time"/>
	

	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> f_ca_NSCC_SM </ci>
			<apply>
          <divide/>
				<cn cellml:units="dimensionless">1.0</cn>
				<apply>
            <plus/>
					<cn cellml:units="dimensionless">1.0</cn>
					<apply>
              <power/>
						<apply>
                <divide/>
							<ci>Ca_i</ci>
							<cn cellml:units="millimolar">0.0002</cn>
							
						</apply>
						<cn cellml:units="dimensionless">-4.0</cn>
					</apply>
				</apply>
			</apply>	
		</apply>
      <!-- closure of equation that calculates fca_NSCC-->	
	</math>
	
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
        <eq/>		
			<ci> rach_NSCC_SM </ci>
			<apply>
          <divide/>
				<cn cellml:units="dimensionless">1.0</cn>
				<apply>
            <plus/>
					<cn cellml:units="dimensionless">1.0</cn>
					<apply>
              <divide/>
						<cn cellml:units="millimolar">0.01</cn>
						<ci>Ach</ci>
					</apply>
				</apply>
			</apply>	
		</apply>
      <!-- closure of equation that calculates rach_NSCC-->	
	</math>
	
	<!-- EQUATION OF I_NSCC-->
	<math xmlns="http://www.w3.org/1998/Math/MathML">   	
		<apply>
        <eq/>		
			<ci>I_NSCC_SM </ci>
			<apply>
          <times/>     			
				<ci> G_max_NSCC_SM </ci>
				<ci> m_NSCC_SM </ci>
				<ci> f_ca_NSCC_SM </ci>
				<ci> rach_NSCC_SM </ci>
				<apply>
            <minus/>
					<ci> Vm_SM </ci>
					<ci> E_NSCC </ci>
				</apply>
			</apply>		
		</apply>	
	</math> <!-- closure ofEquation I_NSCC -->
</component> <!-- closure of I_NSCC-->


<!--Sodium-->
<connection>
	<map_components component_2="I_Na_SM" component_1="m_Na_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="m_Na_SM" variable_1="m_Na_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_Na" variable_1="T_correction_Na"/>
</connection>
<connection>
	<map_components component_2="I_Na_SM" component_1="h_Na_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="h_Na_SM" variable_1="h_Na_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_Na" variable_1="T_correction_Na"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="I_Na_SM"/>
	<map_variables variable_2="Na_i" variable_1="Na_i"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="I_Na_SM" variable_1="I_Na_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="RToF" variable_1="RToF"/>
	<map_variables variable_2="Na_o" variable_1="Na_o"/>
	<map_variables variable_2="T_correction_Na" variable_1="T_correction_Na"/>
</connection>

<!--Ltype-->
<connection>
	<map_components component_2="I_Ltype_SM" component_1="d_Ltype_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="d_Ltype_SM" variable_1="d_Ltype_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_Ca" variable_1="T_correction_Ca"/>
</connection>
<connection>
	<map_components component_2="I_Ltype_SM" component_1="f_Ltype_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="f_Ltype_SM" variable_1="f_Ltype_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_Ca" variable_1="T_correction_Ca"/>
</connection>
<connection>
	<map_components component_2="I_Ltype_SM" component_1="f_ca_Ltype_SM"/>
	<map_variables variable_2="f_ca_Ltype_SM" variable_1="f_ca_Ltype_SM"/>
	<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_Ca" variable_1="T_correction_Ca"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="I_Ltype_SM"/>
	<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="I_Ltype_SM" variable_1="I_Ltype_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="RToF" variable_1="RToF"/>
	<map_variables variable_2="Ca_o" variable_1="Ca_o"/>
	<map_variables variable_2="T_correction_Ca" variable_1="T_correction_Ca"/>
</connection>
	
<!--LVA-->
<connection>
	<map_components component_2="I_LVA_SM" component_1="d_LVA_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="d_LVA_SM" variable_1="d_LVA_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_Ca" variable_1="T_correction_Ca"/>
</connection>
<connection>
	<map_components component_2="I_LVA_SM" component_1="f_LVA_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="f_LVA_SM" variable_1="f_LVA_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_Ca" variable_1="T_correction_Ca"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="I_LVA_SM"/>
	<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="I_LVA_SM" variable_1="I_LVA_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="RToF" variable_1="RToF"/>
	<map_variables variable_2="Ca_o" variable_1="Ca_o"/>
	<map_variables variable_2="T_correction_Ca" variable_1="T_correction_Ca"/>
</connection>

<!--kr-->
<connection>
	<map_components component_2="I_kr_SM" component_1="xr1_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="xr1_SM" variable_1="xr1_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_K" variable_1="T_correction_K"/>
</connection>
<connection>
	<map_components component_2="I_kr_SM" component_1="xr2_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="xr2_SM" variable_1="xr2_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_K" variable_1="T_correction_K"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="I_kr_SM"/>
	<map_variables variable_2="K_i" variable_1="K_i"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="I_kr_SM" variable_1="I_kr_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="RToF" variable_1="RToF"/>
	<map_variables variable_2="K_o" variable_1="K_o"/>
	<map_variables variable_2="T_correction_K" variable_1="T_correction_K"/>
</connection>


<!--katype-->
<connection>
	<map_components component_2="I_ka_SM" component_1="xa1_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="xa1_SM" variable_1="xa1_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_K" variable_1="T_correction_K"/>
</connection>
<connection>
	<map_components component_2="I_ka_SM" component_1="xa2_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="xa2_SM" variable_1="xa2_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="T_correction_K" variable_1="T_correction_K"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="I_ka_SM"/>
	<map_variables variable_2="K_i" variable_1="K_i"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="I_ka_SM" variable_1="I_ka_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="RToF" variable_1="RToF"/>
	<map_variables variable_2="K_o" variable_1="K_o"/>
	<map_variables variable_2="T_correction_K" variable_1="T_correction_K"/>
</connection>
<!--BK-->
<connection>
	<map_components component_2="I_BK_SM" component_1="d_BK_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
	<map_variables variable_2="d_BK_SM" variable_1="d_BK_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="I_BK_SM"/>
	<map_variables variable_2="K_i" variable_1="K_i"/>
	<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="I_BK_SM" variable_1="I_BK_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="RToF" variable_1="RToF"/>
	<map_variables variable_2="K_o" variable_1="K_o"/>
	<map_variables variable_2="T_correction_BK" variable_1="T_correction_BK"/>
</connection>

<!--kb-->
<connection>
	<map_components component_2="SM_Membrane" component_1="I_bk_SM"/>
	<map_variables variable_2="K_i" variable_1="K_i"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="I_bk_SM" variable_1="I_bk_SM"/>
	<map_variables variable_2="RToF" variable_1="RToF"/>
	<map_variables variable_2="K_o" variable_1="K_o"/>
</connection>

<!--NSCC-->
<connection>
	<map_components component_2="I_NSCC_SM" component_1="m_NSCC_SM"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="m_NSCC_SM" variable_1="m_NSCC_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="I_NSCC_SM"/>
	<map_variables variable_2="K_i" variable_1="K_i"/>
	<map_variables variable_2="Na_i" variable_1="Na_i"/>
	<map_variables variable_2="Vm_SM" variable_1="Vm_SM"/>
	<map_variables variable_2="I_NSCC_SM" variable_1="I_NSCC_SM"/>
	<map_variables variable_2="time" variable_1="time"/>
	<map_variables variable_2="RToF" variable_1="RToF"/>
	<map_variables variable_2="K_o" variable_1="K_o"/>
	<map_variables variable_2="Na_o" variable_1="Na_o"/>
	<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
	<map_variables variable_2="Ach" variable_1="Ach"/>
</connection>

<!--Jsr-->
<connection>
	<map_components component_2="SM_Membrane" component_1="J_CaSR_SM"/>
	<map_variables variable_2="J_CaSR_SM" variable_1="J_CaSR_SM"/>
	<map_variables variable_2="Ca_i" variable_1="Ca_i"/>
</connection>

<connection>
	<map_components component_2="Environment" component_1="SM_Membrane"/>
	<map_variables variable_2="Ca_o" variable_1="Ca_o"/>
	<map_variables variable_2="Na_o" variable_1="Na_o"/>
	<map_variables variable_2="K_o" variable_1="K_o"/>
	<map_variables variable_2="Cl_o" variable_1="Cl_o"/>
	<map_variables variable_2="RToF" variable_1="RToF"/>
	<map_variables variable_2="FoRT" variable_1="FoRT"/>
	<map_variables variable_2="F" variable_1="F"/>
	<map_variables variable_2="T_correction_Ca" variable_1="T_correction_Ca"/>
	<map_variables variable_2="T_correction_Na" variable_1="T_correction_Na"/>
	<map_variables variable_2="T_correction_K" variable_1="T_correction_K"/>
	<map_variables variable_2="T_correction_BK" variable_1="T_correction_BK"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="Time"/>
	<map_variables variable_2="time" variable_1="time"/>
</connection>

<connection>
	<map_components component_2="I_stim" component_1="Gap_junction"/>
	<map_variables variable_2="Gcouple" variable_1="Gcouple"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="I_stim"/>
	<map_variables variable_2="I_stim" variable_1="I_stim"/>
</connection>
<connection>
        <map_components component_2="Time" component_1="I_stim"/>
        <map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
	<map_components component_2="SM_Membrane" component_1="Neural_input"/>
	<map_variables variable_2="Ach" variable_1="Ach"/>
</connection>

<group>
    <relationship_ref relationship="encapsulation"/>
	<component_ref component="SM_Membrane">
		<component_ref component="I_Na_SM">
			<component_ref component="m_Na_SM"/>
			<component_ref component="h_Na_SM"/>
		</component_ref>
		<component_ref component="I_Ltype_SM">
			<component_ref component="d_Ltype_SM"/>
			<component_ref component="f_Ltype_SM"/>
			<component_ref component="f_ca_Ltype_SM"/>
		</component_ref>
		<component_ref component="I_LVA_SM">
			<component_ref component="d_LVA_SM"/>
			<component_ref component="f_LVA_SM"/>
		</component_ref>
		<component_ref component="I_BK_SM">
			<component_ref component="d_BK_SM"/>
		</component_ref>
		<component_ref component="I_kr_SM">
			<component_ref component="xr1_SM"/>
			<component_ref component="xr2_SM"/>
		</component_ref>
		<component_ref component="I_ka_SM">
			<component_ref component="xa1_SM"/>
			<component_ref component="xa2_SM"/>
		</component_ref>
		<component_ref component="I_bk_SM">
		</component_ref>
		<component_ref component="I_NSCC_SM">
			<component_ref component="m_NSCC_SM"/>
		</component_ref>
		<component_ref component="J_CaSR_SM">
		</component_ref>
	</component_ref>	
	
</group>
<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#0977f4f9-8049-4ecd-9f95-6023c52b434e">
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>gastric</rdf:li>
    <rdf:li>smooth muscle</rdf:li>
    <rdf:li>ICC</rdf:li>
    <rdf:li>interstitial cells of Cajal</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#9c73ca47-a776-4599-82ae-23a4b1dd802a">
    <rdf:li rdf:resource="rdf:#67dd54c1-ccb7-4726-b4bd-c0232ed28ce7"/>
    <rdf:li rdf:resource="rdf:#be6ff2ad-5443-4700-9c60-9433a6bfffbe"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#d06b5f7f-e46a-401a-bbf2-d69424890530">
    <ns8:y-variable rdf:resource="rdf:#e1d2227d-8b8a-4b23-b694-e52a6053d80b"/>
    <ns8:colour>#ff0000</ns8:colour>
    <ns8:line-type>1</ns8:line-type>
    <ns8:x-variable rdf:resource="rdf:#f50a3400-9738-4f96-9b23-65427674645e"/>
    <ns8:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
  </rdf:Description>
  <rdf:Description rdf:about="#AP_IBK">
    <ns8:traces rdf:resource="rdf:#babad30c-9b84-488f-bad2-832dd8384a08"/>
    <ns8:background-colour>#ffffff</ns8:background-colour>
    <ns8:y-label>Current (pA)</ns8:y-label>
    <ns8:x-label>Time (ms)</ns8:x-label>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>National University of Singapore</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#465ac8c0-65aa-4b36-9ce5-19496bf4ee41"/>
    <dcterms:created rdf:resource="rdf:#3a0a205d-8b8c-45aa-9515-aa1d6a9313a6"/>
    <dc:creator rdf:resource="rdf:#a3458722-cda9-4cc2-a526-d3a0d823bce4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bd83516b-83b5-4054-add0-4e7296a7eac7">
    <ns8:y-variable rdf:resource="rdf:#edd91ee3-b525-45f0-a51c-b5a844eee24a"/>
    <ns8:colour>#ff0000</ns8:colour>
    <ns8:line-type>1</ns8:line-type>
    <ns8:x-variable rdf:resource="rdf:#cfc4353a-3d53-48b2-b579-63ba5efe4391"/>
    <ns8:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
  </rdf:Description>
  <rdf:Description rdf:about="#ICC_stimulus">
    <cmeta:comment rdf:resource="rdf:#d619ded7-61bf-4782-a969-278737f1b968"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c9fc9463-9563-474b-bea1-b5f01eb24231">
    <dc:title>Annals of Biomedical Engineering</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3b423a0f-79ec-4a38-9525-0b100af214a3">
    <dcterms:W3CDTF>2007-09-21</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5780474b-949d-4b3c-a511-7aee58f16c73">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#time"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#be6ff2ad-5443-4700-9c60-9433a6bfffbe">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#05f35ef9-1276-465a-a7f3-80be404f54bf"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#40499aa9-5125-4199-98f5-9c3abf0c658c">
    <bqs:Pubmed_id>17486452</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#68dd4ee3-e8c0-488e-a88b-f737ad1a390f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8b39b093-6388-46da-b00c-61c7135fabda">
    <ns8:y-variable rdf:resource="rdf:#5c30c764-0404-49ca-ba7b-e2d4d8bcb216"/>
    <ns8:colour>#ff0000</ns8:colour>
    <ns8:line-type>1</ns8:line-type>
    <ns8:x-variable rdf:resource="rdf:#5780474b-949d-4b3c-a511-7aee58f16c73"/>
    <ns8:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cfc4353a-3d53-48b2-b579-63ba5efe4391">
    <ns8:minimum>19.5e3</ns8:minimum>
    <ns8:variable rdf:resource="#time"/>
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:maximum>32.0e3</ns8:maximum>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3744fa51-fd80-4fcc-8a60-cd445689464d">
    <rdf:rest rdf:resource="http://www.w3.org/2002/07/owl#Nothing"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#bd83516b-83b5-4054-add0-4e7296a7eac7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#891064af-f187-4e4f-be0e-5a835b73d781">
    <ns8:y-variable rdf:resource="rdf:#e8fa17e6-137b-4847-9268-d53dab888759"/>
    <ns8:colour>#ff0000</ns8:colour>
    <ns8:line-type>1</ns8:line-type>
    <ns8:x-variable rdf:resource="rdf:#e30f3c77-0c82-4882-b0f7-aeb636eb7978"/>
    <ns8:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3ff557bf-3531-4dd7-88c9-0f73aa64d4c6">
    <rdf:rest rdf:resource="http://www.w3.org/2002/07/owl#Nothing"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#d06b5f7f-e46a-401a-bbf2-d69424890530"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#60065d76-13e7-4b34-9218-f877acc48268">
    <dc:creator rdf:resource="rdf:#821c4f4a-1adb-4e61-a60e-1169a5135605"/>
    <rdf:value>This is a CellML description of all the ion channels and intracellular processes contained  in the Corrias and Buist 2007 quantitative model of gastric smooth muscle cellular activation. This encoding of the model includes a stimulus current derived from a psuedo ICC action potential.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#31a19f66-726c-4efd-8247-2d2c2e3c9af6">
    <rdf:rest rdf:resource="http://www.w3.org/2002/07/owl#Nothing"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#8b39b093-6388-46da-b00c-61c7135fabda"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b248ac63-5eaa-4cef-97f3-79445cf5a456">
    <vCard:FN>Alberto Corrias</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0f25d5f0-64c7-409f-8fae-2dae5a22ff6e">
    <vCard:Given>Alberto</vCard:Given>
    <vCard:Family>Corrias</vCard:Family>
    <vCard:Other/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#79382410-660e-4232-9bd7-e695d271d4fd">
    <ns8:y-variable rdf:resource="rdf:#780d4684-55b5-4240-8366-d82fb0fd50ef"/>
    <ns8:colour>#ff0000</ns8:colour>
    <ns8:line-type>1</ns8:line-type>
    <ns8:x-variable rdf:resource="rdf:#494d5d0b-e0ac-4a5c-bf42-95498b143858"/>
    <ns8:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#98c6dafd-7bfd-4ac1-a333-6cc503859bab">
    <dcterms:modified rdf:resource="rdf:#3b423a0f-79ec-4a38-9525-0b100af214a3"/>
    <rdf:value>
				The value of Q10 for Potassium channels (Q10K) was corrected from the value in the published paper (1.5). The correct value to be used is 1.365. 
			</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#b248ac63-5eaa-4cef-97f3-79445cf5a456"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e8fa17e6-137b-4847-9268-d53dab888759">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#IBK"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#05f35ef9-1276-465a-a7f3-80be404f54bf">
    <vCard:Given>Martin</vCard:Given>
    <vCard:Family>Buist</vCard:Family>
    <vCard:Other>L</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fdb520dd-e7a0-4c28-a63e-06868cc155dc">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#IKr"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c9b76ea-df8c-4690-8e99-15dd239a053e">
    <rdf:rest rdf:resource="http://www.w3.org/2002/07/owl#Nothing"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#79382410-660e-4232-9bd7-e695d271d4fd"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e30f3c77-0c82-4882-b0f7-aeb636eb7978">
    <ns8:minimum>19.5e3</ns8:minimum>
    <ns8:variable rdf:resource="#time"/>
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:maximum>32.0e3</ns8:maximum>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3ae7c3df-235b-46e0-a611-f938f3663fa9">
    <rdf:rest rdf:resource="http://www.w3.org/2002/07/owl#Nothing"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#786fa03d-6007-41c2-9d79-6ea11cfbc7cc"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#edd91ee3-b525-45f0-a51c-b5a844eee24a">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#Cai"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a813d857-f369-4b3a-8538-f126229b923d">
    <rdf:rest rdf:resource="rdf:#3ae7c3df-235b-46e0-a611-f938f3663fa9"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#e7cbda8f-ce08-44ac-a842-71a5839c1ba6"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b64687f1-df70-45a3-bde2-995560fe9676">
    <ns8:y-variable rdf:resource="rdf:#fdb520dd-e7a0-4c28-a63e-06868cc155dc"/>
    <ns8:colour>#ff0000</ns8:colour>
    <ns8:line-type>1</ns8:line-type>
    <ns8:x-variable rdf:resource="rdf:#2ae59ac7-7f72-4529-8e91-dad48220e767"/>
    <ns8:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#363b3051-197b-4921-9044-b647d3962d17">
    <dcterms:modified rdf:resource="rdf:#95b01ee3-8735-4157-9f1a-4e9709d4ceec"/>
    <rdf:value>
                          Tidying up the electrical stimulus protocol and adding simulation and graphing metadata for my CellMLSimulator software.
                        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#66a104b0-f5d1-4b94-ba92-a6559d3c47a5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d619ded7-61bf-4782-a969-278737f1b968">
    <rdf:value>
                Here we define a periodic stimulus based on a psuedo ICC action potential. This stimulus protocol is used to test the model under a realistic stimulus protocol, but ideally such a stimulus would originate directly from an ICC cellular electrophysiology model. This stimulus protocol simply applies the psuedo ICC stimulus a fixed number of times.
              </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#67dd54c1-ccb7-4726-b4bd-c0232ed28ce7">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#f258a9e7-dd94-4003-be62-c6c8b11e928c"/>
  </rdf:Description>
  <rdf:Description rdf:about="#AP_INa">
    <ns8:traces rdf:resource="rdf:#3ff557bf-3531-4dd7-88c9-0f73aa64d4c6"/>
    <ns8:background-colour>#ffffff</ns8:background-colour>
    <ns8:y-label>Current (pA)</ns8:y-label>
    <ns8:x-label>Time (ms)</ns8:x-label>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a30c58cf-a734-420d-bf6a-07e20f2d6402">
    <dcterms:W3CDTF>2007</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#71181681-48cc-4661-8404-4875b11baed7">
    <vCard:Orgname>National University of Singapore</vCard:Orgname>
    <vCard:Orgunit>Division of BioEngineering</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f50a3400-9738-4f96-9b23-65427674645e">
    <ns8:minimum>19.5e3</ns8:minimum>
    <ns8:variable rdf:resource="#time"/>
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:maximum>32.0e3</ns8:maximum>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c30c764-0404-49ca-ba7b-e2d4d8bcb216">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#Vm"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#494d5d0b-e0ac-4a5c-bf42-95498b143858">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#time"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3a0a205d-8b8c-45aa-9515-aa1d6a9313a6">
    <dcterms:W3CDTF>2007-09-10T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#AP_Cai">
    <ns8:traces rdf:resource="rdf:#3744fa51-fd80-4fcc-8a60-cd445689464d"/>
    <ns8:background-colour>#ffffff</ns8:background-colour>
    <ns8:y-label>[Ca]i (mM)</ns8:y-label>
    <ns8:x-label>Time (ms)</ns8:x-label>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#16085fa2-8288-44d2-899c-dd4d67643b3f">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="#simulation">
    <ns7:iterationMethod>newton</ns7:iterationMethod>
    <ns7:linearSolver>dense</ns7:linearSolver>
    <ns7:boundIntervals rdf:resource="rdf:#587c163c-1e64-40f7-96da-11ec185b87c1"/>
    <ns7:multistepMethod>bdf</ns7:multistepMethod>
    <ns7:simulationName>SM with psuedo ICC stimulus</ns7:simulationName>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e6938a60-dfbe-4b05-96ac-df6a52e59dc9">
    <ns8:minimum>19.5e3</ns8:minimum>
    <ns8:variable rdf:resource="#time"/>
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:maximum>32.0e3</ns8:maximum>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#786fa03d-6007-41c2-9d79-6ea11cfbc7cc">
    <ns8:line-type>1</ns8:line-type>
    <ns8:x-variable rdf:resource="rdf:#e6938a60-dfbe-4b05-96ac-df6a52e59dc9"/>
    <ns8:label>ILVA</ns8:label>
    <ns8:colour>#0000ff</ns8:colour>
    <ns8:y-variable rdf:resource="rdf:#1400a0f5-06f8-4dc5-9c72-f5b80e0bade6"/>
    <ns8:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
  </rdf:Description>
  <rdf:Description rdf:about="#AP_IKr">
    <ns8:traces rdf:resource="rdf:#6d23387e-153f-40ac-bd1a-c88d838118c1"/>
    <ns8:background-colour>#ffffff</ns8:background-colour>
    <ns8:y-label>Current (pA)</ns8:y-label>
    <ns8:x-label>Time (ms)</ns8:x-label>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a3458722-cda9-4cc2-a526-d3a0d823bce4">
    <vCard:ORG rdf:resource="rdf:#71181681-48cc-4661-8404-4875b11baed7"/>
    <vCard:EMAIL rdf:resource="rdf:#6b908c78-8a41-4810-adef-6f306693eda3"/>
    <vCard:N rdf:resource="rdf:#0f25d5f0-64c7-409f-8fae-2dae5a22ff6e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#465ac8c0-65aa-4b36-9ce5-19496bf4ee41">
    <dc:creator rdf:resource="rdf:#16085fa2-8288-44d2-899c-dd4d67643b3f"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1c47d531-2935-4f3f-8369-1819b231c45d">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#IKA"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2ae59ac7-7f72-4529-8e91-dad48220e767">
    <ns8:minimum>19.5e3</ns8:minimum>
    <ns8:variable rdf:resource="#time"/>
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:maximum>32.0e3</ns8:maximum>
  </rdf:Description>
  <rdf:Description rdf:about="#membrane_potential_graph">
    <ns8:traces rdf:resource="rdf:#31a19f66-726c-4efd-8247-2d2c2e3c9af6"/>
    <ns8:background-colour>#ffffff</ns8:background-colour>
    <ns8:y-label>Membrane potential (mV)</ns8:y-label>
    <ns8:x-label>Time (ms)</ns8:x-label>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#66a104b0-f5d1-4b94-ba92-a6559d3c47a5">
    <vCard:FN>David Nickerson</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d024d030-6aa0-4791-a460-5601273d61a2">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#0977f4f9-8049-4ecd-9f95-6023c52b434e"/>
  </rdf:Description>
  <rdf:Description rdf:about="#AP_IKA">
    <ns8:traces rdf:resource="rdf:#c765e9a6-0749-49e7-9b69-0921e5311b32"/>
    <ns8:background-colour>#ffffff</ns8:background-colour>
    <ns8:y-label>Current (pA)</ns8:y-label>
    <ns8:x-label>Time (ms)</ns8:x-label>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4aa260cf-b0c1-4f79-9598-e0599a79e73a">
    <ns8:y-variable rdf:resource="rdf:#1c47d531-2935-4f3f-8369-1819b231c45d"/>
    <ns8:colour>#ff0000</ns8:colour>
    <ns8:line-type>1</ns8:line-type>
    <ns8:x-variable rdf:resource="rdf:#50c656cb-4f1c-47fb-ad91-da6a4202d2dc"/>
    <ns8:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
  </rdf:Description>
  <rdf:Description rdf:about="#corrias_model_2007">
    <ns7:simulation rdf:resource="#simulation"/>
    <dc:title>
			Corrias and Buist's 2007 quantitative model of gastric smooth muscle cellular electrophysiology.
		</dc:title>
    <ns8:graph rdf:resource="#AP_Cai"/>
    <ns8:graph rdf:resource="#AP_IBK"/>
    <ns8:graph rdf:resource="#AP_ICa"/>
    <ns8:graph rdf:resource="#AP_IKA"/>
    <ns8:graph rdf:resource="#AP_IKr"/>
    <ns8:graph rdf:resource="#AP_INa"/>
    <ns8:graph rdf:resource="#membrane_potential_graph"/>
    <ns8:graph rdf:resource="#stimulus_current_graph"/>
    <bqs:reference rdf:resource="rdf:#40499aa9-5125-4199-98f5-9c3abf0c658c"/>
    <bqs:reference rdf:resource="rdf:#69ba8aa1-3ba1-44f1-bf62-6b86ec7125c9"/>
    <cmeta:comment rdf:resource="rdf:#60065d76-13e7-4b34-9218-f877acc48268"/>
    <cmeta:modification rdf:resource="rdf:#363b3051-197b-4921-9044-b647d3962d17"/>
    <cmeta:modification rdf:resource="rdf:#98c6dafd-7bfd-4ac1-a333-6cc503859bab"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#95b01ee3-8735-4157-9f1a-4e9709d4ceec">
    <dcterms:W3CDTF>2007-09-20</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fff3272f-9f07-404a-a8af-8bef0c1f1e54">
    <ns7:boundVariable rdf:resource="#time"/>
    <ns7:tabulationStepSize>
                                10.0
                              </ns7:tabulationStepSize>
    <ns7:endingValue>
                                60.0e3
                              </ns7:endingValue>
    <ns7:maximumStepSize>
                                1
                              </ns7:maximumStepSize>
    <ns7:startingValue>
                                0.0
                              </ns7:startingValue>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2da2986e-14cd-4fec-94d1-bd392592cbfc">
    <ns8:minimum>19.5e3</ns8:minimum>
    <ns8:variable rdf:resource="#time"/>
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:maximum>32.0e3</ns8:maximum>
  </rdf:Description>
  <rdf:Description rdf:about="#AP_ICa">
    <ns8:traces rdf:resource="rdf:#a813d857-f369-4b3a-8538-f126229b923d"/>
    <ns8:background-colour>#ffffff</ns8:background-colour>
    <ns8:y-label>Current (pA)</ns8:y-label>
    <ns8:x-label>Time (ms)</ns8:x-label>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1400a0f5-06f8-4dc5-9c72-f5b80e0bade6">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#ILVA"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#50c656cb-4f1c-47fb-ad91-da6a4202d2dc">
    <ns8:minimum>19.5e3</ns8:minimum>
    <ns8:variable rdf:resource="#time"/>
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:maximum>32.0e3</ns8:maximum>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#69ba8aa1-3ba1-44f1-bf62-6b86ec7125c9">
    <dc:subject rdf:resource="rdf:#d024d030-6aa0-4791-a460-5601273d61a2"/>
  </rdf:Description>
  <rdf:Description rdf:about="#stimulus_current_graph">
    <ns8:traces rdf:resource="rdf:#5c9b76ea-df8c-4690-8e99-15dd239a053e"/>
    <ns8:background-colour>#ffffff</ns8:background-colour>
    <ns8:y-label>I_stim (pA)</ns8:y-label>
    <ns8:x-label>Time (ms)</ns8:x-label>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6d23387e-153f-40ac-bd1a-c88d838118c1">
    <rdf:rest rdf:resource="http://www.w3.org/2002/07/owl#Nothing"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#b64687f1-df70-45a3-bde2-995560fe9676"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e1d2227d-8b8a-4b23-b694-e52a6053d80b">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#INa"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#821c4f4a-1adb-4e61-a60e-1169a5135605">
    <vCard:FN>Alberto Corrias</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#780d4684-55b5-4240-8366-d82fb0fd50ef">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#I_stim"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#babad30c-9b84-488f-bad2-832dd8384a08">
    <rdf:rest rdf:resource="http://www.w3.org/2002/07/owl#Nothing"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#891064af-f187-4e4f-be0e-5a835b73d781"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#587c163c-1e64-40f7-96da-11ec185b87c1">
    <rdf:rest rdf:resource="http://www.w3.org/2002/07/owl#Nothing"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#fff3272f-9f07-404a-a8af-8bef0c1f1e54"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c765e9a6-0749-49e7-9b69-0921e5311b32">
    <rdf:rest rdf:resource="http://www.w3.org/2002/07/owl#Nothing"/>
    <rdf:type rdf:resource="http://www.w3.org/1999/02/22-rdf-syntax-ns#List"/>
    <rdf:first rdf:resource="rdf:#4aa260cf-b0c1-4f79-9598-e0599a79e73a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f258a9e7-dd94-4003-be62-c6c8b11e928c">
    <vCard:Given>Alberto</vCard:Given>
    <vCard:Family>Corrias</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6bcafcd5-6f07-4000-99c9-08dbcec29294">
    <ns7:simulation rdf:resource="#simulation"/>
    <ns8:variable rdf:resource="#ICaL"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e7cbda8f-ce08-44ac-a842-71a5839c1ba6">
    <ns8:line-type>1</ns8:line-type>
    <ns8:x-variable rdf:resource="rdf:#2da2986e-14cd-4fec-94d1-bd392592cbfc"/>
    <ns8:label>ICaL</ns8:label>
    <ns8:colour>#ff0000</ns8:colour>
    <ns8:y-variable rdf:resource="rdf:#6bcafcd5-6f07-4000-99c9-08dbcec29294"/>
    <ns8:type rdf:resource="http://www.cellml.org/metadata/graphs/1.0#line"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6b908c78-8a41-4810-adef-6f306693eda3">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>albertocorrias@nus.edu.sg</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#68dd4ee3-e8c0-488e-a88b-f737ad1a390f">
    <dc:creator rdf:resource="rdf:#9c73ca47-a776-4599-82ae-23a4b1dd802a"/>
    <dc:title>A quantitative model of gastric smooth muscle cellular activation</dc:title>
    <bqs:volume>35</bqs:volume>
    <bqs:first_page>1595</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#c9fc9463-9563-474b-bea1-b5f01eb24231"/>
    <dcterms:issued rdf:resource="rdf:#a30c58cf-a734-420d-bf6a-07e20f2d6402"/>
    <bqs:last_page>1607</bqs:last_page>
  </rdf:Description>
</rdf:RDF>
</model>