- Author:
- neilstephen2001 <npar276@aucklanduni.ac.nz>
- Date:
- 2021-01-25 20:32:57+13:00
- Desc:
- modifying composite annotations
- Permanent Source URI:
- http://models.cellml.org/workspace/7f1/rawfile/8f4c90b98ea0253dc01be6f1950c133cec18ca5e/devries_sherman_2000/model/devries_sherman_2000.cellml
<?xml version="1.0" encoding="UTF-8"?>
<model xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:semsim="http://www.bhi.washington.edu/SemSim#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="devries_sherman_2000" cmeta:id="devries_2000">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Enhancement Of Emergent Bursting In Pancreatic Beta-Cells</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>This model has been validated and is known to run in OpenCell and COR to recreate the published results (figure 1). The paper by Devries and Sherman describes several different parameter sets including single cell and two cell sets, and also describes a deterministic and a stochastic model. This CellML model uses equations 1-3 and the parameters given in table 1 for the single cell deterministic model. For this particular model lambda = 0.9. Note there is a unit inconsistency in the equation for V, but our attempts to balance the units for this equation resulted in a model which no longer worked.</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>ABSTRACT: Secretion of insulin by electrically coupled populations of pancreatic beta -cells is governed by bursting electrical activity. Isolated beta -cells, however, exhibit atypical bursting or continuous spike activity. We study bursting as an emergent property of the population, focussing on interactions among the subclass of spiking cells. These are modelled by equipping the fast subsystem with a saddle-node-loop bifurcation, which makes it monostable. Such cells can only spike tonically or remain silent when isolated, but can be induced to burst with weak diffusive coupling. With stronger coupling, the cells revert to tonic spiking. We demonstrate that the addition of noise dramatically increases, via a phenomenon like stochastic resonance, the coupling range over which bursting is seen.</para>
<para>The original paper reference is cited below:</para>
<para>
Channel Sharing in Pancreatic Beta-Cells Revisited: Enhancement of Emergent Bursting by Noise, Gerda De Vries and Arthur Sherman, 2000,
<emphasis>The Journal Of Theoretical Biology</emphasis>
, 207, 513-530.
<ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=11093836&dopt=Abstract">PubMed ID: 11093836</ulink>
</para>
<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>diagram of the model</title>
</objectinfo>
<imagedata fileref="devries_2000.png" />
</imageobject>
</mediaobject>
<caption>A schematic representation of the four transmembrane currents captured by the De Vries and Sherman 2000 pancreatic beta-cell model.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<units name="nanoS">
<unit prefix="nano" units="siemens" />
</units>
<units name="millivolt">
<unit prefix="milli" units="volt" />
</units>
<units name="millisecond">
<unit prefix="milli" units="second" />
</units>
<units name="picoA">
<unit prefix="pico" units="ampere" />
</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="-64.0" name="V" public_interface="out" units="millivolt" />
<variable cmeta:id="membrane.tau" initial_value="20.0" name="tau" public_interface="out" units="millisecond" />
<variable cmeta:id="membrane.time" name="time" public_interface="in" units="millisecond" />
<variable cmeta:id="membrane.i_K" name="i_K" public_interface="in" units="picoA" />
<variable cmeta:id="membrane.i_K_ATP" name="i_K_ATP" public_interface="in" units="picoA" />
<variable cmeta:id="membrane.i_Ca" name="i_Ca" public_interface="in" units="picoA" />
<variable cmeta:id="membrane.i_s" name="i_s" public_interface="in" units="picoA" />
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="membrane_voltage_diff_eq">
<eq />
<apply>
<diff />
<bvar>
<ci>time</ci>
</bvar>
<ci>V</ci>
</apply>
<apply>
<divide />
<apply>
<minus />
<apply>
<plus />
<ci>i_Ca</ci>
<ci>i_K</ci>
<ci>i_K_ATP</ci>
<ci>i_s</ci>
</apply>
</apply>
<ci>tau</ci>
</apply>
</apply>
</math>
</component>
<component name="calcium_current">
<variable cmeta:id="calcium_current_i_Ca" name="i_Ca" public_interface="out" units="picoA" />
<variable cmeta:id="calcium_current.g_Ca" initial_value="3.6" name="g_Ca" units="nanoS" />
<variable cmeta:id="calcium_current.V_Ca" initial_value="25.0" name="V_Ca" units="millivolt" />
<variable cmeta:id="calcium_current.time" name="time" public_interface="in" units="millisecond" />
<variable cmeta:id="calcium_current.V" name="V" public_interface="in" private_interface="out" units="millivolt" />
<variable cmeta:id="calcium_current.m_infinity" name="m_infinity" private_interface="in" units="dimensionless" />
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_Ca_calculation">
<eq />
<ci>i_Ca</ci>
<apply>
<times />
<ci>g_Ca</ci>
<ci>m_infinity</ci>
<apply>
<minus />
<ci>V</ci>
<ci>V_Ca</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_current_m_gate">
<variable cmeta:id="calcium_current_m_gate.m_infinity" name="m_infinity" public_interface="out" units="dimensionless" />
<variable cmeta:id="calcium_current_m_gate.V_m" initial_value="-20.0" name="V_m" units="millivolt" />
<variable cmeta:id="calcium_current_m_gate.theta_m" initial_value="12.0" name="theta_m" units="millivolt" />
<variable cmeta:id="calcium_current_m_gate.V" name="V" public_interface="in" units="millivolt" />
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="m_infinity_calculation">
<eq />
<ci>m_infinity</ci>
<apply>
<divide />
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.0</cn>
<apply>
<plus />
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.0</cn>
<apply>
<exp />
<apply>
<divide />
<apply>
<minus />
<ci>V_m</ci>
<ci>V</ci>
</apply>
<ci>theta_m</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="potassium_current">
<variable cmeta:id="potassium_current_i_K" name="i_K" public_interface="out" units="picoA" />
<variable cmeta:id="potassium_current.V_K" initial_value="-75.0" name="V_K" public_interface="out" units="millivolt" />
<variable cmeta:id="potassium_current.g_K" initial_value="10.0" name="g_K" units="nanoS" />
<variable cmeta:id="potassium_current.time" name="time" public_interface="in" private_interface="out" units="millisecond" />
<variable cmeta:id="potassium_current.V" name="V" public_interface="in" private_interface="out" units="millivolt" />
<variable cmeta:id="potassium_current.tau" name="tau" public_interface="in" private_interface="out" units="millisecond" />
<variable cmeta:id="potassium_current.n" name="n" private_interface="in" units="dimensionless" />
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_K_calculation">
<eq />
<ci>i_K</ci>
<apply>
<times />
<ci>g_K</ci>
<ci>n</ci>
<apply>
<minus />
<ci>V</ci>
<ci>V_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="potassium_current_n_gate">
<variable cmeta:id="potassium_current_n_gate.n" initial_value="0.01" name="n" public_interface="out" units="dimensionless" />
<variable cmeta:id="potassium_current_n_gate.n_infinity" name="n_infinity" units="dimensionless" />
<variable cmeta:id="potassium_current_n_gate.V_n" initial_value="-17.0" name="V_n" units="millivolt" />
<variable cmeta:id="potassium_current_n_gate.theta_n" initial_value="5.6" name="theta_n" units="millivolt" />
<variable cmeta:id="potassium_current_n_gate.lambda" initial_value="0.9" name="lambda" units="dimensionless" />
<variable cmeta:id="potassium_current_n_gate.V" name="V" public_interface="in" units="millivolt" />
<variable cmeta:id="potassium_current_n_gate.tau" name="tau" public_interface="in" units="millisecond" />
<variable cmeta:id="potassium_current_n_gate.time" name="time" public_interface="in" units="millisecond" />
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="n_diff_eq">
<eq />
<apply>
<diff />
<bvar>
<ci>time</ci>
</bvar>
<ci>n</ci>
</apply>
<apply>
<divide />
<apply>
<times />
<ci>lambda</ci>
<apply>
<minus />
<ci>n_infinity</ci>
<ci>n</ci>
</apply>
</apply>
<ci>tau</ci>
</apply>
</apply>
<apply id="n_infinity_calculation">
<eq />
<ci>n_infinity</ci>
<apply>
<divide />
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.0</cn>
<apply>
<plus />
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.0</cn>
<apply>
<exp />
<apply>
<divide />
<apply>
<minus />
<ci>V_n</ci>
<ci>V</ci>
</apply>
<ci>theta_n</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_current">
<variable cmeta:id="slow_current_i_s" name="i_s" public_interface="out" units="picoA" />
<variable cmeta:id="slow_current.g_s" initial_value="4.0" name="g_s" units="nanoS" />
<variable cmeta:id="slow_current.time" name="time" public_interface="in" private_interface="out" units="millisecond" />
<variable cmeta:id="slow_current.V" name="V" public_interface="in" private_interface="out" units="millivolt" />
<variable cmeta:id="slow_current.V_K" name="V_K" public_interface="in" units="millivolt" />
<variable cmeta:id="slow_current.s" name="s" private_interface="in" units="dimensionless" />
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_s_calculation">
<eq />
<ci>i_s</ci>
<apply>
<times />
<ci>g_s</ci>
<ci>s</ci>
<apply>
<minus />
<ci>V</ci>
<ci>V_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="slow_current_s_gate">
<variable cmeta:id="slow_current_s_gate.s" initial_value="0.01" name="s" public_interface="out" units="dimensionless" />
<variable cmeta:id="slow_current_s_gate.s_infinity" name="s_infinity" units="dimensionless" />
<variable cmeta:id="slow_current_s_gate.V_s" initial_value="-22.0" name="V_s" units="millivolt" />
<variable cmeta:id="slow_current_s_gate.theta_s" initial_value="8.0" name="theta_s" units="millivolt" />
<variable cmeta:id="slow_current_s_gate.tau_s" initial_value="20000.0" name="tau_s" units="millisecond" />
<variable cmeta:id="slow_current_s_gate.V" name="V" public_interface="in" units="millivolt" />
<variable cmeta:id="slow_current_s_gate.time" name="time" public_interface="in" units="millisecond" />
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="s_diff_eq">
<eq />
<apply>
<diff />
<bvar>
<ci>time</ci>
</bvar>
<ci>s</ci>
</apply>
<apply>
<divide />
<apply>
<minus />
<ci>s_infinity</ci>
<ci>s</ci>
</apply>
<ci>tau_s</ci>
</apply>
</apply>
<apply id="s_infinity_calculation">
<eq />
<ci>s_infinity</ci>
<apply>
<divide />
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.0</cn>
<apply>
<plus />
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1.0</cn>
<apply>
<exp />
<apply>
<divide />
<apply>
<minus />
<ci>V_s</ci>
<ci>V</ci>
</apply>
<ci>theta_s</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="ATP_sensitive_potassium_current">
<variable cmeta:id="ATP_sensitive_potassium_current_i_K_ATP" name="i_K_ATP" public_interface="out" units="picoA" />
<variable cmeta:id="ATP_sensitive_potassium_current.g_K_ATP" initial_value="1.2" name="g_K_ATP" units="nanoS" />
<variable cmeta:id="ATP_sensitive_potassium_current.p" initial_value="0.5" name="p" units="dimensionless" />
<variable cmeta:id="ATP_sensitive_potassium_current.time" name="time" public_interface="in" units="millisecond" />
<variable cmeta:id="ATP_sensitive_potassium_current.V" name="V" public_interface="in" units="millivolt" />
<variable cmeta:id="ATP_sensitive_potassium_current.V_K" name="V_K" public_interface="in" units="millivolt" />
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply id="i_K_ATP_calculation">
<eq />
<ci>i_K_ATP</ci>
<apply>
<times />
<ci>g_K_ATP</ci>
<ci>p</ci>
<apply>
<minus />
<ci>V</ci>
<ci>V_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="encapsulation" />
<component_ref component="slow_current">
<component_ref component="slow_current_s_gate" />
</component_ref>
<component_ref component="potassium_current">
<component_ref component="potassium_current_n_gate" />
</component_ref>
<component_ref component="calcium_current">
<component_ref component="calcium_current_m_gate" />
</component_ref>
</group>
<connection>
<map_components component_1="environment" component_2="potassium_current" />
<map_variables variable_1="time" variable_2="time" />
</connection>
<connection>
<map_components component_1="environment" component_2="membrane" />
<map_variables variable_1="time" variable_2="time" />
</connection>
<connection>
<map_components component_1="membrane" component_2="potassium_current" />
<map_variables variable_1="V" variable_2="V" />
<map_variables variable_1="tau" variable_2="tau" />
</connection>
<connection>
<map_components component_1="potassium_current_n_gate" component_2="potassium_current" />
<map_variables variable_1="n" variable_2="n" />
</connection>
<connection>
<map_components component_1="potassium_current" component_2="ATP_sensitive_potassium_current" />
<map_variables variable_1="V_K" variable_2="V_K" />
</connection>
<connection>
<map_components component_1="membrane" component_2="slow_current" />
<map_variables variable_1="V" variable_2="V" />
</connection>
<connection>
<map_components component_1="potassium_current" component_2="potassium_current_n_gate" />
<map_variables variable_1="tau" variable_2="tau" />
<map_variables variable_1="V" variable_2="V" />
<map_variables variable_1="time" variable_2="time" />
</connection>
<connection>
<map_components component_1="slow_current" component_2="membrane" />
<map_variables variable_1="i_s" variable_2="i_s" />
</connection>
<connection>
<map_components component_1="environment" component_2="slow_current" />
<map_variables variable_1="time" variable_2="time" />
</connection>
<connection>
<map_components component_1="membrane" component_2="calcium_current" />
<map_variables variable_1="V" variable_2="V" />
</connection>
<connection>
<map_components component_1="membrane" component_2="ATP_sensitive_potassium_current" />
<map_variables variable_1="V" variable_2="V" />
</connection>
<connection>
<map_components component_1="potassium_current" component_2="membrane" />
<map_variables variable_1="i_K" variable_2="i_K" />
</connection>
<connection>
<map_components component_1="potassium_current" component_2="slow_current" />
<map_variables variable_1="V_K" variable_2="V_K" />
</connection>
<connection>
<map_components component_1="calcium_current" component_2="calcium_current_m_gate" />
<map_variables variable_1="V" variable_2="V" />
</connection>
<connection>
<map_components component_1="environment" component_2="ATP_sensitive_potassium_current" />
<map_variables variable_1="time" variable_2="time" />
</connection>
<connection>
<map_components component_1="slow_current_s_gate" component_2="slow_current" />
<map_variables variable_1="s" variable_2="s" />
</connection>
<connection>
<map_components component_1="ATP_sensitive_potassium_current" component_2="membrane" />
<map_variables variable_1="i_K_ATP" variable_2="i_K_ATP" />
</connection>
<connection>
<map_components component_1="calcium_current_m_gate" component_2="calcium_current" />
<map_variables variable_1="m_infinity" variable_2="m_infinity" />
</connection>
<connection>
<map_components component_1="environment" component_2="calcium_current" />
<map_variables variable_1="time" variable_2="time" />
</connection>
<connection>
<map_components component_1="slow_current" component_2="slow_current_s_gate" />
<map_variables variable_1="time" variable_2="time" />
<map_variables variable_1="V" variable_2="V" />
</connection>
<connection>
<map_components component_1="calcium_current" component_2="membrane" />
<map_variables variable_1="i_Ca" variable_2="i_Ca" />
</connection>
</model>
