Location: Metabolic Component Library @ e0e2c28e2dfd / examples / Koenig2012_Gamma.cellml

Author:
Tommy Yu <tommy.yu@auckland.ac.nz>
Date:
2015-02-12 19:28:32+13:00
Desc:
Merging to remove the extra head - Since the title is changed and in the interest of maintaining the current structure we keep it as is.
Permanent Source URI:
https://models.cellml.org/w/matthiaskoenig/MetabolicComponentLibrary/rawfile/e0e2c28e2dfdfe37ea2cae3390a6ae26c28a3393/examples/Koenig2012_Gamma.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 3/20/2013 at 2:00:40 at AM using:

COR (0.9.31.1409)
Copyright 2002-2013 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="Gamma" cmeta:id="Gamma" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.1#">
	<import xmlns="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="Koenig2012_HormoneGlucoseResponseDecreasing.cellml">
		<component component_ref="HormoneGlucoseResponseDecreasing" name="EpinephrineResponseDecreasing"/>
	</import>
	<import xmlns="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="Koenig2012_HormoneGlucoseResponseDecreasing.cellml">
		<component component_ref="HormoneGlucoseResponseDecreasing" name="GlucagonResponseDecreasing"/>
	</import>
	
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:nodeID="n1">
			<boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#">
				<rdf:Description rdf:nodeID="n2"/>
			</boundIntervals>
		</rdf:Description>
		<rdf:Description rdf:nodeID="n3">
			<endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">20
      </endingValue>
		</rdf:Description>
		<rdf:Description rdf:about="#Gamma">
			<simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#">
				<rdf:Description rdf:nodeID="n1"/>
			</simulation>
		</rdf:Description>
		<rdf:Description rdf:nodeID="n2">
			<rdf:rest>
				<rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/>
			</rdf:rest>
			<rdf:first>
				<rdf:Description rdf:nodeID="n3"/>
			</rdf:first>
		</rdf:Description>
	</rdf:RDF>   <units name="mM">
      <unit units="mole" prefix="milli"/>
      <unit units="liter" exponent="-1"/>
   </units>
   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="mM_per_second">
      <unit units="mM"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="pM">
      <unit units="mole" prefix="pico"/>
      <unit units="liter" exponent="-1"/>
   </units>
   <component name="Gamma">
      <variable name="glc" units="mM" public_interface="in" private_interface="out"/>
      <variable name="gamma" units="dimensionless" public_interface="out"/>
      <variable name="ins" units="pM" public_interface="out" private_interface="in"/>
      <variable name="glu" units="pM" public_interface="out" private_interface="in"/>
      <variable name="epi" units="pM" public_interface="out" private_interface="in"/>
      <variable name="K_val" units="dimensionless"/>
      <variable name="epi_f" units="dimensionless"/>
      <variable name="K_ins" units="pM"/>
      <variable name="K_glu" units="pM"/>
      <variable name="K_epi" units="pM"/>
      <variable name="ins_max" units="pM"/>
      <variable name="ins_base" units="pM"/>
      <variable name="glu_max" units="pM"/>
      <variable name="glu_base" units="pM"/>
      <variable name="epi_max" units="pM"/>
      <variable name="epi_base" units="pM"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>ins_max</ci>
            <cn cellml:units="pM">818.9</cn>
         </apply>
         <apply>
            <eq/>
            <ci>ins_base</ci>
            <cn cellml:units="pM">0</cn>
         </apply>
         <apply>
            <eq/>
            <ci>glu_max</ci>
            <cn cellml:units="pM">190</cn>
         </apply>
         <apply>
            <eq/>
            <ci>glu_base</ci>
            <cn cellml:units="pM">37.9</cn>
         </apply>
         <apply>
            <eq/>
            <ci>epi_max</ci>
            <cn cellml:units="pM">190</cn>
         </apply>
         <apply>
            <eq/>
            <ci>epi_base</ci>
            <cn cellml:units="pM">37.9</cn>
         </apply>
         <apply>
            <eq/>
            <ci>K_val</ci>
            <cn cellml:units="dimensionless">0.1</cn>
         </apply>
         <apply>
            <eq/>
            <ci>epi_f</ci>
            <cn cellml:units="dimensionless">0.8</cn>
         </apply>
         <apply>
            <eq/>
            <ci>K_ins</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>ins_max</ci>
                  <ci>ins_base</ci>
               </apply>
               <ci>K_val</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K_glu</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>ins_max</ci>
                  <ci>ins_base</ci>
               </apply>
               <ci>K_val</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>K_epi</ci>
            <apply>
               <times/>
               <apply>
                  <minus/>
                  <ci>ins_max</ci>
                  <ci>ins_base</ci>
               </apply>
               <ci>K_val</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>gamma</ci>
            <piecewise>
               <piece>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.5</cn>
                     <apply>
                        <plus/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <ci>ins</ci>
                              <apply>
                                 <plus/>
                                 <ci>ins</ci>
                                 <ci>K_ins</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <divide/>
                           <ci>glu</ci>
                           <apply>
                              <plus/>
                              <ci>glu</ci>
                              <ci>K_glu</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <geq/>
                     <apply>
                        <divide/>
                        <ci>glu</ci>
                        <apply>
                           <plus/>
                           <ci>glu</ci>
                           <ci>K_glu</ci>
                        </apply>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>epi_f</ci>
                           <ci>epi</ci>
                        </apply>
                        <apply>
                           <plus/>
                           <ci>epi</ci>
                           <ci>K_epi</ci>
                        </apply>
                     </apply>
                  </apply>
               </piece>
               <otherwise>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">0.5</cn>
                     <apply>
                        <plus/>
                        <apply>
                           <minus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <ci>ins</ci>
                              <apply>
                                 <plus/>
                                 <ci>ins</ci>
                                 <ci>K_ins</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>epi_f</ci>
                              <ci>epi</ci>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>epi</ci>
                              <ci>K_epi</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="InsulinGlucoseResponse">
      <variable name="ins" units="pM" public_interface="out" private_interface="in"/>
      <variable name="glc" units="mM" public_interface="in" private_interface="out"/>
      <variable name="ins_max" units="pM" private_interface="out"/>
      <variable name="ins_base" units="pM" private_interface="out"/>
      <variable name="ins_km" units="pM" private_interface="out"/>
      <variable name="ins_n" units="dimensionless" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>ins_max</ci>
            <cn cellml:units="pM">818.9</cn>
         </apply>
         <apply>
            <eq/>
            <ci>ins_base</ci>
            <cn cellml:units="pM">0</cn>
         </apply>
         <apply>
            <eq/>
            <ci>ins_km</ci>
            <cn cellml:units="pM">8.6</cn>
         </apply>
         <apply>
            <eq/>
            <ci>ins_n</ci>
            <cn cellml:units="dimensionless">4.2</cn>
         </apply>
      </math>
   </component>
   <component name="GlucagonGlucoseResponse">
      <variable name="glu" units="pM" public_interface="out" private_interface="in"/>
      <variable name="glc" units="mM" public_interface="in" private_interface="out"/>
      <variable name="glu_max" units="pM" private_interface="out"/>
      <variable name="glu_base" units="pM" private_interface="out"/>
      <variable name="glu_km" units="pM" private_interface="out"/>
      <variable name="glu_n" units="dimensionless" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>glu_max</ci>
            <cn cellml:units="pM">190</cn>
         </apply>
         <apply>
            <eq/>
            <ci>glu_base</ci>
            <cn cellml:units="pM">37.9</cn>
         </apply>
         <apply>
            <eq/>
            <ci>glu_km</ci>
            <cn cellml:units="pM">3.01</cn>
         </apply>
         <apply>
            <eq/>
            <ci>glu_n</ci>
            <cn cellml:units="dimensionless">6.4</cn>
         </apply>
      </math>
   </component>
   <component name="EpinephrineGlucoseResponse">
      <variable name="epi" units="pM" public_interface="out" private_interface="in"/>
      <variable name="glc" units="mM" public_interface="in" private_interface="out"/>
      <variable name="epi_max" units="pM" private_interface="out"/>
      <variable name="epi_base" units="pM" private_interface="out"/>
      <variable name="epi_km" units="pM" private_interface="out"/>
      <variable name="epi_n" units="dimensionless" private_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>epi_max</ci>
            <cn cellml:units="pM">6090</cn>
         </apply>
         <apply>
            <eq/>
            <ci>epi_base</ci>
            <cn cellml:units="pM">100</cn>
         </apply>
         <apply>
            <eq/>
            <ci>epi_km</ci>
            <cn cellml:units="pM">3.1</cn>
         </apply>
         <apply>
            <eq/>
            <ci>epi_n</ci>
            <cn cellml:units="dimensionless">8.4</cn>
         </apply>
      </math>
   </component>
   <component name="HormoneGlucoseResponseIncreasing">
      <variable name="h" units="pM" public_interface="out"/>
      <variable name="glc" units="mM" public_interface="in"/>
      <variable name="h_max" units="pM" public_interface="in"/>
      <variable name="h_base" units="pM" public_interface="in"/>
      <variable name="h_km" units="pM" public_interface="in"/>
      <variable name="h_n" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>h</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <ci>h_max</ci>
                     <ci>h_base</ci>
                  </apply>
                  <apply>
                     <power/>
                     <ci>glc</ci>
                     <ci>h_n</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>glc</ci>
                     <ci>h_n</ci>
                  </apply>
                  <apply>
                     <power/>
                     <ci>h_km</ci>
                     <ci>h_n</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="InsulinGlucoseResponse">
         <component_ref component="HormoneGlucoseResponseIncreasing"/>
      </component_ref>
      <component_ref component="GlucagonGlucoseResponse">
         <component_ref component="GlucagonResponseDecreasing"/>
      </component_ref>
      <component_ref component="EpinephrineGlucoseResponse">
         <component_ref component="EpinephrineResponseDecreasing"/>
      </component_ref>
      <component_ref component="Gamma">
         <component_ref component="InsulinGlucoseResponse"/>
         <component_ref component="GlucagonGlucoseResponse"/>
         <component_ref component="EpinephrineGlucoseResponse"/>
      </component_ref>
   </group>
   <connection>
      <map_components component_1="Gamma" component_2="InsulinGlucoseResponse"/>
      <map_variables variable_1="glc" variable_2="glc"/>
      <map_variables variable_1="ins" variable_2="ins"/>
   </connection>
   <connection>
      <map_components component_1="Gamma" component_2="GlucagonGlucoseResponse"/>
      <map_variables variable_1="glc" variable_2="glc"/>
      <map_variables variable_1="glu" variable_2="glu"/>
   </connection>
   <connection>
      <map_components component_1="Gamma" component_2="EpinephrineGlucoseResponse"/>
      <map_variables variable_1="glc" variable_2="glc"/>
      <map_variables variable_1="epi" variable_2="epi"/>
   </connection>
   <connection>
      <map_components component_1="GlucagonGlucoseResponse" component_2="GlucagonResponseDecreasing"/>
      <map_variables variable_1="glu" variable_2="h"/>
      <map_variables variable_1="glc" variable_2="glc"/>
      <map_variables variable_1="glu_max" variable_2="h_max"/>
      <map_variables variable_1="glu_base" variable_2="h_base"/>
      <map_variables variable_1="glu_km" variable_2="h_km"/>
      <map_variables variable_1="glu_n" variable_2="h_n"/>
   </connection>
   <connection>
      <map_components component_1="EpinephrineGlucoseResponse" component_2="EpinephrineResponseDecreasing"/>
      <map_variables variable_1="epi" variable_2="h"/>
      <map_variables variable_1="glc" variable_2="glc"/>
      <map_variables variable_1="epi_max" variable_2="h_max"/>
      <map_variables variable_1="epi_base" variable_2="h_base"/>
      <map_variables variable_1="epi_km" variable_2="h_km"/>
      <map_variables variable_1="epi_n" variable_2="h_n"/>
   </connection>
   <connection>
      <map_components component_1="InsulinGlucoseResponse" component_2="HormoneGlucoseResponseIncreasing"/>
      <map_variables variable_1="ins" variable_2="h"/>
      <map_variables variable_1="glc" variable_2="glc"/>
      <map_variables variable_1="ins_max" variable_2="h_max"/>
      <map_variables variable_1="ins_base" variable_2="h_base"/>
      <map_variables variable_1="ins_km" variable_2="h_km"/>
      <map_variables variable_1="ins_n" variable_2="h_n"/>
   </connection>
</model>