Location: Activation of the Liver Glycogen Phosphorylase by Ca2+ Oscillations: a Theoretical Study @ 9cc2d56056c8 / gall_2000_take2.cellml

Author:
dave <dave@bcs.co.nz>
Date:
2013-01-04 12:32:48+13:00
Desc:
Equation 8 from original paper was wrong in model; correctly define units for all variables.
Permanent Source URI:
https://models.physiomeproject.org/workspace/35d/rawfile/9cc2d56056c85c516e53c1d590fe34916a8b4a48/gall_2000_take2.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 8/16/2011 at 10:43:13 at AM using:

COR (0.9.31.1409)
Copyright 2002-2011 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="gall_2000" cmeta:id="gall_2000" 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#">
   <units name="minute">
      <unit units="second" multiplier="60"/>
   </units>
   <units name="per_minute">
      <unit units="minute" exponent="-1"/>
   </units>
   <units name="millimole">
      <unit units="mole" prefix="milli"/>
   </units>
   <units name="micromole">
      <unit units="mole" prefix="micro"/>
   </units>
   <units name="micromolar">
      <unit units="micromole"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <component name="environment">
      <variable name="time" units="minute" initial_value="0" public_interface="out"/>
   </component>
   <component name="Pha">
      <variable name="Pha" units="dimensionless" initial_value="0.098"/>
      <variable name="V_1" units="dimensionless" public_interface="in"/>
      <variable name="Z" units="dimensionless" public_interface="in"/>
      <variable name="K_1" units="dimensionless" public_interface="in"/>
      <variable name="V_M2" units="per_minute" initial_value="0.6"/>
      <variable name="alpha" units="dimensionless" initial_value="9"/>
      <variable name="Glc" units="micromolar" initial_value="10000"/>
      <variable name="K_2" units="dimensionless" initial_value="0.2"/>
      <variable name="K_a1" units="micromolar" initial_value="10000"/>
      <variable name="K_a2" units="micromolar" initial_value="10000"/>
      <variable name="time" units="minute" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Pha</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>V_1</ci>
                  <ci>Z</ci>
                  <apply>
                     <divide/>
                     <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>Pha</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <apply>
                              <times/>
                              <ci>K_1</ci>
                              <ci>Z</ci>
                           </apply>
                           <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <ci>Pha</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>V_M2</ci>
                     <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                           <times/>
                           <ci>alpha</ci>
                           <ci>Glc</ci>
                           <apply>
                              <plus/>
                              <ci>K_a1</ci>
                              <ci>Glc</ci>
                           </apply>
                        </apply>
                     </apply>
                     <ci>Pha</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <apply>
                        <divide/>
                        <ci>K_2</ci>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <ci>Glc</ci>
                              <ci>K_a2</ci>
                           </apply>
                        </apply>
                     </apply>
                     <ci>Pha</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="V_1">
      <variable name="V_1" units="dimensionless" public_interface="out"/>
      <variable name="V_M1" units="per_minute" initial_value="1.5"/>
      <variable name="gamma" units="dimensionless" initial_value="9"/>
      <variable name="Z" units="dimensionless" public_interface="in"/>
      <variable name="K_a5" units="micromolar" initial_value="0.5"/>
      <variable name="time" units="minute" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_1</ci>
            <apply>
               <times/>
               <ci>V_M1</ci>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <times/>
                     <ci>gamma</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <divide/>
                           <apply>
                              <power/>
                              <ci>Z</ci>
                              <cn cellml:units="dimensionless">4</cn>
                           </apply>
                           <apply>
                              <power/>
                              <ci>K_a5</ci>
                              <cn cellml:units="dimensionless">4</cn>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Z</ci>
                           <cn cellml:units="dimensionless">4</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="K_1">
      <variable name="K_1" units="dimensionless" public_interface="out"/>
      <variable name="K_11" units="dimensionless" initial_value="0.1"/>
      <variable name="Z" units="dimensionless" public_interface="in"/>
      <variable name="K_a6" units="micromolar" initial_value="0.5"/>
      <variable name="time" units="minute" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>K_1</ci>
            <apply>
               <divide/>
               <ci>K_11</ci>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <power/>
                        <ci>Z</ci>
                        <cn cellml:units="dimensionless">4</cn>
                     </apply>
                     <apply>
                        <power/>
                        <ci>K_a6</ci>
                        <cn cellml:units="dimensionless">4</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="Z">
      <variable name="Z" units="dimensionless" public_interface="out"/>
      <variable name="A" units="micromolar" initial_value="0.3"/>
      <variable name="B" units="micromolar" initial_value="0.25"/>
      <variable name="v" units="per_minute" initial_value="1"/>
      <variable name="time" units="minute" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Z</ci>
            <apply>
               <plus/>
               <ci>A</ci>
               <apply>
                  <times/>
                  <ci>B</ci>
                  <apply>
                     <sin/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <pi/>
                        <ci>v</ci>
                        <ci>time</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <connection>
      <map_components component_1="Pha" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="V_1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="K_1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Z" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="Pha" component_2="Z"/>
      <map_variables variable_1="Z" variable_2="Z"/>
   </connection>
   <connection>
      <map_components component_1="Pha" component_2="V_1"/>
      <map_variables variable_1="V_1" variable_2="V_1"/>
   </connection>
   <connection>
      <map_components component_1="K_1" component_2="Pha"/>
      <map_variables variable_1="K_1" variable_2="K_1"/>
   </connection>
   <connection>
      <map_components component_1="K_1" component_2="Z"/>
      <map_variables variable_1="Z" variable_2="Z"/>
   </connection>
   <connection>
      <map_components component_1="V_1" component_2="Z"/>
      <map_variables variable_1="Z" variable_2="Z"/>
   </connection>
</model>