On the stability of a model of testosterone dynamics
Catherine
Lloyd
Auckland Bioengineering Institute, The University of Auckland
Model Status
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.
Model Structure
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).
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.
model diagram
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.
The original paper reference is cited below:
On the stability of a model of testosterone dynamics, German Enciso and Eduardo D. Sontag, 2004,
Journal of Mathematical Biology, 49, 627-634. PubMed ID: 15480670
R
LHRH
luteinizing hormone releasing hormone
$\frac{d R}{d \mathrm{time}}=\mathrm{f\_T}-\mathrm{b1}R\mathrm{f\_T}=\frac{\mathrm{b1}\mathrm{b2}\mathrm{b3}T}{\mathrm{g1}\mathrm{g2}}$
L
LH
luteinizing hormone
$\frac{d L}{d \mathrm{time}}=\mathrm{g1}R-\mathrm{b2}L$
T
testosterone
$\frac{d T}{d \mathrm{time}}=\mathrm{g2}L-\mathrm{b3}T$
testosterone
luteinizing hormone
hypothalamus
pituitary
testis
endocrine
The University of Auckland, Auckland Bioengineering Institute
2007-09-02T11:55:07+12:00
On the stability of a model of testosterone dynamics
49
627
634
2004-12-00 00:00
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.
Catherine
Lloyd
May
This is a CellML description of Enciso and Sontag's 2004 mathematical model of testosterone dynamics.
Catherine Lloyd
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.
c.lloyd@auckland.ac.nz
Journal of Mathematical Biology
Catherine Lloyd
Eduardo
Sontag
D
2007-07-18T00:00:00+00:00
Enciso and Sontag's 2004 mathematical model of testosterone dynamics.
German
Enciso
The University of Auckland
Auckland Bioengineering Institute
Catherine
Lloyd
May
15480670
keyword