- Author:
- Hanne <Hanne@hanne-nielsens-macbook.local>
- Date:
- 2010-04-14 12:19:26+12:00
- Desc:
- Added images in ai and svg format, added clickable SVG to session file
- Permanent Source URI:
- https://models.cellml.org/workspace/mears_sheppard_atwater_rojas_bertram_sherman_1997/rawfile/27df2ef6b23052ad1c0659a2b70432829e7fd742/mears_sheppard_atwater_rojas_bertram_sherman_1997.cellml
<?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, University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>This CellML model runs in both PCEnv 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 </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>