|Potassium starvation in yeast: mechanisms of homeostasis revealed by mathematical modeling.|
|Jump to Full Text|
|PMID: 22737060 Owner: NLM Status: MEDLINE|
|The intrinsic ability of cells to adapt to a wide range of environmental conditions is a fundamental process required for survival. Potassium is the most abundant cation in living cells and is required for essential cellular processes, including the regulation of cell volume, pH and protein synthesis. Yeast cells can grow from low micromolar to molar potassium concentrations and utilize sophisticated control mechanisms to keep the internal potassium concentration in a viable range. We developed a mathematical model for Saccharomyces cerevisiae to explore the complex interplay between biophysical forces and molecular regulation facilitating potassium homeostasis. By using a novel inference method ("the reverse tracking algorithm") we predicted and then verified experimentally that the main regulators under conditions of potassium starvation are proton fluxes responding to changes of potassium concentrations. In contrast to the prevailing view, we show that regulation of the main potassium transport systems (Trk1,2 and Nha1) in the plasma membrane is not sufficient to achieve homeostasis.|
|Matthias Kahm; Clara Navarrete; Vicent Llopis-Torregrosa; Rito Herrera; Lina Barreto; Lynne Yenush; Joaquin Ariño; Jose Ramos; Maik Kschischo|
|Type: Journal Article; Research Support, Non-U.S. Gov't; Validation Studies Date: 2012-06-21|
|Title: PLoS computational biology Volume: 8 ISSN: 1553-7358 ISO Abbreviation: PLoS Comput. Biol. Publication Date: 2012|
|Created Date: 2012-06-27 Completed Date: 2012-10-25 Revised Date: 2013-07-12|
Medline Journal Info:
|Nlm Unique ID: 101238922 Medline TA: PLoS Comput Biol Country: United States|
|Languages: eng Pagination: e1002548 Citation Subset: IM|
|Department of Mathematics and Technology, RheinAhrCampus, University of Applied Sciences, Koblenz, Remagen, Germany.|
|APA/MLA Format Download EndNote Download BibTex|
Cation Transport Proteins / metabolism
Cell Membrane / metabolism
Potassium / metabolism*
Proton-Translocating ATPases / genetics, metabolism
Saccharomyces cerevisiae / genetics, metabolism*
Saccharomyces cerevisiae Proteins / genetics, metabolism
Sodium-Hydrogen Antiporter / metabolism
|0/Cation Transport Proteins; 0/NHA1 protein, S cerevisiae; 0/Saccharomyces cerevisiae Proteins; 0/Sodium-Hydrogen Antiporter; 0/TRK2 protein, S cerevisiae; 136956-54-2/TRK1 protein, S cerevisiae; 7440-09-7/Potassium; EC 3.6.1.-/PMA1 protein, S cerevisiae; EC 220.127.116.11/Proton-Translocating ATPases|
Journal ID (nlm-ta): PLoS Comput Biol
Journal ID (iso-abbrev): PLoS Comput. Biol
Journal ID (publisher-id): plos
Journal ID (pmc): ploscomp
Publisher: Public Library of Science, San Francisco, USA
Kahm et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Received Day: 5 Month: 12 Year: 2011
Accepted Day: 23 Month: 4 Year: 2012
collection publication date: Month: 6 Year: 2012
Print publication date: Month: 6 Year: 2012
Electronic publication date: Day: 21 Month: 6 Year: 2012
Volume: 8 Issue: 6
E-location ID: e1002548
PubMed Id: 22737060
Publisher Id: PCOMPBIOL-D-11-01834
|Potassium Starvation in Yeast: Mechanisms of Homeostasis Revealed by Mathematical Modeling Alternate Title:Potassium Starvation in Yeast|
1Department of Mathematics and Technology, RheinAhrCampus, University of Applied Sciences, Koblenz, Remagen, Germany
2Department of Microbiology, Campus de Rabanales, University of Córdoba, Córdoba, Spain
3Instituto de Biologia Molecular y Celular de Plantas UPV-CSIC, Ciudad Politécnica de la Innovación, Universidad Politécnica de Valencia, Valencia, Spain
4Institut de Biotecnologia I Biomedicina & Department of Biochemistry and Molecular Biology, Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Barcelona, Spain
|Stockholm University, Sweden
|Correspondence: * E-mail: firstname.lastname@example.org
Contributed by footnote: Conceived and designed the experiments: M. Kahm, L. Yenush, J. Ariño, J. Ramos, M. Kschischo. Performed the experiments: C. Navarrete, V. Llopis-Torregrosa, R. Herrera, L. Barreto, J. Ramos. Analyzed the data: M. Kahm, J. Ariño, J. Ramos, M. Kschischo. Contributed reagents/materials/analysis tools: M. Kahm, L. Yenush, J. Ariño, J. Ramos, M. Kschischo. Wrote the paper: M. Kahm, L. Yenush, J. Ariño, J. Ramos, M. Kschischo. Developed the mathematical model: M. Kahm, M. Kschischo.
Potassium is an essential cation required for many cellular processes including the regulation of cell volume, intracellular pH, protein synthesis, activation of enzymes, and maintenance of the plasma membrane potential –. In their natural environment, most cell types have to accumulate intracellular potassium against a strong concentration gradient. Animal cells utilize the energy stored in ATP to directly pump potassium ions into the cell via the / ATPase. This enzyme is absent in most fungi and plants , which have developed alternative mechanisms to control the intracellular potassium concentration. Saccharomyces cerevisiae (S.c.) cells can grow in media with a potassium concentration ranging from to . Despite extensive knowledge about the identity and function of most potassium transporters in this organism , a systems level understanding of the interplay and regulation of the various transport pathways is still lacking.
In S.c., uptake of potassium across the plasma membrane is driven by the membrane potential, which itself is generated by proton pumping via the -ATPase, Pma1 , . The high affinity and high velocity transporter, Trk1, is the major uptake system for potassium. The expression levels of the other Trk protein, Trk2, are low, compared to Trk1, and therefore considered of minor importance , . A low affinity uptake observed by electrophysiological techniques in trk1,2 double mutants has been attributed to the putative calcium blocked channel Nsc1, though the gene responsible for this transport activity has not been found yet , . Efflux of potassium is strongly pH-dependent and coupled to sodium toxicity. The antiporter Nha1 extrudes or ions in exchange for protons under acidic environmental conditions and contributes to the continuous cyclic flux of potassium ions across the plasma membrane and to pH regulation , . It is only at higher external pH that potassium or sodium is actively extruded by the Ena1 ATPase –. Another potassium efflux system is the voltage gated channel, Tok1. Electrophysiological studies revealed that Tok1 opens at positive membrane potentials, which do not occur under normal physiological conditions . Potassium is also stored in intracellular compartments, in particular in the vacuole. The effect of intracellular transport is, however, not sufficiently characterized yet , .
Besides protons, a number of other ions are associated with the transport of potassium. The anion bicarbonate was shown to be important for potassium accumulation . Decarboxylation reactions produce carbon dioxide, which is quickly converted to carbonic acid ( ), by carbonic anhydrase. Carbonic acid can either diffuse freely across the cell membrane or dissociate into bicarbonate ( ), and protons. While protons can be extruded via Pma1, the permeability of bicarbonate is very low compared to that of carbonic acid. The resulting accumulation of bicarbonate provides the link to potassium homeostasis; the negative charges carried by bicarbonate can be balanced by potassium cations. In principle, other weak acids could contribute in a similar way to potassium accumulation, but our results below and previous investigations suggest that the bicarbonate reaction plays an important role . Potassium transport is also related to ammonium toxicity . Under low external potassium conditions, ammonium leaks into the cells, presumably via potassium transporters. Toxic concentrations of ammonium are counteracted by increased production and excretion of amino acids .
The maintenance of a minimal potassium concentration requires the orchestration of the different transport systems under the constraints of various thermodynamic forces. In this article, we use a mathematical model in conjunction with a novel inference algorithm (the reverse tracking algorithm) and model-driven experimentation to identify the key transport mechanisms that must be regulated under the conditions of potassium shortage. We show that the activation of the proton pump, Pma1, and the activation of the bicarbonate reaction sequence are the regulators of potassium homeostasis. We also show that potassium homeostasis is an example of non-perfect adaptation: The intracellular potassium concentration depends on the external potassium concentration and is only regulated to keep minimal levels of potassium required for survival. This is different from other homeostatic systems such as osmoregulation , where certain stationary systems characteristics perfectly adapt, irrespective of the external conditions.
To study the response of S.c. cells to an abrupt decrease of external potassium, we performed potassium starvation experiments using and free media. Cells grown in non-limiting potassium ( KCl) were washed with -free YNB medium (YNB without amino acids and ammonium sulphate, Formedium UK, CYN7505 plus 2% glucose, traces of KCl left: 15 , hereafter referred to as Translucent -free medium ) and resuspended in the same medium . The time course for changes in intracellular potassium concentrations for the wild type strain exhibits two different phases (Figure 1A). In the first hour of starvation there is a large net efflux of potassium indicated by the rapid decrease in the intracellular concentration. Loss of potassium slows down in the second phase and the internal concentration slowly approaches a new stationary state (Table 1). Although the cells cannot perfectly adapt to the large concentration gradients they are able to keep a certain amount of potassium required for survival (approx. ). Interestingly, the second phase of potassium loss is slower for the trk1,2 double mutant than for the wild type (wt). This is surprising, because it is believed  that increased uptake of potassium via Trk1 even at very low external potassium concentrations is a major mechanism of potassium homeostasis. Thus, one would have expected the concentration of internal potassium in the trk1,2 mutant to be lower than in the wild type. The time course for the nha1 mutant is not significantly different from the trk1,2 mutant (see also Figure S7 in Text S1).
Multiple signaling pathways modulate the activity of the various transport systems involved in potassium homeostasis –, , , –. However, it is not entirely clear which of these signals are essential to achieve homeostasis and how they are acting under the constraints set by the thermodynamics of ion transport. To study these constraints, we developed a minimalistic mathematical model which incorporates the essential parts known to be important for potassium homeostasis. The model describes the dynamic coupling between the intracellular potassium concentration , internal pH ( ), carbon dioxide concentration , membrane voltage , and cell volume . A complete description of the equations and parameter values is given in the Materials and Methods section and derivations can be found in the Text S1. Here, only the basic model structure is given:(Figure 1B). The model comprises the Trk1,2 system (abstracted as a single system, ), the Nha1 antiporter ( ), and the Tok1 channel ( ). To mimic the joint contribution of other, mainly non-specific transport pathways for potassium (e.g. Nsc1) we added a potassium leak current to the model. The Ena1 ATPase is neglected because it is known to be inactive at the relatively low external pH used in the experiments .
The dynamics of carbon dioxide (Equation (2)) is coupled to the transport fluxes of bicarbonate and carbonic acid . These transport rates are given in the Materials and Methods section (Equations (18–19)) and a detailed derivation of the bicarbonate model  is given in the Text S1. Carbon dioxide is produced in various metabolic processes such as the TCA cycle or pyruvate decarboxylation. It is impossible to model all these processes explicitly, but we incorporate them in the effective metabolic carbon dioxide production flux . This flux is an input to the model and was initially assumed to be constant.
The change in pH (Equation (3)) per change in proton concentration is described by the buffering capacity . In principle, is a function of the internal pH, but due to the combined action of various buffering species  it can be approximated by a constant for a wide range of intracellular pH values. In addition to the proton fluxes via the -ATPase Pma1 ( ) and the Nha1 antiporter ( ) there are many other proton transport pathways in yeast. The corresponding net flux is subsumed in the proton leak current . The effective proton flux originating from the bicarbonate reaction sequence is given by the term , where is the pH-dependent fraction of dissociated carbon dioxide.
The membrane potential (Equation (4)) is modeled as a charge balance equation ( , specific membrane capacitance; , Faraday constant; , surface area of the cell) . We explicitly modeled the charges carried by potassium, total protons ( ) and bicarbonate. The remaining net charges contributing to the membrane potential are subsumed in , which is determined by the initial conditions of the dynamic variables in the model.
The cell volume (Equation (5)) depends on the balance between internal osmotic pressure , external osmotic pressure and turgor pressure . Ion transport processes change the intra- and extracellular solute concentrations and thus have an osmotic effect (Equations (24–26)) in Materials and Methods). The resistance against volume changes is given by the hydraulic permeability parameter .
The concentration and voltage dependent kinetics of all transport systems were described by simple thermodynamic consistent relationships. The driving force for the transport fluxes of ions across the plasma membrane can be written as the difference of the membrane potential and the equilibrium potential . The equilibrium potential depends on the concentrations and stoichiometry of the ions transported, see Equations (12–14) in the Materials and Methods section. For the potassium fluxes in Equation (1) and the proton leak in Equation (3) we assumed linear relations (Ohm's law) of the form between the driving force and the transport flux , or the corresponding electrical current , respectively. For the leak currents and we initially assumed constant conductivity parameters (Equations (9) and (11) in Materials and Methods). The conductivity of the transport proteins Trk1,2, Nha1 and Tok1 was modeled as a function of the membrane voltage, see Equations (6–8) in Materials and Methods.
This minimalistic model captures the essential biophysical and thermodynamic constraints under which control of potassium homeostasis operates. Despite the simplicity of the model, the experimental data was not sufficient to uniquely identify all the parameters. We decided to use this flexibility to explore the parameter space for regions that are consistent with the data and performed extensive parameter scans and sensitivity analysis simulations. However, we were unable to identify a single parameter combination which reproduced the experimental time courses for the wild type strain observed in Figure 1A. In the model, all potassium inside the cell was rapidly and completely lost upon starvation (Figure S4A in Text S1). Based on our model simulations, this believed to be caused by a strong efflux via the Nha1 antiporter driven by the large concentration gradient across the plasma membrane. This model behavior is robust against various model variations, including the incorporation of an intracellular potassium storage mechanism that mimics the contribution of intracellular compartments to potassium retention. Thus, we conclude that further dynamic mechanisms counteracting the strong potassium gradient are essential for homeostasis. Importantly, the model described so far incorporates only the biophysics of transport but does not account for gene regulatory, signal transduction or metabolic events affecting the transporter activity.
The fact that the minimal model is not able to reproduce the experimental time courses for potassium starvation means that there are some unmodeled dynamics that are not captured by the model. Under the working hypothesis that the model covers the major biophysical effects of potassium transport we assumed that there are additional regulatory responses to a shortage of potassium. Available knowledge ,  and data is currently not sufficient to develop exhaustive models for the metabolic, signal transduction and gene regulatory responses to potassium starvation. It is not even clear which of the transporters or other components are activated or deactivated for the maintenance of homeostasis. In engineering terms , neither the regulators nor the signals triggering their action are sufficiently characterized.
To overcome this limitation, we combined our minimal biophysical model with an inference algorithm for unmodeled dynamics. We assumed that the unknown regulatory events modulate the activity of the transport systems or other components in the model. Mathematically this means that a constant parameter in the model might in fact not be constant, but a function of time. For example, the maximum conductivity (see Equation 6) of the Trk1,2 transport system could be influenced by signal transduction events ,  in response to low potassium. Any attempt to explicitly model this regulation by additional equations is hindered by insufficient knowledge of the structure and dynamics of the regulatory networks involved. However, one might recoin the question and ask: “Is there a function such that the given experimental time course of intracellular potassium and the time course predicted by the model are in sufficient agreement?”. If such a function would exist we would regard the modulation of the Trk1,2 transporter as one potential regulatory mechanism and Trk1,2 as a potential regulator of potassium homeostasis. However, there might be another parameter (e.g. ) associated with a transporter or another component in the model for which a time course exists such that experimental data can be reproduced. Our strategy was now to test different parameters and corresponding processes for being potential regulators, see Figure 2A. We define a transporter or any other component in the model to be a potential regulator if a tracking control signal exists which changes the activity of the component in such a way that the experimental time course and stationary data can be reproduced. We refer to this inference approach as the reverse tracking algorithm, a more detailed mathematical explanation is given in Materials and Methods.
We used the reverse tracking algorithm to test the transporters Trk1,2, Ena1, Nha1, Tok1 and Pma1 and the activity of the bicarbonate reaction for being potential regulators of potassium homeostasis and then compared the predicted tracking control signals to experimental observations. There is no tracking control signal for the major uptake system Trk1,2; see Figure S1 in the Text S1. This is in contrast to the prevailing view that increased uptake of potassium via Trk1 is essential for potassium homeostasis under starvation conditions. The loss of potassium after starvation is slower in trk1,2 double mutants (see Figure 1A) than in wild type cells. It was experimentally observed  that these double mutants have a more negative membrane potential than wild type cells under starvation conditions and also when external potassium is plentiful. This stronger membrane potential (see also Figure S3 in the Text S1) counteracts the outwardly directed potassium concentration gradient and thus explains the higher potassium levels after starvation. Taken together, these results show that the uptake of potassium via Trk1,2 is not the primary mechanism to prevent excessive loss of potassium under starvation conditions.
Although we found a tracking control signal for the Nha1 antiporter, we excluded it from our list of potential regulators based on two observations. First, as indicated in Figure 1A, the time course of potassium loss in nha1 mutants is slower than in the wild type and similar to the trk1,2 mutants. Secondly, it was demonstrated that the influence of Nha1 on the internal potassium concentrations decreases with time . This is in contradiction to our predicted tracking signal (Figure S2 in Text S1), which is nonmonotonic in time.
Similarly, the unspecific transport pathways (leak currents) were excluded, because it is not plausible that unspecific transporters are regulated for the specific purpose of potassium homeostasis. This is based on the well founded assumption that all potassium specific transporters are active under our experimental conditions are known  and included in the model. The non-specific cation uptake system NSC1 can be excluded, because our medium contains enough calcium to render NSC1 inactive . The proton flux includes many co-transport mechanisms with nutrients and other molecules. It is thus unlikely, that one of these transport mechanisms is specifically regulated in response to potassium starvation.
The remaining parts in our model are the Pma1 -ATPase and the bicarbonate reaction sequence. For both of them, the reverse tracking approach predicts a rapid burst of activity in response to the rapid removal of external potassium (Figure 2B and 2C). Activation of proton pumping by Pma1 (Figure 2B) hyperpolarizes the plasma membrane, which counteracts the large concentration gradient of potassium and thus limits potassium efflux. An increased reaction flux (see Figure 2C) through the bicarbonate system has a similar effect: The negative charges carried by bicarbonate increase the magnitude of the membrane potential and thereby compensate the potassium gradient.
To test the prediction that Pma1 is activated after potassium starvation, we measured Pma1 activity from crude membrane preparations  using an in vitro method that has been extensively established as a faithful measure of in vivo Pma1 function , , . Indeed, the activity measurements confirm the prediction of the reverse tracking algorithm that Pma1 activity increases rapidly (timescale of 10 minutes) and slowly declines during the first hours of potassium starvation (Figure 2D). Control experiments revealed that Pma1 protein levels do not change under these conditions. Moreover, we also observe, as predicted by the model, that the Pma1 activity is higher in the trk1,2 mutant than in the wild type strain throughout the time course of potassium starvation (Figures 2B and D). To further substantiate that the activation of Pma1 is essential for the response to low potassium, we measured growth for Pma1 mutants pma1–204 and pma1–205 with decreased expression and ATPase activity (33 and 50% of wild type). Figure 3 shows that the ratio of the growth rates at 1 mM and 50 mM external potassium is much lower for the mutant strains than that of the wild type. These results are in line with the recent finding that the brp1 mutant, which is a PMA1 promotor deletion, that leads to decreased Pma1 protein levels, presents markedly decreased growth in low potassium medium and defective rubidium uptake .
The second prediction from the reverse tracking approach was an increased reaction flux for the bicarbonate system (Figure 2C). This prediction is supported experimentally by an increased mRNA expression of the NCE103 gene coding for carbonic anhydrase, the enzyme catalyzing the bicarbonate reaction (Figure 2E). This result was part of a genome-wide transcriptomic analysis, using DNA microarrays, of the response to potassium starvation (0–120 min) to be published elsewhere (Barreto et al., submitted). It was shown earlier that protein and mRNA levels of carbonic anhydrase are highly correlated . A qRT-PCR measurement confirmed the increase in NCE103 expression in wild type cells shifted to free medium. After 60 minutes of potassium starvation, the NCE103 mRNA levels increase more than four-fold ( independent experiments). These results show that activation of both Pma1 and the bicarbonate reaction sequence are essential for the control of internal potassium concentrations. In non-starved cells the expression of NCE103 is higher for the trk1,2 double mutant than for the wild type (single dot in Figure 2E). A confirmatory semi-quantitative RT-PCR measurement using the same RNA sample as in the microarray experiment and one RNA sample from independent cultures yielded a mean expression ratio of ( data points) for trk1,2 relative to the wild type. These results suggest a simple explanation for the reported hyperpolarization of the trk1,2 double mutant : A high activity of the bicarbonate reaction sequence means that many protons and many bicarbonate ions are produced. Together with a more active proton pump (Figures 2B and D), this results in a more negative membrane potential that counteracts the outwardly directed potassium gradient. The consequence is a higher intracellular concentration of potassium (Figure 1A) in trk1,2 double mutants than wild type cells.
Homeostatic control of a cellular function in response to a changing environment is often mediated by a negative feedback loop. A change in the input signal (e.g. the external potassium concentration) is counteracted by this feedback loop in order to keep an essential cellular quantity (e.g. the intracellular potassium concentration) in a range sufficient for the cell's function. One particular type of feedback is integral control, where the control signal is the time integral of the difference between the reference and the actual quantity . Integral control was observed for a number of cellular processes including bacterial chemotaxis ,  and osmoregulation . A characteristic property of integral control is perfect adaptation, where the steady state input is independent of the steady state output. For potassium this would mean, that the same intracellular potassium concentration (output) is approached irrespective of the extracellular potassium concentration (input).
The activation of proton transport by Pma1 and the activation of the bicarbonate system counteracting low external potassium indicate the existence of a negative feedback loop. To further investigate this feedback, we have modified the potassium starvation experiment. As before, cells were grown at external potassium, but now resuspended in media with different external potassium concentrations. The potassium efflux and the stationary internal concentrations are different for the different external concentrations, which is also reflected by the model (Figure 4A). To test whether these stationary intracellular concentrations are characteristic for the external potassium, we grew cells overnight in media with different external potassium concentrations (Figure 4B). When external potassium is plentiful ( ), the internal concentration attains an upper limit of approx. . For low external potassium ( ), the internal concentration is proportional to the external and agrees with the stationary states of Figure 4A. These experiments show, that perfect adaptation by integral control is not a characteristic of potassium homeostasis for low external potassium. The molecular function and characteristics of this feedback have to be further explored.
In summary, we found that direct regulation of potassium transport proteins is not sufficient for the maintenance of viable potassium levels inside the cell. Although the presence of Trk1,2 influences the dynamics of potassium loss under conditions of low potassium, the regulation of their activity is not the main regulatory process. Cells lacking these proteins have higher intracellular potassium concentrations and the loss of potassium after a rapid shift to low external potassium is slower than in wild type cells. The adaptation to low potassium requires a rapid modulation of proton fluxes as a rescue operation via the increased production of bicarbonate and the activation of the -ATPase Pma1 (Figure 5). The observation that the internal steady state potassium concentration is determined by the external concentration indicates, that potassium homeostasis is an example of non-perfect adaptation, excluding the existence of integral control. The detailed sensing and signaling mechanisms remain to be elucidated and currently we cannot distinguish whether changes in internal or external potassium are sensed directly or indirectly, e.g., as changes of the membrane potential.
Although we cannot completely rule out the possibility that other transport systems not considered in the model contribute to homeostasis, we have reason to believe that our model covers the dominant effects required for the maintenance of viable potassium levels under starvation conditions. All experiments were performed in the presence of calcium, which renders the activity of the calcium blocked non-selective cation pathway Nsc1 unlikely. In addition, non-specific transport of potassium is covered in the model by the leak current. The information about potassium storage in intracellular compartments in the literature is limited. To test the influence of intracellular potassium fluxes originating from an intracellular storage mechanism, we added a hypothetical compartment which can release potassium in response to starvation. This modification did not change the qualitative behavior of the model and was not sufficient to explain the slow efflux of potassium and the maintenance of sufficient intracellular potassium after starvation. Thus, we excluded this modification from the model.
Many cation transporters are evolutionarily conserved in other yeast species and even in higher plants –. However, the current knowledge for these organisms is not as detailed. Considering the importance of ion homeostasis for some pathogenic yeasts  and for the growth and development of plants, the question of whether the regulation of proton fluxes plays a similar dominant role as in S.c. is an interesting starting point for future research.
The development of dynamic mathematical models requires a compilation of all parts and processes which could potentially be important for a cellular mechanism under consideration. Other processes believed to be negligible are often lumped together in the parameter values of the model. The decision of which processes to incorporate or to neglect is often hampered by insufficient biological knowledge. Incorporating too many details is impractical and leads to overly complex models with many parameters and little predictive power. On the other extreme are simplistic models which potentially neglect important processes and cannot reproduce the experimental data. We believe that our strategy to start with such a minimal model and to infer unmodeled dynamics with a reverse tracking approach might be of broader interest in systems biology. The reverse tracking algorithm provides (i) candidate points of applications for regulatory signals not explicitly captured by the model and (ii) an estimate of the corresponding time dependent regulatory signal. We emphasize that these potential regulatory signals have to be checked for biological plausibility and have to be validated by experiments. It can be applied when the core model for the process of interest is “underfitted”, i.e. when it can not sufficiently reproduce the experimental data because other regulatory process influence the parameters in the model. Its main advantage is that it can be applied even when an explicit modeling of the processes generating these regulatory inputs is beyond reach. On the other hand, the algorithm can be used as a tool for prioritizing experiments. In combination with experiments, it also may also help to indicate which model extensions are most promising.
The basic structure of the mathematical model is given by Equations (1–5) in the Results section. Here we report the details of the kinetic relationships. Parameter values, initial conditions and derivations are provided in the Text S1. In the following , and denote the Faraday constant, the gas constant and the temperature.
The concentration and voltage dependent transport currents of the passive or secondary active potassium and proton transporters (Equations (6–11) below) are described by the Ohmian relation , where the conductivity is either constant ( and ) or a function of the membrane voltage (Trk1,2, Nha1, Tok1). The voltage dependent conductivity of the transporter was described by the function, , . The parameters and for Tok1 were taken from the literature ; for Trk1,2 we estimated them from electrophysiological data ; see Figure S6 in Text S1. For Nha1 we assumed this voltage-dependent conductivity on the basis of . This approach leads to the following transport kinetics (here ) are related to the mass fluxes by was introduced to correctly incorporate the dependence of the conductance parameters on the surface area of the cell.
The transport current for Pma1 depends on the free energy of ATP hydrolysis and was modeled as 
The model in  for the bicarbonate reaction sequence was supplemented by an effective metabolic carbon dioxide production or consumption flux . The production rate is an input to the model. Note that the production changes with the volume even if the amount of produced CO does not change. This is taken into account by the relationshipdenotes the initial volume and is the volume independent rate, which is input to the model.
The dynamics of CO in Equation (2) depends also on the transport fluxes for carbonic acid and bicarbonateand are permeabilities for carbonic acid and bicarbonate respectively, and is the volume to surface ratio of the cell, see Text S1. The assumption behind Equation (18) is that the flux of carbonic acid is proportional to its concentration gradient. Equation (19) is a Goldman-Hodgkin-Katz flux equation for the electrodiffusive transport of ions across the membrane , .
The fraction of total undissociated carbon dioxideof carbonic anhydrase) enters the rate of proton production in Equation (3), see Text S1 for a detailed derivation. The production fluxes of carbonic acid and bicarbonate and on the production rate of carbon dioxide .
For the dynamics of the cell volume (Equation (5)) we used a model from the literature  which is based on the balance of internal and external osmotic pressure ( and ) and turgor pressure ( ):is the fraction of the non-osmotic cell volume and determines the elasticity of the cell wall. The internal and external concentration of other osmotically active substances are constant (see Text S1). A cell volume simulation for wildtype and trk1,2 mutant can be found in Figure S5 in Text S1.
Equations (1–5) have the form of a differential algebraic control system. Here, denotes the dynamical variables ( ) and Equation (4) for the membrane potential corresponds to the algebraic equation . The scalar input function is given by the external potassium concentration . The solution of this system for given values of the parameters and a given input function is denoted by . Assume now, that we can observe of the components of experimentally. We collect the experimentally observable components in . This can be written as with a matrix with binary elements . For an experimentally observable variable the -th column of has a single entry and . A zero column with indicates that can not experimentally be observed and is thus excluded from .
Assume further, that we have experimental data for certain time points in response to the known input function . Most parameter estimation techniques aim to minimize the squared errorin order to bring the model prediction for a given input close to the experimental data . However, it might be the case that the minimum error is still too large so that the model cannot be regarded as a reasonable description of the data. This could mean that a dynamical process not explicitly accounted for renders at least one component of the parameter vector to be a time dependent function instead of being constant. The reverse tracking algorithm aims (i) to identify, which of the components of are potentially time dependent and (ii) to predict the time course which minimizes the error. Although the unmodeled dynamical process might effect more than one component, we consider for simplicity each component separately and solve the problem separately. Here, denotes the parameter vector with the l-th component excluded. We then regard as a potential regulatory input, if the problem (29) has a solution with a minimum error smaller than a predefined threshold : . There might be more than one potential regulatory input and the decision of which of these are real can only be made from biological considerations or from additional validation experiments. For example, it might be that has a huge magnitude or takes unrealistic values which could be used to exclude from the list of potential regulatory inputs.
Mathematically, problem (29) is an optimal tracking problem, which often can be solved by a feedback control law . This means that the function is updated according to the local error at time . For a scalar we found the integral controller is a least squares spline fit to the experimental data points . Details and suitable parameter values for are provided in the Text S1.
Cells were grown in Translucent -free medium supplemented with the indicated amount of KCl to an OD600 of 0.4–0.6. Intracellular potassium concentrations were measured by atomic emission spectrometry after extracting the cells with acid as previously described . The time course of internal potassium was obtained by growing the cells in KCl, then cells were washed with Translucent -free medium (traces of KCl left: ) and resuspended to the same free medium or containing the indicated KCl concentrations. Apart from the washing procedure the medium contains 2% glucose.
Data for NCE103 expression changes upon potassium starvation was obtained in the context of a genome-wide transcriptomic analysis by DNA microarray (Barreto et al., Manuscript submitted). Microarray data has been deposited at NCBI's Gene Expression Omnibus  and are accessible through GEO Series accession numbers GSE24711 (trk1 trk2 data) and GSE24712 (time-course data). Briefly, wild-type strain BY4741 cells were grown in Translucent medium supplemented with 50 mM KCl to OD 0.8. Cells were centrifuged and resuspended either in fresh Translucent medium with 50 mM KCl or without potassium. Samples (20 ml) were taken at 10, 20, 40, 60 and 120 min by rapid filtration from 4 biological replicates. Total RNA was extracted by using the Ribo PureTM Yeast kit (Ambion) following the manufacturers instructions. cDNA was prepared and indirectly labeled with Cy3 and Cy5. Images with a resolution of 10 were analyzed with the GenePix Pro 6.0 software (Molecular Devices).
Microarray data was confirmed by qRT-PCR using independent RNA samples. To this end, 60 ng of RNA were amplified using oligonucleotides RT_ NCE103_5 ([gene: TCATTACCTGTCGCACTG]) and RT_ NCE103_3 ([gene: CACAAAAGTTACCCCAAAA]) and the QuantiTect SYBR Green PCR Kit (Quiagen).
Cell cultures were grown at in Translucent YNB medium containing KCl to OD660 0.6, then washed with Translucent - free medium and resuspended in the same medium without KCl. At the indicated times, cell samples were pelleted by centrifugation, resuspended in of fresh media (with KCl for t = 0 and without KCl for the remaining samples), incubated for 5 minutes and frozen in liquid nitrogen. For the crude membrane purification, of 3× extraction buffer (0.3 M Tris-HCl pH 8.0, 180 mM KCl, 30 mM EDTA, 6 mM DTT and Protease Inhibitor Cocktail (Roche)) was added to the thawed samples and cells were broken by vortexing in the presence of an equal volume of glass beads. of GTED20 buffer (20% glycerol, 10 mM Tris-HCl pH 7.6, 1 mM EDTA and 1 mM DTT) were added to the crude extract, which was then centrifuged 5 minutes at 2000 rpm. The supernatant was transferred to a new tube and centrifuged 20 minutes at 13000 rpm. The insoluble fraction was resuspended and homogenized in of GTED20. The total amount of protein present was estimated using the Bradford assay (BioRad). The amount of Pma1 present in this protein fraction was estimated by comparing the amount of Pma1 to a protein standard curve separated in SDS-PAGE gels stained with Coomassie Blue. In a microtiter plate, of total protein (which corresponds to of Pma1) were assayed in the presence and absence of a Pma1-specific inhibitor, dietilstylbestrol (DES, final concentration 0.2 mM). The reaction was started by adding of the reaction buffer (50 mM MES-Tris pH 5.7, 5 mM MgSO4, 50 mM , 5 mM Na Azide, 0.3 mM Molybdate, 2 mM ATP) and the plate was incubated for 20 minutes at . The reaction was stopped by adding of detection solution (2% sulphuric acid, 0.5% ammonium molybdate, 0.5% SDS, 0.1% ascorbic acid) and the color was allowed to develop for 5 minutes before reading the absorbance in microplate reader (BioRad) at 750 nm. Residual activity values in the presence of DES were subtracted from the absolute activity values to obtain the Pma1 activity measurements. The results represent the average of at least 4 measurements at each time point and essentially identical results were observed in two separate experiments. Measurements of Pma1 activity are expressed in mmol/min/g Pma1. Error bars represent the standard deviation.
Supporting text containing additional information about the mathematical model.
Click here for additional data file (pcbi.1002548.s001.pdf)
The authors have declared that no competing interests exist.
Maik Kschischo and Matthias Kahm were supported by BMBF grant 0315786C (SysMo2/Translucent 2). Work in Joaquin Ariño laboratory was supported by grants BFU2008-04188-C03-01, BFU2011-30197-C3-01, GEN2006-27748-C2-1-E/SYS (SysMo/Translucent) and EUI2009-04147 (SysMo2/Translucent 2), from the MICINN, Spain. Joaquin Ariño was the recipient of an Ajut 2009SGR-1091 and an ICREA Academia Award (Generalitat de Catalunya). Work in José Ramos laboratory was supported by grants GEN2006-27748-C2-2-E/SYS (SysMo/Translucent), BFU2008-04188-C03-03 and EUI2009-04153 (SysMo2/Translucent 2), (MICINN, Spain). Work in Lynne Yenush laboratory was supported by grants BFU2008-04188-C03-02, BFU2011-30197-C03-03 (MICINN, Spain) and ACOMP/2011/024 (Generalitat Valenciana). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors thank Hana Sychrova, Hella Lichtenberg, Joachim Ernst, Jost Ludwig, Wolfgang Alt and Milan Höfer for many useful discussions.
|1.||Blatt M,Slayman C. Year: 1987Role of active potassium transport in the regulation of cytoplasmic pH by nonanimal cells.Proc Natl Acad Sci USA84273727413472234|
|2.||Rodriguez-Navarro A. Year: 2000Potassium transport in fungi and plants.Biochim Biophys Acta146913010692635|
|3.||Ariño J,Ramos J,Sychrova H. Year: 2010Alkali-metal-cation transport and homeostasis in yeasts.Microbiol Mol Biol Rev749512020197501|
|4.||Merchan S,Pedelini L,Hueso G,Calzada A,Serrano R. Year: 2011Genetic alterations leading to increases in internal potassium concentrations are detrimental for DNA integrity in Saccharomyces cerevisiae.Genes Cells1615216521143561|
|5.||Buch-Pedersen MJ,Rudashevskaya EL,Berner TS,Venema K,Palmgren MG. Year: 2006Potassium as an intrinsic uncoupler of the plasma membrane H+-ATPase.J Biol Chem281382853829217056603|
|6.||Serrano R. Year: 1983In vivo glucose activation of the yeast plasma membrane ATPase.FEBS Lett15611146221943|
|7.||Ko CH,Buckley AM,Gaber RF. Year: 1990TRK2 is required for low affinity K+ transport in Saccha romyces cerevisiae.Genetics1253053122199312|
|8.||Ramos J,Alijo R,Haro R,Rodriguez-Navarro A. Year: 1994TRK2 is not a low-affinity potassium transporter in Saccharomyces cerevisiae.J Bacteriol1762492528282703|
|9.||Bihler H,Slayman CL,Bertl A. Year: 1998NSC1: A novel high-current inward rectifier for cations in the plasma membrane of Saccharomyces cerevisiae.FEBS Lett43259649710251|
|10.||Roberts SK,Fischer M,Dixon GK,Sanders D. Year: 1999Divalent cation block of inward currents and low-affinity K+ uptake in Saccharomyces cerevisiae.J Bacteriol1812912979864342|
|11.||Bañuelos M,Ruiz M,Jiménez A,Souciet JL,Potier S. Year: 2002Role of the Nha1 antiporter in regulating K+ inux in Saccharomyces cerevisiae.Yeast1991511754478|
|12.||Navarrete C,Petrezsélyová S,Barreto L,Martínez J,Zahrádka J. Year: 2010Lack of main K+ uptake systems in S. cerevisiae cells affects yeast cell physiological parameters both in potassium sufficient and limiting conditions.FEMS Yeast Res1050851720491939|
|13.||Haro R,Garciadeblas B,Rodríguez-Navarro A. Year: 1991A novel P-type ATPase from yeast involved in sodium transport.FEBS Lett2911891911657642|
|14.||Márquez JA,Serrano R. Year: 1996Multiple transduction pathways regulate the sodium-extrusion gene PMR2/ENA1 during salt stress in yeast.FEBS Lett38289928612770|
|15.||Ruiz A,Ariño J. Year: 2007Function and regulation of the Saccharomyces cerevisiae ENA sodium ATPase system.Eukaryot Cell62175218317951516|
|16.||Bertl A,Slayman CL,Gradmann D. Year: 1993Gating and conductance in an outward-rectifying K+ channel from the plasma membrane of Saccharomyces cerevisiae.J Membr Biol1321831998492306|
|17.||Martinez-Muñoz GA,Peña A. Year: 2005In situ study of K+ transport into the vacuole of Saccha romyces cerevisiae.Yeast2268970416034802|
|18.||Lopéz R,Enríquez E,Peña A. Year: 1999Effects of weak acids on cation accumulation, ΔpH and ΔΨ in yeast.Yeast1555356210341418|
|19.||Hess DC,Lu W,Rabinowitz JD,Botstein D. Year: 2006Ammonium toxicity and potassium limitation in yeast.PLoS Biol420122023|
|20.||Muzzey D,Gómez-Uribe CA,Mettetal JT,van Oudenaarden A. Year: 2009A systems-level analysis of perfect adaptation in yeast osmoregulation.Cell13816017119596242|
|21.||Petrezsélyová S,Zahradka J,Sychrová H. Year: 2010Saccharomyces cerevisiae BY4741 and W303-1A laboratory strains differ in salt tolerance.Fungal Biol11414415020960970|
|22.||Casado C,Yenush L,Melero C,del Carmen Ruiz M,Serrano R. Year: 2010Regulation of Trk dependent potassium transport by the calcineurin pathway involves the Hal5 kinase.FEBS Lett5842415242020412803|
|23.||Mulet JM,Leube MP,Kron SJ,Rios G,Fink GR. Year: 1999A novel mechanism of ion home ostasis and salt tolerance in yeast: the Hal4 and Hal5 protein kinases modulate the Trk1-Trk2 potassium transporter.Mol Cell Biol193328333710207057|
|24.||Yenush L,Merchan S,Holmes J,Serrano R. Year: 2005pH-responsive, posttranslational regulation of the Trk1 potassium transporter by the type 1-related Ppz1 phosphatase.Mol Cell Biol258683869216166647|
|25.||Boron W,de Weer P. Year: 1976Intracellular pH transients in squid giant axons caused by CO2, NH3 and metabolic inhibitors.J Gen Physiol67911121460|
|26.||Boron WF. Year: 2004Regulation of intracellular pH.Adv Physiol Educ2816017915545345|
|27.||Keener J,Sneyd J. Year: 2004Mathematical PhysiologyNew YorkSpringer|
|28.||Gennemark P,Nordlander B,Hohmann S,Wedlin D. Year: 2006A simple mathematical model of adaption to high osmolarity in yeast.In Silico Biol634|
|29.||Weiss TF. Year: 1996Cellular Biophysics - Transport.Cambridge: Bradford|
|30.||Åström KJ,Murray RM. Year: 2009Feedback SystemsPrinceton and OxfordPrinceton University Press|
|31.||Yenush L,Mulet JM,Ariño J,Serrano R. Year: 2002The Ppz protein phosphatases are key regulators of K+ and pH homeostasis: implications for salt tolerance, cell wall integrity and cell cycle progression.EMBO J2192092911867520|
|32.||Bihler H,Slayman CL,Bertl A. Year: 2002Low-affinity potassium uptake by Saccharomyces cerevisiae is mediated by NSC1, a calcium-blocked non-specific cation channel.Biochim Biophys Acta155810911811779561|
|33.||Serrano R. Year: 1978Characterization of the plasma membrane ATPase of Saccharomyces cerevisiae.Mol Cell Biochem225163154059|
|34.||Portillo F,Serrano R. Year: 1988Dissection of functional domains of the yeast proton-pumping ATPase by directed mutagenesis.EMBO J7179317982901955|
|35.||Vallejo CG,Serrano R. Year: 1989Physiology of mutants with reduced expression of plasma membrane H+-ATPase.Yeast53073192528864|
|36.||Barreto L,Canadell D,Petrezsélyová S,Navarrete C,Maresová L. Year: 2011A genomewide screen for tolerance to cationic drugs reveals genes important for potassium homeostasis in Sac charomyces cerevisiae.Eukaryot Cell101241125021724935|
|37.||Amoroso G,Morell-Avrahov L,Müller D,Klug K,Sültemeyer D. Year: 2005The gene NCE103 (YNL036w) from Saccharomyces cerevisiae encodes a functional carbonic anhydrase and its tran scription is regulated by the concentration of inorganic carbon in the medium.Mol Microbiol5654955815813743|
|38.||Barkai N,Leibler S. Year: 1997Robustness in simple biochemical networks.Nature3879139179202124|
|39.||Yi TM,Huang Y,Simon MI,Doyle J. Year: 2000Robust perfect adaptation in bacterial chemotaxis through integral feedback control.Proc Natl Acad Sci U S A974649465310781070|
|40.||Miranda M,Bashi E,Vylkova S,Edgerton M,Slayman C. Year: 2009Conservation and dispersion of sequence and function in fungal TRK potassium transporters: focus on Candida albicans.FEMS Yeast Res927829219175416|
|41.||Hille B. Year: 2001Ion channels of excitable membranes.Sunderland: Sinauer|
|42.||Johansson I,Blatt M. Year: 2006Interactive domains between pore loops of the yeast K+ channel Tok1 associate with extracellular K+ sensitivity.Biochem J39364565516287426|
|43.||Kuroda T,Bihler H,Bashi E,Slayman C,Rivetta A. Year: 2004Chloride channel function in the yeast Trk-potassium transporters.J Membr Biol19817719215216418|
|44.||Ohgaki R,Nakamura N,Mitsui K,Kanazawa H. Year: 2005Characterization of the ion transport activity of the budding yeast Na+/H+ antiporter, Nha1p, using isolated secretory vesicles.Biochim Biophys Acta171218519615950597|
|45.||Endresen LP,Hall K,Hoye JS,Myrheim J. Year: 2000A theory for the membrane potential of living cells.Eur Biophys J299010310877018|
|46.||Xian B,Dawson DM,de Queiroz MS,Chen J. Year: 2003A continuous asymptotic tracking control strategy for uncertain multi-input nonlinear systems.5257 In: Proceedings of Intelligent Control 2003 IEEE Int Symp; 8 October 2003. doi:10.1109/ISIC.2003.1253913.|
|47.||Edgar R,Domrachev M,Lash AE. Year: 2002Gene expression omnibus: Ncbi gene expression and hybridization array data repository.Nucleic Acids Res3020721011752295|
Table 1 Optical densities during starvation.
|Strain||Time of starvation (hours)|
Optical densities for the wild type and the trk1,2 double mutant corresponding to the potassium starvation experiments of Figure 1A.
Previous Document: Modelling cell polarization driven by synthetic spatially graded Rac activation.
Next Document: Protein scaffolds can enhance the bistability of multisite phosphorylation systems.