- Author:
- Catherine Lloyd <c.lloyd@auckland.ac.nz>
- Date:
- 2010-02-04 12:14:11+13:00
- Desc:
- Added initial conditions, corected units, updated documentation.
- Permanent Source URI:
- http://models.cellml.org/workspace/gupta_aslakson_gurbaxani_vernon_2007/rawfile/1014228fbdb7bdd602fad7ace80577f7c3a15b9d/gupta_aslakson_gurbaxani_vernon_2007_b.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="gupta_model_2007" name="gupta_model_2007">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Inclusion of the glucocorticoid receptor in a hypothalamic pituitary adrenal axis model reveals bistability</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 particular version of the model has a pulsatile stress and aims to recreate the results in figure 5 of the paper. This CellML model runs in both COR and OpenCell. The units have been checked and they are consistent. Unfortunately the CellML model is still not quite reproducing the published results. We will continue to curate the model and hope to get it fully working.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: BACKGROUND: The body's primary stress management system is the hypothalamic pituitary adrenal (HPA) axis. The HPA axis responds to physical and mental challenge to maintain homeostasis in part by controlling the body's cortisol level. Dysregulation of the HPA axis is implicated in numerous stress-related diseases. RESULTS: We developed a structured model of the HPA axis that includes the glucocorticoid receptor (GR). This model incorporates nonlinear kinetics of pituitary GR synthesis. The nonlinear effect arises from the fact that GR homodimerizes after cortisol activation and induces its own synthesis in the pituitary. This homodimerization makes possible two stable steady states (low and high) and one unstable state of cortisol production resulting in bistability of the HPA axis. In this model, low GR concentration represents the normal steady state, and high GR concentration represents a dysregulated steady state. A short stress in the normal steady state produces a small perturbation in the GR concentration that quickly returns to normal levels. Long, repeated stress produces persistent and high GR concentration that does not return to baseline forcing the HPA axis to an alternate steady state. One consequence of increased steady state GR is reduced steady state cortisol, which has been observed in some stress related disorders such as Chronic Fatigue Syndrome (CFS). CONCLUSION: Inclusion of pituitary GR expression resulted in a biologically plausible model of HPA axis bistability and hypocortisolism. High GR concentration enhanced cortisol negative feedback on the hypothalamus and forced the HPA axis into an alternative, low cortisol state. This model can be used to explore mechanisms underlying disorders of the HPA axis.
</para>
<para>
The original paper reference is cited below:
</para>
<para>
Inclusion of the glucocorticoid receptor in a hypothalamic pituitary adrenal axis model reveals bistability, Shakti Gupta, Eric Aslakson, Brian M. Gurbaxani and Suzanne D. Vernon, 2007,<emphasis>Theoretical Biology and Medical Modelling</emphasis>, volume 4, issue 8. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/17300722">PubMed ID: 17300722</ulink>
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="gupta_2007.png"/>
</imageobject>
</mediaobject>
<caption>The HPA axis model has three compartments, namely the hypothalamus, pituitary and adrenals. F is an external stress that triggers CRH (C) secretion by the hypothalamus. In turn, this signals to the pituitary to release ACTH (A) which stimulates the adrenal gland to release cortisol (O). Cortisol binds to the glucocorticoid receptor (R) in the pituitary and has a negative feedback effect on the release of CRH and ACTH. In addition, the cortisol-receptor complex has a positive feedback effect, promoting further receptor synthesis.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<units name="picomolar">
<unit units="mole" prefix="pico"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="nanomolar">
<unit units="mole" prefix="nano"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="nanomolar2">
<unit units="nanomolar" exponent="2"/>
</units>
<units name="nanomolar2_per_hour">
<unit units="nanomolar" exponent="2"/>
<unit units="hour" exponent="-1"/>
</units>
<units name="hour_per_nanomolar">
<unit units="hour" />
<unit units="nanomolar" exponent="-1"/>
</units>
<units name="hour">
<unit units="second" multiplier="3600.0"/>
</units>
<units name="first_order_rate_constant">
<unit units="hour" exponent="-1"/>
</units>
<component name="environment">
<variable units="hour" public_interface="out" name="time"/>
</component>
<component name="c" cmeta:id="c">
<variable units="picomolar" public_interface="out" initial_value="6.5" name="c"/>
<variable units="dimensionless" name="f"/>
<variable units="first_order_rate_constant" public_interface="in" name="kcd"/>
<variable units="dimensionless" public_interface="in" name="ki1"/>
<variable units="nanomolar" public_interface="in" name="o"/>
<variable units="hour" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci> c </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="nanomolar2_per_hour"> 1.0 </cn>
<apply>
<divide/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> f </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="nanomolar"> 1.0 </cn>
<apply>
<divide/>
<ci> o </ci>
<ci> ki1 </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kcd </ci>
<ci> c </ci>
</apply>
</apply>
</apply>
<apply id="f_calculation">
<eq/>
<ci> f </ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless"> 0.168 </cn>
<apply>
<or/>
<apply>
<and/>
<apply>
<gt/>
<ci> time </ci>
<cn cellml:units="hour"> 0.0 </cn>
</apply>
<apply>
<lt/>
<ci> time </ci>
<cn cellml:units="hour"> 2.0 </cn>
</apply>
</apply>
<apply>
<and/>
<apply>
<gt/>
<ci> time </ci>
<cn cellml:units="hour"> 8.0 </cn>
</apply>
<apply>
<lt/>
<ci> time </ci>
<cn cellml:units="hour"> 10.0 </cn>
</apply>
</apply>
<apply>
<and/>
<apply>
<gt/>
<ci> time </ci>
<cn cellml:units="hour"> 16.0 </cn>
</apply>
<apply>
<lt/>
<ci> time </ci>
<cn cellml:units="hour"> 18.0 </cn>
</apply>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless"> 0.0 </cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="a" cmeta:id="a">
<variable units="picomolar" public_interface="out" initial_value="1.2" name="a"/>
<variable units="first_order_rate_constant" public_interface="in" name="kad"/>
<variable units="dimensionless" public_interface="in" name="ki2"/>
<variable units="picomolar" public_interface="in" name="c"/>
<variable units="nanomolar" public_interface="in" name="o"/>
<variable units="dimensionless" public_interface="in" name="r"/>
<variable units="hour" 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/>
<apply>
<divide/>
<ci> c </ci>
<apply>
<times/>
<cn cellml:units="hour_per_nanomolar"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="nanomolar"> 1.0 </cn>
<apply>
<divide/>
<apply>
<times/>
<ci> o </ci>
<ci> r </ci>
</apply>
<ci> ki2 </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> kad </ci>
<ci> a </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="r" cmeta:id="r">
<variable units="dimensionless" public_interface="out" initial_value="0.08" name="r"/>
<variable units="first_order_rate_constant" public_interface="in" name="kcr"/>
<variable units="first_order_rate_constant" public_interface="in" name="krd"/>
<variable units="nanomolar2" public_interface="in" name="k"/>
<variable units="nanomolar" public_interface="in" name="o"/>
<variable units="hour" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci> r </ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<power/>
<apply>
<times/>
<ci> o </ci>
<ci> r </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<cn cellml:units="hour"> 1.0 </cn>
<apply>
<plus/>
<ci> k </ci>
<apply>
<power/>
<apply>
<times/>
<ci> o </ci>
<ci> r </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
<ci> kcr </ci>
</apply>
<apply>
<times/>
<ci> krd </ci>
<ci> r </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="o" cmeta:id="o">
<variable units="nanomolar" public_interface="out" initial_value="120.0" name="o"/>
<variable units="picomolar" public_interface="in" name="a"/>
<variable units="hour" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci> o </ci>
</apply>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant"> 1.0 </cn>
<apply>
<minus/>
<ci> a </ci>
<ci> o </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="reaction_constants" cmeta:id="reaction_constants">
<variable units="dimensionless" public_interface="out" name="ki1" initial_value="0.1"/>
<variable units="dimensionless" public_interface="out" name="ki2" initial_value="0.1"/>
<variable units="first_order_rate_constant" public_interface="out" name="kcd" initial_value="1.0"/>
<variable units="first_order_rate_constant" public_interface="out" name="kad" initial_value="10.0"/>
<variable units="nanomolar2" public_interface="out" name="k" initial_value="0.001"/>
<variable units="first_order_rate_constant" public_interface="out" name="kcr" initial_value="0.05"/>
<variable units="first_order_rate_constant" public_interface="out" name="krd" initial_value="0.9"/>
</component>
<connection>
<map_components component_2="environment" component_1="c"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<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="r"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="o"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="o" component_1="c"/>
<map_variables variable_2="o" variable_1="o"/>
</connection>
<connection>
<map_components component_2="o" component_1="a"/>
<map_variables variable_2="o" variable_1="o"/>
<map_variables variable_2="a" variable_1="a"/>
</connection>
<connection>
<map_components component_2="c" component_1="a"/>
<map_variables variable_2="c" variable_1="c"/>
</connection>
<connection>
<map_components component_2="r" component_1="a"/>
<map_variables variable_2="r" variable_1="r"/>
</connection>
<connection>
<map_components component_2="o" component_1="r"/>
<map_variables variable_2="o" variable_1="o"/>
</connection>
<connection>
<map_components component_2="reaction_constants" component_1="c"/>
<map_variables variable_2="kcd" variable_1="kcd"/>
<map_variables variable_2="ki1" variable_1="ki1"/>
</connection>
<connection>
<map_components component_2="reaction_constants" component_1="a"/>
<map_variables variable_2="kad" variable_1="kad"/>
<map_variables variable_2="ki2" variable_1="ki2"/>
</connection>
<connection>
<map_components component_2="reaction_constants" component_1="r"/>
<map_variables variable_2="kcr" variable_1="kcr"/>
<map_variables variable_2="krd" variable_1="krd"/>
<map_variables variable_2="k" variable_1="k"/>
</connection>
<rdf:RDF>
<rdf:Bag rdf:about="rdf:#d1123dfb-6743-4b24-868c-c2e3f1cbd7ec">
<rdf:li>glucocorticoid receptor</rdf:li>
<rdf:li>hypothalamus</rdf:li>
<rdf:li>cortisol</rdf:li>
<rdf:li>HPA axis</rdf:li>
</rdf:Bag>
<rdf:Seq rdf:about="rdf:#0c94cd55-b84e-48b5-b64a-ecbb9b6e664d">
<rdf:li rdf:resource="rdf:#2133e783-d498-4326-adb0-0c526b6648b2"/>
<rdf:li rdf:resource="rdf:#c13040af-4b78-440e-a21c-491479bd3ac8"/>
<rdf:li rdf:resource="rdf:#04332ed9-b7fc-45cb-8160-58f84655fd78"/>
<rdf:li rdf:resource="rdf:#ddbedd5e-182f-4b77-a29e-8ffcfcc8ae7f"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#e52791d8-eaa2-4f23-9fe1-2a1bbc1bfa38">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value rdf:resource="rdf:#d1123dfb-6743-4b24-868c-c2e3f1cbd7ec"/>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
<cmeta:comment rdf:resource="rdf:#d11dc723-cda3-4ba9-93b7-54ba78752831"/>
<dcterms:created rdf:resource="rdf:#b95c0e29-a762-4816-b88f-1d8187b8a5ec"/>
<dc:creator rdf:resource="rdf:#ebd3dd29-5544-4e80-837f-f95cb6ae87a9"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ebd3dd29-5544-4e80-837f-f95cb6ae87a9">
<vCard:ORG rdf:resource="rdf:#299f34f9-7669-4124-9115-bdafe7746bfd"/>
<vCard:EMAIL rdf:resource="rdf:#b7dae754-abdc-4387-9524-2e9e74c64f84"/>
<vCard:N rdf:resource="rdf:#caf02fff-30d4-48fd-88b3-1d7a8ca4c70a"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0ed4bc8c-3721-4636-8ddd-036e3a2ed7d5">
<vCard:Given>Shakti</vCard:Given>
<vCard:Family>Gupta</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d003903d-f3da-4082-b252-2c01aa1544e4">
<dc:creator rdf:resource="rdf:#b6c7d725-637d-4666-8243-8bdb263c3697"/>
<rdf:value>This is the CellML description of Gupta et al's 2007 hypothalamic pituitary adrenal axis model. The value of f - the stress input - has been defined such that it pulses to give the results in Figure 5 of the original publication.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="#gupta_model_2007">
<dc:title>
Gupta et al's 2007 hypothalamic pituitary adrenal axis model.
</dc:title>
<cmeta:bio_entity>hypothalamic pituitary adrenal axis</cmeta:bio_entity>
<cmeta:comment rdf:resource="rdf:#d003903d-f3da-4082-b252-2c01aa1544e4"/>
<bqs:reference rdf:resource="rdf:#577334c1-90d9-4b8b-88bf-179e96e4ca6f"/>
<bqs:reference rdf:resource="rdf:#afc253ef-8a4e-46c2-9009-975bf5448b0e"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#922de718-26b4-4e27-8e02-76fb12d635a1">
<vCard:Given>Brian</vCard:Given>
<vCard:Family>Gurbaxani</vCard:Family>
<vCard:Other>M</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="#c">
<dcterms:alternative>CRH</dcterms:alternative>
<dcterms:alternative>corticotropin releasing hormone</dcterms:alternative>
<dc:title>c</dc:title>
</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="rdf:#caf02fff-30d4-48fd-88b3-1d7a8ca4c70a">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b95c0e29-a762-4816-b88f-1d8187b8a5ec">
<dcterms:W3CDTF>2007-06-01T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="#o">
<dcterms:alternative>cortisol</dcterms:alternative>
<dc:title>o</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#04332ed9-b7fc-45cb-8160-58f84655fd78">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#922de718-26b4-4e27-8e02-76fb12d635a1"/>
</rdf:Description>
<rdf:Description rdf:about="#r">
<dcterms:alternative>GR</dcterms:alternative>
<dcterms:alternative>glucocorticoid receptor</dcterms:alternative>
<dc:title>r</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b6c7d725-637d-4666-8243-8bdb263c3697">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ddbedd5e-182f-4b77-a29e-8ffcfcc8ae7f">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#564f25aa-227c-4c75-8af1-20688ed8efa2"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b7dae754-abdc-4387-9524-2e9e74c64f84">
<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:#38c842cb-7134-4653-ac5f-07ee51ba2566">
<dcterms:W3CDTF>2007-02-14</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7a4a77d3-592f-4502-a9e7-e1b95bccd1b9">
<vCard:FN/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#577334c1-90d9-4b8b-88bf-179e96e4ca6f">
<dc:subject rdf:resource="rdf:#e52791d8-eaa2-4f23-9fe1-2a1bbc1bfa38"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#c13040af-4b78-440e-a21c-491479bd3ac8">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#6afd6f58-126b-4392-ae46-1ba85e6cc378"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#299f34f9-7669-4124-9115-bdafe7746bfd">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7de815e8-125b-491b-903a-830fdf4720ee">
<dc:creator rdf:resource="rdf:#0c94cd55-b84e-48b5-b64a-ecbb9b6e664d"/>
<dc:title>Inclusion of the glucocorticoid receptor in a hypothalamic pituitary adrenal axis model reveals bistability</dc:title>
<bqs:volume>4</bqs:volume>
<bqs:first_page/>
<bqs:Journal rdf:resource="rdf:#aafe6b4b-8eb9-4e29-94bd-857cc1c4ece6"/>
<dcterms:issued rdf:resource="rdf:#38c842cb-7134-4653-ac5f-07ee51ba2566"/>
<bqs:last_page/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2133e783-d498-4326-adb0-0c526b6648b2">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#0ed4bc8c-3721-4636-8ddd-036e3a2ed7d5"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6afd6f58-126b-4392-ae46-1ba85e6cc378">
<vCard:Given>Eric</vCard:Given>
<vCard:Family>Aslakson</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#564f25aa-227c-4c75-8af1-20688ed8efa2">
<vCard:Given>Suzanne</vCard:Given>
<vCard:Family>Vernon</vCard:Family>
<vCard:Other>D</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#aafe6b4b-8eb9-4e29-94bd-857cc1c4ece6">
<dc:title>Theoretical Biology and Medical Modelling</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d11dc723-cda3-4ba9-93b7-54ba78752831">
<dc:creator rdf:resource="rdf:#7a4a77d3-592f-4502-a9e7-e1b95bccd1b9"/>
<rdf:value/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#afc253ef-8a4e-46c2-9009-975bf5448b0e">
<bqs:Pubmed_id>17300722</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#7de815e8-125b-491b-903a-830fdf4720ee"/>
</rdf:Description>
</rdf:RDF>
</model>
