Location: Snyed, Wetton, Charles, Sanderson, 1995 @ 3913942bd086 / Sneyd_s_Model.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2009-09-29 15:38:56+13:00
Desc:
Added files to workspace for the first time.
Permanent Source URI:
http://models.cellml.org/workspace/sneyd_wetton_charles_sanderson_1995/rawfile/3913942bd08624a4107673042804900e38cf35a6/Sneyd_s_Model.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 15/09/2009 at 8:43:35 at p.m. using:

COR (0.9.31.1319)
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 name="wliu052_part_iv_project" cmeta:id="wliu052_part_iv_project" 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>
      <section id="sec_status">
        <title>Model Status</title>
        <para>
	This CellML model is the model which was used to produce the original results in the paper. The units have been checked and are consistent and this model is known to be valid CellML.
          </para>
      </section>
      <sect1 id="sec_structure">
        <title>Model Structure</title>

        <para>
        The present of IP3 in the cytosol induces the calcium release from ER by binding with the IP3 receptor on the ER surface, the calcium concentration is also depends on the ER membrane calcium pump acitivity which pump Ca2+ back into ER and the leakage of Ca2+ from ER.
        </para>
	<para>
	This model of the calcium oscillation in induced by intracellular IP3 concentration is based on the original paper cited below:
	</para>
	<para>
	Intercellular calcium waves mediated by diffusion of inostol trisphosphate: a two dimensional model. 7.	J. Sneyd, B.T.R. Wetton, A.C. Charles and M.J. Sanderson. 1995 <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/7611375">PubMed ID: 7611375</ulink>
	</para>

<informalfigure float="0" id="fig_reaction_diagram">
	<mediaobject>
		<imageobject>
			<objectinfo>
			<title>Diagram of model</title>
			</objectinfo>
			<imagedata fileref="Sneyd_s_Model.png"/>
		</imageobject>
	</mediaobject>
<caption>A schematic diagram of the calcium oscillation induced by intracellular IP3 is described by the mathematical model.</caption>
</informalfigure>

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

   <units name="micromolar">
      <unit units="mole" prefix="micro"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="flux">
      <unit units="micromolar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_sec_square">
      <unit units="second" exponent="-2"/>
   </units>
   <units name="micromolar_per_second">
      <unit units="micromolar"/>
      <unit units="per_second"/>
   </units>
   <units name="micromolar_square_per_second">
      <unit units="micromolar" exponent="2"/>
      <unit units="per_second"/>
   </units>
   <component name="enviroment">
      <variable name="time" units="second" public_interface="out"/>
   </component>
   <component name="P">
      <variable name="P" units="micromolar" initial_value="0" public_interface="out"/>
      <variable name="V_p" units="per_second" initial_value="0.08"/>
      <variable name="k_p" units="micromolar" initial_value="1"/>
      <variable name="IPR_3_flux" units="flux" initial_value="0.72"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P</ci>
            </apply>
            <piecewise>
               <piece>
                  <apply>
                     <minus/>
                     <ci>IPR_3_flux</ci>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>V_p</ci>
                           <ci>P</ci>
                           <ci>k_p</ci>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>k_p</ci>
                           <ci>P</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <leq/>
                     <ci>time</ci>
                     <cn cellml:units="second">15</cn>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>V_p</ci>
                        </apply>
                        <ci>P</ci>
                        <ci>k_p</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>k_p</ci>
                        <ci>P</ci>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="c">
      <variable name="c" units="micromolar" initial_value="0.3" public_interface="out"/>
      <variable name="J_flux" units="flux" public_interface="in"/>
      <variable name="J_pump" units="flux" public_interface="in"/>
      <variable name="J_leak" units="flux" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>c</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <ci>J_flux</ci>
                  <ci>J_pump</ci>
               </apply>
               <ci>J_leak</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="J_flux">
      <variable name="J_flux" units="flux" public_interface="out"/>
      <variable name="k_flux" units="micromolar_per_second" initial_value="3"/>
      <variable name="mu" units="dimensionless" public_interface="in"/>
      <variable name="h" units="dimensionless" public_interface="in"/>
      <variable name="b" units="dimensionless" initial_value="0.11"/>
      <variable name="c" units="micromolar" public_interface="in"/>
      <variable name="k_1" units="micromolar" initial_value="0.7"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_flux</ci>
            <apply>
               <times/>
               <ci>k_flux</ci>
               <ci>mu</ci>
               <ci>h</ci>
               <apply>
                  <plus/>
                  <ci>b</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <ci>b</ci>
                        </apply>
                        <ci>c</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>k_1</ci>
                        <ci>c</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="J_pump">
      <variable name="J_pump" units="flux" public_interface="out"/>
      <variable name="c" units="micromolar" public_interface="in"/>
      <variable name="gamma" units="micromolar_per_second" initial_value="1"/>
      <variable name="k_gamma" units="micromolar" initial_value="0.27"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_pump</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>gamma</ci>
                  <apply>
                     <power/>
                     <ci>c</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>k_gamma</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>c</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="J_leak">
      <variable name="J_leak" units="flux" public_interface="out"/>
      <variable name="beta" units="flux" initial_value="0.15"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>J_leak</ci>
            <ci>beta</ci>
         </apply>
      </math>
   </component>
   <component name="mu">
      <variable name="mu" units="dimensionless" public_interface="out"/>
      <variable name="P" units="micromolar" public_interface="in"/>
      <variable name="k_mu" units="micromolar" initial_value="0.01"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>mu</ci>
            <apply>
               <divide/>
               <apply>
                  <power/>
                  <ci>P</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>k_mu</ci>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>P</ci>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="h">
      <variable name="h" units="dimensionless" initial_value="1" public_interface="out"/>
      <variable name="c" units="micromolar" public_interface="in"/>
      <variable name="k_2" units="micromolar" initial_value="0.7"/>
      <variable name="tau_h" units="second" initial_value="0.2"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <apply>
                        <power/>
                        <ci>k_2</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <power/>
                           <ci>k_2</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <power/>
                           <ci>c</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <ci>h</ci>
               </apply>
               <ci>tau_h</ci>
            </apply>
         </apply>
      </math>
   </component>
   <connection>
      <map_components component_1="enviroment" component_2="P"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="enviroment" component_2="c"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="c" component_2="J_flux"/>
      <map_variables variable_1="J_flux" variable_2="J_flux"/>
      <map_variables variable_1="c" variable_2="c"/>
   </connection>
   <connection>
      <map_components component_1="c" component_2="J_pump"/>
      <map_variables variable_1="J_pump" variable_2="J_pump"/>
      <map_variables variable_1="c" variable_2="c"/>
   </connection>
   <connection>
      <map_components component_1="c" component_2="J_leak"/>
      <map_variables variable_1="J_leak" variable_2="J_leak"/>
   </connection>
   <connection>
      <map_components component_1="mu" component_2="J_flux"/>
      <map_variables variable_1="mu" variable_2="mu"/>
   </connection>
   <connection>
      <map_components component_1="h" component_2="J_flux"/>
      <map_variables variable_1="h" variable_2="h"/>
   </connection>
   <connection>
      <map_components component_1="mu" component_2="P"/>
      <map_variables variable_1="P" variable_2="P"/>
   </connection>
   <connection>
      <map_components component_1="h" component_2="c"/>
      <map_variables variable_1="c" variable_2="c"/>
   </connection>
   <connection>
      <map_components component_1="h" component_2="enviroment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
</model>