Location: Gonzalezheydrich, Steingard, Kohane, 1994 @ 129bd19c2178 / gonzalezheydrich_steingard_kohane_1994.cellml

Author:
Hanne <Hanne@hanne-nielsens-macbook.local>
Date:
2009-12-07 13:01:45+13:00
Desc:
Added iamges in ai and svg format
Permanent Source URI:
https://models.physiomeproject.org/workspace/gonzalezheydrich_steingard_kohane_1994/rawfile/129bd19c2178da6f5a77cb0252156f4b88aa7693/gonzalezheydrich_steingard_kohane_1994.cellml

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

CREATED :  28th June 2007

LAST MODIFIED : 28th June 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 Gonzalez-Heydrich et al.'s 1994 mathematical model of the HPA axis.

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="gonzalezheydrich_model_1994" name="gonzalezheydrich_steingard_kohane_model_1994_version01">

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

<article>
  <articleinfo>
  <title>A computer simulation of the hypothalamic-pituitary-adrenal axis</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 hypothalamic-pituitary-adrenal axis (HPA axis) regulates the secretion of cortisol through a complex network of interacting components.  Although our knowledge about each of these separate components is increasing, our understanding of how these components are combined into a network, and how this network functions, is still limited.  In an attempt to address this issue, mathematical models of the HPA axis are being developed, and represent an informative, valuable, exploratory tool.
</para>

<para>
The current model described here in CellML is based on the HPA axis model published by Gonzalez-Heydrich <emphasis>et al.</emphasis> (1994) (the complete original paper reference is cited below).  The mathematical model simulates components of the HPA axis involved in the continuous secretion and elimination of cortisol, adrenocorticotropin (ACTH), and corticotropin releasing hormone (CRH).  The relationships between these components were modelled based on current knowledge of their functioning, and model parameters were derived from experimental results in the published literature.  The presence of CRH and arginine-vasopressin (AVP), which are produced by the hypothalamus, stimulate the pituitary to secrete ACTH.  In turn, ACTH triggers the release of cortisol from the adrenal gland.  Cortisol has a direct negative feedback effect on the pituitary, and also on the hypothalamus, such that the secretion and synthesis of CRH, AVP and cortisol are inhibited.
</para>

<para>
A computer simulation of the hypothalamic-pituitary-adrenal axis, J. Gonzalez-Heydrich, R.J. Steingard and I Kohane, 1994, <emphasis>Eighteenth Annual Proceedings of the Symposium for Computer Applications in Medical Care</emphasis>, pp1010.  (The Annual Symposium on Computer Applications in Medical Care currently does not have a website, but the text of this article can be found online <ulink url="http://www.chip.org/projects/hpa/paper-hpa.html">here</ulink>).  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=7949852&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 7949852</ulink>
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="gonzalezheidrich_1994.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the self-regulatory system for corticol secretion.  CRH represents corticotropin-releasing hormone, and ACTH is adrenocorticotropic hormone.</caption>
</informalfigure>

<para>
Although this model is now more than 10 years old, it is still cited in the literature as one of the few mathematical models of the HPA axis which considers the dynamics of several hormones simultaneously, and doesn't focus solely on cortisol kinetics.
</para>

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



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

  <units name="minute">
    <unit units="second" multiplier="60.0"/>
  </units>

  <units name="flux">
    <unit units="mcg_ml"/>
    <unit units="minute" exponent="-1"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="minute" exponent="-1"/>
  </units>
  
  
  <component name="environment">
    <variable units="minute" public_interface="out" name="time"/>
  </component>
 
  
  <component name="cortisol"> 
    <variable units="mcg_ml" public_interface="out" name="cortisol" initial_value="0.0"/> 
    
    <variable units="first_order_rate_constant" name="k1" initial_value="5E5"/>
    <variable units="flux" name="k2" initial_value="0.01"/>
    <variable units="first_order_rate_constant" name="k3" initial_value="0.01"/>
    
    <variable units="mcg_ml" public_interface="in" name="ACTH"/>
    <variable units="minute" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> cortisol </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/> 
            <apply>
              <times/>
              <ci> k1 </ci>
              <ci> ACTH </ci>
            </apply>
            <ci> k2 </ci>
          </apply>
          <apply>
            <times/>
            <ci> k3 </ci>
            <ci> cortisol </ci>
          </apply>
        </apply>  
      </apply> 
    </math> 
  </component>
  
  <component cmeta:id="ACTH" name="ACTH">
     
    <variable units="mcg_ml" public_interface="out" name="ACTH" initial_value="0.0"/> 
    
    <variable units="first_order_rate_constant" name="k4" initial_value="10.0"/>
    <variable units="flux" name="k5" initial_value="4E-3"/>
    <variable units="first_order_rate_constant" name="k6" initial_value="0.035"/>
    <variable units="mcg_ml" name="Kd" initial_value="0.004"/>
    <variable units="dimensionless" name="Imax" initial_value="0.99"/>
    
    <variable units="mcg_ml" public_interface="in" name="CRH"/>
    <variable units="mcg_ml" public_interface="in" name="cortisol"/>
    <variable units="minute" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> ACTH </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/> 
            <apply>
              <times/>
              <ci> k4 </ci>
              <ci> CRH </ci>
            </apply>
            <ci> k5 </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> k6 </ci>
              <ci> ACTH </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <times/> 
                  <ci> k4 </ci>
                  <ci> CRH </ci>
                </apply>
                <ci> k5 </ci>
              </apply>
              <apply>
                <divide/>
                <apply>
                  <times/> 
                  <ci> Imax </ci>
                  <ci> cortisol </ci>
                </apply>
                <apply>
                  <plus/> 
                  <ci> Kd </ci>
                  <ci> cortisol </ci>
                </apply>
              </apply>  
            </apply> 
          </apply>
        </apply>
      </apply> 
    </math> 
  </component>

  <component cmeta:id="CRH" name="CRH">
     
    <variable units="mcg_ml" public_interface="out" name="CRH" initial_value="50.0"/> 
    
    <variable units="flux" name="k7" initial_value="1E-6"/>
    <variable units="first_order_rate_constant" name="k8" initial_value="0.01"/>
    <variable units="flux" name="pulse" initial_value="50.0"/>
   
    <variable units="minute" public_interface="in" name="time"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> CRH </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/> 
            <ci> pulse </ci>
            <ci> k7 </ci>
          </apply>
          <apply>
            <times/>
            <ci> k8 </ci>
            <ci> CRH </ci>
          </apply>
        </apply>
      </apply> 
    </math> 
  </component>


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

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

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

  <connection>
    <map_components component_2="ACTH" component_1="cortisol"/>
    <map_variables variable_2="cortisol" variable_1="cortisol"/>
    <map_variables variable_2="ACTH" variable_1="ACTH"/>
  </connection>

  <connection>
    <map_components component_2="ACTH" component_1="CRH"/>
    <map_variables variable_2="CRH" variable_1="CRH"/>
  </connection>

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#6aa7b954-3a56-40d5-bacd-0bb38097d15c">
    <rdf:li>endocrine</rdf:li>
    <rdf:li>hpa axis</rdf:li>
    <rdf:li>HPA axis</rdf:li>
    <rdf:li>cortisol</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#f9322b6e-1a17-4ece-8ab9-1d4e14f61430">
    <rdf:li rdf:resource="rdf:#67157005-608c-41d5-aab4-91ed55dfb71a"/>
    <rdf:li rdf:resource="rdf:#ca485a46-242a-4c7b-8882-4ca51ddb40bf"/>
    <rdf:li rdf:resource="rdf:#27dd82cb-f828-48d0-ad84-30350a18e333"/>
  </rdf:Seq>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#6e662e83-1fe0-4265-ae31-67c043257324"/>
    <dcterms:created rdf:resource="rdf:#daa2eb00-6237-46c5-afe0-423132846a42"/>
    <dc:creator rdf:resource="rdf:#ef11c371-0207-4132-98b2-5f0bd5b82875"/>
    <cmeta:modification rdf:resource="rdf:#200e2fae-7578-4f7a-b686-0fd48de3a7f7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8ad2b215-9367-412d-bb82-3e0793aaa00e">
    <vCard:N rdf:resource="rdf:#31f2ec2e-7b59-4472-952a-ba80c7af078c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#068b2416-0341-4959-898c-6d844493199b">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#11e9c07d-1856-4a00-a724-25eb5ee0c80b">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#6aa7b954-3a56-40d5-bacd-0bb38097d15c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ef11c371-0207-4132-98b2-5f0bd5b82875">
    <vCard:ORG rdf:resource="rdf:#feb3c480-5fce-4c75-a3e4-3d45851817df"/>
    <vCard:EMAIL rdf:resource="rdf:#fb49b7e2-e5be-461c-a4d1-1d7f41b61912"/>
    <vCard:N rdf:resource="rdf:#e08e0ada-599e-4e3b-aae1-0bec3ec5748c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#176ac537-a60d-4f4c-a36d-ffa46e104b9e">
    <bqs:Pubmed_id>7949852</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#773c6b80-7324-444a-bb73-173e2dfe0e4e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8d824ed5-f8d1-4995-801a-1a20920f0cce">
    <vCard:Given>R</vCard:Given>
    <vCard:Family>Steingard</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9601adc2-b041-454e-aed0-8728b733f7f8">
    <vCard:Given>J</vCard:Given>
    <vCard:Family>Gonzalez-Heydrich</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f2a7b556-82ee-415f-8ba9-c534ddc21cbe">
    <vCard:Given>I</vCard:Given>
    <vCard:Family>Kohane</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#200e2fae-7578-4f7a-b686-0fd48de3a7f7">
    <dcterms:modified rdf:resource="rdf:#305d706c-30bd-45cb-b84e-99c78afbab60"/>
    <rdf:value>The model has been checked in COR too.  The units have been corrected and are consistent.  The model runs in PCEnv but I'm not sure if it recreates the original results because the paper doesn't include any!  Also, I think there may be a typo in the paper's equation to calculate ACTH - it looks a bit odd.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#8ad2b215-9367-412d-bb82-3e0793aaa00e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#27dd82cb-f828-48d0-ad84-30350a18e333">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#f2a7b556-82ee-415f-8ba9-c534ddc21cbe"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#31f2ec2e-7b59-4472-952a-ba80c7af078c">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6e662e83-1fe0-4265-ae31-67c043257324">
    <dc:creator rdf:resource="rdf:#068b2416-0341-4959-898c-6d844493199b"/>
    <rdf:value>The model has been checked in COR too.  The units have been corrected and are consistent.  The model runs in PCEnv but I'm not sure if it recreates the original results because the paper doesn't include any!  Also, I think there may be a typo in the paper's equation to calculate ACTH - it looks a bit odd.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7864a706-fd96-4f71-91a4-67cc59e0e72a">
    <dcterms:W3CDTF>1994-00-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e08e0ada-599e-4e3b-aae1-0bec3ec5748c">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ca485a46-242a-4c7b-8882-4ca51ddb40bf">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#8d824ed5-f8d1-4995-801a-1a20920f0cce"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#feb3c480-5fce-4c75-a3e4-3d45851817df">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#305d706c-30bd-45cb-b84e-99c78afbab60">
    <dcterms:W3CDTF>2007-09-02T16:04:16+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bb39b022-dd12-4927-8371-57d7e5e27782">
    <dc:title>Eighteenth Annual Proceedings of the Symposium for Computer Applications in Medical Care</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#daa2eb00-6237-46c5-afe0-423132846a42">
    <dcterms:W3CDTF>2007-06-29T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fb49b7e2-e5be-461c-a4d1-1d7f41b61912">
    <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:#25ad4c5b-b17a-494a-b31b-2ad9c72b655f">
    <dc:subject rdf:resource="rdf:#11e9c07d-1856-4a00-a724-25eb5ee0c80b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#773c6b80-7324-444a-bb73-173e2dfe0e4e">
    <dc:creator rdf:resource="rdf:#f9322b6e-1a17-4ece-8ab9-1d4e14f61430"/>
    <dc:title>A computer simulation of the hypothalamic-pituitary-adrenal axis</dc:title>
    <bqs:volume/>
    <bqs:first_page>1010</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#bb39b022-dd12-4927-8371-57d7e5e27782"/>
    <dcterms:issued rdf:resource="rdf:#7864a706-fd96-4f71-91a4-67cc59e0e72a"/>
    <bqs:last_page/>
  </rdf:Description>
  <rdf:Description rdf:about="#gonzalezheydrich_model_1994">
    <dc:title>
        Gonzalez-Heydrich et al.'s 1994 mathematical model of the HPA axis.
      </dc:title>
    <cmeta:bio_entity>cortisol</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#54258794-b4a4-4d71-929c-166791385a26"/>
    <bqs:reference rdf:resource="rdf:#176ac537-a60d-4f4c-a36d-ffa46e104b9e"/>
    <bqs:reference rdf:resource="rdf:#25ad4c5b-b17a-494a-b31b-2ad9c72b655f"/>
  </rdf:Description>
  <rdf:Description rdf:about="/CRH">
    <dcterms:alternative>corticotropin releasing hormone</dcterms:alternative>
    <dc:title>CRH</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6218adcc-ef4f-4937-8e92-2c5fa605c4f6">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="/ACTH">
    <dcterms:alternative>acetylcholine</dcterms:alternative>
    <dc:title>ACTH</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#67157005-608c-41d5-aab4-91ed55dfb71a">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#9601adc2-b041-454e-aed0-8728b733f7f8"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#54258794-b4a4-4d71-929c-166791385a26">
    <dc:creator rdf:resource="rdf:#6218adcc-ef4f-4937-8e92-2c5fa605c4f6"/>
    <rdf:value>This is the CellML description of Gonzalez-Heydrich et al.'s 1994 mathematical model of the HPA axis.</rdf:value>
  </rdf:Description>
</rdf:RDF>
</model>