Location: Wang, Huang, Huang, 2006 @ 32281c9fc6ad / Wang_s_Model.cellml

Author:
Wei Liu <devnull@localhost>
Date:
2009-09-18 11:26:16+12:00
Desc:
Inital Model
Permanent Source URI:
https://models.cellml.org/workspace/wang_huang_huang_2006/rawfile/32281c9fc6add9698fc26075a63f5796f7697fb7/Wang_s_Model.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 21/08/2009 at 9:19:47 at a.m. using:

COR (0.9.31.1309)
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="Calcium_Model_2" cmeta:id="Calcium_Model_2" 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>
        Abstract:
        </para>
        <para>
        A quantitative kinetic model is proposed to simulate the ATP-induced intracellular Ca(2+) oscillations. The quantitative effect of ATP concentration upon the oscillations was successfully simulated. Our simulation results support previous experimental explanations that the Ca(2+) oscillations are mainly due to interaction of Ca(2+) release from the endoplasmic reticulum (ER) and the ATP-dependent Ca(2+) pump back into the ER, and the oscillations are prolonged by extracellular Ca(2+) entry that maintains the constant Ca(2+) supplies to its intracellular stores. The model is also able to simulate the sudden disappearance phenomenon of the Ca(2+) oscillations observed in some cell types by taking into account of the biphasic characteristic of the Ca(2+) release from the endoplasmic reticulum (ER). Moreover, the model simulation results for the Ca(2+) oscillations characteristics such as duration, peak [Ca(2+)](cyt), and average interval, etc., lead to prediction of some possible factors responsible for the variations of Ca(2+) oscillations in different types of cells.
        </para>

        <para>

        </para>
	<para>
	This model of the calcium oscillation induced by ATP binding P2Y receptor activity is based on the original paper cited below:
	</para>
	<para>
	A quantitative kinetic model for ATP-induced intracellular Ca2+ oscillations. . Wang, X. Huang and W. Huang. 2006 <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/17188305">PubMed ID: 17188305</ulink>
	</para>

<informalfigure float="0" id="fig_reaction_diagram">
	<mediaobject>
		<imageobject>
			<objectinfo>
			<title>Diagram of model</title>
			</objectinfo>
			<imagedata fileref="Wang_s_Model.png"/>
		</imageobject>
	</mediaobject>
<caption>A schematic diagram of the calcium oscillation induced by extracellular ATP binding on surface P2Y receptor 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="per_micromolar">
      <unit units="micromolar" exponent="-1"/>
   </units>
   <units name="per_micromole">
      <unit units="mole" prefix="micro" 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_micromolar_per_second">
      <unit units="micromolar" exponent="-1"/>
      <unit units="per_second"/>
   </units>
   <units name="micromolar_square_per_second">
      <unit units="micromolar" exponent="2"/>
      <unit units="per_second"/>
   </units>
   <component name="environment">
      <variable name="time" units="second" public_interface="out"/>
   </component>
   <component name="Constants">
      <variable name="k_0" units="flux" initial_value="1e-4" public_interface="out"/>
      <variable name="k_1" units="per_second" initial_value="3.4" public_interface="out"/>
      <variable name="k_2" units="per_second" initial_value="4" public_interface="out"/>
      <variable name="k_3" units="per_second" initial_value="4.5" public_interface="out"/>
      <variable name="k_4" units="per_second" initial_value="1.2" public_interface="out"/>
      <variable name="k_5" units="per_second" initial_value="0.12" public_interface="out"/>
      <variable name="k_6" units="per_second" initial_value="14" public_interface="out"/>
      <variable name="k_7" units="per_second" initial_value="2" public_interface="out"/>
      <variable name="k_8" units="flux" initial_value="10.5" public_interface="out"/>
      <variable name="k_9" units="flux" initial_value="0.6" public_interface="out"/>
      <variable name="k_10" units="flux" initial_value="3" public_interface="out"/>
      <variable name="k_11" units="flux" initial_value="0.26" public_interface="out"/>
      <variable name="C_PLC_T" units="micromolar" initial_value="0.01" public_interface="out"/>
      <variable name="K_D" units="micromolar" initial_value="0.01" public_interface="out"/>
      <variable name="K_P" units="micromolar" initial_value="0.004" public_interface="out"/>
      <variable name="K_R" units="micromolar" initial_value="0.2" public_interface="out"/>
      <variable name="K_G" units="micromolar" initial_value="0.025" public_interface="out"/>
      <variable name="K_S" units="micromolar" initial_value="0.025" public_interface="out"/>
      <variable name="K_ER" units="micromolar" initial_value="0.075" public_interface="out"/>
      <variable name="K_C1" units="micromolar" initial_value="1" public_interface="out"/>
      <variable name="K_C2" units="micromolar" initial_value="2" public_interface="out"/>
      <variable name="beta" units="dimensionless" initial_value="0.05" public_interface="out"/>
      <variable name="lambda" units="dimensionless" initial_value="0.001" public_interface="out"/>
      <variable name="rho" units="dimensionless" initial_value="0.2" public_interface="out"/>
      <variable name="n" units="dimensionless" initial_value="4" public_interface="out"/>
      <variable name="m" units="dimensionless" initial_value="2" public_interface="out"/>
      <variable name="w" units="dimensionless" initial_value="3" public_interface="out"/>
   </component>
   <component name="R_values">
      <variable name="R_APLC" units="dimensionless" public_interface="out"/>
      <variable name="R_PKC" units="dimensionless" public_interface="out"/>
      <variable name="R_G" units="dimensionless" public_interface="out"/>
      <variable name="R_DG" units="dimensionless" public_interface="out"/>
      <variable name="R_IP_3" units="dimensionless" public_interface="out"/>
      <variable name="R_Cyt1" units="dimensionless" public_interface="out"/>
      <variable name="R_Cyt2" units="dimensionless" public_interface="out"/>
      <variable name="R_ER" units="dimensionless" public_interface="out"/>
      <variable name="APLC" units="micromolar" public_interface="in"/>
      <variable name="DG" units="micromolar" public_interface="in"/>
      <variable name="C_cyt" units="micromolar" public_interface="in"/>
      <variable name="G" units="micromolar" public_interface="in"/>
      <variable name="IP_3" units="micromolar" public_interface="in"/>
      <variable name="C_ER" units="micromolar" public_interface="in"/>
      <variable name="K_P" units="micromolar" public_interface="in"/>
      <variable name="K_D" units="micromolar" public_interface="in"/>
      <variable name="K_R" units="micromolar" public_interface="in"/>
      <variable name="K_G" units="micromolar" public_interface="in"/>
      <variable name="K_S" units="micromolar" public_interface="in"/>
      <variable name="K_C1" units="micromolar" public_interface="in"/>
      <variable name="K_C2" units="micromolar" public_interface="in"/>
      <variable name="K_ER" units="micromolar" public_interface="in"/>
      <variable name="n" units="dimensionless" public_interface="in"/>
      <variable name="m" units="dimensionless" public_interface="in"/>
      <variable name="w" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>R_APLC</ci>
            <apply>
               <divide/>
               <ci>APLC</ci>
               <apply>
                  <plus/>
                  <ci>K_P</ci>
                  <ci>APLC</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>R_PKC</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>DG</ci>
                     <apply>
                        <plus/>
                        <ci>K_D</ci>
                        <ci>DG</ci>
                     </apply>
                  </apply>
                  <ci>C_cyt</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>K_R</ci>
                  <ci>C_cyt</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>R_G</ci>
            <apply>
               <divide/>
               <apply>
                  <power/>
                  <ci>G</ci>
                  <ci>n</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>K_G</ci>
                     <ci>n</ci>
                  </apply>
                  <apply>
                     <power/>
                     <ci>G</ci>
                     <ci>n</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>R_DG</ci>
            <apply>
               <divide/>
               <apply>
                  <power/>
                  <ci>DG</ci>
                  <ci>m</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>K_D</ci>
                     <ci>m</ci>
                  </apply>
                  <apply>
                     <power/>
                     <ci>DG</ci>
                     <ci>m</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>R_IP_3</ci>
            <apply>
               <divide/>
               <apply>
                  <power/>
                  <ci>IP_3</ci>
                  <cn cellml:units="dimensionless">3</cn>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>K_S</ci>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>IP_3</ci>
                     <cn cellml:units="dimensionless">3</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>R_Cyt1</ci>
            <apply>
               <divide/>
               <ci>C_cyt</ci>
               <apply>
                  <plus/>
                  <ci>K_C1</ci>
                  <ci>C_cyt</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>R_Cyt2</ci>
            <apply>
               <divide/>
               <ci>C_cyt</ci>
               <apply>
                  <plus/>
                  <ci>K_C2</ci>
                  <ci>C_cyt</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>R_ER</ci>
            <apply>
               <divide/>
               <apply>
                  <power/>
                  <ci>C_ER</ci>
                  <ci>w</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>K_ER</ci>
                     <ci>w</ci>
                  </apply>
                  <apply>
                     <power/>
                     <ci>C_ER</ci>
                     <ci>w</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="G_GTP">
      <variable name="G" units="micromolar" initial_value="0.001" public_interface="out"/>
      <variable name="k_0" units="flux" public_interface="in"/>
      <variable name="k_1" units="per_second" public_interface="in"/>
      <variable name="k_2" units="per_second" public_interface="in"/>
      <variable name="k_3" units="per_second" public_interface="in"/>
      <variable name="R_APLC" units="dimensionless" public_interface="in"/>
      <variable name="R_PKC" units="dimensionless" 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>G</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <ci>k_0</ci>
                     <apply>
                        <times/>
                        <ci>k_1</ci>
                        <ci>G</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_2</ci>
                     <ci>R_APLC</ci>
                     <ci>G</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_3</ci>
                  <ci>R_PKC</ci>
                  <ci>G</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="APLC">
      <variable name="APLC" units="micromolar" initial_value="0.001" public_interface="out"/>
      <variable name="PLC" units="micromolar"/>
      <variable name="C_PLC_T" units="micromolar" public_interface="in"/>
      <variable name="k_4" units="per_second" public_interface="in"/>
      <variable name="k_5" units="per_second" public_interface="in"/>
      <variable name="R_G" units="dimensionless" public_interface="in"/>
      <variable name="R_DG" units="dimensionless" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>PLC</ci>
            <apply>
               <minus/>
               <ci>C_PLC_T</ci>
               <ci>APLC</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>APLC</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_4</ci>
                  <ci>R_G</ci>
                  <ci>R_DG</ci>
                  <ci>PLC</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_5</ci>
                  <ci>APLC</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="IP_3">
      <variable name="IP_3" units="micromolar" initial_value="0.001" public_interface="out"/>
      <variable name="k_6" units="per_second" public_interface="in"/>
      <variable name="k_7" units="per_second" public_interface="in"/>
      <variable name="APLC" units="micromolar" 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>IP_3</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_6</ci>
                  <ci>APLC</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_7</ci>
                  <ci>IP_3</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="DG">
      <variable name="DG" units="micromolar" public_interface="out"/>
      <variable name="IP_3" units="micromolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>DG</ci>
            <ci>IP_3</ci>
         </apply>
      </math>
   </component>
   <component name="C_cyt">
      <variable name="C_cyt" units="micromolar" initial_value="0.2" public_interface="out"/>
      <variable name="beta" units="dimensionless" public_interface="in"/>
      <variable name="rho" units="dimensionless" public_interface="in"/>
      <variable name="k_8" units="flux" public_interface="in"/>
      <variable name="k_9" units="flux" public_interface="in"/>
      <variable name="k_10" units="flux" public_interface="in"/>
      <variable name="k_11" units="flux" public_interface="in"/>
      <variable name="R_IP_3" units="dimensionless" public_interface="in"/>
      <variable name="R_ER" units="dimensionless" public_interface="in"/>
      <variable name="R_Cyt1" units="dimensionless" public_interface="in"/>
      <variable name="R_Cyt2" units="dimensionless" 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_cyt</ci>
            </apply>
            <apply>
               <times/>
               <ci>beta</ci>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>rho</ci>
                        <apply>
                           <minus/>
                           <apply>
                              <times/>
                              <ci>k_8</ci>
                              <ci>R_IP_3</ci>
                              <ci>R_ER</ci>
                           </apply>
                           <apply>
                              <times/>
                              <ci>k_9</ci>
                              <ci>R_Cyt1</ci>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>k_10</ci>
                        <ci>R_Cyt2</ci>
                     </apply>
                  </apply>
                  <ci>k_11</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="C_ER">
      <variable name="C_ER" units="micromolar" initial_value="1" public_interface="out"/>
      <variable name="lambda" units="dimensionless" public_interface="in"/>
      <variable name="k_8" units="flux" public_interface="in"/>
      <variable name="k_9" units="flux" public_interface="in"/>
      <variable name="R_IP_3" units="dimensionless" public_interface="in"/>
      <variable name="R_ER" units="dimensionless" public_interface="in"/>
      <variable name="R_Cyt1" units="dimensionless" 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_ER</ci>
            </apply>
            <apply>
               <times/>
               <ci>lambda</ci>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <ci>k_8</ci>
                     </apply>
                     <ci>R_IP_3</ci>
                     <ci>R_ER</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>k_9</ci>
                     <ci>R_Cyt1</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <connection>
      <map_components component_1="environment" component_2="G_GTP"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="APLC"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="IP_3"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="C_cyt"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="C_ER"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="R_values" component_2="Constants"/>
      <map_variables variable_1="K_P" variable_2="K_P"/>
      <map_variables variable_1="K_D" variable_2="K_D"/>
      <map_variables variable_1="K_G" variable_2="K_G"/>
      <map_variables variable_1="K_R" variable_2="K_R"/>
      <map_variables variable_1="K_S" variable_2="K_S"/>
      <map_variables variable_1="K_ER" variable_2="K_ER"/>
      <map_variables variable_1="K_C1" variable_2="K_C1"/>
      <map_variables variable_1="K_C2" variable_2="K_C2"/>
      <map_variables variable_1="n" variable_2="n"/>
      <map_variables variable_1="m" variable_2="m"/>
      <map_variables variable_1="w" variable_2="w"/>
   </connection>
   <connection>
      <map_components component_1="R_values" component_2="G_GTP"/>
      <map_variables variable_1="G" variable_2="G"/>
      <map_variables variable_1="R_APLC" variable_2="R_APLC"/>
      <map_variables variable_1="R_PKC" variable_2="R_PKC"/>
   </connection>
   <connection>
      <map_components component_1="R_values" component_2="APLC"/>
      <map_variables variable_1="APLC" variable_2="APLC"/>
      <map_variables variable_1="R_G" variable_2="R_G"/>
      <map_variables variable_1="R_DG" variable_2="R_DG"/>
   </connection>
   <connection>
      <map_components component_1="R_values" component_2="IP_3"/>
      <map_variables variable_1="IP_3" variable_2="IP_3"/>
   </connection>
   <connection>
      <map_components component_1="R_values" component_2="DG"/>
      <map_variables variable_1="DG" variable_2="DG"/>
   </connection>
   <connection>
      <map_components component_1="R_values" component_2="C_cyt"/>
      <map_variables variable_1="C_cyt" variable_2="C_cyt"/>
      <map_variables variable_1="R_IP_3" variable_2="R_IP_3"/>
      <map_variables variable_1="R_ER" variable_2="R_ER"/>
      <map_variables variable_1="R_Cyt1" variable_2="R_Cyt1"/>
      <map_variables variable_1="R_Cyt2" variable_2="R_Cyt2"/>
   </connection>
   <connection>
      <map_components component_1="R_values" component_2="C_ER"/>
      <map_variables variable_1="C_ER" variable_2="C_ER"/>
      <map_variables variable_1="R_IP_3" variable_2="R_IP_3"/>
      <map_variables variable_1="R_ER" variable_2="R_ER"/>
      <map_variables variable_1="R_Cyt1" variable_2="R_Cyt1"/>
   </connection>
   <connection>
      <map_components component_1="G_GTP" component_2="Constants"/>
      <map_variables variable_1="k_0" variable_2="k_0"/>
      <map_variables variable_1="k_1" variable_2="k_1"/>
      <map_variables variable_1="k_2" variable_2="k_2"/>
      <map_variables variable_1="k_3" variable_2="k_3"/>
   </connection>
   <connection>
      <map_components component_1="APLC" component_2="Constants"/>
      <map_variables variable_1="C_PLC_T" variable_2="C_PLC_T"/>
      <map_variables variable_1="k_4" variable_2="k_4"/>
      <map_variables variable_1="k_5" variable_2="k_5"/>
   </connection>
   <connection>
      <map_components component_1="APLC" component_2="IP_3"/>
      <map_variables variable_1="APLC" variable_2="APLC"/>
   </connection>
   <connection>
      <map_components component_1="IP_3" component_2="Constants"/>
      <map_variables variable_1="k_6" variable_2="k_6"/>
      <map_variables variable_1="k_7" variable_2="k_7"/>
   </connection>
   <connection>
      <map_components component_1="IP_3" component_2="DG"/>
      <map_variables variable_1="IP_3" variable_2="IP_3"/>
   </connection>
   <connection>
      <map_components component_1="C_cyt" component_2="Constants"/>
      <map_variables variable_1="beta" variable_2="beta"/>
      <map_variables variable_1="rho" variable_2="rho"/>
      <map_variables variable_1="k_8" variable_2="k_8"/>
      <map_variables variable_1="k_9" variable_2="k_9"/>
      <map_variables variable_1="k_10" variable_2="k_10"/>
      <map_variables variable_1="k_11" variable_2="k_11"/>
   </connection>
   <connection>
      <map_components component_1="C_ER" component_2="Constants"/>
      <map_variables variable_1="lambda" variable_2="lambda"/>
      <map_variables variable_1="k_8" variable_2="k_8"/>
      <map_variables variable_1="k_9" variable_2="k_9"/>
   </connection>
</model>