Rendering of the source text

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  cartwright_model_1986.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 Cartwright and Husain's mathematical model for the control of testosterone 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="cartwright_model_1986" name="cartwright_husain_model_1986_version01">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>A model for the control of testosterone 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 model can be solved and can reproduce one oscillation (which is true to the published results,) but is unable to produce further oscillations.
          </para>
<para>
ValidateCellML verifies this model as valid CellML, although unit inconsistencies are detected.
</para>
  </section>
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
Blood testosterone levels in human males fluctuate over 2-3 hour periods.  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).  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.  These models can be divided into two different classes: 
</para>

<itemizedlist>
  <listitem>
            <para>Those which assume the existence of a <emphasis>neural clock</emphasis>; a pulse-generator which forces the hypothalamic secretion of LHRH.</para>
          </listitem>
  <listitem>
            <para>And those which assume there is no external pulse input, but explain the hormone pulsatility through interactions between the components of the system which produce <emphasis>feedback oscillations</emphasis>.</para>
          </listitem>
</itemizedlist>  

<para>
In the paper described here, Cartwright and Husain propose mathematical model incorporating time delays and negative feedback loops to describe the kinetics of the hormone cascade in the secretion of testosterone.  Their model falls into the second category of models described above.  The model has three components; the hypothalamus, pituitary and testes, and the rate of LH release from the pituitary depends on the the LHRH concentration, and the rate of testosterone synthesis in the testes depends on the LH concentration (as summarised in the figure below).
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="cartwright_1986.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>
A model for the control of testosterone secretion, Mark Cartwright and Masud Husain, 1986, 
            <emphasis>Journal of Theoretical Biology</emphasis>
           123, 239-250.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=3306160&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 3306160</ulink>
</para>

<para>
It should be noted that in its current form, the CellML description of the this model is unable to perfectly capture the simulation results of the published model, this is due to the time delays which are difficult to describe in the CellML code.
</para>

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





  <units name="minute">
    <unit units="second" multiplier="60.0"/>
  </units>

  <units name="hour">
    <unit units="second" multiplier="3600.0"/>
  </units>
  
  <units name="pg_ml">
    <unit units="gram" prefix="pico"/>
    <unit units="litre" prefix="milli" exponent="-1"/>
  </units>

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

  <units name="ng_ml_min">
    <unit units="ng_ml"/>
    <unit units="minute" exponent="-1"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="minute" exponent="-1"/>
  </units>

  
  <component name="environment">
    <variable units="hour" public_interface="out" name="time"/>
  </component>
  
  <component name="R" cmeta:id="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="pg_ml" public_interface="out" name="R" initial_value="0.0"/>
    
    <variable units="ng_ml" name="L_" initial_value="30.0"/>
    <variable units="ng_ml" name="T_" initial_value="8.0"/>
    <variable units="first_order_rate_constant" name="dR" initial_value="0.10"/>
    <variable units="ng_ml_min" name="rR" initial_value="0.1"/>
    <variable units="dimensionless" name="H"/>
    <variable units="dimensionless" name="x"/>
    
    <variable units="ng_ml" public_interface="in" name="L"/>
    <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/>
          <apply>
            <times/>
            <ci> rR </ci>
            <ci> H </ci>
          </apply>
          <apply>
            <times/>
            <ci> dR </ci>
            <ci> R </ci>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> x </ci> 
        <apply>
          <minus/>
          <cn cellml:units="dimensionless"> 2.0 </cn>
		  <apply>
            <plus/>
            <apply>
              <divide/>
              <ci> L </ci>
              <ci> L_ </ci>
			</apply>
            <apply>
              <divide/>
              <ci> T </ci>
              <ci> T_ </ci>
            </apply>
          </apply>
        </apply>
      </apply>

      <apply>
        <eq/> 
        <ci> H </ci> 
        <piecewise>
          <piece>
            <cn cellml:units="dimensionless"> 0.0 </cn>
            <apply>
              <lt/>
              <ci> x </ci>
              <cn cellml:units="dimensionless"> 0.0 </cn>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="dimensionless"> 0.5 </cn>
            <apply>
              <eq/>
              <ci> x </ci>
              <cn cellml:units="dimensionless"> 0.0 </cn>
            </apply>
          </piece>
          <piece>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <gt/>
              <ci> x </ci>
              <cn cellml:units="dimensionless"> 0.0 </cn>
            </apply>
          </piece>
        </piecewise>
      </apply>      
    </math>
  </component>

  <component name="L" cmeta:id="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="10.0"/>
    
    <variable units="first_order_rate_constant" name="dL" initial_value="0.015"/>
    <variable units="first_order_rate_constant" name="rL" initial_value="5.0"/>
    
    <variable units="pg_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> rL </ci>
            <ci> R </ci>
          </apply>
          <apply>
            <times/>
            <ci> dL </ci>
            <ci> L </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="T" cmeta:id="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="12.0"/>
    
    <variable units="first_order_rate_constant" name="dT" initial_value="0.023"/>
    <variable units="first_order_rate_constant" name="rT" initial_value="0.01"/>
    
    <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> rT </ci>
            <ci> L </ci>
          </apply>
          <apply>
            <times/>
            <ci> dT </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"/>
    <map_variables variable_2="L" variable_1="L"/>
  </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:#9936f484-8ba4-4dc6-b93c-de00042113fc">
    <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:#27e629ef-04ad-482e-b33d-577d847ae60d">
    <rdf:li rdf:resource="rdf:#48194ee7-cfee-4f72-a12e-1004473a57b4"/>
    <rdf:li rdf:resource="rdf:#d978cf74-1dfa-46e6-83ef-dbc41a51db3e"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#1f0ee94a-4131-4f14-8db0-6990eaf843db">
    <vCard:ORG rdf:resource="rdf:#5a49b4aa-24bd-4f8c-940c-dc459973e787"/>
    <vCard:EMAIL rdf:resource="rdf:#3b0ff7ef-daa6-44fb-9e49-f6cd48e3bb7f"/>
    <vCard:N rdf:resource="rdf:#7a3bc8ff-c54e-47f5-8918-25c297a67b75"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#084d6044-2ce4-4a5c-9a0f-7501c4abf427">
    <dc:creator rdf:resource="rdf:#f705fd3a-e999-458b-b4c1-1e5810f0f393"/>
    <rdf:value>The model has been curated - the model runs in PCEnv and COR and the units have been checked and are consistent.  The model produces one oscillation which is true to the published results.  However, I can't get it to repeatedly oscillate - I suspect because we need to include the time delays - which CellML can't handle at the moment.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#285840be-a3cc-401f-a037-ea7c13391f2e">
    <vCard:Given>Masud</vCard:Given>
    <vCard:Family>Husain</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#084d6044-2ce4-4a5c-9a0f-7501c4abf427"/>
    <dcterms:created rdf:resource="rdf:#89342f75-a680-4e83-8095-ce3ebd41edbe"/>
    <dc:creator rdf:resource="rdf:#1f0ee94a-4131-4f14-8db0-6990eaf843db"/>
    <cmeta:modification rdf:resource="rdf:#264eb2e6-bdf6-437d-be1f-0108b3a16b76"/>
    <cmeta:modification rdf:resource="rdf:#3fe3b4e8-0bd2-4b6a-b9bc-47c76feee8cd"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#39eab287-1254-4a44-a25e-724c7ad55de8">
    <dc:title>Journal of Theoretical Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#89342f75-a680-4e83-8095-ce3ebd41edbe">
    <dcterms:W3CDTF>2007-07-18T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f705fd3a-e999-458b-b4c1-1e5810f0f393">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4b4403c8-1bd1-4efe-8dc9-6fcb1065428e">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aded85e8-7a1b-4d69-80fa-cf1495b5c239">
    <vCard:N rdf:resource="rdf:#b295a485-18e0-4a94-aa67-9fcd33e5191f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#48194ee7-cfee-4f72-a12e-1004473a57b4">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#b85150af-e402-438e-8fee-a0008429ac0a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f6acbd35-7d42-4e1f-81d0-2dee9f6ceeb9">
    <dcterms:W3CDTF>2009-05-29T12:38:58+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b295a485-18e0-4a94-aa67-9fcd33e5191f">
    <vCard:Given>James</vCard:Given>
    <vCard:Family>Lawson</vCard:Family>
    <vCard:Other>Richard</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d978cf74-1dfa-46e6-83ef-dbc41a51db3e">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#285840be-a3cc-401f-a037-ea7c13391f2e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#264eb2e6-bdf6-437d-be1f-0108b3a16b76">
    <dcterms:modified rdf:resource="rdf:#f6acbd35-7d42-4e1f-81d0-2dee9f6ceeb9"/>
    <rdf:value>updated curation status, removed reference link in documentation</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#aded85e8-7a1b-4d69-80fa-cf1495b5c239"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f5452353-08a4-495a-b5a5-dcc286234780">
    <dcterms:W3CDTF>2007-09-01T16:48:56+12:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#69bf6d70-bef9-4fa9-a4cd-5300efe00880">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#627af779-1201-4385-baa8-1b9857aeebf3">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#9936f484-8ba4-4dc6-b93c-de00042113fc"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b85150af-e402-438e-8fee-a0008429ac0a">
    <vCard:Given>Mark</vCard:Given>
    <vCard:Family>Cartwright</vCard:Family>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3fe3b4e8-0bd2-4b6a-b9bc-47c76feee8cd">
    <dcterms:modified rdf:resource="rdf:#f5452353-08a4-495a-b5a5-dcc286234780"/>
    <rdf:value>The model has been curated - the model runs in PCEnv and COR and the units have been checked and are consistent.  The model produces one oscillation which is true to the published results.  However, I can't get it to repeatedly oscillate - I suspect because we need to include the time delays - which CellML can't handle at the moment.</rdf:value>
    <cmeta:modifier rdf:resource="rdf:#3b54fbb0-63bb-46e6-895a-afa8348abe48"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3b0ff7ef-daa6-44fb-9e49-f6cd48e3bb7f">
    <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:#737ecdc8-9057-407d-b742-2898bb7246dc">
    <dcterms:W3CDTF>1986-11-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d008054f-e2e4-4def-9001-4d3470586fb5">
    <dc:subject rdf:resource="rdf:#627af779-1201-4385-baa8-1b9857aeebf3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5a49b4aa-24bd-4f8c-940c-dc459973e787">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ba43c5fc-c64c-4990-83fe-42a9ca8e910b">
    <dc:creator rdf:resource="rdf:#27e629ef-04ad-482e-b33d-577d847ae60d"/>
    <dc:title>A model for the control of testosterone secretion</dc:title>
    <bqs:volume>123</bqs:volume>
    <bqs:first_page>239</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#39eab287-1254-4a44-a25e-724c7ad55de8"/>
    <dcterms:issued rdf:resource="rdf:#737ecdc8-9057-407d-b742-2898bb7246dc"/>
    <bqs:last_page>250</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="#cartwright_model_1986">
    <dc:title>
        Cartwright and Husain's mathematical model for the control of testosterone secretion. 
      </dc:title>
    <cmeta:comment rdf:resource="rdf:#a5c4bd95-7ceb-41d6-b8f0-e700863933f3"/>
    <bqs:reference rdf:resource="rdf:#4f5a4bfa-6d6e-40d1-a556-e73eb1a8a62a"/>
    <bqs:reference rdf:resource="rdf:#d008054f-e2e4-4def-9001-4d3470586fb5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3b54fbb0-63bb-46e6-895a-afa8348abe48">
    <vCard:N rdf:resource="rdf:#4b4403c8-1bd1-4efe-8dc9-6fcb1065428e"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a5c4bd95-7ceb-41d6-b8f0-e700863933f3">
    <dc:creator rdf:resource="rdf:#69bf6d70-bef9-4fa9-a4cd-5300efe00880"/>
    <rdf:value>This is a CellML description of Cartwright and Husain's mathematical model for the control of testosterone secretion.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7a3bc8ff-c54e-47f5-8918-25c297a67b75">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4f5a4bfa-6d6e-40d1-a556-e73eb1a8a62a">
    <bqs:Pubmed_id>3306160</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#ba43c5fc-c64c-4990-83fe-42a9ca8e910b"/>
  </rdf:Description>
</rdf:RDF>
</model>