Location: Smith, 1980 @ a7736186bc03 / smith_1980.cellml

Author:
Geoff Nunns <gnunns1@jhu.edu>
Date:
2010-03-23 12:53:30+13:00
Desc:
Units inconsistencies have been corrected, and the model now runs to recreate published results.
Permanent Source URI:
http://models.cellml.org/workspace/smith_1980/rawfile/a7736186bc0380ebe3c83593333fcfd11cddfad8/smith_1980.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 23/03/2010 at 12:39:00 at p.m. using:

COR (0.9.31.1351)
Copyright 2002-2010 Dr Alan Garny
http://cor.physiol.ox.ac.uk/ - cor@physiol.ox.ac.uk

CellML 1.0 was used to generate this model
http://www.cellml.org/
-->
<model 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:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="smith_model_1980_version01" cmeta:id="smith_model_1980_version01" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Hypothalamic regulation of pituitary secretion of luteinizing hormone. II. Feedback control of gonadotropin 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 is valid CellML and runs in PCEnv, OpenCell and COR. The units have been checked and are consistent. The model itself runs to reproduce published results within a margin of error, which was judged to be acceptable, as the model authors explicitly state that the figures reproduced in the paper are only representative of the actual results. The CellML model is very small, and is unlikely to contain errors, so the current version can be considered a faithful reproduction of the published model. 
          </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
In human males the secretion of luteinizing hormone (LH) from the anterior pituitary gland is principally controlled by luteinizing hormone releasing hormone (LHRH) which, itself, is secreted by the hypothalamus.  LH in turn influences the secretion of testosterone by the testes, which then has a feedback effect on secretion of its precursor hormones LHRH and LH. 
</para>
				<para>
In the paper described here, William Smith proposes a simple mathematical model incorporating a negative feedback loop to describe the kinetics of the hormone cascade in the secretion of testosterone (as summarised in the figure below).
</para>
				<informalfigure float="0" id="fig_reaction_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>model diagram</title>
							</objectinfo>
							<imagedata fileref="smith_1980.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the mathematical model representing the male hypothalamic-pituitary-gonadal system.  The model has three main components: the hypothalamus which secretes luteinizing hormone releasing hormone (LHRH), the pituitary which responds to LHRH by secreting luteinizing hormone (LH), and the gonads (or 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/pu11363366022540/?p=3e13cfd941d04c37b73466066437cd18&amp;pi=2">Hypothalamic regulation of pituitary secretion of luteinizing hormone. II. Feedback control of gonadotropin secretion</ulink>, William R. Smith, 1980, <ulink url="http://www.springerlink.com/content/119979/">
						<emphasis>Bulletin of Mathematical Biology</emphasis>
					</ulink>, 42, 57-78.  (A <ulink url="http://www.springerlink.com/content/pu11363366022540/fulltext.pdf">PDF</ulink> version of the article is available to journal subscribers on the <emphasis>Bulletin of Mathematical Biology</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=6986927&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 6986927</ulink>
				</para>
			</sect1>
		</article>
	</documentation>
	<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>
	<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>
	<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>
	<rdf:RDF>
		<rdf:Bag rdf:about="rdf:#ee9b163e-6d7c-468d-9f30-8e26f59db0c3">
			<rdf:li>testis</rdf:li>
			<rdf:li>testosterone</rdf:li>
			<rdf:li>luteinizing hormone</rdf:li>
			<rdf:li>endocrine</rdf:li>
			<rdf:li>hypothalamus</rdf:li>
			<rdf:li>pituitary</rdf:li>
		</rdf:Bag>
		<rdf:Seq rdf:about="rdf:#b7c21e02-fd1e-461b-b356-fd002596e37e">
			<rdf:li rdf:resource="rdf:#9ad4fa77-19f7-494c-9df5-455f230723ed"/>
		</rdf:Seq>
		<rdf:Description rdf:about="rdf:#26258ce6-dbc8-4fc5-9c24-b823f294e7e0">
			<dcterms:modified rdf:resource="rdf:#b6bdc0e4-297d-4643-b948-1caa7bfc5d1e"/>
			<rdf:value>I've added initial values and have checked and corrected the units such that they are now consistent.</rdf:value>
			<cmeta:modifier xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#f4753eff-31bc-48bc-a939-ade3db8a3768"/>
		</rdf:Description>
		<rdf:Description rdf:about="">
			<dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
			<cmeta:comment xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#a370189b-7aee-4249-bb79-981ee248016f"/>
			<dcterms:created rdf:resource="rdf:#f45acd09-359b-4e34-b337-8e3b9ef7137f"/>
			<dc:creator rdf:resource="rdf:#f79cd3cf-a7b3-42fc-a1b2-1443eb805e87"/>
			<cmeta:modification xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#26258ce6-dbc8-4fc5-9c24-b823f294e7e0"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#63bd3c41-e8e2-4e52-8dcd-b4c2fa147500">
			<bqs:Pubmed_id>6986927</bqs:Pubmed_id>
			<bqs:JournalArticle rdf:resource="rdf:#f7ae5e2e-8685-4d4d-ab3a-ea21e57215ef"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#f45acd09-359b-4e34-b337-8e3b9ef7137f">
			<dcterms:W3CDTF>2007-07-18T00:00:00+00:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#b7d9ef22-1065-43e5-8391-f15f0c1de290">
			<bqs:subject_type>keyword</bqs:subject_type>
			<rdf:value rdf:resource="rdf:#ee9b163e-6d7c-468d-9f30-8e26f59db0c3"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#b6bdc0e4-297d-4643-b948-1caa7bfc5d1e">
			<dcterms:W3CDTF>2007-09-06T09:34:21+12:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#9ad4fa77-19f7-494c-9df5-455f230723ed">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#4fbdb684-ce08-48d7-b5a0-b804e6774d38"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#a370189b-7aee-4249-bb79-981ee248016f">
			<dc:creator rdf:resource="rdf:#61340964-d357-43ef-a333-273bcfc57293"/>
			<rdf:value>This model runs in PCEnv to replicate the published results.  The units have been checked and are consistent.</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#a2f33434-bbd5-43d0-8ca8-f2104a1f9afb">
			<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:#4288147f-af4d-4e89-8d96-5203ab00fdf4">
			<dc:subject rdf:resource="rdf:#b7d9ef22-1065-43e5-8391-f15f0c1de290"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#cc37bafa-43f6-4e07-ab3e-fd5e93d1b6fe">
			<vCard:Given>Catherine</vCard:Given>
			<vCard:Family>Lloyd</vCard:Family>
			<vCard:Other>May</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#ed99eaa2-229a-4153-9e0d-ef69c8145b19">
			<dcterms:W3CDTF>1980-00-00 00:00</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#cc97b9ae-0519-4ffd-b657-6189fb12cee3">
			<vCard:Given>Catherine</vCard:Given>
			<vCard:Family>Lloyd</vCard:Family>
			<vCard:Other>May</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#ad2ef8fd-ca45-4e88-9bbb-1f95244e718b">
			<dc:title>Bulletin of Mathematical Biology</dc:title>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#23578751-dd48-4aa9-a929-de8f5d2a1ba6">
			<dc:creator rdf:resource="rdf:#1b1261a6-0a99-4c18-ad83-8adedb82d75f"/>
			<rdf:value>This is a CellML description of Smith's 1980 mathematical model for the hypothalamic regulation of pituitary secretion of luteinizing hormone - feedback control of gonadotropin secretion.</rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#4fbdb684-ce08-48d7-b5a0-b804e6774d38">
			<vCard:Given>William</vCard:Given>
			<vCard:Family>Smith</vCard:Family>
			<vCard:Other>R</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#61340964-d357-43ef-a333-273bcfc57293">
			<vCard:FN>Catherine Lloyd</vCard:FN>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#f7ae5e2e-8685-4d4d-ab3a-ea21e57215ef">
			<dc:creator rdf:resource="rdf:#b7c21e02-fd1e-461b-b356-fd002596e37e"/>
			<dc:title>Hypothalamic regulation of pituitary secretion of luteinizing hormone. II. Feedback control of gonadotropin secretion</dc:title>
			<bqs:volume>42</bqs:volume>
			<bqs:first_page>57</bqs:first_page>
			<bqs:Journal rdf:resource="rdf:#ad2ef8fd-ca45-4e88-9bbb-1f95244e718b"/>
			<dcterms:issued rdf:resource="rdf:#ed99eaa2-229a-4153-9e0d-ef69c8145b19"/>
			<bqs:last_page>78</bqs:last_page>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1b1261a6-0a99-4c18-ad83-8adedb82d75f">
			<vCard:FN>Catherine Lloyd</vCard:FN>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#f4753eff-31bc-48bc-a939-ade3db8a3768">
			<vCard:N rdf:resource="rdf:#cc97b9ae-0519-4ffd-b657-6189fb12cee3"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#f79cd3cf-a7b3-42fc-a1b2-1443eb805e87">
			<vCard:ORG rdf:resource="rdf:#3acc86f6-db41-40ec-9b44-cc226f6e19b1"/>
			<vCard:EMAIL rdf:resource="rdf:#a2f33434-bbd5-43d0-8ca8-f2104a1f9afb"/>
			<vCard:N rdf:resource="rdf:#cc37bafa-43f6-4e07-ab3e-fd5e93d1b6fe"/>
		</rdf:Description>
		<rdf:Description rdf:about="#smith_model_1980">
			<dc:title>
        Smith's 1980 mathematical model for the hypothalamic regulation of pituitary secretion of luteinizing hormone - feedback control of gonadotropin secretion.  
      </dc:title>
			<cmeta:comment xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#23578751-dd48-4aa9-a929-de8f5d2a1ba6"/>
			<bqs:reference rdf:resource="rdf:#4288147f-af4d-4e89-8d96-5203ab00fdf4"/>
			<bqs:reference rdf:resource="rdf:#63bd3c41-e8e2-4e52-8dcd-b4c2fa147500"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#3acc86f6-db41-40ec-9b44-cc226f6e19b1">
			<vCard:Orgname>The University of Auckland</vCard:Orgname>
			<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
		</rdf:Description>
	</rdf:RDF>   <units xmlns="http://www.cellml.org/cellml/1.0#" name="hour">
      <unit units="second" multiplier="3600"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="ng_ml">
      <unit units="gram" prefix="nano"/>
      <unit units="litre" prefix="milli" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="ng_ml_hr">
      <unit units="gram" prefix="nano"/>
      <unit units="litre" prefix="milli" exponent="-1"/>
      <unit units="hour" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="first_order_rate_constant">
      <unit units="hour" exponent="-1"/>
   </units>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="environment">
      <variable name="time" units="hour" public_interface="out"/>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="R">
      <variable name="R" units="ng_ml" initial_value="0.5" public_interface="out"/>
      <variable name="h" units="first_order_rate_constant" initial_value="12"/>
      <variable name="c" units="ng_ml_hr" initial_value="100"/>
      <variable name="b1" units="first_order_rate_constant" initial_value="1.29"/>
      <variable name="H" units="dimensionless"/>
      <variable name="x" units="ng_ml"/>
      <variable name="T" units="ng_ml" public_interface="in"/>
      <variable name="time" units="hour" public_interface="in"/>
      <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/>
                  <apply>
                     <minus/>
                     <ci>c</ci>
                     <apply>
                        <times/>
                        <ci>h</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>H</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>b1</ci>
                  <ci>R</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>x</ci>
            <apply>
               <minus/>
               <ci>T</ci>
               <apply>
                  <divide/>
                  <ci>c</ci>
                  <ci>h</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>H</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <gt/>
                     <ci>x</ci>
                     <cn cellml:units="ng_ml">0</cn>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="dimensionless">0</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="L">
      <variable name="L" units="ng_ml" initial_value="22" public_interface="out"/>
      <variable name="g1" units="first_order_rate_constant" initial_value="10"/>
      <variable name="b2" units="first_order_rate_constant" initial_value="0.97"/>
      <variable name="R" units="ng_ml" public_interface="in"/>
      <variable name="time" units="hour" public_interface="in"/>
      <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 xmlns="http://www.cellml.org/cellml/1.0#" name="T">
      <variable name="T" units="ng_ml" initial_value="15" public_interface="out"/>
      <variable name="g2" units="first_order_rate_constant" initial_value="0.7"/>
      <variable name="b3" units="first_order_rate_constant" initial_value="1.39"/>
      <variable name="L" units="ng_ml" public_interface="in"/>
      <variable name="time" units="hour" public_interface="in"/>
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="R" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="L" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="T" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="R" component_2="T"/>
      <map_variables variable_1="T" variable_2="T"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="R" component_2="L"/>
      <map_variables variable_1="R" variable_2="R"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="T" component_2="L"/>
      <map_variables variable_1="L" variable_2="L"/>
   </connection>
</model>