An Optimization Formulation for Characterization of Pulsatile Cortisol Secretion
Total Page:16
File Type:pdf, Size:1020Kb
An optimization formulation for characterization of pulsatile cortisol secretion The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation Faghih, Rose T. et al. "An optimization formulation for characterization of pulsatile cortisol secretion." Frontiers in Neuroscience 9 (August 2015): 228 © 2015 Faghih, Dahleh and Brown As Published http://dx.doi.org/10.3389/fnins.2015.00228 Publisher Frontiers Research Foundation Version Final published version Citable link http://hdl.handle.net/1721.1/112222 Terms of Use Creative Commons Attribution 4.0 International License Detailed Terms http://creativecommons.org/licenses/by/4.0/ ORIGINAL RESEARCH published: 11 August 2015 doi: 10.3389/fnins.2015.00228 An optimization formulation for characterization of pulsatile cortisol secretion Rose T. Faghih 1, 2, 3, 4*, Munther A. Dahleh 1, 4, 5, 6 and Emery N. Brown 2, 3, 7, 8 1 Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, USA, 2 Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, Cambridge, MA, USA, 3 Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital, Boston, MA, USA, 4 Laboratory for Information and Decision Systems, Massachusetts Institute of Technology, Cambridge, MA, USA, 5 Engineering Systems Division, Massachusetts Institute of Technology, Cambridge, MA, USA, 6 Institute for Data, Systems, and Society, Massachusetts Institute of Technology, Cambridge, MA, USA, 7 Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA, USA, 8 Department of Anesthesia, Harvard Medical School, Boston, MA, USA Edited by: Jason Ritt, Cortisol is released to relay information to cells to regulate metabolism and reaction Boston University, USA to stress and inflammation. In particular, cortisol is released in the form of pulsatile Reviewed by: signals. This low-energy method of signaling seems to be more efficient than continuous Alireza Mousavi, signaling. We hypothesize that there is a controller in the anterior pituitary that leads to Brunel University, UK Jose Luis Contreras-Vidal, pulsatile release of cortisol, and propose a mathematical formulation for such controller, University of Houston, USA which leads to impulse control as opposed to continuous control. We postulate that this *Correspondence: controller is minimizing the number of secretory events that result in cortisol secretion, Rose T. Faghih, Neuroscience Statistics Research which is a way of minimizing the energy required for cortisol secretion; this controller Laboratory, Massachusetts Institute of maintains the blood cortisol levels within a specific circadian range while complying Technology, 77 Massachusetts with the first order dynamics underlying cortisol secretion. We use an ℓ -norm cost Avenue, 46-6057 Cambridge, 0 MA, USA function for this controller, and solve a reweighed ℓ1-norm minimization algorithm for [email protected] obtaining the solution to this optimization problem. We use four examples to illustrate the performance of this approach: (i) a toy problem that achieves impulse control, (ii) two Specialty section: This article was submitted to examples that achieve physiologically plausible pulsatile cortisol release, (iii) an example Neuroprosthetics, where the number of pulses is not within the physiologically plausible range for healthy a section of the journal Frontiers in Neuroscience subjects while the cortisol levels are within the desired range. This novel approach results in impulse control where the impulses and the obtained blood cortisol levels Received: 29 March 2015 Accepted: 11 May 2015 have a circadian rhythm and an ultradian rhythm that are in agreement with the known Published: 11 August 2015 physiology of cortisol secretion. The proposed formulation is a first step in developing Citation: intermittent controllers for curing cortisol deficiency. This type of bio-inspired pulse Faghih RT, Dahleh MA and Brown EN (2015) An optimization formulation for controllers can be employed for designing non-continuous controllers in brain-machine characterization of pulsatile cortisol interface design for neuroscience applications. secretion. Front. Neurosci. 9:228. doi: 10.3389/fnins.2015.00228 Keywords: pulsatile control, cortisol secretion, endocrine control, mathematical modeling, circadian rhythm Frontiers in Neuroscience | www.frontiersin.org 1 August 2015 | Volume 9 | Article 228 Faghih et al. Characterization of pulsatile cortisol secretion 1. Introduction of ACTH did not activate cortisol secretion (Spiga et al., 2011). Moreover, studies on sheep in which the hypothalamus has Many hormones that have been well-investigated appear to been disconnected from the pituitary suggest that pulsatile input be released in pulses (Stavreva et al., 2009); for example, from hypothalamic secretagogues (e.g., CRH or vasopressin) is cortisol, gonadal steroids, and insulin are released in a pulsatile not necessary for the ultradian rhythm in cortisol secretion or manner (Veldhuis, 2008). Pulsatility is a physiological way for pulsatile cortisol secretion and pulsatile cortisol secretion is of increasing hormone concentrations rapidly and sending still maintained (Walker et al., 2010a). Hence, pulsatile cortisol distinct signaling information to target cells (Veldhuis, 2008). release is controlled by the dynamics in the anterior pituitary. Ultradian pulsatile hormone secretion allows for encoding Since pulsatile cortisol release seems to be more efficient than information via both amplitude and frequency modulation and is continuous signaling, it might be the case that the anterior a way of frequency encoding (Lightman and Conway-Campbell, pituitary is solving an optimal control problem. 2010; Walker et al., 2010b). Pulsatile signaling permits target We postulate that there is a controller in the anterior pituitary receptor recovery, rapid changes in hormone concentration, that controls the pulsatile secretion of cortisol and the ultradian and greater control, and is also more efficient than continuous rhythm of the pulses via the negative feedback effect of cortisol signaling (Walker et al., 2010b). The mechanism underlying on the anterior pituitary. Hence, by considering the known the generation of hormone pulses and why this method of physiology of the HPA axis, we shall formulate an optimization signaling is chosen by the body over continuous signaling problem that achieves impulse control. In optimal control theory, is not known. Since the transcriptional program prompted impulse control is a special case of bang-bang control, in which by hormone pulses is considerably different from constant an action leads to instantaneous changes in the states of the hormone treatment (Stavreva et al., 2009), it is crucial to system (Sethi and Thompson, 2000). Impulse control occurs understand the physiology underlying pulsatile hormone release. when there is not an upper bound on the control variable and Hormone pulsatility underlies multiple physiological processes. an infinite control is exerted on a state variable in order to cause For example, (i) cortisol oscillations have crucial effects on a finite jump (Sethi and Thompson, 2000). Minimizing an ℓ0- target cell gene expression and glucocorticoids receptor function norm cost function can achieve impulse control and we use a (McMaster et al., 2011; Walker et al., 2012). (ii) Some psychiatric reweighed ℓ1-norm formulation as a relaxation to the ℓ0-norm and metabolic diseases are associated with changes in cortisol to solve the proposed optimization formulation. Moreover, we pulsatility (Walker et al., 2010a). (iii) When the same amount of consider the first-order dynamics underlying cortisol synthesis corticosterone is administered by constant infusion rather than and the circadian amplitude constraints on the cortisol levels a pulsatile infusion, it results in a noticeably reduced ACTH when formulating the optimization problem. response to stress (Lightman and Conway-Campbell, 2010). In this study, we investigate pulsatile release of cortisol and propose 2. Methods a novel mathematical formulation that characterizes pulsatile cortisol secretion. We propose a physiologically plausible optimization problem Cortisol is released from the adrenal glands in pulses in for cortisol secretion by making the following assumptions: (1) response to pulsatile release of ACTH. CRH induces the release of Cortisol levels can be described by first-order kinetics for cortisol ACTH. In return, cortisol has a negative feedback effect on ACTH synthesis in the adrenal glands, cortisol infusion to the blood, and CRH release at the pituitary and hypothalamic levels. The and cortisol clearance by the liver described in Brown et al. timing and amplitudes of cortisol pulses vary throughout the day (2001), Faghih (2010), and Faghih et al. (2011, 2014). (2) There where the amplitude variations are due to the circadian rhythm is a time-varying cortisol demand [h(t)] that should be satisfied underlying cortisol release with periods of 12 and 24 h (Faghih throughout the day, which is a function of the circadian rhythm. et al., 2011), and the variations in the timing of cortisol pulses (3) There is a time-varying upper bound on the cortisol level result from the ultradian rhythm underlying cortisol release. [q(t)] that is a function of the upper bound on the cortisol level Between 15 and 22 secretory pulses of cortisol are expected over