Parathyroid hormone temporal effects on bone formation and resorption
Catherine
Lloyd
Auckland Bioengineering Institute, The University of Auckland
Model Status
This CellML model runs in OpenCell to recreate the published results. The units have been checked and they are consistent. PTH is pulsed every 6 hours and the model needs to be run for some time before the oscillations stabilise. The CellML model will also run in COR but due to the time unit being hours it is not really suitable for full length simulations.
Model Structure
ABSTRACT: Parathyroid hormone (PTH) paradoxically causes net bone loss (resorption) when administered in a continuous fashion, and net bone formation (deposition) when administered intermittently. Currently no pharmacological formulations are available to promote bone formation, as needed for the treatment of osteoporosis. The paradoxical behavior of PTH confuses endocrinologists, thus, a model bone resorption or deposition dependent on the timing of PTH administration would de-mystify this behavior and provide the basis for logical drug formulation. We developed a mathematical model that accounts for net bone loss with continuous PTH administration and net bone formation with intermittent PTH administration, based on the differential effects of PTH on the osteoblastic and osteoclastic populations of cells. Bone, being a major reservoir of body calcium, is under the hormonal control of PTH. The overall effect of PTH is to raise plasma levels of calcium, partly through bone resorption. Osteoclasts resorb bone and liberate calcium, but they lack receptors for PTH. The preosteoblastic precursors and preosteoblasts possess receptors for PTH, upon which the hormone induces differentiation from the precursor to preosteoblast and from the preosteoblast to the osteoblast. The osteoblasts generate IL-6; IL-6 stimulates preosteoclasts to differentiate into osteoclasts. We developed a mathematical model for the differentiation of osteoblastic and osteoclastic populations in bone, using a delay time of 1 hour for differentiation of preosteoblastic precursors into preosteoblasts and 2 hours for the differentiation of preosteoblasts into osteoblasts. The ratio of the number of osteoblasts to osteoclasts indicates the net effect of PTH on bone resorption and deposition; the timing of events producing the maximum ratio would induce net bone deposition. When PTH is pulsed with a frequency of every hour, the preosteoblastic population rises and decreases in nearly a symmetric pattern, with 3.9 peaks every 24 hours, and 4.0 peaks every 24 hours when PTH is administered every 6 hours. Thus, the preosteoblast and osteoblast frequency depends more on the nearly constant value of the PTH, rather than on the frequency of the PTH pulsations. Increasing the time delay gradually increases the mean value for the number of osteoblasts. The osteoblastic population oscillates for all intermittent administrations of PTH and even when the PTH infusion is constant. The maximum ratio of osteoblasts to osteoclasts occurs when PTH is administered in pulses of every 6 hours. The delay features in the model bear most of the responsibility for the occurrence of these oscillations, because without the delay and in the presence of constant PTH infusions, no oscillations occur. However, with a delay, under constant PTH infusions, the model generates oscillations. The osteoblast oscillations express limit cycle behavior. Phase plane analysis show simple and complex attractors. Subsequent to a disturbance in the number of osteoblasts, the osteoblasts quickly regain their oscillatory behavior and cycle back to the original attractor, typical of limit cycle behavior. Further, because the model was constructed with dissipative and nonlinear features, one would expect ensuing oscillations to show limit cycle behavior. The results from our model, increased bone deposition with intermittent PTH administration and increased bone resorption with constant PTH administration, conforms with experimental observations and with an accepted explanation for osteoporosis.
model diagram
Schematic diagram of the effect of PTH on the development of osteoblasts. PTH binds to receptors on the preosteoblast precursors and stimulates their transition to preosteoblasts. However, PTH binding to these preosteoblasts inhibits their differentiation into osteoblasts, and IL-6 (which is secreted by the osteoblasts) is believed to enhance this inhibitory effect. The osteoblasts then differentiate into osteocytes at a rate which is dependent on their number. The IL-6 produced by the osteoblasts also stimulates the differentiation of preosteoclasts to osteoclasts. Osteoclasts become senescent at a rate dependent on their number.
The original paper reference is cited below:
Parathyroid hormone temporal effects on bone formation and resorption, Martin H. Kroll, 2000, Bulletin of Mathematical Biology, 62, 163-187. PubMed ID: 10824426
Y
number of preosteoblasts
$\frac{d Y}{d \mathrm{time}}=\mathrm{k1}\mathrm{C1}\frac{P}{K+P}C-\mathrm{k2}\mathrm{C2}(1.0-\frac{P}{K+P})Y+\mathrm{ky}Y$
X
number of osteoblasts
$\frac{d X}{d \mathrm{time}}=\mathrm{k2}\mathrm{C2}(1.0-\frac{P}{K+P})Y-\mathrm{k3}X$
P
PTH
parathyroid hormone concentration
$\frac{d P}{d \mathrm{time}}=\begin{cases}10.0-\mathrm{k4}P & \text{if $(\mathrm{time}\ge 0.0)\land (\mathrm{time}< 6.0)$}\\ -(\mathrm{k4}P) & \text{if $(\mathrm{time}\ge 6.0)\land (\mathrm{time}< 12.0)$}\\ 10.0-\mathrm{k4}P & \text{if $(\mathrm{time}\ge 12.0)\land (\mathrm{time}< 18.0)$}\\ -(\mathrm{k4}P) & \text{if $(\mathrm{time}\ge 18.0)\land (\mathrm{time}< 24.0)$}\\ 10.0-\mathrm{k4}P & \text{if $(\mathrm{time}\ge 24.0)\land (\mathrm{time}< 30.0)$}\\ -(\mathrm{k4}P) & \text{if $(\mathrm{time}\ge 30.0)\land (\mathrm{time}< 36.0)$}\\ 10.0-\mathrm{k4}P & \text{if $(\mathrm{time}\ge 36.0)\land (\mathrm{time}< 42.0)$}\\ -(\mathrm{k4}P) & \text{if $(\mathrm{time}\ge 42.0)\land (\mathrm{time}< 48.0)$}\\ 10.0-\mathrm{k4}P & \text{if $(\mathrm{time}\ge 48.0)\land (\mathrm{time}< 54.0)$}\\ -(\mathrm{k4}P) & \text{if $(\mathrm{time}\ge 54.0)\land (\mathrm{time}< 60.0)$}\\ 10.0-\mathrm{k4}P & \text{if $(\mathrm{time}\ge 60.0)\land (\mathrm{time}< 66.0)$}\\ -(\mathrm{k4}P) & \text{if $(\mathrm{time}\ge 66.0)\land (\mathrm{time}< 72.0)$}\\ 10.0-\mathrm{k4}P & \text{if $(\mathrm{time}\ge 72.0)\land (\mathrm{time}< 78.0)$}\\ -(\mathrm{k4}P) & \text{if $(\mathrm{time}\ge 78.0)\land (\mathrm{time}< 84.0)$}\\ 10.0-\mathrm{k4}P & \text{if $(\mathrm{time}\ge 84.0)\land (\mathrm{time}< 90.0)$}\\ -(\mathrm{k4}P) & \text{if $(\mathrm{time}\ge 90.0)\land (\mathrm{time}< 96.0)$}\end{cases}$
Z
number of osteoclasts
$\frac{d Z}{d \mathrm{time}}=\mathrm{k5}\mathrm{C3}\frac{\mathrm{IL6}}{\mathrm{K2}+\mathrm{IL6}}-\mathrm{k6}Z$
IL6
interleukin-6 concentration
$\frac{d \mathrm{IL6}}{d \mathrm{time}}=0.1X-10.0\mathrm{IL6}$
parathyroid hormone
osteoporosis
osteoclast
osteoblast
endocrine
2007-07-24T00:00:00+00:00
c.lloyd@auckland.ac.nz
Parathyroid hormone temporal effects on bone formation and resorption
62
163
187
Catherine
Lloyd
May
The University of Auckland
The Bioengineering Institute
Bulletin of Mathematical Biology
This is a CellML description of Martin Kroll's 2000 mathematical model of the temporal effects of parathyroid hormone on bone formation and resorption.
Martin Kroll's 2000 mathematical model of the temporal effects of parathyroid hormone on bone formation and resorption.
Martin
Kroll
H
Catherine Lloyd
10824426
keyword
2000-01-00 00:00
The University of Auckland, Auckland Bioengineering Institute