Personal tools
You are here: Home Model Repository
 
Views

The Phantom Burster Model for Pancreatic Beta Cells

The CellML code.

XML source view

<model xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns="http://www.cellml.org/cellml/1.0#" cmeta:id="bertram_previte_sherman_kinard_satin_2000_version01" name="bertram_previte_sherman_kinard_satin_2000_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  
<articleinfo>
  
<title>The Phantom Burster Model For Pancreatic Beta-Cells</title>
  
<author>
    
<firstname>James</firstname>
          
<surname>Lawson</surname>
    
<affiliation>
      
<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    
</affiliation>
  
</author>
</articleinfo>
  
<section id="sec_status">
    
<title>Model Status</title>
    
<para>
            This model has been rebuilt according to the author's original XPPAUT code, which can be found at http://www.math.fsu.edu/~bertram/software/islet/BJ_04a.ode . The model is still unable to produce the appropriate output, however. This file is known to run in PCEnv.
          
</para>
  
</section>
  
<sect1 id="sec_structure">
<title>Model Structure</title>

<para>
Pancreatic beta-cells have been the subject of both experimental and theoretical studies for several decades.  One reason for this interest has been the essential role beta-cells play in glucose homeostasis - they are the only source of insulin that most cells require in order to take up and metabolise glucose, and impairment of beta-cell function contributes to diabetes.  A major focus of theoretical work has been beta-cell dynamics, especially in the form of bursting electrical activity.  The bursts consist of active phases of Ca
<superscript>2+</superscript>-carrying action potentials alternating with silent phases of repolarisation and are accompanied by oscillations in cytosolic Ca<superscript>2+</superscript>, which drive pulses of insulin secretion.  
</para>

<para>
Experimentally, electrical activity in beta-cells is studied in two distinct preparations: islets of Langerhans, which are microorgans containing thousands of endocrine cells, and isolated cells.  Pancreatic beta-cells exhibit bursting oscillations with a wide range of periods.  Whereas periods in isolated cells are generally either a few seconds or a few minutes, in intact islets of Langerhans they are intermediate (10-60 seconds).  In their 2000 publication, Richard Bertram, Joseph Previte, Arthur Sherman, Tracie A. Kinard and Leslie S. Satin develop a mathematical model for beta-cell electrical activity capable of generating this wide range of bursting oscillations.  Unlike previously published models, bursting is driven by the interaction of two slow processes (I
<subscript>s1</subscript> and I<subscript>s2</subscript> in <xref linkend="fig_cell_diagram" /> below), one with a relatively small time constant (1-5 seconds) and the other with a much larger time constant (1-2 minutes).  Bursting on the intermediate time scale is generated without the need for a slow process having an intermediate time constant, hence phantom bursting.  This mathematical model has been translated into a CellML description which can be downloaded in various formats as described in <xref linkend="sec_download_this_model" />.
</para>

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

<para>
<ulink url="http://www.biophysj.org/cgi/content/abstract/79/6/2880">The Phantom Burster Model for Pancreatic beta-Cells</ulink>, Richard Bertram, Joseph Previte, Arthur Sherman, Tracie A. Kinard and Leslie S. Satin, 2000, <ulink url="http://www.biophysj.org/">
            
<emphasis>Biophysical Journal</emphasis>
          
</ulink>, 79, 2880-2892.  (<ulink url="http://www.biophysj.org/cgi/content/full/79/6/2880">Full text</ulink> and <ulink url="http://www.biophysj.org/cgi/reprint/79/6/2880.pdf">PDF</ulink> versions of the article are available for Journal Members on the Biophysical Journal website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=11106596&amp;dopt=Abstract">PubMed ID: 11106596</ulink>  
</para>


<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  
<imageobject>
    
<objectinfo>
      
<title>cell schematic for the model</title>
    
</objectinfo>
    
<imagedata fileref="../images/bertram_model_2000/bertram_2000.png" />
  
</imageobject>
</mediaobject>
<caption>Schematic diagram of the pancreatic beta-cell plasma membrane showing the ionic currents captured by the phantom burster model.</caption>
</informalfigure>

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


  
<units name="minute">
    
<unit units="second" multiplier="60.0" />
  
</units>
  
<units name="femtoA">
    
<unit units="ampere" prefix="femto" />
  
</units>
  
<units name="femtoF">
    
<unit units="farad" prefix="femto" />
  
</units>
  
<units name="millivolt">
    
<unit units="volt" prefix="milli" />
  
</units>
  
<units name="picoS">
    
<unit units="siemens" prefix="pico" />
  
</units>
  
<component name="environment">
    
<variable units="millisecond" public_interface="out" name="time" />
  
</component>
  
<component name="membrane">
    
<variable units="millisecond" public_interface="in" name="time" />
    
<variable units="femtoF" name="Cm" initial_value="4524" />
    
<variable units="millivolt" public_interface="out" name="V" initial_value="-43" />
    
<variable units="femtoA" public_interface="in" name="ICa" />
    
<variable units="femtoA" public_interface="in" name="IK" />
    
<variable units="femtoA" public_interface="in" name="Il" />
    
<variable units="femtoA" public_interface="in" name="Is1" />
    
<variable units="femtoA" public_interface="in" name="Is2" />
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="membrane d(V)/d(time)">
      
<apply>
        
<eq />
        
<apply>
          
<diff />
          
<bvar>
            
<ci>time</ci>
          
</bvar>
          
<ci>V</ci>
        
</apply>
        
<apply>
          
<divide />
          
<apply>
            
<minus />
            
<apply>
              
<plus />
              
<ci>ICa</ci>
              
<ci>IK</ci>
              
<ci>Il</ci>
              
<ci>Is1</ci>
              
<ci>Is2</ci>
            
</apply>
          
</apply>
          
<ci>Cm</ci>
        
</apply>
      
</apply>
    
</math>
  
</component>
  
<component name="Ca_current">
    
<variable units="millisecond" public_interface="in" name="time" />
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="minf">
      
<apply>
        
<eq />
        
<ci>minf</ci>
        
<apply>
          
<divide />
          
<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless">1</cn>
          
<apply>
            
<plus />
            
<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless">1</cn>
            
<apply>
              
<exp />
              
<apply>
                
<divide />
                
<apply>
                  
<minus />
                  
<ci>Vm</ci>
                  
<ci>V</ci>
                
</apply>
                
<ci>sm</ci>
              
</apply>
            
</apply>
          
</apply>
        
</apply>
      
</apply>
    
</math>
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="ICa">
      
<apply>
        
<eq />
        
<ci>ICa</ci>
        
<apply>
          
<times />
          
<ci>gCa</ci>
          
<ci>minf</ci>
          
<apply>
            
<minus />
            
<ci>V</ci>
            
<ci>VCa</ci>
          
</apply>
        
</apply>
      
</apply>
    
</math>
    
<variable units="millivolt" public_interface="in" name="V" />
    
<variable units="millivolt" name="Vm" initial_value="-22" />
    
<variable units="millivolt" name="VCa" initial_value="100" />
    
<variable units="picoS" name="gCa" initial_value="280" />
    
<variable units="dimensionless" name="minf" />
    
<variable units="millivolt" name="sm" initial_value="7.5" />
    
<variable units="femtoA" public_interface="out" name="ICa" />
  
</component>
  
<component name="rapid_K_current">
    
<variable units="millisecond" public_interface="in" name="time" />
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="IK">
      
<apply>
        
<eq />
        
<ci>IK</ci>
        
<apply>
          
<times />
          
<ci>gK</ci>
          
<ci>n</ci>
          
<apply>
            
<minus />
            
<ci>V</ci>
            
<ci>VK</ci>
          
</apply>
        
</apply>
      
</apply>
    
</math>
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="d(n)/d(time)">
      
<apply>
        
<eq />
        
<apply>
          
<diff />
          
<bvar>
            
<ci>time</ci>
          
</bvar>
          
<ci>n</ci>
        
</apply>
        
<apply>
          
<divide />
          
<apply>
            
<times />
            
<ci>lambda</ci>
            
<apply>
              
<minus />
              
<ci>ninf</ci>
              
<ci>n</ci>
            
</apply>
          
</apply>
          
<apply>
            
<times />
            
<ci>taun</ci>
            
<ci>V</ci>
          
</apply>
        
</apply>
      
</apply>
    
</math>
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="ninf">
      
<apply>
        
<eq />
        
<ci>ninf</ci>
        
<apply>
          
<divide />
          
<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless">1</cn>
          
<apply>
            
<plus />
            
<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" 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>
    
</math>
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="taun">
      
<apply>
        
<eq />
        
<ci>taun</ci>
        
<apply>
          
<divide />
          
<ci>tnbar</ci>
          
<apply>
            
<plus />
            
<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" 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>
    
<variable units="femtoA" public_interface="out" name="IK" />
    
<variable units="millivolt" public_interface="out" name="VK" initial_value="-80" />
    
<variable units="millivolt" public_interface="in" name="V" />
    
<variable units="picoS" name="gK" initial_value="1300" />
    
<variable units="millivolt" name="n" initial_value="0.03" />
    
<variable units="dimensionless" name="lambda" initial_value="1.1" />
    
<variable units="dimensionless" name="tnbar" initial_value="9.09" />
    
<variable units="millivolt" name="Vn" initial_value="-9" />
    
<variable units="millivolt" name="sn" initial_value="10" />
    
<variable units="dimensionless" name="taun" />
    
<variable units="dimensionless" name="ninf" />
  
</component>
  
<component name="slow_K_current">
    
<variable units="millisecond" public_interface="in" name="time" />
    
<variable units="femtoA" public_interface="out" name="Is1" />
    
<variable units="picoS" name="gs1" initial_value="20" />
    
<variable units="dimensionless" name="s1" initial_value="0.1" />
    
<variable units="millivolt" public_interface="in" name="V" />
    
<variable units="millivolt" public_interface="in" name="VK" />
    
<variable units="dimensionless" name="s1inf" />
    
<variable units="millivolt" name="Vs1" initial_value="-40" />
    
<variable units="dimensionless" name="ss1" initial_value="0.5" />
    
<variable units="dimensionless" name="taus1" initial_value="1000" />
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="Is1">
      
<apply>
        
<eq />
        
<ci>Is1</ci>
        
<apply>
          
<times />
          
<ci>gs1</ci>
          
<ci>s1</ci>
          
<apply>
            
<minus />
            
<ci>V</ci>
            
<ci>VK</ci>
          
</apply>
        
</apply>
      
</apply>
    
</math>
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="s1inf">
      
<apply>
        
<eq />
        
<ci>s1inf</ci>
        
<apply>
          
<divide />
          
<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless">1</cn>
          
<apply>
            
<plus />
            
<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless">1</cn>
            
<apply>
              
<exp />
              
<apply>
                
<divide />
                
<apply>
                  
<minus />
                  
<ci>Vs1</ci>
                  
<ci>V</ci>
                
</apply>
                
<ci>ss1</ci>
              
</apply>
            
</apply>
          
</apply>
        
</apply>
      
</apply>
    
</math>
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="d(s1)/d(time)">
      
<apply>
        
<eq />
        
<apply>
          
<diff />
          
<bvar>
            
<ci>time</ci>
          
</bvar>
          
<ci>s1</ci>
        
</apply>
        
<apply>
          
<divide />
          
<apply>
            
<minus />
            
<ci>s1inf</ci>
            
<ci>s1</ci>
          
</apply>
          
<ci>taus1</ci>
        
</apply>
      
</apply>
    
</math>
  
</component>
  
<component name="very_slow_K_current">
    
<variable units="millisecond" public_interface="in" name="time" />
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="Is2">
      
<apply>
        
<eq />
        
<ci>Is2</ci>
        
<apply>
          
<times />
          
<ci>gs2</ci>
          
<ci>s2</ci>
          
<apply>
            
<minus />
            
<ci>V</ci>
            
<ci>VK</ci>
          
</apply>
        
</apply>
      
</apply>
    
</math>
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="s2inf">
      
<apply>
        
<eq />
        
<ci>s2inf</ci>
        
<apply>
          
<divide />
          
<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless">1</cn>
          
<apply>
            
<plus />
            
<cn xmlns:cellml="http://www.cellml.org/cellml/1.1#" cellml:units="dimensionless">1</cn>
            
<apply>
              
<exp />
              
<apply>
                
<divide />
                
<apply>
                  
<minus />
                  
<ci>Vs2</ci>
                  
<ci>V</ci>
                
</apply>
                
<ci>ss2</ci>
              
</apply>
            
</apply>
          
</apply>
        
</apply>
      
</apply>
    
</math>
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="d(s2)/d(time)">
      
<apply>
        
<eq />
        
<apply>
          
<diff />
          
<bvar>
            
<ci>time</ci>
          
</bvar>
          
<ci>s2</ci>
        
</apply>
        
<apply>
          
<divide />
          
<apply>
            
<minus />
            
<ci>s2inf</ci>
            
<ci>s2</ci>
          
</apply>
          
<ci>taus2</ci>
        
</apply>
      
</apply>
    
</math>
    
<variable units="millivolt" public_interface="in" name="V" />
    
<variable units="millivolt" public_interface="in" name="VK" />
    
<variable units="millivolt" name="Vs2" initial_value="-42" />
    
<variable units="dimensionless" name="s2" initial_value="0.434" />
    
<variable units="dimensionless" name="s2inf" />
    
<variable units="dimensionless" name="ss2" initial_value="0.4" />
    
<variable units="picoS" name="gs2" initial_value="32" />
    
<variable units="dimensionless" name="taus2" initial_value="120000" />
    
<variable units="femtoA" public_interface="out" name="Is2" />
  
</component>
  
<component name="leak_current">
    
<variable units="millisecond" public_interface="in" name="time" />
    
<variable units="femtoA" public_interface="out" name="Il" />
    
<variable units="picoS" name="gl" initial_value="25" />
    
<math xmlns="http://www.w3.org/1998/Math/MathML" id="Il">
      
<apply>
        
<eq />
        
<ci>Il</ci>
        
<apply>
          
<times />
          
<ci>gl</ci>
          
<apply>
            
<minus />
            
<ci>V</ci>
            
<ci>Vl</ci>
          
</apply>
        
</apply>
      
</apply>
    
</math>
    
<variable units="millivolt" name="Vl" initial_value="-40" />
    
<variable units="millivolt" public_interface="in" name="V" />
  
</component>
  
<units name="millisecond">
    
<unit units="second" prefix="milli" />
  
</units>
  
<connection>
    
<map_components component_2="membrane" component_1="environment" />
    
<map_variables variable_2="time" variable_1="time" />
  
</connection>
  
<connection>
    
<map_components component_2="membrane" component_1="Ca_current" />
    
<map_variables variable_2="ICa" variable_1="ICa" />
  
</connection>
  
<connection>
    
<map_components component_2="membrane" component_1="rapid_K_current" />
    
<map_variables variable_2="IK" variable_1="IK" />
  
</connection>
  
<connection>
    
<map_components component_2="membrane" component_1="slow_K_current" />
    
<map_variables variable_2="Is1" variable_1="Is1" />
  
</connection>
  
<connection>
    
<map_components component_2="membrane" component_1="very_slow_K_current" />
    
<map_variables variable_2="Is2" variable_1="Is2" />
  
</connection>
  
<connection>
    
<map_components component_2="membrane" component_1="leak_current" />
    
<map_variables variable_2="Il" variable_1="Il" />
  
</connection>
  
<connection>
    
<map_components component_2="Ca_current" component_1="environment" />
    
<map_variables variable_2="time" variable_1="time" />
  
</connection>
  
<connection>
    
<map_components component_2="Ca_current" component_1="membrane" />
    
<map_variables variable_2="V" variable_1="V" />
  
</connection>
  
<connection>
    
<map_components component_2="rapid_K_current" component_1="environment" />
    
<map_variables variable_2="time" variable_1="time" />
  
</connection>
  
<connection>
    
<map_components component_2="rapid_K_current" component_1="membrane" />
    
<map_variables variable_2="V" variable_1="V" />
  
</connection>
  
<connection>
    
<map_components component_2="slow_K_current" component_1="environment" />
    
<map_variables variable_2="time" variable_1="time" />
  
</connection>
  
<connection>
    
<map_components component_2="slow_K_current" component_1="membrane" />
    
<map_variables variable_2="V" variable_1="V" />
  
</connection>
  
<connection>
    
<map_components component_2="slow_K_current" component_1="rapid_K_current" />
    
<map_variables variable_2="VK" variable_1="VK" />
  
</connection>
  
<connection>
    
<map_components component_2="very_slow_K_current" component_1="environment" />
    
<map_variables variable_2="time" variable_1="time" />
  
</connection>
  
<connection>
    
<map_components component_2="very_slow_K_current" component_1="membrane" />
    
<map_variables variable_2="V" variable_1="V" />
  
</connection>
  
<connection>
    
<map_components component_2="very_slow_K_current" component_1="rapid_K_current" />
    
<map_variables variable_2="VK" variable_1="VK" />
  
</connection>
  
<connection>
    
<map_components component_2="leak_current" component_1="environment" />
    
<map_variables variable_2="time" variable_1="time" />
  
</connection>
  
<connection>
    
<map_components component_2="leak_current" component_1="membrane" />
    
<map_variables variable_2="V" variable_1="V" />
  
</connection>


  


<rdf:RDF>
  
<rdf:Seq rdf:about="rdf:#6637b1a4-d8eb-4f9a-b85f-69d1b36f2000">
    
<rdf:li rdf:resource="rdf:#35030caf-b595-49bf-8d36-6f08a0ae0760" />
    
<rdf:li rdf:resource="rdf:#b182afc4-5f51-40d2-9d3e-470af7fc09ce" />
    
<rdf:li rdf:resource="rdf:#df17b235-545a-48d7-a023-0e3e9a2cc876" />
    
<rdf:li rdf:resource="rdf:#f55923fe-5f94-49c1-ba6b-ec70b2ec8c37" />
    
<rdf:li rdf:resource="rdf:#e59cb78c-d1cf-4837-b200-2fe37a586c91" />
  
</rdf:Seq>
  
<rdf:Description rdf:about="rdf:#7ccadafc-8e12-4344-a927-c57d81ca16ef">
    
<dcterms:W3CDTF>2000-12</dcterms:W3CDTF>
  
</rdf:Description>
  
<rdf:Description rdf:about="rdf:#b205fcbd-6643-41f4-b451-9d5e34e45134">
    
<vCard:FN>Catherine Lloyd</vCard:FN>
  
</rdf:Description>
  
<rdf:Description rdf:about="">
    
<dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    
<cmeta:comment rdf:resource="rdf:#47af5cd5-9c00-4fda-ac93-d47f15cd02c9" />
    
<dcterms:created rdf:resource="rdf:#b292721c-0356-46e8-ad68-2b9594fa3c28" />
    
<dc:creator rdf:resource="rdf:#2eee98ba-4648-4ae4-91b0-3eb277147f7a" />
    
<cmeta:modification rdf:resource="rdf:#73e9e190-ac5b-463c-be82-ffc62ca74630" />
    
<cmeta:modification rdf:resource="rdf:#f4225862-d951-4cdf-bd9c-ca17d407fd11" />
  
</rdf:Description>
  
<rdf:Description rdf:about="rdf:#d06465f1-e6e2-43a4-a4b9-357ab90d210e">
    
<vCard:Given>Richard</vCard:Given>
    
<vCard:Family>Bertram</vCard:Family>
  
</rdf:Description>
  
<rdf:Description rdf:about="rdf:#47af5cd5-9c00-4fda-ac93-d47f15cd02c9">
    
<dc:creator rdf:resource="rdf:#04ee5046-2c1a-4d01-b955-40bf86c577cd" />
    
<rdf:value>This model has been rebuilt according to the author's original XPPAUT code, which can be found at http://www.math.fsu.edu/~bertram/software/islet/BJ_04a.ode . Unfortunately the model is still unable to produce the correct output. This file is known to run in PCEnv.</rdf:value>
  
</rdf:Description>
  
<rdf:Description rdf:about="rdf:#f4225862-d951-4cdf-bd9c-ca17d407fd11">
    
<dcterms:modified rdf:resource="rdf:#41b0cfc0-f5f9-41d4-aff4-90d30c1c2c8f" />
    
<rdf:value>
          Added more metadata.
        
</rdf:value>
    
<cmeta:modifier rdf:resource="rdf:#3ca11ded-a806-4f6e-872c-99fbf0f58600" />
  
</rdf:Description>
  
<rdf:Description rdf:about="rdf:#a049c1f5-2c5a-4be0-96b8-d394a9fd496a">
    
<vCard:Given>Autumn</vCard:Given>
    
<vCard:Family>Cuellar</vCard:Family>
    
<vCard:Other>A</vCard:Other>
  
</rdf:Description>
  
<rdf:Description rdf:about="rdf:#0640428e-01f5-49b8-8db5-f0a8b27f0698">
    
<dc:creator rdf:resource="rdf:#b205fcbd-6643-41f4-b451-9d5e34e45134" />
    
<rdf:value>This is the CellML description of Bertram et al's 2000 phantom burster           model for pancreatic beta-cells.</rdf:value>
  
</rdf:Description>
  
<rdf:Description rdf:about="rdf:#b672a99d-8d3d-475d-b5f3-1ef0b1e09fea">
    
<vCard:Given>James</vCard:Given>
    
<vCard:Family>Lawson</vCard:Family>
    
<vCard:Other>Richard</vCard:Other>
  
</rdf:Description>
  
<rdf:D