Membrane localization of scaffold proteins promotes graded signaling in the yeast MAP kinase cascade

Simulating BRAFV600E-MEK-ERK signalling dynamics in response to vertical inhibition treatment strategies

In vertical inhibition treatment strategies, multiple components of an intracellular pathway are simultaneously inhibited. Vertical inhibition of the BRAFV600E-MEK-ERK signalling pathway is a standard of care for treating BRAFV600E-mutated melanoma where two targeted cancer drugs, a BRAFV600E-inhibitor, and a MEK inhibitor, are administered in combination. Targeted therapies have been linked to early onsets of drug resistance, and thus treatment strategies of higher complexities and lower doses have been proposed as alternatives to current clinical strategies. However, finding optimal complex, low-dose treatment strategies is a challenge, as it is possible to design more treatment strategies than are feasibly testable in experimental settings. To quantitatively address this challenge, we develop a mathematical model of BRAFV600E-MEK-ERK signalling dynamics in response to combinations of the BRAFV600E-inhibitor dabrafenib (DBF), the MEK inhibitor trametinib (TMT), and the ERK-inhibitor SCH772984 (SCH). From a model of the BRAFV600E-MEK-ERK pathway, and a set of molecular-level drug-protein interactions, we extract a system of chemical reactions that is parameterised by in vitro data and converted to a system of ordinary differential equations (ODEs) using the law of mass action. The ODEs are solved numerically to produce simulations of how pathway-component concentrations change over time in response to different treatment strategies, i.e., inhibitor combinations and doses. The model can thus be used to limit the search space for effective treatment strategies that target the BRAFV600E-MEK-ERK pathway and warrant further experimental investigation. The results demonstrate that DBF and DBF-TMT-SCH therapies show marked sensitivity to BRAFV600E concentrations in silico, whilst TMT and SCH monotherapies do not.

A guide to ERK dynamics, part 1: mechanisms and models

Extracellular signal-regulated kinase (ERK) has long been studied as a key driver of both essential cellular processes and disease. A persistent question has been how this single pathway is able to direct multiple cell behaviors, including growth, proliferation, and death. Modern biosensor studies have revealed that the temporal pattern of ERK activity is highly variable and heterogeneous, and critically, that these dynamic differences modulate cell fate. This two-part review discusses the current understanding of dynamic activity in the ERK pathway, how it regulates cellular decisions, and how these cell fates lead to tissue regulation and pathology. In part 1, we cover the optogenetic and live-cell imaging technologies that first revealed the dynamic nature of ERK, as well as current challenges in biosensor data analysis. We also discuss advances in mathematical models for the mechanisms of ERK dynamics, including receptor-level regulation, negative feedback, cooperativity, and paracrine signaling. While hurdles still remain, it is clear that higher temporal and spatial resolution provide mechanistic insights into pathway circuitry. Exciting new algorithms and advanced computational tools enable quantitative measurements of single-cell ERK activation, which in turn inform better models of pathway behavior. However, the fact that current models still cannot fully recapitulate the diversity of ERK responses calls for a deeper understanding of network structure and signal transduction in general.

A gene cluster with positive and negative elements controls bistability and hysteresis of the Crippled versus Normal growth in the fungus Podospora anserina

The Crippled Growth (CG) cell degeneration of the model ascomycete Podospora anserina (strain S) is controlled by a prion-like element and has been linked to the self-activation of the PaMpk1 MAP kinase cascade. Here, we report on the identification of the "86-11" locus containing twelve genes, ten of which are involved either in setting up the self-activation loop of CG or in inhibiting this loop, as demonstrated by targeted gene deletion. Interestingly, deletion of the whole locus results only in the elimination of CG and in no detectable additional physiological defect. Sequence comparison shows that these ten genes belong to four different families, each one endowed with a specific activity: two encode factors activating the loop, a third one encodes a factor crucial for inhibition of the loop and the fourth one participates in inhibiting the loop in a pathway parallel to the one controlled by the previously described PDC1 gene. Intriguingly, a very distant homologue of this "86-11" locus is present at the syntenic position in Podospora comata (strain T) that do not present Crippled Growth. Introgression of the P. comata strain T locus in P. anserina strain S and the P. anserina strain S in P. comata strain T showed that both drive CG in the P. anserina strain S genetic background, but not in the genetic background of strain P. comata T, indicating that genetic determinants outside the twelve-gene locus are responsible for lack of CG in P. comata strain T. Our data question the role of this twelve-gene locus in the physiology of P. anserina.

Optimal ligand discrimination by asymmetric dimerization and turnover of interferon receptors

In multicellular organisms, antiviral defense mechanisms evoke a reliable collective immune response despite the noisy nature of biochemical communication between tissue cells. A molecular hub of this response, the interferon I receptor (IFNAR), discriminates between ligand types by their affinity regardless of concentration. To understand how ligand type can be decoded robustly by a single receptor, we frame ligand discrimination as an information-theoretic problem and systematically compare the major classes of receptor architectures: allosteric, homodimerizing, and heterodimerizing. We demonstrate that asymmetric heterodimers achieve the best discrimination power over the entire physiological range of local ligand concentrations. This design enables sensing of ligand presence and type, and it buffers against moderate concentration fluctuations. In addition, receptor turnover, which drives the receptor system out of thermodynamic equilibrium, allows alignment of activation points for ligands of different affinities and thereby makes ligand discrimination practically independent of concentration. IFNAR exhibits this optimal architecture, and our findings thus suggest that this specialized receptor can robustly decode digital messages carried by its different ligands.

Activation of Cph1 causes ß(1,3)-glucan unmasking in Candida albicans and attenuates virulence in mice in a neutrophil-dependent manner

Masking the immunogenic cell wall epitope ß(1,3)-glucan under an outer layer of mannosylated glycoproteins is an important virulence factor deployed by Candida albicans during infection. Consequently, increased ß(1,3)-glucan exposure (unmasking) reveals C. albicans to the host's immune system and attenuates its virulence. We have previously shown that activation of the Cek1 MAPK pathway via expression of a hyperactive allele of an upstream kinase (STE11ΔN467) induced unmasking. It also increased survival of mice in a murine disseminated candidiasis model and attenuated kidney fungal burden by ≥33 fold. In this communication, we utilized cyclophosphamide-induced immunosuppression to test if the clearance of the unmasked STE11ΔN467 mutant was dependent on the host immune system. Suppression of the immune response by cyclophosphamide reduced the attenuation in fungal burden caused by the STE11ΔN467 allele. Moreover, specific depletion of neutrophils via 1A8 antibody treatment also reduced STE11ΔN467-dependent fungal burden attenuation, but to a lesser extent than cyclophosphamide, demonstrating an important role for neutrophils in mediating fungal clearance of unmasked STE11ΔN467 cells. In an effort to understand the mechanism by which Ste11ΔN467 causes unmasking, transcriptomics were used to reveal that several components in the Cek1 MAPK pathway were upregulated, including the transcription factor CPH1 and the cell wall sensor DFI1. In this report we show that a cph1ΔΔ mutation restored ß(1,3)-glucan exposure to wild-type levels in the STE11ΔN467 strain, confirming that Cph1 is the transcription factor mediating Ste11ΔN467-induced unmasking. Furthermore, Cph1 is shown to induce a positive feedback loop that increases Cek1 activation. In addition, full unmasking by STE11ΔN467 is dependent on the upstream cell wall sensor DFI1. However, while deletion of DFI1 significantly reduced Ste11ΔN467-induced unmasking, it did not impact activation of the downstream kinase Cek1. Thus, it appears that once stimulated by Ste11ΔN467, Dfi1 activates a parallel signaling pathway that is involved in Ste11ΔN467-induced unmasking.

Characterization of the Intramolecular Interactions and Regulatory Mechanisms of the Scaffold Protein Tks4

The scaffold protein Tks4 is a member of the p47^phox-related organizer superfamily. It plays a key role in cell motility by being essential for the formation of podosomes and invadopodia. In addition, Tks4 is involved in the epidermal growth factor (EGF) signaling pathway, in which EGF induces the translocation of Tks4 from the cytoplasm to the plasma membrane. The evolutionarily-related protein p47^phox and Tks4 share many similarities in their N-terminal region: a phosphoinositide-binding PX domain is followed by two SH3 domains (so called “tandem SH3”) and a proline-rich region (PRR). In p47^phox, the PRR is followed by a relatively short, disordered C-terminal tail region containing multiple phosphorylation sites. These play a key role in the regulation of the protein. In Tks4, the PRR is followed by a third and a fourth SH3 domain connected by a long (~420 residues) unstructured region. In p47^phox, the tandem SH3 domain binds the PRR while the first SH3 domain interacts with the PX domain, thereby preventing its binding to the membrane. Based on the conserved structural features of p47^phox and Tks4 and the fact that an intramolecular interaction between the third SH3 and the PX domains of Tks4 has already been reported, we hypothesized that Tks4 is similarly regulated by autoinhibition. In this study, we showed, via fluorescence-based titrations, MST, ITC, and SAXS measurements, that the tandem SH3 domain of Tks4 binds the PRR and that the PX domain interacts with the third SH3 domain. We also investigated a phosphomimicking Thr-to-Glu point mutation in the PRR as a possible regulator of intramolecular interactions. Phosphatidylinositol-3-phosphate (PtdIns(3)P) was identified as the main binding partner of the PX domain via lipid-binding assays. In truncated Tks4 fragments, the presence of the tandem SH3, together with the PRR, reduced PtdIns(3)P binding, while the presence of the third SH3 domain led to complete inhibition.

Signal-mediated localization of Candida albicans pheromone response pathway components

 

A MAPK cascade consists of three kinases, (MEKK, MEK and MAPK), that are sequentially activated in response to a stimulus and serve to transmit signals. In C. albicans and in yeast, an MAPK cascade is linked to the pheromone pathway through a scaffold protein (Cst5 and Ste5, respectively). Cst5 is much shorter and lacks key domains compared to Ste5, so in C. albicans, other elements, in particular the MEKK Ste11, play key roles in controlling the associations and localizations of network components.

Abstract

Candida albicans opaque cells release pheromones to stimulate cells of opposite mating type to activate their pheromone response pathway. Although this fungal pathogen shares orthologous proteins involved in the process with Saccharomyces cerevisiae, the pathway in each organism has unique characteristics. We have used GFP-tagged fusion proteins to investigate the localization of the scaffold protein Cst5, as well as the MAP kinases Cek1 and Cek2, during pheromone response in C. albicans. In wild-type cells, pheromone treatment directed Cst5-GFP to surface puncta concentrated at the tips of mating projections. These puncta failed to form in cells defective in either the Gα or β subunits. However, they still formed in response to pheromone in cells missing Ste11, but with the puncta distributed around the cell periphery in the absence of mating projections. These puncta were absent from hst7Δ/Δ cells, but could be detected in the ste11Δ/Δ hst7Δ/Δ double mutant. Cek2-GFP showed a strong nuclear localization late in the response, consistent with a role in adaptation, while Cek1-GFP showed a weaker, but early increase in nuclear localization after pheromone treatment. Activation loop phosphorylation of both Cek1 and Cek2 required the presence of Ste11. In contrast to Cek2-GFP, which showed no localization signal in ste11Δ/Δ cells, Cek1-GFP showed enhanced nuclear localization that was pheromone independent in the ste11Δ/Δ mutant. The results are consistent with CaSte11 facilitating Hst7-mediated MAP kinase phosphorylation and also playing a potentially critical role in both MAP kinase and Cst5 scaffold localization.

How cells determine the number of polarity sites

The diversity of cell morphologies arises, in part, through regulation of cell polarity by Rho-family GTPases. A poorly understood but fundamental question concerns the regulatory mechanisms by which different cells generate different numbers of polarity sites. Mass-conserved activator-substrate (MCAS) models that describe polarity circuits develop multiple initial polarity sites, but then those sites engage in competition, leaving a single winner. Theoretical analyses predicted that competition would slow dramatically as GTPase concentrations at different polarity sites increase toward a 'saturation point', allowing polarity sites to coexist. Here, we test this prediction using budding yeast cells, and confirm that increasing the amount of key polarity proteins results in multiple polarity sites and simultaneous budding. Further, we elucidate a novel design principle whereby cells can switch from competition to equalization among polarity sites. These findings provide insight into how cells with diverse morphologies may determine the number of polarity sites.

Phosphoproteome Response to Dithiothreitol Reveals Unique Versus Shared Features of Saccharomyces cerevisiae Stress Responses

To cope with sudden changes in the external environment, the budding yeast Saccharomyces cerevisiae orchestrates a multifaceted response that spans many levels of physiology. Several studies have interrogated the transcriptome response to endoplasmic reticulum (ER) stress and the role of regulators such as the Ire1 kinase and Hac1 transcription factors. However, less is known about responses to ER stress at other levels of physiology. Here, we used quantitative phosphoproteomics and computational network inference to uncover the yeast phosphoproteome response to the reducing agent dithiothreitol (DTT) and the upstream signaling network that controls it. We profiled wild-type cells and mutants lacking IRE1 or MAPK kinases MKK1 and MKK2, before and at various times after DTT treatment. In addition to revealing downstream targets of these kinases, our inference approach predicted new regulators in the DTT response, including cell-cycle regulator Cdc28 and osmotic-response kinase Rck2, which we validated computationally. Our results also revealed similarities and surprising differences in responses to different stress conditions, especially in the response of protein kinase A targets. These results have implications for the breadth of signaling programs that can give rise to common stress response signatures.

Mechanistic insights into actin-driven polarity site movement in yeast

Directed cell growth or migration are critical for the development and function of many eukaryotic cells. These cells develop a dynamic "front" (also called "polarity site") that can change direction. Polarity establishment involves autocatalytic accumulation of polarity regulators, including the conserved Rho-family GTPase Cdc42, but the mechanisms underlying polarity reorientation remain poorly understood. The tractable model yeast, Saccharomyces cerevisiae, relocates its polarity site when searching for mating partners. Relocation requires polymerized actin, and is thought to involve actin-mediated vesicle traffic to the polarity site. In this study, we provide a quantitative characterization of spontaneous polarity site movement as a search process and use a mechanistic computational model that combines polarity protein biochemical interactions with vesicle trafficking to probe how various processes might affect polarity site movement. Our findings identify two previously documented features of yeast vesicle traffic as being particularly relevant to such movement: tight spatial focusing of exocytosis enhances the directional persistence of movement, and association of Cdc42-directed GTPase-Activating Proteins with secretory vesicles increases the distance moved. Furthermore, we suggest that variation in the rate of exocytosis beyond simple Poisson dynamics may be needed to fully account for the characteristics of polarity site movement in vivo.

Ratiometric GPCR signaling enables directional sensing in yeast

Accurate detection of extracellular chemical gradients is essential for many cellular behaviors. Gradient sensing is challenging for small cells, which can experience little difference in ligand concentrations on the up-gradient and down-gradient sides of the cell. Nevertheless, the tiny cells of the yeast Saccharomyces cerevisiae reliably decode gradients of extracellular pheromones to find their mates. By imaging the behavior of polarity factors and pheromone receptors, we quantified the accuracy of initial polarization during mating encounters. We found that cells bias the orientation of initial polarity up-gradient, even though they have unevenly distributed receptors. Uneven receptor density means that the gradient of ligand-bound receptors does not accurately reflect the external pheromone gradient. Nevertheless, yeast cells appear to avoid being misled by responding to the fraction of occupied receptors rather than simply the concentration of ligand-bound receptors. Such ratiometric sensing also serves to amplify the gradient of active G protein. However, this process is quite error-prone, and initial errors are corrected during a subsequent indecisive phase in which polarity clusters exhibit erratic mobile behavior.

Cells use surface receptors to decode spatial information from chemical gradients, but accurate decoding is hampered by small cell size and the presence of molecular noise. This study shows that yeast cells decode pheromone gradients by measuring the local ratio of bound to unbound receptors. This mechanism corrects for uneven receptor density at the surface and amplifies the gradient transmitted to downstream components.

Negative Feedback Phosphorylation of Gγ Subunit Ste18 and the Ste5 Scaffold Synergistically Regulates MAPK Activation in Yeast

Heterotrimeric G proteins (Gαβγ) are essential transducers in G protein signaling systems in all eukaryotes. In yeast, G protein signaling differentially activates mitogen-activated protein kinases (MAPKs)—Fus3 and Kss1—a phenomenon controlled by plasma membrane (PM) association of the scaffold protein Ste5. Here, we show that phosphorylation of the yeast Gγ subunit (Ste18), together with Fus3 docking on Ste5, controls the rate and stability of Ste5/PM association. Disruption of either element alone by point mutation has mild but reciprocal effects on MAPK activation. Disabling both elements results in ultra-fast and stable bulk Ste5/PM localization and Fus3 activation that is 6 times faster and 4 times more amplified compared to wild-type cells. These results further resolve the mechanism by which MAPK negative feedback phosphorylation controls pathway activation and provides compelling evidence that Gγ subunits can serve as intrinsic regulators of G protein signaling.

Gene networks that compensate for crosstalk with crosstalk

Crosstalk is a major challenge to engineering sophisticated synthetic gene networks. A common approach is to insulate signal-transduction pathways by minimizing molecular-level crosstalk between endogenous and synthetic genetic components, but this strategy can be difficult to apply in the context of complex, natural gene networks and unknown interactions. Here, we show that synthetic gene networks can be engineered to compensate for crosstalk by integrating pathway signals, rather than by pathway insulation. We demonstrate this principle using reactive oxygen species (ROS)-responsive gene circuits in Escherichia coli that exhibit concentration-dependent crosstalk with non-cognate ROS. We quantitatively map the degree of crosstalk and design gene circuits that introduce compensatory crosstalk at the gene network level. The resulting gene network exhibits reduced crosstalk in the sensing of the two different ROS. Our results suggest that simple network motifs that compensate for pathway crosstalk can be used by biological networks to accurately interpret environmental signals.

Crosstalk between genetic circuits is a major challenge for engineering sophisticated networks. Here the authors design networks that compensate for crosstalk by integrating, not insulating, pathways.

Sharing of Phosphatases Promotes Response Plasticity in Phosphorylation Cascades

Sharing of positive or negative regulators between multiple targets is frequently observed in cellular signaling cascades. For instance, phosphatase sharing between multiple kinases is ubiquitous within the MAPK pathway. Here we investigate how such phosphatase sharing could shape robustness and evolvability of the phosphorylation cascade. Through modeling and evolutionary simulations, we demonstrate that 1) phosphatase sharing dramatically increases robustness of a bistable MAPK response, and 2) phosphatase-sharing cascades evolve faster than nonsharing cascades. This faster evolution is particularly pronounced when evolving from a monostable toward a bistable phenotype, whereas the transition speed of a population from a bistable to monostable response is not affected by phosphatase sharing. This property may enable the phosphatase-sharing design to adapt better in a changing environment. Analysis of the respective mutational landscapes reveal that phosphatase sharing reduces the number of limiting mutations required for transition from monostable to bistable responses, hence facilitating a faster transition to such response types. Taken together, using MAPK cascade as an example, our study offers a general theoretical framework to explore robustness and evolutionary plasticity of signal transduction cascades.

MAPK modulation of yeast pheromone signaling output and the role of phosphorylation sites in the scaffold protein Ste5

Mitogen-activated protein kinases (MAPKs) mediate numerous eukaryotic signaling responses. They also can modulate their own signaling output via positive or negative feedback loops. In the yeast pheromone response pathway, the MAPK Fus3 triggers negative feedback that dampens its own activity. One target of this feedback is Ste5, a scaffold protein that promotes Fus3 activation. Binding of Fus3 to a docking motif (D motif) in Ste5 causes signal dampening, which was proposed to involve a central cluster of phosphorylation sites in Ste5. Here, we reanalyzed the role of these central sites. Contrary to prior claims, phosphorylation-mimicking mutations at these sites did not impair signaling. Also, the hyperactive signaling previously observed when these sites were mutated to nonphosphorylatable residues arose from their replacement with valine residues and was not observed with other substitutes. Instead, a cluster of N-terminal sites in Ste5, not the central sites, is required for the rapid dampening of initial responses. Further results suggest that the role of the Fus3 D motif is most simply explained by a tethering effect that promotes Ste5 phosphorylation, rather than an allosteric effect proposed to regulate Fus3 activity. These findings substantially revise our understanding of how MAPK feedback attenuates scaffold-mediated signaling in this model pathway.

Phenotypic instability in fungi

Fungi are prone to phenotypic instability, that is, the vegetative phase of these organisms, be they yeasts or molds, undergoes frequent switching between two or more behaviors, often with different morphologies, but also sometime having different physiologies without any obvious morphological outcome. In the context of industrial utilization of fungi, this can have a negative impact on the maintenance of strains and/or on their productivity. Instabilities have been shown to result from various mechanisms, either genetic or epigenetic. This chapter will review different types of instabilities and discuss some lesser-known ones, mostly in filamentous fungi, while it will direct readers to additional literature in the case of well-known phenomena such as the amyloid prions or fungal senescence. It will present in depth the "white/opaque" switch of Candida albicans and the "crippled growth" degeneration of the model fungus Podospora anserina. These are two of the most thoroughly studied epigenetic phenotypic switches. I will also discuss the "sectors" presented by many filamentous ascomycetes, for which a prion-based model exists but is not demonstrated. Finally, I will also describe intriguing examples of phenotypic instability for which an explanation has yet to be provided.

Single-cell profiling screen identifies microtubule-dependent reduction of variability in signaling

Populations of isogenic cells often respond coherently to signals, despite differences in protein abundance and cell state. Previously, we uncovered processes in the Saccharomyces cerevisiae pheromone response system (PRS) that reduced cell-to-cell variability in signal strength and cellular response. Here, we screened 1,141 non-essential genes to identify 50 "variability genes". Most had distinct, separable effects on strength and variability of the PRS, defining these quantities as genetically distinct "axes" of system behavior. Three genes affected cytoplasmic microtubule function: BIM1, GIM2, and GIM4 We used genetic and chemical perturbations to show that, without microtubules, PRS output is reduced but variability is unaffected, while, when microtubules are present but their function is perturbed, output is sometimes lowered, but its variability is always high. The increased variability caused by microtubule perturbations required the PRS MAP kinase Fus3 and a process at or upstream of Ste5, the membrane-localized scaffold to which Fus3 must bind to be activated. Visualization of Ste5 localization dynamics demonstrated that perturbing microtubules destabilized Ste5 at the membrane signaling site. The fact that such microtubule perturbations cause aberrant fate and polarity decisions in mammals suggests that microtubule-dependent signal stabilization might also operate throughout metazoans.

Hill kinetics as a noise filter: the role of transcription factor autoregulation in gene cascades

An intuition based on deterministic models of chemical kinetics is that population heterogeneity of transcription factor levels in cells is transmitted unchanged downstream to the target genes. We use a stochastic model of a two-gene cascade with a self-regulating upstream gene to show that, counter to the intuition, there is no simple mapping (bimodal to bimodal, unimodal to unimodal) between the shapes of the distributions of transcription factor numbers and target protein numbers in cells. Due to the presence of the two regulations, the system contains two nonlinear transfer functions, defined by the Hill kinetics of transcription factor binding. The transfer function of the regulator can "interfere" with the transfer function of the target, converting the bimodal input into a unimodal output or vice versa. We show that this effect can be predicted by a geometric construction. As an example application of the method, we present a case study of a system of several downstream genes of different sensitivities, controlled by a common transcription factor which also regulates its own transcription. We show that a single regulator can induce qualitatively different patterns (binary or graded) of responses to a signal in different downstream genes, depending on whether the sensitivity regions of the transfer functions of the upstream and downstream genes overlap or not. Alternatively, the same model can be interpreted as describing a single downstream gene that has different sensitivities in different cell lines due to mutations. Our model shows, therefore, a possible kinetic mechanism by which different genes can interpret the same biological signal in a different manner.

Analysis of the thresholds for transcriptional activation by the yeast MAP kinases Fus3 and Kss1

Signaling in the pheromone response pathway of budding yeast activates two distinct MAP kinases (MAPKs), Fus3 and Kss1. Either MAPK alone can mediate pheromone-induced transcription, but it has been unclear to what degree each one contributes to transcriptional output in wild-type cells. Here, we report that transcription reflects the ratio of active to inactive MAPK, and not simply the level of active MAPK. For Kss1 the majority of MAPK molecules must be converted to the active form, whereas for Fus3 only a small minority must be activated. These different activation thresholds reflect two opposing effects of each MAPK, in which the inactive forms inhibit transcription, whereas the active forms promote transcription. Moreover, negative feedback from Fus3 limits activation of Kss1 so that it does not meet its required threshold in wild-type cells but does so only when hyperactivated in cells lacking Fus3. The results suggest that the normal transcriptional response involves asymmetric contributions from the two MAPKs, in which pheromone signaling reduces the negative effect of Kss1 while increasing the positive effect of Fus3. These findings reveal new functional distinctions between these MAPKs, and help illuminate how inhibitory functions shape positive pathway outputs in both pheromone and filamentation pathways.

Heat stress promotes longevity in budding yeast by relaxing the confinement of age-promoting factors in the mother cell

Although individuals of many species inexorably age, a number of observations established that the rate of aging is modulated in response to a variety of mild stresses. Here, we investigated how heat stress promotes longevity in yeast. We show that upon growth at higher temperature, yeast cells relax the retention of DNA circles, which act as aging factors in the mother cell. The enhanced frequency at which circles redistribute to daughter cells was not due to changes of anaphase duration or nuclear shape but solely to the downregulation of the diffusion barrier in the nuclear envelope. This effect depended on the PKA and Tor1 pathways, downstream of stress-response kinase Pkc1. Inhibition of these responses restored barrier function and circle retention and abrogated the effect of heat stress on longevity. Our data indicate that redistribution of aging factors from aged cells to their progeny can be a mechanism for modulating longevity.

Epigenetic control of pheromone MAPK signaling determines sexual fecundity in Candida albicans

Several pathogenic Candida species are capable of heritable and reversible switching between two epigenetic states, "white" and "opaque." In Candida albicans, white cells are essentially sterile, whereas opaque cells are mating-proficient. Here, we interrogate the mechanism by which the white-opaque switch regulates sexual fecundity and identify four genes in the pheromone MAPK pathway that are expressed at significantly higher levels in opaque cells than in white cells. These genes encode the β subunit of the G-protein complex (STE4), the pheromone MAPK scaffold (CST5), and the two terminal MAP kinases (CEK1/CEK2). To define the contribution of each factor to mating, C. albicans white cells were reverse-engineered to express elevated, opaque-like levels of these factors, either singly or in combination. We show that white cells co-overexpressing STE4, CST5, and CEK2 undergo mating four orders of magnitude more efficiently than control white cells and at a frequency approaching that of opaque cells. Moreover, engineered white cells recapitulate the transcriptional and morphological responses of opaque cells to pheromone. These results therefore reveal multiple bottlenecks in pheromone MAPK signaling in white cells and that alleviation of these bottlenecks enables efficient mating by these "sterile" cell types. Taken together, our findings establish that differential expression of several MAPK factors underlies the epigenetic control of mating in C. albicans We also discuss how fitness advantages could have driven the evolution of a toggle switch to regulate sexual reproduction in pathogenic Candida species.

Robust Synthetic Circuits for Two-Dimensional Control of Gene Expression in Yeast

Cellular phenotypes are the result of complex interactions between many components. Understanding and predicting the system level properties of the resulting networks requires the development of perturbation tools that can simultaneously and independently modulate multiple cellular variables. Here, we develop synthetic modules that use different arrangements of two transcriptional regulators to achieve either concurrent and independent control of the expression of two genes, or decoupled control of the mean and variance of a single gene. These modules constitute powerful tools to probe the quantitative attributes of network wiring and function.

Spatial and temporal signal processing and decision making by MAPK pathways

Mitogen-activated protein kinase (MAPK) pathways are conserved from yeast to man and regulate a variety of cellular processes, including proliferation and differentiation. Recent developments show how MAPK pathways perform exquisite spatial and temporal signal processing and underscores the importance of studying the dynamics of signaling pathways to understand their physiological response. The importance of dynamic mechanisms that process input signals into graded downstream responses has been demonstrated in the pheromone-induced and osmotic stress-induced MAPK pathways in yeast and in the mammalian extracellular signal-regulated kinase MAPK pathway. Particularly, recent studies in the yeast pheromone response have shown how positive feedback generates switches, negative feedback enables gradient detection, and coherent feedforward regulation underlies cellular memory. More generally, a new wave of quantitative single-cell studies has begun to elucidate how signaling dynamics determine cell physiology and represents a paradigm shift from descriptive to predictive biology.

Single-cell dynamics and variability of MAPK activity in a yeast differentiation pathway

In response to pheromones, yeast cells activate a MAPK pathway to direct processes important for mating, including gene induction, cell-cycle arrest, and polarized cell growth. Although a variety of assays have been able to elucidate signaling activities at multiple steps in the pathway, measurements of MAPK activity during the pheromone response have remained elusive, and our understanding of single-cell signaling behavior is incomplete. Using a yeast-optimized FRET-based mammalian Erk-activity reporter to monitor Fus3 and Kss1 activity in live yeast cells, we demonstrate that overall mating MAPK activity exhibits distinct temporal dynamics, rapid reversibility, and a graded dose dependence around the K_D of the receptor, where phenotypic transitions occur. The complex dose response was found to be largely a consequence of two feedbacks involving cyclin-mediated scaffold phosphorylation and Fus3 autoregulation. Distinct cell cycle-dependent response patterns comprised a large portion of the cell-to-cell variability at each dose, constituting the major source of extrinsic noise in coupling activity to downstream gene-expression responses. Additionally, we found diverse spatial MAPK activity patterns to emerge over time in cells undergoing default, gradient, and true mating responses. Furthermore, ramping up and rapid loss of activity were closely associated with zygote formation in mating-cell pairs, supporting a role for elevated MAPK activity in successful cell fusion and morphogenic reorganization. Altogether, these findings present a detailed view of spatiotemporal MAPK activity during the pheromone response, elucidating its role in mediating complex long-term developmental fates in a unicellular differentiation system.

Sensory input attenuation allows predictive sexual response in yeast

Animals are known to adjust their sexual behaviour depending on mate competition. Here we report similar regulation for mating behaviour in a sexual unicellular eukaryote, the budding yeast Saccharomyces cerevisiae. We demonstrate that pheromone-based communication between the two mating types, coupled to input attenuation by recipient cells, enables yeast to robustly monitor relative mate abundance (sex ratio) within a mixed population and to adjust their commitment to sexual reproduction in proportion to their estimated chances of successful mating. The mechanism of sex-ratio sensing relies on the diffusible peptidase Bar1, which is known to degrade the pheromone signal produced by mating partners. We further show that such a response to sexual competition within a population can optimize the fitness trade-off between the costs and benefits of mating response induction. Our study thus provides an adaptive explanation for the known molecular mechanism of pheromone degradation in yeast.

Cells of the yeast Saccharomyces cerevisiae can mate with other cells of opposite mating type. Here, the authors show that the combination of a pheromone and a pheromone-degrading enzyme allows yeast cells to monitor relative mate abundance within a population and adjust their commitment to sexual reproduction.

Switch-like Transitions Insulate Network Motifs to Modularize Biological Networks

Cellular decisions are made by complex networks that are difficult to analyze. Although it is common to analyze smaller sub-networks known as network motifs, it is unclear whether this is valid, because these motifs are embedded in complex larger networks. Here, we address the general question of modularity by examining the S. cerevisiae pheromone response. We demonstrate that the feedforward motif controlling the cell-cycle inhibitor Far1 is insulated from cell-cycle dynamics by the positive feedback switch that drives reentry to the cell cycle. Before cells switch on positive feedback, the feedforward motif model predicts the behavior of the larger network. Conversely, after the switch, the feedforward motif is dismantled and has no discernable effect on the cell cycle. When insulation is broken, the feedforward motif no longer predicts network behavior. This work illustrates how, despite the interconnectivity of networks, the activity of motifs can be insulated by switches that generate well-defined cellular states.

Synthetic mixed-signal computation in living cells

Living cells implement complex computations on the continuous environmental signals that they encounter. These computations involve both analogue- and digital-like processing of signals to give rise to complex developmental programs, context-dependent behaviours and homeostatic activities. In contrast to natural biological systems, synthetic biological systems have largely focused on either digital or analogue computation separately. Here we integrate analogue and digital computation to implement complex hybrid synthetic genetic programs in living cells. We present a framework for building comparator gene circuits to digitize analogue inputs based on different thresholds. We then demonstrate that comparators can be predictably composed together to build band-pass filters, ternary logic systems and multi-level analogue-to-digital converters. In addition, we interface these analogue-to-digital circuits with other digital gene circuits to enable concentration-dependent logic. We expect that this hybrid computational paradigm will enable new industrial, diagnostic and therapeutic applications with engineered cells.

Digital and analogue gene circuits each have distinct advantages in natural and engineered cells. Here, Rubens et al. engineer synthetic gene circuits that implement mixed-signal digital and analogue computations in living cells.

Insights from synthetic yeasts

Synthetic biology is one of the most exciting strategies for the investigation of living organisms and lies at the intersection of biology and engineering. Originally developed in prokaryotes, the idea of deciphering biological phenomena through building artificial genetic circuits and studying their behaviours has rapidly demonstrated its potential in a broad range of fields in the life sciences. From the assembly of synthetic genomes to the generation of novel biological functions, yeast cells have imposed themselves as the most powerful eukaryotic model for this approach. However, we are only beginning to explore the possibilities of synthetic biology, and the perspectives it offers in a genetically amenable system such as yeasts are endless. Copyright © 2016 John Wiley & Sons, Ltd.

Sensing a bud in the yeast morphogenesis checkpoint: a role for Elm1

Yeast cells know whether or not they have a bud. The kinase Elm1 and the septin cytoskeleton are key transducers of cell shape information.

Bud formation by Saccharomyces cerevisiae must be coordinated with the nuclear cycle to enable successful proliferation. Many environmental stresses temporarily disrupt bud formation, and in such circumstances, the morphogenesis checkpoint halts nuclear division until bud formation can resume. Bud emergence is essential for degradation of the mitotic inhibitor, Swe1. Swe1 is localized to the septin cytoskeleton at the bud neck by the Swe1-binding protein Hsl7. Neck localization of Swe1 is required for Swe1 degradation. Although septins form a ring at the presumptive bud site before bud emergence, Hsl7 is not recruited to the septins until after bud emergence, suggesting that septins and/or Hsl7 respond to a “bud sensor.” Here we show that recruitment of Hsl7 to the septin ring depends on a combination of two septin-binding kinases: Hsl1 and Elm1. We elucidate which domains of these kinases are needed and show that artificial targeting of those domains suffices to recruit Hsl7 to septin rings even in unbudded cells. Moreover, recruitment of Elm1 is responsive to bud emergence. Our findings suggest that Elm1 plays a key role in sensing bud emergence.

Polarity establishment requires localized activation of Cdc42

Positive feedback by localized activation of Cdc42 drives polarity establishment in the budding yeast Saccharomyces cerevisiae.

Establishment of cell polarity in animal and fungal cells involves localization of the conserved Rho-family guanosine triphosphatase, Cdc42, to the cortical region destined to become the “front” of the cell. The high local concentration of active Cdc42 promotes cytoskeletal polarization through various effectors. Cdc42 accumulation at the front is thought to involve positive feedback, and studies in the budding yeast Saccharomyces cerevisiae have suggested distinct positive feedback mechanisms. One class of mechanisms involves localized activation of Cdc42 at the front, whereas another class involves localized delivery of Cdc42 to the front. Here we show that Cdc42 activation must be localized for successful polarity establishment, supporting local activation rather than local delivery as the dominant mechanism in this system.

Mapping regions in Ste5 that support Msn5-dependent and -independent nuclear export

Careful control of the available pool of the MAPK scaffold Ste5 is important for mating-pathway activation and the prevention of inappropriate mating differentiation in haploid Saccharomyces cerevisiae. Ste5 shuttles constitutively through the nucleus, where it is degraded by a ubiquitin-dependent mechanism triggered by G1 CDK phosphorylation. Here we narrow-down regions of Ste5 that mediate nuclear export. Four regions in Ste5 relocalize SV40-TAgNLS-GFP-GFP from nucleus to cytoplasm. One region is N-terminal, dependent on exportin Msn5/Ste21/Kap142, and interacts with Msn5 in 2 hybrid assays independently of mating pheromone, Fus3, Kss1, Ptc1, the NLS/PM, and RING-H2. A second region overlaps the PH domain and Ste11 binding site and 2 others are on the vWA domain and include residues essential for MAPK activation. We find no evidence for dependence on Crm1/Xpo1, despite numerous potential nuclear export sequences (NESs) detected by LocNES and NetNES1.1 predictors. Thus, Msn5 (homolog of human Exportin-5) and one or more exportins or adaptor molecules besides Crm1/Xpo1 may regulate Ste5 through multiple recognition sites.

Different Mechanisms Confer Gradual Control and Memory at Nutrient- and Stress-Regulated Genes in Yeast

Cells respond to environmental stimuli by fine-tuned regulation of gene expression. Here we investigated the dose-dependent modulation of gene expression at high temporal resolution in response to nutrient and stress signals in yeast. The GAL1 activity in cell populations is modulated in a well-defined range of galactose concentrations, correlating with a dynamic change of histone remodeling and RNA polymerase II (RNAPII) association. This behavior is the result of a heterogeneous induction delay caused by decreasing inducer concentrations across the population. Chromatin remodeling appears to be the basis for the dynamic GAL1 expression, because mutants with impaired histone dynamics show severely truncated dose-response profiles. In contrast, the GRE2 promoter operates like a rapid off/on switch in response to increasing osmotic stress, with almost constant expression rates and exclusively temporal regulation of histone remodeling and RNAPII occupancy. The Gal3 inducer and the Hog1 mitogen-activated protein (MAP) kinase seem to determine the different dose-response strategies at the two promoters. Accordingly, GAL1 becomes highly sensitive and dose independent if previously stimulated because of residual Gal3 levels, whereas GRE2 expression diminishes upon repeated stimulation due to acquired stress resistance. Our analysis reveals important differences in the way dynamic signals create dose-sensitive gene expression outputs.

How MAP kinase modules function as robust, yet adaptable, circuits

Genetic and biochemical studies have revealed that the diversity of cell types and developmental patterns evident within the animal kingdom is generated by a handful of conserved, core modules. Core biological modules must be robust, able to maintain functionality despite perturbations, and yet sufficiently adaptable for random mutations to generate phenotypic variation during evolution. Understanding how robust, adaptable modules have influenced the evolution of eukaryotes will inform both evolutionary and synthetic biology. One such system is the MAP kinase module, which consists of a 3-tiered kinase circuit configuration that has been evolutionarily conserved from yeast to man. MAP kinase signal transduction pathways are used across eukaryotic phyla to drive biological functions that are crucial for life. Here we ask the fundamental question, why do MAPK modules follow a conserved 3-tiered topology rather than some other number? Using computational simulations, we identify a fundamental 2-tiered circuit topology that can be readily reconfigured by feedback loops and scaffolds to generate diverse signal outputs. When this 2-kinase circuit is connected to proximal input kinases, a 3-tiered modular configuration is created that is both robust and adaptable, providing a biological circuit that can regulate multiple phenotypes and maintain functionality in an uncertain world. We propose that the 3-tiered signal transduction module has been conserved through positive selection, because it facilitated the generation of phenotypic variation during eukaryotic evolution.

Signal inhibition by a dynamically regulated pool of monophosphorylated MAPK

MAPKs are activated by dual phosphorylation. Much of the MAPK Fus3 is monophosphorylated and acts to inhibit signaling in vivo. Computational models reveal how a kinase scaffold and phosphatase act together to dynamically regulate dual-phosphorylated and monophosphorylated MAPKs and the downstream signal.

Protein kinases regulate a broad array of cellular processes and do so through the phosphorylation of one or more sites within a given substrate. Many protein kinases are themselves regulated through multisite phosphorylation, and the addition or removal of phosphates can occur in a sequential (processive) or a stepwise (distributive) manner. Here we measured the relative abundance of the monophosphorylated and dual-phosphorylated forms of Fus3, a member of the mitogen-activated protein kinase (MAPK) family in yeast. We found that upon activation with pheromone, a substantial proportion of Fus3 accumulates in the monophosphorylated state. Introduction of an additional copy of Fus3 lacking either phosphorylation site leads to dampened signaling. Conversely, cells lacking the dual-specificity phosphatase (msg5Δ) or that are deficient in docking to the MAPK-scaffold (Ste5^ND) accumulate a greater proportion of dual-phosphorylated Fus3. The double mutant exhibits a synergistic, or “synthetic,” supersensitivity to pheromone. Finally, we present a predictive computational model that combines MAPK scaffold and phosphatase activities and is sufficient to account for the observed MAPK profiles. These results indicate that the monophosphorylated and dual-phosphorylated forms of the MAPK act in opposition to one another. Moreover, they reveal a new mechanism by which the MAPK scaffold acts dynamically to regulate signaling.

Evolution of synthetic signaling scaffolds by recombination of modular protein domains

Signaling scaffolds are proteins that interact via modular domains with multiple partners, regulating signaling networks in space and time and providing an ideal platform from which to alter signaling functions. However, to better exploit scaffolds for signaling engineering, it is necessary to understand the full extent of their modularity. We used a directed evolution approach to identify, from a large library of randomly shuffled protein interaction domains, variants capable of rescuing the signaling defect of a yeast strain in which Ste5, the scaffold in the mating pathway, had been deleted. After a single round of selection, we identified multiple synthetic scaffold variants with diverse domain architectures, able to mediate mating pathway activation in a pheromone-dependent manner. The facility with which this signaling network accommodates changes in scaffold architecture suggests that the mating signaling complex does not possess a single, precisely defined geometry into which the scaffold has to fit. These relaxed geometric constraints may facilitate the evolution of signaling networks, as well as their engineering for applications in synthetic biology.

Compartmentalization of a bistable switch enables memory to cross a feedback-driven transition

Cells make accurate decisions in the face of molecular noise and environmental fluctuations by relying not only on present pathway activity, but also on their memory of past signaling dynamics. Once a decision is made, cellular transitions are often rapid and switch-like due to positive feedback loops in the regulatory network. While positive feedback loops are good at promoting switch-like transitions, they are not expected to retain information to inform subsequent decisions. However, this expectation is based on our current understanding of network motifs that accounts for temporal, but not spatial, dynamics. Here, we show how spatial organization of the feedback-driven yeast G1/S switch enables the transmission of memory of past pheromone exposure across this transition. We expect this to be one of many examples where the exquisite spatial organization of the eukaryotic cell enables previously well-characterized network motifs to perform new and unexpected signal processing functions.

MAPK feedback encodes a switch and timer for tunable stress adaptation in yeast

Signaling pathways can behave as switches or rheostats, generating binary or graded responses to a given cell stimulus. We evaluated whether a single signaling pathway can simultaneously encode a switch and a rheostat. We found that the kinase Hog1 mediated a bifurcated cellular response: Activation and commitment to adaptation to osmotic stress are switchlike, whereas protein induction and the resolution of this commitment are graded. Through experimentation, bioinformatics analysis, and computational modeling, we determined that graded recovery is encoded through feedback phosphorylation and a gene induction program that is both temporally staggered and variable across the population. This switch-to-rheostat signaling mechanism represents a versatile stress adaptation system, wherein a broad range of inputs generate an "all-in" response that is later tuned to allow graded recovery of individual cells over time.

Plasma membrane aminoglycerolipid flippase function is required for signaling competence in the yeast mating pheromone response pathway

The class 4 P-type ATPases ("flippases") maintain membrane asymmetry by translocating phosphatidylethanolamine and phosphatidylserine from the outer leaflet to the cytosolic leaflet of the plasma membrane. In Saccharomyces cerevisiae, five related gene products (Dnf1, Dnf2, Dnf3, Drs2, and Neo1) are implicated in flipping of phosphatidylethanolamine, phosphatidylserine, and phosphatidylcholine. In MAT A: cells responding to α-factor, we found that Dnf1, Dnf2, and Dnf3, as well as the flippase-activating protein kinase Fpk1, localize at the projection ("shmoo") tip where polarized growth is occurring and where Ste5 (the central scaffold protein of the pheromone-initiated MAPK cascade) is recruited. Although viable, a MAT A: dnf1∆ dnf2∆ dnf3∆ triple mutant exhibited a marked decrease in its ability to respond to α-factor, which we could attribute to pronounced reduction in Ste5 stability resulting from an elevated rate of its Cln2⋅Cdc28-initiated degradation. Similarly, a MAT A: dnf1∆ dnf3∆ drs2∆ triple mutant also displayed marked reduction in its ability to respond to α-factor, which we could attribute to inefficient recruitment of Ste5 to the plasma membrane due to severe mislocalization of the cellular phosphatidylinositol 4-phosphate and phosphatidylinositol 4,5-bisphosphate pools. Thus proper remodeling of plasma membrane aminoglycerolipids and phosphoinositides is necessary for efficient recruitment, stability, and function of the pheromone signaling apparatus.

Cdc42p-interacting protein Bem4p regulates the filamentous-growth mitogen-activated protein kinase pathway

The ubiquitous Rho (Ras homology) GTPase Cdc42p can function in different settings to regulate cell polarity and cellular signaling. How Cdc42p and other proteins are directed to function in a particular context remains unclear. We show that the Cdc42p-interacting protein Bem4p regulates the mitogen-activated protein kinase (MAPK) pathway that controls filamentous growth in Saccharomyces cerevisiae. Bem4p controlled the filamentous-growth pathway but not other MAPK pathways (mating or high-osmolarity glycerol response [HOG]) that also require Cdc42p and other shared components. Bem4p associated with the plasma membrane (PM) protein, Sho1p, to regulate MAPK activity and cell polarization under nutrient-limiting conditions that favor filamentous growth. Bem4p also interacted with the major activator of Cdc42p, the guanine nucleotide exchange factor (GEF) Cdc24p, which we show also regulates the filamentous-growth pathway. Bem4p interacted with the pleckstrin homology (PH) domain of Cdc24p, which functions in an autoinhibitory capacity, and was required, along with other pathway regulators, to maintain Cdc24p at polarized sites during filamentous growth. Bem4p also interacted with the MAPK kinase kinase (MAPKKK) Ste11p. Thus, Bem4p is a new regulator of the filamentous-growth MAPK pathway and binds to general proteins, like Cdc42p and Ste11p, to promote a pathway-specific response.

Inhibitory GEF phosphorylation provides negative feedback in the yeast polarity circuit

Cell polarity is critical for the form and function of many cell types. During polarity establishment, cells define a cortical "front" that behaves differently from the rest of the cortex. The front accumulates high levels of the active form of a polarity-determining Rho-family GTPase (Cdc42, Rac, or Rop) that then orients cytoskeletal elements through various effectors to generate the polarized morphology appropriate to the particular cell type [1, 2]. GTPase accumulation is thought to involve positive feedback, such that active GTPase promotes further delivery and/or activation of more GTPase in its vicinity [3]. Recent studies suggest that once a front forms, the concentration of polarity factors at the front can increase and decrease periodically, first clustering the factors at the cortex and then dispersing them back to the cytoplasm [4-7]. Such oscillatory behavior implies the presence of negative feedback in the polarity circuit [8], but the mechanism of negative feedback was not known. Here we show that, in the budding yeast Saccharomyces cerevisiae, the catalytic activity of the Cdc42-directed GEF is inhibited by Cdc42-stimulated effector kinases, thus providing negative feedback. We further show that replacing the GEF with a phosphosite mutant GEF abolishes oscillations and leads to the accumulation of excess GTP-Cdc42 and other polarity factors at the front. These findings reveal a mechanism for negative feedback and suggest that the function of negative feedback via GEF inhibition is to buffer the level of Cdc42 at the polarity site.

Scaffolding proteins: not such innocent bystanders

Sequential transfer of information from one enzyme to the next within the confines of a protein kinase scaffold enhances signal transduction. Though frequently considered to be inert organizational elements, two recent reports implicate kinase-scaffolding proteins as active participants in signal relay.

Synthetic biology tools for programming gene expression without nutritional perturbations in Saccharomyces cerevisiae

A conditional gene expression system that is fast-acting, is tunable and achieves single-gene specificity was recently developed for yeast. A gene placed directly downstream of a modified GAL1 promoter containing six Zif268 binding sequences (with single nucleotide spacing) was shown to be selectively inducible in the presence of β-estradiol, so long as cells express the artificial transcription factor, Z₃EV (a fusion of the Zif268 DNA binding domain, the ligand binding domain of the human estrogen receptor and viral protein 16). We show the strength of Z₃EV-responsive promoters can be modified using straightforward design principles. By moving Zif268 binding sites toward the transcription start site, expression output can be nearly doubled. Despite the reported requirement of estrogen receptor dimerization for hormone-dependent activation, a single binding site suffices for target gene activation. Target gene expression levels correlate with promoter binding site copy number and we engineer a set of inducible promoter chassis with different input–output characteristics. Finally, the coupling between inducer identity and gene activation is flexible: the ligand specificity of Z₃EV can be re-programmed to respond to a non-hormone small molecule with only five amino acid substitutions in the human estrogen receptor domain, which may prove useful for industrial applications.

Identification of the molecular mechanisms for cell-fate selection in budding yeast through mathematical modeling

The specification and maintenance of cell fates is essential to the development of multicellular organisms. However, the precise molecular mechanisms in cell fate selection are, to our knowledge, poorly understood due to the complexity of multiple interconnected pathways. In this study, model-based quantitative analysis is used to explore how to maintain distinguished cell fates between cell-cycle commitment and mating arrest in budding yeast. We develop a full mathematical model of an interlinked regulatory network based on the available experimental data. By theoretically defining the Start transition point, the model is able to reproduce many experimental observations of the dynamical behaviors in wild-type cells as well as in Ste5-8A and Far1-S87A mutants. Furthermore, we demonstrate that a moderate ratio between Cln1/2→Far1 inhibition and Cln1/2→Ste5 inhibition is required to ensure a successful switch between different cell fates. We also show that the different ratios of the mutual Cln1/2 and Far1 inhibition determine the different cell fates. In addition, based on a new, definition of network entropy, we find that the Start point in wild-type cells coincides with the system's point of maximum entropy. This result indicates that Start is a transition point in the network entropy. Therefore, we theoretically explain the Start point from a network dynamics standpoint. Moreover, we analyze the biological bistablity of our model through bifurcation analysis. We find that the Cln1/2 and Cln3 production rates and the nonlinearity of SBF regulation on Cln1/2 production are potential determinants for irreversible entry into a new cell fate. Finally, the quantitative computations further reveal that high specificity and fidelity of the cell-cycle and mating pathways can guarantee specific cell-fate selection. These findings show that quantitative analysis and simulations with a mathematical model are useful tools for understanding the molecular mechanisms in cell-fate decisions.

Multisite phosphorylation networks as signal processors for Cdk1

The order and timing of cell-cycle events is controlled by changing substrate specificity and different activity thresholds of cyclin-dependent kinases (CDKs). However, it is not understood how a single protein kinase can trigger hundreds of switches in a sufficiently time-resolved fashion. We show that cyclin-Cdk1-Cks1-dependent phosphorylation of multisite targets in Saccharomyces cerevisiae is controlled by key substrate parameters including distances between phosphorylation sites, distribution of serines and threonines as phosphoacceptors and positioning of cyclin-docking motifs. The component mediating the key interactions in this process is Cks1, the phosphoadaptor subunit of the cyclin-Cdk1-Cks1 complex. We propose that variation of these parameters within networks of phosphorylation sites in different targets provides a wide range of possibilities for differential amplification of Cdk1 signals, thus providing a mechanism to generate a wide range of thresholds in the cell cycle.

Genetically engineered transvestites reveal novel mating genes in budding yeast

Haploid budding yeast has two mating types, defined by the alleles of the MAT locus, MATa and MATα. Two haploid cells of opposite mating types mate by signaling to each other using reciprocal pheromones and receptors, polarizing and growing toward each other, and eventually fusing to form a single diploid cell. The pheromones and receptors are necessary and sufficient to define a mating type, but other mating-type-specific proteins make mating more efficient. We examined the role of these proteins by genetically engineering "transvestite" cells that swap the pheromone, pheromone receptor, and pheromone processing factors of one mating type for another. These cells mate with each other, but their mating is inefficient. By characterizing their mating defects and examining their transcriptomes, we found Afb1 (a-factor barrier), a novel MATα-specific protein that interferes with a-factor, the pheromone secreted by MATa cells. Strong pheromone secretion is essential for efficient mating, and the weak mating of transvestites can be improved by boosting their pheromone production. Synthetic biology can characterize the factors that control efficiency in biological processes. In yeast, selection for increased mating efficiency is likely to have continually boosted pheromone levels and the ability to discriminate between partners who make more and less pheromone. This discrimination comes at a cost: weak mating in situations where all potential partners make less pheromone.

Msn2 coordinates a stoichiometric gene expression program

BACKGROUND

Many cellular processes operate in an "analog" regime in which the magnitude of the response is precisely tailored to the intensity of the stimulus. In order to maintain the coherence of such responses, the cell must provide for proportional expression of multiple target genes across a wide dynamic range of induction states. Our understanding of the strategies used to achieve graded gene regulation is limited.

RESULTS

In this work, we document a relationship between stress-responsive gene expression and the transcription factor Msn2 that is graded over a large range of Msn2 concentrations. We use computational modeling and in vivo and in vitro analyses to dissect the roots of this relationship. Our studies reveal a simple and general strategy based on noncooperative low-affinity interactions between Msn2 and its cognate binding sites as well as competition over a large number of Msn2 binding sites in the genome relative to the number of Msn2 molecules.

CONCLUSIONS

In addition to enabling precise tuning of gene expression to the state of the environment, this strategy ensures colinear activation of target genes, allowing for stoichiometric expression of large groups of genes without extensive promoter tuning. Furthermore, such a strategy enables precise modulation of the activity of any given promoter by addition of binding sites without altering the qualitative relationship between different genes in a regulon. This feature renders a given regulon highly "evolvable."