- Author:
- Tessa Paris <tpar054@aucklanduni.ac.nz>
- Date:
- 2009-12-17 10:14:39+13:00
- Desc:
- Created the model using the authors original code and the CellML model by Bertram et. al. 1995.
The model runs in pcenv and COR though some dimensional inconsistancies are noted. The model can produce output that matches the figures in the paper.
- Permanent Source URI:
- http://models.cellml.org/workspace/mears_sheppard_atwater_rojas_bertram_sherman_1997/rawfile/73fcb91fcdd19b3e729736c64b930740fed614a3/mears_sheppard_atwater_rojas_bertram_sherman_1997.cellml
<?xml version="1.0" encoding="utf-8"?>
<!--
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 xml:base="file:///C:/Documents%20and%20Settings/tpar054/Desktop/Repository%20models/Mears_1997/Mears_1997.cellml" name="Mears_1997_Low_Glucose_Concentration" cmeta:id="Mears_1997_Low_Glucose_Concentration" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Evidence That Calcium Release-activated Current Mediates the Biphasic Electrical Activity of Mouse Pancreatic β-Cells</title>
<author>
<firstname>Tessa</firstname>
<surname>Paris</surname>
<affiliation>
<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>This model works in both pcenv and COR and it can produce figures similar to in the paper, namely figures 2 and 6. </para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>The electrical response of pancreatic β-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 (I CRAN) that activates as intracellular calcium stores empty during exposure to basal glucose (Bertram et al., 1995). The stores refill during Phase 1, deactivating I CRAN 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 β-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 β-cel.</para>
<para>
The complete original paper reference is cited below:
</para>
<para>Evidence That Calcium Release-activated Current Mediates the Biphasic Electrical Activity of Mouse Pancreatic β-Cells, Authors, year, <emphasis>Journal name</emphasis>, edition, firstpage-lastpage. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/9002424">PubMed ID: 9002424</ulink></para>
<informalfigure float="0" id="______">
<mediaobject>
<imageobject>
<objectinfo>
<title>Image title goes here</title>
</objectinfo>
<imagedata fileref="author_year.png"/>
</imageobject>
</mediaobject>
<caption>Image caption goes here</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-snytax-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>Your email address (keep in mind this will be freely available)</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="author_model_year">
<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 β-Cells</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">
<rdf:value>
— — — Acute insulin response — Mathematical model — Glucose
<rdf:li>Islet of Langerhans</rdf:li>
<rdf:li>Stimulus-secretion coupling</rdf:li>
<rdf:li>Endoplasmic reticulum</rdf:li>
<rdf:li>Glucose</rdf:li>
<rdf:li>Acute insulin response</rdf:li>
</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">
Richard Bertram, Paul Smolen, Arthur Sherman, David Mears, Illani Atwater, Franz Martin and Bernat Soria, 1995,
<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>David</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>Norman</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>Illani</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>Betrtram</vCard:Family>
<vCard:Given>Richard</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>Arthur</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 β-Cells</dc:title>
<dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:ParseType="Resource">
<dcterms:W3CDTF>1996-08-15</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 xmlns="http://www.cellml.org/cellml/1.0#" name="millisecond">
<unit units="second" prefix="milli"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="per_millisecond">
<unit units="millisecond" exponent="-1"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="micromolar_per_millisecond">
<unit units="micromolar"/>
<unit units="millisecond" exponent="-1"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="micromolar_per_millisecond_picoA">
<unit units="micromolar"/>
<unit units="millisecond" exponent="-1"/>
<unit units="picoA" exponent="-1"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="millivolt">
<unit units="volt" prefix="milli"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="micromolar">
<unit units="mole" prefix="micro"/>
<unit units="litre" exponent="-1"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="micrometer_3">
<unit units="meter" prefix="micro" exponent="3"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="picoS">
<unit units="siemens" prefix="pico"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="picoA">
<unit units="ampere" prefix="pico"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="micromolar_per_picoA">
<unit units="micromolar"/>
<unit units="ampere" prefix="pico" exponent="-1"/>
</units>
<units xmlns="http://www.cellml.org/cellml/1.0#" name="femtoF">
<unit units="farad" prefix="femto"/>
</units>
<component xmlns="http://www.cellml.org/cellml/1.0#" name="environment">
<variable name="time" units="millisecond" public_interface="out"/>
</component>
<component xmlns="http://www.cellml.org/cellml/1.0#" name="membrane">
<variable name="V" units="millivolt" initial_value="-61" public_interface="out"/>
<variable name="Cm" units="femtoF" initial_value="6158"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="i_K" units="picoA" public_interface="in"/>
<variable name="i_K_Ca" units="picoA" public_interface="in"/>
<variable name="i_K_ATP" units="picoA" public_interface="in"/>
<variable name="i_CRAC" units="picoA" public_interface="in"/>
<variable name="i_Ca" units="picoA" public_interface="in"/>
<variable name="i_leak" units="picoA" public_interface="in"/>
<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 xmlns="http://www.cellml.org/cellml/1.0#" name="K_current">
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="i_K" units="picoA" public_interface="out"/>
<variable name="V_K" units="millivolt" initial_value="-70" public_interface="out"/>
<variable name="g_K" units="picoS" initial_value="3900"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="n" units="dimensionless" private_interface="in"/>
<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 xmlns="http://www.cellml.org/cellml/1.0#" name="K_channel_n_gate">
<variable name="n" units="dimensionless" initial_value="0.0005" public_interface="out"/>
<variable name="n_infinity" units="dimensionless"/>
<variable name="tau_n" units="millisecond"/>
<variable name="Vn" units="millivolt" initial_value="-15"/>
<variable name="Sn" units="millivolt" initial_value="6"/>
<variable name="lambda_n" units="dimensionless" initial_value="1.85"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<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 cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn 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 cellml:units="millisecond">9.09</cn>
<apply>
<plus/>
<cn 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 xmlns="http://www.cellml.org/cellml/1.0#" name="K_ATP_current">
<variable name="i_K_ATP" units="picoA" public_interface="out"/>
<variable name="g_K_ATP" units="picoS" initial_value="150"/>
<variable name="V_K" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<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 xmlns="http://www.cellml.org/cellml/1.0#" name="fast_Ca_current">
<variable name="i_Ca_f" units="picoA" public_interface="out"/>
<variable name="V_Ca" units="millivolt" initial_value="100" public_interface="out"/>
<variable name="g_Ca_f" units="picoS" initial_value="810"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="m_f_infinity" units="dimensionless" private_interface="in"/>
<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 xmlns="http://www.cellml.org/cellml/1.0#" name="fast_Ca_channel_m_gate">
<variable name="m_f_infinity" units="dimensionless" public_interface="out"/>
<variable name="Vm_f" units="millivolt" initial_value="-20"/>
<variable name="Sm_f" units="millivolt" initial_value="7.5"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>m_f_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn 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 xmlns="http://www.cellml.org/cellml/1.0#" name="slow_Ca_current">
<variable name="i_Ca_s" units="picoA" public_interface="out"/>
<variable name="g_Ca_s" units="picoS" initial_value="510"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="V_Ca" units="millivolt" public_interface="in"/>
<variable name="m_s_infinity" units="dimensionless" private_interface="in"/>
<variable name="jm" units="dimensionless" private_interface="in"/>
<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 cellml:units="dimensionless">1</cn>
<ci>jm</ci>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>V_Ca</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component xmlns="http://www.cellml.org/cellml/1.0#" name="slow_Ca_channel_m_gate">
<variable name="m_s_infinity" units="dimensionless" public_interface="out"/>
<variable name="Vm_s" units="millivolt" initial_value="-16"/>
<variable name="Sm_s" units="millivolt" initial_value="10"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>m_s_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn 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 xmlns="http://www.cellml.org/cellml/1.0#" name="slow_Ca_channel_j_gate">
<variable name="jm" units="dimensionless" initial_value="0.12" public_interface="out"/>
<variable name="jm_infinity" units="dimensionless"/>
<variable name="Vj" units="millivolt" initial_value="-53"/>
<variable name="tau_j" units="millisecond"/>
<variable name="Sj" units="millivolt" initial_value="2"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>jm_infinity</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn 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 cellml:units="millisecond">50000</cn>
<apply>
<plus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
<ci>Vj</ci>
</apply>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>Vj</ci>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">4</cn>
</apply>
</apply>
</apply>
</apply>
<cn 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 xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_current_total">
<variable name="i_Ca" units="picoA" public_interface="out"/>
<variable name="i_Ca_f" units="picoA" public_interface="in"/>
<variable name="i_Ca_s" units="picoA" public_interface="in"/>
<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 xmlns="http://www.cellml.org/cellml/1.0#" name="K_Ca_current">
<variable name="i_K_Ca" units="picoA" public_interface="out"/>
<variable name="g_K_Ca" units="picoS" initial_value="1200"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="V_K" units="millivolt" public_interface="in"/>
<variable name="Ca_i" units="micromolar" public_interface="in"/>
<variable name="kdkca" units="micromolar" initial_value="0.55"/>
<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 cellml:units="dimensionless">5</cn>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">5</cn>
</apply>
<apply>
<power/>
<ci>kdkca</ci>
<cn cellml:units="dimensionless">5</cn>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>V_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component xmlns="http://www.cellml.org/cellml/1.0#" name="CRAC_current">
<variable name="i_CRAC" units="picoA" public_interface="out"/>
<variable name="g_CRAC" units="picoS" initial_value="75"/>
<variable name="V_CRAC" units="millivolt" initial_value="0" public_interface="out"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<variable name="Ca_er" units="micromolar" public_interface="in" private_interface="out"/>
<variable name="r_infinity" units="dimensionless" private_interface="in"/>
<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 xmlns="http://www.cellml.org/cellml/1.0#" name="CRAC_r_gate">
<variable name="r_infinity" units="dimensionless" public_interface="out"/>
<variable name="V" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Ca_er" units="micromolar" public_interface="in"/>
<variable name="Ca_er_bar" units="micromolar" initial_value="40"/>
<variable name="sloper" units="micromolar" initial_value="3"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>r_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn 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 xmlns="http://www.cellml.org/cellml/1.0#" name="leak_current">
<variable name="i_leak" units="picoA" public_interface="out"/>
<variable name="g_leak" units="picoS" initial_value="0"/>
<variable name="V_CRAC" units="millivolt" public_interface="in"/>
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="V" units="millivolt" public_interface="in" private_interface="out"/>
<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 xmlns="http://www.cellml.org/cellml/1.0#" name="ER_parameters">
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="Ca_i" units="micromolar" public_interface="in"/>
<variable name="J_er_p" units="micromolar_per_millisecond" public_interface="out"/>
<variable name="IP3" units="micromolar" initial_value="0"/>
<variable name="kerp" units="micromolar" initial_value="0.09"/>
<variable name="verp" units="micromolar_per_millisecond" initial_value="0.24"/>
<variable name="dact" units="micromolar" initial_value="0.35"/>
<variable name="dinh" units="micromolar" initial_value="0.4"/>
<variable name="dip3" units="micromolar" initial_value="0.2"/>
<variable name="a_infinity" units="dimensionless"/>
<variable name="b_infinity" units="dimensionless"/>
<variable name="h_infinity" units="dimensionless"/>
<variable name="O" units="per_millisecond" public_interface="out"/>
<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 cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>kerp</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>a_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn 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 cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn 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 cellml:units="dimensionless">3</cn>
</apply>
<apply>
<power/>
<ci>b_infinity</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<apply>
<power/>
<ci>h_infinity</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<cn cellml:units="per_millisecond">1</cn>
</apply>
</apply>
</math>
</component>
<component xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_equations">
<variable name="time" units="millisecond" public_interface="in" private_interface="out"/>
<variable name="Ca_er" units="micromolar" initial_value="60" public_interface="out"/>
<variable name="Ca_i" units="micromolar" initial_value="0.11" public_interface="out" private_interface="out"/>
<variable name="J_er_tot" units="micromolar_per_millisecond" public_interface="out"/>
<variable name="J_er_p" units="micromolar_per_millisecond" private_interface="in"/>
<variable name="J_er_IP3" units="micromolar_per_millisecond"/>
<variable name="J_er_leak" units="micromolar_per_millisecond"/>
<variable name="J_mem_tot" units="micromolar_per_millisecond" public_interface="in"/>
<variable name="O" units="per_millisecond" private_interface="in"/>
<variable name="perl" units="per_millisecond" initial_value="0.003"/>
<variable name="lambda_er" units="dimensionless" initial_value="250"/>
<variable name="sigma_er" units="dimensionless" initial_value="1"/>
<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 xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_membrane_flux">
<variable name="time" units="millisecond" public_interface="in"/>
<variable name="kmp" units="micromolar" initial_value="0.35"/>
<variable name="vmp" units="micromolar" initial_value="0.08"/>
<variable name="gamma" units="micromolar_per_picoA" initial_value="0.000003607"/>
<variable name="i_Ca" units="picoA" public_interface="in"/>
<variable name="Ca_i" units="micromolar" public_interface="in"/>
<variable name="J_mem_tot" units="micromolar_per_millisecond" public_interface="out"/>
<variable name="Jmp" units="micromolar"/>
<variable name="f" units="per_millisecond" initial_value="0.01"/>
<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 cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>Ca_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>kmp</ci>
<cn 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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="membrane" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="K_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="K_ATP_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="fast_Ca_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="slow_Ca_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="K_Ca_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="CRAC_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="Ca_equations" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="Ca_membrane_flux" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="leak_current" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="fast_Ca_current" component_2="membrane"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="slow_Ca_current" component_2="membrane"/>
<map_variables variable_1="V" variable_2="V"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="K_current" component_2="K_ATP_current"/>
<map_variables variable_1="V_K" variable_2="V_K"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="K_current" component_2="K_Ca_current"/>
<map_variables variable_1="V_K" variable_2="V_K"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="fast_Ca_current" component_2="slow_Ca_current"/>
<map_variables variable_1="V_Ca" variable_2="V_Ca"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="Ca_current_total" component_2="membrane"/>
<map_variables variable_1="i_Ca" variable_2="i_Ca"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="Ca_current_total" component_2="Ca_membrane_flux"/>
<map_variables variable_1="i_Ca" variable_2="i_Ca"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="CRAC_current" component_2="Ca_equations"/>
<map_variables variable_1="Ca_er" variable_2="Ca_er"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="K_Ca_current" component_2="Ca_equations"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<map_components component_1="CRAC_current" component_2="leak_current"/>
<map_variables variable_1="V_CRAC" variable_2="V_CRAC"/>
</connection>
<connection xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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 xmlns="http://www.cellml.org/cellml/1.0#">
<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>
