Quantitative analysis of glycerol accumulation, glycolysis and growth under hyper osmotic stress.

We provide an integrated dynamic view on a eukaryotic osmolyte system, linking signaling with regulation of gene expression, metabolic control and growth. Adaptation to osmotic changes enables cells to adjust cellular activity and turgor pressure to an altered environment. The yeast Saccharomyces cerevisiae adapts to hyperosmotic stress by activating the HOG signaling cascade, which controls glycerol accumulation. The Hog1 kinase stimulates transcription of genes encoding enzymes required for glycerol production (Gpd1, Gpp2) and glycerol import (Stl1) and activates a regulatory enzyme in glycolysis (Pfk26/27). In addition, glycerol outflow is prevented by closure of the Fps1 glycerol facilitator. In order to better understand the contributions to glycerol accumulation of these different mechanisms and how redox and energy metabolism as well as biomass production are maintained under such conditions we collected an extensive dataset. Over a period of 180 min after hyperosmotic shock we monitored in wild type and different mutant cells the concentrations of key metabolites and proteins relevant for osmoadaptation. The dataset was used to parameterize an ODE model that reproduces the generated data very well. A detailed computational analysis using time-dependent response coefficients showed that Pfk26/27 contributes to rerouting glycolytic flux towards lower glycolysis. The transient growth arrest following hyperosmotic shock further adds to redirecting almost all glycolytic flux from biomass towards glycerol production. Osmoadaptation is robust to loss of individual adaptation pathways because of the existence and upregulation of alternative routes of glycerol accumulation. For instance, the Stl1 glycerol importer contributes to glycerol accumulation in a mutant with diminished glycerol production capacity. In addition, our observations suggest a role for trehalose accumulation in osmoadaptation and that Hog1 probably directly contributes to the regulation of the Fps1 glycerol facilitator. Taken together, we elucidated how different metabolic adaptation mechanisms cooperate and provide hypotheses for further experimental studies.

Network modeling of the transcriptional effects of copy number aberrations in glioblastoma.

DNA copy number aberrations (CNAs) are a hallmark of cancer genomes. However, little is known about how such changes affect global gene expression. We develop a modeling framework, EPoC (Endogenous Perturbation analysis of Cancer), to (1) detect disease-driving CNAs and their effect on target mRNA expression, and to (2) stratify cancer patients into long- and short-term survivors. Our method constructs causal network models of gene expression by combining genome-wide DNA- and RNA-level data. Prognostic scores are obtained from a singular value decomposition of the networks. By applying EPoC to glioblastoma data from The Cancer Genome Atlas consortium, we demonstrate that the resulting network models contain known disease-relevant hub genes, reveal interesting candidate hubs, and uncover predictors of patient survival. Targeted validations in four glioblastoma cell lines support selected predictions, and implicate the p53-interacting protein Necdin in suppressing glioblastoma cell growth. We conclude that large-scale network modeling of the effects of CNAs on gene expression may provide insights into the biology of human cancer. Free software in MATLAB and R is provided.

Robustness and fragility in the yeast high osmolarity glycerol (HOG) signal-transduction pathway.

Cellular signalling networks integrate environmental stimuli with the information on cellular status. These networks must be robust against stochastic fluctuations in stimuli as well as in the amounts of signalling components. Here, we challenge the yeast HOG signal-transduction pathway with systematic perturbations in components' expression levels under various external conditions in search for nodes of fragility. We observe a substantially higher frequency of fragile nodes in this signal-transduction pathway than that has been observed for other cellular processes. These fragilities disperse without any clear pattern over biochemical functions or location in pathway topology and they are largely independent of pathway activation by external stimuli. However, the strongest toxicities are caused by pathway hyperactivation. In silico analysis highlights the impact of model structure on in silico robustness, and suggests complex formation and scaffolding as important contributors to the observed fragility patterns. Thus, in vivo robustness data can be used to discriminate and improve mathematical models.

Exploring the impact of osmoadaptation on glycolysis using time-varying response-coefficients.

We present a model of osmoadaptation in S. cerevisiae based on existing experimental and theoretical work. In order to investigate the impact of osmoadaptation on glycolysis, this model focuses on the interactions between glycolysis and osmoadaptation, namely the production of glycerol and its influence on flux towards pyruvate. Evaluation of this model shows that, depending on initial relations between glycerol and pyruvate production, the increased glycerol production can have a substantial negative effect on the pyruvate production rate. Existing experimental data and a detailed analysis of the model lead to the suggestion of an interaction between activated Hog1 and activators of glycolysis such as Pfk26.

Integrative model of the response of yeast to osmotic shock.

Integration of experimental studies with mathematical modeling allows insight into systems properties, prediction of perturbation effects and generation of hypotheses for further research. We present a comprehensive mathematical description of the cellular response of yeast to hyperosmotic shock. The model integrates a biochemical reaction network comprising receptor stimulation, mitogen-activated protein kinase cascade dynamics, activation of gene expression and adaptation of cellular metabolism with a thermodynamic description of volume regulation and osmotic pressure. Simulations agree well with experimental results obtained under different stress conditions or with specific mutants. The model is predictive since it suggests previously unrecognized features of the system with respect to osmolyte accumulation and feedback control, as confirmed with experiments. The mathematical description presented is a valuable tool for future studies on osmoregulation in yeast and-with appropriate modifications-other organisms. It also serves as a starting point for a comprehensive description of cellular signaling.

A microfluidic system in combination with optical tweezers for analyzing rapid and reversible cytological alterations in single cells upon environmental changes.

We report on the development of an experimental platform where epi-fluorescence microscopy and optical tweezers are combined with a microfluidic system to enable the analysis of rapid cytological responses in single cells. The microfluidic system allows two different media to be merged in a Y-shaped channel. Microscale channel dimensions ensure purely laminar flow and, as a result, an environmental gradient can be created between the two media. Optical tweezers are used to move a single trapped cell repeatedly between the different environments. The cell is monitored continuously by fluorescence microscopy during the experiment. In a first experiment on yeast (Saccharomyces cerevisiae) we observed changes in cell volume as the cell was moved between environments with different osmolarity. This demonstrated that the platform allowed analysis of cytological alterations on a time scale shorter than 0.2 s. In a second experiment we observed the spatial migration of the Yap1p transcription factor fused to GFP as a cell was moved from an environment of low to high oxidative capacity. The system is universal allowing the response to numerous environmental changes to be studied on the sub second time scale in a variety of model cells. We intend to use the platform to study how the age of cells, their progression through the cell cycle, or their genetic landscape, alter their capacity (kinetics and amplitude) to respond to environmental changes.

Yeast osmoregulation.

Osmoregulation is the active control of the cellular water balance and encompasses homeostatic mechanisms crucial for life. The osmoregulatory system in the yeast Saccharomyces cerevisiae is particularly well understood. Key to yeast osmoregulation is the production and accumulation of the compatible solute glycerol, which is partly controlled by the high osmolarity glycerol (HOG) signaling system. Genetic analyses combined with studies on protein-protein interactions have revealed the wiring scheme of the HOG signaling network, a branched mitogen-activated protein (MAP) kinase (MAPK) pathway that eventually converges on the MAPK Hog1. Hog1 is activated following cell shrinking and controls posttranscriptional processes in the cytosol as well as gene expression in the nucleus. HOG pathway activity can easily and rapidly be controlled experimentally by extracellular stimuli, and signaling and adaptation can be separated by a system of forced adaptation. This makes yeast osmoregulation suitable for studies on system properties of signaling and cellular adaptation via mathematical modeling. Computational simulations and parallel quantitative time course experimentation on different levels of the regulatory system have provided a stepping stone toward a holistic understanding, revealing how the HOG pathway can combine rigorous feedback control with maintenance of signaling competence. The abundant tools make yeast a suitable model for an integrated analysis of cellular osmoregulation. Maintenance of the cellular water balance is fundamental for life. All cells, even those in multicellular organisms with an organism-wide osmoregulation, have the ability to actively control their water balance. Osmoregulation encompasses homeostatic processes that maintain an appropriate intracellular environment for biochemical processes as well as turgor of cells and organism. In the laboratory, the osmoregulatory system is studied most conveniently as a response to osmotic shock, causing rapid and dramatic changes in the extracellular water activity. Those rapid changes mediate either water efflux (hyperosmotic shock), and hence cell shrinkage, or influx (hypoosmotic shock), causing cell swelling. The yeast S. cerevisiae, as a free-living organism experiencing both slow and rapid changes in extracellular water activity, has proven a suitable and genetically tractable experimental system in studying the underlying signaling pathways and regulatory processes governing osmoregulation. Although far from complete, the present picture of yeast osmoregulation is both extensive and detailed (de Nadal et al., 2002; Hohmann, 2002; Klipp et al., 2005). Simulations using mathematical models combined with time course measurements of different molecular processes in signaling and adaptation have allowed elucidation of the first system properties on the yeast osmoregulatory network.

Initiation of the transcriptional response to hyperosmotic shock correlates with the potential for volume recovery.

The control of activity and localization of transcription factors is critical for appropriate transcriptional responses. In eukaryotes, signal transduction components such as mitogen-activated protein kinase (MAPK) shuttle into the nucleus to activate transcription. It is not known in detail how different amounts of nuclear MAPK over time affect the transcriptional response. In the present study, we aimed to address this issue by studying the high osmolarity glycerol (HOG) system in Saccharomyces cerevisiae. We employed a conditional osmotic system, which changes the period of the MAPK Hog1 signal independent of the initial stress level. We determined the dynamics of the Hog1 nuclear localization and cell volume by single-cell analysis in well-controlled microfluidics systems and compared the responses with the global transcriptional output of cell populations. We discovered that the onset of the initial transcriptional response correlates with the potential of cells for rapid adaptation; cells that are capable of recovering quickly initiate the transcriptional responses immediately, whereas cells that require longer time to adapt also respond later. This is reflected by Hog1 nuclear localization, Hog1 promoter association and the transcriptional response, but not Hog1 phosphorylation, suggesting that a presently uncharacterized rapid adaptive mechanism precedes the Hog1 nuclear response. Furthermore, we found that the period of Hog1 nuclear residence affects the amplitude of the transcriptional response rather than the spectrum of responsive genes.

A simple mathematical model of adaptation to high osmolarity in yeast.

We present a simple ordinary differential equation (ODE) model of the adaptive response to an osmotic shock in the yeast Saccharomyces cerevisiae. The model consists of two main components. First, a biophysical model describing how the cell volume and the turgor pressure are affected by varying extra-cellular osmolarity. The second component describes how the cell controls the biophysical system in order to keep turgor pressure, or equivalently volume, constant. This is done by adjusting the glycerol production and the glycerol outflow from the cell. The complete model consists of 4 ODEs, 3 algebraic equations and 10 parameters. The parameters are constrained from various literature sources and estimated from new and previously published absolute time series data on intra-cellular and total glycerol. The qualitative behaviour of the model has been successfully tested on data from other genetically modified strains as well as data for different input signals. Compared to a previous detailed model of osmoregulation, the main strength of our model is its lower complexity, contributing to a better understanding of osmoregulation by focusing on relationships which are obscured in the more detailed model. Besides, the low complexity makes it possible to obtain more reliable parameter estimates.

Conditional osmotic stress in yeast: a system to study transport through aquaglyceroporins and osmostress signaling.

The accumulation and transport of solutes are hallmarks of osmoadaptation. In this study we have employed the inability of the Saccharomyces cerevisiae gpd1Delta gpd2Delta mutant both to produce glycerol and to adapt to high osmolarity to study solute transport through aquaglyceroporins and the control of osmostress-induced signaling. High levels of different polyols, including glycerol, inhibited growth of the gpd1Delta gpd2Delta mutant. This growth inhibition was suppressed by expression of the hyperactive allele Fps1-Delta1 of the osmogated yeast aquaglyceroporin, Fps1. The degree of suppression correlated with the relative rate of transport of the different polyols tested. Transport studies in secretory vesicles confirmed that Fps1-Delta1 transports polyols at increased rates compared with wild type Fps1. Importantly, wild type Fps1 and Fps1-Delta1 showed similarly low permeability for water. The growth defect on polyols in the gpd1Delta gpd2Delta mutant was also suppressed by expression of a heterologous aquaglyceroporin, rat AQP9. We surmised that this suppression was due to polyol influx, causing the cells to passively adapt to the stress. Indeed, when aquaglyceroporin-expressing gpd1Delta gpd2Delta mutants were treated with glycerol, xylitol, or sorbitol, the osmosensing HOG pathway was activated, and the period of activation correlated with the apparent rate of polyol uptake. This observation supports the notion that deactivation of the HOG pathway is closely coupled to osmotic adaptation. Taken together, our "conditional" osmotic stress system facilitates studies on aquaglyceroporin function and reveals features of the osmosensing and signaling system.

Anaerobicity prepares Saccharomyces cerevisiae cells for faster adaptation to osmotic shock.

Yeast cells adapt to hyperosmotic shock by accumulating glycerol and altering expression of hundreds of genes. This transcriptional response of Saccharomyces cerevisiae to osmotic shock encompasses genes whose products are implicated in protection from oxidative damage. We addressed the question of whether osmotic shock caused oxidative stress. Osmotic shock did not result in the generation of detectable levels of reactive oxygen species (ROS). To preclude any generation of ROS, osmotic shock treatments were performed in anaerobic cultures. Global gene expression response profiles were compared by employing a novel two-dimensional cluster analysis. The transcriptional profiles following osmotic shock under anaerobic and aerobic conditions were qualitatively very similar. In particular, it appeared that expression of the oxidative stress genes was stimulated upon osmotic shock even if there was no apparent need for their function. Interestingly, cells adapted to osmotic shock much more rapidly under anaerobiosis, and the signaling as well as the transcriptional response was clearly attenuated under these conditions. This more rapid adaptation is due to an enhanced glycerol production capacity in anaerobic cells, which is caused by the need for glycerol production in redox balancing. Artificially enhanced glycerol production led to an attenuated response even under aerobic conditions. These observations demonstrate the crucial role of glycerol accumulation and turgor recovery in determining the period of osmotic shock-induced signaling and the profile of cellular adaptation to osmotic shock.