HepG2 cells undergo regulatory volume decrease by mechanically induced efflux of water and solutes

This study challenges the conventional belief that animal cell membranes lack a significant hydrostatic gradient, particularly under anisotonic conditions, as demonstrated in the human hepatoma cell line HepG2. The Boyle van't Hoff (BvH) relation describes volumetric equilibration to anisotonic conditions for many cells. However, the BvH relation is simple and does not include many cellular components such as the cytoskeleton and actin cortex, mechanosensitive channels, and ion pumps. Here we present alternative models that account for mechanical resistance to volumetric expansion, solute leakage, and active ion pumping. We found the BvH relation works well to describe hypertonic volume equilibration but not hypotonic volume equilibration. After anisotonic exposure and return isotonic conditions cell volumes were smaller than their initial isotonic volume, indicating solutes had leaked out of the cell during swelling. Finally, we observed HepG2 cells undergo regulatory volume decrease at both 20 °C and 4 °C, indicating regulatory volume decrease to be a relatively passive phenomenon and not driven by ion pumps. We determined the turgor-leak model, which accounts for mechanical resistance and solute leakage, best fits the observations found in the suite of experiments performed, while other models were rejected.

Meta-analysis of the Boyle van 't Hoff relation: Turgor and leak models explain non-ideal volume equilibrium

There has been much recent attention paid to the interaction of cell volume, its regulation, and the molecular biology of the cell. Cells are generally assumed to behave as linear osmometers, with their water volume linearly proportionate to the inverse of osmotic pressure as described by the Boyle van 't Hoff (BvH) relation. This study evaluates the generality of this and other long-standing assumptions about cell responses to anisotonic conditions. We present alternative models that account for osmoregulation including mechanical resistance to volumetric expansion (the turgor model) and ion-osmolyte leakage (the leak model). To evaluate the generality of the BvH relation and determine the suitability of alternative models, we performed a comprehensive survey of the literature and a careful analysis of the resulting data, and then we used these data to compare among models. We identified 137 articles published from 1964 to 2019 spanning 14 animal species and 26 cell types and determined the BvH relation is not an appropriate general model but is adequate when restricted to the hypertonic region. In contrast, models that account for either mechanical resistance or ion-osmolyte leakage fit well to almost all collected data. The leak model has fitted parameters that are in the same range as the current literature estimate, while the turgor model typically requires an elastic modulus value of one or multiple orders of magnitude larger than literature values. However, confirmation of the underlying mechanism of osmotic regulation is required at the cell-specific level and cannot be assumed a priori.

Modeling Calcium Signaling in S. cerevisiae Highlights the Role and Regulation of the Calmodulin-Calcineurin Pathway in Response to Hypotonic Shock

Calcium homeostasis and signaling processes in Saccharomyces cerevisiae, as well as in any eukaryotic organism, depend on various transporters and channels located on both the plasma and intracellular membranes. The activity of these proteins is regulated by a number of feedback mechanisms that act through the calmodulin-calcineurin pathway. When exposed to hypotonic shock (HTS), yeast cells respond with an increased cytosolic calcium transient, which seems to be conditioned by the opening of stretch-activated channels. To better understand the role of each channel and transporter involved in the generation and recovery of the calcium transient—and of their feedback regulations—we defined and analyzed a mathematical model of the calcium signaling response to HTS in yeast cells. The model was validated by comparing the simulation outcomes with calcium concentration variations before and during the HTS response, which were observed experimentally in both wild-type and mutant strains. Our results show that calcium normally enters the cell through the High Affinity Calcium influx System and mechanosensitive channels. The increase of the plasma membrane tension, caused by HTS, boosts the opening probability of mechanosensitive channels. This event causes a sudden calcium pulse that is rapidly dissipated by the activity of the vacuolar transporter Pmc1. According to model simulations, the role of another vacuolar transporter, Vcx1, is instead marginal, unless calcineurin is inhibited or removed. Our results also suggest that the mechanosensitive channels are subject to a calcium-dependent feedback inhibition, possibly involving calmodulin. Noteworthy, the model predictions are in accordance with literature results concerning some aspects of calcium homeostasis and signaling that were not specifically addressed within the model itself, suggesting that it actually depicts all the main cellular components and interactions that constitute the HTS calcium pathway, and thus can correctly reproduce the shaping of the calcium signature by calmodulin- and calcineurin-dependent complex regulations. The model predictions also allowed to provide an interpretation of different regulatory schemes involved in calcium handling in both wild-type and mutants yeast strains. The model could be easily extended to represent different calcium signals in other eukaryotic cells.

Control of nuclear size by osmotic forces in Schizosaccharomyces pombe

The size of the nucleus scales robustly with cell size so that the nuclear-to-cell volume ratio (N/C ratio) is maintained during cell growth in many cell types. The mechanism responsible for this scaling remains mysterious. Previous studies have established that the N/C ratio is not determined by DNA amount but is instead influenced by factors such as nuclear envelope mechanics and nuclear transport. Here, we developed a quantitative model for nuclear size control based upon colloid osmotic pressure and tested key predictions in the fission yeast Schizosaccharomyces pombe. This model posits that the N/C ratio is determined by the numbers of macromolecules in the nucleoplasm and cytoplasm. Osmotic shift experiments showed that the fission yeast nucleus behaves as an ideal osmometer whose volume is primarily dictated by osmotic forces. Inhibition of nuclear export caused accumulation of macromolecules in the nucleoplasm, leading to nuclear swelling. We further demonstrated that the N/C ratio is maintained by a homeostasis mechanism based upon synthesis of macromolecules during growth. These studies demonstrate the functions of colloid osmotic pressure in intracellular organization and size control.

Microfluidic measurement of individual cell membrane water permeability

This paper reports a microfluidic lab-on-chip for dynamic particle sizing and real time individual cell membrane permeability measurements. To achieve this, the device measures the impedance change of individual cells or particles at up to ten time points after mixing with different media, e.g. dimethyl sulfoxide or DI water, from separate inlets. These measurements are enabled by ten gold electrode pairs spread across a 20 mm long microchannel. The device measures impedance values within 0.26 s after mixing with other media, has a detection throughput of 150 samples/second, measures impedance values at all ten electrodes at this rate, and allows tracking of individual cell volume changes caused by cell osmosis in anisosmotic fluids over a 1.3 s postmixing timespan, facilitating accurate individual cell estimates of water permeability. The design and testing were performed using yeast cells (Saccharomyces cerevisiae). The relationship between volume and impedance in both polystyrene calibration beads as well as the volume-osmolality relationship in yeast were demonstrated. Moreover, we present the first noninvasive and non-optically-based water permeability measurements in individual cells.

Yeast two- and three-species hybrids and high-sugar fermentation

The dominating strains of most sugar-based natural and industrial fermentations either belong to Saccharomyces cerevisiae and Saccharomyces uvarum or are their chimeric derivatives. Osmotolerance is an essential trait of these strains for industrial applications in which typically high concentrations of sugars are used. As the ability of the cells to cope with the hyperosmotic stress is under polygenic control, significant improvement can be expected from concerted modification of the activity of multiple genes or from creating new genomes harbouring positive alleles of strains of two or more species. In this review, the application of the methods of intergeneric and interspecies hybridization to fitness improvement of strains used under high-sugar fermentation conditions is discussed. By protoplast fusion and heterospecific mating, hybrids can be obtained that outperform the parental strains in certain technological parameters including osmotolerance. Spontaneous postzygotic genome evolution during mitotic propagation (GARMi) and meiosis after the breakdown of the sterility barrier by loss of MAT heterozygosity (GARMe) can be exploited for further improvement. Both processes result in derivatives of chimeric genomes, some of which can be superior both to the parental strains and to the hybrid. Three-species hybridization represents further perspectives.

The yeast osmostress response is carbon source dependent

Adaptation to altered osmotic conditions is a fundamental property of living cells and has been studied in detail in the yeast Saccharomyces cerevisiae. Yeast cells accumulate glycerol as compatible solute, controlled at different levels by the High Osmolarity Glycerol (HOG) response pathway. Up to now, essentially all osmostress studies in yeast have been performed with glucose as carbon and energy source, which is metabolised by glycolysis with glycerol as a by-product. Here we investigated the response of yeast to osmotic stress when yeast is respiring ethanol as carbon and energy source. Remarkably, yeast cells do not accumulate glycerol under these conditions and it appears that trehalose may partly take over the role as compatible solute. The HOG pathway is activated in very much the same way as during growth on glucose and is also required for osmotic adaptation. Slower volume recovery was observed in ethanol-grown cells as compared to glucose-grown cells. Dependence on key regulators as well as the global gene expression profile were similar in many ways to those previously observed in glucose-grown cells. However, there are indications that cells re-arrange redox-metabolism when respiration is hampered under osmostress, a feature that could not be observed in glucose-grown cells.

Systems Level Analysis of the Yeast Osmo-Stat

Adaptation is an important property of living organisms enabling them to cope with environmental stress and maintaining homeostasis. Adaptation is mediated by signaling pathways responding to different stimuli. Those signaling pathways might communicate in order to orchestrate the cellular response to multiple simultaneous stimuli, a phenomenon called crosstalk. Here, we investigate possible mechanisms of crosstalk between the High Osmolarity Glycerol (HOG) and the Cell Wall Integrity (CWI) pathways in yeast, which mediate adaptation to hyper- and hypo-osmotic challenges, respectively. We combine ensemble modeling with experimental investigations to test in quantitative terms different hypotheses about the crosstalk of the HOG and the CWI pathways. Our analyses indicate that for the conditions studied i) the CWI pathway activation employs an adaptive mechanism with a variable volume-dependent threshold, in contrast to the HOG pathway, whose activation relies on a fixed volume-dependent threshold, ii) there is no or little direct crosstalk between the HOG and CWI pathways, and iii) its mainly the HOG alone mediating adaptation of cellular osmotic pressure for both hyper- as well as hypo-osmotic stress. Thus, by iteratively combining mathematical modeling with experimentation we achieved a better understanding of regulatory mechanisms of yeast osmo-homeostasis and formulated new hypotheses about osmo-sensing.

Global stability and exact solution of an arbitrary-solute nonlinear cellular mass transport system

The prediction of the cellular state as a function of extracellular concentrations and temperatures has been of interest to physiologists for nearly a century. One of the most widely used models in the field is one where mass flux is linearly proportional to the concentration difference across the membrane. These fluxes define a nonlinear differential equation system for the intracellular state, which when coupled with appropriate initial conditions, define the intracellular state as a function of the extracellular concentrations of both permeating and nonpermeating solutes. Here we take advantage of a reparametrization scheme to extend existing stability results to a more general setting and to a develop analytical solutions to this model for an arbitrary number of extracellular solutes.

Mechanosensitive channels Msy1 and Msy2 are required for maintaining organelle integrity upon hypoosmotic shock in Schizosaccharomyces pombe

The mechanosensitive channels, Mys1 and Msy2, in fission yeast are localized in the endoplasmic reticulum membrane and control cytoplasmic Ca(2+) levels in the hypoosmotic response. We here investigated changes in organellar structures with hypoosmotic shock using transmission electron microscopy. While msy1(-) and msy2(-) single mutant cells developed a number of swollen vacuoles following hypoosmotic shock, similar to wild-type cells, msy1(-) msy2(-) double mutant cells only had two abnormally large vacuoles and cracks between the inner and outer nuclear membranes. These results suggest that Msy1 and Msy2 may be involved in maintaining vacuole integrity and protecting the nuclear envelope upon hypoosmotic shock and also that these two channels are functionally complementary.

Nested autoinhibitory feedbacks alter the resistance of homeostatic adaptive biochemical networks

Negative feedback control is a ubiquitous feature of biochemical systems, as is time delay between a signal and its response. Negative feedback in conjunction with time delay can lead to oscillations. In a cellular context, it might be beneficial to mitigate oscillatory behaviour to avoid recurring stress situations. This can be achieved by increasing the distance between the parameters of the system and certain thresholds, beyond which oscillations occur. This distance has been termed resistance. Here, we prove that in a generic three-dimensional negative feedback system the resistance of the system is modified by nested autoinhibitory feedbacks. Our system features negative feedbacks through both input-inhibition as well as output-activation, a signalling component with mass conservation and perfect adaptation. We show that these features render the system applicable to biological data, exemplified by the high osmolarity glycerol system in yeast and the mammalian p53 system. Output-activation is better supported by data than input-inhibition and also shows distinguished properties with respect to the system's stimulus. Our general approach might be useful in designing synthetic systems in which oscillations can be tuned by synthetic autoinhibitory feedbacks.

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.

Osmotic changes are common environmental challenges for cells, even in multi-cellular organisms, having led to sophisticated adaptation mechanisms. In order to adapt to hyperosmotic stress, yeast cells accumulate glycerol. This is achieved by short-term responses involving metabolic and transmembrane transport changes as well as long-term transcriptional responses. By integrating experimentation and simulation of a mathematical model we resolve the quantitative and temporal characteristics of different processes contributing to glycerol accumulation. We show that osmoadaptation prioritizes the redox and energy balance in glycolysis while rerouting flux from biomass to glycerol production. We further show that the glycerol accumulation network provides osmoadaptation with robustness by compensating for the loss of certain nodes and with the flexibility necessary for responding to different stress situations. Finally we provide novel insight into the roles of transport processes in glycerol accumulation and evidence that trehalose may play a role in yeast osmoadaptation. The present work provides for the first time an integrated dynamic view on a eukaryotic osmolyte system and links signaling with regulation of gene expression and metabolic control.

Modelling reveals novel roles of two parallel signalling pathways and homeostatic feedbacks in yeast

Ensemble modelling is used to study the yeast high osmolarity glycerol (HOG) pathway, a prototype for eukaryotic mitogen-activated kinase signalling systems. The best fit model provides new insights into the function of this system, some of which are then experimentally validated.

The main mechanism for osmo-adaptation is a fast and transient non-transcriptional Hog1-mediated activation of glycerol production.

The transcriptional response rather serves to maintain an increased steady-state glycerol production with low steady-state Hog1 activity after adaptation.

A fast negative feedback of activated Hog1 on the upstream signalling branches serves to stabilise the adaptation response by preventing oscillatory behaviour.

Two parallel redundant signalling branches elicit a more robust and swifter adaptation than a single branch alone, at least for low osmotic shock. This notion could be corroborated by dedicated measurements of single-cell volume recovery for the wild-type and single-branch mutants.

The high osmolarity glycerol (HOG) pathway in yeast serves as a prototype signalling system for eukaryotes. We used an unprecedented amount of data to parameterise 192 models capturing different hypotheses about molecular mechanisms underlying osmo-adaptation and selected a best approximating model. This model implied novel mechanisms regulating osmo-adaptation in yeast. The model suggested that (i) the main mechanism for osmo-adaptation is a fast and transient non-transcriptional Hog1-mediated activation of glycerol production, (ii) the transcriptional response serves to maintain an increased steady-state glycerol production with low steady-state Hog1 activity, and (iii) fast negative feedbacks of activated Hog1 on upstream signalling branches serves to stabilise adaptation response. The best approximating model also indicated that homoeostatic adaptive systems with two parallel redundant signalling branches show a more robust and faster response than single-branch systems. We corroborated this notion to a large extent by dedicated measurements of volume recovery in single cells. Our study also demonstrates that systematically testing a model ensemble against data has the potential to achieve a better and unbiased understanding of molecular mechanisms.

Quantification of cell volume changes upon hyperosmotic stress in Saccharomyces cerevisiae

Cell volume is a biophysical property, which is of great importance for quantitative characterisations of biological processes, such as osmotic adaptation. It also is a crucial parameter in the most common type of mathematical description of cellular behaviour-ordinary differential equation (ODE) models, e.g. the integrative model of the osmotic stress response in baker's yeast (E. Klipp, B. Nordlander, R. Kruger, P. Gennemark and S. Hohmann, Nat. Biotechnol., 2005, 23, 975-982). Until recently only rough estimates of this value were available. In this study we measured the mean volume of more than 300 individual yeast cells (Saccharomyces cerevisiae). We quantitatively characterised the dependence between the relative cell volume and the concentration of osmoticum in the cell surrounding. We also followed the recovery of the cellular volume over time, as well as the influence of increased external osmolarity on the nuclear volume. We found that cell shrinkage caused by shifts in the external osmolarity is proportional to the stress intensity only up to 1000 mM NaCl. At this concentration the yeast cells shrink to approximately 55% of their unstressed volume and this volume is maintained even in the case of further osmolarity increase. We observed that returning to the initial, unstressed volume takes more than 45 minutes for stress concentrations exceeding 100 mM NaCl and that only cells treated with the latter concentration are able to fully regain their initial size within the course of the experiment. We postulate that the cytoplasm plays a protective role for the nucleus by buffering the changes in volume caused by external osmolarity shifts. In conclusion, we quantitatively characterised the dynamics of cell volume changes caused by hyperosmotic stress, providing an accurate description of a biophysical cell property, which is crucial for precise mathematical simulations of cellular processes.

Biophysical properties of Saccharomyces cerevisiae and their relationship with HOG pathway activation

Parameterized models of biophysical and mechanical cell properties are important for predictive mathematical modeling of cellular processes. The concepts of turgor, cell wall elasticity, osmotically active volume, and intracellular osmolarity have been investigated for decades, but a consistent rigorous parameterization of these concepts is lacking. Here, we subjected several data sets of minimum volume measurements in yeast obtained after hyper-osmotic shock to a thermodynamic modeling framework. We estimated parameters for several relevant biophysical cell properties and tested alternative hypotheses about these concepts using a model discrimination approach. In accordance with previous reports, we estimated an average initial turgor of 0.6 ± 0.2 MPa and found that turgor becomes negligible at a relative volume of 93.3 ± 6.3% corresponding to an osmotic shock of 0.4 ± 0.2 Osm/l. At high stress levels (4 Osm/l), plasmolysis may occur. We found that the volumetric elastic modulus, a measure of cell wall elasticity, is 14.3 ± 10.4 MPa. Our model discrimination analysis suggests that other thermodynamic quantities affecting the intracellular water potential, for example the matrix potential, can be neglected under physiological conditions. The parameterized turgor models showed that activation of the osmosensing high osmolarity glycerol (HOG) signaling pathway correlates with turgor loss in a 1:1 relationship. This finding suggests that mechanical properties of the membrane trigger HOG pathway activation, which can be represented and quantitatively modeled by turgor.