Rendering of the source text

<?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="cartwright_model_1986" name="cartwright_model_1986">

<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>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This CellML model runs in OpenCell and can reproduce one oscillation (which is true to the published results) but is unable to produce further oscillations. This is probably because there are time delays in the origianl model which as yet can not be described in CellML.  The units have been checked and they are consistent with just small equivalencies in pico and nano grams per ml. The CellML model also runs in COR but the long duration of the simulation is not ideal for this tool.
          </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>
ABSTRACT: We produce here a model to explain the control of testosterone secretion. In this model the hypothalamic secretion of the hormone LHRH (luteinizing hormone releasing hormone) is controlled by a combination of local testosterone concentration and of the local concentration of the pituitary hormone LH (luteinizing hormone). Since LHRH stimulates the release of LH, and LH in turn stimulates the release of testosterone, the three hormones constitute a three-component "feedback" network. We show how this model is able to account for the pulsatility of the release of these three hormones. Furthermore, the model is consistent with results obtained from a wide range of experimental manipulations of the system. For example, it accounts for the changes observed in hormone release patterns after castration. In particular, it follows that no "neural clock", or "neural pulse-generator", is required to force the system into pulsatile behaviour.
</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 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/pubmed/3306160">PubMed ID: 3306160</ulink>
</para>

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





  <units name="minute">
    <unit units="second" multiplier="60.0"/>
  </units>
  
  <units name="pg_per_ml">
    <unit units="gram" prefix="pico"/>
    <unit units="litre" prefix="milli" exponent="-1"/>
  </units>

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

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

  
  <component name="environment">
    <variable units="minute" 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_per_ml" public_interface="out" name="R" initial_value="0.0"/>
    
    <variable units="ng_per_ml" name="L_" initial_value="30.0"/>
    <variable units="ng_per_ml" name="T_" initial_value="8.0"/>
    <variable units="first_order_rate_constant" name="dR" initial_value="0.10"/>
    <variable units="ng_per_ml_min" name="rR" initial_value="0.1"/>
    <variable units="dimensionless" name="H"/>
    <variable units="dimensionless" name="x"/>
    
    <variable units="ng_per_ml" public_interface="in" name="L"/>
    <variable units="ng_per_ml" public_interface="in" name="T"/>
    <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>
          <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_per_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_per_ml" public_interface="in" name="R"/>
    <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> 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_per_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_per_ml" public_interface="in" name="L"/>
    <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> 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, 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:#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">
    <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>Auckland 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>