Rendering of the source text

<?xml version="1.0"?>
<!--
This CellML file was generated on 17/12/2009 at 9:56:26 at a.m. using:

COR (0.9.31.1333)
Copyright 2002-2009 Dr Alan Garny
http://cor.physiol.ox.ac.uk/ - cor@physiol.ox.ac.uk

CellML 1.0 was used to generate this model
http://www.cellml.org/
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xml:base="" cmeta:id="mears_1997" name="mears_1997">
	
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Evidence That Calcium Release-activated Current Mediates the Biphasic Electrical Activity of Mouse Pancreatic </title>
				<author>
					<firstname>Tessa</firstname>
					<surname>Paris</surname>
					<affiliation>
						<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>This CellML model runs in both OpenCell and COR to produce simulation graphs which are similar to those in the original published paper, namely figures 2 and 6.</para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>ABSTRACT:  The electrical response of pancreatic beta-cells to step increases in glucose concentration is biphasic, consisting of a prolonged depolarization with action potentials (Phase 1) followed by membrane potential oscillations known as bursts. We have proposed that the Phase 1 response results from the combined depolarizing influences of potassium channel closure and an inward, nonselective cation current (ICRAN) that activates as intracellular calcium stores empty during exposure to basal glucose (Bertram et al., 1995). The stores refill during Phase 1, deactivating ICRAN and allowing steady-state bursting to commence. We support this hypothesis with additional simulations and experimental results indicating that Phase 1 duration is sensitive to the filling state of intracellular calcium stores. First, the duration of the Phase 1 transient increases with duration of prior exposure to basal (2.8 mM) glucose, reflecting the increased time required to fill calcium stores that have been emptying for longer periods. Second, Phase 1 duration is reduced when islets are exposed to elevated K+ to refill calcium stores in the presence of basal glucose. Third, when extracellular calcium is removed during the basal glucose exposure to reduce calcium influx into the stores, Phase 1 duration increases. Finally, no Phase 1 is observed following hyperpolarization of the beta-cell membrane with diazoxide in the continued presence of 11 mm glucose, a condition in which intracellular calcium stores remain full. Application of carbachol to empty calcium stores during basal glucose exposure did not increase Phase 1 duration as the model predicts. Despite this discrepancy, the good agreement between most of the experimental results and the model predictions provides evidence that a calcium release-activated current mediates the Phase 1 electrical response of the pancreatic beta-cell.</para>
				<para>
The original paper reference is cited below:
</para>
				<para>Evidence That Calcium Release-activated Current Mediates the Biphasic Electrical Activity of Mouse Pancreatic <emphasis>Journal name</emphasis>, edition, firstpage-lastpage.  <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/9002424">PubMed ID: 9002424</ulink></para>
				
				<informalfigure>
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>Diagram</title>
							</objectinfo>
							<imagedata fileref="mears_1997.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the model.</caption>
				</informalfigure>
			</sect1>
		</article>
	</documentation>
	
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">

		<rdf:Description rdf:about="">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
				<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<vCard:Family>Paris</vCard:Family>
					<vCard:Given>Tessa</vCard:Given>
					<vCard:Other>Hazel</vCard:Other>
				</vCard:N>
				<vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<rdf:value>tpar054@aucklanduni.ac.nz</rdf:value>
					<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
				</vCard:EMAIL>
				<vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<vCard:Orgname>The University of Auckland</vCard:Orgname>
					<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
				</vCard:ORG>
			</dc:creator>
			<dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
				<dcterms:W3CDTF>2009-12-16</dcterms:W3CDTF>
			</dcterms:created>
		</rdf:Description>
		
		<rdf:Description rdf:about="#mears_1997">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Evidence That Calcium Release-activated Current Mediates the Biphasic Electrical Activity of Mouse Pancreatic </dc:title>
			
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
            <dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
            <bqs:subject_type>keyword</bqs:subject_type>
            <rdf:value>
            <rdf:Bag>
              <rdf:li>Islet of Langerhans</rdf:li>
			  <rdf:li>pancreatic beta cells</rdf:li>
			  <rdf:li>calcium dynamics</rdf:li>
			  <rdf:li>electrophysiology</rdf:li>
            </rdf:Bag>
            </rdf:value>
            </dc:subject>
            </bqs:reference>
			
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
        <bqs:Pubmed_id>9002424</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/">
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Mears</vCard:Family>
                    <vCard:Given>D</vCard:Given>
			      </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Sheppard</vCard:Family>
                    <vCard:Given>N</vCard:Given>
					<vCard:Other>F</vCard:Other>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Atwater</vCard:Family>
                    <vCard:Given>I</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Rojas</vCard:Family>
                    <vCard:Given>E</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Bertram</vCard:Family>
                    <vCard:Given>R</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
                    <vCard:Family>Sherman</vCard:Family>
                    <vCard:Given>A</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
            Evidence That Calcium Release-activated Current Mediates the Biphasic Electrical Activity of Mouse Pancreatic Beta-Cells
          </dc:title>
          <dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
            <dcterms:W3CDTF>1997-01-01</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Journal of Membrane Biology</dc:title>
          </bqs:Journal>
		  <bqs:volume>155</bqs:volume>
          <bqs:first_page>47</bqs:first_page>
          <bqs:last_page>59</bqs:last_page>
        </bqs:JournalArticle>
      </bqs:reference>
    </rdf:Description>
  </rdf:RDF>
  
  
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"><rdf:Description rdf:about="#mears_1997"><simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n1"/></simulation></rdf:Description><rdf:Description rdf:nodeID="n2"><rdf:rest><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/></rdf:rest><rdf:first><rdf:Description rdf:nodeID="n3"/></rdf:first></rdf:Description><rdf:Description rdf:nodeID="n3"><endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">60000</endingValue><pointDensity xmlns="http://www.cellml.org/metadata/simulation/1.0#nonstandard-" xml:lang="en">180000</pointDensity></rdf:Description><rdf:Description rdf:nodeID="n1"><boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n2"/></boundIntervals></rdf:Description></rdf:RDF>   <units name="millisecond">
      <unit prefix="milli" units="second"/>
   </units>
   <units name="per_millisecond">
      <unit exponent="-1" units="millisecond"/>
   </units>
   <units name="micromolar_per_millisecond">
      <unit units="micromolar"/>
      <unit exponent="-1" units="millisecond"/>
   </units>
   <units name="micromolar_per_millisecond_picoA">
      <unit units="micromolar"/>
      <unit exponent="-1" units="millisecond"/>
      <unit exponent="-1" units="picoA"/>
   </units>
   <units name="millivolt">
      <unit prefix="milli" units="volt"/>
   </units>
   <units name="micromolar">
      <unit prefix="micro" units="mole"/>
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="micrometer_3">
      <unit exponent="3" prefix="micro" units="meter"/>
   </units>
   <units name="picoS">
      <unit prefix="pico" units="siemens"/>
   </units>
   <units name="picoA">
      <unit prefix="pico" units="ampere"/>
   </units>
   <units name="micromolar_per_picoA">
      <unit units="micromolar"/>
      <unit exponent="-1" prefix="pico" units="ampere"/>
   </units>
   <units name="femtoF">
      <unit prefix="femto" units="farad"/>
   </units>
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="millisecond"/>
   </component>
   <component name="membrane">
      <variable cmeta:id="membrane_V" initial_value="-61" name="V" public_interface="out" units="millivolt"/>
      <variable initial_value="6158" name="Cm" units="femtoF"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="i_K" public_interface="in" units="picoA"/>
      <variable name="i_K_Ca" public_interface="in" units="picoA"/>
      <variable name="i_K_ATP" public_interface="in" units="picoA"/>
      <variable name="i_CRAC" public_interface="in" units="picoA"/>
      <variable name="i_Ca" public_interface="in" units="picoA"/>
      <variable name="i_leak" public_interface="in" units="picoA"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>i_Ca</ci>
                     <ci>i_K</ci>
                     <ci>i_K_ATP</ci>
                     <ci>i_K_Ca</ci>
                     <ci>i_CRAC</ci>
                     <ci>i_leak</ci>
                  </apply>
               </apply>
               <ci>Cm</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="K_current">
      <variable name="time" private_interface="out" public_interface="in" units="millisecond"/>
      <variable cmeta:id="K_current_i_K" name="i_K" public_interface="out" units="picoA"/>
      <variable initial_value="-70" name="V_K" public_interface="out" units="millivolt"/>
      <variable initial_value="3900" name="g_K" units="picoS"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="n" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K</ci>
            <apply>
               <times/>
               <ci>g_K</ci>
               <ci>n</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="K_channel_n_gate">
      <variable initial_value="0.0005" name="n" public_interface="out" units="dimensionless"/>
      <variable name="n_infinity" units="dimensionless"/>
      <variable name="tau_n" units="millisecond"/>
      <variable initial_value="-15" name="Vn" units="millivolt"/>
      <variable initial_value="6" name="Sn" units="millivolt"/>
      <variable initial_value="1.85" name="lambda_n" units="dimensionless"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>n</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>lambda_n</ci>
                  <apply>
                     <minus/>
                     <ci>n_infinity</ci>
                     <ci>n</ci>
                  </apply>
               </apply>
               <ci>tau_n</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>n_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>Vn</ci>
                           <ci>V</ci>
                        </apply>
                        <ci>Sn</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_n</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">9.09</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>V</ci>
                           <ci>Vn</ci>
                        </apply>
                        <ci>Sn</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="K_ATP_current">
      <variable cmeta:id="K_ATP_current_i_K_ATP" name="i_K_ATP" public_interface="out" units="picoA"/>
      <variable initial_value="150" name="g_K_ATP" units="picoS"/>
      <variable name="V_K" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K_ATP</ci>
            <apply>
               <times/>
               <ci>g_K_ATP</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_Ca_current">
      <variable cmeta:id="fast_Ca_current_i_Ca_f" name="i_Ca_f" public_interface="out" units="picoA"/>
      <variable initial_value="100" name="V_Ca" public_interface="out" units="millivolt"/>
      <variable initial_value="810" name="g_Ca_f" units="picoS"/>
      <variable name="time" private_interface="out" public_interface="in" units="millisecond"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="m_f_infinity" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca_f</ci>
            <apply>
               <times/>
               <ci>g_Ca_f</ci>
               <ci>m_f_infinity</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="fast_Ca_channel_m_gate">
      <variable name="m_f_infinity" public_interface="out" units="dimensionless"/>
      <variable initial_value="-20" name="Vm_f" units="millivolt"/>
      <variable initial_value="7.5" name="Sm_f" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>m_f_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>Vm_f</ci>
                           <ci>V</ci>
                        </apply>
                        <ci>Sm_f</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_Ca_current">
      <variable cmeta:id="slow_Ca_current_i_Ca_s" name="i_Ca_s" public_interface="out" units="picoA"/>
      <variable initial_value="510" name="g_Ca_s" units="picoS"/>
      <variable name="time" private_interface="out" public_interface="in" units="millisecond"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="V_Ca" public_interface="in" units="millivolt"/>
      <variable name="m_s_infinity" private_interface="in" units="dimensionless"/>
      <variable name="jm" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca_s</ci>
            <apply>
               <times/>
               <ci>g_Ca_s</ci>
               <ci>m_s_infinity</ci>
               <apply>
                  <minus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <ci>jm</ci>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_Ca</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_Ca_channel_m_gate">
      <variable name="m_s_infinity" public_interface="out" units="dimensionless"/>
      <variable initial_value="-16" name="Vm_s" units="millivolt"/>
      <variable initial_value="10" name="Sm_s" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>m_s_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>Vm_s</ci>
                           <ci>V</ci>
                        </apply>
                        <ci>Sm_s</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="slow_Ca_channel_j_gate">
      <variable initial_value="0.12" name="jm" public_interface="out" units="dimensionless"/>
      <variable name="jm_infinity" units="dimensionless"/>
      <variable initial_value="-53" name="Vj" units="millivolt"/>
      <variable name="tau_j" units="millisecond"/>
      <variable initial_value="2" name="Sj" units="millivolt"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>jm_infinity</ci>
            <apply>
               <minus/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <plus/>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <ci>Vj</ci>
                           </apply>
                           <ci>Sj</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>tau_j</ci>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">50000</cn>
                  <apply>
                     <plus/>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>V</ci>
                              <ci>Vj</ci>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
                        </apply>
                     </apply>
                     <apply>
                        <exp/>
                        <apply>
                           <divide/>
                           <apply>
                              <minus/>
                              <ci>Vj</ci>
                              <ci>V</ci>
                           </apply>
                           <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millivolt">4</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="millisecond">1500</cn>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>jm</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>jm_infinity</ci>
                  <ci>jm</ci>
               </apply>
               <ci>tau_j</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_current_total">
      <variable name="i_Ca" public_interface="out" units="picoA"/>
      <variable name="i_Ca_f" public_interface="in" units="picoA"/>
      <variable name="i_Ca_s" public_interface="in" units="picoA"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca</ci>
            <apply>
               <plus/>
               <ci>i_Ca_f</ci>
               <ci>i_Ca_s</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="K_Ca_current">
      <variable cmeta:id="K_Ca_current_i_K_Ca" name="i_K_Ca" public_interface="out" units="picoA"/>
      <variable initial_value="1200" name="g_K_Ca" units="picoS"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="V_K" public_interface="in" units="millivolt"/>
      <variable name="Ca_i" public_interface="in" units="micromolar"/>
      <variable initial_value="0.55" name="kdkca" units="micromolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K_Ca</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>g_K_Ca</ci>
                     <apply>
                        <power/>
                        <ci>Ca_i</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">5</cn>
                     </apply>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <power/>
                        <ci>Ca_i</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">5</cn>
                     </apply>
                     <apply>
                        <power/>
                        <ci>kdkca</ci>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">5</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="CRAC_current">
      <variable cmeta:id="CRAC_current_i_CRAC" name="i_CRAC" public_interface="out" units="picoA"/>
      <variable initial_value="75" name="g_CRAC" units="picoS"/>
      <variable initial_value="0" name="V_CRAC" public_interface="out" units="millivolt"/>
      <variable name="time" private_interface="out" public_interface="in" units="millisecond"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <variable name="Ca_er" private_interface="out" public_interface="in" units="micromolar"/>
      <variable name="r_infinity" private_interface="in" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_CRAC</ci>
            <apply>
               <times/>
               <ci>g_CRAC</ci>
               <ci>r_infinity</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_CRAC</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="CRAC_r_gate">
      <variable name="r_infinity" public_interface="out" units="dimensionless"/>
      <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="Ca_er" public_interface="in" units="micromolar"/>
      <variable initial_value="40" name="Ca_er_bar" units="micromolar"/>
      <variable initial_value="3" name="sloper" units="micromolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>r_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <minus/>
                           <ci>Ca_er</ci>
                           <ci>Ca_er_bar</ci>
                        </apply>
                        <ci>sloper</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="leak_current">
      <variable name="i_leak" public_interface="out" units="picoA"/>
      <variable initial_value="0" name="g_leak" units="picoS"/>
      <variable name="V_CRAC" public_interface="in" units="millivolt"/>
      <variable name="time" private_interface="out" public_interface="in" units="millisecond"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_leak</ci>
            <apply>
               <times/>
               <ci>g_leak</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>V_CRAC</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="ER_parameters">
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable name="Ca_i" public_interface="in" units="micromolar"/>
      <variable cmeta:id="ER_parameters_J_er_p" name="J_er_p" public_interface="out" units="micromolar_per_millisecond"/>
      <variable initial_value="0" name="IP3" units="micromolar"/>
      <variable initial_value="0.09" name="kerp" units="micromolar"/>
      <variable initial_value="0.24" name="verp" units="micromolar_per_millisecond"/>
      <variable initial_value="0.35" name="dact" units="micromolar"/>
      <variable initial_value="0.4" name="dinh" units="micromolar"/>
      <variable initial_value="0.2" name="dip3" units="micromolar"/>
      <variable name="a_infinity" units="dimensionless"/>
      <variable name="b_infinity" units="dimensionless"/>
      <variable name="h_infinity" units="dimensionless"/>
      <variable name="O" public_interface="out" units="per_millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_er_p</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>verp</ci>
                  <apply>
                     <power/>
                     <ci>Ca_i</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Ca_i</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>kerp</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>a_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>dact</ci>
                     <ci>Ca_i</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>b_infinity</ci>
            <apply>
               <divide/>
               <ci>IP3</ci>
               <apply>
                  <plus/>
                  <ci>IP3</ci>
                  <ci>dip3</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>h_infinity</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>Ca_i</ci>
                     <ci>dinh</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>O</ci>
            <apply>
               <times/>
               <apply>
                  <power/>
                  <ci>a_infinity</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
               </apply>
               <apply>
                  <power/>
                  <ci>b_infinity</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
               </apply>
               <apply>
                  <power/>
                  <ci>h_infinity</ci>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">3</cn>
               </apply>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_millisecond">1</cn>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_equations">
      <variable name="time" private_interface="out" public_interface="in" units="millisecond"/>
      <variable initial_value="60" name="Ca_er" public_interface="out" units="micromolar"/>
      <variable initial_value="0.11" name="Ca_i" private_interface="out" public_interface="out" units="micromolar"/>
      <variable name="J_er_tot" public_interface="out" units="micromolar_per_millisecond"/>
      <variable name="J_er_p" private_interface="in" units="micromolar_per_millisecond"/>
      <variable name="J_er_IP3" units="micromolar_per_millisecond"/>
      <variable name="J_er_leak" units="micromolar_per_millisecond"/>
      <variable name="J_mem_tot" public_interface="in" units="micromolar_per_millisecond"/>
      <variable name="O" private_interface="in" units="per_millisecond"/>
      <variable initial_value="0.003" name="perl" units="per_millisecond"/>
      <variable initial_value="250" name="lambda_er" units="dimensionless"/>
      <variable initial_value="1" name="sigma_er" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_er_tot</ci>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <ci>J_er_leak</ci>
                  <ci>J_er_IP3</ci>
               </apply>
               <ci>J_er_p</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_er_leak</ci>
            <apply>
               <times/>
               <ci>perl</ci>
               <apply>
                  <minus/>
                  <ci>Ca_er</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_er_IP3</ci>
            <apply>
               <times/>
               <ci>O</ci>
               <apply>
                  <minus/>
                  <ci>Ca_er</ci>
                  <ci>Ca_i</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_er</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>J_er_tot</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>lambda_er</ci>
                  <ci>sigma_er</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_i</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <divide/>
                  <ci>J_er_tot</ci>
                  <ci>lambda_er</ci>
               </apply>
               <ci>J_mem_tot</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Ca_membrane_flux">
      <variable name="time" public_interface="in" units="millisecond"/>
      <variable initial_value="0.35" name="kmp" units="micromolar"/>
      <variable initial_value="0.08" name="vmp" units="micromolar"/>
      <variable initial_value="0.000003607" name="gamma" units="micromolar_per_picoA"/>
      <variable name="i_Ca" public_interface="in" units="picoA"/>
      <variable name="Ca_i" public_interface="in" units="micromolar"/>
      <variable name="J_mem_tot" public_interface="out" units="micromolar_per_millisecond"/>
      <variable name="Jmp" units="micromolar"/>
      <variable initial_value="0.01" name="f" units="per_millisecond"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Jmp</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>vmp</ci>
                  <apply>
                     <power/>
                     <ci>Ca_i</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Ca_i</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>kmp</ci>
                     <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>J_mem_tot</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>f</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>gamma</ci>
                     <ci>i_Ca</ci>
                  </apply>
                  <ci>Jmp</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="K_current">
            <component_ref component="K_channel_n_gate"/>
         </component_ref>
         <component_ref component="K_ATP_current"/>
         <component_ref component="fast_Ca_current">
            <component_ref component="fast_Ca_channel_m_gate"/>
         </component_ref>
         <component_ref component="slow_Ca_current">
            <component_ref component="slow_Ca_channel_m_gate"/>
            <component_ref component="slow_Ca_channel_j_gate"/>
         </component_ref>
         <component_ref component="K_Ca_current"/>
         <component_ref component="CRAC_current">
            <component_ref component="CRAC_r_gate"/>
         </component_ref>
         <component_ref component="ER_parameters"/>
         <component_ref component="Ca_equations"/>
         <component_ref component="Ca_membrane_flux"/>
      </component_ref>
   </group>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="K_current">
         <component_ref component="K_channel_n_gate"/>
      </component_ref>
      <component_ref component="fast_Ca_current">
         <component_ref component="fast_Ca_channel_m_gate"/>
      </component_ref>
      <component_ref component="slow_Ca_current">
         <component_ref component="slow_Ca_channel_m_gate"/>
         <component_ref component="slow_Ca_channel_j_gate"/>
      </component_ref>
      <component_ref component="CRAC_current">
         <component_ref component="CRAC_r_gate"/>
      </component_ref>
      <component_ref component="Ca_equations">
         <component_ref component="ER_parameters"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="membrane" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="K_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="K_ATP_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="fast_Ca_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="slow_Ca_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="K_Ca_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="CRAC_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Ca_equations" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Ca_membrane_flux" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="leak_current" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="K_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_K" variable_2="i_K"/>
   </connection>
   <connection>
      <map_components component_1="K_ATP_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_K_ATP" variable_2="i_K_ATP"/>
   </connection>
   <connection>
      <map_components component_1="fast_Ca_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="slow_Ca_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="K_Ca_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_K_Ca" variable_2="i_K_Ca"/>
   </connection>
   <connection>
      <map_components component_1="CRAC_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_CRAC" variable_2="i_CRAC"/>
   </connection>
   <connection>
      <map_components component_1="leak_current" component_2="membrane"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_leak" variable_2="i_leak"/>
   </connection>
   <connection>
      <map_components component_1="K_current" component_2="K_ATP_current"/>
      <map_variables variable_1="V_K" variable_2="V_K"/>
   </connection>
   <connection>
      <map_components component_1="K_current" component_2="K_Ca_current"/>
      <map_variables variable_1="V_K" variable_2="V_K"/>
   </connection>
   <connection>
      <map_components component_1="fast_Ca_current" component_2="slow_Ca_current"/>
      <map_variables variable_1="V_Ca" variable_2="V_Ca"/>
   </connection>
   <connection>
      <map_components component_1="slow_Ca_current" component_2="Ca_current_total"/>
      <map_variables variable_1="i_Ca_s" variable_2="i_Ca_s"/>
   </connection>
   <connection>
      <map_components component_1="fast_Ca_current" component_2="Ca_current_total"/>
      <map_variables variable_1="i_Ca_f" variable_2="i_Ca_f"/>
   </connection>
   <connection>
      <map_components component_1="Ca_current_total" component_2="membrane"/>
      <map_variables variable_1="i_Ca" variable_2="i_Ca"/>
   </connection>
   <connection>
      <map_components component_1="Ca_current_total" component_2="Ca_membrane_flux"/>
      <map_variables variable_1="i_Ca" variable_2="i_Ca"/>
   </connection>
   <connection>
      <map_components component_1="CRAC_current" component_2="Ca_equations"/>
      <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
   </connection>
   <connection>
      <map_components component_1="K_Ca_current" component_2="Ca_equations"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
   </connection>
   <connection>
      <map_components component_1="CRAC_current" component_2="leak_current"/>
      <map_variables variable_1="V_CRAC" variable_2="V_CRAC"/>
   </connection>
   <connection>
      <map_components component_1="Ca_equations" component_2="Ca_membrane_flux"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="J_mem_tot" variable_2="J_mem_tot"/>
   </connection>
   <connection>
      <map_components component_1="ER_parameters" component_2="Ca_equations"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
      <map_variables variable_1="J_er_p" variable_2="J_er_p"/>
      <map_variables variable_1="O" variable_2="O"/>
   </connection>
   <connection>
      <map_components component_1="K_current" component_2="K_channel_n_gate"/>
      <map_variables variable_1="n" variable_2="n"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="fast_Ca_current" component_2="fast_Ca_channel_m_gate"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="m_f_infinity" variable_2="m_f_infinity"/>
   </connection>
   <connection>
      <map_components component_1="slow_Ca_current" component_2="slow_Ca_channel_m_gate"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="m_s_infinity" variable_2="m_s_infinity"/>
   </connection>
   <connection>
      <map_components component_1="slow_Ca_current" component_2="slow_Ca_channel_j_gate"/>
      <map_variables variable_1="jm" variable_2="jm"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>
   <connection>
      <map_components component_1="CRAC_current" component_2="CRAC_r_gate"/>
      <map_variables variable_1="r_infinity" variable_2="r_infinity"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
   </connection>
</model>