 Research article
 Open Access
 Published:
Modeling mutant phenotypes and oscillatory dynamics in the Saccharomyces cerevisiae cAMPPKA pathway
BMC Systems Biology volume 7, Article number: 40 (2013)
Abstract
Background
The cyclic AMPProtein Kinase A (cAMPPKA) pathway is an evolutionarily conserved signal transduction mechanism that regulates cellular growth and differentiation in animals and fungi. We present a mathematical model that recapitulates the shortterm and longterm dynamics of this pathway in the budding yeast, Saccharomyces cerevisiae. Our model is aimed at recapitulating the dynamics of cAMP signaling for wildtype cells as well as single (pde1 Δ and pde2 Δ) and double (pde1 Δ pde2 Δ) phosphodiesterase mutants.
Results
Our model focuses on PKAmediated negative feedback on the activity of phosphodiesterases and the Ras branch of the cAMPPKA pathway. We show that both of these types of negative feedback are required to reproduce the wildtype signaling behavior that occurs on both short and long time scales, as well as the the observed responses of phosphodiesterase mutants. A novel feature of our model is that, for a wide range of parameters, it predicts that intracellular cAMP concentrations should exhibit decaying oscillatory dynamics in their approach to steady state following glucose stimulation. Experimental measurements of cAMP levels in two genetic backgrounds of S. cerevisiae confirmed the presence of decaying cAMP oscillations as predicted by the model.
Conclusions
Our model of the cAMPPKA pathway provides new insights into how yeast respond to alterations in their nutrient environment. Because the model has both predictive and explanatory power it will serve as a foundation for future mathematical and experimental studies of this important signaling network.
Background
In eukaryotic cells, the cyclic adenosine monophosphate (cAMP) – Protein Kinase A (PKA) pathway plays a central role in mediating diverse biological responses such as growth, development, and cell differentiation [1, 2]. Stimuli such as hormones, neurotransmitters, nutrients, and physiological stress agents trigger signaling cascades that collectively mount a response through cAMPmediated PKA signaling [3–9]. In the budding yeast, Saccharomyces cerevisiae, activation of PKA in response to essential nutrients and fermentable carbon sources is directed by intracellular levels of cAMP [10]. Synthesis of cAMP from ATP is catalyzed by the enzyme adenylate cyclase (Cyr1) and governed by two different Gprotein systems (Figure 1) [10–12]. Ras2, a small GTPbinding protein, in its GTPbound (active) state stimulates adenylate cyclase, causing a rapid increase in intracellular cAMP levels (Figure 1, top right) [11]. Ras2 activity is positively regulated by the guanine nucleotide exchange factors (GEFs) Cdc25 and Sdc25 as a function of intracellular glucose levels [13–20]. Conversely, GTPase activating proteins, Ira1 and Ira2, down regulate Ras2 activity by stimulating the hydrolysis of GTP [21, 22]. In parallel to Ras, a second Gprotein pathway, involving the proteins Gpr1, Gpa2, and Rgs2, responds to extracellular levels of glucose and increases adenylate cyclase activity (Figure 1, top left). Gpr1 is a membrane bound Gprotein coupled receptor that activates the G α protein Gpa2 in response to extracellular glucose levels. The GTP bound form of Gpa2 in turn stimulates adenylate cyclase activity [12, 20, 23, 24]. cAMP synthesis through this pathway is negatively controlled by Rgs2, the GTPase activating protein of Gpa2 [25, 26]. An increase in cellular cAMP levels upon glucose induction activates cAMPdependent PKA. PKA is a holoenzyme that includes both regulatory (Bcy1) and catalytic subunits (Tpk1, Tpk2 and Tpk3) [27, 28]. cAMP binding to the regulatory subunits leads to the release of the catalytic PKA subunits ([29]; Figure 1, bottom) which are then free to interact with downstream targets such as metabolic enzymes, transcription factors and other kinases [30–33].
Downstream of the adenylate cyclase, cAMP levels are modified by the action of phosphodiesterases, enzymes that catalyze the conversion of cAMP to AMP, thus preventing cAMP accumulation in the cell [34, 35]. In yeast there are two phosphodiesterases, the lowaffinity phosphodiesterase Pde1 and the higher affinity Pde2. Pde2 is a class II phosphodiesterase and shares homology with other eukaryotic phosphodiesterases [36], while Pde1 homologues have been identified in only a small number of eukaryotes [35]. The highaffinity Pde2 is thought to play the key role in maintaining cAMP at steady state levels [37]. The role of Pde1 is less well understood, though genetic studies in both S. cerevisiae and the distantly related yeast Candida albicans suggests that Pde1 is a target of PKA, and that PKAmediated feedback on Pde1 is an important component of maintaining tight regulation of cAMP levels [35, 38].
In this study we present a mathematical model of the cAMPPKA pathway in Saccharomyces cerevisiae. The primary goal of this model is to explain the observed response of wildtype cells and PDE mutants (pde1 Δ, pde2 Δ, and pde1 Δ pde2 Δ) to glucose stimulation. Wildtype cells exhibit a rapid, transient increase in cAMP levels followed by a quick return to a new steady state; this response typically occurs on a time scale of 6090 seconds [35]. pde2 Δ mutants have approximately wildtype dynamics, while pde1 Δ mutants exhibit elevated peak cAMP levels and a slow return to a steady state that is higher than that of wildtype cells. Surprisingly, the double phosphodiesterase mutant, pde1 Δ pde2 Δ, shows no response to glucose stimulus, and maintains cAMP concentrations at essentially a prestimulus steady state [35]. The question we posed was whether it was possible to mathematically reconcile all of these observed dynamics with currently known genetic and biochemical interactions. Our model focuses on PKAmediated negative feedback interactions, and suggests that feedback on both phosphodiesterases and on the Ras branch of the pathway are important for reproducing the dynamics observed in wildtype cells and PDE mutants.
While a number of studies have proposed mathematical models for the cAMPPKA network in yeast [39–41], none of these previous models is capable of recapitulating the observed dynamics of both wildtype cells and the double phosphodiesterase mutant (pde1 Δ pde2 Δ). Our proposed model is thus unique in this aspect. Another notable feature of our model is that, for a wide range of parameters, it predicts the presence of decaying cAMP oscillations following glucose stimulus. Oscillatory behaviors were not an explicit goal of our initial modeling effort, and thus we treated this as a novel model prediction. We tested this prediction by measuring the cAMP response to glucose stimulus in two genetically diverse yeast strains (S288c and Σ 1278b). Our experiments reveal variation in oscillatory behaviors between these genetic backgrounds, and we demonstrate that our model, in turn, is capable of reproducing the diversity of observed oscillatory patterns. Finally, our core model can also be easily extended to generate sustained cAMP oscillations that occur on longer time scales (tens of minutes; [42]). The results presented here thus demonstrate both the descriptive and predictive power of the model we have developed, and suggest that this model can serve as a basis for future experimental and modeling studies of this important eukaryotic signal transduction pathway.
Summary of motivating experiments
There is a large body of experimental work that provides information on genetic and biochemical interactions relevant to the cAMPPKA pathway. The strategy we adopted in this study was to use this body of work as a basis for constraining interactions in our model. In particular our model was developed to understand two specific experimental results: the shortterm dynamics of cAMP following glucose stimulation reported by Ma et al. [35]; and the longterm dynamics of cAMP under various stress levels reported by GarmendiaTorres et al. [42].
Shortterm cAMP dynamics
Ma et al. [35] investigated the shortterm behavior of cAMP signaling following glucose stimulation in cells that had previously been grown under carbonsource limitation. These data therefore represents the dynamics of cAMP signaling as cells transition from starved to fed states, on a time scale of seconds to minutes. Ma et al.’s analysis focused particularly on PKAmediated negative feedback involving the phosphodiesterases Pde1 and Pde2. Figure 2 summarizes the five dynamic patterns we consider here, based on Ma et al.’s experiments involving wildtype cells and various phosphodiesterase mutants (compare to Ma et al. Figures 2a and 3a). These cases are:
Case 1: Wildtype (wt; blue line) – the concentration of cAMP in wt cells rapidly increases following glucose stimulation. After reaching a peak level, the concentration of cAMP declines to a new steady state that is higher than its initial concentration. We refer to this as the wildtype transient response.
Case 2: Pde1 knockout (pde1 Δ; red line) – the concentration of cAMP exhibits a much larger and longer transient response, about four times as large as the transient response exhibited by wildtype.
Case 3: Pde2 knockout (pde2 Δ; green line) – the concentration of cAMP exhibits similar dynamics to wildtype.
Case 4: Pde1 phosphorylation mutant (pde1^{ala 152}; pink line) – the cAMP transient response is about twice as large as that exhibited by wildtype.
Case 5: Double PDE knockout (pde1 Δ pde2 Δ; black line) – glucose stimulation has no effect on the concentration of cAMP present in the cytosol. The initial value in this case is much larger then the wildtype initial value.
These results are surprising in two ways: 1) the dynamics observed in the pde2 Δ mutant are almost the same as those observed in wildtype cells, while the response in the pde1 Δ mutant is drastically different, despite the fact that Pde2 is thought to be the higher affinity phosphodiesterase; and 2) in the double mutant (pde1 Δ pde2 Δ) glucose stimulation has no effect on the concentration of cAMP, indicating a synergistic interaction between the PDE nullalleles.
Longterm cAMP dynamics
GarmendiaTorres et al. [42] focused on the longterm behavior of cAMP signaling in yeast exposed to stress conditions. They predicted sustained oscillations, with a period of 510 minutes, in the concentration of cAMP at intermediate stress levels. This prediction was motivated by oscillations observed in the nucleocytoplasmic shuttling of the transcriptional activator Msn2. Msn2 is negatively regulated by the cAMPPKA pathway through phosphorylation of Msn2 by PKA; this phosphorylation promotes export of Msn2 from the nucleus [43]. GarmendiaTorres et al. predicted that the nucleocytoplasmic shuttling of Msn2 that they observed was driven by oscillations in the concentration of cAMP, and they constructed a mathematical model to explore the proposed dynamics. In their model, high levels of stress drive the concentration of cAMP to a low steady state, and at low stress levels the cAMP levels go to a high steady state. Only under intermediate levels of stress did cAMP oscillations occur. This oscillatory behavior in their model was driven by PKAmediated negative feedback on Ras ·GTP.
Mechanisms of PKAmediated negative feedback on cAMP signaling
Both sets of observations described above suggest that PKA plays an important role in attenuating cAMP signaling via negative feedback interactions. PKAmediated negative feedback on cAMP signaling has been proposed to occur at two levels – via enhancement of phosphodiesterase activity and via attenuation of Ras signaling.
PKA interactions with Pde1 are supported by experiments by Ma et al. that showed that mutagenesis of a putative PKA phosphorylation site in Pde1 (pde1^{ala 152}) causes a dramatic increase in cAMP accumulation (Case 4, Figure 2). There is less direct evidence for PKA interactions with Pde2 but other class II phosphodiesterases are known to be targets of PKA [36]. Hu et al. [44] showed that hyperactive PKA activity leads to elevated Pde2 levels, and that this resulted from increased protein stability of Pde2.
The observation that cAMP signaling did not significantly change in the pde2 Δ mutant might seem to suggest that Pde2 plays no role in regulating cAMP immediately following glucose stimulus. However, the cAMP transient in the pde1 Δ mutant, and the synergistic effect observed in the pde1 Δ pde2 Δ double mutant, argue for a role for Pde2, though with a slower activation rate as compared to Pde1. However, while we favor the hypothesis of PKA feedback on both PDEs, we can not rule out the possibility that the transient cAMP dynamics observed in the pde1 Δ mutant might involve other feedback mechanisms.
With respect to Ras signaling, one or more elements of this pathway may be targets of PKAmediated feedback. For example, Ira1 and Ira2 both have predicted PKA phosphorylation sites [12, 45]. Similarly, PKA hyperphosphorylates Cdc25, causing its dissociation from the adenylate cyclase/Ras2 complex. Cdc25 may also directly bind to the adenylate cyclase to facilitate membrane anchoring [46–48]. Relocalization of Cdc25 to the cytoplasm attenuates the glucose response of adenylate cyclase, inhibiting synthesis of cAMP [49, 50]. Moreover, membrane dissociated Cdc25 can no longer activate Ras2, triggering an increase in Ras ·GDP levels [4]. Ma et al. suggested that the likely explanation for the pattern they observed in the pde1 Δ pde2 Δ double mutant (no change in cAMP levels after glucose stimulation) resulted from constitutively high PKA feedback inhibition in the absence of the phosphodiesterases.
In summary, both studies described above suggest that PKA down regulates cAMP signaling by increasing the rate of PDE mediated decay and/or by decreasing the rate of cAMP production via interactions with the Ras effectors. Thus a key goal of our modeling effort was to explore PKA feedback on phosphodiesterase activity and the Ras pathway. In our model PKA interacts with both PDEs to increase their activity, though we assume PKA has a higher affinity for Pde1 than for Pde2. Similarly, we model negative feedback on Ras ·GTP through PKA phosphorylation of the RasGTPases Ira1 and Ira2. These interactions are highlighted by the red arrows in Figure 1.
Model
To fully describe the key reactions of the cAMPPKA pathway as depicted in Figure 1 requires a system of nine timedependent differential equations to model the concentrations of: 1) Gpa2 ·GTP (active Gpa2); 2) Cdc25; 3) Ras ·GTPases Ira1 and Ira2; 4) Ras ·GTP (active Ras2); 5) adenylate cyclase activity; 6) activated Pde1; 7) activated Pde2; 8) cAMP; and 9) active PKA (free catalytic subunits). The full ninedimensional system is described in detail in Additional file 1: Supplementary model (Equations (S1)(S9)).
For the purposes of exposition, we first describe the motivation behind the equations that represent PKA mediated feeedback. We then describe a simplification of the model to a fourdimensional system that replicates the key dynamical behaviors of cAMP signaling as described above.
Key reactions for PKAmediated feedback
PKA feedback (red arrows in Figure 1) in our model occurs in three ways: (i) PKA phosphorylates Pde1, enhancing Pde1 breakdown of cAMP; (ii) PKA phosphorylates Pde2, enhancing Pde2 breakdown of cAMP; and (iii) PKA phosphorylates the Ras ·GTPases Ira1 and Ira2. In Cases (i) and (ii), PKA inhibits the concentration of cAMP by enhancing the breakdown of cAMP. In Case (iii) PKA inhibits the production of cAMP by inactivating Ras ·GTP.
We assume that there are much greater concentrations of Pde1, Pde2, Ira1 and Ira2 in the cytoplasm than of PKA. We also assume that Pde1, Pde2 and Ira1/2 compete with each other for activation by PKA. We further assume that Ira1 and Ira2 can be simultaneously activated by PKA, thus we treat them as a single variable, Ira. These assumptions are taken in order to understand the Ma et al. cases. The three competitive reactions are specified as follows:
Here the asterisk (^{∗}) indicates the activated or phosphorylated form of the enzyme. Here and below, PKA refers the catalytic subunits (Tpk1, Tpk2, or Tpk3). Equations (1)–(3), below, model the velocities of these reactions using an extension of MichaelisMenten kinetics [51] that allows for competition between Pde1, Pde2 and Ira for PKA.
In order for the model to replicate the dynamics observed by Ma et al. [35] (Figure 2) we impose the following conditions:
Condition (a) The following inequalities must hold:
Condition (b) In comparing analogous reactions of Pde1 and Pde2, the reactions of Pde2 are uniformly slower.
Condition (c) PKA rapidly phosphorylates Ira.
By invoking Condition (a), the effect of Pde2 is negligible when Pde1 is active. This assumption is justified by the observation that the dynamics of cAMP signaling following glucose stimulus are essentially unchanged in the pde2 Δ mutant. By invoking Condition (b) we can model both Pde1 and Pde2 with a single equation (for details and derivation see the Additional file 1: Supplementary Model). Condition (b) is justified by the large transient cAMP signal observed in the pde1 Δ mutant. Neither Condition (a) nor (b) are required to reproduce the behavior of the wildtype or mutant strains, but are mathematically convenient for exploring the model analytically as described below.
Condition (c) reflects negative feedback on the Ras branch of the pathway. We impose Condition (c) so that even when PKA’s phosphorylation of Ira is slowed by the presence of either Pde1 or Pde2, by Condition (a), the effect of PKA on Ira can not be ignored. In the double mutant case (pde1 Δ pde2 Δ) the loss of competitive inhibition of PKA greatly enhances PKA’s effect on Ira. Thus in the double mutant case PKA has a strong negative effect on Ras ·GTP. Several studies have provided experimental support for a role of the Ira proteins in negative feedback on cAMP signaling [21, 22] as well as evidence that this feedback is regulated by PKA [52].
Steady state assumptions
The model may be further simplified by assuming that the following four reactions are fast and hence proceed to steady state.

1.
Gpa2 ·GDP $\rightleftharpoons $ Gpa2 ·GTP (reaction 1 in Figure 1).

2.
Activation/inactivation of Cdc25 (reaction 2 in Figure 1).

3.
Activation/inactivation of adenylate cyclase (reaction 5 in Figure 1).

4.
cAMP+[Bcy1::PKA] $\rightleftharpoons $ [cAMP::Bcy1] +PKA (reaction 9 in Figure 1).
Again, these steady state assumptions are made for mathematical simplicity, and are not a requirement to mathematically reproduce wildtype or mutant signaling dynamics. We will only focus on the last of these steady state assumptions, for the rest we refer the reader to the Additional file 1: Supplementary Model. We assume that PKA activates rapidly, that is we assume that the concentration of active PKA can be treated as being at steady state. This assumption was taken for mathematical simplicity of the model, it is not a necessary condition to numerically replicate the dynamics observed by either Ma et al. or GarmendiaTorres et al. In our model PKA feedback plays a central role, thus we will examine this steady state condition in detail. In the cAMPPKA pathway PKA is activated by cAMP in the following manner: four cAMP molecules bind to two regulatory subunits (Bcy1) and release two catalytic subunits, creating the active form of PKA.
When modeling this reaction we can simplify this in two ways: 1) by assuming that the concentration of active PKA does not approach its maximum; and 2) by approximating the 4 to 2 ratio by a 2 to 1 ratio. Applying this simplification we can model the change in PKA activity as proportional to cAMP squared (x^{2}); thus in Equations (4a)–(4d) (below) x^{2} represents PKA feedback (for full details see the Additional file 1: Supplementary Model).
Simplified model
Combining feedback Conditions (a), (b), and (c) with the four steady state assumptions, and by modeling both forms of PDE with a single case dependent variable, p, we reduce the ninedimensional model shown in Figure 1 to the simpler fourdimensional model depicted in Figure 3. The model in Figure 3 is described by the following system of equations:
In Equations (4a)–(4d) and in Figure 3 the concentration of Ras ·GTP is represented by variable r; the concentration of Ras ·GDP is represented by (1−r); the concentrations of Ira1 and Ira2 are represented by z; the variable p represents the concentration of Pde1 in Cases 1, 3, and 4 (wt, pde2 Δ, and pde1^{ala 152}) and the concentration of Pde2 in Case 2 (pde1 Δ); the concentration of cAMP is represented by x; as stated above x^{2} represents the concentration of active PKA; and time is given by variable t. All variables have been nondimensionalized. The constants in Equations (4a)–(4d) are nondimensional composites of dimensional parameters of the full system, derived in the Additional file 1: Supplementary Model. Intuitively, A is the activation rate of Ras ·GTP catalyzed by Cdc25; B is the inactivation rate of Ras ·GTP catalyzed by Ira1 and Ira2; Γ_{1} is the affinity of both Cdc25 and Ira1/2 for Ras ·GTP; N is the reaction rate of Ira1 and Ira2; M is the quantitative expression of Condition (b); C is the production rate of cAMP due to the basal activity of adenylate cyclase; G accounts for the glucose, normalized so that G=1 after a glucose stimulus is applied; D_{0} represents a “basal” decay rate of cAMP in the absence of activated Pde; D represents enhanced decay due to PKA feedback on Pde modified by Γ, the affinity of activated Pde for cAMP.
Results
In this section we examine analytical, numerical, and experimental results motivated by our model of the cAMPPKA pathway. We show that the simplified version of the model, given by Equations (4a)–(4d), can adequately replicate the shortterm dynamics of cAMP, reported by Ma et al. [35] for wildtype cells and all three phosphodiesterase mutants (pde1 Δ, pde2 Δ, pde1 Δ pde2 Δ). For the wildtype case we analyze the behavior of the system as it approaches steady state and find that for a wide range of parameters the model predicts that cAMP levels should exhibit decaying oscillations. Since oscillatory behavior was not an explict input into our modeling effort, we considered this a novel prediction, and undertook an experimental validation of the model using time series measurements of cAMP signaling for two diploid yeast strains. As predicted by the mathematical model we observed decaying cAMP oscillations following glucose stimulation. We then extended the model to consider longertime scale, sustained cAMP oscillations (510 minute periods) in response to nutrient stress, as predicted by GarmendiaTorres et al. [42]. We show that the core model can also reproduce sustained oscillations on this longer time scale, and we discuss the model parameters and the corresponding biochemical interactions, that are required to generate such sustained oscillations.
Shortterm dynamics
Our first task is to demonstrate that Equations (4a)–(4d) can adequately replicate the dynamics of cAMP signaling immediately following glucose stimulation, for both wildtype cells and mutants lacking either one or both phosphodiesterases. Figure 4 shows the modeled concentration of cAMP as a function of time, following glucose stimulation of wildtype and PDE mutants. This outcome was generated using the parameter choices given in Tables 1 and 2. The dynamical patterns generated by the model are a good match to those illustrated in Figure 2. The dimensional concentrations in the figure were obtained by multiplying x in Equations (4a)–(4d), the dimensionless variable representing cAMP concentrations, by 24.95 fmol·10^{−6} cells; and multiplying time, given by dimensionless variable t, by.038 minutes. These scale factors emerge from the discussion in the Additional file 1: Supplementary Model. Initial conditions were the steadystate concentrations that occur under glucose starved conditions.
Note that the numbers in columns of Table 1 that represent the wildtype and pde2 Δ backgrounds are identical. Because of Condition (a) above, when Pde1 is present, the effect of Pde2 is negligible. Thus, the values in these columns of the table characterize properties of Pde1. Similarly, the values in the pde1 Δ column characterize properties of Pde2. In the case of pde1^{ala 152} both the activation and inactivation rate of Pde1 decrease, but the affinity of Pde1 for cAMP (Γ) does not change. The value of B is increased in the double PDE knockout because, by Condition (a), when both forms of PDE are eliminated the effect of PKA on Ira1/2 is greatly enhanced. As mentioned above, the dimensionless parameters in Equations (4a)–(4d) are composites of dimensional parameters. In choosing dimensional parameters we used values given by GarmendiaTorres et al. [42] when available. Moreover, the parameters determined by fitting the model also lie in biologically appropriate ranges. We refer to the Additional file 1: Supplementary Model for the details of the fitting process.
An intuitive explanation for the behavior of wildtype cells is that the transient cAMP peak is due to the delay between the increase of adenylate cyclase activity and PKA mediated enhancement of phosphodiesterase activity. Activation of adenylate cyclase activity increases cAMP concentrations, leading to the release of PKA catalytic subunits. PKA in turn activates the phosphodiesterases, which enhances the decay of cAMP and brings the concentration of cAMP down to its steady state level. In wild type cells, we propose that Pde1 is the primary effector of PKA feedback. In the pde1 Δ and pde1^{ala 152} cases the slower activation of Pde2 accounts for the greater transient observed in these cases.
Oscillatory approaches to steady state
For a wide range of parameters the model described by Equations (4a)–(4d) suggests that cAMP concentration should exhibit decaying oscillations as it approaches its steady state value following glucose stimulus. Here we describe the conditions required to generate this behavior.
Decaying cAMP oscillations via feedback on PDEs
PKA feedback through the PDEs alone is sufficient to generate decaying cAMP oscillations following glucose stimulus. To see this we examine the case when Ras is either completely in the Ras ·GTP or Ras ·GDP state. When A is sufficiently large, or sufficiently small, we can reduce the model by making the approximation that the concentration of Ras stays completely either in the GTPbound (r_{ s s }≈1), or GDPbound form (r_{ s s }≈0). This reduces our model, Equations (4a)–(4d), to a twobytwo system, representing change in concentration of Pde (p) and cAMP (x):
where C_{0}=1+C when Ras ·GTP is the dominant Ras form (r_{ s s }≈1) or C_{0}=C when Ras ·GDP dominates (r_{ s s }≈0). In these cases we can find an explicit expression for when the concentration of cAMP exhibits decaying oscillation as it approaches steady state.
Claim 1
Solutions of (5a) and (5b) exhibit an oscillatory approach to steadystate if and only if
where x_{ s s } is the equilibrium solution of (5a, 5b); i.e.,
Condition (6) is derived by computing the Jacobian matrix of Equations (5a) and (5b) at the equilibrium point and determining when the eigenvalues have a nonzero imaginary part. Because of the implicit definition of x_{ s s }, it is rather difficult to apply directly. Therefore we refer to Figures 5a and 5b, in which the overlapping shaded regions show the ranges of D_{0} and D for which the inequality (6) is satisfied, for several values of Γ and for M assumed equal to.01. Figure 5a shows the parameter range for C_{0}=1+C=1.044, the case when r_{ s s }≈1, and Figure 5b shows the parameter range for C_{0}=C=.044, the case when r_{ s s }≈0. Observe that the value of Γ has a large effect on when oscillation occurs: when Γ is small (i.e. the PDEs have high affinity for cAMP), oscillations occur over a large parameter range, but as Γ increases, the parameter range decreases significantly. Finally, because of the factor of M on the righthand side of (6), it is difficult to satisfy the inequality if M is very small; thus decaying oscillations are less likely for the pde1 Δ and pde1^{ala 152} mutants.
General conditions for decaying oscillations
We now seek to understand the parameters that cause decaying cAMP oscillations when Ras activation is free to vary. To do so we numerically examine when the Jacobian matrix determined by the four–by–four system has a complex conjugate pair of eigenvalues with negative real parts. We summarize the results here. We find that as parameters Γ_{1}, N and the ratio $\frac{A}{B}$ increase we are more likely to observe decaying oscillations. We also find that as the parameters that enhance the decay of cAMP are increased, that is as parameters D_{0} and D are increased or as M decreases decaying oscillations are less likely. For a more detailed analysis we refer the reader to the Additional file 1: Supplementary Model.
Experimental validation of predicted cAMP oscillations
Our model makes the novel prediction that for many choices of parameter values, cAMP levels should exhibit decaying oscillations towards a steadystate following glucose stimulus. In order to test this hypothesis we measured the dynamics of the cAMP response in diploid cells of two strains of S. cerevisiae – S288c and Σ 1278b. S288c is considered the standard “reference” genome for yeast studies while Σ 1278b is commonly used for studies of developmental pathways in yeast [53, 54].
For each strain we monitored cAMP levels for eight to twelve minutes following a glucose stimulus. Two typical experimental time series are illustrated in Figure 6. There are cAMP oscillations in both genetic backgrounds but the quantitative features of the oscillations appear to be strain dependent. Σ 1278b exhibits a classical form of decaying oscillations characterized by a large cAMP peak immediately after the stimulus, and dampening oscillations towards a new steady state. In the S288c background, by contrast, the oscillations appear to be delayed, occur with a lower amplitude and do not decay as rapidly. The Σ 1278b strain has a cAMP peak approximately twice as large as that of S288c.
There are various ways in which we could fit the dynamics reported here. We could assume that Ras ·GTP is saturated, as in Ma et al. Cases 1–4 and seek to fit only the two–by–two system, given by Equations (5a) and (5b), to this data. Alternatively we could fit all four equations, (4a)–(4d), to this data. Since this experiment was done in a similar way to the Ma et al. [35] experiment, we seek to fit the data to parameters relating to only Pde feedback, that is parameters D, Γ, and M. The data is taken in increments of fifteen seconds for the first two minutes, thirty seconds for the next five minutes and every minute for the last five minutes. We used a least squares approach to fit the parameters of the mathematical model to the observed data. Because the density of samples is higher in the initial two minutes, the fitting more closely resembles the data over this interval.
We can approximate the dynamics observed for both Σ 1278b and S288c using the model described by Equations (4a)–(4d) (Figure 6A). Table 3 shows the parameter values used to fit the oscillations in Σ 1278b and S288c compared to the parameter values used to fit Ma et al. wild type. For both cases the model requires Pde’s affinity for cAMP, $\frac{1}{\Gamma}$, to be much greater than Pde’s affinity for cAMP in the Ma et al. case. In our fitting we note that the most significant difference between the two strains is in the decay of cAMP with respect to Pde (parameter D).
PDE mutants in the Σ 1278b background
The Σ 1278b wildtype strain exhibited a considerably larger cAMP peak and more pronounced oscillations than S288c (Figure 6), or the the W303 backround used by Ma et al. We therefore undertook additional experiments to examine the behavior of the pde1 Δ and pde2 Δ knockout mutants in Σ 1278b. Typical experimental time series for the PDE mutants in the Σ 1278b background are shown in Figure 7. These experiments qualitatively agree with the dynamics reported by Ma et al. for the W303 background: 1) the pde2 Δ mutant has very similar cAMP dynamics to the wildtype background; and 2) the pde1 Δ mutant reaches a much higher cAMP maximum. One notable difference between the W303 and Σ 1278b pde1 Δ mutants is that cAMP levels in the Σ 1278b background reach a steadystate at or near their maximum, while the W303 pde1 Δ mutant shows a more transient response with cAMP levels decaying towards a steady state significantly lower than the maximum. Because of this observed difference in the cAMP signaling behavior, when fitting model parameters for the pde1 Δ and pde2 Δ mutants we relaxed Condition (b), meaning that we no longer assumed a uniformly slower rate of activation/inactivation for Pde2. Model parameter estimates for the Σ 1278b PDE mutants are given in Table 4. Comparing the parameters of the pde1 Δ mutants in W303 (Table 1) and Σ 1278b (Table 4), we see that the nondimensional parameter M, representing the rate at which phosphodiesterase activity reaches its steady state, increases relative to the wildtype in the Σ 1278b background, reflecting the lack of a transient cAMP response in this strain. The difference between the backgrounds in terms of steady state cAMP levels following stimulus is primarily reflected in the relative change in the parameter D.
Longterm dynamics
If we consider our model over time scales longer than five minutes, our model the cAMPPKA pathway predicts that cAMP levels may either go to a unique steady state or experience sustained oscillations. GarmendiaTorres et al. [42] predicted sustained oscillations in the concentration of cAMP at intermediate stress levels. Since stress affects our model through parameter A (the activation rate of Ras ·GTP), we choose A as our bifurcation parameter. At intermediate stress levels, for example A=.014, the longterm concentration of cAMP experiences sustained oscillations, as seen by the red curve in Figure 8. When stress levels are low, for example A=1.4, the longterm concentration of cAMP remains at a high steadystate value, as seen by the blue line in Figure 8. When stress levels are high, for example A=.005, the longterm concentration of cAMP remains at a low steadystate value, as seen by the black line in Figure 8.
An intuitive biological explanation of this behavior is that at intermediate stress levels Ras (r) alternates between active (GTPbound, r_{ ss }≈1) and inactive (GDPbound, r_{ ss }≈0) states. This oscillation in Ras states causes oscillations in the concentration of cAMP. When stress levels are low, for example in the glucose fed state, Ras is almost completely in the Ras ·GTP state, as in Ma et al. [35] Cases 1–4, forcing the longterm concentration of cAMP to steady state. Similarly, when stress levels are high, Ras is primarily in the Ras ·GDP state, as in Ma et al. [35] Case 5, forcing the longterm concentration of cAMP to a low steady state.
Analysis of sustained oscillations
Analysis of our model suggests that sustained oscillations are the result of PKA feedback through the Ras pathway (variables r and z). In contrast, negative feedback through the PDEs (variable p) stabilizes the system. This analysis was done by examining the stability of the system by evaluating the Jacobian matrix at the equilibrium. We find that our system loses stability through a Hopf bifurcation; thus, indicating that the system experiences sustained oscillations (for details see Additional file 1: Supplementary Model). This instability, leading to sustained oscillations, is a result of Ras varying between its active (GTP) and inactive (GDP) states. Activation/inactivation of Ras is controlled by parameters A, B, and Γ_{ 1 } in our model. We find that as the value of Γ_{ 1 } increases our system is less likely to experience sustained oscillations. Oscillations are also more likely as the ratio $\frac{A}{B}$ increases. In regards to the period and amplitude of oscillations, we observe that as both forms of PKA feedback (through Ras and PDE) slow down, that is as the values of parameters N and M decrease, the length of the period of oscillation increases. Finally, we conclude that as the parameters that enhance the decay of cAMP are decreased, that is as parameters D_{ 0 } and D are decreased and parameter Γ is increased, oscillations are more likely to have longer periods and larger amplitudes.
Discussion
The model presented here is, to the best of our knowledge, the first analytical model of the yeast cAMPPKA pathway that is capable of recapitulating both the shortterm and longterm dynamics of cAMP signaling in wildetype cells. The model we propose is also capable of replicating cAMP dynamics that have been observed single and double phosphodiesterase mutants. In addition to capturing behaviors previously described our model makes new predictions about cAMP oscillations. Our work also highlights variation in cAMP signaling between genetic backgrounds.
cAMP oscillations
Analysis of a nondimensional system derived from our model suggests that a substantial fraction of the model parameter space for the wildtype signaling network should result in decaying cAMP oscillations following glucose stimulus. While long timescale oscillations have previously been predicted [39, 42], shortterm oscillations, on the scale of minutes, have not been noted in the yeast literature. To test our predictions we experimentally measured the cAMP response in two different diploid strains of S. cerevisiae and observed oscillations as predicted by the model. While ours is the first study to document short timescale cAMP oscillations in yeast, such oscillations have been noted in other systems including mammals, amphibians, and slime molds [55–57]. For example, Dyachock et al. [55] reported short timescale glucose induced oscillations in mouse pancreatic βcells. They showed that these oscillations correlated with pulsatile insulin release. Similar oscillations are observed in response to stimulation of such cells by the peptide hormone GLP1 [58]. The slime mold Dictyostelium also exhibits short timescale cAMP oscillations [57]. Similar to the model we present here, Maeda et al. [57] suggested that positive feedback by PKA on phosphodiesterase activity and PKA negative feedback on adenylate cyclase activity are critical features of the molecular circuit that generates these oscillations in Dictyostelium. Whether cAMP oscillations in these diverse eukaryotic groups are due to conserved interactions in the cAMP pathway or whether such dynamics represent a convergent feature of cAMP signaling is a key question for future comparative studies.
Our model is also able to generate longertime scale cAMP oscillations. Such oscillatory dynamics should be reflected in PKA activity. Recent studies [59, 60] have demonstrated that PKA activation plays an important role in regulating oscillatory patterns of nuclear localization of stressresponsive transcription factors such as Msn2. Different dynamical patterns of Msn2 regulation in turn can lead to qualitatively different expression outputs of downstream targets [59].
Comparison to other models
There have been a number of recent studies that have proposed mathematical models for the cAMPPKA network in yeast, including Gonze et al. [39], Cazzaniga et al. [40], and Williamson et al. [41]. A limitation shared by all three of these models is that they can not recapitulate the dynamical behaviors of the double phosphodiesterase mutant (pde1 Δ pde2 Δ). This can be seen by setting the variables representing the phosphodiesterases to zero in any of these models (i.e. representing the double knockout); the result is a rapid accumulation of cAMP. This contrary to the observed data reported by Ma et al. [35]. Our model uses competitive inhibition between Pde1, Pde2 and the GTPase activating proteins Ira1 and/or Ira2 to explain the double PDE knockout case; thus, competitive inhibition between Pde1, Pde2, and Ira1/2 is what sets our model apart from these previous models.
Gonze et al. [39] propose a stochastic model, that is a direct extension of the GarmendiaTorres et al. model [42], to explore the effects of the cAMPPKA pathway on the transcription factor Msn2. This model only focuses on the long term dynamics and can not be used to understand any of the PDE mutant cases proposed by Ma et al. Cazzaniga and colleagues [40] employed a stochastic modeling approach to study the yeast cAMPPKA network. Though they do explore the effects of different PDE activities, their model does not explicitly account for the mutant cases described by Ma et al. and again they can not replicate the double phosphodiesterase knockout. Also they do not consider the long term dynamics proposed by GarmendiaTorres et al. [42]. The deterministic model presented by Williamson et al. [41] is the closest in approach to ours. Their model is able to recapitulate the relevant dynamics of the Pde1 and Pde2 single mutants but does not account for the behavior of the PDE double mutant. As we mentioned above this is significant since the double knockout case is the most surprising case and the hardest case to explain. They also only examined the short term dynamics of cAMP, where as our model accounts for both the long and short term dynamics of cAMP.
Conclusions
The mathematical model we present here replicates the dynamical behavior of wildtype and mutant cells and leads us to novel predictions of oscillatory cAMP dynamics in S. cerevisiae. Experiments motivated by this modeling effort confirmed these predictions and also highlighted variation in oscillatory phenomena among yeast isolates. The mathematical model, in turn, is capable of reproducing the different signaling phenotypes we observed and allows us to predict which parts of the network are most likely to contribute to signaling variation. Because it has both predictive and explanatory power this model can serve as a foundation for future mathematical and experimental explorations of cAMP signaling. A future challenge for both experimental and modeling studies of cAMP signaling will be to determine whether the grossly similar oscillatory dynamics observed across eukaryotes is due to a common mechanistic basis or whether this represents an example of convergent evolution.
Methods
Strains and cAMP assay
Diploid strains of the laboratory backgrounds Σ 1278b (G85 [‘wildtype’, ura3 Δ 0, his3 Δ 0::hisG], G30025 [‘pde1 Δ’, pde1::bcKanMX4, ura3 Δ 0, his3 Δ 0::hisG], G30025 [‘pde2 Δ’, pde2::bcKanMX4, ura3 Δ 0, his3 Δ 0::hisG], gifts of T. Galitski [61]) and S288c (BY4743, gift of D. Lew) were used to study cAMP dynamics. Cells were prepared following a protocol adapted from Paiardi et al. [17]. Approximately 2×10^{9} cells/ml were grown in rich media (1% yeast extract, 2% peptone and 2% glucose) at 30°C and transferred to SC media with 0.1% glucose and 3% gycerol for overnight incubation at 30°C. Cells were collected, washed and incubated in 25 mM MES buffer, pH 6 (Boston Bio Products, Worcester, MA) for 30 minutes at 30°C before glucose induction. Glucose was added to 2 ml of cell culture(approx. 2 ×10^{8} cells/ml) to a final concentration of 100 mM and 25 mM. After glucose induction cells were collected over a 12minute time course in 15 second increments for the first minute and at 30 second increments for the rest of the time course. 3 ×10^{6} cells were collected at each time point, fixed in 300 μ l nbutanolsaturated 1M formic acid (Jiang et.al., 1998) and frozen at 80°C. Cells were lysed using four freezethaw cycles and freezedried in a Speed Vac concentrator (Savant Instruments, Farmingdale, NY). Cell extracts were used to determine cAMP concentration using a cAMP Biotrak Enzyme immunoassay kit (Amersham, GE Healthcare) and a SpectraMax spectrophotometer (Molecular Devices, Sunnyvale, CA).
Model fitting
For the purposes of model fitting, where available we used estimates of key biochemical parameters as reported in [42]. In cases where such parameters were not avaiable, we used a method of least squares, implemented in the subroutine lmdif.f in the Fortran MINPACK library [62]. This subroutine minimizes the sum of the squares of m nonlinear functions in n variables by a modification of the LevenbergMarquardt algorithm. We wrote a C computer program that calculates the appropriate functions that are used by this algorithm for the model described above. Additional details of the model fitting procedure are provided in the Additional file 1: Supplemental materials.
References
 1.
Fimia GM, SassoneCorsi P: Cyclic AMP signalling. J Cell Sci. 2001, 114 (Pt 11): 19711972.
 2.
Beavo JA, Brunton LL: Cyclic nucleotide research – still expanding after half a century. Nat Rev Mol Cell Biol. 2002, 3 (9): 710718. 10.1038/nrm911.
 3.
Holz GG, Kang G, Harbeck M, Roe MW, Chepurny OG: Cell physiology of cAMP sensor Epac. J Physiol. 2006, 577 (Pt 1): 515. [http://dx.doi.org/10.1113/jphysiol.2006.119644]
 4.
Rall TW, Sutherland EW: Formation of a cyclic adenine ribonucleotide by tissue particles. J Biol Chem. 1958, 232 (2): 10651076.
 5.
Carlone DL, Richards JS: Functional interactions, phosphorylation, and levels of 3’,5’cyclic adenosine monophosphateregulatory element binding protein and steroidogenic factor1 mediate hormoneregulated and constitutive expression of aromatase in gonadal cells. Mol Endocrinol. 1997, 11 (3): 292304. 10.1210/me.11.3.292.
 6.
Buck LB: Information coding in the vertebrate olfactory system. Annu Rev Neurosci. 1996, 19: 517544. 10.1146/annurev.ne.19.030196.002505. [http://dx.doi.org/10.1146/annurev.ne.19.030196.002505] 10.1146/annurev.ne.19.030196.002505
 7.
Walsh DA, Perkins JP, Krebs EG: An adenosine 3’,5’monophosphatedependant protein kinase from rabbit skeletal muscle. J Biol Chem. 1968, 243 (13): 37633765.
 8.
Hong S, Bang S, Hyun S, Kang J, Jeong K, Paik D, Chung J, Kim J: cAMP signalling in mushroom bodies modulates temperature preference behaviour in Drosophila. Nature. 2008, 454 (7205): 771775. [http://dx.doi.org/10.1038/nature07090]
 9.
Haynes RC, Sutherland EW, Rall TW: The role of cyclic adenylic acid in hormone action. Recent Prog Horm Res. 1960, 16: 121138.
 10.
Thevelein JM, de Winde JH: Novel sensing mechanisms and targets for the cAMPprotein kinase A pathway in the yeast Saccharomyces cerevisiae. Mol Microbiol. 1999, 33 (5): 904918. 10.1046/j.13652958.1999.01538.x.
 11.
Toda T, Uno I, Ishikawa T, Powers S, Kataoka T, Broek D, Cameron S, Broach J, Matsumoto K, Wigler M: In yeast, RAS proteins are controlling elements of adenylate cyclase. Cell. 1985, 40: 2736. 10.1016/00928674(85)903058.
 12.
Kraakman L, Lemaire K, Ma P, Teunissen AW, Donaton MC, Dijck PV, Winderickx J, de Winde JH, Thevelein JM: A Saccharomyces cerevisiae Gprotein coupled receptor, Gpr1, is specifically required for glucose activation of the cAMP pathway during the transition to growth on glucose. Mol Microbiol. 1999, 32 (5): 10021012. 10.1046/j.13652958.1999.01413.x.
 13.
van Aelst, Jans AW, Thevelein JM: Involvement of the CDC25 gene product in the signal transmission pathway of the glucoseinduced RASmediated cAMP signal in the yeast Saccharomyces cerevisiae. J Gen Microbiol. 1991, 137 (2): 341349. 10.1099/002212871372341.
 14.
Broek D, Toda T, Michaeli T, Levin L, Birchmeier C, Zoller M, Powers S, Wigler M: The S. cerevisiae CDC25 gene product regulates the RAS/adenylate cyclase pathway. Cell. 1987, 48 (5): 789799. 10.1016/00928674(87)900766.
 15.
Créchet JB, Poullet P, Mistou MY, Parmeggiani A, Camonis J, BoyMarcotte E, Damak F, Jacquet M: Enhancement of the GDPGTP exchange of RAS proteins by the carboxylterminal domain of SCD25. Science. 1990, 248 (4957): 866868. 10.1126/science.2188363.
 16.
BoyMarcotte E, Ikonomi P, Jacquet M: SDC25, a dispensable Ras guanine nucleotide exchange factor of Saccharomyces cerevisiae differs from CDC25 by its regulation. Mol Biol Cell. 1996, 7 (4): 529539.
 17.
Paiardi C, Belotti F, Colombo S, Tisi R, Martegani E: The large Nterminal domain of Cdc25 protein of the yeast Saccharomyces cerevisiae is required for glucoseinduced Ras2 activation. FEMS Yeast Res. 2007, 7 (8): 12701275. 10.1111/j.15671364.2007.00300.x. [http://dx.doi.org/10.1111/j.15671364.2007.00300.x] 10.1111/j.15671364.2007.00300.x
 18.
Colombo S, Ma P, Cauwenberg L, Winderickx J, Crauwels M, Teunissen A, Nauwelaers D, de Winde JH, Gorwa MF, Colavizza D, Thevelein JM: Involvement of distinct Gproteins, Gpa2 and Ras, in glucose and intracellular acidificationinduced cAMP signalling in the yeast Saccharomyces cerevisiae. EMBO J. 1998, 17 (12): 33263341. 10.1093/emboj/17.12.3326. [http://dx.doi.org/10.1093/emboj/17.12.3326] 10.1093/emboj/17.12.3326
 19.
Colombo S, Ronchetti D, Thevelein JM, Winderickx J, Martegani E: Activation state of the Ras2 protein and glucoseinduced signaling in Saccharomyces cerevisiae. J Biol Chem. 2004, 279 (45): 4671546722. 10.1074/jbc.M405136200. [http://dx.doi.org/10.1074/jbc.M405136200] 10.1074/jbc.M405136200
 20.
Thevelein JM, Bonini B, Castermans D, Haesendonckx S, Kriel J, Louwet W, Thayumanavan P, Popova Y, RubioTexeira M, Schepers W, Vandormael P, Zeebroeck GV, Verhaert P, Versele M, Voordeckers K: Novel mechanisms in nutrient activation of the yeast protein kinase A pathway. Acta Microbiol Immunol Hung. 2008, 55 (2): 7589. 10.1556/AMicr.55.2008.2.1. [http://dx.doi.org/10.1556/AMicr.55.2008.2.1] 10.1556/AMicr.55.2008.2.1
 21.
Tanaka K, Matsumoto K, TohE A: IRA1, an inhibitory regulator of the RAScyclic AMP pathway in Saccharomyces cerevisiae. Mol Cell Biol. 1989, 9 (2): 757768.
 22.
Tanaka K, Nakafuku M, Tamanoi F, Kaziro Y, Matsumoto K, Tohe A: IRA2, a second gene of Saccharomyces cerevisiae that encodes a protein with a domain homologous to mammalian Ras GTPaseactivating protein. Mol Cell Biol. 1990, 10 (8): 43034313.
 23.
Rolland F, Winde JHD, Lemaire K, Boles E, Thevelein JM, Winderickx J: Glucoseinduced cAMP signalling in yeast requires both a Gprotein coupled receptor system for extracellular glucose detection and a separable hexose kinasedependent sensing process. Mol Microbiol. 2000, 38 (2): 348358. 10.1046/j.13652958.2000.02125.x.
 24.
Santangelo GM: Glucose signaling in Saccharomyces cerevisiae. Microbiol Mol Biol Rev. 2006, 70: 253282. 10.1128/MMBR.70.1.253282.2006. [http://dx.doi.org/10.1128/MMBR.70.1.253282.2006] 10.1128/MMBR.70.1.253282.2006
 25.
Versele M, de Winde JH, Thevelein JM: A novel regulator of G protein signalling in yeast, Rgs2, downregulates glucoseactivation of the cAMP pathway through direct inhibition of Gpa2. EMBO J. 1999, 18 (20): 55775591. 10.1093/emboj/18.20.5577. [http://dx.doi.org/10.1093/emboj/18.20.5577] 10.1093/emboj/18.20.5577
 26.
Kehrl J, Sinnarajah S: RGS2: a multifunctional regulator of Gprotein signaling. Int J Biochem Cell Biol. 2002, 34 (5): 432438. 10.1016/S13572725(01)001418.
 27.
Toda T, Cameron S, Sass P, Zoller M, Scott JD, McMullen B, Hurwitz M, Krebs EG, Wigler M: Cloning and characterization of BCY1, a locus encoding a regulatory subunit of the cyclic AMPdependent protein kinase in Saccharomyces cerevisiae. Mol Cell Biol. 1987, 7 (4): 13711377.
 28.
Toda T, Cameron S, Sass P, Zoller M, Wigler M: Three different genes in S.cerevisiae encode the catalytic subunits of the cAMPdependent protein kinase. Cell. 1987, 50 (2): 277287. 10.1016/00928674(87)902236.
 29.
Gill GN, Garren LD: Role of the receptor in the mechanism of action of adenosine 3’:5’cyclic monophosphate. Proc Natl Acad Sci U S A. 1971, 68 (4): 786790. 10.1073/pnas.68.4.786.
 30.
Proud CG, Rylatt DB, Yeaman SJ, Cohen P: Amino acid sequences at the two sites on glycogen synthetase phosphorylated by cyclic AMPdependent protein kinase and their dephosphorylation by protein phosphataseIII. FEBS Lett. 1977, 80 (2): 435442. 10.1016/00145793(77)804936.
 31.
Fang X, Yu SX, Lu Y, Bast RCJr, Woodgett JR, Mills GB: Phosphorylation and inactivation of glycogen synthase kinase 3 by protein kinase A. Proc Natl Acad Sci U S A. 2000, 97 (22): 119601165. 10.1073/pnas.220413597.
 32.
Taylor SS, Buechler JA, Yonemoto W: cAMPdependent protein kinase: framework for a diverse family of regulatory enzymes. Annu Rev Biochem. 1990, 59: 9711005. 10.1146/annurev.bi.59.070190.004543. [http://dx.doi.org/10.1146/annurev.bi.59.070190.004543] 10.1146/annurev.bi.59.070190.004543
 33.
Francis SH, Corbin JD: Structure and function of cyclic nucleotidedependent protein kinases. Annu Rev Physiol. 1994, 56: 237272. 10.1146/annurev.ph.56.030194.001321. [http://dx.doi.org/10.1146/annurev.ph.56.030194.001321] 10.1146/annurev.ph.56.030194.001321
 34.
Mitsuzawa H: Responsiveness to exogenous cAMP of a Saccharomyces cerevisiae strain conferred by naturally occurring alleles of PDE1 and PDE2. Genetics. 1993, 135 (2): 321326.
 35.
Ma P, Wera S, Dijck PV, Thevelein JM: The PDE1encoded lowaffinity phosphodiesterase in the yeast Saccharomyces cerevisiae has a specific function in controlling agonistinduced cAMP signaling. Mol Biol Cell. 1999, 10: 91104.
 36.
Conti M, Beavo J: Biochemistry and physiology of cyclic nucleotide phosphodiesterases: essential components in cyclic nucleotide signaling. Annu Rev Biochem. 2007, 76: 481511. 10.1146/annurev.biochem.76.060305.150444. [http://dx.doi.org/10.1146/annurev.biochem.76.060305.150444] 10.1146/annurev.biochem.76.060305.150444
 37.
Park J, Grant CM, Dawes IW: The highaffinity cAMP phosphodiesterase of Saccharomyces cerevisiae is the major determinant of cAMP levels in stationary phase: involvement of different branches of the Rascyclic AMP pathway in stress responses. Biochem Biophys Res Commun. 2005, 327: 311319. 10.1016/j.bbrc.2004.12.019. [http://dx.doi.org/10.1016/j.bbrc.2004.12.019] 10.1016/j.bbrc.2004.12.019
 38.
Wilson D, Fiori A, Brucker KD, Dijck PV, Stateva L: Candida albicans Pde1p and Gpa2p comprise a regulatory module mediating agonistinduced cAMP signalling and environmental adaptation. Fungal Genet Biol. 2010, 47 (9): 742752. 10.1016/j.fgb.2010.06.006.
 39.
Gonze D, Goldbeter MJA: Stochastic modelling of nucleocytoplasmic oscillations of the transcription factor Msn2 in yeast. J R Soc Interface. 2008, 5 (Suppl 1): S95S109. [http://dx.doi.org/10.1098/rsif.2008.0141.focus]
 40.
Cazzaniga P, Pescini D, Besozzi D, Mauri G, Colombo S, Martegani E: Modeling and stochastic simulation of the Ras/cAMP/PKA pathway in the yeast Saccharomyces cerevisiae evidences a key regulatory function for intracellular guanine nucleotides pools. J Biotechnol. 2008, 133 (3): 377385. 10.1016/j.jbiotec.2007.09.019. [http://dx.doi.org/10.1016/j.jbiotec.2007.09.019] 10.1016/j.jbiotec.2007.09.019
 41.
Williamson T, Schwartz J, Kell DB, Stateva L: Deterministic mathematical models of the cAMP pathway in Saccharomyces cerevisiae. BMC Syst Biol. 2009, 3: 7010.1186/17520509370. [http://dx.doi.org/10.1186/17520509370] 10.1186/17520509370
 42.
GarmendiaTorres C, Goldbeter A, Jacquet M: Nucleocytoplasmic oscillations of the yeast transcription factor Msn2: evidence for periodic PKA activation. Current Biol. 2007, 17: 10441049. 10.1016/j.cub.2007.05.032.
 43.
Görner W, Durchschlag E, MartinezPastor MT, Estruch F, Ammerer G, Hamilton B, Ruis H, Schüller C: Nuclear localization of the C2H2 zinc finger protein Msn2p is regulated by stress and protein kinase A activity. Genes Dev. 1998, 12 (4): 586597. 10.1101/gad.12.4.586. [http://genesdev.cshlp.org/content/12/4/586.long] 10.1101/gad.12.4.586
 44.
Hu Y, Liu E, Bai X, Zhang A: The localization and concentration of the PDE2encoded highaffinity cAMP phosphodiesterase is regulated by cAMPdependent protein kinase A in the yeast Saccharomyces cerevisiae. FEMS Yeast Res. 2010, 10 (2): 177187. 10.1111/j.15671364.2009.00598.x.
 45.
Lacal JC, McCormick F: The Ras Superfamily of GTPases. 1993, CRC Press
 46.
Pardo LA, Lazo PS, Ramos S: Activation of adenylate cyclase in Cdc25 mutants of Saccharomyces cerevisiae. FEBS Lett. 1993, 319 (3): 237243. 10.1016/00145793(93)805548.
 47.
Garreau H, Geymonat M, Renault G, Jacquet M: Membraneanchoring domains of Cdc25p, a Saccharomyces cerevisiae Ras exchange factor. Biol Cell. 1996, 86 (23): 93102. 10.1016/02484900(96)84771X.
 48.
Mintzer KA, Field J: The SH3 domain of the S. cerevisiae Cdc25p binds adenylyl cyclase and facilitates Ras regulation of cAMP signalling. Cell Signal. 1999, 11 (2): 127135. 10.1016/S08986568(98)000448.
 49.
Gross E, Goldberg D, Levitzki A: Phosphorylation of the S. cerevisiae Cdc25 in response to glucose results in its dissociation from Ras. Nature. 1992, 360 (6406): 762765. 10.1038/360762a0. [http://dx.doi.org/10.1038/360762a0] 10.1038/360762a0
 50.
Gross A, Winograd S, Marbach I, Levitzki A: The Nterminal half of Cdc25 is essential for processing glucose signaling in Saccharomyces cerevisiae. Biochemistry. 1999, 38 (40): 1325213262. 10.1021/bi9909849.
 51.
Segel I H: Enzyme Kinetics: Behavior and Analysis of Rapid Equilibrium and SteadyState Enzyme Systems. 1993, New York: John Wiley and Sons Inc.
 52.
Belotti F, Tisi R, Paiardi C, Rigamonti M, Groppi S, Martegani E: Localization of Ras signaling complex in budding yeast. Biochim Biophys Acta. 2012, 1823 (7): 12081216. 10.1016/j.bbamcr.2012.04.016.
 53.
Gimeno CJ, Ljungdahl PO, Styles CA, Fink GR: Unipolar cell divisions in the yeast S. cerevisiae lead to filamentous growth: regulation by starvation and RAS. Cell. 1992, 68 (6): 10771090. 10.1016/00928674(92)90079R.
 54.
Lorenz MC, Heitman J: Yeast pseudohyphal growth is regulated by GPA2, a G protein alpha homolog. EMBO J. 1997, 16 (23): 70087018. 10.1093/emboj/16.23.7008. [http://dx.doi.org/10.1093/emboj/16.23.7008] 10.1093/emboj/16.23.7008
 55.
Dyachok O, IdevallHagren O, Sagetorp J, Tian G, Wuttke A, Arrieumerlou C, Akusjarvi G, Gylfe E, Tengholm A: Glucoseinduced cyclic AMP oscillations regulate pulsatile insulin secretion. Cell Metab. 2008, 8: 2637. 10.1016/j.cmet.2008.06.003. [http://dx.doi.org/10.1016/j.cmet.2008.06.003] 10.1016/j.cmet.2008.06.003
 56.
Gorbunova Y, Spitzer NC: Dynamic interactions of cyclic AMP transients and spontaneous Ca(2+) spikes. Nature. 2002, 418 (6893): 9396. 10.1038/nature00835. [http://dx.doi.org/10.1038/nature00835] 10.1038/nature00835
 57.
Maeda M, Lu S, Shaulsky G, Miyazaki Y, Kuwayama H, Tanaka Y, Kuspa A, Loomis WF: Periodic signaling controlled by an oscillatory circuit that includes protein kinases ERK2 and PKA. Science. 2004, 304 (5672): 875878. 10.1126/science.1094647. [http://dx.doi.org/10.1126/science.1094647] 10.1126/science.1094647
 58.
Dyachok O, Sagetorp J, Isakov Y, Tengholm A: cAMP oscillations restrict protein kinase A redistribution in insulinsecreting cells. Biochem Soc Trans. 2006, 34 (Pt 4): 498501. [http://dx.doi.org/10.1042/BST0340498]
 59.
Hao N, O’Shea EK: Signaldependent dynamics of transcription factor translocation controls gene expression. Nat Struct Mol Biol. 2012, 19: 3139.
 60.
Cai L, Dalal CK, Elowitz MB: Frequencymodulated nuclear localization bursts coordinate gene regulation. Nature. 2008, 455 (7212): 485490. 10.1038/nature07292.
 61.
Drees BL, Thorsson V, Carter GW, Rives AW, Raymond MZ, AvilaCampillo I, Shannon P, Galitski T: Derivation of genetic interaction networks from quantitative phenotype data. Genome Biol. 2005, 6 (4): R3810.1186/gb200564r38. [http://dx.doi.org/10.1186/gb200564r38] 10.1186/gb200564r38
 62.
Moré JJ, Sorensen DC, Hillstrom KE, Garbow BS: The MINPACK Project. Sources and Development of Mathematical Software. Edited by: Cowell WJ, Cowell WJ. 1984, Englewood Cliffs: PrenticeHall
 63.
Jorgensen P, Nishikawa JL, Breitkreutz BJ, Tyers M: Systematic identification of pathways that couple cell growth and division in yeast. Science. 2002, 297 (5580): 395400. 10.1126/science.1070850.
Acknowledgements
We acknowledge Joseph Heitman, David McCandlish, and two anonymous reviewers for helpful comments on the manuscript. Timothy Galitski, Institute for Systems Biology provided the Σ 1278b strains used in this study. This work was supported by an NIH award to PMM and the Duke Center for Systems Biology (P50GM08188301).
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
KG carried out the mathematical modeling. OK conducted biochemical measurements of cAMP signaling. DGS contributed to the mathematical modeling. PMM contributed to design of the study and carried out analysis of the experimental data. All authors participated in the drafting of the manuscript, and have read and approved the final manuscript.
Electronic supplementary material
Additional file 1:Supplementary model[35, 42, 62, 63].(PDF 86 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Gonzales, K., Kayıkçı, Ö., Schaeffer, D.G. et al. Modeling mutant phenotypes and oscillatory dynamics in the Saccharomyces cerevisiae cAMPPKA pathway. BMC Syst Biol 7, 40 (2013). https://doi.org/10.1186/17520509740
Received:
Accepted:
Published:
Keywords
 Yeast
 Signal transduction
 Second messenger
 Genetic variation