- Author:
- Hanne <Hanne@hanne-nielsens-macbook.local>
- Date:
- 2009-11-25 13:36:24+13:00
- Desc:
- Added images in ai svg and png format, removed non pub med references
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/jelic_cupic_kolaranic_2005/rawfile/ee2918cd3ca8af928506f167113463b1fecc03e5/jelic_cupic_kolaranic_2005.cellml
<?xml version='1.0' encoding='utf-8'?>
<!-- FILE : jelic_model_2005.xml
CREATED : 15th June 2007
LAST MODIFIED : 15th 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 Jelic et al.'s 2005 hypothalamic pituitary adrenal axis model.
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="jelic_model_2005" name="jelic_cupic_kolaranic_model_2005_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>A Hypothalamic-Pituitary-Adrenal System Model</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 (HPA) axis represents a neuroendocrine system which plays an essential role in maintaining body homeostasis in response to stress. These stresses can be physical (e.g. infection) or psychological (e.g. fear), and both activate the hypothalamus to secrete corticotropin releasing hormone (CRH). In turn, the presence of CRH in the circulation stimulates the pituitary to release adrenocorticotropic hormone (ACTH) into the blood, where it travels to the adrenal glands and induces the secretion of cortisol from the adrenal cortex. Cortisol has a negative feedback effect on the hypothalamus and pituitary such that the secretion of CRH and ACTH are inhibited.
</para>
<para>
The HPA system operates as a complex feedforward and feedback control system. There are three main pathways of activation for this system:
</para>
<itemizedlist>
<listitem>
<para>The internal biological clock, the hypothalamic supraiasmatic nucleus (SCN) with circadian oscillations in hormone levels;</para>
</listitem>
<listitem>
<para>Ultradian oscillations due to pulsatile patterns of HPA system hormone secretion and negative feedback loops; and</para>
</listitem>
<listitem>
<para>Stressful events (as described above) can also activate the HPA system and additional cortisol secretion occurs on top of the circadian and ultradian hormone oscillations.</para>
</listitem>
</itemizedlist>
<para>
Direct experimental assessment of complex biological systems is generally very difficult, and our understanding of these systems has been improved by the development of mathematical models. Several mathematical models describing the HPA system have been published (most recently by Gupta <emphasis>et al.</emphasis> (2007), and this model can be found in the CellML model repository along with other HPA models). The mathematical model described here in CellML (see figure below) is a four-dimensional non-linear differential equation model which can be used to analyse plasma cortisol levels in humans.
</para>
<para>
The complete original paper reference is cited below:
</para>
<para>
Mathematical modeling of the hypothalamic-pituitary-adrenal system activity., Smiljana Jelic, Zeljko Cupic, and Ljiljana Kolar-Anic, 2005, <emphasis>Mathematical Biosciences</emphasis>, 197, pp173-187. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&cmd=Retrieve&dopt=AbstractPlus&list_uids=16112688&query_hl=1&itool=pubmed_docsum">PubMed ID: 16112688</ulink>
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="jelic_2005.png"/>
</imageobject>
</mediaobject>
<caption>A schematic diagram of the hypothalamic-pituitary-adrenal system. Corticotrophin-releasing hormone (CRH) and arginin-vasopressin (AVP) stimulate adrenocorticotropin (ACTH) secretion from the pituitary, followed by cortisol secretion from the adrenal cortex. GR represents a glucocorticoid receptor and MR denotes a mineralocorticoid receptor. + represents a postive feedback loop while - represents a negative feedback loop.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<component name="environment">
<variable units="second" public_interface="out" name="time"/>
</component>
<component cmeta:id="a" name="a">
<variable units="dimensionless" public_interface="out" name="a" initial_value="7.68E-8"/>
<variable units="dimensionless" public_interface="in" name="alpha"/>
<variable units="dimensionless" public_interface="in" name="beta"/>
<variable units="dimensionless" public_interface="in" name="KP"/>
<variable units="dimensionless" public_interface="in" name="g"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci> a </ci>
</apply>
<apply>
<minus/>
<ci> KP </ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> alpha </ci>
<ci> beta </ci>
</apply>
<ci> a </ci>
</apply>
<apply>
<times/>
<ci> a </ci>
<apply>
<power/>
<ci> g </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="g" name="g">
<variable units="dimensionless" public_interface="out" name="g" initial_value="3.43E-8"/>
<variable units="dimensionless" public_interface="in" name="alpha"/>
<variable units="dimensionless" public_interface="in" name="gamma"/>
<variable units="dimensionless" public_interface="in" name="L"/>
<variable units="dimensionless" public_interface="in" name="a"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci> g </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> alpha </ci>
</apply>
<ci> a </ci>
</apply>
<apply>
<times/>
<ci> a </ci>
<apply>
<power/>
<ci> g </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci> L </ci>
<apply>
<times/>
<ci> gamma </ci>
<ci> g </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="reaction_constants" name="reaction_constants">
<variable units="dimensionless" public_interface="out" name="alpha"/>
<variable units="dimensionless" public_interface="out" name="beta"/>
<variable units="dimensionless" public_interface="out" name="gamma"/>
<variable units="dimensionless" public_interface="out" name="KP"/>
<variable units="dimensionless" public_interface="out" name="L"/>
<variable units="dimensionless" name="k2" initial_value="6.0E-4"/>
<variable units="dimensionless" name="k3" initial_value="0.0000048"/>
<variable units="dimensionless" name="k6" initial_value="0.000891"/>
<variable units="dimensionless" name="k7" initial_value="0.006831"/>
<variable units="dimensionless" name="k0" initial_value="8.7831E-11"/>
<variable units="dimensionless" name="k4" initial_value="2.1E12"/>
<variable units="dimensionless" name="km" initial_value="6.9001E-14"/>
<variable units="dimensionless" name="K"/>
<variable units="dimensionless" name="P"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci> alpha </ci>
<apply>
<divide/>
<ci> k3 </ci>
<ci> k2 </ci>
</apply>
</apply>
<apply>
<eq/>
<ci> beta </ci>
<apply>
<divide/>
<ci> k6 </ci>
<ci> k2 </ci>
</apply>
</apply>
<apply>
<eq/>
<ci> gamma </ci>
<apply>
<divide/>
<ci> k7 </ci>
<ci> k2 </ci>
</apply>
</apply>
<apply>
<eq/>
<ci> KP </ci>
<apply>
<times/>
<ci> K </ci>
<ci> P </ci>
</apply>
</apply>
<apply>
<eq/>
<ci> K </ci>
<apply>
<power/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<power/>
<ci> k0 </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> k4 </ci>
</apply>
<apply>
<power/>
<ci> k2 </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
<cn cellml:units="dimensionless"> 0.5 </cn>
</apply>
</apply>
<apply>
<eq/>
<ci> P </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 0.735 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.1563 </cn>
<apply>
<sin/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 3.141592653 </cn>
<ci> time </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 24.0 </cn>
<cn cellml:units="dimensionless"> 3600.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.3127 </cn>
<apply>
<abs/>
<apply>
<sin/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.141592653 </cn>
<ci> time </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 24.0 </cn>
<cn cellml:units="dimensionless"> 3600.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> L </ci>
<apply>
<power/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<power/>
<ci> km </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> k4 </ci>
</apply>
<apply>
<power/>
<ci> k2 </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
<cn cellml:units="dimensionless"> 0.5 </cn>
</apply>
</apply>
</math>
</component>
<connection>
<map_components component_2="environment" component_1="a"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<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="reaction_constants"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="g" component_1="a"/>
<map_variables variable_2="g" variable_1="g"/>
<map_variables variable_2="a" variable_1="a"/>
</connection>
<connection>
<map_components component_2="reaction_constants" component_1="a"/>
<map_variables variable_2="alpha" variable_1="alpha"/>
<map_variables variable_2="beta" variable_1="beta"/>
<map_variables variable_2="KP" variable_1="KP"/>
</connection>
<connection>
<map_components component_2="reaction_constants" component_1="g"/>
<map_variables variable_2="alpha" variable_1="alpha"/>
<map_variables variable_2="gamma" variable_1="gamma"/>
<map_variables variable_2="L" variable_1="L"/>
</connection>
<rdf:RDF>
<rdf:Bag rdf:about="rdf:#2e4d4ce4-6ae2-4061-bcf5-b2ef34531aff">
<rdf:li>endocrine</rdf:li>
<rdf:li>hpa axis</rdf:li>
<rdf:li>HPA axis</rdf:li>
<rdf:li>hypothalamus</rdf:li>
<rdf:li>cortisol</rdf:li>
</rdf:Bag>
<rdf:Seq rdf:about="rdf:#b1112dc7-e29b-46eb-a7df-f53c39f74743">
<rdf:li rdf:resource="rdf:#bd414160-67ab-46a4-8758-ca5babf67bc6"/>
<rdf:li rdf:resource="rdf:#8543587a-38f3-4d4c-a6ad-1ae3953632c2"/>
<rdf:li rdf:resource="rdf:#4d396b19-197d-4bb5-886f-e5448827352e"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#f58e6aa6-2b7b-4ac6-a009-77ada78baa68">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
<cmeta:comment rdf:resource="rdf:#67b53216-9e4e-4ed6-96a1-c124f83b1702"/>
<dcterms:created rdf:resource="rdf:#8e59f600-e339-4e25-b182-2329b0eed66a"/>
<dc:creator rdf:resource="rdf:#09b1e34e-4b56-4f51-8ac1-d37b6bef445a"/>
<cmeta:modification rdf:resource="rdf:#dd414d54-bb0c-41bc-a297-f9151606fb95"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#43c6d79f-c4d2-459a-8cd6-8f86280ab104">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#da1eb32e-2dc8-48f7-b748-8c9a5f240b35">
<dcterms:W3CDTF>2005-00-00 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ca677a9f-3125-4f4a-b992-99d9851754db">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#67b53216-9e4e-4ed6-96a1-c124f83b1702">
<dc:creator rdf:resource="rdf:#43c6d79f-c4d2-459a-8cd6-8f86280ab104"/>
<rdf:value>The model has now been checked in COR too. I've added a complicated pulse stimulus equation to the model - but I'm not sure it's working properly! Can PCEnv handle "sin"?
I've added initially values from the paper and have also made every varible (except time) dimensionless.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="/a">
<dcterms:alternative>ACTH</dcterms:alternative>
<dcterms:alternative>adrenocorticotropic hormone</dcterms:alternative>
<dc:title>a</dc:title>
</rdf:Description>
<rdf:Description rdf:about="/g">
<dcterms:alternative>cortisol</dcterms:alternative>
<dc:title>o</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8e59f600-e339-4e25-b182-2329b0eed66a">
<dcterms:W3CDTF>2007-06-15T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#633f05ae-1918-4673-842b-9016079fba7f">
<dc:creator rdf:resource="rdf:#b1112dc7-e29b-46eb-a7df-f53c39f74743"/>
<dc:title>Mathematical modeling of the hypothalamic-pituitary-adrenal system activity</dc:title>
<bqs:volume>197</bqs:volume>
<bqs:first_page>173</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#1b108538-9188-49c6-874a-f4198660c7ae"/>
<dcterms:issued rdf:resource="rdf:#da1eb32e-2dc8-48f7-b748-8c9a5f240b35"/>
<bqs:last_page>187</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#09b1e34e-4b56-4f51-8ac1-d37b6bef445a">
<vCard:ORG rdf:resource="rdf:#f58e6aa6-2b7b-4ac6-a009-77ada78baa68"/>
<vCard:EMAIL rdf:resource="rdf:#8f48c3b2-a7ae-468a-bdaa-25bbae42e8b1"/>
<vCard:N rdf:resource="rdf:#fd43087c-75d5-413e-b592-15972c11a718"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fd43087c-75d5-413e-b592-15972c11a718">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#bd414160-67ab-46a4-8758-ca5babf67bc6">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#e563758e-7d71-4220-a5d3-20403e3a7213"/>
</rdf:Description>
<rdf:Description rdf:about="#jelic_model_2005">
<dc:title>
Jelic et al's 2005 hypothalamic pituitary adrenal axis model.
</dc:title>
<cmeta:bio_entity>hypothalamic pituitary adrenal axis</cmeta:bio_entity>
<cmeta:comment rdf:resource="rdf:#09ad4983-41f9-45bf-a05d-70646bd83cd2"/>
<bqs:reference rdf:resource="rdf:#9848fe8c-27fc-4472-9b78-cce4ecfe58ee"/>
<bqs:reference rdf:resource="rdf:#d100751a-fd7a-4d33-94a5-0bceba99b70a"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a40f38cd-1dd5-4029-80a1-bb8ed596e7cc">
<vCard:Given>Ljiljana</vCard:Given>
<vCard:Family>Kolar-Anic</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c8d32cc3-51fb-418f-952d-bbc98c024d88">
<vCard:Given>Zeljko</vCard:Given>
<vCard:Family>Cupic</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8543587a-38f3-4d4c-a6ad-1ae3953632c2">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#c8d32cc3-51fb-418f-952d-bbc98c024d88"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d100751a-fd7a-4d33-94a5-0bceba99b70a">
<bqs:Pubmed_id>16112688</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#633f05ae-1918-4673-842b-9016079fba7f"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#09ad4983-41f9-45bf-a05d-70646bd83cd2">
<dc:creator rdf:resource="rdf:#761e2d67-2343-4106-9fbe-0f6427f13163"/>
<rdf:value>This is the CellML description of Jelic et al's 2005 hypothalamic pituitary adrenal axis model.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d95174cd-f8f8-445c-939f-784969bdc067">
<dcterms:W3CDTF>2007-09-03T07:48:26+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3fc83b18-f1a0-456a-a674-0203e2571610">
<vCard:N rdf:resource="rdf:#ca677a9f-3125-4f4a-b992-99d9851754db"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4d396b19-197d-4bb5-886f-e5448827352e">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#a40f38cd-1dd5-4029-80a1-bb8ed596e7cc"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#dd414d54-bb0c-41bc-a297-f9151606fb95">
<dcterms:modified rdf:resource="rdf:#d95174cd-f8f8-445c-939f-784969bdc067"/>
<rdf:value>The model has now been checked in COR too. I've added a complicated pulse stimulus equation to the model - but I'm not sure it's working properly! Can PCEnv handle "sin"?
I've added initially values from the paper and have also made every varible (except time) dimensionless.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#3fc83b18-f1a0-456a-a674-0203e2571610"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8f48c3b2-a7ae-468a-bdaa-25bbae42e8b1">
<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:#e563758e-7d71-4220-a5d3-20403e3a7213">
<vCard:Given>Smiljana</vCard:Given>
<vCard:Family>Jelic</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b963de07-ee98-4d4a-a54c-2ce4b25a212a">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value rdf:resource="rdf:#2e4d4ce4-6ae2-4061-bcf5-b2ef34531aff"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1b108538-9188-49c6-874a-f4198660c7ae">
<dc:title>Mathematical Biosciences</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#761e2d67-2343-4106-9fbe-0f6427f13163">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9848fe8c-27fc-4472-9b78-cce4ecfe58ee">
<dc:subject rdf:resource="rdf:#b963de07-ee98-4d4a-a54c-2ce4b25a212a"/>
</rdf:Description>
</rdf:RDF>
</model>
