Major depressive disorder and bistability in an HPA-CNS toggle switch

Ben Ron Mizrachi; Avichai Tendler; Omer Karin; Tomer Milo; Dafna Haran; Avi Mayo; Uri Alon

doi:10.1371/journal.pcbi.1011645

Abstract

Major depressive disorder (MDD) is the most common psychiatric disorder. It has a complex and heterogeneous etiology. Most treatments take weeks to show effects and work well only for a fraction of the patients. Thus, new concepts are needed to understand MDD and its dynamics. One of the strong correlates of MDD is increased activity and dysregulation of the hypothalamic-pituitary-adrenal (HPA) axis which produces the stress hormone cortisol. Existing mathematical models of the HPA axis describe its operation on the scale of hours, and thus are unable to explore the dynamic on the scale of weeks that characterizes many aspects of MDD. Here, we propose a mathematical model of MDD on the scale of weeks, a timescale provided by the growth of the HPA hormone glands under control of HPA hormones. We add to this the mutual inhibition of the HPA axis and the hippocampus and other regions of the central nervous system (CNS) that forms a toggle switch. The model shows bistability between euthymic and depressed states, with a slow timescale of weeks in its dynamics. It explains why prolonged but not acute stress can trigger a self-sustaining depressive episode that persists even after the stress is removed. The model explains the weeks timescale for drugs to take effect, as well as the dysregulation of the HPA axis in MDD, based on gland mass changes. This understanding of MDD dynamics may help to guide strategies for treatment.

Author summary

Major depressive disorder (MDD) is a prevalent psychiatric condition whose mechanisms are not fully understood. Existing treatments often take weeks to be effective and only benefit a portion of patients, necessitating a fresh perspective on understanding MDD. Research has shown a strong association between MDD and the dysregulation of the human stress hormone pathway. However, previous mathematical models of the stress hormone pathway have only explored its dynamics over a few hours or days, not capturing the weeks-long timescale relevant to MDD. To address this, we present a new mathematical model that incorporates the growth of the stress hormone glands, and their mutually inhibitory interactions with the brain. This model reveals how the brain and the stress axis toggle between normal and depressed states, with slow dynamics over weeks. It also explains why prolonged stress—but not brief stress—can trigger a sustained depressive episode. The model further sheds light on why antidepressant drugs take weeks to take effect. By comprehending the dynamics of depression, this model may provide valuable insights for guiding treatment strategies in the future.

Citation: Ron Mizrachi B, Tendler A, Karin O, Milo T, Haran D, Mayo A, et al. (2023) Major depressive disorder and bistability in an HPA-CNS toggle switch. PLoS Comput Biol 19(12): e1011645. https://doi.org/10.1371/journal.pcbi.1011645

Editor: James R. Faeder, University of Pittsburgh, UNITED STATES

Received: May 18, 2023; Accepted: November 1, 2023; Published: December 6, 2023

Copyright: © 2023 Ron Mizrachi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Code files, animations, generated by Mathematica, of the phase portraits under changing u, a or b parameters can be found here: https://github.com/benron20/MDD-and-bistability-in-the-HPA-stress-axis.

Funding: This project was funded by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No 856487). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Major depressive disorder (MDD) is the most common psychiatric syndrome, affecting about 20% of the world population at some point in life [1]. There were 21 million recorded cases of depression in adults in 2020 in the USA, about 8% of the adult population [2]. MDD in adolescents and children has risen sharply in the past decade in many countries [3]. People with MDD suffer from depressed mood or sadness, anhedonia, and may suffer from sleeping and eating disorders. The etiology of MDD is complex and heterogeneous, and multiple genetic and environmental factors contribute to its incidence [4].

Treatment for MDD includes medication, psychotherapy, physical exercise, and electroconvulsive therapy (ECT). Treatment is usually effective for only a fraction of the patients [5,6] and often takes several weeks to show effects [7]. In the past five decades few novel therapeutics for MDD have emerged that go beyond the monoamine theory of depression first developed in the 1960s [8].

Among the endocrine, inflammatory, metabolic and neuronal factors that have been explored, the stress hormone cortisol is one of the most significant correlates of MDD [9]. Stress is associated with MDD, and cortisol is higher than controls in most studies of people with MDD [10]. High cortisol is not only a readout of MDD—it can also be causal for depression. This is seen in Cushing’s syndrome in which a tumor induces high levels of cortisol, resulting in depression in over 50% of patients [11–13]. Similarly, prolonged high-dose glucocorticoid treatment, such as the cortisol analogues dexamethasone (DEX) or prednisone taken for a few weeks or longer, is associated with depression in a dose-dependent manner [14].

Furthermore, the endocrine system that secretes cortisol, the hypothalamic-pituitary-adrenal (HPA) axis, is dysregulated in many patients with MDD [15]. Enlarged adrenal glands in suicide victims and stressed animals were one of the clues that led to characterization of the HPA axis [16]. HPA dysregulation in patients with MDD is evident in enlarged adrenal glands [17,18], blunted adrenocorticotropic hormone (ACTH) response in corticotropin-releasing-hormone (CRH) tests [19] and elevated ACTH response in DEX-CRH tests [20,21].

MDD also involves dysfunctions of the central nervous system (CNS) such as impaired memory and hippocampal atrophy [22–24]. The CNS and HPA have mutual interactions involved in MDD. The hippocampus and the prefrontal cortex have inhibitory effects on HPA activation [25–27]. Conversely, high levels of cortisol can have inhibitory effects on these CNS regions. Cortisol is sensed by glucocorticoid receptors (GR) in the hippocampus and prefrontal cortex, causing reduced neurogenesis, synaptic arbor reduction, neuronal death [28,29], and a long-term decrease in CNS volume [30]. This interplay between the CNS and HPA is believed to be mechanistically important for MDD [26].

There are several open questions concerning the origin and dynamics of MDD. How does prolonged but not acute stress trigger depression? Once triggered, how does depression persist even after the stressor is gone? What is the origin of the timescale of weeks or longer of each depressive episode as well as the timescale for drugs to take effect? What is the origin of the difference between individuals susceptible to MDD and those not susceptible, in which a given stress triggers depression in the former and not the latter? Understanding these questions can help to guide MDD research and treatment.

One way to explore these open questions is to use mathematical modeling [31,32]. Some of the previous mathematical models of MDD did not include specific physiological mechanisms. They show recurrence of depressive episodes [33]. More mechanistic models based on the HPA axis describe the dynamics of hormones and their receptors on the scale of hours [34–43]. Some show the interplay between the CNS and HPA [39], and present bistability dynamics as an explanation of stress vulnerability and MDD [36,41,42,44]. These HPA models do not discuss the timescale of weeks which is seen in many aspects of MDD.

To address the dynamics of MDD on the scale of weeks, we consider the interplay between the HPA and the CNS within a minimal mathematical model. The model supplies the relevant timescale of weeks by including the slow growth of HPA glands under control of HPA hormones, an approach recently developed in our group [45–47]. We add to this the mutual inhibition by the CNS to obtain a new model which describes the mechanism of depression as a bistable system with a slow dynamical timescale of weeks. The model can explain why prolonged but not acute stress triggers depression, why treatment can take weeks, and why the HPA axis is dysregulated in MDD. It also reveals a hysteresis property that makes recovery from depression more difficult than entry into it.

Results

Toggle switch model of the HPA—CNS system

We present a minimal mathematical model of the interaction between the HPA axis and the CNS processes that inhibit it (Fig 1A). The model has two variables: the functional mass of the cortisol-secreting cells in the adrenal cortex, A, and the inhibitory activity of the CNS, h. The inhibitory activity h is inspired by the hippocampus, which inhibits the hypothalamic inputs to the HPA axis [26]. Our approach is agnostic, however, to the exact nature of the CNS variable h, which could comprise in addition to the hippocampus, also other regions including the prefrontal cortex [27,48].

Download:

Fig 1. Mathematical model for the interactions between the HPA axis and its inhibitory CNS regions.

A) Schematic of HPA and CNS interactions in the model (icons taken from Biorender [52] according to their Academic License Terms [53]. B) The model can be summarized as a toggle switch between CNS inhibitory activity h and the adrenal mass, A. CNS inhibits the hypothalamic input and thus leads to a decrease in A. A releases cortisol, which inhibits CNS through the glucocorticoid receptor GR. C) Nullclines show three fixed points in susceptible individuals. Euthymic and depressed states are marked in black dots. D) In non-susceptible individuals there is only one euthymic fixed point. C and D used identical parameters except stress input u = 0.75 and u = 0.4 respectively. An analogous shift between 3 and 1 fixed points can occur when changing the CNS system parameters to represent genetic variants. Figure created with Biorender.

https://doi.org/10.1371/journal.pcbi.1011645.g001

The main idea is that cortisol produced by the HPA axis acts through the glucocorticoid receptors (GR) [30,49] to inhibit the inhibitory CNS activity h [26]. This activity, in turn, inhibits the HPA axis. The mutual inhibition creates a toggle-switch circuit: A and h inhibit each other (Fig 1B). Since the adrenal cortex functional mass A changes over weeks, a slow timescale is introduced to the dynamics of the system.

The model can be compactly understood by the nullcline method, an approach which separates the toggle switch into two arms [50,51]. We first treat each arm independently, and then combine them to find the fixed points of the system, given by the crossing points of the nullclines.

The conclusions are robust in the sense that they depend only on the qualitative shapes of the nullclines and not on their precise mathematical forms. This is useful for MDD in which many of the physiological mechanisms are not yet fully characterized. In addition to the qualitative approach, we describe in Methods a specific mathematical model based on a model for the HPA axis [45] together with equations for the CNS inhibitory activity h. The Methods section details the equations and parameters used to simulate the dynamics.

The first nullcline describes the inhibition of the HPA axis by the CNS. It is defined by the adrenal cortex mass steady state, , as a function of the inhibitory activity h. h inhibits the HPA axis, since a large h reduces CRH [26] the inducer of ACTH which is the growth factor for the adrenal cortex A. If one clamps h at a given value, A changes its functional mass over the timescale of weeks to reach a steady state A_st. The higher h, the lower the HPA activity and hence the smaller A_st. This provides a decreasing curve plotted in the phase plane whose axes are A and h (Fig 1C, red line).

In the model described in Methods, solving Eq (18) gives a specific form for this nullcline, so that A drops with h. Note for future reference that adrenal mass A rises with the stress input u.

The second arm of the toggle switch, in which A inhibits h, is given by steady state CNS activity h_st as a function of A. The steady state is the value of h at which . Due to the effects of cortisol, h declines with A. Since cortisol inhibits the activity of relevant brain regions such as the hippocampus and the prefrontal cortex through the GR receptor, which has a cooperative (step like) activity dependence on cortisol, we assume that h goes from a high level h_high to a low level h_low in a step-like manner (Fig 1C, black line) when cortisol produced by A crosses the dissociation constant k_d of the GR receptor.

The HPA-CNS model gives the steady state solution , which is a step-like function due to the high cooperativity of the GR receptor with a Hill coefficient of about 3.

The qualitative shapes of the nullclines are more general than the specific model equations and suffice for many of the conclusions below.

Model shows bistability between euthymic and depressed states

There are two main scenarios for the model nullclines, which correspond to individuals susceptible and resilient to MDD. At the points where the nullclines intersect, the system is at steady state, because both and .

In the first scenario, the two nullclines intersect at three fixed points (Fig 1C). One fixed point has low adrenal cortex mass, A, and high CNS inhibition, h, and can be considered as the euthymic or healthy state. A second fixed point has high adrenal mass A and low CNS inhibition h, and represents the depressed state, with high cortisol and an enlarged adrenal cortex. In between, there is an unstable fixed point. This system is bistable: depending on initial conditions, it reaches either the healthy or the MDD state. There is a basin of attraction for each state. The boundary between these basins is a curve called the separatrix which goes through the unstable fixed point.

In the second scenario, the parameters of the nullclines are such that there is only a single intersection point (Fig 1D). This is the healthy state with low A—the MDD state with high A no longer exists. Such individuals are resilient with respect to major depressive disorder. All initial conditions eventually flow to the healthy fixed point.

Susceptibility to MDD is partly due to genetics [11,14]. To conceptualize this susceptibility using the model we note that most MDD-related mutations are not in HPA-related genes (except for genes related to GR), and thus we may assume that the A nullcline is similar between individuals. On the other hand, the h nullcline may depend on many genetic variants that affect the CNS. Multi-gene variations might push this nullcline up or down, and contribute to the hereditary component of susceptibility or resilience to MDD [54]. We analyzed which parameters changes most sensitively change the number of fixed points (Table 1).

Download:

Table 1. HPA-CNS model nominal parameter values and sensitivities.

https://doi.org/10.1371/journal.pcbi.1011645.t001

Prolonged but not acute stress can trigger self-sustained depression

The model can help understand the differential effects of prolonged and acute stress. Stress can be modeled by an increase in hypothalamic input u. This shifts the A nullcline up and to the right—because at a given level of CNS inhibition, high stress leads to an enlarged adrenal relative to the low stress condition. It also shifts the h nullcline to the left—because higher stress is more likely to activate GR receptors in the CNS. In the bistable case, a large enough shift can abolish the healthy fixed point, leaving MDD as the only possible steady-state (Fig 2A1).

Download:

Fig 2. Dynamics of the model show how prolonged but not acute stress can lead to depression.

A1) Stress shifts the nullclines, destroying the euthymic fixed point. A2) Dynamics under brief stress show recovery. A3) Dynamics under prolonged stress reach the depressed state on the timescale of weeks and remain there even after the stress is removed. B1) exiting depression requires stress to be reduced below baseline due to a hysteresis effect. B2) brief relief is not enough to exit depression B3) recovery on the timescale of weeks upon a prolonged stress relief. C) bifurcation plot of adrenal steady state size versus stress input u shows hysteresis. The transition from euthymic state to depressed state happens for prolonged stress levels of u = u₁. To return from the depressed state to the euthymic state, one must decrease stress below u = u₂<u₁.

https://doi.org/10.1371/journal.pcbi.1011645.g002

If stress lasts less than a few weeks, A starts growing but does not reach the separatrix. After stress is over, the nullcline shifts back to produce a healthy fixed point again, and A returns to its healthy size within weeks. No long-term depression develops (Fig 2A2).

If stress however lasts longer than a few weeks, A can grow so much that it crosses the original (pre-stress) separatrix. Now even if stress is over, and the healthy fixed point is reinstated, A continues to flow to the MDD fixed point because it lies within its basin of attraction. Once it reaches the MDD fixed point, the system stays there even if stress returns to its previous low value. Therefore, the MDD fixed point is self-sustaining, with low inhibition that fuels high cortisol, which in turn reduces inhibition in a vicious cycle (Fig 2A3).

This explains why prolonged but not acute stress can cause MDD in susceptible individuals, and why MDD persists even after the triggering stressor is removed. Animations of these dynamics are provided in a link in Methods.

Stress must be reduced below baseline to overcome depression due to a hysteresis effect

According to the model, MDD can be relieved by sufficiently reducing stress input. If one reduces stress input u, the A nullcline shifts down and to the left, and the h nullcline shifts to the right. This can abolish the MDD fixed point, leaving the euthymic state as the only possibility (Fig 2B1).

As before, this shift to low stress must last long enough (weeks) for the dynamics to cross the separatrix and converge to the euthymic fixed point. Otherwise the system returns to the MDD fixed point once stress is heightened (Fig 2B2 and 2B3).

Notably, the stress input must be reduced below the previous normal level in order to restore the healthy fixed point. This effect is known as hysteresis, and is caused because the positive feedback in the MDD state requires a larger than normal reduction of stress to compensate [51].

This hysteresis can be seen in a bifurcation plot, showing the steady state of A versus stress u. Starting in the healthy state, stress must cross a threshold u₁ to enter the MDD state of high A. However, starting in the MDD state, stress must drop below a lower threshold u₂<u₁ in order to return to the healthy state of low A (Fig 2D).

Hysteresis may explain why in about 30% of Cushing patients with depression, depression does not abate after the tumor is treated and cortisol level decreases [11–13]. We predict these to be patients with pre-existing susceptibility to MDD, namely nullclines that cross at two stable fixed points. Cushing syndrome causes the transition to the depressed state, which remains occupied even after treatment.

Model explains dysregulation of endocrine tests in MDD

The dysregulation of the HPA axis in MDD is characterized by two main endocrine tests, the CRH and DEX-CRH tests. One explanation for the dysregulation involves impaired GR function [29]. Here we offer a different explanation, in which GR is normal and the dysregulation is caused by enlarged glands in MDD.

The CRH test involves injecting CRH and measuring the ACTH response. Patients with MDD often show blunted ACTH responses in the CRH test relative to controls [19].

The present model explains this effect. Cortisol inhibits hypothalamic activity through MR and GR receptors, whereas it inhibits the pituitary activity through GR receptors only [49,55]. As a result, the model predicts that at the MDD fixed point, both the pituitary and the adrenal are larger than their euthymic size, and that the adrenal cortex grows more than the pituitary at steady state.

To model the CRH test, we add external CRH to the model in a pulse that lasts 30 min, corresponding to the lifetime of synthetic CRH (normal and synthetic CRH half-lives are about 4 min and 43 min respectively [45])

Control individuals are modeled with baseline stress levels and normal steady state adrenal size A. The CRH pulse causes an ACTH and cortisol response in both control and depressed simulations (Fig 3A1 and 3A2). However, the response is blunted in an individual with MDD, modeled with HPA glands that correspond to the MDD fixed point, namely enlarged adrenal and pituitary and high cortisol. The blunting is due to the repressive effect of high endogenous cortisol on the pituitary ACTH secretion that is not compensated by sufficient pituitary gland size.

Download:

Fig 3. Model explains CRH and DEX/CRH test aberrations in MDD based on enlarged gland sizes.

A1) Simulation shows blunted ACTH levels in CRH test in depressed patients compared to controls. A2) MDD patients show blunted ACTH in the CRH test (100μg CRH) compared to controls, adapted from Bardeleben at el [19]. B1) Simulation shows elevated ACTH levels in the DEX/CRH test in depressed patients compared to controls. B2) MDD patients show elevated ACTH in the DEX/CRH test (1.5 mg DEX and 100μg CRH 15 hours later), adapted from Heuser at el [21].

https://doi.org/10.1371/journal.pcbi.1011645.g003

A similar analysis explains the elevated ACTH in DEX-CRH test in depressed patients [20,21]. In this test, the long-lived cortisol analogue dexamethasone (DEX) is given at night. DEX suppresses the HPA axis. The next day, ACTH is measured after a CRH injection. Individuals with MDD show an elevated ACTH response. The model explains this due to their enlarged pituitary corticotroph mass P. The presence of DEX masks the higher cortisol found in MDD patients, and thus ACTH measured in the test corresponds primarily to the pituitary secretion capacity. MDD individuals in the model have larger P, which causes more ACTH to be secreted per CRH unit than controls (Fig 3B1, 3B2).

We note that these explanations of the tests do not require assuming that the GR receptor pathway is different in people with MDD and in controls.

Model explains timescale of weeks for antidepressant action

We now consider the effects of changes in the h nullcline. These can presumably be caused by drugs or behavioral interventions that affect the CNS. A sufficient rise in h should abrogate the MDD fixed point (Fig 4A1) in susceptible individuals. Now the healthy point is the only possible steady state.

Download:

Fig 4. Dynamics of the model may explain the timescale of weeks for antidepressant action.

A1) treatments that enhance CNS inhibitory activity can abolish the depressed fixed point. A2) brief treatment is not effective. A3) prolonged treatment of weeks or more is effective even after treatment ends.

https://doi.org/10.1371/journal.pcbi.1011645.g004

The slow timescale of the HPA glands, of weeks, indicates that the intervention must last at least a few weeks, otherwise the separatrix will not be crossed, and the dynamics will return to the MDD point (Fig 4A1-4A3).

This may explain why SSRIs and other medications for MDD take a few weeks to show effects [56,57]. Interestingly, SSRIs also activate the HPA axis on short timescales [58]. The model suggests they can be effective in the long term because they boost hippocampal health on the slow timescale of weeks [59]; such effects might cause the shift in the h nullcline that produces a therapeutic effect.

Discussion

We present a simple mathematical model for depression based on a CNS-HPA toggle switch which controls euthymic and depressed states on the timescale of weeks. This slow timescale is due to the turnover rate of the HPA glands. The timescale of the model explains why acute stress usually doesn’t cause depression whereas prolonged stress can. It also explains why drug treatment takes weeks to affect patients. The model shows bistability, in which patients enter a self-sustaining depressed state that lasts even after stress is removed. This is due to hysteresis in which entry into the depressed state is easier than exit from this state. In addition, the model provides a mechanism for blunted CRH and elevated DEX-CRH tests in depression. This mechanism is based on HPA gland growth and differs from standard explanations based on GR receptor dysregulation.

One class of MDD drugs that seems to work without the delay of weeks is ketamine-based treatments, which often show effects within a day [60,61]. If one attempts to understand such drugs within the framework of the model, one may assume that they work in the CNS downstream of cortisol. They mitigate the negative effect of cortisol on behavior, cognition and mood, as indicated in experiments on mice models [61]. Another (nonexclusive) possibility is that ketamine affects the HPA axis directly, lowering cortisol for a given stressor, as found in male mice [62], and in a human case study in which DEX suppression tracked ketamine in a treatment-resistant patient within days [63].

The present model supplies a mechanistic explanation for the dysregulation of the HPA axis observed in MDD—namely blunted CRH and elevated DEX-CRH tests. This explanation differs from the standard explanation in which the GR receptor feedback is weaker in patients with MDD [29]. In the CRH test, ACTH is blunted because the high cortisol in MDD inhibits ACTH production in the pituitary. A subtle point is that this effect depends on the known distribution of cortisol receptors: GR in the pituitary and GR and MR in the hypothalamus. This causes the pituitary to grow in MDD but to a lesser extent than the growth of the adrenal. If both glands grew by the same factor, the CRH test would not be blunted because the large pituitary mass would compensate.

One upshot of this blunted response is that beta-endorphin secretion from the pituitary, which is thought to be proportional to ACTH secretion since both are produced from the same polypeptide POMC [64], is decreased in MDD. Beta endorphin is a euphoric and analgesic hormone. Its blunting may thus contribute to anhedonia, because a given HPA stimulus—a given CRH pulse—causes lower levels of beta endorphin. This blunted response may contribute to the dysphoria and pain sensitivity in MDD [65,66].

In the DEX-CRH test, the extrinsic DEX acts as a cortisol analogue and masks the high endogenous cortisol level in MDD. The test exposes the functional mass of the pituitary, which is enlarged in MDD in our model. As a result, cortisol is elevated in this test in MDD relative to controls.

The model also presents a picture of susceptibility to MDD. Individuals vary in terms of their molecular/physiological parameters. This is described by the nullclines of the model. Individuals in whom the nullclines intersect three times have the potential to be stuck in MDD. Individuals for whom the nullclines intersect only once don’t have the potential for MDD, and instead show a return to euthymic state even after prolonged stress. The model may, in principle, offer a way to determine susceptibility, but this would require experimental approaches to quantify the interaction strengths of the CNS and the HPA. One approach may be quantifying emotional reactivity in terms of HPA dynamics [67,68].

In summary, the HPA-CNS toggle switch model offers an understanding of aspects of MDD and their dynamics on the scale of weeks. The model is based on qualitative physiological facts; future research on the specific mechanisms can help make it more precise. The present study rationalizes susceptibility to MDD in terms of bistability of the HPA-CNS feedback loop, explains why prolonged but not acute stress can trigger MDD, why many interventions take weeks to work and how HPA dysregulation occurs in MDD. This understanding may provide concepts to guide the development of treatment strategies for MDD.

Material and methods

HPA-CNS toggle switch model

We model the HPA axis using the model of Karin et al [45]. We denote the input signal to the HPA system as u [69]. This input is inhibited by the CNS activity h [25–27], a novel feature of the present model. The concentrations of the three hormones CRH, ACTH, and cortisol are denoted x₁, x₂, x₃. Their fast time scale dynamics are described by (1–3): (1) (2) (3) Where P is the total pituitary corticotroph functional mass and A is the functional mass of the cortisol secreting cells in the adrenal cortex. Cortisol inhibits hypothalamic activity through MR and GR receptors, whereas it inhibits the pituitary activity through GR receptors only [33,38]. We assume that the effects of MR and GR are multiplicative and use the same MR and GR expressions as done in Karin et al [45]. MR inhibitory activity is non cooperative so that , whereas GR inhibitory activity is Hill like with cooperativity coefficient 3 [45,70], so that . We note that at physiological cortisol levels, MR receptors tend to be nearly saturated whereas GR receptors are hardly activated [39,71], namely K_MR<<x₃<<K_GR. In this case we can approximate the MR and GR expressions—MR(x₃)≈x₃, GR(x₃)≈1. This provides a simpler set of equations for hormone levels (4–6): (4) (5) (6)

Next we add the equations for the functional mass of the glands. The growth factor of P is x₁ [72,73] and the growth factor of A is x₂ [74,75], as detailed in Karin et al [45]. Their dynamic is on a slow time scale of weeks, as a result of the pituitary and adrenal cortex turnover rates (7–8), as given by Karin et al [45] (7) (8)

The main new feature of the present model is the variable h that represents CNS inhibitory activity. We posit a slow time scale of weeks for h in Eq (9) by setting its removal rate appropriately. We note that this choice does not affect the conclusions: a timescale of hours or days for h would leave the qualitative results unchanged, and only have mild effects on the dynamics. The relevant CNS regions, hippocampus and frontal cortex, express both high-affinity MR and low-affinity GR receptors for cortisol. The GR generally has inhibitory effects. As mentioned above, MR receptors are nearly saturated by cortisol in its physiological range [71]. The GR receptors have a cooperative response to cortisol which is approximately step-like. We therefore, for simplicity, model the inhibitory effect of cortisol on the CNS as a step function, defined by a+bΘ(x₃>T). Using a Hill-type function instead of a step function leaves the conclusions unchanged but makes the analytical solution much more difficult. When cortisol level is lower than the threshold T, only MR receptors are activated, causing CNS inhibition by a factor of a. But when cortisol level exceeds the threshold T, both MR and GR are activated, inhibiting CNS by a factor of a+b.

(9)

Our HPA-CNS model is the set of equations (Eqs 1,2,3,7,8,9) which describe the dynamics of hormones and glands in the HPA-CNS system.

In order to demonstrate the general features of the model and nullclines, we first use a simplified version (Eqs 4,5,6,7,8,9). When considering compensation mechanisms where the GR plays an important role, we use the full HPA-CNS model Eq (Eqs 1,2,3,7,8,9).

The nominal parameter values for the model are presented in Table 1. These parameters provide 3 fixed points, and provide hormone levels equal to 1. The fold-change of each parameter required to bifurcate to a single fixed point are also shown. Some parameters cause a bifurcation both when raised and lowered, as represented by two fold-change values in the table.

The equations have two timescales—Eqs (4–6) have a timescale of hours, whereas Eqs (7–9) have a timescale of weeks. Solving the fast equations at quasi-steady-state (qss) gives the fast time levels of the hormones x₁, x₂, x₃ (10–12): (10) (11) (12) where one assumes that gland sizes are constant on the hour timescale.

Solving Eqs 7 and 8 at steady state shows that the slow time steady state of x₁ and x₂ is quite simple, and independent on the input u and most other model parameters, due to the properties of integral feedback in these equations [45].

(13)

(14)

We normalize our equations , so the baseline level is now . Eqs 13 and 14 show that upon a step change in input u, the change in hormone levels x₁ and x₂ is transient and the pituitary and adrenal glands grow to compensate to return CRH and ACTH to their original baseline levels within weeks. This growth and compensation can be calculated by substituting in Eqs (5–6). From Eq (5) we find that the adrenal gland steady-state sizes rises linearly with steady-state cortisol A_st ∝ x₃. On the other hand, from Eq (6) we find that the pituitary’s steady state size is independent on cortisol level, . Below we will show that if we solve the original Eqs (1–3) and consider GR inhibitive effect on the pituitary, we find a compensation mechanism from both adrenal and pituitary glands.

Now, if we substitute in Eq (4), we find that the steady state of cortisol is

This makes biological sense by having cortisol proportional to its hypothalamic input signal u. To obtain the HPA-CNS model slow time steady state solution, we substitute fast time scale solutions (10–12) into the slow Eqs (7–9): (15) (16) (17)

In order to demonstrate this slow steady state solution space in a two-dimensional phase plane, we use the pituitary’s steady state size, , removing one of the three equations and allowing a 2D analysis. This allows us to solve and to find the slow time scale steady state of the system, and to understand the shape of h and A nullclines (18–19): (18) (19)

These two nullclines can generate a bistable system as seen in Fig 1.

Now, we explore the compensation mechanisms of the system, using the steady state solution of (2–3) and (7–8). This compensation can be estimated by substituting in Eqs (2–3) steady state solutions. We find that A_st ∝ x₃ whereas . Both the pituitary and adrenal cortex grow in order to compensate for high CRH and ACTH levels, but the adrenal cortex grows more than the pituitary.

Interventions and drugs that affect the CNS are modeled by increasing CNS inhibitory activity by a factor of D: (20)

Animation of phase portrait under parameter changes

Animations, generated by Mathematica, of the phase portraits under changing u, a or b parameters can be found here: https://github.com/benron20/MDD-and-bistability-in-the-HPA-stress-axis.

Supporting information

The Supporting Information, including all code files and animations can be found here:

https://github.com/benron20/MDD-and-bistability-in-the-HPA-stress-axis.

S1 File. Hysteresis.

https://doi.org/10.1371/journal.pcbi.1011645.s001

(PY)

S2 File. CRH test graph.

https://doi.org/10.1371/journal.pcbi.1011645.s002

(PY)

S3 File. DEX/CRH test graph.

https://doi.org/10.1371/journal.pcbi.1011645.s003

(PY)

S4 File. CRH test data.

https://doi.org/10.1371/journal.pcbi.1011645.s004

(PY)

S5 File. DEX/CRH test data.

https://doi.org/10.1371/journal.pcbi.1011645.s005

(PY)

S6 File. HPA phase portrait.

https://doi.org/10.1371/journal.pcbi.1011645.s006

(IPYNB)

S7 File. Dynamic of brief stress.

https://doi.org/10.1371/journal.pcbi.1011645.s007

(PY)

S8 File. Dynamic of prolonged stress.

https://doi.org/10.1371/journal.pcbi.1011645.s008

(PY)

S9 File. Dynamic of brief stress reduction.

https://doi.org/10.1371/journal.pcbi.1011645.s009

(PY)

S10 File. Dynamic of prolonged stress reduction.

https://doi.org/10.1371/journal.pcbi.1011645.s010

(PY)

S11 File. Dynamic of short treatment.

https://doi.org/10.1371/journal.pcbi.1011645.s011

(PY)

S12 File. Dynamic of prolonged treatment.

https://doi.org/10.1371/journal.pcbi.1011645.s012

(PY)

References

1. Gutiérrez-Rojas L, Porras-Segovia A, Dunne H, Andrade-González N, Cervilla JA. Prevalence and correlates of major depressive disorder: a systematic review. Braz J Psychiatry. 2020 Aug 3;42(6):657–72. pmid:32756809
- View Article
- PubMed/NCBI
- Google Scholar
2. Major Depression [Internet]. National Institute of Mental Health (NIMH). [cited 2023 Jan 29]. Available from: https://www.nimh.nih.gov/health/statistics/major-depression.
3. Thapar A, Eyre O, Patel V, Brent D. Depression in young people. The Lancet. 2022 Aug 20;400(10352):617–31.
- View Article
- Google Scholar
4. Wong ML, Licinio J. Research and treatment approaches to depression. Nat Rev Neurosci. 2001 May;2(5):343–51. pmid:11331918
- View Article
- PubMed/NCBI
- Google Scholar
5. Moncrieff J, Cooper RE, Stockmann T, Amendola S, Hengartner MP, Horowitz MA. The serotonin theory of depression: a systematic umbrella review of the evidence. Mol Psychiatry. 2022 Jul 20;1–14. pmid:35854107
- View Article
- PubMed/NCBI
- Google Scholar
6. Karrouri R, Hammani Z, Benjelloun R, Otheman Y. Major depressive disorder: Validated treatments and future challenges. World J Clin Cases. 2021 Nov 6;9(31):9350–67. pmid:34877271
- View Article
- PubMed/NCBI
- Google Scholar
7. Rayner L, Price A, Evans A, Valsraj K, Hotopf M, Higginson IJ. Antidepressants for the treatment of depression in palliative care: systematic review and meta-analysis. Palliat Med. 2011 Jan 1;25(1):36–51. pmid:20935027
- View Article
- PubMed/NCBI
- Google Scholar
8. Schildkraut JJ. The catecholamine hypothesis of affective disorders: a review of supporting evidence. Am J Psychiatry. 1965 Nov;122(5):509–22. pmid:5319766
- View Article
- PubMed/NCBI
- Google Scholar
9. Kennis M, Gerritsen L, van Dalen M, Williams A, Cuijpers P, Bockting C. Prospective biomarkers of major depressive disorder: a systematic review and meta-analysis. Mol Psychiatry. 2020 Feb;25(2):321–38. pmid:31745238
- View Article
- PubMed/NCBI
- Google Scholar
10. Knorr U, Vinberg M, Kessing LV, Wetterslev J. Salivary cortisol in depressed patients versus control persons: A systematic review and meta-analysis. Psychoneuroendocrinology. 2010 Oct 1;35(9):1275–86. pmid:20447770
- View Article
- PubMed/NCBI
- Google Scholar
11. Sonino N, Fava GA. Psychiatric Disorders Associated with Cushing’s Syndrome. CNS Drugs. 2001 May 1;15(5):361–73.
- View Article
- Google Scholar
12. Sonino N, Fava GA, Belluardo P, Girelli ME, Boscaro M. Course of Depression in Cushing’s Syndrome: Response to Treatment and Comparison with Graves’ Disease. Horm Res Paediatr. 1993;39(5–6):202–6. pmid:8314204
- View Article
- PubMed/NCBI
- Google Scholar
13. Starkman MN, Schteingart DE, Schork MA. Depressed mood and other psychiatric manifestations of Cushing’s syndrome: relationship to hormone levels. Psychosom Med. 1981 Feb;43(1):3–18. pmid:6259680
- View Article
- PubMed/NCBI
- Google Scholar
14. Judd LL, Schettler PJ, Brown ES, Wolkowitz OM, Sternberg EM, Bender BG, et al. Adverse Consequences of Glucocorticoid Medication: Psychological, Cognitive, and Behavioral Effects. Am J Psychiatry. 2014 Oct;171(10):1045–51. pmid:25272344
- View Article
- PubMed/NCBI
- Google Scholar
15. Pariante CM, Lightman SL. The HPA axis in major depression: classical theories and new developments. Trends Neurosci. 2008 Sep;31(9):464–8. pmid:18675469
- View Article
- PubMed/NCBI
- Google Scholar
16. Miller WL. The Hypothalamic-Pituitary-Adrenal Axis: A Brief History. Horm Res Paediatr. 2018;89(4):212–23.
- View Article
- Google Scholar
17. Rubin RT, Phillips JJ, Sadow TF, McCracken JT. Adrenal gland volume in major depression. Increase during the depressive episode and decrease with successful treatment. Arch Gen Psychiatry. 1995 Mar;52(3):213–8. pmid:7872849
- View Article
- PubMed/NCBI
- Google Scholar
18. Rubin RT, Phillips JJ, McCracken JT, Sadow TF. Adrenal gland volume in major depression: Relationship to basal and stimulated pituitary-adrenal cortical axis function. Biol Psychiatry. 1996 Jul 15;40(2):89–97. pmid:8793041
- View Article
- PubMed/NCBI
- Google Scholar
19. Bardeleben U von Stalla GK, Müller OA, Holsboer F. Blunting of ACTH response to human CRH in depressed patients is avoided by metyrapone pretreatment. Biol Psychiatry. 1988 Nov 1;24(7):782–6. pmid:2852510
- View Article
- PubMed/NCBI
- Google Scholar
20. Kunugi H, Ida I, Owashi T, Kimura M, Inoue Y, Nakagawa S, et al. Assessment of the Dexamethasone/CRH Test as a State-Dependent Marker for Hypothalamic-Pituitary-Adrenal (HPA) Axis Abnormalities in Major Depressive Episode: A Multicenter Study. Neuropsychopharmacology. 2006 Jan;31(1):212–20. pmid:16123748
- View Article
- PubMed/NCBI
- Google Scholar
21. Heuser I, Yassouridis A, Holsboer F. The combined dexamethasone/CRH test: A refined laboratory test for psychiatric disorders. J Psychiatr Res. 1994 Jul 1;28(4):341–56. pmid:7877114
- View Article
- PubMed/NCBI
- Google Scholar
22. Sheline YI, Gado MH, Kraemer HC. Untreated Depression and Hippocampal Volume Loss. Am J Psychiatry. 2003 Aug;160(8):1516–8. pmid:12900317
- View Article
- PubMed/NCBI
- Google Scholar
23. Depression, antidepressants, and the shrinking hippocampus | PNAS [Internet]. [cited 2023 Jan 29]. Available from: https://www.pnas.org/doi/full/10.1073/pnas.231475998.
24. Bear MF, Abraham WC. Long-Term Depression in Hippocampus. Annu Rev Neurosci. 1996;19(1):437–62. pmid:8833450
- View Article
- PubMed/NCBI
- Google Scholar
25. Sahay A, Hen R. Adult hippocampal neurogenesis in depression. Nat Neurosci. 2007 Sep;10(9):1110–5. pmid:17726477
- View Article
- PubMed/NCBI
- Google Scholar
26. Jacobson L, Sapolsky R. The Role of the Hippocampus in Feedback Regulation of the Hypothalamic-Pituitary-Adrenocortical Axis*. Endocr Rev. 1991 May 1;12(2):118–34. pmid:2070776
- View Article
- PubMed/NCBI
- Google Scholar
27. Weidenfeld J, Ovadia H, Weidenfeld J, Ovadia H. The Role of the Amygdala in Regulating the Hypothalamic-Pituitary-Adrenal Axis [Internet]. The Amygdala—Where Emotions Shape Perception, Learning and Memories. IntechOpen; 2017 [cited 2023 Feb 19]. Available from: https://www.intechopen.com/chapters/54565.
- View Article
- Google Scholar
28. Sapolsky RM. The possibility of neurotoxicity in the hippocampus in major depression: a primer on neuron death. Biol Psychiatry. 2000 Oct 15;48(8):755–65. pmid:11063972
- View Article
- PubMed/NCBI
- Google Scholar
29. Anacker C, Zunszain PA, Carvalho LA, Pariante CM. The glucocorticoid receptor: Pivot of depression and of antidepressant treatment? Psychoneuroendocrinology. 2011 Apr;36(3):415–25. pmid:20399565
- View Article
- PubMed/NCBI
- Google Scholar
30. Gregrowicz J, Rogalska J. Glucocorticoid receptors function in the pathophysiology of brain hypoxia. Postępy Hig Med Dośw. 2019 Dec 31;73:838–49.
- View Article
- Google Scholar
31. Alon U. An Introduction to Systems Biology: Design Principles of Biological Circuits. CRC Press; 2019. 343 p.
32. Keener J, Sneyd J. Mathematical Physiology: I: Cellular Physiology. Springer Science & Business Media; 2008. 1067 p.
- View Article
- Google Scholar
33. Demic S, Cheng S. Modeling the Dynamics of Disease States in Depression. PLoS ONE. 2014 Oct 17;9(10):e110358. pmid:25330102
- View Article
- PubMed/NCBI
- Google Scholar
34. Kyrylov V, Severyanova LA, Vieira A. Modeling robust oscillatory behavior of the hypothalamic-pituitary-adrenal axis. IEEE Trans Biomed Eng. 2005 Dec;52(12):1977–83. pmid:16366221
- View Article
- PubMed/NCBI
- Google Scholar
35. Mathematical modeling of the hypothalamic-pituitary-adrenal system activity—PubMed [Internet]. [cited 2023 May 16]. Available from: https://pubmed.ncbi.nlm.nih.gov/16112688/.
36. Gupta S, Aslakson E, Gurbaxani BM, Vernon SD. Inclusion of the glucocorticoid receptor in a hypothalamic pituitary adrenal axis model reveals bistability. Theor Biol Med Model. 2007 Feb 14;4:8. pmid:17300722
- View Article
- PubMed/NCBI
- Google Scholar
37. Vinther F, Andersen M, Ottesen JT. The minimal model of the hypothalamic-pituitary-adrenal axis. J Math Biol. 2011 Oct;63(4):663–90. pmid:21107577
- View Article
- PubMed/NCBI
- Google Scholar
38. Rankin J, Walker JJ, Windle R, Lightman SL, Terry JR. Characterizing Dynamic Interactions between Ultradian Glucocorticoid Rhythmicity and Acute Stress Using the Phase Response Curve. PLOS ONE. 2012;7(2):e30978. pmid:22363526
- View Article
- PubMed/NCBI
- Google Scholar
39. Andersen M, Vinther F, Ottesen JT. Mathematical modeling of the hypothalamic–pituitary–adrenal gland (HPA) axis, including hippocampal mechanisms. Math Biosci. 2013 Nov 1;246(1):122–38. pmid:24012602
- View Article
- PubMed/NCBI
- Google Scholar
40. Savić D, Knežević G, Opačić G. A mathematical model of stress reaction: Individual differences in threshold and duration. Psychobiology. 2000 Dec 1;28(4):581–92.
- View Article
- Google Scholar
41. Kim LU D’Orsogna MR, Chou T. Onset, timing, and exposure therapy of stress disorders: mechanistic insight from a mathematical model of oscillating neuroendocrine dynamics. Biol Direct. 2016 Mar 25;11(1):13. pmid:27013324
- View Article
- PubMed/NCBI
- Google Scholar
42. Kim LU D’Orsogna MR, Chou T. Perturbing the Hypothalamic–Pituitary–Adrenal Axis: A Mathematical Model for Interpreting PTSD Assessment Tests. Comput Psychiatry Camb Mass. 2018 Feb 1;2:28–49.
- View Article
- Google Scholar
43. Bangsgaard EO, Ottesen JT. Patient specific modeling of the HPA axis related to clinical diagnosis of depression. Math Biosci. 2017 May 1;287:24–35. pmid:27816534
- View Article
- PubMed/NCBI
- Google Scholar
44. Cheng X, Dorsogna M, Chou T. Mathematical modeling of depressive disorders: Circadian driving, bistability and dynamical transitions. Comput Struct Biotechnol J. 2020 Nov 1;19. pmid:33510869
- View Article
- PubMed/NCBI
- Google Scholar
45. Karin O, Raz M, Tendler A, Bar A, Korem Kohanim Y, Milo T, et al. A new model for the HPA axis explains dysregulation of stress hormones on the timescale of weeks. Mol Syst Biol. 2020 Jul 16;16(7):e9510. pmid:32672906
- View Article
- PubMed/NCBI
- Google Scholar
46. Tendler A, Bar A, Mendelsohn-Cohen N, Karin O, Korem Kohanim Y, Maimon L, et al. Hormone seasonality in medical records suggests circannual endocrine circuits. Proc Natl Acad Sci. 2021 Feb 16;118(7):e2003926118. pmid:33531344
- View Article
- PubMed/NCBI
- Google Scholar
47. Hair cortisol as a biomarker for mood in bipolar disorder: a pilot study [Internet]. 2023 [cited 2023 Apr 2]. Available from: https://www.researchsquare.com.
- View Article
- Google Scholar
48. Radley JJ, Arias CM, Sawchenko PE. Regional Differentiation of the Medial Prefrontal Cortex in Regulating Adaptive Responses to Acute Emotional Stress. J Neurosci. 2006 Dec 13;26(50):12967–76. pmid:17167086
- View Article
- PubMed/NCBI
- Google Scholar
49. Zhe D, Fang H, Yuxiu S. Expressions of Hippocampal Mineralocorticoid Receptor (MR) and Glucocorticoid Receptor (GR) in the Single-Prolonged Stress-Rats. Acta Histochem Cytochem. 2008 Aug 28;41(4):89–95. pmid:18787639
- View Article
- PubMed/NCBI
- Google Scholar
50. Bar A, Karin O, Mayo A, Ben-Zvi D, Alon U. Rules for body fat interventions based on an operating point mechanism. iScience. 2023 Feb 17;26(2):106047. pmid:36818281
- View Article
- PubMed/NCBI
- Google Scholar
51. Korem Kohanim Y, Milo T, Raz M, Karin O, Bar A, Mayo A, et al. Dynamics of thyroid diseases and thyroid-axis gland masses. Mol Syst Biol. 2022 Aug;18(8):e10919. pmid:35938225
- View Article
- PubMed/NCBI
- Google Scholar
52. BioRender [Internet]. [cited 2023 Nov 8]. Available from: https://app.biorender.com/.
53. BioRender—Academic License [Internet]. [cited 2023 Nov 8]. Available from: https://www.biorender.com/academic-license.
54. Arnett MG, Muglia LM, Laryea G, Muglia LJ. Genetic Approaches to Hypothalamic-Pituitary-Adrenal Axis Regulation. Neuropsychopharmacology. 2016 Jan;41(1):245–60. pmid:26189452
- View Article
- PubMed/NCBI
- Google Scholar
55. Briassoulis G, Damjanovic S, Xekouki P, Lefebvre H, Stratakis CA. THE GLUCOCORTICOID RECEPTOR AND ITS EXPRESSION IN THE ANTERIOR PITUITARY AND THE ADRENAL CORTEX: A SOURCE OF VARIATION IN HYPOTHALAMIC-PITUITARY-ADRENAL AXIS FUNCTION; IMPLICATIONS FOR PITUITARY AND ADRENAL TUMORS. Endocr Pract Off J Am Coll Endocrinol Am Assoc Clin Endocrinol. 2011;17(6):941–8. pmid:21742609
- View Article
- PubMed/NCBI
- Google Scholar
56. Machado-Vieira R, Baumann J, Wheeler-Castillo C, Latov D, Henter ID, Salvadore G, et al. The Timing of Antidepressant Effects: A Comparison of Diverse Pharmacological and Somatic Treatments. Pharmaceuticals. 2010 Jan 6;3(1):19–41. pmid:27713241
- View Article
- PubMed/NCBI
- Google Scholar
57. Malhi GS, Bell E, Morris G, Hamilton A. The delay in response to antidepressant therapy: A window of opportunity? Aust N Z J Psychiatry. 2020 Feb 1;54(2):127–9.
- View Article
- Google Scholar
58. Nadeem HS, Attenburrow MJ, Cowen PJ. Comparison of the Effects of Citalopram and Escitalopram on 5-Ht-Mediated Neuroendocrine Responses. Neuropsychopharmacology. 2004 Sep;29(9):1699–703. pmid:15127082
- View Article
- PubMed/NCBI
- Google Scholar
59. Alenina N, Klempin F. The role of serotonin in adult hippocampal neurogenesis. Behav Brain Res. 2015 Jan 15;277:49–57. pmid:25125239
- View Article
- PubMed/NCBI
- Google Scholar
60. Murrough JW, Iosifescu DV, Chang LC, Al Jurdi RK, Green CE, Perez AM, et al. Antidepressant Efficacy of Ketamine in Treatment-Resistant Major Depression: A Two-Site Randomized Controlled Trial. Am J Psychiatry. 2013 Oct;170(10):1134–42. pmid:23982301
- View Article
- PubMed/NCBI
- Google Scholar
61. Sowa J, Kusek M, Bobula B, Hess G, Tokarski K. Ketamine Administration Reverses Corticosterone-Induced Alterations in Excitatory and Inhibitory Transmission in the Rat Dorsal Raphe Nucleus. Neural Plast. 2019 Aug 15;2019:3219490. pmid:31511771
- View Article
- PubMed/NCBI
- Google Scholar
62. Johnston CJ, Fitzgerald PJ, Gewarges JS, Watson BO, Spencer-Segal JL. Ketamine decreases HPA axis reactivity to a novel stressor in male but not female mice [Internet]. bioRxiv; 2021 [cited 2023 Feb 19]. p. 2021.06.29.450387. Available from: https://www.biorxiv.org/content/10.1101/2021.06.29.450387v1.
- View Article
- Google Scholar
63. Ostroff R, Kothari JS. Reversal of Non-Suppression of Cortisol Levels in a Patient With Refractory Depression Receiving Ketamine. Am J Psychiatry. 2015 Jan;172(1):95–6. pmid:25553498
- View Article
- PubMed/NCBI
- Google Scholar
64. Harno E, Gali Ramamoorthy T, Coll AP, White A. POMC: The Physiological Power of Hormone Processing. Physiol Rev. 2018 Oct 1;98(4):2381–430. pmid:30156493
- View Article
- PubMed/NCBI
- Google Scholar
65. Goodwin GM. Depression and associated physical diseases and symptoms. Dialogues Clin Neurosci. 2006 Jun;8(2):259–65. pmid:16889110
- View Article
- PubMed/NCBI
- Google Scholar
66. Thompson T, Correll CU, Gallop K, Vancampfort D, Stubbs B. Is Pain Perception Altered in People With Depression? A Systematic Review and Meta-Analysis of Experimental Pain Research. J Pain. 2016 Dec;17(12):1257–72. pmid:27589910
- View Article
- PubMed/NCBI
- Google Scholar
67. Bylsma L, Morris B, Rottenberg J. A meta-analysis of emotional reactivity in major depressive disorder. Clin Psychol Rev. 2008 May 1;28:676–91. pmid:18006196
- View Article
- PubMed/NCBI
- Google Scholar
68. Kao K, Tuladhar CT, Meyer JS, Tarullo AR. Emotion regulation moderates the association between parent and child hair cortisol concentrations. Dev Psychobiol. 2019 Nov;61(7):1064–78. pmid:30953366
- View Article
- PubMed/NCBI
- Google Scholar
69. Tsigos C, Chrousos GP. Hypothalamic–pituitary–adrenal axis, neuroendocrine factors and stress. J Psychosom Res. 2002 Oct 1;53(4):865–71. pmid:12377295
- View Article
- PubMed/NCBI
- Google Scholar
70. Arterbery A, Fergus D, Fogarty E, Mayberry J, Deitcher D, Kraus W, et al. Evolution of ligand specificity in vertebrate corticosteroid receptors. BMC Evol Biol. 2011 Jan 14;11:14. pmid:21232159
- View Article
- PubMed/NCBI
- Google Scholar
71. Koning ASCAM Buurstede JC, van Weert LTCM Meijer OC. Glucocorticoid and Mineralocorticoid Receptors in the Brain: A Transcriptional Perspective. J Endocr Soc. 2019 Jul 24;3(10):1917–30. pmid:31598572
- View Article
- PubMed/NCBI
- Google Scholar
72. Chronic Administration of Corticotropin-Releasing Factor Increases Pituitary Corticotroph Number* | Endocrinology | Oxford Academic [Internet]. [cited 2023 Sep 10]. Available from: https://academic.oup.com/endo/article/120/1/381/2540834.
73. The mitogenic effects of corticotrophin-releasing factor on the anterior pituitary gland of the rat in: Journal of Endocrinology Volume 118 Issue 2 (1988) [Internet]. [cited 2023 Sep 10]. Available from: https://joe.bioscientifica.com/view/journals/joe/118/2/joe_118_2_012.xml.
- View Article
- Google Scholar
74. Swann HG. THE PITUITARY-ADRENOCORTICAL RELATIONSHIP. Physiol Rev. 1940 Oct 1;20(4):493–521.
- View Article
- Google Scholar
75. Lotfi CFP, de Mendonca POR. Comparative Effect of ACTH and Related Peptides on Proliferation and Growth of Rat Adrenal Gland. Front Endocrinol [Internet]. 2016 [cited 2023 Sep 10];7. Available from: https://www.frontiersin.org/articles/10.3389/fendo.2016.00039. pmid:27242663
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Gutiérrez-Rojas L, Porras-Segovia A, Dunne H, Andrade-González N, Cervilla JA. Prevalence and correlates of major depressive disorder: a systematic review. Braz J Psychiatry. 2020 Aug 3;42(6):657–72. pmid:32756809
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Major Depression [Internet]. National Institute of Mental Health (NIMH). [cited 2023 Jan 29]. Available from: https://www.nimh.nih.gov/health/statistics/major-depression.

[ref3] 3. Thapar A, Eyre O, Patel V, Brent D. Depression in young people. The Lancet. 2022 Aug 20;400(10352):617–31.
View Article
Google Scholar

[7] View Article

[8] Google Scholar

[ref4] 4. Wong ML, Licinio J. Research and treatment approaches to depression. Nat Rev Neurosci. 2001 May;2(5):343–51. pmid:11331918
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref5] 5. Moncrieff J, Cooper RE, Stockmann T, Amendola S, Hengartner MP, Horowitz MA. The serotonin theory of depression: a systematic umbrella review of the evidence. Mol Psychiatry. 2022 Jul 20;1–14. pmid:35854107
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref6] 6. Karrouri R, Hammani Z, Benjelloun R, Otheman Y. Major depressive disorder: Validated treatments and future challenges. World J Clin Cases. 2021 Nov 6;9(31):9350–67. pmid:34877271
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref7] 7. Rayner L, Price A, Evans A, Valsraj K, Hotopf M, Higginson IJ. Antidepressants for the treatment of depression in palliative care: systematic review and meta-analysis. Palliat Med. 2011 Jan 1;25(1):36–51. pmid:20935027
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref8] 8. Schildkraut JJ. The catecholamine hypothesis of affective disorders: a review of supporting evidence. Am J Psychiatry. 1965 Nov;122(5):509–22. pmid:5319766
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref9] 9. Kennis M, Gerritsen L, van Dalen M, Williams A, Cuijpers P, Bockting C. Prospective biomarkers of major depressive disorder: a systematic review and meta-analysis. Mol Psychiatry. 2020 Feb;25(2):321–38. pmid:31745238
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref10] 10. Knorr U, Vinberg M, Kessing LV, Wetterslev J. Salivary cortisol in depressed patients versus control persons: A systematic review and meta-analysis. Psychoneuroendocrinology. 2010 Oct 1;35(9):1275–86. pmid:20447770
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref11] 11. Sonino N, Fava GA. Psychiatric Disorders Associated with Cushing’s Syndrome. CNS Drugs. 2001 May 1;15(5):361–73.
View Article
Google Scholar

[38] View Article

[39] Google Scholar

[ref12] 12. Sonino N, Fava GA, Belluardo P, Girelli ME, Boscaro M. Course of Depression in Cushing’s Syndrome: Response to Treatment and Comparison with Graves’ Disease. Horm Res Paediatr. 1993;39(5–6):202–6. pmid:8314204
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref13] 13. Starkman MN, Schteingart DE, Schork MA. Depressed mood and other psychiatric manifestations of Cushing’s syndrome: relationship to hormone levels. Psychosom Med. 1981 Feb;43(1):3–18. pmid:6259680
View Article
PubMed/NCBI
Google Scholar

[45] View Article

[46] PubMed/NCBI

[47] Google Scholar

[ref14] 14. Judd LL, Schettler PJ, Brown ES, Wolkowitz OM, Sternberg EM, Bender BG, et al. Adverse Consequences of Glucocorticoid Medication: Psychological, Cognitive, and Behavioral Effects. Am J Psychiatry. 2014 Oct;171(10):1045–51. pmid:25272344
View Article
PubMed/NCBI
Google Scholar

[49] View Article

[50] PubMed/NCBI

[51] Google Scholar

[ref15] 15. Pariante CM, Lightman SL. The HPA axis in major depression: classical theories and new developments. Trends Neurosci. 2008 Sep;31(9):464–8. pmid:18675469
View Article
PubMed/NCBI
Google Scholar

[53] View Article

[54] PubMed/NCBI

[55] Google Scholar

[ref16] 16. Miller WL. The Hypothalamic-Pituitary-Adrenal Axis: A Brief History. Horm Res Paediatr. 2018;89(4):212–23.
View Article
Google Scholar

[57] View Article

[58] Google Scholar

[ref17] 17. Rubin RT, Phillips JJ, Sadow TF, McCracken JT. Adrenal gland volume in major depression. Increase during the depressive episode and decrease with successful treatment. Arch Gen Psychiatry. 1995 Mar;52(3):213–8. pmid:7872849
View Article
PubMed/NCBI
Google Scholar

[60] View Article

[61] PubMed/NCBI

[62] Google Scholar

[ref18] 18. Rubin RT, Phillips JJ, McCracken JT, Sadow TF. Adrenal gland volume in major depression: Relationship to basal and stimulated pituitary-adrenal cortical axis function. Biol Psychiatry. 1996 Jul 15;40(2):89–97. pmid:8793041
View Article
PubMed/NCBI
Google Scholar

[64] View Article

[65] PubMed/NCBI

[66] Google Scholar

[ref19] 19. Bardeleben U von Stalla GK, Müller OA, Holsboer F. Blunting of ACTH response to human CRH in depressed patients is avoided by metyrapone pretreatment. Biol Psychiatry. 1988 Nov 1;24(7):782–6. pmid:2852510
View Article
PubMed/NCBI
Google Scholar

[68] View Article

[69] PubMed/NCBI

[70] Google Scholar

[ref20] 20. Kunugi H, Ida I, Owashi T, Kimura M, Inoue Y, Nakagawa S, et al. Assessment of the Dexamethasone/CRH Test as a State-Dependent Marker for Hypothalamic-Pituitary-Adrenal (HPA) Axis Abnormalities in Major Depressive Episode: A Multicenter Study. Neuropsychopharmacology. 2006 Jan;31(1):212–20. pmid:16123748
View Article
PubMed/NCBI
Google Scholar

[72] View Article

[73] PubMed/NCBI

[74] Google Scholar

[ref21] 21. Heuser I, Yassouridis A, Holsboer F. The combined dexamethasone/CRH test: A refined laboratory test for psychiatric disorders. J Psychiatr Res. 1994 Jul 1;28(4):341–56. pmid:7877114
View Article
PubMed/NCBI
Google Scholar

[76] View Article

[77] PubMed/NCBI

[78] Google Scholar

[ref22] 22. Sheline YI, Gado MH, Kraemer HC. Untreated Depression and Hippocampal Volume Loss. Am J Psychiatry. 2003 Aug;160(8):1516–8. pmid:12900317
View Article
PubMed/NCBI
Google Scholar

[80] View Article

[81] PubMed/NCBI

[82] Google Scholar

[ref23] 23. Depression, antidepressants, and the shrinking hippocampus | PNAS [Internet]. [cited 2023 Jan 29]. Available from: https://www.pnas.org/doi/full/10.1073/pnas.231475998.

[ref24] 24. Bear MF, Abraham WC. Long-Term Depression in Hippocampus. Annu Rev Neurosci. 1996;19(1):437–62. pmid:8833450
View Article
PubMed/NCBI
Google Scholar

[85] View Article

[86] PubMed/NCBI

[87] Google Scholar

[ref25] 25. Sahay A, Hen R. Adult hippocampal neurogenesis in depression. Nat Neurosci. 2007 Sep;10(9):1110–5. pmid:17726477
View Article
PubMed/NCBI
Google Scholar

[89] View Article

[90] PubMed/NCBI

[91] Google Scholar

[ref26] 26. Jacobson L, Sapolsky R. The Role of the Hippocampus in Feedback Regulation of the Hypothalamic-Pituitary-Adrenocortical Axis*. Endocr Rev. 1991 May 1;12(2):118–34. pmid:2070776
View Article
PubMed/NCBI
Google Scholar

[93] View Article

[94] PubMed/NCBI

[95] Google Scholar

[ref27] 27. Weidenfeld J, Ovadia H, Weidenfeld J, Ovadia H. The Role of the Amygdala in Regulating the Hypothalamic-Pituitary-Adrenal Axis [Internet]. The Amygdala—Where Emotions Shape Perception, Learning and Memories. IntechOpen; 2017 [cited 2023 Feb 19]. Available from: https://www.intechopen.com/chapters/54565.
View Article
Google Scholar

[97] View Article

[98] Google Scholar

[ref28] 28. Sapolsky RM. The possibility of neurotoxicity in the hippocampus in major depression: a primer on neuron death. Biol Psychiatry. 2000 Oct 15;48(8):755–65. pmid:11063972
View Article
PubMed/NCBI
Google Scholar

[100] View Article

[101] PubMed/NCBI

[102] Google Scholar

[ref29] 29. Anacker C, Zunszain PA, Carvalho LA, Pariante CM. The glucocorticoid receptor: Pivot of depression and of antidepressant treatment? Psychoneuroendocrinology. 2011 Apr;36(3):415–25. pmid:20399565
View Article
PubMed/NCBI
Google Scholar

[104] View Article

[105] PubMed/NCBI

[106] Google Scholar

[ref30] 30. Gregrowicz J, Rogalska J. Glucocorticoid receptors function in the pathophysiology of brain hypoxia. Postępy Hig Med Dośw. 2019 Dec 31;73:838–49.
View Article
Google Scholar

[108] View Article

[109] Google Scholar

[ref31] 31. Alon U. An Introduction to Systems Biology: Design Principles of Biological Circuits. CRC Press; 2019. 343 p.

[ref32] 32. Keener J, Sneyd J. Mathematical Physiology: I: Cellular Physiology. Springer Science & Business Media; 2008. 1067 p.
View Article
Google Scholar

[112] View Article

[113] Google Scholar

[ref33] 33. Demic S, Cheng S. Modeling the Dynamics of Disease States in Depression. PLoS ONE. 2014 Oct 17;9(10):e110358. pmid:25330102
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref34] 34. Kyrylov V, Severyanova LA, Vieira A. Modeling robust oscillatory behavior of the hypothalamic-pituitary-adrenal axis. IEEE Trans Biomed Eng. 2005 Dec;52(12):1977–83. pmid:16366221
View Article
PubMed/NCBI
Google Scholar

[119] View Article

[120] PubMed/NCBI

[121] Google Scholar

[ref35] 35. Mathematical modeling of the hypothalamic-pituitary-adrenal system activity—PubMed [Internet]. [cited 2023 May 16]. Available from: https://pubmed.ncbi.nlm.nih.gov/16112688/.

[ref36] 36. Gupta S, Aslakson E, Gurbaxani BM, Vernon SD. Inclusion of the glucocorticoid receptor in a hypothalamic pituitary adrenal axis model reveals bistability. Theor Biol Med Model. 2007 Feb 14;4:8. pmid:17300722
View Article
PubMed/NCBI
Google Scholar

[124] View Article

[125] PubMed/NCBI

[126] Google Scholar

[ref37] 37. Vinther F, Andersen M, Ottesen JT. The minimal model of the hypothalamic-pituitary-adrenal axis. J Math Biol. 2011 Oct;63(4):663–90. pmid:21107577
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref38] 38. Rankin J, Walker JJ, Windle R, Lightman SL, Terry JR. Characterizing Dynamic Interactions between Ultradian Glucocorticoid Rhythmicity and Acute Stress Using the Phase Response Curve. PLOS ONE. 2012;7(2):e30978. pmid:22363526
View Article
PubMed/NCBI
Google Scholar

[132] View Article

[133] PubMed/NCBI

[134] Google Scholar

[ref39] 39. Andersen M, Vinther F, Ottesen JT. Mathematical modeling of the hypothalamic–pituitary–adrenal gland (HPA) axis, including hippocampal mechanisms. Math Biosci. 2013 Nov 1;246(1):122–38. pmid:24012602
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref40] 40. Savić D, Knežević G, Opačić G. A mathematical model of stress reaction: Individual differences in threshold and duration. Psychobiology. 2000 Dec 1;28(4):581–92.
View Article
Google Scholar

[140] View Article

[141] Google Scholar

[ref41] 41. Kim LU D’Orsogna MR, Chou T. Onset, timing, and exposure therapy of stress disorders: mechanistic insight from a mathematical model of oscillating neuroendocrine dynamics. Biol Direct. 2016 Mar 25;11(1):13. pmid:27013324
View Article
PubMed/NCBI
Google Scholar

[143] View Article

[144] PubMed/NCBI

[145] Google Scholar

[ref42] 42. Kim LU D’Orsogna MR, Chou T. Perturbing the Hypothalamic–Pituitary–Adrenal Axis: A Mathematical Model for Interpreting PTSD Assessment Tests. Comput Psychiatry Camb Mass. 2018 Feb 1;2:28–49.
View Article
Google Scholar

[147] View Article

[148] Google Scholar

[ref43] 43. Bangsgaard EO, Ottesen JT. Patient specific modeling of the HPA axis related to clinical diagnosis of depression. Math Biosci. 2017 May 1;287:24–35. pmid:27816534
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref44] 44. Cheng X, Dorsogna M, Chou T. Mathematical modeling of depressive disorders: Circadian driving, bistability and dynamical transitions. Comput Struct Biotechnol J. 2020 Nov 1;19. pmid:33510869
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

[ref45] 45. Karin O, Raz M, Tendler A, Bar A, Korem Kohanim Y, Milo T, et al. A new model for the HPA axis explains dysregulation of stress hormones on the timescale of weeks. Mol Syst Biol. 2020 Jul 16;16(7):e9510. pmid:32672906
View Article
PubMed/NCBI
Google Scholar

[158] View Article

[159] PubMed/NCBI

[160] Google Scholar

[ref46] 46. Tendler A, Bar A, Mendelsohn-Cohen N, Karin O, Korem Kohanim Y, Maimon L, et al. Hormone seasonality in medical records suggests circannual endocrine circuits. Proc Natl Acad Sci. 2021 Feb 16;118(7):e2003926118. pmid:33531344
View Article
PubMed/NCBI
Google Scholar

[162] View Article

[163] PubMed/NCBI

[164] Google Scholar

[ref47] 47. Hair cortisol as a biomarker for mood in bipolar disorder: a pilot study [Internet]. 2023 [cited 2023 Apr 2]. Available from: https://www.researchsquare.com.
View Article
Google Scholar

[166] View Article

[167] Google Scholar

[ref48] 48. Radley JJ, Arias CM, Sawchenko PE. Regional Differentiation of the Medial Prefrontal Cortex in Regulating Adaptive Responses to Acute Emotional Stress. J Neurosci. 2006 Dec 13;26(50):12967–76. pmid:17167086
View Article
PubMed/NCBI
Google Scholar

[169] View Article

[170] PubMed/NCBI

[171] Google Scholar

[ref49] 49. Zhe D, Fang H, Yuxiu S. Expressions of Hippocampal Mineralocorticoid Receptor (MR) and Glucocorticoid Receptor (GR) in the Single-Prolonged Stress-Rats. Acta Histochem Cytochem. 2008 Aug 28;41(4):89–95. pmid:18787639
View Article
PubMed/NCBI
Google Scholar

[173] View Article

[174] PubMed/NCBI

[175] Google Scholar

[ref50] 50. Bar A, Karin O, Mayo A, Ben-Zvi D, Alon U. Rules for body fat interventions based on an operating point mechanism. iScience. 2023 Feb 17;26(2):106047. pmid:36818281
View Article
PubMed/NCBI
Google Scholar

[177] View Article

[178] PubMed/NCBI

[179] Google Scholar

[ref51] 51. Korem Kohanim Y, Milo T, Raz M, Karin O, Bar A, Mayo A, et al. Dynamics of thyroid diseases and thyroid-axis gland masses. Mol Syst Biol. 2022 Aug;18(8):e10919. pmid:35938225
View Article
PubMed/NCBI
Google Scholar

[181] View Article

[182] PubMed/NCBI

[183] Google Scholar

[ref52] 52. BioRender [Internet]. [cited 2023 Nov 8]. Available from: https://app.biorender.com/.

[ref53] 53. BioRender—Academic License [Internet]. [cited 2023 Nov 8]. Available from: https://www.biorender.com/academic-license.

[ref54] 54. Arnett MG, Muglia LM, Laryea G, Muglia LJ. Genetic Approaches to Hypothalamic-Pituitary-Adrenal Axis Regulation. Neuropsychopharmacology. 2016 Jan;41(1):245–60. pmid:26189452
View Article
PubMed/NCBI
Google Scholar

[187] View Article

[188] PubMed/NCBI

[189] Google Scholar

[ref55] 55. Briassoulis G, Damjanovic S, Xekouki P, Lefebvre H, Stratakis CA. THE GLUCOCORTICOID RECEPTOR AND ITS EXPRESSION IN THE ANTERIOR PITUITARY AND THE ADRENAL CORTEX: A SOURCE OF VARIATION IN HYPOTHALAMIC-PITUITARY-ADRENAL AXIS FUNCTION; IMPLICATIONS FOR PITUITARY AND ADRENAL TUMORS. Endocr Pract Off J Am Coll Endocrinol Am Assoc Clin Endocrinol. 2011;17(6):941–8. pmid:21742609
View Article
PubMed/NCBI
Google Scholar

[191] View Article

[192] PubMed/NCBI

[193] Google Scholar

[ref56] 56. Machado-Vieira R, Baumann J, Wheeler-Castillo C, Latov D, Henter ID, Salvadore G, et al. The Timing of Antidepressant Effects: A Comparison of Diverse Pharmacological and Somatic Treatments. Pharmaceuticals. 2010 Jan 6;3(1):19–41. pmid:27713241
View Article
PubMed/NCBI
Google Scholar

[195] View Article

[196] PubMed/NCBI

[197] Google Scholar

[ref57] 57. Malhi GS, Bell E, Morris G, Hamilton A. The delay in response to antidepressant therapy: A window of opportunity? Aust N Z J Psychiatry. 2020 Feb 1;54(2):127–9.
View Article
Google Scholar

[199] View Article

[200] Google Scholar

[ref58] 58. Nadeem HS, Attenburrow MJ, Cowen PJ. Comparison of the Effects of Citalopram and Escitalopram on 5-Ht-Mediated Neuroendocrine Responses. Neuropsychopharmacology. 2004 Sep;29(9):1699–703. pmid:15127082
View Article
PubMed/NCBI
Google Scholar

[202] View Article

[203] PubMed/NCBI

[204] Google Scholar

[ref59] 59. Alenina N, Klempin F. The role of serotonin in adult hippocampal neurogenesis. Behav Brain Res. 2015 Jan 15;277:49–57. pmid:25125239
View Article
PubMed/NCBI
Google Scholar

[206] View Article

[207] PubMed/NCBI

[208] Google Scholar

[ref60] 60. Murrough JW, Iosifescu DV, Chang LC, Al Jurdi RK, Green CE, Perez AM, et al. Antidepressant Efficacy of Ketamine in Treatment-Resistant Major Depression: A Two-Site Randomized Controlled Trial. Am J Psychiatry. 2013 Oct;170(10):1134–42. pmid:23982301
View Article
PubMed/NCBI
Google Scholar

[210] View Article

[211] PubMed/NCBI

[212] Google Scholar

[ref61] 61. Sowa J, Kusek M, Bobula B, Hess G, Tokarski K. Ketamine Administration Reverses Corticosterone-Induced Alterations in Excitatory and Inhibitory Transmission in the Rat Dorsal Raphe Nucleus. Neural Plast. 2019 Aug 15;2019:3219490. pmid:31511771
View Article
PubMed/NCBI
Google Scholar

[214] View Article

[215] PubMed/NCBI

[216] Google Scholar

[ref62] 62. Johnston CJ, Fitzgerald PJ, Gewarges JS, Watson BO, Spencer-Segal JL. Ketamine decreases HPA axis reactivity to a novel stressor in male but not female mice [Internet]. bioRxiv; 2021 [cited 2023 Feb 19]. p. 2021.06.29.450387. Available from: https://www.biorxiv.org/content/10.1101/2021.06.29.450387v1.
View Article
Google Scholar

[218] View Article

[219] Google Scholar

[ref63] 63. Ostroff R, Kothari JS. Reversal of Non-Suppression of Cortisol Levels in a Patient With Refractory Depression Receiving Ketamine. Am J Psychiatry. 2015 Jan;172(1):95–6. pmid:25553498
View Article
PubMed/NCBI
Google Scholar

[221] View Article

[222] PubMed/NCBI

[223] Google Scholar

[ref64] 64. Harno E, Gali Ramamoorthy T, Coll AP, White A. POMC: The Physiological Power of Hormone Processing. Physiol Rev. 2018 Oct 1;98(4):2381–430. pmid:30156493
View Article
PubMed/NCBI
Google Scholar

[225] View Article

[226] PubMed/NCBI

[227] Google Scholar

[ref65] 65. Goodwin GM. Depression and associated physical diseases and symptoms. Dialogues Clin Neurosci. 2006 Jun;8(2):259–65. pmid:16889110
View Article
PubMed/NCBI
Google Scholar

[229] View Article

[230] PubMed/NCBI

[231] Google Scholar

[ref66] 66. Thompson T, Correll CU, Gallop K, Vancampfort D, Stubbs B. Is Pain Perception Altered in People With Depression? A Systematic Review and Meta-Analysis of Experimental Pain Research. J Pain. 2016 Dec;17(12):1257–72. pmid:27589910
View Article
PubMed/NCBI
Google Scholar

[233] View Article

[234] PubMed/NCBI

[235] Google Scholar

[ref67] 67. Bylsma L, Morris B, Rottenberg J. A meta-analysis of emotional reactivity in major depressive disorder. Clin Psychol Rev. 2008 May 1;28:676–91. pmid:18006196
View Article
PubMed/NCBI
Google Scholar

[237] View Article

[238] PubMed/NCBI

[239] Google Scholar

[ref68] 68. Kao K, Tuladhar CT, Meyer JS, Tarullo AR. Emotion regulation moderates the association between parent and child hair cortisol concentrations. Dev Psychobiol. 2019 Nov;61(7):1064–78. pmid:30953366
View Article
PubMed/NCBI
Google Scholar

[241] View Article

[242] PubMed/NCBI

[243] Google Scholar

[ref69] 69. Tsigos C, Chrousos GP. Hypothalamic–pituitary–adrenal axis, neuroendocrine factors and stress. J Psychosom Res. 2002 Oct 1;53(4):865–71. pmid:12377295
View Article
PubMed/NCBI
Google Scholar

[245] View Article

[246] PubMed/NCBI

[247] Google Scholar

[ref70] 70. Arterbery A, Fergus D, Fogarty E, Mayberry J, Deitcher D, Kraus W, et al. Evolution of ligand specificity in vertebrate corticosteroid receptors. BMC Evol Biol. 2011 Jan 14;11:14. pmid:21232159
View Article
PubMed/NCBI
Google Scholar

[249] View Article

[250] PubMed/NCBI

[251] Google Scholar

[ref71] 71. Koning ASCAM Buurstede JC, van Weert LTCM Meijer OC. Glucocorticoid and Mineralocorticoid Receptors in the Brain: A Transcriptional Perspective. J Endocr Soc. 2019 Jul 24;3(10):1917–30. pmid:31598572
View Article
PubMed/NCBI
Google Scholar

[253] View Article

[254] PubMed/NCBI

[255] Google Scholar

[ref72] 72. Chronic Administration of Corticotropin-Releasing Factor Increases Pituitary Corticotroph Number* | Endocrinology | Oxford Academic [Internet]. [cited 2023 Sep 10]. Available from: https://academic.oup.com/endo/article/120/1/381/2540834.

[ref73] 73. The mitogenic effects of corticotrophin-releasing factor on the anterior pituitary gland of the rat in: Journal of Endocrinology Volume 118 Issue 2 (1988) [Internet]. [cited 2023 Sep 10]. Available from: https://joe.bioscientifica.com/view/journals/joe/118/2/joe_118_2_012.xml.
View Article
Google Scholar

[258] View Article

[259] Google Scholar

[ref74] 74. Swann HG. THE PITUITARY-ADRENOCORTICAL RELATIONSHIP. Physiol Rev. 1940 Oct 1;20(4):493–521.
View Article
Google Scholar

[261] View Article

[262] Google Scholar

[ref75] 75. Lotfi CFP, de Mendonca POR. Comparative Effect of ACTH and Related Peptides on Proliferation and Growth of Rat Adrenal Gland. Front Endocrinol [Internet]. 2016 [cited 2023 Sep 10];7. Available from: https://www.frontiersin.org/articles/10.3389/fendo.2016.00039. pmid:27242663
View Article
PubMed/NCBI
Google Scholar

[264] View Article

[265] PubMed/NCBI

[266] Google Scholar

Figures

Abstract

Author summary

Introduction

Results

Toggle switch model of the HPA—CNS system

Model shows bistability between euthymic and depressed states

Prolonged but not acute stress can trigger self-sustained depression

Stress must be reduced below baseline to overcome depression due to a hysteresis effect

Model explains dysregulation of endocrine tests in MDD

Model explains timescale of weeks for antidepressant action

Discussion

Material and methods

HPA-CNS toggle switch model

Animation of phase portrait under parameter changes

Supporting information

S1 File. Hysteresis.

S2 File. CRH test graph.

S3 File. DEX/CRH test graph.

S4 File. CRH test data.

S5 File. DEX/CRH test data.

S6 File. HPA phase portrait.

S7 File. Dynamic of brief stress.

S8 File. Dynamic of prolonged stress.

S9 File. Dynamic of brief stress reduction.

S10 File. Dynamic of prolonged stress reduction.

S11 File. Dynamic of short treatment.

S12 File. Dynamic of prolonged treatment.

References