Location: Lenbury, Pornsawad, 2005 @ b31eb49ec36d / lenbury_pornsawad_2005.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2009-11-25 13:48:56+13:00
Desc:
Added images in ai and svg fomrat, removed non pub med reference
Permanent Source URI:
http://models.cellml.org/workspace/lenbury_pornsawad_2005/rawfile/b31eb49ec36d0404b6d8a7f48de6e57941046ca4/lenbury_pornsawad_2005.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  lenbury_model_2005.xml

CREATED :  12th July 2007

LAST MODIFIED : 12th July 2007

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.1 Specification.

DESCRIPTION :  This file contains a CellML description of Lenbury and Pornsawad's 2005 delay-differential equation model of the feedback-controlled hypothalamus-pituitary-adrenal axis in humans. 

CHANGES:  
  
--><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="lenbury_model_2005" name="lenbury_pornsawad_model_2005_version01">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>A delay-differential equation model of the feedback-controlled hypothalamus-pituitary-adrenal axis in humans</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
The hypothalamus-pituitary-adrenal (HPA) axis is an essential stress-responsive component which initiates and co-ordinates physiological responses to both internal and external stresses, including disease and sleep deprivation.  These stimuli are initially detected by the central nervous system, and are fed into the HPA axis via the hypothalamus.  The hypothalamus responds by synthesising and secreting corticotropin-releasing hormone (CRH), which acts directly on the anterior pituitary to stimulate the production and release of corticotropin (ACTH).  ACTH binds to a specific cell surface receptor expressed by the cells of the adrenal cortex and cortisol is released as a result of this ACTH-binding.  Finally, cortisol can have negative feedback effects on various components of the HPA axis, including CRH release from the hypothalamus and ACTH secretion from the anterior pituitary.  All of these steps in the hormone cascade system are summarised schematically in the diagram below.  
</para>

<para>
In the paper described here, Lenbury and Pornsawad propose mathematical model incorporating time delays and negative feedback loops to describe the kinetics of the hormone cascade in the HPA axis.  They build further on the earlier HPA axis model of Lenbury and Pacheenburawana (published in 1991, and also described in CellML in the model repository), presenting a physiologically realistic, but relatively simple, mathematical representation of the hormone cascade.
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="lenbury_2005.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the hypothalamus-pituitary-adrenal (HPA) axis.  Stimulatory and inhibitory paths are indicated by the arrows and + or - signs respectively.  CRH represents corticotropin-releasing hormone and ACTH represents corticotropin.</caption>
</informalfigure>

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

<para>
A delay-differential equation model of the feedback-controlled hypothalamus-pituitary-adrenal axis in humans, Yongwimon Lenbury and Pornsarp Pornsawad, 2005, 
            <emphasis>Mathematical Medicine and Biology</emphasis>, 22, 15-33.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=15716298&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 15716298</ulink>
</para>

<para>
It should be noted that in its current form, the CellML description of the this model is unable to perfectly capture the simulation results of the published model, this is due to the time delays which are difficult to describe in the CellML code.
</para>

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




  <component name="environment">
    <variable units="dimensionless" public_interface="out" name="time"/>
  </component>
  
  <component cmeta:id="x" name="x">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="x">
        <dc:title>x</dc:title>
        <dcterms:alternative>CRH</dcterms:alternative>
        <dcterms:alternative>corticotropin-releasing hormone</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="dimensionless" public_interface="out" name="x"/>
    
    <variable units="dimensionless" name="y_delay2"/>
    
    <variable units="dimensionless" public_interface="in" name="y"/>
    <variable units="dimensionless" public_interface="in" name="beta"/>
    <variable units="dimensionless" public_interface="in" name="tau2"/>
    <variable units="dimensionless" public_interface="in" name="delta1"/>
    <variable units="dimensionless" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> x </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> delta1 </ci>
            <apply>
              <exp/>
              <apply>
                <times/>
                <ci> beta </ci>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <times/>
                    <apply>
                      <power/>
                      <ci> y_delay2 </ci>
                      <cn cellml:units="dimensionless"> 2.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> delta1 </ci>
            <ci> x </ci>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> y_delay2 </ci> 
        <apply>
          <times/>
          <ci> y </ci>
          <apply>
            <minus/>
            <ci> time </ci>
            <ci> tau2 </ci>
          </apply>
        </apply>
      </apply>      
    </math>
  </component>

  <component cmeta:id="y" name="y">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="y">
        <dc:title>y</dc:title>
        <dcterms:alternative>ACTH</dcterms:alternative>
        <dcterms:alternative>corticotropin</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="dimensionless" public_interface="out" name="y"/>

    <variable units="dimensionless" name="x_delay"/>
    <variable units="dimensionless" name="z_delay"/>
    
    <variable units="dimensionless" public_interface="in" name="x"/>
    <variable units="dimensionless" public_interface="in" name="z"/>
    <variable units="dimensionless" public_interface="in" name="beta"/>
    <variable units="dimensionless" public_interface="in" name="tau1"/>
    <variable units="dimensionless" public_interface="in" name="tau2"/>
    <variable units="dimensionless" public_interface="in" name="delta2"/>
    <variable units="dimensionless" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> y </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> delta2 </ci>
            <ci> x_delay </ci>
            <apply>
              <exp/>
              <apply>
                <times/>
                <ci> beta </ci>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <apply>
                    <times/>
                    <apply>
                      <power/>
                      <ci> z_delay </ci>
                      <cn cellml:units="dimensionless"> 2.0 </cn>
                    </apply>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> delta2 </ci>
            <ci> y </ci>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> x_delay </ci> 
        <apply>
          <times/>
          <ci> x </ci>
          <apply>
            <minus/>
            <ci> time </ci>
            <ci> tau1 </ci>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> z_delay </ci> 
        <apply>
          <times/>
          <ci> z </ci>
          <apply>
            <minus/>
            <ci> time </ci>
            <ci> tau2 </ci>
          </apply>
        </apply>
      </apply>     
    </math>
  </component>

  <component cmeta:id="z" name="z">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="z">
        <dc:title>z</dc:title>
        <dcterms:alternative>cortisol</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="dimensionless" public_interface="out" name="z"/>
    
    <variable units="dimensionless" name="y_delay1"/>

    <variable units="dimensionless" public_interface="in" name="y"/>
    <variable units="dimensionless" public_interface="in" name="tau1"/>
    <variable units="dimensionless" public_interface="in" name="delta3"/>
    <variable units="dimensionless" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> z </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> delta3 </ci>
            <ci> y_delay1 </ci>
          </apply>
          <apply>
            <times/>
            <ci> delta3 </ci>
            <ci> z </ci>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> y_delay1 </ci> 
        <apply>
          <times/>
          <ci> y </ci>
          <apply>
            <minus/>
            <ci> time </ci>
            <ci> tau1 </ci>
          </apply>
        </apply>
      </apply> 
    </math>
  </component>

  <component name="model_parameters"> 
    <variable units="dimensionless" public_interface="out" name="delta1" initial_value="0.5"/>
    <variable units="dimensionless" public_interface="out" name="delta2" initial_value="0.38"/>
    <variable units="dimensionless" public_interface="out" name="delta3" initial_value="0.6"/>
    <variable units="dimensionless" public_interface="out" name="beta" initial_value="1.091"/>
    <variable units="dimensionless" public_interface="out" name="tau1" initial_value="0.5"/>
    <variable units="dimensionless" public_interface="out" name="tau2" initial_value="0.77"/>
  </component>
  
  <connection>
    <map_components component_2="environment" component_1="x"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="y"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="z"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="y" component_1="x"/>
    <map_variables variable_2="x" variable_1="x"/>
    <map_variables variable_2="y" variable_1="y"/>
  </connection> 

  <connection>
    <map_components component_2="model_parameters" component_1="x"/>
    <map_variables variable_2="delta1" variable_1="delta1"/>
    <map_variables variable_2="tau2" variable_1="tau2"/>
    <map_variables variable_2="beta" variable_1="beta"/>
  </connection>

  <connection>
    <map_components component_2="z" component_1="y"/>
    <map_variables variable_2="y" variable_1="y"/>
    <map_variables variable_2="z" variable_1="z"/>
  </connection> 

  <connection>
    <map_components component_2="model_parameters" component_1="y"/>
    <map_variables variable_2="delta2" variable_1="delta2"/>
    <map_variables variable_2="tau1" variable_1="tau1"/>
    <map_variables variable_2="beta" variable_1="beta"/>
    <map_variables variable_2="tau2" variable_1="tau2"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="z"/>
    <map_variables variable_2="delta3" variable_1="delta3"/>
    <map_variables variable_2="tau1" variable_1="tau1"/>
  </connection>

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#8dfca5d4-902a-4cdb-a7b0-5a08b7fc9111">
    <rdf:li>cortisol</rdf:li>
    <rdf:li>endocrine</rdf:li>
    <rdf:li>adrenal cortex</rdf:li>
    <rdf:li>HPA axis</rdf:li>
    <rdf:li>hypothalamus</rdf:li>
    <rdf:li>hpa axis</rdf:li>
    <rdf:li>anterior pituitary</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#76567ea3-72c2-4f0f-9011-de826e133640">
    <rdf:li rdf:resource="rdf:#59fb4ec7-2e8b-4755-8d42-23c8c9d03d15"/>
    <rdf:li rdf:resource="rdf:#dc4692f6-6c75-43d4-9e9a-c138ee34140a"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#18b4cd84-64ad-4ced-a441-ca22a7fb5ceb">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bc9f0708-49c1-4da2-98a5-45b72a037b3b">
    <dc:creator rdf:resource="rdf:#4ba55f1e-2498-4772-9355-7c91aac3deb5"/>
    <rdf:value>This is a CellML description of Lenbury and Pornsawad's 2005 delay-differential equation model of the feedback-controlled hypothalamus-pituitary-adrenal axis in humans.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dc4692f6-6c75-43d4-9e9a-c138ee34140a">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#c5ed21ca-3ed5-4c81-baaa-40ff8356642e"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#3a32f1da-5f5f-4c3b-8a88-f22113484d27"/>
    <dcterms:created rdf:resource="rdf:#2bc08a7a-85df-49d9-b7a4-d1966123ebe3"/>
    <dc:creator rdf:resource="rdf:#9f4df56a-c154-43a5-b66a-a50f3e8d88b9"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6a2bad9f-0983-4d48-a61c-e15799caf551">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2bc08a7a-85df-49d9-b7a4-d1966123ebe3">
    <dcterms:W3CDTF>2007-07-13T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cdf6316c-6c14-4d0a-b0c6-3986a78cfa42">
    <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:#08ebe658-70db-4037-99d5-5c8dc7964d9e">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#8dfca5d4-902a-4cdb-a7b0-5a08b7fc9111"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9f4df56a-c154-43a5-b66a-a50f3e8d88b9">
    <vCard:ORG rdf:resource="rdf:#18b4cd84-64ad-4ced-a441-ca22a7fb5ceb"/>
    <vCard:EMAIL rdf:resource="rdf:#cdf6316c-6c14-4d0a-b0c6-3986a78cfa42"/>
    <vCard:N rdf:resource="rdf:#9b9920c0-f5b4-4aa2-97a0-c6ccc2cccc2a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c5ed21ca-3ed5-4c81-baaa-40ff8356642e">
    <vCard:Given>Pornsarp</vCard:Given>
    <vCard:Family>Pornsawad</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3a32f1da-5f5f-4c3b-8a88-f22113484d27">
    <dc:creator rdf:resource="rdf:#6a2bad9f-0983-4d48-a61c-e15799caf551"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3e4d0902-594e-4072-9fb0-ecb6af59bacc">
    <dcterms:W3CDTF>2005-03</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4ba55f1e-2498-4772-9355-7c91aac3deb5">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f5e7fdec-d063-437c-b17f-77ec4360ec00">
    <vCard:Given>Yongwimon</vCard:Given>
    <vCard:Family>Lenbury</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#90808be9-f152-481f-a389-fbc5a88d6b4b">
    <dc:title>Mathematical Medicine and Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1feacdb7-a0ae-42f2-a637-1623f4888d9b">
    <dc:creator rdf:resource="rdf:#76567ea3-72c2-4f0f-9011-de826e133640"/>
    <dc:title>A delay-differential equation model of the feedback-controlled hypothalamus-pituitary-adrenal axis in humans</dc:title>
    <bqs:volume>22</bqs:volume>
    <bqs:first_page>15</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#90808be9-f152-481f-a389-fbc5a88d6b4b"/>
    <dcterms:issued rdf:resource="rdf:#3e4d0902-594e-4072-9fb0-ecb6af59bacc"/>
    <bqs:last_page>33</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="#lenbury_model_2005">
    <dc:title>
        Lenbury and Pornsawad's 2005 delay-differential equation model of the feedback-controlled hypothalamus-pituitary-adrenal axis in humans.
      </dc:title>
    <cmeta:comment rdf:resource="rdf:#bc9f0708-49c1-4da2-98a5-45b72a037b3b"/>
    <bqs:reference rdf:resource="rdf:#69ff8280-9a4e-44ba-8943-821bb02b1cd0"/>
    <bqs:reference rdf:resource="rdf:#8f4de9c1-da1c-450a-b05c-848bc6aec772"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9b9920c0-f5b4-4aa2-97a0-c6ccc2cccc2a">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#69ff8280-9a4e-44ba-8943-821bb02b1cd0">
    <bqs:Pubmed_id>15716298</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#1feacdb7-a0ae-42f2-a637-1623f4888d9b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#59fb4ec7-2e8b-4755-8d42-23c8c9d03d15">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#f5e7fdec-d063-437c-b17f-77ec4360ec00"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8f4de9c1-da1c-450a-b05c-848bc6aec772">
    <dc:subject rdf:resource="rdf:#08ebe658-70db-4037-99d5-5c8dc7964d9e"/>
  </rdf:Description>
</rdf:RDF>
</model>