- Author:
- Geoff Nunns <gnunns1@jhu.edu>
- Date:
- 2010-03-30 12:37:16+13:00
- Desc:
- Updated model status with abstract, fixed ulink errors and corrected typos in RDF.
- Permanent Source URI:
- https://models.cellml.org/workspace/macgregor_leng_2005/rawfile/b41ac05c216730316f3ab4b522bd40219a2a4694/macgregor_leng_2005_variant01.cellml
<?xml version='1.0' encoding='utf-8'?>
<!-- FILE : macgregor_model_2005.xml
CREATED : 19th June 2007
LAST MODIFIED : 19th 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 MacGregor and Leng's 2005 mathematical model of hypothalamic control of growth hormone secretion.
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="macgregor_model_2005" name="macgregor_leng_model_2005_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Modelling the Hypothalamic Control of Growth Hormone Secretion</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This particular CellML model describes the full (pituitary + hypothalamic) model of the published article. The units have been checked and are consistent, and the model runs in PCEnv to replicate the published results (as shown in figure 12 of the paper).
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
Abstract: Here, we construct a mathematical model of the hypothalamic systems that control the secretion of growth hormone (GH). The work extends a recent model of the pituitary GH system, adding representations of the hypothalamic GH-releasing hormone (GHRH) and somatostatin neurones, each modelled as a single synchronised unit. An unpatterned stochastic input drives the GHRH neurones generating pulses of GHRH release that trigger GH pulses. Delayed feedback from GH results in increased somatostatin release, which inhibits both GH secretion and GHRH release, producing an overall pattern of 3-h pulses of GH secretion that is very similar to the secretory profile observed in male rats. Rather than directly stimulating somatostatin release, GH feedback triggers a priming effect, increasing releasable stores of somatostatin. Varying this priming effect to reduce the effect of GH can reproduce the less pulsatile form of GH release observed in the female rat. The model behaviour is tested by comparison with experimental observations with a range of different experimental protocols involving GHRH injections and somatostatin and GH infusion.
</para>
<para>
The complete original paper reference is cited below:
</para>
<para>
Modelling the hypothalamic control of growth hormone secretion, D. J. MacGregor and G. Leng, 2005, <emphasis>Journal of Neuroendocrinology</emphasis>, volume 17, 788-803. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&cmd=Retrieve&dopt=AbstractPlus&list_uids=16280026&query_hl=1&itool=pubmed_docsum">PubMed ID: 16280026</ulink>
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="macgregor_2005.png"/>
</imageobject>
</mediaobject>
<caption>A schematic diagram of the different components and connections in the growth hormone secretory system. Somatostatin (SOS) inhibits the secretion of both growth hormone-releasing hormone (GHRH) and growth hormone (GH). After a delay, GH feeds back to prime the pool of releasable somatostatin, resulting in the bulk release of somatostatin, which in turn allows further GHRH release and and another burst of GH.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<units name="nanomolar">
<unit units="mole" prefix="nano"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="minute">
<unit units="second" multiplier="60.0"/>
</units>
<units name="first_order_rate_constant">
<unit units="minute" exponent="-1"/>
</units>
<units name="second_order_rate_constant">
<unit units="nanomolar" exponent="-1"/>
<unit units="minute" exponent="-1"/>
</units>
<units name="flux">
<unit units="nanomolar"/>
<unit units="minute" exponent="-1"/>
</units>
<units name="per_nanomolar">
<unit units="nanomolar" exponent="-1"/>
</units>
<component name="environment">
<variable units="minute" public_interface="out" name="time"/>
</component>
<component cmeta:id="r" name="r">
<variable units="nanomolar" public_interface="out" name="r" initial_value="0.0"/>
<variable units="dimensionless" name="ract"/>
<variable units="dimensionless" name="hr"/>
<variable units="flux" name="Ir"/>
<variable units="flux" name="Ir2"/>
<variable units="first_order_rate_constant" name="k6" initial_value="5.0"/>
<variable units="dimensionless" name="n1" initial_value="4.0"/>
<variable units="nanomolar" name="theta_1" initial_value="1.0"/>
<variable units="dimensionless" name="n3" initial_value="5.0"/>
<variable units="dimensionless" name="theta_3" initial_value="30.0"/>
<variable units="nanomolar" public_interface="in" name="s"/>
<variable units="dimensionless" public_interface="in" name="h_delta1"/>
<variable units="dimensionless" public_interface="in" name="j1"/>
<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> r </ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci> ract </ci>
<ci> hr </ci>
<ci> Ir </ci>
</apply>
<apply>
<times/>
<ci> k6 </ci>
<ci> r </ci>
</apply>
</apply>
<ci> Ir2 </ci>
</apply>
</apply>
<apply id="ract_calculation">
<eq/>
<ci> ract </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<apply>
<power/>
<ci> s </ci>
<ci> n1 </ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci> s </ci>
<ci> n1 </ci>
</apply>
<apply>
<power/>
<ci> theta_1 </ci>
<ci> n1 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="hr_calculation">
<eq/>
<ci> hr </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<apply>
<power/>
<ci> h_delta1 </ci>
<ci> n3 </ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci> h_delta1 </ci>
<ci> n3 </ci>
</apply>
<apply>
<power/>
<apply>
<times/>
<ci> j1 </ci>
<ci> theta_3 </ci>
</apply>
<ci> n3 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> Ir </ci>
<piecewise>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 0.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 90.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 91.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 92.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 93.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 113.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 114.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 115.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 116.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 136.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 137.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 138.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 139.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 159.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 160.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 161.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 162.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 252.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 253.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 254.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 255.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 275.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 276.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 277.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 278.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 298.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 299.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 300.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 301.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 321.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 322.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 323.0 </cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="flux"> 0.0 </cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci> Ir2 </ci>
<piecewise>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 0.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 5.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 1000.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 6.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 7.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 8.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 12.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 1000.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 13.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 14.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 15.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 21.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 1000.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 22.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 23.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 24.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 204.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 1000.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 205.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 206.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 207.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 217.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 1000.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 218.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 219.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 220.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 227.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 1000.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 228.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 229.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 230.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 310.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 1000.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 311.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 312.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 313.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 321.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 1000.0 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="minute"> 322.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 323.0 </cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="flux"> 0.0 </cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component cmeta:id="s" name="s">
<variable units="nanomolar" public_interface="out" name="s" initial_value="0.0"/>
<variable units="flux" name="Is2" initial_value="50.0"/>
<variable units="flux" public_interface="in" name="Is"/>
<variable units="first_order_rate_constant" public_interface="in" name="k7"/>
<variable units="dimensionless" public_interface="in" name="g"/>
<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> s </ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci> Is </ci>
<ci> g </ci>
</apply>
<apply>
<times/>
<ci> k7 </ci>
<ci> s </ci>
</apply>
</apply>
<ci> Is2 </ci>
</apply>
</apply>
</math>
</component>
<component cmeta:id="se" name="se">
<variable units="nanomolar" public_interface="out" name="se" initial_value="0.0"/>
<variable units="flux" name="Ise" initial_value="10.0"/>
<variable units="first_order_rate_constant" public_interface="in" name="k7"/>
<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> se </ci>
</apply>
<apply>
<minus/>
<ci> Ise </ci>
<apply>
<times/>
<ci> k7 </ci>
<ci> se </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="f" name="f">
<variable units="dimensionless" public_interface="out" name="f" initial_value="0.3"/>
<variable units="second_order_rate_constant" name="k1" initial_value="0.1"/>
<variable units="first_order_rate_constant" name="k2" initial_value="0.002"/>
<variable units="first_order_rate_constant" name="k3" initial_value="0.018"/>
<variable units="dimensionless" name="phi_b_s"/>
<variable units="dimensionless" name="sb" initial_value="0.029"/>
<variable units="dimensionless" name="delta_b" initial_value="0.3"/>
<variable units="nanomolar" public_interface="in" name="r"/>
<variable units="nanomolar" public_interface="in" name="c"/>
<variable units="nanomolar" public_interface="in" name="s"/>
<variable units="nanomolar" public_interface="in" name="se"/>
<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> f </ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci> k1 </ci>
<apply>
<plus/>
<ci> r </ci>
<ci> c </ci>
</apply>
<ci> f </ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci> k2 </ci>
<apply>
<times/>
<ci> k3 </ci>
<ci> phi_b_s </ci>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> f </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="phi_b_s_calculation">
<eq/>
<ci> phi_b_s </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<log/>
<apply>
<times/>
<cn cellml:units="per_nanomolar"> 1.0 </cn>
<apply>
<plus/>
<ci> s </ci>
<ci> se </ci>
</apply>
</apply>
</apply>
<ci> sb </ci>
</apply>
</apply>
<ci> delta_b </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="h" name="h">
<variable units="nanomolar" public_interface="out" name="h" initial_value="0.0"/>
<variable units="first_order_rate_constant" name="k4" initial_value="0.5"/>
<variable units="first_order_rate_constant" name="k5" initial_value="71.0"/>
<variable units="dimensionless" name="phi_r_s"/>
<variable units="dimensionless" name="sr" initial_value="-0.56"/>
<variable units="dimensionless" name="delta_r" initial_value="0.2"/>
<variable units="dimensionless" public_interface="in" name="f"/>
<variable units="first_order_rate_constant" public_interface="in" name="k8"/>
<variable units="dimensionless" public_interface="in" name="j1"/>
<variable units="nanomolar" public_interface="in" name="r"/>
<variable units="nanomolar" public_interface="in" name="c"/>
<variable units="nanomolar" public_interface="in" name="s"/>
<variable units="nanomolar" public_interface="in" name="se"/>
<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> h </ci>
</apply>
<apply>
<times/>
<ci> j1 </ci>
<apply>
<minus/>
<apply>
<times/>
<apply>
<plus/>
<ci> k4 </ci>
<apply>
<times/>
<ci> k5 </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> phi_r_s </ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci> r </ci>
<ci> c </ci>
</apply>
<ci> f </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k8 </ci>
<ci> h </ci>
</apply>
</apply>
</apply>
</apply>
<apply id="phi_r_s_calculation">
<eq/>
<ci> phi_r_s </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<log/>
<apply>
<times/>
<cn cellml:units="per_nanomolar"> 1.0 </cn>
<apply>
<plus/>
<ci> s </ci>
<ci> se </ci>
</apply>
</apply>
</apply>
<ci> sr </ci>
</apply>
</apply>
<ci> delta_r </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="g" name="g">
<variable units="dimensionless" public_interface="out" name="g" initial_value="2.0"/>
<variable units="first_order_rate_constant" name="g1" initial_value="1.0"/>
<variable units="dimensionless" name="gmax" initial_value="5.0"/>
<variable units="per_nanomolar" name="g2" initial_value="0.008"/>
<variable units="dimensionless" name="hact"/>
<variable units="dimensionless" name="n2" initial_value="2.0"/>
<variable units="dimensionless" name="theta_2" initial_value="30.0"/>
<variable units="dimensionless" public_interface="in" name="h_delta"/>
<variable units="flux" public_interface="in" name="Is"/>
<variable units="dimensionless" public_interface="in" name="j1"/>
<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> g </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> g1 </ci>
<ci> hact </ci>
<apply>
<divide/>
<apply>
<minus/>
<ci> gmax </ci>
<ci> g </ci>
</apply>
<ci> gmax </ci>
</apply>
</apply>
<apply>
<times/>
<ci> g2 </ci>
<ci> Is </ci>
<ci> g </ci>
</apply>
</apply>
</apply>
<apply id="hact_calculation">
<eq/>
<ci> hact </ci>
<apply>
<divide/>
<apply>
<power/>
<ci> h_delta </ci>
<ci> n2 </ci>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci> h_delta </ci>
<ci> n2 </ci>
</apply>
<apply>
<power/>
<apply>
<times/>
<ci> j1 </ci>
<ci> theta_2 </ci>
</apply>
<ci> n2 </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="hh" name="hh">
<variable units="nanomolar" public_interface="out" name="hh" initial_value="0.0"/>
<variable units="flux" name="Ih" initial_value="50.0"/>
<variable units="dimensionless" public_interface="in" name="j1"/>
<variable units="first_order_rate_constant" public_interface="in" name="k8"/>
<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> hh </ci>
</apply>
<apply>
<times/>
<ci> j1 </ci>
<apply>
<minus/>
<ci> Ih </ci>
<apply>
<times/>
<ci> k8 </ci>
<ci> hh </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="model_parameters">
<variable units="dimensionless" public_interface="out" name="h_delta"/>
<variable units="dimensionless" public_interface="out" name="h_delta1"/>
<variable units="dimensionless" public_interface="out" name="j1" initial_value="10"/>
<variable units="first_order_rate_constant" public_interface="out" name="k7" initial_value="5.0"/>
<variable units="flux" public_interface="out" name="Is"/>
<variable units="first_order_rate_constant" public_interface="out" name="k8" initial_value="0.07"/>
<variable units="nanomolar" public_interface="out" name="c" initial_value="0.01"/>
<variable units="per_nanomolar" name="delta" initial_value="60.0"/>
<variable units="per_nanomolar" name="delta1" initial_value="15.0"/>
<variable units="nanomolar" public_interface="in" name="h"/>
<variable units="nanomolar" public_interface="in" name="hh"/>
<variable units="minute" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci> h_delta </ci>
<apply>
<times/>
<apply>
<plus/>
<ci> h </ci>
<ci> hh </ci>
</apply>
<ci> delta </ci>
</apply>
</apply>
<apply>
<eq/>
<ci> h_delta1 </ci>
<apply>
<times/>
<apply>
<plus/>
<ci> h </ci>
<ci> hh </ci>
</apply>
<ci> delta1 </ci>
</apply>
</apply>
<apply>
<eq/>
<ci> Is </ci>
<piecewise>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<gt/>
<ci> time </ci>
<cn cellml:units="minute"> 0.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 90.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<gt/>
<ci> time </ci>
<cn cellml:units="minute"> 90.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 180.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 10.0 </cn>
<apply>
<and/>
<apply>
<gt/>
<ci> time </ci>
<cn cellml:units="minute"> 180.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 270.0 </cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="flux"> 0.0 </cn>
<apply>
<and/>
<apply>
<gt/>
<ci> time </ci>
<cn cellml:units="minute"> 270.0 </cn>
</apply>
<apply>
<leq/>
<ci> time </ci>
<cn cellml:units="minute"> 360.0 </cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="flux"> 0.0 </cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<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="s"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="se"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="f"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="h"/>
<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="hh"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="model_parameters"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="s" component_1="r"/>
<map_variables variable_2="s" variable_1="s"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="r"/>
<map_variables variable_2="h_delta1" variable_1="h_delta1"/>
<map_variables variable_2="j1" variable_1="j1"/>
</connection>
<connection>
<map_components component_2="g" component_1="s"/>
<map_variables variable_2="g" variable_1="g"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="s"/>
<map_variables variable_2="Is" variable_1="Is"/>
<map_variables variable_2="k7" variable_1="k7"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="se"/>
<map_variables variable_2="k7" variable_1="k7"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="f"/>
<map_variables variable_2="c" variable_1="c"/>
</connection>
<connection>
<map_components component_2="r" component_1="f"/>
<map_variables variable_2="r" variable_1="r"/>
</connection>
<connection>
<map_components component_2="s" component_1="f"/>
<map_variables variable_2="s" variable_1="s"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="h"/>
<map_variables variable_2="c" variable_1="c"/>
<map_variables variable_2="j1" variable_1="j1"/>
<map_variables variable_2="k8" variable_1="k8"/>
<map_variables variable_2="h" variable_1="h"/>
</connection>
<connection>
<map_components component_2="r" component_1="h"/>
<map_variables variable_2="r" variable_1="r"/>
</connection>
<connection>
<map_components component_2="s" component_1="h"/>
<map_variables variable_2="s" variable_1="s"/>
</connection>
<connection>
<map_components component_2="se" component_1="h"/>
<map_variables variable_2="se" variable_1="se"/>
</connection>
<connection>
<map_components component_2="se" component_1="f"/>
<map_variables variable_2="se" variable_1="se"/>
</connection>
<connection>
<map_components component_2="f" component_1="h"/>
<map_variables variable_2="f" variable_1="f"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="g"/>
<map_variables variable_2="Is" variable_1="Is"/>
<map_variables variable_2="j1" variable_1="j1"/>
<map_variables variable_2="h_delta" variable_1="h_delta"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="hh"/>
<map_variables variable_2="j1" variable_1="j1"/>
<map_variables variable_2="k8" variable_1="k8"/>
<map_variables variable_2="hh" variable_1="hh"/>
</connection>
<rdf:RDF>
<rdf:Bag rdf:about="rdf:#21108480-0177-4661-aa67-a32d8cc46324">
<rdf:li>growth hormone</rdf:li>
<rdf:li>hypothalamus</rdf:li>
<rdf:li>somatostatin</rdf:li>
<rdf:li>growth hormone releasing hormone</rdf:li>
</rdf:Bag>
<rdf:Seq rdf:about="rdf:#c9b84d74-3bc4-4c84-aa34-bf92fb3d6e2a">
<rdf:li rdf:resource="rdf:#defe8854-2305-4768-abe1-366b4c1c1e41"/>
<rdf:li rdf:resource="rdf:#5d83f7df-7c6e-43e0-84df-2cb20acf993b"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#a1307b00-a49b-4b59-bd31-ff8cff4c7e15">
<vCard:ORG rdf:resource="rdf:#4dc2f283-a7b9-413c-9822-c96c644a2874"/>
<vCard:EMAIL rdf:resource="rdf:#0074e35a-7320-417c-8ffd-da04aa3f74ca"/>
<vCard:N rdf:resource="rdf:#7973be78-1364-49f9-8098-5a9724200a00"/>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
<cmeta:comment rdf:resource="rdf:#790fdbd3-9a41-4af0-bea7-89f7be75237c"/>
<dcterms:created rdf:resource="rdf:#53397ee2-b1f5-417f-8d93-961225d7d1ad"/>
<dc:creator rdf:resource="rdf:#a1307b00-a49b-4b59-bd31-ff8cff4c7e15"/>
<cmeta:modification rdf:resource="rdf:#fb240565-8ff6-43f8-8a7a-1e56437f76fd"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#defe8854-2305-4768-abe1-366b4c1c1e41">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#39694593-b277-45a7-b193-2a5668a1a2cb"/>
</rdf:Description>
<rdf:Description rdf:about="#hh">
<dcterms:alternative>human growth hormone</dcterms:alternative>
<dc:title>hh</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4dc2f283-a7b9-413c-9822-c96c644a2874">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#be2f82ba-dcbc-4733-94ee-0ecd2fddbdd5">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="#se">
<dcterms:alternative>exogenous somatostatin</dcterms:alternative>
<dc:title>se</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5d83f7df-7c6e-43e0-84df-2cb20acf993b">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#f52889a3-8e66-4864-aef2-88c1f43b54e7"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fb240565-8ff6-43f8-8a7a-1e56437f76fd">
<dcterms:modified rdf:resource="rdf:#97a1988a-9346-4cb6-8b05-c569e944e8fc"/>
<rdf:value>I've checked the model in COR and have corrected all the units. I've also added some piecewise equations to stimulate the model (as set out in the figure legends of the paper).</rdf:value>
<cmeta:modifier rdf:resource="rdf:#738722ff-ed9d-472c-a085-f3433a6fb047"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5cfc0d4a-77e7-4e58-ac08-ca282aab57db">
<bqs:Pubmed_id>16280026</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#9ae82fcb-97af-4380-a8a3-0d71d815b744"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#73ab6072-f902-4e94-8cc9-ae0e3dd3f257">
<dcterms:W3CDTF>2005-00-00 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#53397ee2-b1f5-417f-8d93-961225d7d1ad">
<dcterms:W3CDTF>2007-05-28T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#97a1988a-9346-4cb6-8b05-c569e944e8fc">
<dcterms:W3CDTF>2007-09-04T15:52:40+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="#f">
<dcterms:alternative>density of free GHRH receptors</dcterms:alternative>
<dc:title>f</dc:title>
</rdf:Description>
<rdf:Description rdf:about="#g">
<dcterms:alternative>releasable store of somatostatin</dcterms:alternative>
<dc:title>g</dc:title>
</rdf:Description>
<rdf:Description rdf:about="#macgregor_model_2005">
<dc:title>
MacGregor and Leng's 2005 mathematical model of hypothalamic control of growth hormone secretion.
</dc:title>
<cmeta:bio_entity>hypothalamus</cmeta:bio_entity>
<cmeta:comment rdf:resource="rdf:#a8dfbb63-684f-4412-8aad-4afa2b5b47f7"/>
<bqs:reference rdf:resource="rdf:#5cfc0d4a-77e7-4e58-ac08-ca282aab57db"/>
<bqs:reference rdf:resource="rdf:#7b422dfe-2fa6-4552-81e9-1b7af3b5a536"/>
</rdf:Description>
<rdf:Description rdf:about="#h">
<dcterms:alternative>GH</dcterms:alternative>
<dcterms:alternative>growth hormone</dcterms:alternative>
<dc:title>h</dc:title>
</rdf:Description>
<rdf:Description rdf:about="#r">
<dcterms:alternative>GHRH</dcterms:alternative>
<dcterms:alternative>growth hormone releasing hormone</dcterms:alternative>
<dc:title>r</dc:title>
</rdf:Description>
<rdf:Description rdf:about="#s">
<dcterms:alternative>endogenous somatostatin</dcterms:alternative>
<dc:title>s</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7973be78-1364-49f9-8098-5a9724200a00">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#26b2d58c-fe42-4106-9f62-9caec4ebd5de">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#abd20f21-286b-40ad-9ab6-f9ac642238b6">
<dc:title>Journal of Neuroendocrinology</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#790fdbd3-9a41-4af0-bea7-89f7be75237c">
<dc:creator rdf:resource="rdf:#26b2d58c-fe42-4106-9f62-9caec4ebd5de"/>
<rdf:value>The model now runs in PCEnv to replicate the published results (as shown in figure 12). The only difference is my stimulus was hand written by a piecewise equation where as the model authors used some sort of a random pulse generator... I probably need to make the piecewise equation carry on for longer to get this result.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9ae82fcb-97af-4380-a8a3-0d71d815b744">
<dc:creator rdf:resource="rdf:#c9b84d74-3bc4-4c84-aa34-bf92fb3d6e2a"/>
<dc:title>Modelling the hypothalamic control of growth hormone secretion.</dc:title>
<bqs:volume>17</bqs:volume>
<bqs:first_page>788</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#abd20f21-286b-40ad-9ab6-f9ac642238b6"/>
<dcterms:issued rdf:resource="rdf:#73ab6072-f902-4e94-8cc9-ae0e3dd3f257"/>
<bqs:last_page>803</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0074e35a-7320-417c-8ffd-da04aa3f74ca">
<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:#39694593-b277-45a7-b193-2a5668a1a2cb">
<vCard:Given>D</vCard:Given>
<vCard:Family>MacGregor</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7b422dfe-2fa6-4552-81e9-1b7af3b5a536">
<dc:subject rdf:resource="rdf:#5151e189-41bc-42ea-bcce-66b4ba368648"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#738722ff-ed9d-472c-a085-f3433a6fb047">
<vCard:N rdf:resource="rdf:#be2f82ba-dcbc-4733-94ee-0ecd2fddbdd5"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f0d77927-6fdc-4fb1-8d74-1a4fd5df0de0">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5151e189-41bc-42ea-bcce-66b4ba368648">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value rdf:resource="rdf:#21108480-0177-4661-aa67-a32d8cc46324"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a8dfbb63-684f-4412-8aad-4afa2b5b47f7">
<dc:creator rdf:resource="rdf:#f0d77927-6fdc-4fb1-8d74-1a4fd5df0de0"/>
<rdf:value>This is the CellML description of MacGregor and Leng's 2005 mathematical model of hypothalamic control of growth hormone secretion.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f52889a3-8e66-4864-aef2-88c1f43b54e7">
<vCard:Given>G</vCard:Given>
<vCard:Family>Leng</vCard:Family>
</rdf:Description>
</rdf:RDF>
</model>
