Model-based therapeutic correction of hypothalamic-pituitary-adrenal axis dysfunction

Advancing Research and Treatment: An Overview of Clinical Trials in Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) and Future Perspectives

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a chronic, debilitating, and multi-faceted illness. Heterogenous onset and clinical presentation with additional comorbidities make it difficult to diagnose, characterize, and successfully treat. Current treatment guidelines focus on symptom management, but with no clear target or causative mechanism, remission rates are low, and fewer than 5% of patients return to their pre-morbid activity levels. Therefore, there is an urgent need to undertake robust clinical trials to identify effective treatments. This review synthesizes insights from clinical trials exploring pharmacological interventions and dietary supplements targeting immunological, metabolic, gastrointestinal, neurological, and neuroendocrine dysfunction in ME/CFS patients which require further exploration. Additionally, the trialling of alternative interventions in ME/CFS based on reported efficacy in the treatment of illnesses with overlapping symptomology is also discussed. Finally, we provide important considerations and make recommendations, focusing on outcome measures, to ensure the execution of future high-quality clinical trials to establish clinical efficacy of evidence-based interventions that are needed for adoption in clinical practice.

Identification of an Unmet Medical Need: Height of Depression, Hypersomnia, and Sleep Apnea Positively Correlate With the Level of Fatigue in Patients With Immune Thrombocytopenia

INTRODUCTION

Immune thrombocytopenia (ITP) is a rare chronic disease, frequently accompanied by fatigue, which is an important comorbidity associated with this disease. Patients experience difficulties in managing their daily activities and a reduction in their overall quality of life (QoL). The causes of fatigue in ITP are not clarified yet, and underlying causes seem to be multifactorial. The development of fatigue may not solely be influenced by a decrease in platelet count but also by unknown factors as well as psychological reasons.

METHODS

This prospective, multicenter, exploratory, pilot study aimed to investigate which parameters contribute to the occurrence of fatigue in patients with ITP. Adult patients with ITP and with or without fatigue who visited the study center for their regular appointments were asked to complete questionnaires pertaining to patient-reported outcome measures regarding bleeding symptoms, depression, sleep apnea, and hypersomnia. Blood tests included platelet count as well as different parameters like vitamin D.

RESULTS

A total of 36 patients (100%; 27 females (75%) and nine males (25%)) with primary ITP, with a median age of 46.5 years (range 19‑83 years) were analyzed. The median duration of ITP was 4.5 years (min‑max 0-21). Approximately one-third of patients (29.4%; 10/34 patients) had no comorbidities. The two most frequently used current treatment options were "watch-and-wait" (38.9%; 14/36 patients) and "avatrombopag" (30.6%; 11/36 patients); eight patients (22.2%; 8/36 patients) needed rescue therapy with corticosteroids. There was a statistically negative correlation between fatigue and year of diagnosis (r=-0.41, p=0.014). Results indicated no statistically significant relationship between fatigue and age or differences in fatigue between the genders. Ferritin predicted fatigue with statistical significance. Platelet count was not correlated with the level of fatigue. A significant correlation was obvious between fatigue, depression, and obstructive sleep apnea syndrome (OSAS) as well as sleep-related problems (p<0.01).

DISCUSSION

Patient characteristics were comparable to that of other studies. The level of fatigue negatively impacts the lives of patients with ITP. Age and gender were not correlated with fatigue in ITP, which is in line with other reports. Interestingly, the fatigue level was higher in patients presenting with additional depression and poor sleeping quality due to, e.g., hypersomnia, which seems common. Fatigue levels seem independent from thrombocyte levels, which were reported elsewhere.

CONCLUSION

Patients diagnosed with ITP several years ago cope with their condition better than patients with a more recent diagnosis, who have higher levels of fatigue. Concurrent depression, hypersomnia, and sleep apnea are important underestimated factors, which do have a negative effect on the QoL of patients with ITP. We were able to show that patients with ITP might face an unmet medical need in terms of delayed diagnosis and supportive therapy. To our knowledge, this is the first report on combined findings of depression, hypersomnia, and sleep apnea in patients with ITP.

Targeting the Arginine Vasopressin V1b Receptor System and Stress Response in Depression and Other Neuropsychiatric Disorders

A healthy stress response is critical for good mental and overall health and promotes neuronal growth and adaptation, but the intricately balanced biological mechanisms that facilitate a stress response can also result in predisposition to disease when that equilibrium is disrupted. The hypothalamic-pituitary-adrenal (HPA) axis neuroendocrine system plays a critical role in the body's response and adaptation to stress, and vasopressinergic regulation of the HPA axis is critical to maintaining system responsiveness during chronic stress. However, exposure to repeated or excessive physical or emotional stress or trauma can shift the body's stress response equilibrium to a "new normal" underpinned by enduring changes in HPA axis function. Exposure to early life stress due to adverse childhood experiences can also lead to lasting neurobiological changes, including in HPA axis function. HPA axis impairment in patients with depression is considered among the most reliable findings in biological psychiatry, and chronic stress has been shown to play a major role in the pathogenesis and onset of depression and other neuropsychiatric disorders. Modulating HPA axis activity, for example via targeted antagonism of the vasopressin V_1b receptor, is a promising approach for patients with depression and other neuropsychiatric disorders associated with HPA axis impairment. Despite favorable preclinical indications in animal models, demonstration of clinical efficacy for the treatment of depressive disorders by targeting HPA axis dysfunction has been challenging, possibly due to the heterogeneity and syndromal nature of depressive disorders. Measures of HPA axis function, such as elevated cortisol levels, may be useful biomarkers for identifying patients who may benefit from treatments that modulate HPA axis activity. Utilizing clinical biomarkers to identify subsets of patients with impaired HPA axis function who may benefit is a promising next step in fine-tuning HPA axis activity via targeted antagonism of the V_1b receptor.

Perspective: Drawing on Findings From Critical Illness to Explain Myalgic Encephalomyelitis/Chronic Fatigue Syndrome

We propose an initial explanation for how myalgic encephalomyelitis / chronic fatigue syndrome (ME/CFS) could originate and perpetuate by drawing on findings from critical illness research. Specifically, we combine emerging findings regarding (a) hypoperfusion and endotheliopathy, and (b) intestinal injury in these illnesses with our previously published hypothesis about the role of (c) pituitary suppression, and (d) low thyroid hormone function associated with redox imbalance in ME/CFS. Moreover, we describe interlinkages between these pathophysiological mechanisms as well as “vicious cycles” involving cytokines and inflammation that may contribute to explain the chronic nature of these illnesses. This paper summarizes and expands on our previous publications about the relevance of findings from critical illness for ME/CFS. New knowledge on diagnostics, prognostics and treatment strategies could be gained through active collaboration between critical illness and ME/CFS researchers, which could lead to improved outcomes for both conditions.

Innervated adrenomedullary microphysiological system to model nicotine and opioid exposure

Transition to extrauterine life results in a surge of catecholamines necessary for increased cardiovascular, respiratory, and metabolic activity. Mechanisms mediating adrenomedullary catecholamine release are poorly understood. Important mechanistic insight is provided by newborns delivered by cesarean section or subjected to prenatal nicotine or opioid exposure, demonstrating impaired release of adrenomedullary catecholamines. To investigate mechanisms regulating adrenomedullary innervation, we developed compartmentalized 3D microphysiological systems (MPS) by exploiting GelPins, capillary pressure barriers between cell-laden hydrogels. The MPS comprises discrete cultures of adrenal chromaffin cells and preganglionic sympathetic neurons within a contiguous bioengineered microtissue. Using this model, we demonstrate that adrenal chromaffin innervation plays a critical role in hypoxia-mediated catecholamine release. Opioids and nicotine were shown to affect adrenal chromaffin cell response to a reduced oxygen environment, but neurogenic control mechanisms remained intact. GelPin containing MPS represent an inexpensive and highly adaptable approach to study innervated organ systems and improve drug screening platforms.

Modeling Neuroimmune Interactions in Human Subjects and Animal Models to Predict Subtype-Specific Multidrug Treatments for Gulf War Illness

Gulf War Illness (GWI) is a persistent chronic neuroinflammatory illness exacerbated by external stressors and characterized by fatigue, musculoskeletal pain, cognitive, and neurological problems linked to underlying immunological dysfunction for which there is no known treatment. As the immune system and the brain communicate through several signaling pathways, including the hypothalamic–pituitary–adrenal (HPA) axis, it underlies many of the behavioral and physiological responses to stressors via blood-borne mediators, such as cytokines, chemokines, and hormones. Signaling by these molecules is mediated by the semipermeable blood–brain barrier (BBB) made up of a monocellular layer forming an integral part of the neuroimmune axis. BBB permeability can be altered and even diminished by both external factors (e.g., chemical agents) and internal conditions (e.g., acute or chronic stress, or cross-signaling from the hypothalamic–pituitary–gonadal (HPG) axis). Such a complex network of regulatory interactions that possess feed-forward and feedback connections can have multiple response dynamics that may include several stable homeostatic states beyond normal health. Here we compare immune and hormone measures in the blood of human clinical samples and mouse models of Gulf War Illness (GWI) subtyped by exposure to traumatic stress for subtyping this complex illness. We do this via constructing a detailed logic model of HPA–HPG–Immune regulatory behavior that also considers signaling pathways across the BBB to neuronal–glial interactions within the brain. We apply conditional interactions to model the effects of changes in BBB permeability. Several stable states are identified in the system beyond typical health. Following alignment of the human and mouse blood profiles in the context of the model, mouse brain sample measures were used to infer the neuroinflammatory state in human GWI and perform treatment simulations using a genetic algorithm to optimize the Monte Carlo simulations of the putative treatment strategies aimed at returning the ill system back to health. We identify several ideal multi-intervention strategies and potential drug candidates that may be used to treat chronic neuroinflammation in GWI.

Theory: Treatments for Prolonged ICU Patients May Provide New Therapeutic Avenues for Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS)

We here provide an overview of treatment trials for prolonged intensive care unit (ICU) patients and theorize about their relevance for potential treatment of myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS). Specifically, these treatment trials generally target: (a) the correction of suppressed endocrine axes, notably through a "reactivation" of the pituitary gland's pulsatile secretion of tropic hormones, or (b) the interruption of the "vicious circle" between inflammation, oxidative and nitrosative stress (O&NS), and low thyroid hormone function. There are significant parallels in the treatment trials for prolonged critical illness and ME/CFS; this is consistent with the hypothesis of an overlap in the mechanisms that prevent recovery in both conditions. Early successes in the simultaneous reactivation of pulsatile pituitary secretions in ICU patients-and the resulting positive metabolic effects-could indicate an avenue for treating ME/CFS. The therapeutic effects of thyroid hormones-including in mitigating O&NS and inflammation and in stimulating the adreno-cortical axis-also merit further studies. Collaborative research projects should further investigate the lessons from treatment trials for prolonged critical illness for solving ME/CFS.

Hypothesis: Mechanisms That Prevent Recovery in Prolonged ICU Patients Also Underlie Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS)

Here the hypothesis is advanced that maladaptive mechanisms that prevent recovery in some intensive care unit (ICU) patients may also underlie Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS). Specifically, these mechanisms are: (a) suppression of the pituitary gland's pulsatile secretion of tropic hormones, and (b) a "vicious circle" between inflammation, oxidative and nitrosative stress (O&NS), and low thyroid hormone function. This hypothesis should be investigated through collaborative research projects.

Timescales of Human Hair Cortisol Dynamics

Cortisol is a major human stress hormone, secreted within minutes of acute stress. Cortisol also has slower patterns of variation: a strong circadian rhythm and a seasonal rhythm. However, longitudinal cortisol dynamics in healthy individuals over timescales of months has rarely been studied. Here, we measured longitudinal cortisol in 55 healthy participants using 12 cm of hair, which provides a retrospective measurement over one year. Individuals showed (non-seasonal) fluctuations averaging about 22% around their baseline. Fourier analysis reveals dominant slow frequencies with periods of months to a year. These frequencies can be explained by a mathematical model of the hormonal cascade that controls cortisol, the HPA axis, when including the slow timescales of tissue turnover of the glands. Measuring these dynamics is important for understanding disorders in which cortisol secretion is impaired over months, such as mood disorders, and to test models of cortisol feedback control.

A new model for the HPA axis explains dysregulation of stress hormones on the timescale of weeks

Stress activates a complex network of hormones known as the hypothalamic-pituitary-adrenal (HPA) axis. The HPA axis is dysregulated in chronic stress and psychiatric disorders, but the origin of this dysregulation is unclear and cannot be explained by current HPA models. To address this, we developed a mathematical model for the HPA axis that incorporates changes in the total functional mass of the HPA hormone-secreting glands. The mass changes are caused by HPA hormones which act as growth factors for the glands in the axis. We find that the HPA axis shows the property of dynamical compensation, where gland masses adjust over weeks to buffer variation in physiological parameters. These mass changes explain the experimental findings on dysregulation of cortisol and ACTH dynamics in alcoholism, anorexia, and postpartum. Dysregulation occurs for a wide range of parameters and is exacerbated by impaired glucocorticoid receptor (GR) feedback, providing an explanation for the implication of GR in mood disorders. These findings suggest that gland-mass dynamics may play an important role in the pathophysiology of stress-related disorders.

Leveraging Prior Knowledge of Endocrine Immune Regulation in the Therapeutically Relevant Phenotyping of Women With Chronic Fatigue Syndrome

PURPOSE

The complex and varied presentation of myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) has made it difficult to diagnose, study, and treat. Its symptoms and likely etiology involve multiple components of endocrine and immune regulation, including the hypothalamic-pituitary-adrenal axis, the hypothalamic-pituitary-gonadal axis, and their interactive oversight of immune function. We propose that the persistence of ME/CFS may involve changes in the regulatory interactions across these physiological axes. We also propose that the robustness of this new pathogenic equilibrium may at least in part explain the limited success of conventional single-target therapies.

METHODS

A comprehensive model was constructed of female endocrine-immune signaling consisting of 28 markers linked by 214 documented regulatory interactions. This detailed model was then constrained to adhere to experimental measurements in a subset of 17 candidate immune markers measured in peripheral blood of patients with ME/CFS and healthy control subjects before, during, and after a maximal exercise challenge. A set of 26 competing numerical models satisfied these data to within 5% error.

FINDINGS

Mechanistically informed predictions of endocrine and immune markers that were either unmeasured or exhibited high subject-to-subject variability pointed to possible context-specific overexpression in ME/CFS at rest of corticotropin-releasing hormone, chemokine (C-X-C motif) ligand 8, estrogen, follicle-stimulating hormone (FSH), gonadotropin-releasing hormone 1, interleukin (IL)-23, and luteinizing hormone, and underexpression of adrenocorticotropic hormone, cortisol, interferon-γ, IL-10, IL-17, and IL-1α. Simulations of rintatolimod and rituximab treatment predicted a shift in the repertoire of available endocrine-immune regulatory regimens. Rintatolimod was predicted to make available substantial remission in a significant subset of subjects, in particular those with low levels of IL-1α, IL-17, and cortisol; intermediate levels of progesterone and FSH; and high estrogen levels. Rituximab treatment was predicted to support partial remission in a smaller subset of patients with ME/CFS, specifically those with low norepinephrine, IL-1α, chemokine (C-X-C motif) ligand 8, and cortisol levels; intermediate FSH and gonadotropin-releasing hormone 1 levels; and elevated expression of tumor necrosis factor-α, luteinizing hormone, IL-12, and B-cell activation.

IMPLICATIONS

Applying a rigorous filter of known signaling mechanisms to experimentally measured immune marker expression in ME/CFS has highlighted potential new context-specific markers of illness. These novel endocrine and immune markers may offer useful candidates in delineating new subtypes of ME/CFS and may inform on refinements to the inclusion criteria and instrumentation of new and ongoing trials involving rintatolimod and rituximab treatment protocols.

High-fidelity discrete modeling of the HPA axis: a study of regulatory plasticity in biology

BACKGROUND

The hypothalamic-pituitary-adrenal (HPA) axis is a central regulator of stress response and its dysfunction has been associated with a broad range of complex illnesses including Gulf War Illness (GWI) and Chronic Fatigue Syndrome (CFS). Though classical mathematical approaches have been used to model HPA function in isolation, its broad regulatory interactions with immune and central nervous function are such that the biological fidelity of simulations is undermined by the limited availability of reliable parameter estimates.

METHOD

Here we introduce and apply a generalized discrete formalism to recover multiple stable regulatory programs of the HPA axis using little more than connectivity between physiological components. This simple discrete model captures cyclic attractors such as the circadian rhythm by applying generic constraints to a minimal parameter set; this is distinct from Ordinary Differential Equation (ODE) models, which require broad and precise parameter sets. Parameter tuning is accomplished by decomposition of the overall regulatory network into isolated sub-networks that support cyclic attractors. Network behavior is simulated using a novel asynchronous updating scheme that enforces priority with memory within and between physiological compartments.

RESULTS

Consistent with much more complex conventional models of the HPA axis, this parsimonious framework supports two cyclic attractors, governed by higher and lower levels of cortisol respectively. Importantly, results suggest that stress may remodel the stability landscape of this system, favoring migration from one stable circadian cycle to the other. Access to each regime is dependent on HPA axis tone, captured here by the tunable parameters of the multi-valued logic. Likewise, an idealized glucocorticoid receptor blocker alters the regulatory topology such that maintenance of persistently low cortisol levels is rendered unstable, favoring a return to normal circadian oscillation in both cortisol and glucocorticoid receptor expression.

CONCLUSION

These results emphasize the significance of regulatory connectivity alone and how regulatory plasticity may be explored using simple discrete logic and minimal data compared to conventional methods.

Old muscle in young body: an aphorism describing the Chronic Fatigue Syndrome

The chronic fatigue syndrome (CFS) otherwise known as myalgic encephalomyelitis (ME), is a debilitating syndrome whose identification is very complex due to lack of precise diagnostic criteria. This pathology begins with limitations in duration and intensity of exercise and rapid onset of pain during physical activity. Its etiology is unknown, and symptoms are not limited to the muscles. Epidemiology is rather difficult to delimit, even if it affects mainly young (20-40 years), female subjects. The results of muscular research show some peculiarities that can justify what has been observed in vivo. In particular, 1. presence of oxidative damage of lipid component of biological membranes and DNA not compensated by the increase of the scavenger activity; 2. Excitation-Contraction (E-C) alteration with modification of Ca2+ transport; 3. passage from slow to fast fiber phenotype; 4. inability to increase glucose uptake; 5. presence of mitochondrial dysfunction; and 6. genes expressed differentially (particularly those involved in energy production). The skeletal muscles of CFS / ME patients show a significant alteration of the oxidative balance due to mitochondrial alteration and of the fiber phenotype composition as shown in sarcopenic muscles of the elderly. Vice versa, the muscle catabolism does not appear to be involved in the onset of this syndrome. The data support the hypothesis that patients with CFS are subjected to some of the problems typical for muscle aging, which is probably related to disorders of muscle protein synthesis and biogenesis of mitochondria. Patients with CFS can benefit from an appropriate training program because no evidence suggests that physical exercise worsens symptoms. Type, intensity and duration of any physical activity that activates muscle contraction (including Electrical Stimulation) require further investigation even if it is known that non-exhaustive physical activity decreases painful symptomatology.

Exploring the Diagnostic Potential of Immune Biomarker Co-expression in Gulf War Illness

Complex disorders like Gulf War illness (GWI) often defy diagnosis on the basis of a single biomarker and may only be distinguishable by considering the co-expression of multiple markers measured in response to a challenge. We demonstrate the practical application of such an approach using an example where blood was collected from 26 GWI, 13 healthy control subjects, and 9 unhealthy controls with chronic fatigue at three points during a graded exercise challenge. A 3-way multivariate projection model based on 12 markers of endocrine and immune function was constructed using a training set of n = 10 GWI and n = 11 healthy controls. These groups were separated almost completely on the basis of two co-expression patterns. In a separate test set these same features allowed for discrimination of new GWI subjects (n = 16) from unhealthy (n = 9) and healthy control subjects with a sensitivity of 70% and a specificity of 90%.

Corticosterone oscillations during mania induction in the lateral hypothalamic kindled rat-Experimental observations and mathematical modeling

Changes in the hypothalamic-pituitary-adrenal (HPA) axis activity constitute a key component of bipolar mania, but the extent and nature of these alterations are not fully understood. We use here the lateral hypothalamic-kindled (LHK) rat model to deliberately induce an acute manic-like episode and measure serum corticosterone concentrations to assess changes in HPA axis activity. A mathematical model is developed to succinctly describe the entwined biochemical transformations that underlay the HPA axis and emulate by numerical simulations the considerable increase in serum corticosterone concentration induced by LHK. Synergistic combination of the LHK rat model and dynamical systems theory allows us to quantitatively characterize changes in HPA axis activity under controlled induction of acute manic-like states and provides a framework to study in silico how the dynamic integration of neurochemical transformations underlying the HPA axis is disrupted in these states.

Onset, timing, and exposure therapy of stress disorders: mechanistic insight from a mathematical model of oscillating neuroendocrine dynamics

BACKGROUND

The hypothalamic-pituitary-adrenal (HPA) axis is a neuroendocrine system that regulates numerous physiological processes. Disruptions in the activity of the HPA axis are correlated with stress-related diseases such as post-traumatic stress disorder (PTSD) and major depressive disorder. In this paper, we characterize "normal" and "diseased" states of the HPA axis as basins of attraction of a dynamical system describing the inhibition of peptide hormones such as corticotropin-releasing hormone (CRH) and adrenocorticotropic hormone (ACTH) by circulating glucocorticoids such as cortisol (CORT).

RESULTS

In addition to including key physiological features such as ultradian oscillations in cortisol levels and self-upregulation of CRH neuron activity, our model distinguishes the relatively slow process of cortisol-mediated CRH biosynthesis from rapid trans-synaptic effects that regulate the CRH secretion process. We show that the slow component of the negative feedback allows external stress-induced reversible transitions between "normal" and "diseased" states in novel intensity-, duration-, and timing-dependent ways.

CONCLUSION

Our two-step negative feedback model suggests a mechanism whereby exposure therapy of stress disorders such as PTSD may act to normalize downstream dysregulation of the HPA axis. Our analysis provides a causative rationale for improving treatments and guiding the design of new protocols.

Correcting hypothalamic-pituitary-adrenal axis dysfunction using observer-based explicit nonlinear model predictive control

The hypothalamic-pituitary-adrenal (HPA) axis is critical in maintaining homeostasis under physical and psychological stress by modulating cortisol levels in the body. Dysregulation of cortisol levels is linked to numerous stress-related disorders. In this paper, an automated treatment methodology is proposed, employing a variant of nonlinear model predictive control (NMPC), called explicit MPC (EMPC). The controller is informed by an unknown input observer (UIO), which estimates various hormonal levels in the HPA axis system in conjunction with the magnitude of the stress applied on the body, based on measured concentrations of adreno-corticotropic hormones (ACTH). The proposed closed-loop control strategy is tested on multiple in silico patients and the effectiveness of the controller performance is demonstrated.

Modeling the hypothalamus-pituitary-adrenal axis: A review and extension

Multiple models of the hypothalamus-pituitary-adrenal (HPA) axis have been developed to characterize the oscillations seen in the hormone concentrations and to examine HPA axis dysfunction. We reviewed the existing models, then replicated and compared five of them by finding their correspondence to a dataset consisting of ACTH and cortisol concentrations of 17 healthy individuals. We found that existing models use different feedback mechanisms, vary in the level of details and complexities, and offer inconsistent conclusions. None of the models fit the validation dataset well. Therefore, we re-calibrated the best performing model using partial calibration and extended the model by adding individual fixed effects and an exogenous circadian function. Our estimated parameters reduced the mean absolute percent error significantly and offer a validated reference model that can be used in diverse applications. Our analysis suggests that the circadian and ultradian cycles are not created endogenously by the HPA axis feedbacks, which is consistent with the recent literature on the circadian clock and HPA axis.

Achieving Remission in Gulf War Illness: A Simulation-Based Approach to Treatment Design

Gulf War Illness (GWI) is a chronic multi-symptom disorder affecting up to one-third of the 700,000 returning veterans of the 1991 Persian Gulf War and for which there is no known cure. GWI symptoms span several of the body’s principal regulatory systems and include debilitating fatigue, severe musculoskeletal pain, cognitive and neurological problems. Using computational models, our group reported previously that GWI might be perpetuated at least in part by natural homeostatic regulation of the neuroendocrine-immune network. In this work, we attempt to harness these regulatory dynamics to identify treatment courses that might produce lasting remission. Towards this we apply a combinatorial optimization scheme to the Monte Carlo simulation of a discrete ternary logic model that represents combined hypothalamic-pituitary-adrenal (HPA), gonadal (HPG), and immune system regulation in males. In this work we found that no single intervention target allowed a robust return to normal homeostatic control. All combined interventions leading to a predicted remission involved an initial inhibition of Th1 inflammatory cytokines (Th1Cyt) followed by a subsequent inhibition of glucocorticoid receptor function (GR). These first two intervention events alone ended in stable and lasting return to the normal regulatory control in 40% of the simulated cases. Applying a second cycle of this combined treatment improved this predicted remission rate to 2 out of 3 simulated subjects (63%). These results suggest that in a complex illness such as GWI, a multi-tiered intervention strategy that formally accounts for regulatory dynamics may be required to reset neuroendocrine-immune homeostasis and support extended remission.

A role for homeostatic drive in the perpetuation of complex chronic illness: Gulf War Illness and chronic fatigue syndrome

A key component in the body's stress response, the hypothalamic-pituitary-adrenal (HPA) axis orchestrates changes across a broad range of major biological systems. Its dysfunction has been associated with numerous chronic diseases including Gulf War Illness (GWI) and chronic fatigue syndrome (CFS). Though tightly coupled with other components of endocrine and immune function, few models of HPA function account for these interactions. Here we extend conventional models of HPA function by including feed-forward and feedback interaction with sex hormone regulation and immune response. We use this multi-axis model to explore the role of homeostatic regulation in perpetuating chronic conditions, specifically GWI and CFS. An important obstacle in building these models across regulatory systems remains the scarcity of detailed human in vivo kinetic data as its collection can present significant health risks to subjects. We circumvented this using a discrete logic representation based solely on literature of physiological and biochemical connectivity to provide a qualitative description of system behavior. This connectivity model linked molecular variables across the HPA axis, hypothalamic-pituitary-gonadal (HPG) axis in men and women, as well as a simple immune network. Inclusion of these interactions produced multiple alternate homeostatic states and sexually dimorphic responses. Experimental data for endocrine-immune markers measured in male GWI subjects showed the greatest alignment with predictions of a naturally occurring alternate steady state presenting with hypercortisolism, low testosterone and a shift towards a Th1 immune response. In female CFS subjects, expression of these markers aligned with an alternate homeostatic state displaying hypocortisolism, high estradiol, and a shift towards an anti-inflammatory Th2 activation. These results support a role for homeostatic drive in perpetuating dysfunctional cortisol levels through persistent interaction with the immune system and HPG axis. Though coarse, these models may nonetheless support the design of robust treatments that might exploit these regulatory regimes.

Differentiation between genomic and non-genomic feedback controls yields an HPA axis model featuring hypercortisolism as an irreversible bistable switch

Background

The hypothalamic-pituitary-adrenal axis (HPA axis) is a major part of the neuroendocrine system responsible for the regulation of the response to physical or mental stress and for the control of the synthesis of the stress hormone cortisol. Dysfunctions of the HPA axis characterized by either low (hypocortisolism) or increased (hypercortisolism) cortisol levels are implicated in various pathological conditions. Their understanding and therapeutic correction may be supported by mathematical modeling and simulation of the HPA axis.

Methods

Mass action and Michaelis Menten enzyme kinetics were used to provide a mechanistic description of the feedback mechanisms within the pituitary gland cells by which cortisol inhibits its own production. A separation of the nucleus from the cytoplasm by compartments enabled a differentiation between slow genomic and fast non-genomic processes. The model in parts was trained against time resolved ACTH stress response data from an in vitro cell culture of murine AtT-20 pituitary tumor cells and analyzed by bifurcation discovery tools.

Results

A recently found pituitary gland cell membrane receptor that mediates rapid non-genomic actions of glucocorticoids has been incorporated into our model of the HPA axis. As a consequence of the distinction between genomic and non-genomic feedback processes our model possesses an extended dynamic repertoire in comparison to existing HPA models. In particular, our model exhibits limit cycle oscillations and bistable behavior associated to hypocortisolism but also features a (second) bistable switch which captures irreversible transitions in hypercortisolism to elevated cortisol levels.

Conclusions

Model predictive control and inverse bifurcation analysis have been previously applied in the simulation-based design of therapeutic strategies for the correction of hypocortisolism. Given the HPA model extension presented in this paper, these techniques may also be used in the study of hypercortisolism. As an example, we show how sparsity enforcing penalization may suggest network interventions that allow the return from elevated cortisol levels back to nominal ones.

A review of hypothalamic-pituitary-adrenal axis function in chronic fatigue syndrome

Hypothalamic-pituitary-adrenal (HPA) axis dysfunction has been found in a high proportion of chronic fatigue syndrome (CFS) patients and includes enhanced corticosteroid-induced negative feedback, basal hypocortisolism, attenuated diurnal variation, and a reduced responsivity to challenge. A putative causal role for genetic profile, childhood trauma, and oxidative stress has been considered. In addition, the impact of gender is demonstrated by the increased frequency of HPA axis dysregulation in females. Despite the temporal relationship, it is not yet established whether the endocrine dysregulation is causal, consequent, or an epiphenomenon of the disorder. Nonetheless, given the interindividual variation in the effectiveness of existing biological and psychological treatments, the need for novel treatment strategies such as those which target the HPA axis is clear.

Transportation decreases the pulse frequency of growth hormone in the blood of prepubertal male calves

Both the mean concentration and the pulse pattern of growth hormone (GH) in the blood are important for the metabolism and body growth of calves. Transportation is reported to decrease blood GH concentrations in prepubertal male calves. However, the effect of transportation on GH pulsatility remains unknown. Because transportation is important in moving these calves from calf-production farms to markets or fattening farms, we tested whether transportation decreases their GH pulse frequency. Five calves were subjected to transportation by trucking (transport group), while five were left in their shed (non-transport group). Both groups were subsequently subjected to frequent blood sampling at 15-min intervals for 5 h. In the transport group, the cortisol concentrations increased in the first hour (P < 0.05) but significantly decreased thereafter (P < 0.05) to lower than those of the non-transport group. During the 5-hour study period, the transport group displayed a similar mean GH concentration relative to the non-transport group, but displayed a delayed first GH pulse, and a lower number of GH pulses than the non-transport group (P < 0.05). Hence, transportation is suggested to decrease GH pulse frequency under abnormal cortisol states, presumably suppressing metabolism and body growth in prepubertal male calves.

Biomarkers for chronic fatigue

Fatigue that persists for 6 months or more is termed chronic fatigue. Chronic fatigue (CF) in combination with a minimum of 4 of 8 symptoms and the absence of diseases that could explain these symptoms, constitute the case definition for chronic fatigue syndrome/myalgic encephalomyelitis (CFS/ME). Inflammation, immune system activation, autonomic dysfunction, impaired functioning in the hypothalamic-pituitary-adrenal axis, and neuroendocrine dysregulation have all been suggested as root causes of fatigue. The identification of objective markers consistently associated with CFS/ME is an important goal in relation to diagnosis and treatment, as the current case definitions are based entirely on physical signs and symptoms. This review is focused on the recent literature related to biomarkers for fatigue associated with CFS/ME and, for comparison, those associated with other diseases. These markers are distributed across several of the body's core regulatory systems. A complex construct of symptoms emerges from alterations and/or dysfunctions in the nervous, endocrine and immune systems. We propose that new insight will depend on our ability to develop and deploy an integrative profiling of CFS/ME pathogenesis at the molecular level. Until such a molecular signature is obtained efforts to develop effective treatments will continue to be severely limited.

Inverse problems from biomedicine: inference of putative disease mechanisms and robust therapeutic strategies

Many complex diseases that are difficult to treat cannot be mapped onto a single cause, but arise from the interplay of multiple contributing factors. In the study of such diseases, it is becoming apparent that therapeutic strategies targeting a single protein or metabolite are often not efficacious. Rather, a systems perspective describing the interaction of physiological components is needed. In this paper, we demonstrate via examples of disease models the kind of inverse problems that arise from the need to infer disease mechanisms and/or therapeutic strategies. We identify the challenges that arise, in particular the need to devise strategies that are robust against variable physiological states and parametric uncertainties.

Cortisol suppression and hearing thresholds in tinnitus after low-dose dexamethasone challenge

BACKGROUND

Tinnitus is a frequent, debilitating hearing disorder associated with severe emotional and psychological suffering. Although a link between stress and tinnitus has been widely recognized, the empirical evidence is scant. Our aims were to test for dysregulation of the stress-related hypothalamus-pituitary adrenal (HPA) axis in tinnitus and to examine ear sensitivity variations with cortisol manipulation.

METHODS

Twenty-one tinnitus participants and 21 controls comparable in age, education, and overall health status but without tinnitus underwent basal cortisol assessments on three non-consecutive days and took 0.5 mg of dexamethasone (DEX) at 23:00 on the first day. Cortisol levels were measured hourly the next morning. Detection and discomfort hearing thresholds were measured before and after dexamethasone suppression test.

RESULTS

Both groups displayed similar basal cortisol levels, but tinnitus participants showed stronger and longer-lasting cortisol suppression after DEX administration. Suppression was unrelated to hearing loss. Discomfort threshold was lower after cortisol suppression in tinnitus ears.

CONCLUSIONS

Our findings suggest heightened glucocorticoid sensitivity in tinnitus in terms of an abnormally strong glucocorticoid receptor (GR)-mediated HPA-axis feedback (despite a normal mineralocorticoid receptor (MR)-mediated tone) and lower tolerance for sound loudness with suppressed cortisol levels. Long-term stress exposure and its deleterious effects therefore constitute an important predisposing factor for, or a significant pathological consequence of, this debilitating hearing disorder.

Modeling cortisol dynamics in the neuro-endocrine axis distinguishes normal, depression, and post-traumatic stress disorder (PTSD) in humans

Cortisol, secreted in the adrenal cortex in response to stress, is an informative biomarker that distinguishes anxiety disorders such as major depression and post-traumatic stress disorder (PTSD) from normal subjects. Yehuda et al. proposed a hypothesis that, in humans, the hypersensitive hypothalamus-pituitary-adrenal (HPA) axis is responsible for the occurrence of differing levels of cortisol in anxiety disorders. Specifically, PTSD subjects have lower cortisol levels during the late subjective night in comparison to normal subjects, and this was assumed to occur due to strong negative feedback loops in the HPA axis. In the present work, to address this hypothesis, we modeled the cortisol dynamics using nonlinear ordinary differential equations and estimated the kinetic parameters of the model to fit the experimental data of three categories, namely, normal, depressed, and PTSD human subjects. We concatenated the subjects (n = 3) in each category and created a model subject (n = 1) without considering the patient-to-patient variability in each case. The parameters of the model for the three categories were simultaneously obtained through global optimization. Bifurcation analysis carried out with the optimized parameters exhibited two supercritical Hopf points and, for the choice of parameters, the oscillations were found to be circadian in nature. The fitted kinetic parameters indicate that PTSD subjects have a strong negative feedback loop and, as a result, the predicted oscillating cortisol levels are extremely low at the nadir in contrast to normal subjects, albeit within the endocrinologic range. We also simulated the phenotypes for each of the categories and, as observed in the clinical data of PTSD patients, the simulated cortisol levels are consistently low at the nadir, and correspondingly the negative feedback was found to be extremely strong. These results from the model support the hypothesis that high stress intensity and strong negative feedback loop may cause hypersensitive neuro-endocrine axis that results in hypocortisolemia in PTSD.

PTSD is an anxiety disorder that occurs among persons exposed to a traumatic event involving life threat and injury. This is a co-morbid psychiatric disorder that occurs along with depression. Cortisol is an informative endocrine biomarker that can distinguish PTSD from other co-morbid disorders. In comparison to normal subjects, hypocortisolemia was observed during the night in PTSD, while hypercortisolemia was observed in depressed subjects. From analyzing the clinical data, Yehuda et al. hypothesized that hypocortisolemia in PTSD was due to the strong negative feedback loop operating in the neuroendocrine axis under severe stress. We complemented this hypothesis by constructing a mathematical model for cortisol dynamics in HPA axis and estimated the kinetic parameters that fitted the cortisol time series obtained from the clinical data of normal, depressed and PTSD patients. The parameters obtained from the simulated phenotypes also strongly support the hypothesis that, due to disruptive negative feedback loops, cortisol levels are different in normal, PTSD and depressed subjects during the night. Importantly, the model predicted the transitions from normal to various diseased states, and these transitions were shown to occur due to changes in the strength of the negative feedback loop and the stress intensity in the neuro-endocrine axis.

Exploring the diagnostic potential of immune biomarker coexpression in Gulf War Illness

Complex disorders like Gulf War Illness (GWI) often defy diagnosis on the basis of a single biomarker and may only be distinguishable by considering the coexpression of multiple markers measured in response to a challenge. We demonstrate the practical application of such an approach using an example where blood was collected from 26 GWI, 13 healthy control subjects, and 9 unhealthy controls with Chronic Fatigue at three points during a graded exercise challenge. A 3-way multivariate projection model based on 12 markers of endocrine and immune function was constructed using a training set of n = 10 GWI and n = 11 healthy controls. These groups were separated almost completely on the basis of two coexpression patterns. In a separate test set these same features allowed for discrimination of new GWI subjects (n = 16) from unhealthy (n = 9) and healthy control subjects with a sensitivity of 70% and a specificity of 90%.

Plasma neuropeptide Y: a biomarker for symptom severity in chronic fatigue syndrome

BACKGROUND

Chronic fatigue syndrome (CFS) is a complex, multi-symptom illness with a multisystem pathogenesis involving alterations in the nervous, endocrine and immune systems.Abnormalities in stress responses have been identified as potential triggers or mediators of CFS symptoms. This study focused on the stress mediator neuropeptide Y (NPY). We hypothesized that NPY would be a useful biomarker for CFS.

METHODS

The CFS patients (n = 93) were from the Chronic Fatigue and Related Disorders Clinic at the University of Miami and met the 1994 case definition of Fukuda and colleagues. Healthy sedentary controls (n = 100)) were from NIH or VA funded studies. Another fatiguing, multi-symptom illness, Gulf War Illness (GWI), was also compared to CFS. We measured NPY in plasma using a radioimmunoassay (RIA). Psychometric measures, available for a subset of CFS patients included: Perceived Stress Scale, Profile of Mood States, ATQ Positive & Negative Self-Talk Scores, the COPE, the Beck Depression Inventory, Fatigue Symptom Inventory, Cognitive Capacity Screening Examination, Medical Outcomes Survey Short Form-36, and the Quality of Life Scale.

RESULTS

Plasma NPY was elevated in CFS subjects, compared to controls (p = .000) and to GWI cases (p = .000). Receiver operating characteristics (ROC) curve analyses indicated that the predictive ability of plasma NPY to distinguish CFS patients from healthy controls and from GWI was significantly better than chance alone. In 42 patients with CFS, plasma NPY had significant correlations (<0.05) with perceived stress, depression, anger/hostility, confusion, negative thoughts, positive thoughts, general health, and cognitive status. In each case the correlation (+ or -) was in the anticipated direction.

CONCLUSIONS

This study is the first in the CFS literature to report that plasma NPY is elevated compared to healthy controls and to a fatigued comparison group, GWI patients. The significant correlations of NPY with stress, negative mood, general health, depression and cognitive function strongly suggest that this peptide be considered as a biomarker to distinguish subsets of CFS.

The minimal model of the hypothalamic-pituitary-adrenal axis

This paper concerns ODE modeling of the hypothalamic-pituitary- adrenal axis (HPA axis) using an analytical and numerical approach, combined with biological knowledge regarding physiological mechanisms and parameters. The three hormones, CRH, ACTH, and cortisol, which interact in the HPA axis are modeled as a system of three coupled, nonlinear differential equations. Experimental data shows the circadian as well as the ultradian rhythm. This paper focuses on the ultradian rhythm. The ultradian rhythm can mathematically be explained by oscillating solutions. Oscillating solutions to an ODE emerges from an unstable fixed point with complex eigenvalues with a positive real parts and a non-zero imaginary parts. The first part of the paper describes the general considerations to be obeyed for a mathematical model of the HPA axis. In this paper we only include the most widely accepted mechanisms that influence the dynamics of the HPA axis, i.e. a negative feedback from cortisol on CRH and ACTH. Therefore we term our model the minimal model. The minimal model, encompasses a wide class of different realizations, obeying only a few physiologically reasonable demands. The results include the existence of a trapping region guaranteeing that concentrations do not become negative or tend to infinity. Furthermore, this treatment guarantees the existence of a unique fixed point. A change in local stability of the fixed point, from stable to unstable, implies a Hopf bifurcation; thereby, oscillating solutions may emerge from the model. Sufficient criteria for local stability of the fixed point, and an easily applicable sufficient criteria guaranteeing global stability of the fixed point, is formulated. If the latter is fulfilled, ultradian rhythm is an impossible outcome of the minimal model and all realizations thereof. The second part of the paper concerns a specific realization of the minimal model in which feedback functions are built explicitly using receptor dynamics. Using physiologically reasonable parameter values, along with the results of the general case, it is demonstrated that un-physiological values of the parameters are needed in order to achieve local instability of the fixed point. Small changes in physiologically relevant parameters cause the system to be globally stable using the analytical criteria. All simulations show a globally stable fixed point, ruling out periodic solutions even when an investigation of the 'worst case parameters' is performed.

Kindling and Oxidative Stress as Contributors to Myalgic Encephalomyelitis/Chronic Fatigue Syndrome

Myalgic Encephalomyelitis/chronic fatigue syndrome (ME/CFS) is one of the more complex illnesses involving multiple systems within the body. Onset of ME/CFS frequently occurs quickly, and many patients report a prior exposure to a viral infection. This debilitating illness can affect the immune, neuroendocrine, autonomic, and neurologic systems. Abnormal biological findings among some patients have included aberrant ion transport and ion channel activity, cortisol deficiency, sympathetic nervous system hyperactivity, EEG spike waves, left ventricular dysfunction in the heart, low natural killer cell cytotoxicity, and a shift from Th1 to Th2 cytokines. We propose that the kindling and oxidative stress theories provide a heuristic template for better understanding the at times conflicting findings regarding the etiology and pathophysiology of this illness.