AJP - Regu Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Regul Integr Comp Physiol 293: R83-R98, 2007. First published April 25, 2007; doi:10.1152/ajpregu.00907.2006
0363-6119/07 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Figures and Tables
Right arrow All Versions of this Article:
293/1/R83    most recent
00907.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via ISI Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Peters, A.
Right arrow Articles by Fehm, H. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Peters, A.
Right arrow Articles by Fehm, H. L.

APPETITE, OBESITY, DIGESTION, AND METABOLISM

The principle of homeostasis in the hypothalamus-pituitary-adrenal system: new insight from positive feedback

A. Peters,1 M. Conrad,1,2 C. Hubold,1 U. Schweiger,3 B. Fischer,2 and H. L. Fehm1

1Department of Internal Medicine I, Endocrinology, 2Institute for Mathematics, and 3Clinic of Psychiatry and Psychotherapy, University of Luebeck, Luebeck, Germany

Submitted 28 December 2006 ; accepted in final form 23 April 2007


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Feedback control, both negative and positive, is a fundamental feature of biological systems. Some of these systems strive to achieve a state of equilibrium or "homeostasis". The major endocrine systems are regulated by negative feedback, a process believed to maintain hormonal levels within a relatively narrow range. Positive feedback is often thought to have a destabilizing effect. Here, we present a "principle of homeostasis," which makes use of both positive and negative feedback loops. To test the hypothesis that this homeostatic concept is valid for the regulation of cortisol, we assessed experimental data in humans with different conditions (gender, obesity, endocrine disorders, medication) and analyzed these data by a novel computational approach. We showed that all obtained data sets were in agreement with the presented concept of homeostasis in the hypothalamus-pituitary-adrenal axis. According to this concept, a homeostatic system can stabilize itself with the help of a positive feedback loop. The brain mineralocorticoid and glucocorticoid receptors—with their known characteristics—fulfill the key functions in the homeostatic concept: binding cortisol with high and low affinities, acting in opposing manners, and mediating feedback effects on cortisol. This study supports the interaction between positive and negative feedback loops in the hypothalamus-pituitary-adrenal system and in this way sheds new light on the function of dual receptor regulation. Current knowledge suggests that this principle of homeostasis could also apply to other biological systems.

mathematical modeling; mineralocorticoid receptor; glucocorticoid receptor; obesity


HORMONES, NEUROPEPTIDES, NEUROTRANSMITTERS, and growth factors help the organism to achieve a state of equilibrium, which is referred to as homeostasis (18). These messengers (e.g., glucocorticoids, neurotrophins) typically bind as ligands to at least two different kinds of receptors that display different affinities and in this way convey their stabilizing signals. The principle, however, that underlies such homeostatic processes and that may explain why there are two receptors, is a long-standing mystery. Therefore, we assessed experimental data and analyzed these by a novel computational approach, showing that the data were in agreement with a fundamental homeostatic principle in the hypothalamic-pituitary-adrenal (HPA) system. In addition, we provide evidence from the literature that this concept could apply to other biological systems.

Usually, equilibria are thought to be driven by negative feedback only, for example, a deviated pendulum suffers a retrieving force to its stationary state. Biological systems, however, often display the typical components: a ligand, a low-affinity receptor, and a high-affinity receptor. The presence of such components suggests a general "principle of homeostasis" which includes positive feedback in particular ranges of ligand concentrations. As a recent discovery, positive feedback has been shown to be a ubiquitous signal transduction motif that allows systems to convert graded inputs into decisive, all-or-none outputs (2, 13). Thus, focusing on positive feedback already has shed light on another class of processes, that is, nonhomeostatic "multistable" processes like cellular on-off switches. Here we elaborate a basic model of homeostasis and introduce positive feedback loops as key mechanisms of homeostasis (Fig. 1, A and B; see also Box 1).
Box 1 The principle of biological homeostasis

A, B, and C are molecules: one ligand, two receptors.

Rule 1: The ligand A binds with high affinity to a receptor B, but with low affinity to a receptor C.

Rule 2: The two ligand-bound receptors act in opposing manners.

Rule 3: Ligand-bound B receptors increase* while ligand-bound C receptors decrease{dagger} the concentration of A. These interactions result in a homeostatic state of A.

*Positive feedback; {dagger}negative feedback

 


Figure 1
View larger version (38K):
[in this window]
[in a new window]

 
Fig. 1. Schematic representation of the "principle of homeostasis". A: biological homeostatic system stabilizes itself by producing ligands, A, which exert both positive feedback via high-affinity receptors, B, and negative feedback via low-affinity receptors, C. B: efficacies of the B and C receptors depend on the concentration of A (top); the combined feedback effect of both receptors (bottom) also depends on the concentration of A. The feedback effect is stimulatory at low concentrations of A, but inhibitory at high concentrations of A. C: flowchart of the hypothalamo-pituitary-adrenal (HPA) system. Five parameters are indicated (gray), which describe the efficacies of action at the respective sites. The brain-pituitary system releases ACTH, which stimulates the adrenal glands to secrete cortisol. In a long feedback loop, cortisol exerts positive feedback via mineralocorticoid receptor (MR) and negative via glucocorticoid receptor (GR) on the brain-pituitary system. In a short feedback loop, glucocorticoids bind to GR in the adrenal cortex, and in so doing, suppress adrenal glucocorticoid release (63, 75). External corticotropin-releasing-hormone (CRH) administration stimulates the brain-pituitary system to release ACTH. D: brain MR and GR efficacies depend on cortisol concentrations. We regard the brain MR as the major representative member of the nuclear receptor superfamily, which binds cortisol with high affinity. For the sake of simplicity, we attribute the effects of all such high-affinity cortisol receptors to MR. Accordingly, we consider the brain GR as a major representative member binding cortisol with low affinity. The maximal efficacies MRE maxbrain and GRE maxbrain are indicated. These efficacies represent the integrated actions of all MRs and GRs localized in different brain regions (e.g., amygdala, hippocampus, hypothalamus, and pituitary). Thus, computing such global MR and GR efficacies, which describe the actions on plasma ACTH, does not allow localizing the specific site of action within the brain-pituitary system. We designate the ratio of global MRE maxbrain and GRE maxbrain as "MR:GR balance". MREC50 denotes the molar concentration of cortisol, which produces 50% of the maximum possible MR-mediated response, that is, the change of ACTH in time. GREC50 is the respective cortisol concentration for GR-mediated changes of ACTH in time. The summation curve shows the combined feedback effect of both receptors on the brain-pituitary system, which also depends on cortisol concentrations. The brain setpoint of the HPA system (the zero indicated by the vertical arrow) is characterized by equilibrium between positive and negative feedback.

 
Which components of the HPA system fulfill the functions of A, B, and C in the general rules (Box 1)? The brain-pituitary system via ACTH stimulates the adrenal glands, increasing cortisol (Fig. 1C). In a feedback loop, cortisol acts on the brain by binding as a ligand to cerebral glucocorticoid receptors: at low concentrations, cortisol binds to the mineralocorticoid receptor (MR), and at higher concentrations, cortisol binds to the glucocorticoid receptor (GR) (Fig. 1D). MR and GR descend from a gene duplication deep in the vertebrate lineage and function as dual cortisol sensors (14). Starting in the year 1968 with the milestone paper of Bruce McEwen and colleagues (67), a large number of researchers have since managed to characterize the two brain receptor subtypes both biochemically and functionally. High densities of MR were found in the hippocampus, dentate gyrus, lateral septum, and amygdala, while abundant GR were detected in all brain regions (74).

It is known that the two corticosteroid receptors MR and GR operate in an opposing manner on several processes such as synaptic transmission (28), synaptic plasticity (4, 52, 53, 76, 77, 92), and cell survival (25). It is unknown, however, whether brain MR and GR also operate counteractively in the central feedback on the HPA system. Some scientists have referred to pharmacological studies (3, 31, 35, 47, 73, 81, 89, 99) and have concluded that MR and GR both exert inhibitory feedback in the HPA system. Ron de Kloet and colleagues (28) emphasized the "proactive" mode of MR and the "reactive" mode of GR and that the balance between MR and GR is essential for cell homeostasis, mental performance, and health. Nevertheless, he did not bring up positive feedback in this context. Considering the opposing nature of MR and GR in other biological processes, we assume here that brain MR and GR make use of positive and negative feedback in the HPA system, respectively. We, therefore, formulate the homeostatic rules with respect to the HPA system by replacing A, B, and C by cortisol, MR, and GR (Box 2).
Box 2 The principle of homeostasis in the HPA system

Rule 1: Cortisol binds at low concentrations to the brain MR and only at high concentrations to the brain GR.

Rule 2: Activated brain MR and GR operate in an opposing manner.

Rule 3: Cortisol raises* its own serum concentration via activated brain MR, while it reduces{dagger} it via activated brain GR.

These interactions result in a homeostatic state of cortisol.

*Positive feedback; {dagger}negative feedback

 

To test the hypothesis that the "principle of homeostasis" applies to the regulation of cortisol, we applied a computational approach to experimental data (Fig. 2) :


Figure 2
View larger version (29K):
[in this window]
[in a new window]

 
Fig. 2. Flowchart illustrating the working steps used in this paper. A: we formulated the principle of homeostasis using two differential equations with known constants and up to five unknown parameters (gray). B: we obtained experimental data, here exemplified by data from a CRH challenge test in a healthy man (case C.H.). C: we identified the unknown parameters (gray) for each individual (e.g., Case C.H.) by analyzing his/her experimental data. D: with these identified parameters, we solved the differential equations: the solution constitutes a predicted time-dependent concentration of ACTH and cortisol (solid line). We compared the theoretical prediction obtained (solid line) with the experimental data (data dots), computed the prediction error, and tested the hypothesis (significant prediction error).

 
In order to test the hypothesis over a broad scope, we assessed experimental data under various physiological, pharmacological, and pathological conditions. We obtained data from healthy women, healthy men, and obese men during a corticotropin-releasing hormone (CRH) challenge test (Fig. 3, A and B); we also referred to data that we had already published earlier, obtained in healthy subjects during MR blockade (Fig. 3C) (96), as well as in patients with Addison's disease and in adrenalectomized patients with Cushing's disease during a cortisol-infusion challenge (Fig. 3D) (38).


Figure 3
Figure 3
View larger version (50K):
[in this window]
[in a new window]

 
Fig. 3. Experimental data and theoretical predictions of ACTH and cortisol concentrations during various challenge tests. In all conditions, the experimental data (left) are similar to the theoretical predictions (right). A: CRH challenge in healthy women and men. Inset, right: Feedback-analysis showing the dose-response curves for these women and men: the change of ACTH in time (fmol·l–1·s–1) depends on serum cortisol (nmol/l). These functions are computed using the differential equations on the basis of the experimental data. The shaded area indicates the range of cortisol concentrations that is covered by experimental observations. The curves depicted outside this range are extrapolated with the help of data-based parameter identification and the differential equations. B: CRH challenge in obese and healthy men. C: treatment with canrenoate or placebo in healthy men. Inset, right: feedback-analysis illustrates that the biphasic dose-response curve (cortisol -> change in ACTH in time) is shifted to the right. With the right shift of the biphasic curve, the brain setpoint of the HPA system (vertical arrow) is altered, and homeostasis is achieved only at higher cortisol concentrations. D: continuous hydrocortisone infusion (50 mg/h, 0–120 min) in patients with Addison's disease or in patients with Cushing's disease/ADX. Inset, right: note the broad range of cortisol concentrations that is covered by experimental observations (shaded area). Error bars indicate SE.

 

    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Differential equations. To mathematically analyze experimental data sets, we established a system involving two differential equations. The rationale for the choice of terms in the equations is based on the following aspects: Basically, the equations represent the three rules (Box 2). We used a limited number of extensions that take into account specific physiological features of the HPA system, that is, ACTH as an intermediary signal posed in the brain-to-adrenal pathway, a clearance of hormones from blood, and a peripheral feedback at the adrenal level. As commented in the discussion, we deliberately abstained from adding extra compartments (regions of the brain) or extraregulatory processes (e.g., heterodimerization). For the denotations used in the equations, see Table 1:

Formula 1(1)


View this table:
[in this window]
[in a new window]

 
Table 1. Functions, parameters, and constants used in the differential equations

 
The term [dACTH(t)/dt] denotes the derivative of plasma ACTH concentrations over a time t. This change of ACTH in time depends on 1) ACTH and its first-order clearance c1 from plasma (see Table 1), 2) a positive MR-mediated effect, which is a function of CORT(t{tau}1) ({tau}1 indicates a delay time) and the pharmacokinetic parameters MRFormula 1 MREC50, and r, 3) a negative central GR effect, which is also a function of CORT(t{tau}1) and the respective pharmacokinetic parameters GRFormula 1GREC50, and r, and finally on an external disturbance e1(t), which is a typical gamma distribution (21) acting on ACTH, for example, an injection of CRH (with the parameter CRHFormula 1

Formula 2(2)

The term [dCORT(t)/dt] denotes the derivative of serum cortisol concentrations over time t. This change in cortisol over time depends on 1) cortisol and its first-order clearance c2 from plasma, 2) a positive ACTH effect, which is a function of ACTH(t{tau}2) ({tau}2 indicates a delay time) and the parameter ACTHFormula 23) a negative peripheral GR effect, which is a function of CORT(t{tau}1) and the parameter GRFormula 2 and finally on an external disturbance e2(t) acting on serum cortisol, for example, an infusion of cortisol (with the pharmacokinetic parameter s).

Equations 1 and 2 are "delay differential equations" (9). In these equations, ACTH and cortisol function as "integral controllers" (34). In the forward pathway, there is first a delay time {tau}2 (1.2 min) until ACTH has an effect on the derivative of serum cortisol concentration over time [dCORT(t)/dt]. Second, the full effect of ACTH on serum cortisol concentration per se is carried out time-lagged by this integral controller. This delayed effect individually arises from the constants and identified parameters in the differential equations. To estimate the total time lag in the forward pathway, we assessed the peak-to-peak phase shift between the ACTH and cortisol curve and found the lag to be 25.0 min (Fig. 2D). In the feedback pathway, there are similar time lags: a delay time {tau}1 (2.0 min) until cortisol influences the derivative of plasma ACTH concentrations [dACTH(t)/dt], and another time lag until cortisol measurably affects plasma ACTH concentrations per se. To estimate the total time lag in the brain corticosteroid feedback, we assessed the peak-to-nadir phase shift between the cortisol and ACTH curve and found the lag to be 38.0 min (Supplemental Figure S3 for this article can be found online at the American Journal of Physiology-Regulatory, Integrative and Comparative Physiology website.). The differential equations were solved numerically using MATLAB 7.0.

Data collection. We performed CRH challenge-tests in 10 normal men, 10 normal women, and 7 obese men. Exclusion criteria were chronic or acute physical and mental illness, alcohol or drug abuse, smoking, competitive sports, exceptional physical or psychological stress, and current medication of any kind. (Group characteristics are listed in Supplemental Table 1, which can be found online at the American Journal of Physiology–Regulatory, Integrative and Comparative Physiology website). All data collection procedures were reviewed and approved by the local ethics committee, and all subjects provided written informed consent.

The participants were instructed to abstain from alcohol, not to perform any kind of exhausting physical activity, and to go to bed no later than 11:00 PM on the day preceding the study. On the days of the experiments, subjects arrived at noontime in the medical research unit and received a standardized meal. A cannula was inserted into an antecubital vein at 12:00 PM. Afterward, each participant remained fasted except for the free intake of water. Three baseline blood samples for determining plasma ACTH and serum cortisol were collected from 3:30 PM until 4:00 PM. After the baseline period, CRH (CRH Ferring, Ferring Arzneimittel, Kiel, Germany) was intravenously injected at 4:00 PM at a dose of 1 µg/kg adapted to body mass. Blood samples were subsequently drawn at short intervals of 5 min until 5:30 PM and at intervals of 15 min until the end of the test at 8:00 PM.

All blood samples were immediately centrifuged and the supernatants were stored at –24°C until assaying. Cortisol serum concentrations were determined by electrochemiluminescence immunoassay in conjunction with the Elecsys 2010 (Roche Molecular Biochemical Diagnostics, Mannheim, Germany; intraassay CV < 1.3%; interassay CV < 1.5%), and plasma ACTH concentrations were measured by an immunoluminometric two-step assay (Lumitest ACTH, Brahms Diagnostics, Berlin, Germany; intraassay CV < 6.0%; interassay CV < 8.0%). Corticosteroid-binding globulin (CBG) was measured by a highly sensitive RIA (IBL Immuno-Biological Laboratories, Hamburg, Germany; intraassay CV < 3.8% and interassay CV < 5.0%).

The experimental methods that we had used in the study on MR blockade and the study on patients with patients with Addison's and Cushing's disease are described in the supplementary methods, which can be found online at the American Journal of Physiology–Regulatory, Integrative and Comparative Physiology website.

Parameter identification. The parameter identification method is a mathematical procedure for obtaining optimal parameters from the comparison of experimental data and a system of differential equations. The parameter identification procedure repeats three steps until an optimal parameter set has been found. We begin with a set of n parameters on an n-dimensional grid.

In the first step of the algorithm, approximate concentration profiles ACTH(t) and CORT(t) are created. For this purpose, conventional numeric "delay differential-equation-solvers" are employed (9, 46). In a second step, the algorithm compares the differences between the obtained concentration profiles with the actual experimental ACTH and cortisol data from a test subject. This difference provides a measure for the prediction error. The difference is calculated using the root-mean-square procedure. The third step of the algorithm modifies the parameter values so that the prediction error decreases. Here, we employed what we refer to as direct exploratory methods of "numeric optimization" (39).

These three steps are repeated until the difference can no longer be reduced. Then we carry out this procedure for all points of the n-dimensional grid. The outcome of the parameter identification method is the parameter set that provides the lowest prediction error. A prediction error <0.05 is considered as high agreement between experimental data and theoretical prediction. In the sense of the used error measure, an error <0.20 can be regarded as satisfactory, a larger error as unacceptable agreement. We apply this procedure for every test individual separately. The parameter identification was done by MATLAB 7.0.

Statistical analysis. To compare data among different groups of subjects, we used the following statistical tools: Dependent and independent Student's t-tests were used for comparing characteristics and identified parameters between subject groups. One sample t-test was employed to test whether the mean of one variable differs from a constant. ANOVA was performed for comparing parameters among three groups. P values <0.05 were considered to be significant. Values are expressed as means ± SE. The analysis was done with SPSS 12.0.


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Experiments in healthy humans. When CRH was administered to healthy subjects (characteristics found online in supplementary Table 1), ACTH profiles in men and women were similar, although women had slightly higher cortisol concentrations in the recovery phase (Fig. 3A, left). The ACTH curves reflect the combined effects of CRH and cortisol on the brain-pituitary system. The cortisol profiles reflect the combined effects of ACTH and cortisol on the adrenal gland. In women, cortisol concentrations were higher than in men, suggesting that the brain corticosteroid feedback was different between the sexes. When the parameter identification procedure was applied to these data, the parameter MRFormula 2 was found to be higher (2.3 ± 0.1 vs. 2.0 ± 0.1 fmol·l–1·s–1; P < 0.05) and GRFormula 2 lower (2.8 ± 0.1 vs. 3.2 ± 0.1 fmol·l–1·s–1; P < 0.05) in healthy women than it was in healthy men (Table 2). To interpret these results easier, we calculated the ratio between MRFormula 2 and GRFormula 2 as an estimate of the MR:GR balance; the balance was larger in healthy women (1:1.22) compared with healthy men (1:1.60). Note that MRFormula 2 quantifies a positive and GRFormula 2 a negative feedback (see Eq. 1).


View this table:
[in this window]
[in a new window]

 
Table 2. Results from the parameter identification method

 
When the identified parameters were used to reconstruct the ACTH and cortisol profiles, these theoretic predictions were similar to the observed data (prediction error < 0.05; Fig. 3A, right; Table 2). When analyzing how the variability of estimated parameters affected the predictions, we found that changes in the central parameters MRFormula 2GRFormula 2 and CRHFormula 2 sensitively influenced the predicted ACTH and cortisol concentrations (Fig. 4, A and B), while changes in the peripheral parameters ACTHFormula 2 and GRFormula 2 mainly affected predicted serum cortisol (Fig. 4, CE). Systematic variation of MRFormula 2 and GRFormula 2 revealed a high robustness of the prediction since there is only one parameter set representing the best model fit (Fig. 5).


Figure 4
Figure 4
View larger version (31K):
[in this window]
[in a new window]

 
Fig. 4. Effect of parameter-modifications on the prediction of ACTH and cortisol concentrations in healthy men. For validation of a model, it is mandatory to perform a sensitivity analysis (33, 34). Results of our sensitivity analysis are illustrated.: A: stepwise increasing the values of the central parameter MRE maxbrain by 125, 150, and 175%, leads to marked elevations of predicted ACTH and cortisol profiles. Thus, the predictions are very sensitive to changes in this parameter. B: using stepwise increased values of the central parameter GRE maxbrain, the predicted profiles of ACTH and cortisol are clearly lowered, indicating sensitivity to changes in this parameter. C: when the parameter CRHexogE is altered, we also find the predicted ACTH and cortisol profiles sensitively influenced. D and E: modifying the peripheral parameters ACTHEadrenal and GREadrenal mainly affects the predicted cortisol profiles in a sensitive manner.

 

Figure 5
View larger version (99K):
[in this window]
[in a new window]

 
Fig. 5. Effect of two-parameter-modifications on the precision of prediction in healthy men. Precision is determined by the deviation between two predictions: 1) the optimal prediction of data obtained in healthy men, 2) the prediction calculated with the modified parameters MRE maxbrain and GRE maxbrain; high precision (red), low precision (blue). If the central parameters MRE maxbrain and GRE maxbrain are slightly modified, predictions markedly deviate from the optimal predictions in healthy men. Note that there is only one parameter set of MRE maxbrain and GRE maxbrain fulfilling the best model fit (dark red).

 
When the identified parameters were used to analyze the dose-response function for the brain corticosteroid feedback, women had a weaker inhibitory GR effect than men (inset in Fig. 3A). In both sexes, the maximum MR-mediated positive feedback effects were evident at a low serum cortisol concentration of 50 nmol/l (typical value at bedtime), whereas at 300 nmol/l (typical value at wake-up time) the counteracting MR- and GR effects were of equal strength, and at 600 nmol/l (typical value during a stressful event), the negative feedback effect of GR predominated. The theoretically determined weaker central feedback inhibition in women reflects the observation that they had higher cortisol compared with men. In this experiment, we found that a model of the HPA system that makes use of positive and negative feedback loops is consistent with the data set.

Pharmacological blockade in healthy humans. We analyzed a data set derived from experiments using the MR antagonist canrenoate (96). When canrenoate was administered intravenously, serum cortisol concentrations were elevated compared with placebo (Fig. 3C, left). Similar results have been obtained in several other studies (3, 31, 35, 47, 99). Looking at the adrenal-to-brain feedback pathway, the fact that higher cortisol concentrations were observed suggests a different MR-GR feedback in the brain. Since canrenoate competes as a competitive antagonist with cortisol for binding to MR and (nonselectively to) GR, we identified the following parameters under canrenoate treatment: the affinities (described by MREC50 and GREC50) but not the efficacies (MRFormula 2 and GRFormula 2(Table 1). When the parameter identification procedure was applied to the experimental data, canrenoate was found to increase the MREC50 by around 7- fold (P < 0.005), but the GREC50 only by around 5-fold (Table 2).

When the identified parameters were used to reconstruct the ACTH and cortisol profiles, these theoretic predictions were similar to the observed data (prediction error < 0.05; Fig. 3C, right, Table 2). When the identified parameters were used to analyze the brain corticosteroid feedback, the pharmacological MR blockade was plausibly found to shift the dose-response curve to the right (Fig. 3C, inset). Thus, we used the computational approach in the pharmacologically blocked HPA system and found results (i.e., increased EC50 values for the corticoid receptors) that are in agreement with what is known about the characteristics of competitive blockade.

To interpret the influences of canrenoate, we refer to the three rules (Box 2). Rule 1: In the presence of competitive blockade (which affects MR and nonselectively also GR), cortisol binds only at medium concentrations to the brain MR and at very high concentrations to the brain GR (Fig. 6, top). Rule 2: Activated brain MR and GR operate in opposing manners. Rule 3: As long as brain MR are predominantly activated, a stimulatory effect prevails that acts to increase ACTH and consequently cortisol concentrations; as long as brain GR are predominantly activated, there is an inhibitory effect working that acts to decrease ACTH and so cortisol concentrations. In this way, the cortisol concentrations strive to a point of equilibrium at which the stimulatory and the inhibitory effects of the brain corticosteroid receptors are of the same magnitude (Fig. 6, bottom). These interactions would result in a homeostatic hypercortisolemic state under competitive blockade with canrenoate. On the basis of this reasoning, we can interpret that a blockade of the corticosteroid receptors elevates the brain setpoint of the HPA system (Fig. 3C inset); when cortisol has to compete with canrenoate for binding at the MR, it can only exert its full (MR-meditated) stimulatory drive in the HPA system at higher cortisol concentrations.


Figure 6
View larger version (19K):
[in this window]
[in a new window]

 
Fig. 6. Top: competitive blockade of the brain corticosteroid receptors and the right-shift of dose-response curves (horizontal arrows). Canrenoate predominantly affects the MR dose-response curve but also the GR dose-response curve in a nonselective manner (normal binding = solid lines; competitive blockade = dashed lines). In this way, the intersection (crossed circle) of the MR and GR receptor-kinetic curves is shifted to higher cortisol concentrations. Bottom: The summation curve of stimulatory MR and inhibitory GR feedback effects is also shifted to the right under competitive blockade with canrenoate. Note that canrenoate shifts the MR dose-response curve closer to the GR dose-response curve and in this way reduces the maximum of the biphasic feedback curve. The zero of the summation curve constitutes the brain setpoint of the HPA system (vertical arrows). At this setpoint, the brain's stimulatory and inhibitory effects on cortisol secretion are of the same magnitude, and the HPA system stabilizes in a homeostatic state. Canrenoate can be viewed to shift the brain setpoint of the HPA system to higher cortisol concentrations.

 
Looking in turn at the brain-to-adrenal forward pathway, the ACTH time courses were similar, but the cortisol concentrations were higher under canrenoate. We found canrenoate administration to decrease GRFormula 2 that is, efficacy of adrenal GR feedback. This finding is in line with the view of other researches suggesting that MR blockade affects the sensitivity of the adrenal gland to ACTH (99).

Genetic MR inactivation within the central nervous system represents another approach for testing our model of homeostasis. Central MR mutant mice show lower serum corticosteroid concentrations than their control littermates (10). We applied these published data in the parameter identification procedure and found a 20% (albeit not complete) reduction in the MRFormula 2 Merely, it remained unclear why the functional loss of MR efficacy was incomplete. Basically, our model of the HPA system also applies to this data set, in which the MR was genetically inactivated in the brain.

Experiments in patients with Addison's and Cushing's disease. We resorted to a data set that we had already published on patients with primary adrenal insufficiency, that is, Addison's disease, and on patients with central Cushing's disease, whereby the latter had had adrenalectomy (ADX) for therapeutic purposes (38). This experimental condition now made it possible to systematically investigate the brain corticosteroid feedback on ACTH almost over the entire range of cortisol concentrations, particularly focusing on the low-concentration range. When hydrocortisone was infused intravenously, the plasma ACTH in patients with Cushing's disease/ADX was found to steeply rise in the beginning and then to continuously fall (Fig. 3D, left), indicating both positive and negative feedback. Recently, we performed a replication study using up-to-date laboratory methods, and confirmed the stimulatory effect on ACTH at low cortisol concentrations (see supplemental Fig. S1 online). In patients with Addison's disease, low cortisol is the result of adrenal insufficiency. However, low cortisol, as it occurred after replacement therapy was withheld for 24 h, may also lead to secondary alterations in the MR:GR balance, as can be expected from experiments in adrenalectomized animals (69). We found that plasma ACTH continuously declined from the beginning of hydrocortisone infusion, thus indicating a dominant negative feedback effect of cortisol. Of note, in a replication study on patients with Addison's disease, we were able to detect an initial dose-dependent rise in plasma ACTH (albeit less pronounced than in patients Cushing's disease) (37). A similar initial and transient rise in plasma ACTH has also been reported in normal humans, when hydrocortisone was infused during oral metyrapone treatment (97). With Cushing's disease, the central MR efficacy is large compared with the MR efficacy seen in patients with Addison's disease/ADX (23.4 ± 6.3 vs. 0.3 ± 0.2 fmol·l–1·s–1; P < 0.005; Table 2). Moreover, in patients with Cushing's disease the "MR:GR balance" was markedly larger (1:1.09) than in normal men (1:1.60).

When the identified parameters were used to reconstruct the ACTH and cortisol profiles, the agreement between experimental data and theoretical prediction was satisfactory (prediction error <0.20; Fig. 3D, right, Table 2). When the identified parameters were used to analyze the brain corticosteroid feedback, central feedback effects over the entire cortisol range were found to be mainly inhibitory in patients with Addison's disease and stimulatory in patients with Cushing's disease (inset in Fig. 3D). Our model, therefore, adequately explains why the homeostatic level of cortisol is pathologically raised in Cushing's disease.

Since some scientists have referred to experiments (3, 31, 35, 47, 73, 81, 89, 99) and concluded that the MR exerted a tonic inhibition at low cortisol concentrations, we tested the alternative system of differential equations that described an inhibitory feedback effect of MR, in addition to the inhibitory feedback of GR. When this alternative system was used to reconstruct the ACTH and cortisol profiles in Cushing's disease/ADX, the observed peak ACTH concentration was found to deviate unacceptably from the prediction (prediction error > 0.60; see supplementary Fig. S2 online). Thus, our model had a wider scope of validity with the data sets examined here than a concept based on a MR tonic inhibition (Table 3).


View this table:
[in this window]
[in a new window]

 
Table 3. Agreement, between two alternative models and different data sets

 
Experiments in obese humans. Is the brain corticosteroid-feedback altered in obesity? When CRH was administered to obese male subjects (characteristics found in supplementary Table 1 online), ACTH responses were twofold higher, and the cortisol responses were slightly higher than in normal weight men (Fig. 3B, left). The obese men revealed a higher MR efficacy (2.4 ± 0.1 vs. 2.0 ± 0.1 fmol·l–1·s–1; P < 0.01) and a lower GR efficacy (2.8 ± 0.2 vs. 3.2 ± 0.1 fmol·l–1·s–1; P < 0.01) compared with the normal controls. In the obese men the "MR:GR balance" was larger (1:1.17) compared with healthy men (1:1.60). The peripheral adrenal parameters were of the same magnitude in the obese and control individuals.

The theoretical prediction shows a high agreement with the experimental observations (prediction error < 0.05; Fig. 3B right, Table 2). In obese men, central feedback inhibition was markedly reduced, especially at the high cortisol concentrations (inset in Fig. 3B). The computational approach allows locating specific derangements in the HPA system in all situations that are characterized by disturbed ACTH/cortisol secretion.

One homeostatic system often conveys an input signal to another related homeostatic system in the organism. For example, energy homoeostasis is closely linked to the HPA system (78). Our analysis shows that with obese men, the function of the GR is decreased and that as a result, the HPA system is less well-braked, that is, inhibited. The two stress systems, that is, the sympathetic nervous system and the HPA system, increase the energy supply of the brain. Recent studies have shown that in adrenalectomized animals the components of the HPA system are highly stimulated and that an exogenous application of glucose can almost completely return the hyperfunctioning HPA system to normal activity (8, 56, 57). Glucose can therefore substitute for the GR-mediated braking function of the HPA system (27). In a challenged organism, the alteration of one homeostatic system (e.g., brain energy metabolism) may strongly affect the homeostasis of another system (e.g., the HPA). Thus, input signals from related homeostatic systems can be influential, but, as we have demonstrated in this paper, they are not required to generate homeostasis.


    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
The question addressed by this study was whether the presented "principle of homeostasis" (Box 2) is in agreement with data on the HPA system. The main finding of the study is that this model applies even for extreme conditions in the HPA system, like Addison's and Cushing's disease.

The discussion about homeostasis is closely linked to the question of how setpoints are generated in biological systems. In this paper, we showed that the HPA system can be described by a "pn-system", that is, one that is driven by positive (p) MR-mediated and negative (n) GR-mediated feedback loops. We had also investigated an alternative "nn-system" that is driven by negative MR and negative GR feedback loops. The pn-systems are able to stabilize themselves with the help of a positive feedback loop to achieve equilibrium at a given setpoint. In contrast, the nn-systems do require a positive signal input into their regulatory loop; otherwise, they strive to equilibrium at zero values. Basically, generation of a setpoint requires a positive signal that originates either from the inside (positive feedback loop) or from the outside of a closed-loop system (external reference input). In humans, the HPA system is the major systemic stress system that normally stabilizes at a setpoint about midnight (ACTH 2 pmol/l, cortisol 50 nmol/l). It is important to maintain the low serum concentrations of cortisol during the slow wave sleep periods in the early night, because cortisol then predominantly binds to hippocampal MR, which is necessary for memory consolidation, and underlying synaptic plasticity in these neuronal networks (95). If the HPA system were an nn-system, it would require a steady positive signal input from outside the system (e.g., from the isocortex) to maintain midnight concentrations of ACTH and cortisol. If such a second system also were an nn-system, a third system would be necessary, and so on. Adding further nn-systems does not solve the problem of how homeostatic setpoints are created. From this theoretical point of view, homeostasis is possible only when at least one pn-system exists in the organism. Our data support the view that the HPA system itself is a required pn-system that allows achieving homeostasis at a given setpoint.

The problem that one experimental data set can be in agreement with two distinct theoretical approaches has often prompted a scientific debate. This issue is often referred to as "inverse problem" (87). Some researchers consider the HPA systems to be characterized by an nn-system, while we present the pn-system here. During the past decades, extensive empirical data accumulated that led researchers to assume an nn-system, which includes MR-mediated tonic inhibition, that is, negative feedback. The majority of these data has been obtained from experiments using more or less specific pharmacological MR blockade in humans (3, 31, 35, 47, 99) and animals (73, 81, 89). As a representative of such studies, we selected one from our own group for a more detailed mathematical investigation. In the original publication on the collected data (96), it was at first concluded that "MR blockade leads to a tonically increased cortisol secretion," suggesting an nn-system regulation of cortisol. When the novel computational approach was used to reinvestigate the data, an nn-model was, in fact, found to fit these experimental data, but in contrast to some expectations, the pn-model, which included MR-mediated positive feedback, also fit the experimental data (Table 3). Such an ambiguity typically arises from an "inverse problem". At this point, we emphasize that "MR-mediated tonic inhibition" is not necessarily the only explanation in this data set, but "MR-mediated positive feedback" is a plausible possibility. The above-mentioned data sets on canrenoate cannot help us to distinguish which of the alternative models is more accurate in making predictions. Neither can the data sets on CRH challenge tests in healthy or obese humans help to make this distinction (Table 3). To solve an "inverse problem," it might be helpful to investigate additional data sets. We therefore selected a data set that particularly focused on regulatory behavior of the HPA system in the low-cortisol concentration range in adrenalectomized patients. The comparison between the two alternative models shows that the pn-model fits with this data set, while an nn-model fails. Thus, we conclude that our model, which uses MR-mediated positive feedback components, has a wider scope of validity compared with a model with pure negative feedback (Table 3).

What are the underlying molecular mechanisms that could mediate positive and negative feedback in the HPA system? First, the brain mineralocorticoid and glucocorticoid receptors fulfill the key functions in the homeostatic rules (Box 2): binding cortisol with high and low affinities, acting in opposing manners, and mediating feedback effects on cortisol. The fact that the pair of MR and GR fulfills these required criteria supports MR and GR as one mechanism underling positive and negative feedback in the HPA system, but we cannot rule out that other molecules are involved. Second, the fact that MR and GR fulfill the criteria (Box 2) does not allow us to infer whether genomic and nongenomic molecular mechanisms underlie the HPA system's feedback. By genomic mechanisms, the corticosteroid receptors modify transcription of responsive genes, either through DNA binding or through protein-protein interactions with other transcription factors (7). Glucocorticoids influence de novo protein synthesis within a short, intermediate, and prolonged timeframe (41, 42). By nongenomic mechanisms, the glucocorticoids can exert fast effects on excitation-secretion coupling (48). Furthermore, MR were found indispensable for rapid nongenomic modulation of hippocampal glutamate transmission by corticosteroids (51). Corticosteroids also exert fast-onset inhibitory effects via membrane receptors (32, 93). There is the notion that the fast glucocorticoid actions, which are mediated by membrane receptors, are an ancient type of sterol/steroid-mediated effect and that these may be the primordial glucocorticoid receptors (26). It is conceivable that brain corticosteroid feedback consists of a combination of genomic and nongenomic MR and GR actions. Whereas nongenomic MR and GR mechanisms modulate the early response, the late response is mainly influenced by genomic mechanisms, with the latter being more important in the homeostatic phase. Of note, the rules underlying our model (Box 2) do not include auxiliary conditions about time. Most relevant for our considerations are the balances assessed in the homeostatic state of equilibrium. Third, MR and GR display heterogenic distribution in the brain, and they may be located on excitatory or inhibitory neurons. Thus, a diversity of specific neuronal pathways and molecular mechanisms that underlies the brain corticosteroid feedback is conceivable. In all, on the basis of the current evidence, we regard MR and GR as the major mechanisms that fulfill the functions in the homeostatic rules and that underlie brain corticosteroid feedback.

Is our model of the HPA system oversimplified? One point of view might be that important control mechanisms escape mention here so that the true complexity of the HPA system is not delved into. We have, in fact, refrained from including a large number of mechanisms that might be relevant in the regulation, for example, CBG, autoregulation of MR and GR and heterodimerization. First, CBG binds cortisol with high affinity but low capacity, so that the extent to which it binds cortisol saturates with midlevel cortisol concentrations. In this way, CBG acts like a functional buffer. When we tentatively included CBG in the model, we found that CBG had minor influence on the binding kinetics of MR and GR (particularity with rising or falling cortisol), but the crucial characteristics of two intersecting receptor-kinetic curves (Fig. 1D, top) were preserved. Second, autoregulation of MR and GR describes a short loop feedback process, by which both cerebral receptors influence their own gene expression. This process affects the efficacies of MR and GR. When we included autoregulation in the model, we observed that the system became more stable, but the essential features of two intersecting receptor-kinetic curves (Fig. 1D, top) remained conserved. Third, heterodimerization of MR and GR mainly affects how the separate receptor effects are integrated (61, 72, 88). When heterodimerization was included in a preliminary model (22), the combined feedback effect of both receptors effects (Fig. 1D, bottom) displayed an enhanced contrast between the positive and negative feedback effects, but again the critical aspects of the regulation remained unaffected. In all, there are several regulatory processes in the HPA system involved that might, of course, modify its regulatory behavior. Since using the presented minimal model of the HPA system led to satisfactory approximation in all of our data sets, we abstained from including extra biological mechanisms into the final model. We are aware that the restrictions proposed here might be the subject of debate, but we feel that the specificities of the model presented are less important than the general basic conception proposed here for homeostasis.

Our findings on the HPA system regulation lead us to speculate that the outlined principle also applies to other biological systems. Almost all known ligands bind to at least two receptors (a selection is shown in Table 4). In most of these dual receptor systems, the first and second rule (Box 2) has been demonstrated (17). In many dual-receptor systems, also the third rule has been established. These systems are the homeostatic ones. As a first example, the brain-derived neurotrophic factor (BDNF) binds to high- and low-affinity receptors and exerts counteracting effects on the survival and apoptosis of cells (Table 4). BDNF acts as an autocrine survival factor in neurons supporting its own release at low BDNF concentrations and inhibiting it at high concentrations. Thus, there is experimental evidence that all three homeostatic rules also apply to BDNF. As a second example, insulin has been shown to act by a dual receptor system with evidence for the existence of all three homeostatic rules (Table 4). The pancreatic beta-cell secretes insulin molecules, which in an autocrine manner bind to high- or low-affinity insulin receptors located on the beta-cell surface. Insulin acts on its own secretion either stimulatory at low insulin concentrations or inhibitory at high concentrations. In this way, insulin can stabilize its basal pancreatic secretion in a self-perpetuating process. As a third example, the three rules may also apply to the regulation of neuronal ATP (79). High- and low-affinity ATP-sensitive potassium channels adjust the excitation of neurons, the neuronal energy provision, and neuronal energy consumption. Here, dual ATP-sensitive receptor systems exert a homoeostatic influence on cellular energy metabolism. In summary, we can, therefore, find evidence for examples from most receptor classes, which suggest that many biological systems have evolved positive feedback loops to create homeostasis.


View this table:
[in this window]
[in a new window]

 
Table 4. Examples of biological systems and respective evidence from the the literature that the rules of homeostasis apply

 
In this paper, we did not intend to show that the outlined model of the HPA-system regulation represents reality or that any other possibilities are excluded. We did demonstrate, however, that this model could fit all of our experimental data in a sensitive and robust manner and that among alternative models, the "principle of homeostasis" was in agreement with the largest data set, representing the widest scope of validity.


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This work was supported by grants (Clinical Research Group KFO-126) from the German Research Foundation.


    ACKNOWLEDGMENTS
 
We are grateful to M. Hallschmid (Institute for Neuroendocrinology, Luebeck), D. Langemann (Institute for Mathematics, Luebeck) and K. Oltmanns (Department of Psychiatry and Psychotherapy, Luebeck) for critically reading the manuscript and giving their most constructive comments. We thank P. Wellhöner and C. Dodt (Department of Internal Medicine I, Luebeck) for their contributions (mineralocorticoid receptor blockade study).


    FOOTNOTES
 

Address for reprint requests and other correspondence: A. Peters, Dept. of Internal Medicine I, Endocrinology, Univ. of Luebeck, 23538 Luebeck, Germany (e-mail: achim.peters{at}uni-luebeck.de)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 

  1. Acheson A, Conover JC, Fandl JP, DeChiara TM, Russell M, Thadani A, Squinto SP, Yancopoulos GD, Lindsay RM. A BDNF autocrine loop in adult sensory neurons prevents cell death. Nature 374: 450–453, 1995.[CrossRef][Medline]
  2. Angeli D, Ferrell JE Jr., Sontag ED. Detection of multistability, bifurcations, and hysteresis in a large class of biological positive-feedback systems. Proc Natl Acad Sci USA 101: 1822–1827, 2004.[Abstract/Free Full Text]
  3. Arvat E, Maccagno B, Giordano R, Pellegrino M, Broglio F, Gianotti L, Maccario M, Camanni F, Ghigo E. Mineralocorticoid receptor blockade by canrenoate increases both spontaneous and stimulated adrenal function in humans. J Clin Endocrinol Metab 86: 3176–3181, 2001.[Abstract/Free Full Text]
  4. Avital A, Segal M, Richter-Levin G. Contrasting roles of corticosteroid receptors in hippocampal plasticity. J Neurosci 26: 9130–9134, 2006.[Abstract/Free Full Text]
  5. Bahrenberg G, Behrmann I, Barthel A, Hekerman P, Heinrich PC, Joost HG, Becker W. Identification of the critical sequence elements in the cytoplasmic domain of leptin receptor isoforms required for Janus kinase/signal transducer and activator of transcription activation by receptor heterodimers. Mol Endocrinol 16: 859–872, 2002.[Abstract/Free Full Text]
  6. Barker PA. p75NTR: A study in contrasts. Cell Death Differ 5: 346–356, 1998.[CrossRef][ISI][Medline]
  7. Beato M, Herrlich P, Schutz G. Steroid hormone receptors: many actors in search of a plot. Cell 83: 851–857, 1995.[CrossRef][ISI][Medline]
  8. Bell ME, Bhatnagar S, Liang J, Soriano L, Nagy TR, Dallman MF. Voluntary sucrose ingestion, like corticosterone replacement, prevents the metabolic deficits of adrenalectomy. J Neuroendocrinol 12: 461–470, 2000.[CrossRef][ISI][Medline]
  9. Bellen A, Zennaro M. Numerical Methods for Delay Differential Equations. New York: Oxford University, 2003.
  10. Berger S, Wolfer DP, Selbach O, Alter H, Erdmann G, Reichardt HM, Chepkova AN, Welzl H, Haas HL, Lipp HP, Schutz G. Loss of the limbic mineralocorticoid receptor impairs behavioral plasticity. Proc Natl Acad Sci USA 103: 195–200, 2006.[Abstract/Free Full Text]
  11. Bibel M, Hoppe E, Barde YA. Biochemical and functional interactions between the neurotrophin receptors trk and p75NTR. EMBO J 18: 616–622, 1999.[CrossRef][ISI][Medline]
  12. Born J, DeKloet ER, Wenz H, Kern W, Fehm HL. Sleep-influence of mr/gr- corticosteronreceptor. Am J Physiol Endocrinol Metab 260: E183–E188, 1991.[Abstract/Free Full Text]
  13. Brandman O, Ferrell JE, Li R, Meyer T. Interlinked fast and slow positive feedback loops drive reliable cell decisions. Science 310: 496–498, 2005.[Abstract/Free Full Text]
  14. Bridgham JT, Carroll SM, Thornton JW. Evolution of hormone-receptor complexity by molecular exploitation. Science 312: 97–101, 2006.[Abstract/Free Full Text]
  15. Buccafusco JJ, Brezenoff HE. Opposing influences on behavior mediated by muscarinic and nicotinic receptors in the rat posterior hypothalamic nucleus. Psychopharmacology (Berl) 67: 249–254, 1980.[CrossRef][Medline]
  16. Calabrese EJ. Adrenergic receptors: biphasic dose responses. Crit Rev Toxicol 31: 523–538, 2001.[CrossRef][ISI][Medline]
  17. Calabrese EJ, Baldwin LA. Peptides and hormesis. Crit Rev Toxicol 33: 355–405, 2003.[ISI][Medline]
  18. Cannon WB. The Wisdom of the Body. New York: W. W. Norton, 1932.
  19. Chao MV, Hempstead BL. p75 and Trk: a two-receptor system. Trends Neurosci 18: 321–326, 1995.[CrossRef][ISI][Medline]
  20. Chen SC, Kochan JP, Campfield LA, Burn P, Smeyne RJ. Splice variants of the OB receptor gene are differentially expressed in brain and peripheral tissues of mice. J Recept Signal Transduct Res 19: 245–266, 1999.[ISI][Medline]
  21. Chow YS, Teicher H. Probability Theory: Independence Interchangeability, Martingales. New York: Springer, 1997.
  22. Conrad M, Hubold C, Fischer B, Peters A. The selfish brain: a new model of the LHPA system. Exp Clin Endocrinol Diabetes 112: 481–482. 2004.
  23. Croce A, Astier H, Recasens M, Vignes M. Opposite effects of alpha 1- and beta-adrenoceptor stimulation on both glutamate- and gamma-aminobutyric acid-mediated spontaneous transmission in cultured rat hippocampal neurons. J Neurosci Res 71: 516–525, 2003.[CrossRef][ISI][Medline]
  24. Crochemore C, Lu J, Wu Y, Liposits Z, Sousa N, Holsboer F, Almeida OF. Direct targeting of hippocampal neurons for apoptosis by glucocorticoids is reversible by mineralocorticoid receptor activation. Mol Psychiatry 10: 790–798, 2005.[CrossRef][ISI][Medline]
  25. Crochemore C, Lu J, Wu Y, Liposits Z, Sousa N, Holsboer F, Almeida OF. Direct targeting of hippocampal neurons for apoptosis by glucocorticoids is reversible by mineralocorticoid receptor activation. Mol Psychiatry 10: 790–798, 2005.[CrossRef][ISI][Medline]
  26. Dallman MF. Fast glucocorticoid actions on brain: Back to the future. Front Neuroendocrinol 26: 103–108, 2005.[CrossRef][ISI][Medline]
  27. Dallman MF, Pecoraro N, Akana SF, la Fleur SE, Gomez F, Houshyar H, Bell ME, Bhatnagar S, Laugero KD, Manalo S. Chronic stress and obesity: a new view of "comfort food". Proc Natl Acad Sci USA 100: 11696–11701, 2003.[Abstract/Free Full Text]
  28. de Kloet ER, Vreugdenhil E, Oitzl MS, Joels M. Brain corticosteroid receptor balance in health and disease. Endocr Rev 19: 269–301, 1998.[Abstract/Free Full Text]
  29. De Meyts P, Wallach B, Christoffersen CT, Urso B, Gronskov K, Latus LJ, Yakushiji F, Ilondo MM, Shymko RM. The insulin-like growth factor-I receptor. Structure, ligand-binding mechanism and signal transduction. Horm Res 42: 152–169, 1994.[ISI][Medline]
  30. Derendorf H, Mollmann H, Barth J, Mollmann C, Tunn S, Krieg M. Pharmacokinetics and oral bioavailability of hydrocortisone. J Clin Pharmacol 31: 473–476, 1991.[Abstract]
  31. Deuschle M, Weber B, Colla M, Muller M, Kniest A, Heuser I. Mineralocorticoid receptor also modulates basal activity of hypothalamus-pituitary-adrenocortical system in humans. Neuroendocrinology 68: 355–360, 1998.[CrossRef][ISI][Medline]
  32. Di S, Malcher-Lopes R, Halmos KC, Tasker JG. Nongenomic glucocorticoid inhibition via endocannabinoid release in the hypothalamus: a fast feedback mechanism. J Neurosci 23: 4850–4857, 2003.[Abstract/Free Full Text]
  33. Distefano JJ, III. A retrospective on the modeling methodology forum. The review process for modeling papers: a revisionist's experiences and perspectives. Am J Physiol Endocrinol Metab 267: E485–E488, 1994.[Free Full Text]
  34. DiStefano JJ, Stubberud AR, Williams IJ. Theory and Problems of Feedback and Control Systems. New York: McGraw-Hill, 1997.
  35. Dodt C, Kern W, Fehm HL, Born J. Antimineralocorticoid canrenoate enhances secretory activity of the hypothalamus-pituitary-adrenocortical (HPA) axis in humans. Neuroendocrinology 58: 570–574, 1993.[ISI][Medline]
  36. During MJ, Leone P, Davis KE, Kerr D, Sherwin RS. Glucose modulates rat substantia nigra GABA release in vivo via ATP-sensitive potassium channels. J Clin Invest 95: 2403–2408, 1995.[ISI][Medline]
  37. Fehm HL, Voigt KH, Kummer G, Lang R, Pfeiffer EF. Differential and integral corticosteroid feedback effects on ACTH secretion in hypoadrenocorticism. J Clin Invest 63: 247–253, 1979.[ISI][Medline]
  38. Fehm HL, Voigt KH, Lang RE, Beinert KE, Kummer GW, Pfeiffer EF. Paradoxical ACTH response to glucocorticoids in Cushing's disease. N Engl J Med 297: 904–907, 1977.[Abstract]
  39. Fletcher R. Practical Methods of Optimization. New York: Wiley, 1999.
  40. Fukuda N, Kubo A, Izumi Y, Soma M, Kanmatsuse K. Characteristics and expression of transforming growth factor-beta receptor subtypes on vascular smooth muscle cells from spontaneously hypertensive rats. J Hypertens 13: 831–837, 1995.[CrossRef][ISI][Medline]
  41. Ginsberg AB, Campeau S, Day HE, Spencer RL. Acute glucocorticoid pretreatment suppresses stress-induced hypothalamic-pituitary-adrenal axis hormone secretion and expression of corticotropin-releasing hormone hnRNA but does not affect c-fos mRNA or fos protein expression in the paraventricular nucleus of the hypothalamus. J Neuroendocrinol 15: 1075–1083, 2003.[CrossRef][ISI][Medline]
  42. Ginsberg AB, Frank MG, Francis AB, Rubin BA, O'Connor KA, Spencer RL. Specific and time-dependent effects of glucocorticoid receptor agonist RU28362 on stress-induced pro-opiomelanocortin hnRNA, c-fos mRNA and zif268 mRNA in the pituitary. J Neuroendocrinol 18: 129–138, 2006.[CrossRef][ISI][Medline]
  43. Goodman LV, Majack RA. Vascular smooth muscle cells express distinct transforming growth factor-beta receptor phenotypes as a function of cell density in culture. J Biol Chem 264: 5241–5244, 1989.[Abstract/Free Full Text]
  44. Goumans MJ, Valdimarsdottir G, Itoh S, Rosendahl A, Sideras P, Ten Dijke P. Balancing the activation state of the endothelium via two distinct TGF-beta type I receptors. EMBO J 21: 1743–1753, 2002.[CrossRef][ISI][Medline]
  45. Gudmundsson A, Carnes M. Pulsatile adrenocorticotropic hormone: an overview. Biol Psychiatry 41: 342–365, 1997.[CrossRef][ISI][Medline]
  46. Hairer E, Nørsett SP, Wanner G. Solving Ordinary Differential Equations I. New York: Springer Verlag, 1993.
  47. Heuser I, Deuschle M, Weber B, Stalla GK, Holsboer F. Increased activity of the hypothalamus-pituitary-adrenal system after treatment with the mineralocorticoid receptor antagonist spironolactone. Psychoneuroendocrinology 25: 513–518, 2000.[CrossRef][ISI][Medline]
  48. Hinz B, Hirschelmann R. Rapid non-genomic feedback effects of glucocorticoids on CRF-induced ACTH secretion in rats. Pharm Res 17: 1273–1277, 2000.[CrossRef][ISI][Medline]
  49. Inagaki N, Gonoi T, Clement JP, Wang CZ, Aguilar-Bryan L, Bryan J, Seino S. A family of sulfonylurea receptors determines the pharmacological properties of ATP-sensitive K+ channels. Neuron 16: 1011–1017, 1996.[CrossRef][ISI][Medline]
  50. Joels M, Hesen W, de Kloet ER. Long-term control of neuronal excitability by corticosteroid hormones. J Steroid Biochem Mol Biol 53: 315–323, 1995.[CrossRef][ISI][Medline]
  51. Karst H, Berger S, Turiault M, Tronche F, Schutz G, Joels M. Mineralocorticoid receptors are indispensable for nongenomic modulation of hippocampal glutamate transmission by corticosterone. Proc Natl Acad Sci USA 102: 19204–19207, 2005.[Abstract/Free Full Text]
  52. Kim JJ, Diamond DM. The stressed hippocampus, synaptic plasticity and lost memories. Nat Rev Neurosci 3: 453–462, 2002.[CrossRef][ISI][Medline]
  53. Korz V, Frey JU. Stress-related modulation of hippocampal long-term potentiation in rats: Involvement of adrenal steroid receptors. J Neurosci 23: 7281–7287, 2003.[Abstract/Free Full Text]
  54. Krozowski ZS, Funder JW. Renal mineralocorticoid receptors and hippocampal corticosterone-binding species have identical intrinsic steroid specificity. Proc Natl Acad Sci USA 80: 6056–6060, 1983.[Abstract/Free Full&n