Location: Li, Kuang, Mason, 2006 @ 7ed85b69dfcd / li_kuang_mason_2006.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2010-04-09 12:36:37+12:00
Desc:
Balanced units. Fixed broken RDF. Updated curation comment.
Permanent Source URI:
https://models.cellml.org/workspace/li_kuang_mason_2006/rawfile/7ed85b69dfcd9aa15a660ad92a4d1a206c321de0/li_kuang_mason_2006.cellml

<?xml version='1.0' encoding='utf-8'?>

<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="li_model_2006" name="li_model_2006">

<documentation xmlns="http://cellml.org/tmp-documentation">

<article>
  <articleinfo>
  <title>Modelling the glucose-insulin regulatory system and ultradian insulin secretory oscillations</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
  <title>Model Status</title>
				<para>
				This CellML model runs in both COR and OpenCell however it does not recreate the published results. Glucose and insulin concentrations do not oscillate over time. This is most likely to be due to the lack of time delays in the CellML model (which are currently impossible to represent in CellML, but are required for equation 1 in the paper). The units are consistent. Parameter values have been taken from table 1 and di=0.06 and Gin=0.54 - as in figure 4.
				</para>
			</section>
  
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
ABSTRACT: In the glucose-insulin regulatory system, ultradian insulin secretory oscillations are observed to have a period of 50-150 min. After pioneering work traced back to the 1960s, several mathematical models have been proposed during the last decade to model these ultradian oscillations as well as the metabolic system producing them. These currently existing models still lack some of the key physiological aspects of the glucose-insulin system. Applying the mass conservation law, we introduce two explicit time delays and propose a more robust alternative model for better understanding the glucose-insulin endocrine metabolic regulatory system and the ultradian insulin secretory oscillations for the cases of continuous enteral nutrition and constant glucose infusion. We compare the simulation profiles obtained from this two time delay model with those from the other existing models. As a result, we notice many unique features of this two delay model. Based on our intensive simulations, we suspect that one of the possibly many causes of ultradian insulin secretion oscillations is the time delay of the insulin secretion stimulated by the elevated glucose concentration.
</para>


<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="li_2006.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the glucose-insulin regulatory system model.  The dotted lines indicate that elevated glucose concentration stimulates insulin synthesis and secretion by the pancreatic beta-cells, and also insulin promotes glucose utilisation in muscle, fat and other tissues.</caption>
</informalfigure>

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

<para>
Modeling the glucose-insulin regulatory system and ultradian insulin secretory oscillations with two explicit time delays, Jiaxu Li, Yang Kuang, and Clinton C. Mason, 2006, Journal of Theoretical Biology, 242, 722-735. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=16712872&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 16712872</ulink>
</para>
</sect1>
</article>
</documentation> 
  
  <units name="minute">
    <unit multiplier="60.0" units="second"/>
  </units>

  <units name="U" base_units="yes"/>

  <units name="mU_min">
    <unit prefix="milli" units="U"/>
    <unit units="minute" exponent="-1"/>
  </units>
  
  <units name="per_l">
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="l_min">
    <unit units="litre"/>
    <unit units="minute" exponent="-1"/>
  </units>

  <units name="l_mU">
    <unit units="litre"/>
    <unit prefix="milli" units="U" exponent="-1"/>
  </units>
  
  <units name="mg_l_min">
    <unit prefix="milli" units="gram"/>
    <unit units="litre" exponent="-1"/>
    <unit units="minute" exponent="-1"/>
  </units>

  <units name="mg_l">
    <unit prefix="milli" units="gram"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="mU_l">
    <unit prefix="milli" units="U"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="mg_min">
    <unit prefix="milli" units="gram"/>
    <unit units="minute" exponent="-1"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="minute" exponent="-1"/>
  </units>
  
  
  <component name="environment">
    <variable name="time" public_interface="out" units="minute"/>
  </component>
 
  <component name="G" cmeta:id="G">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/">
      <rdf:Description rdf:about="#G">
        <dc:title>G</dc:title>
        <dcterms:alternative>glucose</dcterms:alternative>    
      </rdf:Description>
    </rdf:RDF> 
    <variable units="mg_l" public_interface="out" initial_value="1000.0" name="G"/>
    <variable units="litre" public_interface="out" name="Vg" initial_value="10.0"/>
    
	<variable units="mg_l_min" name="Gin" initial_value="0.54"/>
    <variable units="mg_min" name="f2_G"/>
    <variable units="per_l" name="f3_G"/>
    <variable units="mg_min" name="f4_I"/>
    <variable units="mg_min" name="f5_I"/>
    <variable units="mg_l" name="C2" initial_value="144.0"/>
    <variable units="mg_l" name="C3" initial_value="1000.0"/>
    <variable units="mU_l" name="C4" initial_value="80.0"/>
    <variable units="mU_l" name="C5" initial_value="26.0"/>
    <variable units="mg_min" name="U0" initial_value="40.0"/>
    <variable units="mg_min" name="Um" initial_value="940.0"/>
    <variable units="mg_min" name="Ub" initial_value="72.0"/>
    <variable units="dimensionless" name="beta" initial_value="1.77"/>
    <variable units="mg_min" name="Rg" initial_value="180.0"/>
    <variable units="l_mU" name="alpha" initial_value="0.29"/>
    <variable units="minute" name="ti" initial_value="100"/>
    <variable units="l_min" name="E" initial_value="0.2"/>
    <variable units="litre" name="Vp" initial_value="3.0"/>
    <variable units="litre" name="Vi" initial_value="11.0"/>

    <variable units="minute" public_interface="in" name="time"/>
    <variable units="mU_l" public_interface="in" name="I"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="G_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> G </ci>
        </apply>
        <apply>
          <plus/>
          <ci> Gin </ci>
          <apply>
            <times/>
            <cn cellml:units="per_l"> 1.0 </cn>
			<ci> f5_I </ci>
		  </apply>
          <apply>
            <minus/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <cn cellml:units="per_l"> 1.0 </cn>
				<ci> f2_G </ci>
			  </apply>
              <apply>
                <times/>
                <ci> f3_G </ci>
                <ci> f4_I </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    
      <apply id="f2_G_calculation">
        <eq/>
        <ci> f2_G </ci>
        <apply>
          <times/>
          <ci> Ub </ci>
          <apply>
            <minus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
				<cn cellml:units="litre"> 1.0 </cn>
				<apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> G </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> C2 </ci>
                    <ci> Vg </ci>
				  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="f3_G_calculation">
        <eq/>
        <ci> f3_G </ci>
        <apply>
          <divide/>
          <ci> G </ci>
          <apply>
            <times/>
            <ci> C3 </ci>
            <ci> Vg </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="f4_I_calculation">
        <eq/>
        <ci> f4_I </ci>
        <apply>
          <plus/>
          <ci> U0 </ci>
          <apply>
            <divide/>
            <apply>
              <minus/>
              <ci> Um </ci>
              <ci> U0 </ci>
            </apply>
            <apply>
              <plus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <exp/>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <ci> beta </ci>
				  </apply>
                  <apply>
                    <ln/>
                    <apply>
                      <times/>
					  <cn cellml:units="per_l"> 1.0 </cn>
					  <apply>
                      <divide/>
                      <ci> I </ci>
                      <apply>
                        <times/>
                        <ci> C4 </ci>
						<apply>
                          <plus/>
                          <apply>
                            <divide/>
                            <cn cellml:units="dimensionless"> 1.0 </cn>
							<ci> Vi </ci>
						  </apply>
                          <apply>
                            <divide/>
                            <cn cellml:units="dimensionless"> 1.0 </cn>
                            <apply>
                              <times/>
                              <ci> E </ci>
                              <ci> ti </ci>
                            </apply>
                          </apply>
                        </apply>
					  </apply>
                      </apply>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply id="f5_I_calculation">
        <eq/>
        <ci> f5_I </ci>
        <apply>
          <divide/>
          <ci> Rg </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <times/>
                <ci> alpha </ci>
                <apply>
                  <divide/>
                  <ci> I </ci>
                  <apply>
                    <minus/>
                    <apply>
                      <times/>
				      <cn cellml:units="per_l"> 1.0 </cn>
					  <ci> Vp </ci>
					</apply>
                    <apply>
                      <times/>
				      <cn cellml:units="l_mU"> 1.0 </cn>
					  <ci> C5 </ci>
					</apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="I" cmeta:id="I">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/">
      <rdf:Description rdf:about="#I">
        <dc:title>I</dc:title>
        <dcterms:alternative>insulin</dcterms:alternative>    
      </rdf:Description>
    </rdf:RDF> 
    <variable units="mU_l" public_interface="out" initial_value="9000.0" name="I"/>
    
    <variable units="mU_min" name="f1_G"/>
    <variable units="mU_min" name="Rm" initial_value="210.0"/>
    <variable units="mg_l" name="C1" initial_value="2000.0"/>
    <variable units="mg_l" name="a1" initial_value="300.0"/>
    <variable units="first_order_rate_constant" name="di" initial_value="0.06"/>
    
    <variable units="minute" public_interface="in" name="time"/>
    <variable units="litre" public_interface="in" name="Vg"/>
    <variable units="mg_l" public_interface="in" name="G"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="I_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> I </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <cn cellml:units="per_l"> 1.0 </cn>
			<ci> f1_G </ci>
		  </apply>
          <apply>
            <times/>
            <ci> di </ci>
            <ci> I </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="f1_G_calculation">
        <eq/>
        <ci> f1_G </ci>
        <apply>
          <divide/>
          <ci> Rm </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <exp/>
              <apply>
                <divide/>
                <apply>
                  <divide/>
                  <apply>
                    <minus/>
                    <ci> C1 </ci>
                    <ci> G </ci>
                  </apply>
				  <apply>
                    <times/>
                    <cn cellml:units="per_l"> 1.0 </cn>
					<ci> Vg </ci>
				  </apply>
                </apply>
                <ci> a1 </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>


  <connection>
    <map_components component_2="environment" component_1="G"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="I"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_1="G" component_2="I"/>
    <map_variables variable_1="G" variable_2="G"/>
    <map_variables variable_1="I" variable_2="I"/>
	<map_variables variable_2="Vg" variable_1="Vg"/>
  </connection>


<rdf:RDF xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
  <rdf:Bag rdf:about="rdf:#e45ef7f6-d80f-4835-a330-7e8142f493da">
    <rdf:li>pancreas</rdf:li>
    <rdf:li>insulin</rdf:li>
    <rdf:li>glucose</rdf:li>
    <rdf:li>beta cell</rdf:li>
    <rdf:li>ultradian oscillations</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#af427b9c-8615-4631-a35d-303f9fe8018a">
    <rdf:li rdf:resource="rdf:#6f346fb1-e485-4ff2-bee5-7abe791c6f7f"/>
    <rdf:li rdf:resource="rdf:#11553c1c-7d1f-4586-8793-8fcec770e289"/>
    <rdf:li rdf:resource="rdf:#872b95ea-e6b9-4251-ab51-169487e6089a"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#1ce675f0-f19c-440f-8541-19a0339711b4">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#e45ef7f6-d80f-4835-a330-7e8142f493da"/>
  </rdf:Description>
  <rdf:Description rdf:about="#li_model_2006">
    <dc:title>
        Li et al.'s 2006 mathematical model of the glucose-insulin regulatory system and ultradian insulin secretory oscillations with two explicit time delays
      </dc:title>
    <cmeta:bio_entity>pancreas</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#a1035b85-5d58-46d0-bffe-f3eeee8e9dc9"/>
    <bqs:reference rdf:resource="rdf:#70cbbb59-5a0f-4a57-9af7-4e7529e44ff1"/>
    <bqs:reference rdf:resource="rdf:#b1084e00-edec-493d-a4f0-b0b5f8d6426f"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Auckland Bioengineering Institute
      </dc:publisher>
    <dcterms:created rdf:resource="rdf:#c6a77f38-8133-401f-94c1-100ee875f204"/>
    <dc:creator rdf:resource="rdf:#8cc644c9-e254-41dd-9436-784e92a6d2a3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6f346fb1-e485-4ff2-bee5-7abe791c6f7f">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#d66b4de0-a227-4c38-8ec3-f40ad469f2bc"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8cc644c9-e254-41dd-9436-784e92a6d2a3">
    <vCard:ORG rdf:resource="rdf:#b0951346-b994-4900-a796-1d4504c0d2ba"/>
    <vCard:EMAIL rdf:resource="rdf:#72c2b4c4-5fa3-4fc8-a2de-75cd957ab204"/>
    <vCard:N rdf:resource="rdf:#40099f9c-c3f2-48fb-b4fe-2d467a87b81d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#872b95ea-e6b9-4251-ab51-169487e6089a">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#90407fee-53f9-44da-b1ec-273066e8dadf"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#11553c1c-7d1f-4586-8793-8fcec770e289">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#455a1ef7-1081-4424-b78c-24332f60d25c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4f672018-f4e9-4d27-a92e-270f6a8056ba">
    <dc:creator rdf:resource="rdf:#af427b9c-8615-4631-a35d-303f9fe8018a"/>
    <dc:title>
            Modeling the glucose-insulin regulatory system and ultradian insulin secretory oscillations with two explicit time delays
          </dc:title>
    <bqs:volume>242</bqs:volume>
    <bqs:first_page>722</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#6c0dabdc-b0b1-47b5-991b-8075ad61a11a"/>
    <dcterms:issued rdf:resource="rdf:#93774758-841e-4bd5-a48c-a87fdbe59927"/>
    <bqs:last_page>735</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#72c2b4c4-5fa3-4fc8-a2de-75cd957ab204">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b0951346-b994-4900-a796-1d4504c0d2ba">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ca2a58ef-e3e2-46a3-86b5-55bb4df564aa">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#90407fee-53f9-44da-b1ec-273066e8dadf">
    <vCard:Given>Clinton</vCard:Given>
    <vCard:Family>Mason</vCard:Family>
    <vCard:Other>C</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6c0dabdc-b0b1-47b5-991b-8075ad61a11a">
    <dc:title>Journal of Theoretical Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d66b4de0-a227-4c38-8ec3-f40ad469f2bc">
    <vCard:Given>Jiaxu</vCard:Given>
    <vCard:Family>Li</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c6a77f38-8133-401f-94c1-100ee875f204">
    <dcterms:W3CDTF>2007-06-21</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#70cbbb59-5a0f-4a57-9af7-4e7529e44ff1">
    <bqs:Pubmed_id>16712872</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#4f672018-f4e9-4d27-a92e-270f6a8056ba"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#93774758-841e-4bd5-a48c-a87fdbe59927">
    <dcterms:W3CDTF>2006-10-07</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a1035b85-5d58-46d0-bffe-f3eeee8e9dc9">
    <dc:creator rdf:resource="rdf:#ca2a58ef-e3e2-46a3-86b5-55bb4df564aa"/>
    <rdf:value>
          This is the CellML description of Li et al.'s 2006 mathematical model of the glucose-insulin regulatory system and ultradian insulin secretory oscillations.  Note that the two time delays - tau1 and tau2 - included in the original published model have not been included in the current CellML description of the model.  At present, there is no way for the CellML to handle these delays.
        </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#455a1ef7-1081-4424-b78c-24332f60d25c">
    <vCard:Given>Yang</vCard:Given>
    <vCard:Family>Kuang</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#40099f9c-c3f2-48fb-b4fe-2d467a87b81d">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b1084e00-edec-493d-a4f0-b0b5f8d6426f">
    <dc:subject rdf:resource="rdf:#1ce675f0-f19c-440f-8541-19a0339711b4"/>
  </rdf:Description>
</rdf:RDF>
</model>