A systems biology analysis of long and short-term memories of osmotic stress adaptation in fungi
© You et al.; licensee BioMed Central Ltd. 2012
Received: 26 October 2011
Accepted: 25 May 2012
Published: 25 May 2012
Saccharomyces cerevisiae senses hyperosmotic conditions via the HOG signaling network that activates the stress-activated protein kinase, Hog1, and modulates metabolic fluxes and gene expression to generate appropriate adaptive responses. The integral control mechanism by which Hog1 modulates glycerol production remains uncharacterized. An additional Hog1-independent mechanism retains intracellular glycerol for adaptation. Candida albicans also adapts to hyperosmolarity via a HOG signaling network. However, it remains unknown whether Hog1 exerts integral or proportional control over glycerol production in C. albicans.
We combined modeling and experimental approaches to study osmotic stress responses in S. cerevisiae and C. albicans. We propose a simple ordinary differential equation (ODE) model that highlights the integral control that Hog1 exerts over glycerol biosynthesis in these species. If integral control arises from a separation of time scales (i.e. rapid HOG activation of glycerol production capacity which decays slowly under hyperosmotic conditions), then the model predicts that glycerol production rates elevate upon adaptation to a first stress and this makes the cell adapts faster to a second hyperosmotic stress. It appears as if the cell is able to remember the stress history that is longer than the timescale of signal transduction. This is termed the long-term stress memory. Our experimental data verify this. Like S. cerevisiae, C. albicans mimimizes glycerol efflux during adaptation to hyperosmolarity. Also, transient activation of intermediate kinases in the HOG pathway results in a short-term memory in the signaling pathway. This determines the amplitude of Hog1 phosphorylation under a periodic sequence of stress and non-stressed intervals. Our model suggests that the long-term memory also affects the way a cell responds to periodic stress conditions. Hence, during osmohomeostasis, short-term memory is dependent upon long-term memory. This is relevant in the context of fungal responses to dynamic and changing environments.
Our experiments and modeling have provided an example of identifying integral control that arises from time-scale separation in different processes, which is an important functional module in various contexts.
Adaptation to dynamically changing environments is required for any life form to survive. Unicellular microorganisms often adapt to environmental changes by modulating their metabolism and reprogramming their gene expression patterns. Generally, microbes are constantly monitoring their environment via specific receptors, and then adapting their cellular physiology accordingly via changes in gene expression and metabolism that are driven by specific signal transduction networks.
Osmosensing in Saccharomyces cerevisiae
The yeast Saccharomyces cerevisiae is able to adapt to acute changes in extracellular osmolarity. The hyperosmotic signal is amplified and transmitted via the High Osmolarity Glycerol (HOG) signaling network, and this leads to the intracellular accumulation of the osmolyte, glycerol. Increasing intracellular glycerol concentrations decrease the intracellular water potential, restore water influx and consequently restore cell volume and turgor pressure.
Osmosensing in Candida albicans
In contrast to the benign model yeast S. cerevisiae Candida albicans is a major fungal pathogen of humans that causes frequent mucosal infections in otherwise healthy individuals and potentially lethal infections in immunocompromised intensive care patients . HOG signaling is essential for the normal virulence of C. albicans, although osmosensing in C. albicans is not well characterised (Figure 1). Similar to S. cerevisiae, in C. albicans osmotic stress signals are relayed to Hog1 via Pbs2 . However, in C. albicans, osmotic stresses activate Hog1 only through Ssk2, and the Sho1 branch does not contribute significantly to osmosensing . Therefore, Ssk2 appears to be the only MAPKKK that activates the HOG MAPK module under the stress conditions examined so far . It is not yet clear whether the stress signal is relayed via Sln1. A novel osmotic stress signaling pathway might work in parallel with Sln1 to relay osmotic signals to Pbs2 . In contrast, bioinformatic surveys of osmotic signaling components have suggested that key components of the Sho1 branch might be absent in C. albicans (Additional file 1: Figure S1; ). In summary, uncertainties remain regarding the upstream and downstream components of the HOG signaling network in C. albicans.
A systems view of HOG signaling network
Given the importance of HOG signaling components for the virulence of C. albicans, it is important to study how this pathogen responds to osmotic shocks. In this study, we have investigated how the HOG signaling network operates in C. albicans using a model that encompasses the HOG signaling network and downstream adaptive processes.
Some signal transduction pathways returns to basal activity levels once a cell has adapted to environmental change. For instance, chemotaxis in E. coli is achieved through a sensing network whose activity always returns to the basal levels once the cell moves towards the attractant or away from the repellent . Also, the phosphorylation status of the heat shock transcription factor returns to basal levels once C. albicans cells have adapted to elevated ambient temperatures . This type of “perfect adaptation” may be significant in minimizing the potential impact of signaling crosstalk. In particular, the HOG pathway shares common components with the pheromone sensing pathway and the filamentous growth pathway. Perpetual HOG signaling activity might inappropriately activate other pathways.
The integral controller employs the time integral of the error signal K ∫ 0 t e(τ)dτ as the feedback mechanism.
Interestingly, integral control mechanisms have arisen frequently in evolution, presumably because they provide robust means by which systems can adapt effectively to constant environmental stimuli despite stochasticity. Indeed, systems biology studies have revealed a variety of biological processes that employ integral control mechanisms including bacterial chemotaxis , calcium homeostasis  and energy metabolism .
Elegant single cell analyses of nuclear Hog1 enrichment in S. cerevisiae following hyperosmotic shock have revealed that nuclear Hog1 levels return to basal levels once the cell has adapted to hyperosmotic conditions, and therefore that the HOG system displays perfect adaptation . These authors showed that upon exposure to osmotic stress, the time integral of Hog1 activity is linearly related to the increase in extracellular osmotic pressure within a certain range . This verifies that Hog1 exercises integral control over glycerol production in response to a persistent hyperosmolarity. This key discovery offers a starting point for developing simple mathematical models of osmotic stress responses in S. cerevisiae and C. albicans.
Using modeling and experimentation, we demonstrated that such integral control arise through a separation of the time scales between rapid signaling events and the relatively slower downstream adaptive processes. In addition, a hyperosmotic shock triggers long-term changes in the physiology of yeast cells that prepares them for subsequent hyperosmotic shocks. It appears as if a yeast cell can “remember” a previous hyperosmotic shock over a period that is longer than the time scale of signal transduction, and therefore adapts quickly to subsequent hyperosmotic shocks. We term this phenomenon a long-term stress memory. Our model suggests the HOG system also has a short-term memory that is engendered in the activation and inactivation of the intermediate kinases during the signal relay. The short-term memory might be affected by the regulation of glycerol channels. Our work highlights the need to study the properties of a signal transduction network in an appropriate biological context.
Cellular memory and model formulation
An alternative mechanism by which integral control can arise is via differential time scales in biological processes (see Figure 2C for details). Phosphorylated Hog1 activates several parallel mechanisms to increase the intracellular levels of glycerol in S. cerevisiae. Active Hog1 up-regulates the expression of genes encoding glycerol biosynthetic enzymes (GPD1 GPD2 and GPP2). Hog1 also increases the activity of a key enzyme (Pfk26) that diverts metabolic flux towards glycerol production . Furthermore, Hog1 stimulates the production of Stl1, a component of the glycerol/H+ symport system that assimilates glycerol from the growth medium [11, 30]. In principle, if the timescale for the inactivation any of these mechanisms is much longer than the timescale for achieving adaptation to higher osmolarity, then we can safely ignore the inactivation of the adaptive processes during the adaptation process. Since the glycerol channel activities change rapidly in response to external osmotic pressure , the adaptive processes we consider here are only related to glycerol production. Hence, if inactivation of glycerol production is a much slower process than adaptation itself, then the change in glycerol production rate (i.e. the output of the HOG signaling network) is determined by the time integral of phosphorylated Hog1. Indeed, in S. cerevisiae the rate of increase in glycerol biosynthetic enzyme levels is proportional to the time integral of Hog1 phosphorylation . This change in physiological status would be likely to generate a molecular memory whose duration is determined by the rate at which the glycerol production rate reduces to basal levels following a return to hypo-osmotic conditions. Therefore, a significant memory was predicted because the degradation of stable enzymes is relatively slow. A descriptive model detailing Hog1 integral control was experimentally validated . However, the idea that this integral control arises in S. cerevisiae through a separation of time scales between signaling and glycerol production rate changes was not fully recognized. The connection between Hog1 adaptation and integral control was not reported till later .
Whilst upstream components of HOG signaling have diverged significantly in the evolutionarily divergent yeast, C. albicans[19, 20], Hog1 integral control of osmotic stress adaptation might be conserved. Therefore our first aim was to build a simple mathematical model that highlights Hog1 integral control for both S. cerevisiae and C. albicans, and then use it to test alternative mechanisms by which a long-term memory can arise.
In our model we also assume that glycerol productivity does not decrease under hyperosmotic conditions (v7 = 0). Thus, the increase in the activity of G lycerol P roduction M achinery (GPM) will be proportional to the time integral of phosphorylation levels of Hog1. We also note that in the long term in S. cerevisiae, active Hog1 is reported to phosphorylate the aquaglyceroporin to target it for subsequent degradation [11, 32]. However, this regulation has yet to be quantified in S. cerevisiae, and it is not known whether the aquaglyceroporin in C. albicans is regulated in the same way. Hence, the degradation of aquaglyceroporin, despite its potential interest, is not included in our model.
The simple model explains published data
where F nuc is nuclear fluorescence signal, F cell is whole cell fluorescence.
This signal reduces to a small value following cellular adaptation. This suggests that the activity of aquaglyceroporin will return to its resting value b AGP after adaptation. Taken together, the rate of glycerol efflux is determined by and is proportional to the intracellular glycerol concentration G i alone. Because glycerol production P G balances glycerol efflux v 6 , the former is also proportional to intracellular glycerol concentration G i .
The increase in intracellular osmolarity caused by the rise in the intracellular glycerol concentration G i eventually matches the new extracellular osmolarity ∏ NaCl to reduce the signal upon adaptation. Thus, increase in the extracellular osmolarity that triggered the response ∏ NaCl is linearly proportional to the increase in glycerol production rate P G , which is determined by ∫ 0 t Hog 1P(t)dt due to the integral control. Hence, the time integral of the nuclear enrichment of Hog1, which is highly correlated with the phosphorylation of Hog1, has a linear relationship with ∏ NaCl .
Hyperosmotic responses in C. albicans
In response to hyperosmotic conditions, S. cerevisiae closes the aquaglyceroporin channel, Fps1, to prevent glycerol leakage and help to maintain a high intracellular glycerol concentration [27, 31]. Similarly, C. albicans appears to have a high restriction over glycerol efflux by an aquaglyceroporin. Our experimental data show that, in the presence of 1 M NaCl, the increase in total glycerol (i.e. the sum of intra- and extra-cellular glycerol) is mainly due to the induction of intracellular glycerol, and extracellular glycerol (i.e. the difference between total and intracellular glycerol) does not change significantly (Figure 6D). This is consistent with the idea that the aquaglyceroporin closes quickly in response to hyperosmotic stress thereby inhibiting glycerol leakage. This would maximize the efficiency with which cells could accumulate glycerol. However, we note that the experimental data suggest that intracellular glycerol does not plateau upon Hog1 adaptation (Figure 6B6D).
Our measurements of intracellular and total glycerol concentrations revealed relatively high initial intracellular glycerol concentrations in unstressed cells, which were induced about four-fold following a hyperosmotic shock. Our experimental procedures, which examined exponential cells grown in rich media at 30°C, were designed to minimize environmental stress. We would point out that the glycerol assay does not have the sensitivity to measure low concentrations of glycerol accurately . Nevertheless, our data indicate that under our experimental conditions unstressed C. albicans cells contain intracellular glycerol (Figure 6D). In our model we assume that the intracellular glycerol concentration is in equilibrium with the extracellular glycerol concentration pre-stress.
Our experimental data revealed two significant differences between the responses of C. albicans and S. cerevisiae. Firstly, in S. cerevisiae the proportion of phosphorylated Hog1 reaches its maximum value when concentrations of NaCl are greater than or equal to about 0.2 M (Figure 6A), whereas in C. albicans the percentage of phosphorylated Hog1 continues to increase, even in the presence of 0.5 M NaCl (data not shown). The NaCl concentration at which Hog1 phosphorylation amplitude becomes saturated was used to constrain parameters of the HOG signaling pathway (for details see supplemental information, section 2.1). Secondly, in S. cerevisiae GPD1 expression was strongly induced at the protein level and this regulation was mediated by Hog1 . However, in C. albicans GPD1 is expressed at high levels even in the absence of hyperosmotic stress, and Gpd1 is not highly induced under hyperosmotic conditions . However, expression of GPD2 was found to be induced at the mRNA level . This indicates differential regulation of the Gpd isoenzymes by Hog1 in C. albicans.
As mentioned previously, the existence of a long term memory is the key to distinguishing the mode of control that Hog1 exerts upon glycerol production (i.e. proportional, derivative or integral) and to explain its underpinning mechanism (i.e. saturated enzymatic kinetics or the separation of time scales). The existence of a long-term cellular memory would imply a mechanism of Hog1 integral control that arises through the separation of time scales between Hog1 signaling events and glycerol production rate changes.
We then tested whether this molecular memory was due to the retention of a relatively high glycerol synthetic capacity after adaptation to the first hyperosmotic stress. This was achieved by monitoring the intracellular glycerol levels generated by C. albicans cells after repeated exposure to hyperosmotic conditions (Figure 7C). The data showed that intracellular glycerol levels plateau within the first 45 min, which was coincident with the disappearance of Hog1 phosphorylation during the first hyperosmotic shock. Intracellular glycerol levels were then restored within 20 min of the second shock, also coincident with the disappearance of Hog1 phosphorylation. This provided further evidence that C. albicans cells retain high glycerol productivity for some time after the first adaptation. It also strongly supports the idea that the integral control in the HOG system is endowed by the separation of time scales. We note that the total glycerol was induced about 2-fold while it is simulated to increase 4-fold. However this minor discrepancy does not change the main conclusions drawn above.
Duration of long-term memory
Long-term memory impacts on short-term memory
The HOG signaling networks in S. cerevisiae and C. albicans also possesses a “short-term memory” of hyperosmolarity. This results from the time it takes for the signaling components to be inactivated via dephosphorylation/phosphorylation. In response to a fluctuating osmotic signal, if the duration of the “off” phase of the signal is longer than the time required to inactivate the signaling components on HOG pathway, the phosphorylation levels of Hog1 would be predicted to follow the fluctuating signal. However, if the duration of the “off” phase is shorter than the inactivation time of the HOG pathway, the HOG network would remain active into the next cycle. The short-term memory of HOG pathway (i.e. the inactivation time for HOG pathway) typically lasts several minutes [37, 38].
To better understand this phenomenon, we developed a theoretical framework for the basic HOG signaling module, which is represented by a protein that is phosphorylated and dephosphorylated by its kinase and phosphatase, respectively. Employing Michaelis-Menten kinetics, we show that this module is essentially a low-pass filter (supplemental section 2.2). In other words, given a constant phosphatase level, the phosphorylation level of the protein faithfully follows a low-frequency activation signal while integrating a high-frequency signal, as reported previously [28, 37, 38]. We also derived analytically how the frequency of an activation signal would impact upon the phosphorylation level of an enzyme (see supplemental section 2.2 for details of the solution). In this way, our analytic solution of the frequency response has furthered our understanding of the behaviour of this simple module. In addition, we showed that in response to an alternating signal, the amplitude of the output of a network decreases with the number of cascades it contains (supplemental section 2.2). This ability to attenuate high frequency signals would allow cells to filter out random fluctuations (white noise) in environmental osmolarities, whilst maintaining its ability to sense significant environmental changes over longer time scales.
Meanwhile, the HOG signaling network integrates high frequency signals. It is expected that a cell will respond to a high frequency sinusoidal signal in a similar manner as to a constant hyperosmotic stress with the time-average amplitude of the sine signal. In addition, due to the integral control, Hog1 phosphorylation levels might be expected to return to basal levels once the cell has adapted. However, our simulations predict that Hog1 phosphorylation does not return to the initial condition. Instead the model predicts that, due to the quick release of glycerol in the “off” phase of the signal, the cell does not accumulate intracellular glycerol to the same level as it does under a constant signal (Figure 9B), and the glycerol production rate does not reach the same level as that under a constant stress (Figure 9E). This would be similar to the molecular response of mutants that are unable to close the aquaglyceroporin channel, which must constantly activate Hog1 to maintain a high intracellular glycerol concentration. Based on this consideration, we designed an asymmetric square wave signal with the same frequency of 1 rad/min, but where each “on” phase is longer than each “off” phase. The simulations predicted that under these conditions, Hog1 phosphorylation levels will recover when the “off” phase is sufficiently short (see supplement section 2.6 for details). These results suggest that the frequency of the cellular response is dependent upon both the short-term memory in the signaling pathway (i.e. the low-pass filter) and the property of glycerol channels. The results highlight the need to place a model of the signal transduction network into an appropriate context when studying its properties.
Recently, Gennemark and coworkers reported a simple mathematical model of hyperosmotic responses in S. cerevisiae. Their model concisely described the biophysical changes of a hyperosmotically stressed yeast cell, and related those changes to glycerol production. This model is able to reproduce some published data on the changes in intracellular and extracellular glycerol concentrations in response to different physiological perturbations. Instead of integral control, their model assumes a time-delayed proportional control of glycerol production, i.e. glycerol productivity at a particular time is assumed to be proportional to the change in turgor pressure (relative to a constant reference turgor pressure). A consequence of this assumption is that there should be no memory effects in the stress response dynamics. This motivated us to investigate whether a model including integral control would be more appropriate.
Also inspired by the recent finding of Hog1 integral control in S. cerevisiae, we investigated osmotic stress responses in C. albicans by constructing a low granularity mathematical model assuming Hog1 integral control. This model integrates the HOG signal transduction pathway with downstream adaptive processes, and encompasses both hyper- and hypo-osmotic conditions. In this model, a signaling or biochemical pathway was concisely represented by an overall reaction (Figure 5). Hence, the property of the pathway was defined by the associated kinetic parameters. The model was calibrated using in-house generated time-course data of Hog1 phosphorylation plus intracellular and total glycerol levels following exposure to two doses of NaCl. Subsequently, Hog1 integral control was experimentally verified in C. albicans using experiments in which cells were subjected to repeated NaCl insults. The data showed that pre-stressed cells adapted to a second hyperosmotic shock more quickly, as if they remembered the hyperosmotic history. This phenomenon is known as the “long-term stress memory”, because its duration was longer than the time required for signal transduction to take place. Using this model, we reasoned that the cellular memory observed under hyperosmotic conditions is a hallmark of the separation of time scales between HOG signaling and reduction of glycerol production rates. In this case this mechanism enables Hog1 integral control.
The separation of time scales in different biological processes is a universal phenomenon. Often, signaling networks that perceive environmental conditions such as nutrient availability, chemical insults, and mating pheromones are required to function at much faster rates (i.e. seconds) than the ensuing regulatory networks which encompass transcriptional and translational controls (i.e. tens of minutes). It is also common for signals to be converted, via modulation of transcriptional regulators, into the expression of a stable product (i.e. a stable mRNA or protein). Such a separation of time scales could result in an integral control mechanism that would generate an experimentally testable long-term memory. Therefore our experimental approach could be further applied to other signaling networks to examine such a property. In addition, the identification of such a phenomenon allows mathematical simplifications of the regulatory network under study, an important step towards a thorough understanding of signaling networks which may be obscured in exhaustive descriptive models.
The teleological question as to why the yeast cell employs an integral control mechanism remains to be answered. Among the various types of feedback mechanism, proportional control and derivative control mechanisms would allow a system to adapt to a new status that is different from the original condition. However, integral control ensures perfect adaptation of the system, whereby signaling returns to basal levels once adaptation is achieved. This would reduce the potential for cross-talk between signaling pathways. The ability of the cell to restore turgor pressure in the face of an osmotic challenge is crucial for the restoration of growth. Should Hog1 signaling regulate glycerol production via proportional control, a change in osmotic pressure would be reflected by a proportional change in Hog1 activity. In other words, the ability to sense diverse amplitudes of osmotic signal would require high concentrations of HOG signaling proteins, the production of which would demand considerable cellular resources. In contrast, the integral control mechanism converts the error signal into the temporal integral of HOG activity. This avoids the need to express HOG pathway components at high levels. In fact in S. cerevisiae, Hog1 phosphorylation levels become saturated in response to relatively low hyperosmotic signals. In addition, derivative control offers gradual changes in a system, which might be too slow to generate the necessary adaptation to external osmotic challenges. Indeed, derivative control is often found to provide a buffering capacity in cellular systems. For instance, in energy metabolism, phosphorylated molecules, such as phosphocreatine, act as a buffer during high energy demand periods. Phosphocreatine implements a biological derivative control over ATP concentration .
Our modeling further suggests that the frequency response of the HOG signaling pathway can be affected by the regulation over the glycerol channels. This prediction is significant because this frequency response is generally thought to be a property of the short-term memory mediated by the regulation of intermediate protein kinases alone. In addition, our simulations show that the long-term memory may also affect the way a cell responds to a high frequency signal. Our simulations suggest that a cell that has previously adapted to hyperosmolarity retains relatively high intracellular glycerol levels for a period. Consequently, the high intracellular osmolarity counterbalances the hyperosmotic signal that the cell perceives. These results highlight the importance of a systems approach to study the signaling networks.
Our measurements of Hog1 phosphorylation dynamics in C. albicans suggest that, in this pathogen, Hog1 achieves perfect adaptation to a constant hyperosmotic shock through an integral control mechanism that is dependent on separation in time scales between signalling and downstream adaptation. We also predict that Hog1 signalling exhibits a short-term memory that can be revealed by tracking responses to an alternating signal. Furthermore, we predict that this short-term memory depends on the long-term memory. For example, Hog1 phosphorylation levels are predicted to adapt relatively quickly to high frequency signals in cells that have pre-adapted to a constant hyperosmotic stress.
Strains and culture conditions
C. albicans strain CA1655 (ura3::λimm434/ura3::λimm434 RPS1/rps1::pACT1-FLAG-GFP) was used for the Hog1-phosphorylation assays. This strain is CAI-4  containing pACT1-FLAG-GFP  integrated at the RPS1 locus. C. albicans CA1655 expresses a FLAG-tagged version of GFP from the ACT1 promoter and functions as an internal loading control for Western blots. C. albicans NGY152 , which is CAI4 transformed with CIp10 , was used for the glycerol assays. Strains were routinely maintained on YPD-T medium (2 % (w/v) glucose, 2 % (w/v) mycological peptone (Oxoid), 1 % (w/v) yeast extract (Oxoid), 100 mM Tris–HCl, pH 7.4) .
Cell growth and stress experiments
Cells were grown at 30°C with shaking (200 rpm). To ensure consistency in the starting physiology of the cultures used in our experiments, cells were pre-cultured twice overnight in YPD-T, diluted to an OD600 of 0.2 in fresh medium and then grown to mid-log phase (OD600 0.8) before any stress was applied . A single hyperosmotic stress was applied by splitting one large starting culture and diluting the cells back to an OD600 of 0.2 in fresh YPD-T with NaCl added to achieve concentrations of 0, 0.3, 0.5, 0.75 or 1 M. Cells were harvested by centrifugation at selected time points and the pellets flash frozen in liquid nitrogen.
For the experiments involving the imposition of sequential stresses, cells were initially cultured and stressed as described above. At the end of the first stress period, cells were harvested, resuspended in fresh YPD-T with no NaCl and incubated at 30°C for 10 min. NaCl was then added to the cultures and samples were collected after the second stress period. Three biological replicates were prepared for each condition.
Protein extraction and western blotting
To make protein extracts, flash-frozen cell pellets were resuspended in freshly prepared protein lysis buffer (50 mM Tris–HCl pH 7.5, 150 mM NaCl, 0.5 % NP40) containing inhibitors (2 mg/ml Leupeptin, 2 mg/ml Pepstatin, 1 mM PMSF, 2 mM Na3VO4 and 50 mM NaF) and sheared with acid-washed glass beads (425–600 μm, Sigma-Alrich) using a Fast-Prep FP120 machine (Thermo Scientific) (6 × 15 s bursts, Speed 6.5 with 1 min on ice between bursts). Protein extracts were clarified by centrifugation and protein concentration determined using a Bradford assay .
Proteins were separated by SDS-PAGE on NuPAGE®Novex Bis-Tris 4-12 % pre-cast polyacrylamide gels (Invitrogen, Paisley, UK). 15 μg protein was loaded in each lane. Proteins were transferred to PVDF membranes, rinsed in phosphate-buffered saline (PBS) and then blocked in PBS-T containing 10 % BSA (PBS, 0.1 % Tween20, 10 % (w/v) BSA) for 30 min at room temperature. Blots were incubated overnight at 4°C in PBS-T containing 5%BSA containing a 1:2000 dilution of anti-phospho-p38 MAPK (Thr180/Tyr182) antibody (Cell Signaling Technology) and a 1:20000 dilution of anti-FLAG antibody (Sigma-Aldrich, UK). Blots were washed in PBS-T and then incubated for 1 h at room temperature with Anti-rabbit IgG, HRP-linked antibody (Cell Signaling Technology). The blots were washed in PBS-T and the signal was detected using LumiGLO™ (Cell Signaling Technology) or SuperSignal®WestFemto (Thermo Scientific) according to manufacturer’s instructions. Phosphorylated Hog1 and GFP-FLAG levels were visualized and quantified using the FluorChem®FC2 (Alpha Innotech) or the FUSION SL™ (Peqlab) systems and the levels of phosphorylated Hog1 were expressed relative to the internal standard (GFP-FLAG).
Extracellular and total glycerol concentrations were determined in stressed and unstressed cells at various time points. Cultures were grown and the stress applied as described above. The Free Glycerol Determination Kit (Sigma-Aldrich) was used to determine glycerol concentrations according to the manufacturer’s instructions. To measure extracellular glycerol concentrations, 1.5 ml samples of culture were collected and the cells harvested by centrifugation. Supernatants were transferred to a fresh tube and heated at 100°C for 10 min before their glycerol content was determined. Total glycerol concentrations were measured using 1.5 ml of culture after heating at 100°C for 10 min without prior removal of cells. Cell debris was spun down and the supernatant assayed for glycerol content. Intracellular glycerol concentrations were calculated by deducting the extracellular glycerol concentration from the total glycerol concentration.
Availability of supporting data
Supporting data on the following topics are collated in an additional Supplemental data file with the manuscript:
Bioinformatic analysis of the conservation of proteins in HOG signaling network in Candida albicans
· Steady-state response of HOG signaling network
· Frequency response of HOG signaling network
· Model assumptions
· The model
· Model parameterisation
· Frequency response of the system
This paper is dedicated to Jaroslav Stark, our source of wisdom and inspiration, who tragically died in 2010. TY would like to thank Wei Pang for the valuable discussions on model development and Christopher Brackley for the help on the analysis of frequency response. All authors would like to thank Laura Selway for undertaking some Gpd1 time course experimental work, and Despoina Kaloriti, Anna Tillman, Melanie Puttnam, Elahe Radmaneshfar, George Coghill, Janet Quinn and Mauricio Barahona for general discussions. All authors gratefully acknowledge the funding support by the CRISP project (Combinatorial Responses In Stress Pathways) funded by the BBSRC (BB/F00513X/1) under the Systems Approaches to Biological Research (SABR) Initiative. NARG and AJPB were also funded by the European Marie Curie Initial Training Networks, ARIADNE and FINSysB, respectively (PITN-GA-2008-214004). AJPB was also supported by an ERC Advanced Grant (STRIFE: ERC-2009-AdG_20090506).
- Maeda T, Wurgler-Murphy SM, Saito H: A two-component system that regulates an osmosensing MAP kinase cascade in yeast. Nature. 1994, 369: 242-245. 10.1038/369242a0.PubMedView ArticleGoogle Scholar
- Posas F, Wurgler-Murphy SM, Maeda T, Witten EA, Thai TC, Saito H: YeastHOG1MAP kinase cascade is regulated by a multistep phosphorelay mechanism in the SLN1-YPD1-SSK1 “two-component” osmosensor. Cell. 1996, 86: 865-875. 10.1016/S0092-8674(00)80162-2.PubMedView ArticleGoogle Scholar
- Maeda T, Takekawa M, Saito H: Activation of yeast PBS2 MAPKK by MAPKKKs or by binding of an SH3-containing osmosensor. Science. 1995, 269: 554-558. 10.1126/science.7624781.PubMedView ArticleGoogle Scholar
- Posas F, Saito H: Osmotic activation of the HOG MAP kinase pathway via Ste11p MAPKKK: scaffold role of Pbs2p MAPKK. Science. 1997, 276: 1702-1705. 10.1126/science.276.5319.1702.PubMedView ArticleGoogle Scholar
- Hohmann S: Osmotic stress signaling and osmoadaptation in Yeasts. Microbiol Mol Biol Rev. 2002, 66: 300-372. 10.1128/MMBR.66.2.300-372.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Hao N, Behar M, Parnell SC, Torres MP, Borchers CH, Elston TC, Dohlman HG: A systems-biology analysis of feedback inhibition in the Sho1 osmotic-stress-response pathway. Curr Biol. 2007, 17: 659-667. 10.1016/j.cub.2007.02.044.PubMedView ArticleGoogle Scholar
- Hao N, Zeng Y, Elston TC, Dohlman HG: Control of MAPK specificity by feedback phosphorylation of shared adaptor protein Ste50. J Biol Chem. 2008, 283: 33798-33802. 10.1074/jbc.C800179200.PubMedPubMed CentralView ArticleGoogle Scholar
- Sturm OE, Orton R, Grindlay J, Birtwistle M, Vyshemirsky V, Gilbert D, Calder M, Pitt A, Kholodenko B, Kolch W: The mammalian MAPK/ERK pathway exhibits properties of a negative feedback amplifier. Sci Signal. 2010, 3 (153): ra90-10.1126/scisignal.2001212.PubMedView ArticleGoogle Scholar
- Macia J, Regot S, Peeters T, Conde N, Solé R, Posas F: Dynamic signaling in the Hog1 MAPK pathway relies on high basal signal transduction. Sci Signal. 2009, 2: ra13-10.1126/scisignal.2000056.PubMedView ArticleGoogle Scholar
- Sato N, Kawahara H, Toh-e A, Maeda T: Phosphorelay-regulated degradation of the yeast Ssk1p response regulator by the ubiquitin-proteasome system. Mol Cell Biol. 2003, 23 (18): 6662-6671. 10.1128/MCB.23.18.6662-6671.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Hohmann S: Control of high osmolarity signaling in the yeast Saccharomyces cerevisiae. FEBS Lett. 2009, 583: 4025-4029. 10.1016/j.febslet.2009.10.069.PubMedView ArticleGoogle Scholar
- Tatebayashi K, Tanaka K, Yang HY, Yamamoto K, Matsushita Y, Tomida T, Imai M, Saito H: Transmembrane mucins Hkr1 and Msb2 are putative osmosensors in the SHO1 branch of yeast HOG pathway. EMBO J. 2007, 26: 3521-3533. 10.1038/sj.emboj.7601796.PubMedPubMed CentralView ArticleGoogle Scholar
- Takahashi S, Pryciak PM: Identification of novel membrane-binding domains in multiple yeast Cdc42 effectors. Mol Biol Cell. 2007, 18: 4945-4956. 10.1091/mbc.E07-07-0676.PubMedPubMed CentralView ArticleGoogle Scholar
- Tatebayashi K, Yamamoto K, Tanaka K, Tomida T, Maruoka T, Kasukawa E, Saito H: Adaptor functions of Cdc42, Ste50, and Sho1 in the yeast osmoregulatory HOG MAPK pathway. EMBO J. 2006, 25: 3033-3044. 10.1038/sj.emboj.7601192.PubMedPubMed CentralView ArticleGoogle Scholar
- Brown AJP, Haynes K, Gow NAR, Quinn J: Stress responses inCandida. Candida and Candidiasis. Edited by: Clancy CJ, Calderone RA. 2011, ASM Press, Clancy CJ, Calderone RA, in pressGoogle Scholar
- Calderone RA: Candida and Candidiasis. 2002, ASM Press, Washington, D.CGoogle Scholar
- Alonso-Monge R, Navarro-Garcia F, Molero G, Diez-Orejas R, Gustin M, Pla J, Sanchez M, Nombela C: Role of the mitogen-activated protein kinase Hog1p in morphogenesis and virulence ofCandida albicans. J Bacteriol. 1999, 181: 3058-3068.PubMedPubMed CentralGoogle Scholar
- Arana DM, Nombela C, Alonso-Monge R, Pla J: The Pbs2 MAP kinase kinase is essential for the oxidative-stress response in the fungal pathogenCandida albicans. Microbiol. 2005, 151: 1033-1049. 10.1099/mic.0.27723-0.View ArticleGoogle Scholar
- Cheetham J, Smith DA, da Silva DA, Doris KS, Patterson MJ, Bruce CR, Quinn J: A single MAPKKK regulates the Hog1 MAPK pathway in the pathogenic fungusCandida albicans. Mol Bio Cell. 2007, 18: 4603-4614. 10.1091/mbc.E07-06-0581.View ArticleGoogle Scholar
- Nikolaou E, Agrafioti I, Stumpf M, Quinn J, Stansfield I, Brown AJP: Phylogenetic diversity of stress signalling pathways in fungi. BMC Evol Biol. 2009, 9: 44-10.1186/1471-2148-9-44.PubMedPubMed CentralView ArticleGoogle Scholar
- Alon U, Surette MG, Barkai N, Leibler S: Robustness in bacterial chemotaxis. Nature. 1999, 397: 168-171. 10.1038/16483.PubMedView ArticleGoogle Scholar
- Leach MD, Tyc KM, Brown AJP, Klipp E: Modelling the regulation of thermal adaptation inCandida albicans, a major fungal pathogen of humans. PLoS One. 2012, 7: e3246-Google Scholar
- Ma W, Trusina A, El-Samad H, Lim WA, Tang C: Defining network topologies that can achieve biochemical adaptation. Cell. 2009, 138 (4): 760-73. 10.1016/j.cell.2009.06.013.PubMedPubMed CentralView ArticleGoogle Scholar
- Yi TM, Huang Y, Simon MI, Doyle J: Robust perfect adaptation in bacterial chemotaxis through integral feedback control. Proc Natl Acad Sci U S A. 2000, 97: 4649-4653. 10.1073/pnas.97.9.4649.PubMedPubMed CentralView ArticleGoogle Scholar
- El-Samad H, Goff JP, Khammash M: Calcium homeostasis and parturient hypocalcemia: an integral feedback perspective. J Theor Biol. 2002, 214: 17-29. 10.1006/jtbi.2001.2422.PubMedView ArticleGoogle Scholar
- Cloutier M, Wellstead P: The control systems structures of energy metabolism. J R Soc Interface. 2010, 7: 651-665. 10.1098/rsif.2009.0371.PubMedPubMed CentralView ArticleGoogle Scholar
- Muzzey D, Gómez-Uribe CA, Mettetal JT, van Oudenaarden A: A systems-level analysis of perfect adaptation in yeast osmoregulation. Cell. 2009, 138: 160-171. 10.1016/j.cell.2009.04.047.PubMedPubMed CentralView ArticleGoogle Scholar
- Gomez-Uribe C, Verghese GC, Mirny LA: Operating regimes of signaling cycles: statics, dynamics, and noise filtering. PLoS Comp Biol. 2007, 3: 2487-2497.View ArticleGoogle Scholar
- Dihazi H, Kessler R, Eschrich K: High osmolarity glycerol (HOG) pathway-induced phosphorylation and activation of 6-phosphofructo-2-kinase are essential for glycerol accumulation and yeast cell proliferation under hyperosmotic stress. J Biol Chem. 2004, 279 (23): 23961-23968. 10.1074/jbc.M312974200.PubMedView ArticleGoogle Scholar
- Westfall PJ, Patterson JC, Chen RE, Thorner J: Stress resistance and signal fidelity independent of nuclear MAPK function. Proc Natl Acad Sci U S A. 2008, 105: 12212-12217. 10.1073/pnas.0805797105.PubMedPubMed CentralView ArticleGoogle Scholar
- Klipp E, Nordlander B, Krüger R, Gennemark P, Hohmann S: Integrative model of the response of yeast to osmotic shock. Nat Biotech. 2005, 23: 975-982. 10.1038/nbt1114.View ArticleGoogle Scholar
- Mollapour M, Piper PW: Hog1 mitogen-activated protein kinase phosphorylation targets the yeast Fps1 aquaglyceroporin for endocytosis, thereby rendering cells resistant to acetic acid. Mol Cell Biol. 2007, 27: 6446-6456. 10.1128/MCB.02205-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Gennemark P, Nordlander B, Hohmann S, Wedelin D: A simple mathematical model of adaptation to high osmolarity in yeast. In Silico Biol. 2006, 6: 193-214.PubMedGoogle Scholar
- Reed RH, Chudek JA, Foster R, Gadd GM: Osmotic significance of glycerol accumulation in exponentially growing yeasts. Appl Environ Microbiol. 1987, 53: 2119-2123.PubMedPubMed CentralGoogle Scholar
- Zi Z, Liebermeister W, Klipp E: A quantitative study of the Hog1 MAPK response to fluctuating osmotic stress in Saccharomyces cerevisiae. PLoS One. 2010, 5: e9522-10.1371/journal.pone.0009522.PubMedPubMed CentralView ArticleGoogle Scholar
- Enjalbert B, Smith DA, Cornell MJ, Alam I, Nicholls S, Brown AJP, Quinn J: Role of the Hog1 Stress-Activated Protein Kinase in the Global Transcriptional Response to Stress in the Fungal PathogenCandida albicans. Molec. Biol. Cell. 2006, 17: 1018-1032.PubMedPubMed CentralView ArticleGoogle Scholar
- Hersen P, McClean MN, Mahadevan L, Ramanathan S: Signal processing by the HOG MAP kinase pathway. Proc Natl Acad Sci U S A. 2008, 105: 7165-7170. 10.1073/pnas.0710770105.PubMedPubMed CentralView ArticleGoogle Scholar
- Mettetal JT, Muzzey D, Gómez-Uribe C, van Oudenaarden A: The frequency dependence of osmo-adaptation inSaccharomyces cerevisiae. Science. 2008, 319: 482-484. 10.1126/science.1151582.PubMedPubMed CentralView ArticleGoogle Scholar
- Fonzi WA, Irwin MY: Isogenic strain construction and gene mapping inCandida albicans. Genetics. 1993, 134: 717-728.PubMedPubMed CentralGoogle Scholar
- Nicholls S, Leach M, Priest C, Brown AJP: Role of the heat shock transcription factor, Hsf1, in a major fungal pathogen that is obligately associated with warm-blooded animals. Molec. Microbiol. 2009, 74: 844-861. 10.1111/j.1365-2958.2009.06883.x.View ArticleGoogle Scholar
- Walker LA, MacCallum DM, Bertram G, Gow NAR, Odds FC, Brown AJP: Genome-wide analysis ofCandida albicansgene expression patterns during infection of the mammalian kidney. Fungal Gen Biol. 2009, 46: 210-219. 10.1016/j.fgb.2008.10.012.View ArticleGoogle Scholar
- Murad AM, Lee PR, Broadbent ID, Barelle CJ, Brown AJ: CIp10, an efficient and convenient integrating vector forCandida albicans. Yeast. 2000, 16: 325-327. 10.1002/1097-0061(20000315)16:4<325::AID-YEA538>3.0.CO;2-#.PubMedView ArticleGoogle Scholar
- Kaloriti D, Tillmann A, Cook E, Jacobsen MD, You T, Lenardon MD, Ames L, Barahona M, Chandrasekaran K, Coghill G, Goodman D, Gow NAR, Grebogi C, Ho HL, Ingram P, McDonagh A, de Moura APS, Pang W, Puttnam M, Radmaneshfar E, Romano MC, Silk D, Stark J, Stumpf M, Thiel M, Thorne T, Usher J, Yin Z, Haynes K, Brown AJP: Combinatorial stresses kill pathogenicCandidaspecies. Medical Mycol. 2012, 10.3109/13693786.2012.672770.Google Scholar
- Bradford MM: A rapid and sensitive method for the quantification of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976, 72: 248-254. 10.1016/0003-2697(76)90527-3.PubMedView ArticleGoogle Scholar