- Author:
- pmr2.import <nobody@models.cellml.org>
- Date:
- 2009-06-17 13:42:41+12:00
- Desc:
- committing version02 of enciso_sontag_2004
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/enciso_sontag_2004/rawfile/04197cebc58f2d1ef6e7dbae35ab1e6ddc651473/enciso_sontag_2004.cellml
<?xml version='1.0' encoding='utf-8'?>
<!-- FILE : enciso_model_2004.xml
CREATED : 18th July 2007
LAST MODIFIED : 18th 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 Enciso and Sontag's 2004 mathematical model of testosterone dynamics.
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="enciso_model_2004" name="enciso_sontag_model_2004_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>On the stability of a model of testosterone dynamics</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>
Blood testosterone levels in human males are known to oscillate periodically with a period of a few hours. Testosterone secretion from the testes is stimulated by luteinizing hormone (LH) and the production of LH is, in turn, stimulated by the hypothalamic hormone luteinizing hormone releasing hormone (LHRH) (as summarised in the figure below). Interestingly, the observed fluctuations in the blood levels of testosterone can be traced back to the pulsatile release of LHRH. Several mathematical models have been proposed to describe the pulsatile release of these three hormones. Some of these models depend on an introduced delay in one of the variables in order to create oscillations in an otherwise stable model. Examples of such delayed mathematical models include those published by Cartwright and Husain (1980) which describes the control of testosterone secretion, and also the model published by Lenbury and Pornsawad (2005) to characterise feedback control in the hypothalamus-pituitary-adrenal axis in humans. Both of these models have been translated into CellML and can be found in the model repository. However, while delayed models can be biologically realistic (for example there is often a delay between the secretion of a hormone and the binding of that hormone to its receptor and the subsequent response), introducing a delay to a mathematical model often comes at the expense of increased model complexity. (In particular, it should be noted that in its current form, CellML is unable to express time delays).
</para>
<para>
In the paper described here, Enciso and Sontag critique published mathematical models which describe testosterone dynamics. In particular they highlight how models which use time delays to create hormone oscillations can exhibit stability.
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="enciso_2004.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram of the mathematical model representing the luteinizing hormone releasing hormone (LHRH), luteinizing hormone (LH), and testosterone system. The model has three main components: the hypothalamus which secretes LHRH, the pituitary which responds to LHRH by secreting LH, and the testes which, in responds to LH, secrete testosterone. Finally, LH and testosterone feedback on the hypothalamus to regulate LHRH secretion.</caption>
</informalfigure>
<para>
The complete original paper reference is cited below:
</para>
<para>
<ulink url="http://www.springerlink.com/content/4hajnqthj9nyk7lk/">On the stability of a model of testosterone dynamics</ulink>, German Enciso and Eduardo D. Sontag, 2004, <ulink url="http://www.springerlink.com/content/100436/">
<emphasis>Journal of Mathematical Biology</emphasis>
</ulink>, 49, 627-634. (A<ulink url="http://www.springerlink.com/content/4hajnqthj9nyk7lk/fulltext.pdf">PDF</ulink> version of the article is available to journal subscribers on the <emphasis>Journal of Mathematical Biology</emphasis> website.) <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&cmd=Retrieve&dopt=AbstractPlus&list_uids=15480670&query_hl=1&itool=pubmed_docsum">PubMed ID: 15480670</ulink>
</para>
</sect1>
</article>
</documentation>
<units name="hour">
<unit units="second" multiplier="3600.0"/>
</units>
<units name="ng_ml">
<unit units="gram" prefix="nano"/>
<unit units="litre" prefix="milli" exponent="-1"/>
</units>
<units name="ng_ml_hr">
<unit units="ng_ml"/>
<unit units="hour" exponent="-1"/>
</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 cmeta:id="R" name="R">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="R">
<dc:title>R</dc:title>
<dcterms:alternative>LHRH</dcterms:alternative>
<dcterms:alternative>luteinizing hormone releasing hormone</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="ng_ml" public_interface="out" name="R" initial_value="0.1"/>
<variable units="first_order_rate_constant" name="b1" initial_value="1.0"/>
<variable units="dimensionless" name="b2" initial_value="1.0"/>
<variable units="dimensionless" name="b3" initial_value="1.0"/>
<variable units="dimensionless" name="g1" initial_value="10.0"/>
<variable units="dimensionless" name="g2" initial_value="10.0"/>
<variable units="ng_ml_hr" name="f_T"/>
<variable units="ng_ml" public_interface="in" name="T"/>
<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/>
<ci> f_T </ci>
<apply>
<times/>
<ci> b1 </ci>
<ci> R </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> f_T </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> b1 </ci>
<ci> b2 </ci>
<ci> b3 </ci>
<ci> T </ci>
</apply>
<apply>
<times/>
<ci> g1 </ci>
<ci> g2 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="L" name="L">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="L">
<dc:title>L</dc:title>
<dcterms:alternative>LH</dcterms:alternative>
<dcterms:alternative>luteinizing hormone</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="ng_ml" public_interface="out" name="L" initial_value="0.0"/>
<variable units="first_order_rate_constant" name="b2" initial_value="1.0"/>
<variable units="first_order_rate_constant" name="g1" initial_value="10.0"/>
<variable units="ng_ml" 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> L </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> g1 </ci>
<ci> R </ci>
</apply>
<apply>
<times/>
<ci> b2 </ci>
<ci> L </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="T" name="T">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="T">
<dc:title>T</dc:title>
<dcterms:alternative>testosterone</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="ng_ml" public_interface="out" name="T" initial_value="0.0"/>
<variable units="first_order_rate_constant" name="b3" initial_value="1.0"/>
<variable units="first_order_rate_constant" name="g2" initial_value="10.0"/>
<variable units="ng_ml" public_interface="in" name="L"/>
<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> T </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> g2 </ci>
<ci> L </ci>
</apply>
<apply>
<times/>
<ci> b3 </ci>
<ci> T </ci>
</apply>
</apply>
</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="L"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="T"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="L" component_1="R"/>
<map_variables variable_2="R" variable_1="R"/>
</connection>
<connection>
<map_components component_2="T" component_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
</connection>
<connection>
<map_components component_2="L" component_1="T"/>
<map_variables variable_2="L" variable_1="L"/>
</connection>
<rdf:RDF>
<rdf:Bag rdf:about="rdf:#f233d76d-0ca6-4d21-a7f1-b8691c2ec7de">
<rdf:li>testosterone</rdf:li>
<rdf:li>luteinizing hormone</rdf:li>
<rdf:li>hypothalamus</rdf:li>
<rdf:li>pituitary</rdf:li>
<rdf:li>testis</rdf:li>
<rdf:li>endocrine</rdf:li>
</rdf:Bag>
<rdf:Seq rdf:about="rdf:#6d7384ff-aea0-44b0-960e-827ba559c87b">
<rdf:li rdf:resource="rdf:#fb0539ce-323c-4ef9-b473-b578319fc6d9"/>
<rdf:li rdf:resource="rdf:#8e077a3e-4dd6-4a82-adf3-1f38d8cdceb5"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#d614fd1a-0f6c-4b47-b8cc-ae19080b0f0e">
<dc:subject rdf:resource="rdf:#d537c69f-001d-4ee2-87db-71d8463e1671"/>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
<cmeta:comment rdf:resource="rdf:#b7e843f1-f17b-4873-a16f-df929daa6036"/>
<dcterms:created rdf:resource="rdf:#f2840ea7-e126-46f1-90d0-4cb1e5e4ee6d"/>
<dc:creator rdf:resource="rdf:#43424633-8e66-4dc0-9327-18aa474312c8"/>
<cmeta:modification rdf:resource="rdf:#bcd0fff9-7935-49b5-aaca-5f1d49b480eb"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ca3f5590-4c10-4f71-8b60-4c45863575e1">
<dcterms:W3CDTF>2007-09-02T11:55:07+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#395750ec-e22a-4c22-8f9f-41e96e4e5e64">
<dc:creator rdf:resource="rdf:#6d7384ff-aea0-44b0-960e-827ba559c87b"/>
<dc:title>On the stability of a model of testosterone dynamics</dc:title>
<bqs:volume>49</bqs:volume>
<bqs:first_page>627</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#cca55c33-50f4-4561-b872-bc572109e5d0"/>
<dcterms:issued rdf:resource="rdf:#57c844e1-29d8-45dd-b908-4b60dad35382"/>
<bqs:last_page>634</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#43424633-8e66-4dc0-9327-18aa474312c8">
<vCard:ORG rdf:resource="rdf:#6dd7a57f-1fea-4986-8fd9-70f9c27b4f7a"/>
<vCard:EMAIL rdf:resource="rdf:#90e6deff-9629-48a1-967d-cc4b4cfcf1e2"/>
<vCard:N rdf:resource="rdf:#432f2cfa-0962-4698-855e-80eab7d5583a"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#57c844e1-29d8-45dd-b908-4b60dad35382">
<dcterms:W3CDTF>2004-12-00 00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#b7e843f1-f17b-4873-a16f-df929daa6036">
<dc:creator rdf:resource="rdf:#30530678-0b48-4b01-aef6-8a97bd0b4dd8"/>
<rdf:value>I've fixed the units using COR and they are now consistent. The model runs in PCEnv but I can't get it to oscillate - I'm guessing because I can't translate the time delays which are included in the model equations.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#166bf7e8-88ff-4348-ba56-b13098bd3a14">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5aa6836b-57b7-47fc-b543-9eaf5f2f5c81">
<dc:creator rdf:resource="rdf:#15908299-75d5-4cec-905a-14b137a74378"/>
<rdf:value>This is a CellML description of Enciso and Sontag's 2004 mathematical model of testosterone dynamics.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#30530678-0b48-4b01-aef6-8a97bd0b4dd8">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#bcd0fff9-7935-49b5-aaca-5f1d49b480eb">
<dcterms:modified rdf:resource="rdf:#ca3f5590-4c10-4f71-8b60-4c45863575e1"/>
<rdf:value>I've fixed the units using COR and they are now consistent. The model runs in PCEnv but I can't get it to oscillate - I'm guessing because I can't translate the time delays which are included in the model equations.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#9a34946e-1b99-482c-9502-b4ead4969696"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#90e6deff-9629-48a1-967d-cc4b4cfcf1e2">
<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:#cca55c33-50f4-4561-b872-bc572109e5d0">
<dc:title>Journal of Mathematical Biology</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#15908299-75d5-4cec-905a-14b137a74378">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#1a5b9d5c-1cd9-4ad8-8d71-ebfedf8060e4">
<vCard:Given>Eduardo</vCard:Given>
<vCard:Family>Sontag</vCard:Family>
<vCard:Other>D</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f2840ea7-e126-46f1-90d0-4cb1e5e4ee6d">
<dcterms:W3CDTF>2007-07-18T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#8e077a3e-4dd6-4a82-adf3-1f38d8cdceb5">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#1a5b9d5c-1cd9-4ad8-8d71-ebfedf8060e4"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9a34946e-1b99-482c-9502-b4ead4969696">
<vCard:N rdf:resource="rdf:#166bf7e8-88ff-4348-ba56-b13098bd3a14"/>
</rdf:Description>
<rdf:Description rdf:about="#enciso_model_2004">
<dc:title>
Enciso and Sontag's 2004 mathematical model of testosterone dynamics.
</dc:title>
<cmeta:comment rdf:resource="rdf:#5aa6836b-57b7-47fc-b543-9eaf5f2f5c81"/>
<bqs:reference rdf:resource="rdf:#d614fd1a-0f6c-4b47-b8cc-ae19080b0f0e"/>
<bqs:reference rdf:resource="rdf:#e61f84da-402b-435e-b712-20b5d46500f3"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#fb0539ce-323c-4ef9-b473-b578319fc6d9">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#7f0eea60-e115-4c27-be42-196f60ae50a5"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#7f0eea60-e115-4c27-be42-196f60ae50a5">
<vCard:Given>German</vCard:Given>
<vCard:Family>Enciso</vCard:Family>
</rdf:Description>
<rdf:Description rdf:about="rdf:#6dd7a57f-1fea-4986-8fd9-70f9c27b4f7a">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#432f2cfa-0962-4698-855e-80eab7d5583a">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e61f84da-402b-435e-b712-20b5d46500f3">
<bqs:Pubmed_id>15480670</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#395750ec-e22a-4c22-8f9f-41e96e4e5e64"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#d537c69f-001d-4ee2-87db-71d8463e1671">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value rdf:resource="rdf:#f233d76d-0ca6-4d21-a7f1-b8691c2ec7de"/>
</rdf:Description>
</rdf:RDF>
</model>
