TY - JOUR
T1 - Probabilistic recovery of neuroendocrine pulsatile, secretory and kinetic structure
T2 - An alternating discrete and continuous scheme
AU - Chattopadhyay, Somesh
AU - Keenan, Daniel M.
AU - Veldhuis, Johannes D.
PY - 2008
Y1 - 2008
N2 - The brain (hypothalamus) directs hormone secretion by the pituitary gland via burst-like (pulsatile) release of specific peptides at inferentially random times. These pulsatile signals supervise growth, reproduction, lactation, stress adaptations, water balance and immune responses. However, hypothalamic molecules are diluted > 3000-fold in systemic blood, leaving pituitary-hormone pulses as measurable surrogates. The latter (roughly) mirror hypothalamic peptide bursts on a 1:1 basis, albeit being observed in a noisy environment. As a window to the brain, one must accurately recover the pulse (onset) times, and thereby estimate hormone secretion and kinetic parameters (θ ∈ ⊖) without distortion. Based upon limited observed data, one would like to obtain probability statements about underlying pulsatility, secretion and kinetics. Moreover, to be applicable in today's clinical setting, it is important that any such procedure require minimal or no human input. We propose and justify the following method. First, the data (a pituitary hormone concentration time-profile) is "selectively smoothed" by a nonlinear diffusion equation, whose diffusion coefficient is inversely related to the degree of rapid increase. This procedure generates a collection of potential pulse time sets (T). Then, via an algorithm which alternates between a Metropolis algorithm on T and a time-homogeneous diffusion process on ⊖, a compact manifold with boundary, simulation from an appropriately formulated (posterior) probability measure is achieved. The method is applied to recover the underlying structure of brain-pituitary regulation in disease and aging.
AB - The brain (hypothalamus) directs hormone secretion by the pituitary gland via burst-like (pulsatile) release of specific peptides at inferentially random times. These pulsatile signals supervise growth, reproduction, lactation, stress adaptations, water balance and immune responses. However, hypothalamic molecules are diluted > 3000-fold in systemic blood, leaving pituitary-hormone pulses as measurable surrogates. The latter (roughly) mirror hypothalamic peptide bursts on a 1:1 basis, albeit being observed in a noisy environment. As a window to the brain, one must accurately recover the pulse (onset) times, and thereby estimate hormone secretion and kinetic parameters (θ ∈ ⊖) without distortion. Based upon limited observed data, one would like to obtain probability statements about underlying pulsatility, secretion and kinetics. Moreover, to be applicable in today's clinical setting, it is important that any such procedure require minimal or no human input. We propose and justify the following method. First, the data (a pituitary hormone concentration time-profile) is "selectively smoothed" by a nonlinear diffusion equation, whose diffusion coefficient is inversely related to the degree of rapid increase. This procedure generates a collection of potential pulse time sets (T). Then, via an algorithm which alternates between a Metropolis algorithm on T and a time-homogeneous diffusion process on ⊖, a compact manifold with boundary, simulation from an appropriately formulated (posterior) probability measure is achieved. The method is applied to recover the underlying structure of brain-pituitary regulation in disease and aging.
KW - Estimation
KW - Hormonal secretion
KW - Pulse detection
KW - Simulation by diffusion
UR - http://www.scopus.com/inward/record.url?scp=55449095322&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=55449095322&partnerID=8YFLogxK
U2 - 10.1090/S0033-569X-08-01024-4
DO - 10.1090/S0033-569X-08-01024-4
M3 - Article
AN - SCOPUS:55449095322
SN - 0033-569X
VL - 66
SP - 401
EP - 421
JO - Quarterly of Applied Mathematics
JF - Quarterly of Applied Mathematics
IS - 3
ER -