Rendering of the source text

<?xml version="1.0"?>
<!--  FILE :  enciso_sontag_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#" xml:base="" cmeta:id="enciso_sontag_2004" name="enciso_sontag_2004">

<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>Auckland Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>

			<section id="sec_status">
				<title>Model Status</title>
				<para>
           This model is known to run in both OpenCell and COR, although there are no figures in the paper with which to compare the CellML model simulation output. The units have been checked and they are consistent.
          </para>
			</section>

  <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>
On the stability of a model of testosterone dynamics, German Enciso and Eduardo D. Sontag, 2004, 
            <emphasis>Journal of Mathematical Biology</emphasis>, 49, 627-634. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/15480670">PubMed ID: 15480670</ulink>
</para>

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





  <units name="hour">
    <unit multiplier="3600.0" units="second"/>
  </units>

  <units name="ng_ml">
    <unit prefix="nano" units="gram"/>
    <unit exponent="-1" prefix="milli" units="litre"/>
  </units>

  <units name="ng_ml_hr">
    <unit units="ng_ml"/>
    <unit exponent="-1" units="hour"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit exponent="-1" units="hour"/>
  </units>

  
  <component name="environment">
    <variable cmeta:id="environment_time" name="time" public_interface="out" units="hour"/>
  </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 xmlns:dc="http://purl.org/dc/elements/1.1/">R</dc:title>
        <dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">LHRH</dcterms:alternative>
        <dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">luteinizing hormone releasing hormone</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable cmeta:id="R_R" initial_value="0.1" name="R" public_interface="out" units="ng_ml"/>
    
    <variable initial_value="1.0" name="b1" units="first_order_rate_constant"/>   
    <variable initial_value="1.0" name="b2" units="dimensionless"/>
    <variable initial_value="1.0" name="b3" units="dimensionless"/>
	<variable initial_value="10.0" name="g1" units="dimensionless"/>
    <variable initial_value="10.0" name="g2" units="dimensionless"/>
    <variable name="f_T" units="ng_ml_hr"/>
    
    <variable name="T" public_interface="in" units="ng_ml"/>
    <variable name="time" public_interface="in" units="hour"/>
    
    <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 xmlns:dc="http://purl.org/dc/elements/1.1/">L</dc:title>
        <dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">LH</dcterms:alternative>
        <dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">luteinizing hormone</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable cmeta:id="L_L" initial_value="0.0" name="L" public_interface="out" units="ng_ml"/>
    
    <variable initial_value="1.0" name="b2" units="first_order_rate_constant"/>
    <variable initial_value="10.0" name="g1" units="first_order_rate_constant"/>
    
    <variable name="R" public_interface="in" units="ng_ml"/>
    <variable name="time" public_interface="in" units="hour"/>
    
    <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 xmlns:dc="http://purl.org/dc/elements/1.1/">T</dc:title>
        <dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">testosterone</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable cmeta:id="T_T" initial_value="0.0" name="T" public_interface="out" units="ng_ml"/>
    
    <variable initial_value="1.0" name="b3" units="first_order_rate_constant"/>
    <variable initial_value="10.0" name="g2" units="first_order_rate_constant"/>
    
    <variable name="L" public_interface="in" units="ng_ml"/>
    <variable name="time" public_interface="in" units="hour"/>
    
    <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_1="R" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  
  <connection>
    <map_components component_1="L" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>

  <connection>
    <map_components component_1="T" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  
  <connection>
    <map_components component_1="R" component_2="L"/>
    <map_variables variable_1="R" variable_2="R"/>
  </connection>

  <connection>
    <map_components component_1="R" component_2="T"/>
    <map_variables variable_1="T" variable_2="T"/>
  </connection>

  <connection>
    <map_components component_1="T" component_2="L"/>
    <map_variables variable_1="L" variable_2="L"/>
  </connection>

<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <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 xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#d537c69f-001d-4ee2-87db-71d8463e1671"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher xmlns:dc="http://purl.org/dc/elements/1.1/">The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#b7e843f1-f17b-4873-a16f-df929daa6036"/>
    <dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#f2840ea7-e126-46f1-90d0-4cb1e5e4ee6d"/>
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" 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 xmlns:dcterms="http://purl.org/dc/terms/">2007-09-02T11:55:07+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#395750ec-e22a-4c22-8f9f-41e96e4e5e64">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:resource="rdf:#6d7384ff-aea0-44b0-960e-827ba559c87b"/>
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">On the stability of a model of testosterone dynamics</dc:title>
    <bqs:volume xmlns:bqs="http://www.cellml.org/bqs/1.0#">49</bqs:volume>
    <bqs:first_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">627</bqs:first_page>
    <bqs:Journal xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#cca55c33-50f4-4561-b872-bc572109e5d0"/>
    <dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:resource="rdf:#57c844e1-29d8-45dd-b908-4b60dad35382"/>
    <bqs:last_page xmlns:bqs="http://www.cellml.org/bqs/1.0#">634</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#43424633-8e66-4dc0-9327-18aa474312c8">
    <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#6dd7a57f-1fea-4986-8fd9-70f9c27b4f7a"/>
    <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#90e6deff-9629-48a1-967d-cc4b4cfcf1e2"/>
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#432f2cfa-0962-4698-855e-80eab7d5583a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#57c844e1-29d8-45dd-b908-4b60dad35382">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">2004-12-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b7e843f1-f17b-4873-a16f-df929daa6036">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" 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 xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lloyd</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5aa6836b-57b7-47fc-b543-9eaf5f2f5c81">
    <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" 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 xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bcd0fff9-7935-49b5-aaca-5f1d49b480eb">
    <dcterms:modified xmlns:dcterms="http://purl.org/dc/terms/" 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 xmlns:dc="http://purl.org/dc/elements/1.1/">Journal of Mathematical Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#15908299-75d5-4cec-905a-14b137a74378">
    <vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1a5b9d5c-1cd9-4ad8-8d71-ebfedf8060e4">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Eduardo</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Sontag</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">D</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f2840ea7-e126-46f1-90d0-4cb1e5e4ee6d">
    <dcterms:W3CDTF xmlns:dcterms="http://purl.org/dc/terms/">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 xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#1a5b9d5c-1cd9-4ad8-8d71-ebfedf8060e4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9a34946e-1b99-482c-9502-b4ead4969696">
    <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#166bf7e8-88ff-4348-ba56-b13098bd3a14"/>
  </rdf:Description>
  <rdf:Description rdf:about="#enciso_sontag_2004">
    <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
        Enciso and Sontag's 2004 mathematical model of testosterone dynamics.  
      </dc:title>
    <cmeta:comment rdf:resource="rdf:#5aa6836b-57b7-47fc-b543-9eaf5f2f5c81"/>
    <bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#d614fd1a-0f6c-4b47-b8cc-ae19080b0f0e"/>
    <bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" 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 xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:resource="rdf:#7f0eea60-e115-4c27-be42-196f60ae50a5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7f0eea60-e115-4c27-be42-196f60ae50a5">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">German</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Enciso</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6dd7a57f-1fea-4986-8fd9-70f9c27b4f7a">
    <vCard:Orgname xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">The University of Auckland</vCard:Orgname>
    <vCard:Orgunit xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#432f2cfa-0962-4698-855e-80eab7d5583a">
    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Catherine</vCard:Given>
    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Lloyd</vCard:Family>
    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e61f84da-402b-435e-b712-20b5d46500f3">
    <bqs:Pubmed_id xmlns:bqs="http://www.cellml.org/bqs/1.0#">15480670</bqs:Pubmed_id>
    <bqs:JournalArticle xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:resource="rdf:#395750ec-e22a-4c22-8f9f-41e96e4e5e64"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d537c69f-001d-4ee2-87db-71d8463e1671">
    <bqs:subject_type xmlns:bqs="http://www.cellml.org/bqs/1.0#">keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#f233d76d-0ca6-4d21-a7f1-b8691c2ec7de"/>
  </rdf:Description>
</rdf:RDF>
</model>