Report

The Frequency Dependence of Osmo-Adaptation in Saccharomyces cerevisiae

See allHide authors and affiliations

Science  25 Jan 2008:
Vol. 319, Issue 5862, pp. 482-484
DOI: 10.1126/science.1151582

Abstract

The propagation of information through signaling cascades spans a wide range of time scales, including the rapid ligand-receptor interaction and the much slower response of downstream gene expression. To determine which dynamic range dominates a response, we used periodic stimuli to measure the frequency dependence of signal transduction in the osmo-adaptation pathway of Saccharomyces cerevisiae. We applied system identification methods to infer a concise predictive model. We found that the dynamics of the osmo-adaptation response are dominated by a fast-acting negative feedback through the kinase Hog1 that does not require protein synthesis. After large osmotic shocks, an additional, much slower, negative feedback through gene expression allows cells to respond faster to future stimuli.

The mechanisms cells use to sense and respond to environmental changes include complicated systems of biochemical reactions that occur with rates spanning a wide dynamic range. Reactions can be fast, such as association and dissociation between a ligand and its receptor (<1 s), or slow, such as protein synthesis (>103 s). Although a system may comprise hundreds of reactions, often only a few of them dictate the system dynamics. Unfortunately, identification of the dominant processes is often difficult, and many models instead incorporate knowledge of all reactions in the system. Although occasionally successful (14), this exhaustive approach often suffers from missing information, such as unknown interactions or parameters.

Here, we used systems-engineering tools to study how oscillatory signals propagate through a signal transduction cascade, which allowed us to identify and to model concisely the interactions that dominate system dynamics. The cornerstone of this approach is to measure the cascade output in response to input signals oscillating at a range of frequencies (5, 6). By comparing the frequency response of the wild-type network to that of mutants, the molecular underpinnings of network dynamics can be determined. Studies of neural and other physiological systems have used systems theory (6), and control theory has also been applied to cellular networks (714).

We focused on the high-osmolarity glycerol (HOG) mitogen-activated protein kinase (MAPK) cascade in the budding yeast Saccharomyces cerevisiae. This cascade forms a core module of the hyperosmotic shock–response system and is particularly well suited to frequency-response analysis for several reasons. First, both the input (extracellular osmolyte concentration) and output (activity of the MAPK Hog1) of the network are easily measured and manipulated. Second, the molecular components of the network have been well studied, which facilitates connecting dynamic models with molecular events. Finally, the system contains multiple negative-feedback loops that operate on different time scales (4, 15, 16). It is still unclear which negative-feedback loop or loops dominate the signaling dynamics and whether the different feedback loops have distinct biological functions. We determined the properties of the main negative-feedback loops in the HOG network and arrived at a concise predictive model of the signaling dynamics. Furthermore, by analyzing the system's dynamics over a range of osmotic-shock strengths, we begin to understand how the multiple-feedback architecture might be beneficial for osmotic homeostasis in fluctuating environments.

After a hyperosmotic shock, membrane proteins trigger a signal transduction cascade that culminates in the activation of the MAPK Hog1, which is primarily cytoplasmic before the osmoshock (17, 18). When activated, Hog1 accumulates in the nucleus (Fig. 1A), where it activates a broad transcriptional response to osmotic stress (19). Constitutively active phosphatases dephosphorylate and deactivate Hog1, which leads to its export from the nucleus. When osmotic balance is restored, through changes either to the extracellular environment or to the intracellular osmolyte concentration, cascade activity ceases, and the Hog1 nuclear enrichment decreases (Fig. 1A). To estimate the amount of phosphorylated Hog1 in living cells, we simultaneously monitored the cellular localization of Hog1-YFP, a yellow fluorescent protein fused to Hog1, and Nrd1-RFP, a red fluorescent protein fused to a strictly nuclear protein. To quantify Hog1 nuclear localization, we define the function, R(t) = (<YFP>nucleus/<YFP>cell)population, as the ratio (averaged over many cells) of mean YFP pixel intensities in the nucleus and the whole cell [(Fig. 1A), red circles].

Fig. 1.

Enrichment of Hog1 nuclear localization is driven by pulsed salt shocks. (A) Localization of the fusion protein Hog1-YFP and the nuclear marker Nrd1-RFP by fluorescence microscopy. We applied and removed NaCl (0.2 M) as shown by the blue line. The population average translocation response (red circles) was defined by the ratio of average YFP fluorescence in the nucleus to the average whole-cell YFP fluorescence. (B) Oscillations of Hog1-YFP translocation in a population of cells (red circles) in response to square-wave oscillations in the input of extracellular NaCl (blue line).

Cells were periodically shocked in a flow chamber (fig. S1) in which a computer-actuated valve supplied square-wave pulses of medium, with and without 0.2 M NaCl [(Fig. 1B), blue line]. Localization of fluorescent proteins was concurrently imaged [(Fig. 1B), red circles]. Using Fourier analysis (20), we approximated both the input and output signals as sine waves oscillating with a period T0 = 2π/ω and quantified the output signal by the magnitude of the amplitude (A) at the driving frequency A(ω) and the relative phase shift φ(ω) (Fig. 2A) (21). We measured the response of the system to square-wave stimuli with periods ranging from T0 = 2 min to T0 = 128 min (fig. S2). These responses were analyzed to obtain so-called Bode plots (Fig. 2, B and C) (20), representing the frequency-dependent amplitude A(ω) and phase φ(ω).

Fig. 2.

Fourier analysis, model fits, and model predictions of Hog1 nuclear enrichment. (A) Illustration of the input (NaCl concentration, blue line), the network response (Hog1-YFP translocation, red circles), and the sine wave (black line) corresponding to the Fourier component of the response at the driving frequency ω. This Fourier component is described by three parameters: A(ω) (green) representing the amplitude of the oscillations, φ(ω) (brown) representing the phase delay between the input and the response oscillations, and y0 (ω) representing the signal offset. (B) Measurement of the Fourier amplitude A(ω) (green dots) over a range of driving frequencies along with model fit (green line). (C) Phase of the response measured relative to the driving signal (brown dots) along with model prediction of the phase (brown line). (D) Response of the system to a step increase of 0.2 M NaCl compared with the step response predicted by the model. The “low Pbs2” data (boxes) are gathered from the Pbs2 underexpressing mutant strain and were used to generate the model fit [gray line in (B)] and model predictions [gray lines in (C) and (D)].

We used linear-systems theory to develop a predictive model for the response to arbitrary osmotic input signals u(t) (20). We fitted a general second-order linear time–invariant (LTI) model to the data in Fig. 2B and used the extracted parameters and a simple nonlinear element (20, 22) to predict the response to a step input of 0.2 M NaCl. The model accurately predicted both the response amplitude and the time required to return to basal activity [(Fig. 2D), blue circles].

Because our model was not instructed by knowledge of the underlying biology, we sought to explore how it is similar to and different from the canonical molecular model of the hyperosmotic-shock response. Thus, we converted our LTI model into a model that is more readily interpreted in terms of biological processes (20): Math(1) Math This model contains two negative-feedback loops, which act to reduce the difference, (A0ux), between the stimulus A0u and the internal-state variable x. Fig. 3A shows a schematic of the model, and Fig. 3B shows the canonical biochemical mechanisms that regulate the osmotic-shock response. Because of the high degree of correspondence between our derived model and the extensive Hog1 literature, we interpret the state variable, x, as the intracellular osmolyte concentration and y as enrichment of phosphorylated Hog1 above its baseline level. Thus, interpreted in biological terms, our model predicts that one feedback pathway depends on Hog1-induced glycerol accumulation (i.e., βy changes x through the activity of the observable y), whereas a second glycerol-accumulating pathway is Hog1-independent [i.e., α(A0ux) changes x independently of y].

Fig. 3.

Network topology implied by pulsed-input analysis corresponds to biological network. (A) Diagrammatic representation of the mechanistic model shows two linear negative-feedback responses: one dependent on Hog1 activity with strength β and a second independent of Hog1 activity with strength α. The intracellular osmotic pressure was modeled as an integrator, whereas the MAPK signal transduction pathway was modeled by the linear impulse–response function. The output of circumscribed-plus-sign symbols is simply the sum of its inputs. (B) The osmo-adaptation network structure. With the inducement of osmotic stress, cells increase or decrease their export rate of glycerol through the transmembrane protein Fps1, which is modified by Hog1-independent and Hog1-dependent mechanisms. In addition, under high osmotic stresses, active nuclear Hog1 is known to modify the expression of glycerol-producing proteins over longer time scales. ΔP represents the difference between internal and external pressure relative to its optimal value.

To gauge the relative strength of the two feedbacks, we applied the same Fourier technique to a mutant strain expressing a reduced amount of Pbs2, the MAPK kinase that phosphorylates Hog1 (Fig. 3B). Because Hog1 is not activated as highly after a hyperosmotic shock in this strain, we can effectively isolate Hog1-dependent feedback from Hog1-independent feedback by comparing the response dynamics (fig. S4) to those of the wild-type strain. Bode plots for this strain were different from those of the wild-type strain (Fig. 2, B and C). Matching the LTI model to the A(ω) data alone yielded a good fit [(Fig. 2B), black squares], as well as an accurate prediction for the φ(ω) data [(Fig. 2C), black squares] and the step response [(Fig. 2D), boxes]. The latter prediction [(Fig. 2D), gray line] reproduced both the reduced maximum response and the slower response dynamics observed in this strain. This suggests that the Hog1-dependent feedback loop plays a major role in rapidly regulating the osmotic-shock response (20).

We compared our model with known biological details of the hyperosmotic-shock response (4, 15, 16). Yeast regulate their intracellular osmolyte concentration through two parallel mechanisms. In a Hog1-independent manner, the membrane protein Fps1 quickly (<2 min) responds by decreasing the glycerol-export rate (23, 24), thereby leading to glycerol retention. Further, active Hog1 increases expression of the glycerol-producing proteins Gpd1 and Gpp2. This raises the intracellular glycerol level over longer time scales (>30 min) (25).

Although the topology of our derived model corresponds closely to that of the known biological system (Fig. 3, A and B), dynamic differences suggest that the current view of the MAPK's role in osmotic regulation is incomplete. Cells begin to recover from the NaCl pulse within 5 min and are finished responding within 15 min. Both of these time scales are faster than required for gene expression, which is typically greater than 15 min (4). This suggests that both feedback loops in our model control the rapid accumulation of glycerol, consistent with previous reports (23, 24, 26).

Our model suggests that gene expression plays a minimal role in the hyperosmotic-shock response, yet the expression of hundreds of genes changes in response to hyperosmotic shock. We hypothesized that gene expression may be more important as a longer–time scale feedback in this system, so we looked for pulse-to-pulse variability in the response of cells stimulated with periodic pulses of NaCl (fig. S5). Cells were shocked either in the absence (Fig. 4, A to D) or presence (Fig. 4, E to H) of cycloheximide, a small molecule that inhibits protein synthesis. As predicted by the initial data, cells responded very similarly to an initial pulse of osmolyte regardless of their ability to synthesize new proteins (Fig. 4, black lines). Nevertheless, we found that cells stimulated multiple times recovered from each subsequent pulse faster in the absence of cycloheximide and slower in its presence, revealing a longer–time scale component absent from our earlier data. These results suggest that nontranscriptional feedback mediates short–time scale osmolyte accumulation (4, 16, 26, 27), whereas gene expression plays a role in osmolyte production only on longer time scales and for more intense shocks. Accordingly, we found that stronger shocks cause cells to increase their rate of glycerol production (fig. S6) in a manner that depends on gene expression, which permits faster recovery to subsequent fluctuations (28).

Fig. 4.

Gene expression facilitates response to subsequent pulses. Consecutive-pulse responses were compared for cells treated with (A) 0.1 M NaCl, 16-min period (i.e., 8 min at 0.1 M followed by 8 min at 0.0 M); (B) 0.2 M NaCl, 32-min period; (C) 0.35 M NaCl, 45-min period; and (D) 0.5 M NaCl, 60-min period. (E to H) Similarly treated cells also exposed to cycloheximide.

These results demonstrate the promise of applying engineering principles to cellular networks, particularly when predicting the response of the system to dynamic stimuli. In more complex systems, measuring the activity level of all relevant state-space variables could help with determining the effective network structure.

Supporting Online Material

www.sciencemag.org/cgi/content/full/319/5862/482/DC1

Materials and Methods

SOM Text

Figs. S1 to S7

Table S1

References

References and Notes

View Abstract

Stay Connected to Science

Navigate This Article